Elemental Spectra of Cosmic Rays: Two Galactic Sources and the Extragalactic Transition Beyond the Knee ✛ ✎ ✦
Context
The nature of the knee in the spectrum of galactic cosmic rays at ~3 PeV remains one of the central puzzles of particle astrophysics. Historically, it was thought it could be linked to the maximum energy of accelerators, propagation effects, or photonuclear processes. Already in the 1990s, hints emerged that the spectrum is shaped by multiple source populations. The classic idea of Fritz Zwicky (1934) that supernova remnants are the main accelerators remains dominant, but the precise configuration of sources and their contributions across different energy ranges were unclear until the advent of next-generation precision measurements.
Methods
The study used spectral data on protons, helium, and iron from the DAMPE (direct satellite measurements), LHAASO, and IceTop (ground-based extensive air shower arrays) experiments. High-resolution spectroscopy enabled the isolation of individual elements. The fit function was a sum of three components (two galactic and one extragalactic). To describe the spectrum of a single source, a new phenomenological form—'double log-parabola with exponential cutoff' (EDLP)—was proposed, which reproduces the observed flattening before the cutoff. The source population was modeled by convolving the single-source spectrum with a luminosity function, after which parameters were fitted to data using MINUIT (MIGRAD) and MCMC (emcee) methods. Fragmentation losses for heavy nuclei were additionally accounted for using a model dependence of the path length on rigidity.
Results
It is shown that the spectra of hydrogen, helium, and iron in the 40 GeV – 500 PeV range are satisfactorily described by the three-component model. The first galactic component has a convex shape with a low-energy spectral index of 2.76 (the same for all nuclei) and ends after the TeV knee at a rigidity of ~15×Z TeV. The second component dominates up to the PeV knee (3×Z PeV) and is responsible for the main break. The third, extragalactic, component appears at ~100 PeV for protons and becomes dominant above 1 EeV. It was found that the so-called 'ankles'—breaks toward flattening at 500×Z GeV and ~150 TeV (protons) or ~1 PeV (helium)—result from the intersection of different components and do not always obey a rigidity dependence due to differences in the chemical composition of the populations. The formal fit of the helium spectrum requires a high-energy cutoff of the first component at 700 TeV, but this detail is not mandatory when accounting for systematic uncertainties. Importantly, no additional third galactic population is needed to describe data in the 10 PeV – 1 EeV region.
Implications
The proposed interpretation means that the transition from galactic to extragalactic cosmic rays occurs 'seamlessly' just beyond the PeV knee, rather than at the traditionally considered 'ankle' in the all-particle spectrum at 4 EeV. This supports the 'disappointing model' for ultra-high energies and underscores that the Galaxy ceases to be an efficient proton accelerator already at energies of ~100 PeV. The convex shape of the first component may indicate nonlinear effects in acceleration processes (diffusive shock acceleration) or a contribution from local sources in the late stages of supernova evolution.
Future development
Future precision measurements of the spectra of CNO group nuclei (carbon, oxygen) and iron in the knee region by LHAASO and the planned SWGO will allow independent testing of the three-component model predictions and the rigidity dependence of breaks. Joint analysis of anisotropy and spectra in the 0.1–10 PeV range, especially accounting for the dipole phase, will help isolate the contribution of nearby sources. Numerical simulations of acceleration in supernova remnants with realistic nonlinear feedback and variable ambient composition will be an important step toward a theoretical foundation of the observed picture.
Impact
The results will impact cosmic-ray physics, high-energy gamma-ray astronomy (search for galactic PeVatrons), and the interpretation of neutrino telescope data, as the sharp transition to extragalactic protons alters expectations for the diffuse neutrino background.
Next steps
Verification of the predicted dependence of the 'running' transition on nuclear charge using new data on CNO and iron from LHAASO; refinement of fragmentation cross sections of heavy nuclei in the TeV region; modeling of secondary nuclei spectra (boron, beryllium) within the three-component model.
Key open problems
The proposed model directly links the origin of the knee to the exhaustion of galactic accelerator capabilities and the immediate appearance of extragalactic sources, necessitating a revision of traditional ideas about the 'galactic' range up to 1 EeV. Moreover, the convex shape of the first component and its possible connection to nonlinear acceleration at shock fronts in supernova remnants remains a subject of active debate, linking the knee problem to the intricacies of plasma physics.
The knowledge machine's view
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
The study constructs a three-component model of proton, helium, and iron spectra based on data from DAMPE, LHAASO, and IceTop. A new phenomenological form for the single-source spectrum (EDLP) is introduced, accounting for the luminosity function and nuclear fragmentation. It is shown that the transition from galactic to extragalactic cosmic rays occurs just beyond the PeV knee, and the 'ankle' is explained by the crossing of components.
Where its strength is
The flexible parameterization and simultaneous fitting of multiple elements allowed the identification of fine spectral structure and the role of different source populations, providing a basis for testing acceleration and propagation models.
Where it could go next
- It might be interesting to compare the derived source parameters with predictions from acceleration models in supernova remnants and other candidates, such as the work on why there is a knee in cosmic rays.based on: Why Cosmic Rays Have a Knee arXiv:2606.28141
- Related is the analysis of the origin of the Milky Way's neutrino emission; comparing with the expected contribution from galactic sources in your model could clarify their nature.based on: The Milky Way Speaks the Language of Neutrinos arXiv:2607.25966
- It is worth considering how the extragalactic component relates to particle acceleration in active galactic nuclei or other powerful sources discussed in a neighboring work on cosmic guns.based on: Cosmic Cannons: How Black Holes Accelerate Particles to Unimaginable Speeds arXiv:2606.00816
Worth trying
- Test whether the model reproduces the observed proton-to-helium flux ratio in the transition region using higher-resolution data.
- Estimate the contribution of the proposed galactic components to the diffuse gamma-ray background and compare with measurements, e.g., in the Galactic center region.
- Use MCMC methods to explore parameter degeneracies and the robustness of conclusions about the number of components.
Why this may matter
If the proposed interpretation is correct, it changes our understanding of the boundary between galactic and extragalactic cosmic rays and could impact the interpretation of data from neutrino and gamma-ray observatories.
Nearby works in our archive
- The Milky Way Speaks the Language of Neutrinossimilarity 0.64
- Why Cosmic Rays Have a Kneesimilarity 0.63
- Dark Matter Without Fine-Tuningsimilarity 0.63
- Cosmic Cannons: How Black Holes Accelerate Particles to Unimaginable Speedssimilarity 0.62
- Gamma Glow at the Galactic Center: Not Annihilation, but Collisionssimilarity 0.62
- Moon Shadows Cast by Giant Cosmic Rayssimilarity 0.62
The closest work in the archive sits at 0.64 by meaning — the usual level for our corpus.
A proton with an energy of 1 EeV in a typical galactic magnetic field (3 μG) has a Larmor radius of about 300 parsecs — almost half the disk thickness, so the Galaxy cannot confine such particles, and they almost certainly come from beyond its boundaries.
Key numbers
- TeV knee energy for protons: 15 TeV
- PeV knee energy for protons: 3 PeV
- first break (ankle) energy for protons: 500 GeV
- second break energy for protons: 150 TeV
- low-energy spectral index of the first component: 2.76