Tibet Hybrid Experiment Research Articles

A Monte Carlo study to check muon numbers in hadronic interaction models observed by the Tibet hybrid experiment (YAC-II + Tibet-III + MD)

The interpretation of the EAS data relies on the employment of hadronic interaction models, which are subject to theoretical and experimental uncertainties that may hamper composition studies of cosmic rays. To check the reliability of such models at the energies relevant for EAS studies, the predictions of the models can be compared with data from air shower observatories. In this regard, the study of the number of muons becomes extremely useful, since they are sensitive to the hadronic interactions that occur in the early phases of the EAS development. In the paper, we propose a Monte Carlo study to test the number of muons in hadronic interaction models by the hybrid experiment (YAC-II + Tibet-III + MD). For an air-shower event, the Tibet air-shower array (Tibet-III) provides the arrival direction and the air-shower size which are interrelated to primary energy, the Yangbajing Air shower Core detector (YAC-II) array measures the high energy electromagnetic particles in the very forward region so as to obtain the characteristic parameters of air-shower cores, at the same time, the underground MDs record the number of high-energy muons above 1 GeV. Since we can select proton events with high accuracy by YAC-II almost independently of the hadronic interaction models, the accompanying number of muons induced by proton events can be fed out. With the unique settle of YAC-II, our results show that the description of muon numbers in different hadronic interaction models can be well systematics-checked to avoid the ambiguity of the primary cosmic-ray mass composition around the knee energy region by the Tibet hybrid experiment (YAC-II + Tibet-III + MD).

Towards a better estimation of energy and species of primary cosmic rays in the knee region with the Tibet hybrid experiment: Utilization of refined EAS lateral distributions

A new hybrid experiment (YAC-II+Tibet-III+MD) located at Yangbajing has started to improve the capability of explicit measurements of the cosmic-ray components (P, He, CNO, Fe, etc) at the knee energy region since 2014. Considering the different features of air shower development for different primary cosmic-ray nuclei, using a full Monte Carlo simulation, we modified the lateral distribution functions for proton-induced, helium-induced and iron-induced air showers, respectively. The results show that the air shower size of different nuclei obtained by the modified Nishimura-Kamata-Greisen function is consistent with the true shower size within 5% systematic errors. Furthermore, we studied the Extensive Air Shower age parameter determination from the lateral distribution functions of charged particles near the air shower core with the Tibet hybrid experiment. The derived age parameter can be well used to estimate the particle type of the incident primary cosmic rays.

The cosmic ray energy spectrum measured with the new Tibet hybrid experiment

We have upgraded the new Tibet ASgamma experiment in China since 2014 to measure the chemical composition of cosmic rays around the knee. This hybrid experiment consist of an air-shower-core detector array (YAC-II) to detect high energy electromagnetic component, the Tibet air-shower array (Tibet-III) and a large underground water-Cherenkov muon-detector array (MD). We have carried out a detailed air-shower Monte Carlo (MC) simulation to study the performance of the hybrid detectors by using CORSIKA (version 7.5000), which includes EPOS-LHC, QGSJETII-04, SIBYLL2.1 and SIBYLL2.3 hadronic interaction models. The preliminary results of the interaction model checking above 50 TeV energy region are reported in this paper, and the primary proton and helium spectra in the energy range 50 TeV to 1015 eV was derived from YAC-I data and is smoothly connected with direct observation data at lower energies and also with our previously reported works at higher energies within statistical errors. The knee of the (P+He) spectra is located around 400 TeV. The interaction model dependence in deriving the primary (P+He) spectra is found to be small (less than 25% in absolute intensity, 10% in position of the knee), and the composition model dependence is less than 10% in absolute intensity.

Performance of the Tibet hybrid experiment (YAC-II + Tibet-III + MD) to measure the energy spectra of the light primary cosmic rays at energies 50–10,000 TeV

A new hybrid detector system has been constructed by the Tibet ASγ collaboration at Tibet, China, since 2014 to measure the chemical composition of cosmic rays around the knee in the wide energy range. They consist of an air-shower-core detector-grid (YAC-II) to detect high energy electromagnetic component, the Tibet air-shower array (Tibet-III) and a large underground water-Cherenkov muon-detector array (MD). We have carried out a detailed air-shower Monte Carlo (MC) simulation to study the performance of the hybrid detectors by using CORSIKA (version 6.204), which includes QGSJET01c and SIBYLL2.1 hadronic interaction models. Assumed primary cosmic ray models are based on helium poor, helium rich and Gaisser’s fit compositions around the knee. All detector responses are calculated using Geant4 (version 9.5) according to the real detector configurations and the MC events are reconstructed by the same procedure as the experimental data analysis. The energy determination is made by lateral density fitting (LDF) method using modified NKG function and the separation of the light components (proton, helium) is made by means of the artificial neural network (ANN) method and the random forest (RF) method. The systematic errors of the spectra of proton and helium caused by each steps of the analysis procedure are investigated including the dependence of the MC data on the hadronic interaction models and the primary composition models, and the algorithms for the primary mass identification. The systematic errors of the flux to be obtained by the new experiment are summarized as less than 30% in total. Our results show that the new hybrid experiment is powerful enough to study the chemical composition of the cosmic rays, in particular, to obtain the light-component spectra of the primary cosmic rays in 50–10,000TeV energy range overlapping to the direct observation data at low energy side and ground-based indirect observations at high energy side. It is possible in this energy range to find the break points of the power indices of proton and helium (the knee of individual component spectrum) which are basically important parameter for the study of the cosmic-ray origin.

The complex EAS hybrid arrays in Tibet

The Tibet hybrid experiment is an extensive air shower experiment that operated at Yangbajing (E90°31', N30°06'; 4300 m above sea level) in Tibet, China, since 1996. This multi-detector system is used to observe air showers of high energy celestial gamma-ray origin and that of nuclear-component origin. Both observations aim to investigate the origin of cosmic rays using different methods. In this paper, based on the energy spectrum of the nuclear component and its chemical composition around the knee, we would like to discuss whether they could be described by the theoretical picture of the diffusive shock acceleration at SNRs. The consistency of the experimental data with rigidity dependent cutoff spectrum of CRs is examined using Tibet data and KASCADE data. It is found that both experimental data cannot be fully interpreted by such a simple approach. The sharp knee observed by Tibet experiment suggests the contribution of nearby source(s).

Cosmic ray data and their interpretation: about the Tibet Shower Array: Chemical Composition of Cosmic Rays at the Knee Viewed from the Tibet Air Shower Experiment

A brief summary is given about present status of composition measurement of cosmic rays by direct observations and indirect ones which provides key information on the origin of the cosmic rays. The result of Tibet hybrid experiment measuring proton and helium spectra at the knee is presented, in which the decreasing fractions of both of the proton and helium component to the all particle spectrum with increasing energy are observed suggesting that the knee is dominated by heavy nuclei. It is planed to explicitly measure the iron spectrum at the knee by detecting the core structure of the air showers in the next phase of Tibet experiment ‘YAC’ ( Y angbajing A ir shower C ore array).

Primary cosmic-ray energy spectrum around the knee energy region measured by the Tibet hybrid experiment

The Tibet hybrid experiment composed of emulsion chambers, burst detectors and the Tibet-II air-shower array was done at Yangbajing (4,300m above sea level, 606 g/cm2) in Tibet. For 3-year operation, we have observed 177 γ-families (ΣEγ > 20 TeV) accompanying air showers with the shower size Ne > 2 × 105. Using this data set and a neural network method to select proton and helium induced events, the primary proton and helium energy spectra are deduced between 1015 and 1016eV. The proton spectrum, thus obtained, can be expressed by a single power-law function with a differential index of −3.01±0.11 and −3.05±0.12 based on the QGSJET+HD and SIBYLL+HD models, respectively, which are steeper than that extrapolated from the direct observations of −2.74±0.01 in the energy range below 1014 eV. Our result suggests that the main component responsible for making the knee structure of the all-particle spectrum should be composed of nuclei heavier than helium.

Automatic analysis of the emulsion chamber using the image scanner applied to the Tibet hybrid experiment

A new method of emulsion chamber analysis by using an image scanner is developed in order to process a large number of X-ray films used in the Tibet hybrid experiment which is carried out to observe primary proton spectrum in the PeV region. The experimental procedure in this method consists of multiple steps of processing the image of the X-ray films, i.e. the detection of cascade showers induced by high-energy γ-rays, the noise reduction, the spatial reconstruction of γ-families and their energy determination. The quality of the results obtained by the new method is compared with that obtained by the current manual method or Monte Carlo simulations. It is shown that sufficiently high quality of the emulsion chamber analysis is realized by the automatic method both in energy and geometrical resolution.