The Scientific World Journal

Volume 2014 (2014), Article ID 278968, 6 pages

http://dx.doi.org/10.1155/2014/278968

Research Article

Simulating a Measurement of the 2nd Knee in the Cosmic Ray Spectrum with an Atmospheric Fluorescence Telescope Tower Array

Jiali Liu,¹ Qunyu Yang,² Yunxiang Bai,³ and Zhen Cao³

¹Physics Department, Kunming University, Kunming, Yunnan 650214, China

²Physics Department, Yunnan University, Kunming, Yunnan 650214, China

³Astro-Particle Physics Center, Institute of High Energy Physics, CAS, Beijing 100049, China

Received 28 December 2013; Accepted 21 January 2014; Published 9 March 2014

Academic Editors: P. Branchini, V. Flaminio, and A. Takamura

Copyright © 2014 Jiali Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A fluorescence telescope tower array has been designed to measure cosmic rays in the energy range of 10¹⁷–10¹⁸ eV. A full Monte Carlo simulation, including air shower production, light generation and propagation, detector response, electronics, and trigger system, has been developed for that purpose. Using such a simulation tool, the detector configuration, which includes one main tower array and two side-trigger arrays, 24 telescopes in total, has been optimized. The aperture and the event rate have been estimated. Furthermore, the performance of the technique in measuring composition has also been studied.

1. Introduction

The cosmic ray (CR) all-particle spectrum roughly follows a power law, that is, . The index is not the same in different energy regions. It reveals different features as shown in Figure 1, where the flux is multiplied by [1]. One important feature is the knee at  eV, which the spectrum index changes from in the low energy region to . The second important feature is the second knee around  eV, where the spectrum index suddenly changes to . Another important feature is the ankle at about  eV, where the spectrum seems to become flat with again. The other important feature is Greisen-Zatsepin-Kuzmin suppression (called the GZK cutoff) [2] above about  eV, which was observed for the first time by the HiRes experiment [3].

Figure 1: The all-particle cosmic ray energy spectrum multiplied by . The dots represent the results observed by all experiments. The curves show the results calculated following the polygonato model [15]. The contribution of each single element (marked with different number ) and the all-particle spectrum are shown separately.

Measuring these structures precisely will improve our understanding of the origin of CRs. However, because the dynamic range of individual experiments is limited and has its own energy scale and uncertainty, the observed results are not consistent among them. Furthermore, in composition observations, different experiments have different systematic uncertainties since they use different techniques, leading to the inconsistencies among the results as summarized in Figure 2. Therefore, to accurately measure the structures in the FF energy spectrum and the composition of CRs, it is ideal to have experiments being carried out at the same place with sufficient overlap among them. This would allow us to obtain a complete and self-consistent observation over the whole energy range from 100 TeV to several EeV.

Figure 2: Mean distribution as a function of energy [16]. The dots are experimental results. The lines are simulation results for proton and iron predicted by Corsika [17] with the hadronic interaction model QGSJET 01 (solid line), QGSJET II-03 (dashed line), SIBYLL (dot line), and EPOS 1.6 (dot-dashed line) [18, 19].

The project of the Wide Field of View (FOV) Cherenkov/Fluorescence Telescope Array (WFCTA), composed of 24 telescopes, has been designed to reach these goals. One of the physical aims is to observe the energy spectrum and composition of CRs from  eV to  eV, covering both knees [4, 5]. To cover such a wide dynamic range, three observational stages with the same telescopes will be used. The first two stages will use the Cherenkov technique under two different configurations, covering the energy range of  eV– eV and  eV– eV. The first stage will overlap at lower energies with space/balloon borne measurements, which will provide a bridge between space/balloon borne measurements and ground based (low altitude) measurements. The third stage will focus on the energy range of  eV– eV, using the fluorescence technique by reconfigurating the telescopes. The overlap between these three stages will transfer the accurate energy scale of space/balloon borne experiments to the ground observatory and provide cross-calibration in energy. The focus of this paper is to estimate the performance of the fluorescence telescope tower array (the third stage of WFCTA) for the measurement of spectrum and composition around the second knee. Section 2 describes the detector features and the telescope array configuration. Section 3 will give details of the simulation chain. The simulation results on aperture, event rate, and on the performance of composition measurement are presented in Section 4. Section 5 provides summary of results.

2. Detectors

The proposed fluorescence detector array is composed of 24 telescopes, located at the three vertices of a triangle. As suggested by the configuration optimization, the chosen base line is 8 km and the height is 5 km. The three positions are denoted as FD1, FD2, and FD3, respectively. FD1 is configured as a array, covering a range of in azimuth and in elevation starting from . FD2 and FD3 are located at the bottom vertices, which are configured as a array, covering a range of in azimuth and in elevation starting from . Such a configuration is shown in Figure 3.

Figure 3: The configuration of the fluorescence telescope array. The solid circles are the fluorescence telescopes. The lines illustrate the detector FOV along the azimuth.

A light collecting mirror, composed of a 20 hexagonal mirror segments, with a reflectivity of 82% is used for each telescope. The image camera is made of pixels. Each pixel is a 40 mm hexagonal photomultiplier tube (PMT) that has about a FOV. Each PMT is read out by a 50 MHz flash ADC to provide a measurement of the shower waveform signals. A peak finding algorithm has been developed in order to provide an individual channel trigger using a field programmable gate array (FPGA). Three trigger levels are required for each event. The first level trigger is set by requiring the signal-noise ratio to be greater than for the FD1 array and for FD2 and FD3, where is the standard deviation of the total photoelectron noise within a running window of 320 ns. The second level trigger uses adjacent triggered tubes within a box running over a telescope camera to form pattern recognition. The second level trigger is referred to as the telescope trigger. The “track-type” pattern requires at least six triggered pixels forming a straight line. The third level trigger is the event trigger which requires that both the tower array and one side array be triggered. The detailed description of the telescope can be found in [6].

3. Monte Carlo Simulation

In the simulation, cosmic rays are sampled isotropically and uniformly. The impact parameter, (i.e., the distance from FD1 to the shower axis) is limited and is less than 10 km. An spectrum is assumed for the resolution study. 728,555 proton and 854,342 iron primary showers have been generated in the energy range from  eV to  eV. To maintain statistics at high energies, 386,080 proton showers following an spectrum covering the same energy range used in the resolution study are generated for the detector aperture estimation.

3.1. Air Shower Simulation

Given the shower energy and geometry (i.e., zenith angle, azimuth angle, , and core location), the shower longitudinal development is parameterized by a function with three parameters [7]: the maximum position of shower development , the maximum number of charged particles , and the width of the shower , including their energy dependence, fluctuations, and the correlations between them. The number of charged particles at atmospheric depth is calculated as

The shower lateral spread is taken into account in the simulation of photon production.

3.2. Photon Production and Light Propagation

The fluorescence and Cherenkov photons production are considered in detail as explained below.

Ultraviolet fluorescence light is generated as charged shower particles pass through the atmosphere. The fluorescence yield measured in [8] is used. Laterally, fluorescence photons are spread out using the Nishimura-Kamata-Greisen (NKG) function [9]

where is the Moliere radius and is the age of the shower. The normalized function reads as

Cherenkov photons are also generated by charged particles if the particle energy is higher than the threshold energy [10]. Scattered by the atmospheric molecules (Rayleigh scattering) and aerosols (Mie scattering), photons are distributed in all directions according to corresponding phase functions. A standard desert aerosol model [11] with a scale height of 1 km and a horizontal attenuation length of 25 km is assumed in the simulation. The same model is used for estimating the photon flux attenuation due to scattering. A ray-tracing procedure is carried out to trace each photon to the photocathodes of PMTs once the photon arrives to the entrance of a telescope. All detector responses are considered in the ray-tracing procedure. In Figure 4, a typical example of profiles for all kinds of photons that are produced by a shower in the simulation is plotted as a function of slant atmospheric depth along with the shower longitudinal development, the range between two lines being covered by FD1. The detailed description of photon production and propagation can be found in [12, 13] and references therein.

Figure 4: Profiles of all photons along with the shower longitudinal development that were produced by a shower in simulation. The solid curve represents the sum of all photons. The short-dot-dashed curve represents fluorescence photons. The long-dot-dashed curve represents direct Cherenkov components. The dashed curve represents Cherenkov photons scattered by aerosols. The lowest dotted curve represents Cherenkov photons scattered by atmospheric molecules (Rayleigh scattering). The range between two vertical lines is covered by FD1.

3.3. Electronics and Noise

All photons collected by one PMT are distributed in flash ADC bins according to their arrival time. Night sky background (NSB) photons with an average flux of [14] are randomly added to the waveform. The electronic noise with a mean of 1.2 FADC counts is also added to every 20 ns time window with fluctuations. The Trigger algorithm of the three levels as described in Section 2 is used in the simulation. An example of a detected CR event is shown in Figure 5.

Figure 5: A typical air shower seen by the WFCTA fluorescence detectors. Each square marked with a number shows the FOV of each telescope. The open circles in the square represent triggered tubes and the size of each circle is proportional to logarithm of the number of photons. The solid curve crossing the circles represents the projection of a plane containing the shower axis and the detector. Energy (in PeV) and zenith and azimuth angle (in degrees) are displayed at the top.

4. Simulation Results

Using the simulation procedure and the optimized detector configuration described above, the detector aperture and event rate have been estimated.

4.1. Detector Aperture and Event Rate

The detector triggered aperture has been estimated. This is displayed in Figure 6 with solid circles.

Figure 6: Detector aperture in different conditions. Solid circles are the aperture without any cut. Solid triangles are obtained using cut1 described in the text. Asterisks show the results obtained using cut2.

It can be seen in Figure 5 that the tower detector has two vertical edges that could cause an incomplete measurement of images, which may cause errors in estimates of energy and direction. To maintain event reconstruction quality, showers that touch the detector boundaries, with an angle between image center (weighted by signals) and mirror edge smaller than , should be cut. To avoid Cherenkov light contamination, all tubes that have viewing angles smaller than are removed and track lengths smaller than are cut. Furthermore, at the WFCTA site located at a height of 4.3 km above sea-level, some showers with large will hit the ground before reaching the . On the other hand, for some showers having small , the location of maximum shower development will fall outside the upper edge of the tower array FOV. To guarantee the required resolution in the measurement of , the showers’ location must fall inside the detector FOV. There are about 18% of showers with falling outside the detector FOV. We denote these cut conditions by cut1, after applying cut1, the detector aperture is reduced as shown in Figure 6 by solid triangles.

In order to measure the structure of the cosmic ray energy spectrum around with minimum bias, a flat aperture as a function of energy is required. The triggered event density distributions as a function of in different energy regions are shown in Figure 7. According to such an event density distribution, following geometric constraint based on is applied. The events in the intervals 16.6~17.0, 17.0~17.8, and 17.8~18.3 have smaller than 5.5 km, 7.0 km, and 8.0 km, respectively. Cut1 together with geometric constraint is denoted cut2. The application of cut2 will yield a flat aperture as shown in Figure 6 by asterisks. The corresponding event rates are displayed in Figure 8. After all cuts, 14.3 k events with and 162 events with are obtained per year.

Figure 7: Event density distribution as a function of . Different points represent different energies.

Figure 8: Event rate per year in three conditions, respectively.

4.2. Distribution

The distributions as a function of energy are displayed in Figure 9. The two point-curves are for primary proton and iron nuclei. The solid points are triggered events without any cut. The open points are for events after cut1 as described in Section 4.1. It indicates that the maximum detected bias is smaller than . Estimations of the resolution are ongoing.

Figure 9: distribution as a function of energy. Two lines are for primary proton and iron. Solid points are triggered events without any cut. Opened points are for events after cut1 described in Section 4.1.

5. Summary

The WFCTA fluorescence detector array simulation procedure has been completed. Based on the simulation tools, the detector array configuration optimization has been performed. To maintain data reconstruction quality, some cuts have been applied. In order to obtain the flat aperture needed for a precise study of the energy spectrum, further cuts have been applied. Detector apertures and event rates under three conditions have been evaluated. 14.3 k events with and 162 events with are left after all cuts. The distribution as a function of energy has been obtained. There is only a small detection bias both before and after event cuts. Therefore, it will be possible to carry out a detailed study of the CR energy spectrum and composition in the Sub-EeV range using such a fluorescence telescope array.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors appreciate the reviewers for their careful reading of a preliminary version of the paper and their helpful comments and valuable suggestions, which have led to a great improvement in this paper. This work is funded by Kunming University through the Grant of YJL11008.

References

J. R. Horandel, “Cosmic rays from the knee to the second knee: 10¹⁷ to 10¹⁸eV,” Modern Physics Letters A, vol. 22, pp. 1533–1551, 2007.

G. T. Zatsepin and V. A. Kuzmin, “Upper Limit of the Spectrum of Cosmic Rays,” JETP Letters, vol. 4, pp. 78–80, 1966.

D. R. Bergman, “First observation of the Greisen-Zatsepin-Kuzmin suppression,” Physical Review Letters, vol. 100, no. 10, Article ID 101101, 5 pages, 2008.

C. Zhen, “A future project at tibet: the large high altitude air shower observatory (LHAASO),” Chinese Physics C, vol. 34, no. 2, pp. 249–252, 2010.

M. Zha, “Status of the large high altitude air shower observatory project,” Nuclear Instruments and Methods in Physics Research A, vol. 692, pp. 77–82, 2012.

S. S. Zhang, Y. X. Bai, Z. Cao et al., “Properties and performance of two wide field of view Cherenkov/fluorescence telescope array prototypes,” Nuclear Instruments and Methods in Physics Research A, vol. 629, no. 1, pp. 57–65, 2011.

T. Abu-Zayyad, K. Belov, D. J. Bird et al., “A measurement of the average longitudinal development profile of cosmic ray air showers between 1017 and 1018 eV,” Astroparticle Physics, vol. 16, no. 1, pp. 1–11, 2001.

F. Kakimoto, E. C. Loh, M. Nagano, H. Okuno, M. Teshima, and S. Ueno, “A measurement of the air fluorescence yield,” Nuclear Instruments and Methods in Physics Research A, vol. 372, no. 3, pp. 527–533, 1996.

K. Kamata and J. Nishimura, “The lateral and the angular structure functions of electron showers,” Progress of Theoretical Physics Supplements, vol. 6, pp. 93–155, 1958.

R. M. Baltrusaitis, R. Cady, G. L. Cassiday et al., “The Utah Fly's Eye detector,” Nuclear Instruments and Methods in Physics Research A, vol. 240, no. 2, pp. 410–428, 1985.

D. R. Longtin, “A wind dependent desert aerosol model: radiative properties,” Tech. Rep. AFL-TR-88-0112, Air Force Geophysics Laboratory, 1998.

Z. Cao, M. A. Huang, P. Sokolsky, and Y. Hu, “Ultra high energy ν τ detection with a cosmic ray tau neutrino telescope using fluorescence/Cerenkov light technique,” Journal of Physics G, vol. 31, no. 7, pp. 571–582, 2005.

J. L. Liu, S. S. Zhang, Z. Cao, et al., “Simulation of the cosmic ray tau neutrino telescope (CRTNT) experiment,” Journal of Physics G, vol. 36, no. 7, Article ID 075201, 2009.

D. J. Bird, S. C. Corbató, H. Y. Dai, et al., “Evidence for correlated changes in the spectrum and composition of cosmic rays at extremely high energies,” Physics Review Letter, vol. 71, no. 21, pp. 3401–3404, 1993.

J. R. Hörandel, “On the knee in the energy spectrum of cosmic rays,” Astroparticle Physics, vol. 19, no. 2, pp. 193–220, 2003.

J. R. Horandel, “A review of experimental results at the knee,” Journal of Physics: Conference Series, vol. 47, no. 41, 2006.

D. Heck, J. Knapp, J. N. Capdevielle, G. Schatz, and T. Thouw, Corsika: A Monte Carlo Code to Simulate Extensive Air Showers, vol. 6019 of Forschungszentrum Karlsruhe Report, FZKA, 1998.

J. R. Horandel, “On total inelastic cross sections and the average depth of the maximum of extensive air showers,” Journal of Physics G, vol. 29, no. 11, p. 2439, 2002.

J. Blumer, R. Engel, and J. R. Horandel, “Cosmic rays from the knee to the highest energies,” Progress in Particle and Nuclear Physics, vol. 63, no. 2, pp. 293–338, 2009.