The all-particle spectrum of Galactic cosmic rays steepens abruptly at about 3 PeV (the 'knee'), and the prevailing paradigm attributes cosmic rays up to the knee to diffusive shock acceleration in supernova remnants (SNRs). Classical estimates cap ordinary SNR shocks near 0.1-0.3 PeV per unit charge, roughly an order of magnitude below the knee, unless magnetic fields at the shock are strongly amplified; moreover, rigidity-dependent cutoffs across a smooth source population generically produce a smooth rather than sharp knee. This problem asks whether standard SNR shock acceleration can accelerate hadrons to PeV energies, i.e. whether SNRs act as hadronic PeVatrons, and whether the SNR population with rigidity-dependent maximum energies can account quantitatively for the sharpness and composition of the knee. The literature from 2016 to 2026, spanning theory (magnetic-field amplification, time-dependent maximum energy, particle escape) and gamma-ray observations (H.E.S.S., HAWC, MAGIC, LHAASO; the Cassiopeia A limits and the G106.3+2.7 debate), finds the question open: no SNR has been confirmed as a hadronic PeVatron, deep limits on Cas A challenge the paradigm while competing analyses keep it alive, and the origin of the sharp knee is still contested between source and propagation explanations.

Background

Cosmic rays are relativistic charged particles (mostly protons and nuclei) that reach Earth almost isotropically from Galactic sources. Their all-particle flux follows a power law in energy that steepens at about 3 PeV (1 PeV = 10^15 eV), a feature called the 'knee', and flattens again near 3 EeV (the 'ankle'), beyond which an extragalactic component presumably dominates. The standard paradigm holds that Galactic cosmic rays are accelerated by diffusive shock acceleration (DSA) — a first-order Fermi process in which particles repeatedly cross a shock front and gain energy — in the shock waves of supernova remnants (SNRs), the expanding debris of stellar explosions. The paradigm is attractive because the kinetic energy released by Galactic supernovae matches the power needed to maintain the cosmic-ray population, but it must be tested at the highest Galactic energies: an accelerator that reaches PeV energies is called a 'PeVatron'.

The maximum energy EmaxE_{max} attainable by DSA in an SNR of age tt, shock speed vsv_s, and magnetic field BB is limited by the acceleration time and by the finite shock lifetime and size. For 'ordinary' remnant parameters, without invoking any mechanism beyond standard DSA in the compressed interstellar magnetic field, classical estimates give EmaxE_{max} of order Z×100Z \times 100 TeV for a nucleus of charge ZZ — about 0.3 PeV for protons, an order of magnitude below the knee. The route to higher energies is amplification of the magnetic field near the shock by the accelerated particles themselves, most prominently through the non-resonant streaming instability identified by Bell (2004), in which the electric current of escaping accelerated particles excites exponentially growing magnetic turbulence ahead of the shock; whether strongly amplified fields and the resulting PeV energies are a robust outcome of standard SNR evolution, or a fine-tuned possibility restricted to the first decades or centuries of a remnant's life, remains unsettled. Because EmaxE_{max} is expected to decline as the remnant ages and decelerates, most remnants observable today should no longer be PeVatrons — the 'low rate of Galactic pevatrons' — which sharply constrains how much of the PeV cosmic-ray flux the SNR population can supply.

A second, independent tension concerns the shape of the knee. If each nucleus species is cut off at a maximum energy set by its rigidity (momentum per unit electric charge), Emax≈Z×E0E_{max} \approx Z \times E_0, the all-particle spectrum is a superposition of smooth cutoffs and should bend gradually, whereas the measured knee is sharp: the spectral index changes abruptly near 3 PeV, and the LHAASO Collaboration's precise proton spectrum shows a hardening that culminates at 3 PeV followed by a distinct softening. Explanations that attribute the knee to a change in hadronic interaction physics are disfavored by collider data, leaving source-cutoff, source-population, and propagation explanations in competition. Whether an SNR population with rigidity-dependent cutoffs can reproduce the observed sharpness, together with the energy-dependent mass composition across the knee, without fine tuning, is an open quantitative question.

The problem is observationally testable because protons accelerated to PeV energies can be detected indirectly: collisions between accelerated hadrons and ambient gas produce neutral pions that decay into gamma rays, so gamma rays with energies of at least ~100 TeV trace protons of at least ~1 PeV. Current instruments (H.E.S.S., MAGIC, HAWC, Tibet ASgamma, and above all LHAASO with its km2-scale particle arrays) have reached the sensitivity required for this test. The observational record to date is tantalizing but inconclusive: the H.E.S.S. Collaboration has demonstrated that PeV protons are being accelerated in the diffuse source at the Galactic Centre (which is not an SNR); the twelve ultra-high-energy gamma-ray sources reported by LHAASO, including photons up to 1.4 PeV, are mostly associated with pulsar-wind nebulae — nebulae of energetic particles powered by the relativistic wind of a spinning neutron star (pulsar) — rather than SNR shocks; deep LHAASO limits on the young remnant Cassiopeia A challenge the prevailing paradigm that such remnants are major suppliers of PeV cosmic rays; and the promising candidate G106.3+2.7 is disputed between an SNR-shock origin and a pulsar-wind-nebula origin.

Problem Statement

Determine whether shock acceleration in supernova remnants can accelerate hadrons to petaelectronvolt energies and supply the cosmic-ray flux at the knee, specifically: (a) Theoretically, establish whether standard diffusive shock acceleration in SNRs — with realistic shock evolution, magnetic-field amplification driven by the accelerated particles themselves, and particle escape — reaches a proton maximum energy of at least ~3 PeV, or prove a quantitative upper bound below it; the calculation must state explicitly which physical ingredients (e.g. non-resonant field amplification, the very youngest, ejecta-dominated phases, in which the shock is still sweeping up the stellar ejecta rather than the ambient gas) are required and whether they are robust outcomes of standard SNR evolution. (b) Observationally, establish or exclude the existence of hadronic PeVatrons among supernova remnants, through gamma-ray (or neutrino) signatures of pion-decay emission extending to at least ~100 TeV from a source identified with an SNR shock rather than an associated pulsar-wind nebula or molecular cloud (a dense interstellar gas cloud), or through population-level limits that exclude SNRs as the dominant source of PeV cosmic rays. (c) Determine whether a Galactic population of SNR sources with rigidity-dependent cutoffs, together with standard propagation, can reproduce quantitatively the sharpness of the all-particle knee near 3 PeV, the LHAASO proton spectral hardening culminating at 3 PeV with subsequent sharp softening, and the energy-dependent composition across the knee, without fine tuning.

References

  1. Biplab Bijay and Arunava Bhadra, 'Progenitor model of cosmic ray knee', Research in Astronomy and Astrophysics 16, 006 (2016), DOI 10.1088/1674-4527/16/1/006, preprint arXiv:1412.0818, https://doi.org/10.1088/1674-4527/16/1/006 (preprint: https://arxiv.org/abs/1412.0818)
  2. A. R. Bell, 'Turbulent amplification of magnetic field and diffusive shock acceleration of cosmic rays', Monthly Notices of the Royal Astronomical Society 353, 550-558 (2004), DOI 10.1111/j.1365-2966.2004.08097.x, https://doi.org/10.1111/j.1365-2966.2004.08097.x
  3. K. M. Schure, A. R. Bell, L. O'C Drury and A. M. Bykov, 'Diffusive Shock Acceleration and Magnetic Field Amplification', Space Science Reviews 173, 491-519 (2012), DOI 10.1007/s11214-012-9871-7, https://doi.org/10.1007/s11214-012-9871-7
  4. P. Blasi, 'The origin of galactic cosmic rays', The Astronomy and Astrophysics Review 21, 70 (2013), DOI 10.1007/s00159-013-0070-7, https://doi.org/10.1007/s00159-013-0070-7
  5. S. Gabici, C. Evoli, D. Gaggero, P. Lipari, P. Mertsch, E. Orlando, A. Strong and A. Vittino, 'The origin of Galactic cosmic rays: challenges to the standard paradigm', International Journal of Modern Physics D 28, 1930022 (2019), DOI 10.1142/S0218271819300222, preprint arXiv:1903.11584, https://doi.org/10.1142/S0218271819300222 (preprint: https://arxiv.org/abs/1903.11584)
  6. P. Cristofari, S. Gabici, R. Terrier and T. B. Humensky, 'On the search for Galactic supernova remnant PeVatrons with current TeV instruments', Monthly Notices of the Royal Astronomical Society 479, 3415-3421 (2018), DOI 10.1093/mnras/sty1589, preprint arXiv:1803.09728, https://doi.org/10.1093/mnras/sty1589 (preprint: https://arxiv.org/abs/1803.09728)
  7. P. Cristofari, P. Blasi and E. Amato, 'The low rate of Galactic pevatrons', Astroparticle Physics 123, 102492 (2020), DOI 10.1016/j.astropartphys.2020.102492, preprint arXiv:2007.04294, https://doi.org/10.1016/j.astropartphys.2020.102492 (preprint: https://arxiv.org/abs/2007.04294)
  8. P. Cristofari, 'The Hunt for Pevatrons: The Case of Supernova Remnants', Universe 7, 324 (2021), DOI 10.3390/universe7090324, preprint arXiv:2110.07956, https://doi.org/10.3390/universe7090324 (preprint: https://arxiv.org/abs/2110.07956)
  9. H.E.S.S. Collaboration, 'Acceleration of petaelectronvolt protons in the Galactic Centre', Nature 531, 476-479 (2016), DOI 10.1038/nature17147, https://doi.org/10.1038/nature17147
  10. Z. Cao, F. Aharonian, et al. (LHAASO Collaboration), 'Ultrahigh-energy photons up to 1.4 petaelectronvolts from 12 gamma-ray Galactic sources', Nature 594, 33-36 (2021), DOI 10.1038/s41586-021-03498-z, https://doi.org/10.1038/s41586-021-03498-z
  11. A. Albert et al. (HAWC Collaboration), 'HAWC J2227+610 and its association with G106.3+2.7, a new potential Galactic PeVatron', The Astrophysical Journal Letters 896, L29 (2020), DOI 10.3847/2041-8213/ab96cc, preprint arXiv:2005.13699, https://doi.org/10.3847/2041-8213/ab96cc (preprint: https://arxiv.org/abs/2005.13699)
  12. H. Abe et al. (MAGIC Collaboration), 'MAGIC observations provide compelling evidence of the hadronic multi-TeV emission from the putative PeVatron SNR G106.3+2.7', Astronomy & Astrophysics 671, A12 (2023), DOI 10.1051/0004-6361/202244931, preprint arXiv:2211.15321, https://doi.org/10.1051/0004-6361/202244931 (preprint: https://arxiv.org/abs/2211.15321)
  13. R. Alfaro et al. (HAWC Collaboration), 'Testing the Molecular Cloud Paradigm for Ultra-High-Energy Gamma Ray Emission from the Direction of SNR G106.3+2.7', Astronomy & Astrophysics 691, A89 (2024), DOI 10.1051/0004-6361/202451514, preprint arXiv:2407.10729, https://doi.org/10.1051/0004-6361/202451514 (preprint: https://arxiv.org/abs/2407.10729)
  14. Z. Cao et al. (LHAASO Collaboration), 'Does or did the supernova remnant Cassiopeia A operate as a PeVatron?', arXiv:2310.17082 (2023), https://arxiv.org/abs/2310.17082
  15. Z. Cao et al. (LHAASO Collaboration), 'Broadband gamma-ray spectrum of supernova remnant Cassiopeia A', arXiv:2502.04848 (2025), https://arxiv.org/abs/2502.04848
  16. P. Blasi, 'Gamma radiation from cosmic rays escaping a young supernova remnant: The case of Cas A', arXiv:2509.08447 (2025), https://arxiv.org/abs/2509.08447
  17. Bo-Tao Li, Wei Wang and Zhuo Li, 'Particle Acceleration in Cassiopeia A Revealed by Broadband High-Energy Spectrum', arXiv:2603.26326 (2026), https://arxiv.org/abs/2603.26326
  18. P. Cristofari, G. Emery, T. Lubrano di Vavaria, H. Costantini, F. Cassol, M.-S. Carrasco and B. Le Nagat Neher, 'Can a single supernova remnant account for the gamma-ray emission of G106.3+2.7?', Astronomy & Astrophysics 710, A197 (2026), DOI 10.1051/0004-6361/202658874, preprint arXiv:2603.28701, https://doi.org/10.1051/0004-6361/202658874 (preprint: https://arxiv.org/abs/2603.28701)
  19. Z. Cao et al. (LHAASO Collaboration), 'Precise measurements of the cosmic ray proton energy spectrum in the "knee" region', Science Bulletin 70, 4173-4180 (2025), DOI 10.1016/j.scib.2025.10.048, preprint arXiv:2505.14447, https://doi.org/10.1016/j.scib.2025.10.048 (preprint: https://arxiv.org/abs/2505.14447)
  20. C. Prévotat, M. Kachelrieß, S. Koldobskiy, A. Neronov and D. Semikoz, 'Energy dependence of the knee in the cosmic ray spectrum across the Milky Way', Physical Review D 110, 103035 (2024), DOI 10.1103/PhysRevD.110.103035, preprint arXiv:2407.11911, https://doi.org/10.1103/PhysRevD.110.103035 (preprint: https://arxiv.org/abs/2407.11911)
  21. Hua Yue, Lin Nie, Yi-Qing Guo and Hong-Bo Hu, 'Joint constraint on the propagation origin of the cosmic-ray spectral knee from energy spectrum and anisotropy observations', The Astrophysical Journal 998, 264 (2026), DOI 10.3847/1538-4357/ae3d2d, preprint arXiv:2601.17851, https://doi.org/10.3847/1538-4357/ae3d2d (preprint: https://arxiv.org/abs/2601.17851)