Fermi (1949) proposed that Galactic cosmic rays are accelerated stochastically by collisions against moving interstellar magnetic irregularities, and explicitly flagged what he could not solve: the injection of the heavy nuclei observed in the primary radiation (his theory's 'chief difficulty', with pre-acceleration thresholds rising from about 200 MeV for protons to about 300 GeV for iron), the microphysics behind the roughly one-light-year effective scattering mean free path, the value of the proton regeneration factor k, and the rate and energy limits of the mechanism. Most of this bundle has since been resolved by a change of framework: interstellar second-order acceleration was superseded as the primary source by first-order diffusive shock acceleration at supernova-remnant shocks (Bell 1978; Blandford and Ostriker 1978); the heavy-nuclei source composition is quantitatively explained by shock acceleration of gas-phase ions together with dust-grain acceleration and sputtering (Meyer, Drury and Ellison 1997), with fresh massive-star seeds (Binns et al. 2016) and direct evidence of proton acceleration in remnants (Ackermann et al. 2013); the energy budget is supplied by supernova power, dissolving the regeneration question. What survives is a precise joint core: identify the kinetic mechanism and efficiency with which thermal ions, including heavy nuclei, are injected into diffusive shock acceleration - fixing the mass-to-charge- and volatility-dependent source composition - and determine which Galactic source class (supernova remnants, massive-star superbubbles, or black-hole microquasars) accelerates nuclei to the knee near 3 PeV. In the 2025-2026 LHAASO measurements this core remains open: microquasars and a Cygnus superbubble are established PeV accelerators and the proton spectrum shows a new hard PeV component, but its origin and the injection microphysics remain unsettled.

Background

Cosmic rays and the knee. Primary cosmic rays are the relativistic nuclei (mostly hydrogen and helium, with about one percent heavier nuclei such as carbon, oxygen and iron, plus a small electron component) arriving at Earth from astrophysical sources, as distinct from the secondary particles produced when primaries strike the atmosphere. Their intensity follows an approximate power law in energy over ten decades and steepens at the knee near 3×10153\times10^{15} eV. The knee is generally interpreted as a rigidity-dependent effect (rigidity is momentum per unit charge, R=pc/(Ze)R=pc/(Ze)): each nuclear species cuts off at roughly the same rigidity, so heavier species dominate progressively above the knee. Fermi's mechanism and its open issues. Fermi (1949) proposed that cosmic rays are accelerated in interstellar space itself: charged particles spiral along the Galactic magnetic field and are reflected by moving, conducting irregularities of that field, gaining energy stochastically at each encounter. Since the fractional energy change per collision is second order in the irregularities' velocity, ΔE/E∼β2\Delta E/E\sim\beta^2 with β=V/c≈10−4\beta=V/c\approx10^{-4} for streaming speeds V∼30V\sim30 km/s, this is called second-order Fermi acceleration. With a mean interval τ≈1.3\tau\approx1.3 years between effective scattering collisions (an effective mean free path of order one light year) and a residence time T≈6×107T\approx6\times10^{7} years set by nuclear collisions with interstellar gas, the mechanism produces an inverse power-law spectrum with exponent 1+τ/(β2T)1+\tau/(\beta^2 T), which Fermi normalized to the then-observed value 2.9. Because ionization losses dominate below a threshold, the process cannot start from thermal energies: Fermi estimated injection requires about 200 MeV for protons, about 1 GeV for helium nuclei, and - because ionization loss grows with charge squared - roughly 20 GeV for oxygen and 300 GeV for an iron nucleus. He acknowledged he could not explain how the observed heavy nuclei are injected ('I do not know a plausible answer'), suggesting they are 'injected by a totally different mechanism, perhaps as a consequence of the stellar magnetism'; his abstract calls this the theory's chief difficulty. He likewise left at order-of-magnitude level the statistical properties of the magnetic irregularities, the value of the proton reproduction factor kk (the mean number of new above-threshold protons produced per interstellar hadronic collision, argued to sit near unity by a self-stabilization argument rather than calculated), and the rate, upper energy, and energy budget of the mechanism. The shock-acceleration replacement. In 1977-1978 several groups independently proposed what is now called diffusive shock acceleration (DSA, or first-order Fermi acceleration), with canonical treatments by Bell (1978) and Blandford and Ostriker (1978): particles that scatter back and forth across a collisionless shock - for example the blast wave of a supernova remnant (SNR), the diffuse shock produced when a stellar explosion sweeps up interstellar gas - gain energy in a systematic, first-order way at each crossing, and the repeated crossings generically yield a power-law spectrum with an exponent close to the observed one. Because the gain is first order in the shock speed, DSA is orders of magnitude faster than Fermi's interstellar process, whose own parameters imply an energy e-folding time of order 10810^{8} years. The reviews of Drury (1983), Blandford and Eichler (1987), and Blasi (2013) document the resulting paradigm - SNR shocks, powered by supernova kinetic energy, accelerate the bulk of Galactic nuclei up to rigidity 10610^{6}-10710^{7} GV - while recording its weak points. The modern forms of Fermi's open issues. Within DSA, injection reappears as a kinetic question: what fraction of the thermal ions at a shock enters the acceleration process, with what dependence on mass-to-charge ratio A/ZA/Z? The observed source composition - recovered by subtracting the secondaries produced by spallation (fragmentation in interstellar collisions) - shows refractory elements (those depleted into interstellar dust) enhanced relative to volatile ones; Meyer, Drury and Ellison (1997) showed this pattern is reproduced if shocks accelerate gas-phase ions together with dust grains whose sputtering (erosion by impacting particles) feeds refractory ions back into the accelerator. Isotopic measurements with the Advanced Composition Explorer (ACE) satellite's cosmic-ray isotope spectrometer show the accelerated material contains the radioactive clock isotope 60^{60}Fe, indicating it was synthesized within a few million years of acceleration (Binns et al. 2016). The maximum energy is governed by the scattering turbulence itself, which the accelerated particles amplify (Bell 2004), so Fermi's mean-free-path question survives inside the energy-ceiling problem. Hadronic acceleration is diagnosed through gamma rays from the decay of neutral pions produced when accelerated nuclei strike ambient gas - the pion-decay bump detected in SNR spectra by the Fermi satellite's Large Area Telescope (LAT; Ackermann et al. 2013) - and through ultra-high-energy (above 100 TeV) gamma-ray surveys by the High Energy Stereoscopic System (H.E.S.S.) and the Large High Altitude Air Shower Observatory (LHAASO), which have revealed candidate PeVatrons (sources accelerating particles past 101510^{15} eV), including a superbubble (a low-density cavity in interstellar gas blown by the combined winds and supernovae of a massive-star cluster) and microquasars (accreting stellar-mass black holes with relativistic jets), while classical SNRs have not yet been shown to reach PeV energies.

Problem Statement

The question is the one Fermi (1949) left open, reduced to its surviving core as established by the later literature (the resolved components are documented in Background). Stated in its modern, source-faithful form, the problem is to determine, self-consistently and quantitatively, how Galactic cosmic-ray nuclei are injected into, and accelerated to the highest observed energies by, their sources. Concretely: (1) identify the kinetic mechanism by which thermal background ions - including the heavy nuclei observed in the primary cosmic radiation - enter diffusive shock acceleration at collisionless shocks, and predict the injection efficiency and its dependence on mass-to-charge ratio and volatility well enough to reproduce the measured cosmic-ray source composition, including the refractory-versus-volatile enhancement pattern; and (2) determine the maximum energy reachable by each candidate Galactic source class - classical supernova remnants, massive-star-cluster superbubbles, and black-hole microquasars - and thereby establish which class, or combination of classes, supplies the nuclei observed at and above the knee of the all-particle spectrum near 3×10153\times10^{15} eV. The two sub-questions are physically coupled: the same coupled physics (injection, scattering by self-generated turbulence, acceleration, and escape) fixes both the composition and the maximum energy, and both are constrained by the same body of composition, gamma-ray, and charged-particle data.

References

  1. E. Fermi, 'On the Origin of the Cosmic Radiation', Physical Review 75 (8), 1169-1174 (1949), DOI: 10.1103/PhysRev.75.1169, https://doi.org/10.1103/PhysRev.75.1169
  2. E. Fermi, 'Galactic Magnetic Fields and the Origin of Cosmic Radiation', The Astrophysical Journal 119, 1 (1954), DOI: 10.1086/145789, https://doi.org/10.1086/145789
  3. A. R. Bell, 'The acceleration of cosmic rays in shock fronts - I', Monthly Notices of the Royal Astronomical Society 182 (2), 147-156 (1978), DOI: 10.1093/mnras/182.2.147, https://doi.org/10.1093/mnras/182.2.147
  4. R. D. Blandford and J. P. Ostriker, 'Particle acceleration by astrophysical shocks', The Astrophysical Journal 221, L29 (1978), DOI: 10.1086/182658, https://doi.org/10.1086/182658
  5. L. O'C. Drury, 'An introduction to the theory of diffusive shock acceleration of energetic particles in tenuous plasmas', Reports on Progress in Physics 46 (8), 973-1027 (1983), DOI: 10.1088/0034-4885/46/8/002, https://doi.org/10.1088/0034-4885/46/8/002
  6. R. Blandford and D. Eichler, 'Particle acceleration at astrophysical shocks: A theory of cosmic ray origin', Physics Reports 154 (1), 1-75 (1987), DOI: 10.1016/0370-1573(87)90134-7, https://doi.org/10.1016/0370-1573(87)90134-7
  7. I. D. Palmer, 'Transport coefficients of low-energy cosmic rays in interplanetary space', Reviews of Geophysics 20 (2), 335-351 (1982), DOI: 10.1029/RG020i002p00335, https://doi.org/10.1029/RG020i002p00335
  8. J.-P. Meyer, L. O'C. Drury and D. C. Ellison, 'Galactic Cosmic Rays from Supernova Remnants. I. A Cosmic-Ray Composition Controlled by Volatility and Mass-to-Charge Ratio', The Astrophysical Journal 487, 182-196 (1997), DOI: 10.1086/304599 (preprint arXiv:astro-ph/9704267), https://doi.org/10.1086/304599
  9. D. C. Ellison, L. O'C. Drury and J.-P. Meyer, 'Galactic Cosmic Rays from Supernova Remnants. II. Shock Acceleration of Gas and Dust', The Astrophysical Journal 487, 197-217 (1997), DOI: 10.1086/304580 (preprint arXiv:astro-ph/9704293), https://doi.org/10.1086/304580
  10. M. A. Malkov and L. O'C. Drury, 'Nonlinear theory of diffusive acceleration of particles by shock waves', Reports on Progress in Physics 64 (4), 429-481 (2001), DOI: 10.1088/0034-4885/64/4/201, https://doi.org/10.1088/0034-4885/64/4/201
  11. M. Ackermann et al. (Fermi-LAT Collaboration), 'Detection of the Characteristic Pion-Decay Signature in Supernova Remnants', Science 339 (6121), 807-811 (2013), DOI: 10.1126/science.1231160 (preprint arXiv:1302.3307), https://doi.org/10.1126/science.1231160
  12. W. R. Binns, M. H. Israel, E. R. Christian, A. C. Cummings, G. A. de Nolfo, K. A. Lave, R. A. Leske, R. A. Mewaldt, E. C. Stone, T. T. von Rosenvinge and M. E. Wiedenbeck, 'Observation of the 60Fe nucleosynthesis-clock isotope in galactic cosmic rays', Science 352 (6286), 677-680 (2016), DOI: 10.1126/science.aad6004, https://doi.org/10.1126/science.aad6004
  13. 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
  14. A. R. Bell, 'Turbulent amplification of magnetic field and diffusive shock acceleration of cosmic rays', Monthly Notices of the Royal Astronomical Society 353 (2), 550-558 (2004), DOI: 10.1111/j.1365-2966.2004.08097.x, https://doi.org/10.1111/j.1365-2966.2004.08097.x
  15. A. R. Bell, 'Cosmic Ray Origins in Supernova Blast Waves' (2014), arXiv:1412.7294, https://arxiv.org/abs/1412.7294
  16. D. Caprioli, A.-R. Pop and A. Spitkovsky, 'Simulations and Theory of Ion Injection at Non-relativistic Collisionless Shocks', The Astrophysical Journal Letters 798 (2), L28 (2014), DOI: 10.1088/2041-8205/798/2/L28 (preprint arXiv:1409.8291), https://doi.org/10.1088/2041-8205/798/2/L28
  17. A. Abramowski et al. (H.E.S.S. Collaboration), 'Acceleration of petaelectronvolt protons in the Galactic Centre', Nature 531, 476-479 (2016), DOI: 10.1038/nature17147 (preprint arXiv:1603.07730), https://arxiv.org/abs/1603.07730
  18. Z. Cao et al. (LHAASO Collaboration), 'Ultrahigh-energy photons up to 1.4 petaelectronvolts from 12 γ-ray Galactic sources', Nature 594, 33-36 (2021), DOI: 10.1038/s41586-021-03498-z, https://doi.org/10.1038/s41586-021-03498-z
  19. LHAASO Collaboration, 'An ultrahigh-energy γ-ray bubble powered by a super PeVatron', Science Bulletin 69 (4), 449-457 (2024), DOI: 10.1016/j.scib.2023.12.040, https://doi.org/10.1016/j.scib.2023.12.040
  20. Z. Cao et al. (LHAASO Collaboration), 'Ultrahigh-energy gamma-ray emission associated with black hole–jet systems', National Science Review 12 (12), nwaf496 (2025), DOI: 10.1093/nsr/nwaf496 (preprint arXiv:2410.08988), https://doi.org/10.1093/nsr/nwaf496
  21. 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