Test of the hadronic interaction models at around *10 TeV with Tibet EAS core data

Tibet AsγCollaboration

DOI: https://doi.org/10.48550/arXiv.1303.3094

2013-03-13

Abstract:A hybrid experiment has been started by AS{\gamma} collaboration at Tibet, China, since May 2009, that consists of a burst-detector-grid (YAC, Yangbajing Air shower Core array) and the Tibet air-shower array (Tibet-III). The first step of YAC, called YAC-I, contains 16 detector units and observes high energy electromagnetic particles in air-shower cores within several meters from the shower axis, and Tibet-III array measures the total energy and the arrival direction of air showers. YAC-I is used to check hadronic interaction models currently used for air-shower simulations such as QGSJET, SIBYLL, EPOS etc. through the multi-parameter measurement in air-shower cores. In this paper, we used a data set collected from May 1st 2009 through February 23rd 2010 by the YAC-I. The effective live time used for the present analysis is 169.65 days. The preliminary results of the interaction model checking at *10 TeV energy region is reported in this paper.

High Energy Physics - Experiment,High Energy Astrophysical Phenomena

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to test the accuracy of hadronic interaction models, such as QGSJET, SIBYLL and EPOS models, in the simulation of extensive air showers (EAS) in the energy region of about 10 TeV. Specifically, the research team used the data collected by the Yangbajing Air Shower Core Array (YAC - I) and Tibet - III air shower arrays in Tibet, China, to verify whether these models can accurately describe the behavior of secondary particles produced after high - energy cosmic rays collide with atmospheric molecules.

### Main research contents include:

1. **Experimental device**: The YAC - I experiment consists of 16 scintillator detectors. Each detector contains a 3.5 - centimeter - thick lead layer (approximately equal to 7 radiation lengths) and a 40cm x 50cm x 1cm plastic scintillator. These detectors are used to record high - energy electromagnetic particles in the air shower core.

2. **Data collection**: From May 1, 2009 to February 23, 2010, YAC - I collected data with an effective live time of 169.65 days.

3. **Monte Carlo simulation**: Using the CORSIKA simulation software (version 6.204), combined with the QGSJET2 and SIBYLL2.1 hadronic interaction models, as well as two primary composition models (the heavy - element - dominated model HD and the non - linear acceleration model NLA), air shower events were generated. The primary energy range from 1 TeV to infinity and the incident angle from 0 to 60 degrees were considered in the simulation.

4. **Data analysis**: By selecting different thresholds \( N_b \) and event selection conditions, the research team obtained samples in different energy intervals. The distribution of \( N_b / \sum N_b \) was mainly analyzed and compared with the Monte Carlo simulation results.

### Key findings:

- **Consistency of distribution shape**: In all four cases (QGSJET2 + HD, SIBYLL2.1 + HD, QGSJET2 + NLA, SIBYLL2.1 + NLA), the shape of the \( N_b / \sum N_b \) distribution is consistent with the experimental data, indicating that these models may correctly reflect the actual situation in the 10 TeV energy region with an error of about 10%.

- **Difference in absolute intensity**: In all cases, the absolute intensity of the experimental data is generally higher than that of the Monte Carlo simulation results. Among them, the results of the SIBYLL2.1 + NLA model are the closest to the experimental data, while the QGSJET2 + HD model has the largest difference.

- **Influence of composition models**: For a given hadronic interaction model, the NLA composition model performs better than the HD composition model. For a given composition model, the SIBYLL2.1 model performs better than the QGSJET2 model.

- **Influence of helium spectra**: Considering the new results of the PAMELA and CREAM experiments, both the NLA and HD models underestimate the number of helium nuclei, especially in the relevant energy region. If this factor is considered, the results of the QGSJET2 model may be closer to the experimental data.

### Conclusion:

The YAC - I experiment shows the ability and sensitivity to check hadronic interaction models in the 10 TeV energy region. Although there are some differences in absolute intensity, in general, the experimental data and the Monte Carlo simulation results show good consistency in the distribution shape. Further research is underway to more accurately understand the applicability of these models in high - energy cosmic ray physics.