The all-particle cosmic-ray spectrum shows a 'second knee' near 10^17 eV, between the first knee (~3x10^15 eV) and the ankle (~5x10^18 eV). Theory expects this interval to contain the transition from Galactic to extragalactic dominance, which should imprint a characteristic evolution on the primary mass composition: increasingly heavy nuclei up to the iron cutoff of Galactic accelerators, followed by a return to lighter, proton-rich primaries of extragalactic origin. Establishing the actual mass mixture and its energy evolution across 10^16-10^18 eV, at a level that is robust against hadronic-interaction-model and detector systematics and consistent across experiments, remains an open measurement problem. It is the decisive observable for where and how the Galactic-to-extragalactic transition occurs.

Background

Cosmic rays are atomic nuclei (and electrons) accelerated by astrophysical sources and arriving at Earth as an all-particle flux. Above ~10^14 eV the flux is too low for direct satellite or balloon detection, so the mass of the primary particle must be inferred indirectly from extensive air showers (EAS), the cascades of secondary particles the primary initiates in the atmosphere. The most composition-sensitive air-shower observable is XmaxX_\mathrm{max}, the atmospheric depth (in g/cm2^2) at which the shower reaches its maximum particle number: heavier nuclei interact higher up and produce on average shallower showers, so the mean ⟨Xmax⟩\langle X_\mathrm{max}\rangle, its fluctuation σ(Xmax)\sigma(X_\mathrm{max}), and their variation with energy (the 'elongation rate') constrain the primary mass mixture. Translating them into elemental fractions or the mean logarithmic mass ⟨ln⁡A⟩\langle\ln A\rangle requires comparing data with Monte Carlo simulations built on post-LHC hadronic interaction models such as QGSJetII-04, EPOS-LHC, and Sibyll 2.3d, which differ in their predicted ⟨Xmax⟩\langle X_\mathrm{max}\rangle by 10-20 g/cm2^2 - comparable to one mass step; this model dependence is a central systematic in all composition measurements. Non-imaging Cherenkov arrays (KASCADE-Grande, Tunka-133, Yakutsk), muon/electron surface arrays (IceTop, LHAASO-KM2A), and radio arrays (LOFAR, the LOw-Frequency ARray radio telescope) provide complementary, similarly model-dependent observables. The all-particle spectrum steepens at the first knee (~3x10^15 eV), commonly interpreted as the rigidity-dependent (R=pc/(eZ)R = pc/(eZ), scaling as the nuclear charge ZZ) acceleration or confinement limit of Galactic sources ('Peters cycle'), so that heavier nuclei cut off at ZZ-times-higher energies, culminating in an iron cutoff near 10^17 eV where a further steepening - the second knee - is observed by several experiments. At the highest energies, measurements by the Pierre Auger Observatory and Telescope Array (TA) show the composition becoming lighter again near the ankle, and Auger's large-scale anisotropy measurement above 8x10^18 eV indicates extragalactic origin there. The interval 10^17-10^18 eV is therefore expected to hold the transition from Galactic to extragalactic dominance, and the composition is the decisive observable for locating and characterizing it: a rigidity-cutoff picture predicts heavy (iron-group-dominated) composition up to ~10^17 eV, while alternative scenarios - a second Galactic component (e.g. Wolf-Rayet supernovae or reacceleration at the Galactic wind termination shock), or an early appearance of extragalactic primaries - predict different mixtures. The experimental situation in this interval has historically been sparse: fluorescence measurements - imaging the ultraviolet light emitted along the shower track, possible only on clear, moonless nights - have excellent energy resolution but a ~10% duty cycle and limited statistics below 10^18 eV, while ground arrays offer exposure but stronger model dependence. LOFAR radio measurements found a surprisingly large light (p+He) fraction of ~0.8 (99% CL interval [0.38, 0.98]) at 10^17-10^17.5 eV, in tension with a naive iron-dominated Galactic picture at that energy; the NICHE (Non-Imaging CHErenkov Array)/TALE program was motivated specifically by the lack of composition-sensitive measurements in the 10^16-10^18 eV transition.

Problem Statement

Determine the mass composition of primary cosmic rays across the energy interval 101610^{16}-1018 eV10^{18}\ \mathrm{eV}, spanning the second knee (∼1017 eV\sim 10^{17}\ \mathrm{eV}) and the expected transition from Galactic to extragalactic sources: establish the primary mass mixture (elemental or mass-group fractions, or an equivalent characterization such as the energy-dependent mean logarithmic mass ⟨ln⁡A⟩\langle\ln A\rangle with its distribution width) and its evolution with energy over this interval, with uncertainties controlled well enough - across hadronic-interaction-model choices and different measurement techniques - to establish (a) whether the composition indeed becomes heavier up to the second knee and lighter above it, as the rigidity-dependent Galactic-cutoff picture predicts, and (b) at what energy and how sharply the dominant contribution shifts from Galactic to extragalactic sources, or alternatively whether a distinct Galactic component (or an early extragalactic component) is required to explain the observed mixture.

References

  1. Telescope Array Collaboration: R. U. Abbasi et al., 'Cosmic ray mass composition measurement in the energy range from 10^16.5 eV to 10^18.5 eV observed with the TALE hybrid detector', Phys. Rev. D 113, 062003 (2026). DOI: 10.48550/arXiv.2603.14804; https://arxiv.org/abs/2603.14804
  2. K. Fujita (Telescope Array Collaboration), 'Cosmic ray mass composition measurement with the TALE hybrid detector', EPJ Web of Conferences 283, 02009 (2023). DOI: 10.1051/epjconf/202328302009; https://doi.org/10.1051/epjconf/202328302009
  3. Y. Omura, D. R. Bergman, J. F. Krizmanic, K. Nakai, Y. Tsunesada, 'NICHE detector and operations', Journal of Physics: Conference Series 1468, 012097 (2020). DOI: 10.1088/1742-6596/1468/1/012097; https://doi.org/10.1088/1742-6596/1468/1/012097
  4. S. Buitink et al., 'A large light-mass component of cosmic rays at 10^17-10^17.5 electronvolts from radio observations', Nature 531, 70-73 (2016). DOI: 10.1038/nature16976; https://doi.org/10.1038/nature16976 (arXiv:1603.01594)
  5. M. G. Aartsen et al. (IceCube Collaboration), 'Cosmic ray spectrum and composition from PeV to EeV using 3 years of data from IceTop and IceCube', Physical Review D 100, 082002 (2019). DOI: 10.1103/PhysRevD.100.082002; https://doi.org/10.1103/PhysRevD.100.082002 (arXiv:1906.04317)
  6. Z. Cao et al. (LHAASO Collaboration), 'Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A', Physical Review Letters 132, 131002 (2024). DOI: 10.1103/PhysRevLett.132.131002; https://doi.org/10.1103/PhysRevLett.132.131002 (arXiv:2403.10010)
  7. R. U. Abbasi et al. (Telescope Array Collaboration), 'Depth of Ultra High Energy Cosmic Ray Induced Air Shower Maxima Measured by the Telescope Array Black Rock and Long Ridge FADC Fluorescence Detectors and Surface Array in Hybrid Mode', The Astrophysical Journal 858, 76 (2018). DOI: 10.3847/1538-4357/aabad7; https://doi.org/10.3847/1538-4357/aabad7 (arXiv:1801.09784)
  8. R. U. Abbasi et al. (Telescope Array Collaboration), 'The Cosmic-Ray Composition between 2 PeV and 2 EeV Observed with the TALE Detector in Monocular Mode', The Astrophysical Journal 909, 178 (2021). DOI: 10.3847/1538-4357/abdd30; https://doi.org/10.3847/1538-4357/abdd30 (arXiv:2012.10372)
  9. The Pierre Auger Collaboration: A. Abdul Halim et al., 'Constraining the sources of ultra-high-energy cosmic rays across and above the ankle with the spectrum and composition data measured at the Pierre Auger Observatory', Journal of Cosmology and Astroparticle Physics 2023(05), 024 (2023). DOI: 10.1088/1475-7516/2023/05/024; https://doi.org/10.1088/1475-7516/2023/05/024 (arXiv:2211.02857)
  10. The Pierre Auger Collaboration: P. Abreu et al., 'The energy spectrum of cosmic rays beyond the turn-down around 10^17 eV as measured with the surface detector of the Pierre Auger Observatory', The European Physical Journal C 81, 966 (2021). DOI: 10.1140/epjc/s10052-021-09700-w; https://doi.org/10.1140/epjc/s10052-021-09700-w (arXiv:2109.13400)
  11. S. Thoudam, J. P. Rachen, A. van Vliet, A. Achterberg, S. Buitink, H. Falcke, J. R. Hoerandel, 'Cosmic-ray energy spectrum and composition up to the ankle: the case for a second Galactic component', Astronomy & Astrophysics 595, A33 (2016). DOI: 10.1051/0004-6361/201628894; https://doi.org/10.1051/0004-6361/201628894 (arXiv:1605.03111)
  12. F. G. Schroeder, 'News from Cosmic Ray Air Showers (ICRC 2019 - Cosmic Ray Indirect Rapporteur)', published as PoS(ICRC2019)030 (2019). DOI: 10.48550/arXiv.1910.03721; https://arxiv.org/abs/1910.03721
  13. R. Alves Batista, J. Biteau, M. Bustamante, K. Dolag, R. Engel, K. Fang, K.-H. Kampert, D. Kostunin, M. Mostafa, K. Murase, F. Oikonomou, A. V. Olinto, M. I. Panasyuk, G. Sigl, A. M. Taylor, M. Unger, 'Open Questions in Cosmic-Ray Research at Ultrahigh Energies', Frontiers in Astronomy and Space Sciences 6, 23 (2019). DOI: 10.3389/fspas.2019.00023; https://doi.org/10.3389/fspas.2019.00023
  14. J. Albrecht et al., 'The Muon Puzzle in cosmic-ray induced air showers and its connection to the Large Hadron Collider', Astrophysics and Space Science 367, 27 (2022). DOI: 10.1007/s10509-022-04054-5; https://doi.org/10.1007/s10509-022-04054-5 (arXiv:2105.06148)
  15. B.-Q. Qiao, Q. Yuan, Y.-Q. Guo, 'Co-evolution of cosmic ray energy spectra, composition, and anisotropies', Journal of Cosmology and Astroparticle Physics 2026(07), 050 (2026). DOI: 10.1088/1475-7516/2026/07/050; https://doi.org/10.1088/1475-7516/2026/07/050 (arXiv:2506.18118)
  16. K. Fujita (Telescope Array Collaboration), 'Mass Composition Measurement with the TALE Hybrid Detector', Proceedings of the 39th International Cosmic Ray Conference (PoS ICRC2025), 0261 (2025). DOI: 10.22323/1.501.0261; https://doi.org/10.22323/1.501.0261
  17. S. Buitink et al., 'Anomalous Air Showers and What They Reveal About Hadronic Interactions and Cosmic-ray Masses', contribution to 'Advancing Astrophysics with the SKA II' (AASKAII), 2026. DOI: 10.48550/arXiv.2607.00135; https://arxiv.org/abs/2607.00135
  18. E. Roulet, 'Combined fit above 0.1 EeV to the cosmic-ray spectrum and composition measured with the Pierre Auger Observatory', Proceedings of the 39th International Cosmic Ray Conference (PoS ICRC2025), 373 (2025). DOI: 10.48550/arXiv.2507.10806; https://arxiv.org/abs/2507.10806