🔭 astrophysics

Co-evolution of cosmic ray energy spectra, composition, and anisotropies

This paper proposes a four-component model involving Galactic, local, and two extra-galactic source populations to explain the co-evolution of cosmic ray energy spectra, composition, and anisotropies across a wide energy range, attributing specific spectral features like the knee, ankle, and suppression to distinct source characteristics and suggesting a Galactic-to-extra-galactic transition occurring around $10^8$ GeV.

Original authors: Bing-Qiang Qiao, Qiang Yuan, Yi-Qing Guo

Published 2026-05-27

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Bing-Qiang Qiao, Qiang Yuan, Yi-Qing Guo

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Mystery: Where Do Cosmic Rays Come From?

Imagine the universe is a giant, chaotic storm of invisible particles called cosmic rays. These are high-energy atoms (like protons and iron nuclei) zooming through space at nearly the speed of light. For over 100 years, scientists have been trying to figure out where they come from.

The problem is that space is filled with invisible magnetic fields that act like a giant, swirling blender. When these particles travel through space, the magnetic fields twist and turn them, scrambling their paths. By the time they hit Earth, they have forgotten their original address. It's like trying to find out where a specific raindrop fell by looking at a puddle in a hurricane; the direction is all mixed up.

The New Clue: A "Co-Evolution" Dance

This paper suggests a new way to solve the mystery. Instead of looking at just one thing (like how many particles hit us), the authors looked at three things happening at the same time:

The Energy Spectrum: How much energy the particles have.
The Composition: What kind of "stuff" they are made of (light elements like hydrogen vs. heavy elements like iron).
The Anisotropy: The direction they seem to be coming from (even if it's slightly scrambled).

The authors noticed that these three things change together in a very specific, synchronized way as the energy goes up. They call this "co-evolution."

The Analogy: Imagine a band playing a song.

The Energy is the volume.
The Composition is the mix of instruments (drums vs. guitars).
The Anisotropy is the direction the sound is coming from.

The authors realized that as the volume gets louder (higher energy), the mix of instruments changes, and the direction of the sound shifts in a predictable pattern. This pattern tells them that different "bands" (source groups) are taking turns playing the song.

The Solution: A Four-Component Model

To explain these synchronized changes, the authors built a model with four different source groups (like four different bands) contributing to the cosmic ray "music."

The Galactic Background (The "Crowd"):
- Who: A vast, diffuse collection of sources all over our Milky Way galaxy (like supernova remnants).
- Role: They provide the steady background noise. They are responsible for the "Knee" in the spectrum (a drop-off in energy) because they simply run out of power to accelerate protons beyond a certain point.
The Local Neighbor (The "Soloist"):
- Who: A single, relatively nearby source (likely the Geminga pulsar, a dead star spinning nearby).
- Role: Because it's close, it dominates the "music" at medium energies. It causes a sudden "hardening" (a spike in energy) and a "softening" (a drop) around 100 TeV. Crucially, because it's in a specific direction, it causes the "wind" of cosmic rays to shift direction, creating a phase reversal in the data.
Extragalactic Low-Energy (The "Distant Band 1"):
- Who: Sources outside our galaxy.
- Role: As the Galactic sources run out of steam, this group steps in. Their arrival causes a "hardening" around 30 PeV and marks the transition where cosmic rays start coming from outside our galaxy.
Extragalactic High-Energy (The "Distant Band 2"):
- Who: Another group of ultra-powerful sources outside our galaxy.
- Role: They take over at the highest energies. The "Ankle" (a bump in the spectrum) is actually the moment this second group takes over from the first extragalactic group. Finally, this group hits its own acceleration limit, causing the spectrum to cut off at the very highest energies.

How the Puzzle Pieces Fit Together

The paper argues that the complex shapes we see in the data aren't random glitches. They are the result of these four groups overlapping:

The "Knee" (3-4 PeV): This is where the main Galactic sources (Group 1) stop accelerating protons.
The "Second Knee" (200 PeV): This is where the Galactic sources stop accelerating heavy iron nuclei.
The "Ankle" (5 EeV): This isn't just a random bump; it's the moment the first extragalactic group (Group 3) hands the baton to the second extragalactic group (Group 4).
The Composition Changes: As the "bands" switch, the mix of particles changes. When a new group starts playing, it usually starts with lighter particles (protons), making the "average weight" of the cosmic rays drop. When they fade out, the heavier particles from the old group dominate again. This creates a "bump and dip" pattern in the composition data that perfectly matches the energy bumps.

The Directional Shift (Anisotropy)

The paper also explains why the direction of the cosmic rays seems to flip-flop.

At low energies, the Local Neighbor (Geminga) is so close that its "wind" blows toward us from a specific direction.
At medium energies, the Galactic Background takes over, and the wind shifts to point away from the Galactic Center.
At high energies, the Extragalactic sources take over, and the wind shifts again.

The "phase reversal" (where the direction flips 180 degrees) happens exactly when one source group stops dominating and another takes over.

The Bottom Line

The authors conclude that cosmic rays aren't coming from one single place or one single mechanism. Instead, the universe is a relay race. Different groups of sources (Galactic, Local, and two types of Extragalactic) pass the baton to each other as the energy increases.

By watching how the energy, composition, and direction change together, the authors have successfully mapped out this relay race, explaining the complex "knees," "ankles," and "bumps" in the cosmic ray spectrum as the natural result of these four distinct populations co-existing and co-evolving.

Technical Summary: Co-evolution of Cosmic Ray Energy Spectra, Composition, and Anisotropies

Problem Statement
The origin of cosmic rays (CRs) remains a fundamental unresolved problem in astrophysics, primarily due to the deflection of charged particles by magnetic fields, which obscures source locations, and the difficulty in determining mass composition via indirect ground-based measurements. While recent advances in space-borne and ground-based experiments have provided high-precision data on energy spectra, mass composition (mean logarithmic mass, $\langle \ln A \rangle$ ), and large-scale anisotropies, these observables have often been analyzed independently. A critical gap exists in understanding the correlated evolution of these multi-messenger observables across a vast energy range (from $\mathcal{O}(10)$ GeV to $10^{11}$ GeV). Specifically, previous models have struggled to simultaneously explain the fine spectral structures (hardening and softening features), the variation in composition, and the complex phase reversals in anisotropy without invoking multiple distinct source populations.

Methodology
The authors propose a unified four-component model to jointly interpret the energy spectra, composition, and anisotropies of CRs. The model decomposes the total CR flux into four distinct source populations:

Galactic Background Sources (Component A): A distributed population (e.g., supernova remnants) modeled with a spatially dependent propagation (SDP) scenario to account for diffusion variations between the Galactic disk and halo. The injection spectrum is an exponentially cutoff power-law in rigidity.
Galactic Nearby Source (Component B): A single local source, illustrated by the pulsar Geminga, injected as a point source with a specific age and distance. Its contribution is calculated using a Green's function method.
Extragalactic Low-Energy Sources (Component C): A population responsible for the transition from Galactic to extragalactic dominance.
Extragalactic High-Energy Sources (Component D): A second extragalactic population responsible for the highest energy cutoff.

The model employs the DRAGON code to solve the transport equation for Galactic components, incorporating diffusion, convection, and re-acceleration. For extragalactic components, a phenomenological diffusion coefficient is used to fit anisotropy data above EeV energies. The total observed spectrum is the algebraic sum of the four components, while the anisotropy is the vector sum of their individual streaming directions. The model parameters (spectral indices, cutoff rigidities, normalizations) are tuned to reproduce the "hardening-softening" groups observed in the data, the evolution of $\langle \ln A \rangle$ , and the dipole amplitude and phase of anisotropies.

Key Contributions and Results
The study identifies four correlated groups of "hardening-softening" spectral features across the energy spectrum, linking each to specific source transitions and physical limits:

Group 1 ( $\sim 10$ GeV – $100$ TeV): Characterized by a hardening at $\sim 200$ GV and a softening at $\sim 15$ TV. The model attributes this to the superposition of the Galactic background (A) and the nearby source (B). The nearby source dominates the flux and anisotropy direction between $\sim 100$ GeV and $100$ TeV, explaining the observed phase reversal of the dipole anisotropy (shifting from the Galactic Center direction toward the nearby source direction).
Group 2 ( $\sim 100$ TeV – $10$ PeV): Includes the "knee" at $\sim 3\text{--}4$ PeV. This feature arises as the contribution of the nearby source (B) diminishes, allowing the Galactic background (A) to dominate again, until reaching the maximum acceleration energy of protons by Galactic sources.
Group 3 ( $\sim 10$ PeV – $1$ EeV): Features a hardening at $\sim 30$ PeV and the "second knee" at $\sim 200$ PeV. The hardening is caused by the emergence of the first extragalactic component (C), while the second knee corresponds to the maximum acceleration energy of iron nuclei by Galactic sources (A). The transition in anisotropy phase occurs around $1$ EeV as the streaming from component C overtakes component A.
Group 4 ( $\sim 1$ EeV – $100$ EeV): Includes the "ankle" at $\sim 5$ EeV and the spectral suppression at $\sim 50$ EeV. The ankle is interpreted as the transition from extragalactic component C to component D. The highest energy suppression is attributed to the acceleration limit of component D, rather than the GZK process.

Composition and Anisotropy Correlations
The model successfully reproduces the observed "bump-dip" structures in the mean logarithmic mass ( $\langle \ln A \rangle$ ) that correlate with the spectral hardening and softening. The authors argue that as new components appear (typically proton-dominated), the average mass lightens, and as they fade or reach cutoffs, the composition heavyens. Similarly, the large-scale anisotropy amplitude and phase show distinct transitions at the characteristic energies of these source switches. Notably, the model suggests the transition from Galactic to extragalactic origin occurs around $\mathcal{O}(10^8)$ GeV, which is lower than the traditional ankle energy.

Significance and Claims
The paper claims that the co-evolution of spectra, composition, and anisotropies provides strong evidence for a multi-population origin of cosmic rays. By unifying these observables, the four-component model offers a consistent explanation for the complex structures in the CR spectrum without relying on a single homogeneous source population. The authors emphasize that while current measurements of composition and anisotropy at high energies (above tens of PeV) and low energies (below TeV) carry significant uncertainties, the correlated trends observed support the proposed scenario. The work highlights the necessity of future high-precision measurements from both direct (e.g., HERD, VLAST) and indirect detection experiments to further constrain the source populations and validate the specific transition energies and anisotropy phase behaviors predicted by the model.