Fundamentals

The Cosmic Ray Energy Spectrum

The cosmic ray spectrum is one of the most remarkable features in all of physics — a nearly perfect power law spanning 12 orders of magnitude in energy, from 10⁹ eV to beyond 10²⁰ eV. Yet hidden in this smooth decline are subtle features — the "knee," the "ankle," and the "cutoff" — each marking a fundamental transition in cosmic ray physics.

The Power Law

Count cosmic rays hitting Earth's atmosphere and plot the number at each energy. The result is startlingly simple: a power law that stretches across an almost unimaginable range.

📊 The Spectrum Formula

dN/dE ∝ E

where γ (the spectral index) is approximately 2.7 below the knee, 3.1 between knee and ankle, and ~2.6 above the ankle.

This power-law form isn't coincidental — it's a signature of the underlying acceleration mechanisms. Fermi acceleration at shock fronts naturally produces power-law spectra with indices near 2.0-2.3. The observed steeper index (~2.7) results from energy-dependent escape and propagation through the Galaxy.

Flux: From Common to Rare

The power law means lower energies are vastly more common than higher energies. Here's the cosmic ray flux at different energies:

🌊 Cosmic Ray Flux

Energy Flux Detection
10⁹ eV (GeV) ~1000/m²/second Satellite
10¹² eV (TeV) ~1/m²/second Satellite
10¹⁵ eV (PeV) - Knee ~1/m²/year Ground array
10¹⁸ eV (EeV) - Ankle ~1/km²/year Large array
10¹⁹ eV (10 EeV) ~1/km²/decade Giant array
10²⁰ eV (100 EeV) ~1/km²/century Giant array

The flux drops by a factor of 10 billion between GeV and EeV energies. This is why detecting UHECRs requires enormous collection areas — the Pierre Auger Observatory's 3,000 km² is barely sufficient to catch a few events per year above 10²⁰ eV.

The Knee (~3 × 10¹⁵ eV)

Around 3 PeV (3 × 10¹⁵ eV), the spectrum steepens from γ ≈ 2.7 to γ ≈ 3.1. This "knee" was discovered in 1958 by the Moscow State University group.

What Causes the Knee?

Several explanations have been proposed:

Galactic accelerator limit: Supernova remnants — the primary Galactic cosmic ray sources — may reach their maximum acceleration energy around the knee. Above this energy, particles escape the accelerating shock.

Rigidity-dependent composition: Galactic magnetic fields confine particles based on rigidity (E/Z). Protons might begin leaking from the Galaxy at lower energies than heavier nuclei. The knee could mark where protons escape.

Evidence: Measurements show the composition becomes heavier above the knee, consistent with lighter elements dropping out first. A "second knee" around 10¹⁷ eV might mark where iron nuclei begin escaping.

The Second Knee (~10¹⁷ eV)

A subtler steepening occurs around 10¹⁷ eV. This is often interpreted as the endpoint of the Galactic cosmic ray population — the energy above which even heavy nuclei can no longer be confined to the Galaxy.

The region between 10¹⁷ and 10¹⁸ eV is transitional: Galactic cosmic rays are fading while extragalactic cosmic rays begin to dominate.

The Ankle (~5 × 10¹⁸ eV)

Around 5 EeV, the spectrum flattens from γ ≈ 3.3 to γ ≈ 2.6. This "ankle" feature is thought to mark the transition to a predominantly extragalactic cosmic ray population.

Two Interpretations

Dip model: The ankle is caused by energy losses of extragalactic protons interacting with the cosmic microwave background via pair production (p + γ_CMB → p + e⁺ + e⁻). This "pair production dip" creates a hardening in the observed spectrum.

Mixed composition model: The ankle marks where a new, harder extragalactic component overtakes the steeper Galactic component. This requires less fine-tuning but predicts a mixed (heavier) composition.

Current data (particularly from Auger) favor the mixed composition model — the composition at the ankle appears to be mixed rather than pure proton.

The Instep (~13 EeV)

A more recently identified feature: the spectrum shows a slight steepening around 13 EeV before continuing to the suppression. Auger calls this the "instep."

Its origin is unclear — it might relate to propagation effects, source properties, or composition changes.

The GZK Suppression (~5 × 10¹⁹ eV)

Above about 50 EeV, the spectrum drops sharply. This suppression was predicted in 1966 by Kenneth Greisen, Georgiy Zatsepin, and Vadim Kuzmin — before it was observed.

The GZK Mechanism

Ultra-high-energy protons interact with cosmic microwave background (CMB) photons:

p + γ_CMB → Δ⁺ → p + π⁰ (or n + π⁺)

Each interaction costs the proton about 20% of its energy. After traveling ~100 Mpc, a 10²⁰ eV proton has degraded below the threshold.

This creates a "GZK horizon" — particles above the threshold can only come from relatively nearby sources (<100-200 Mpc).

Suppression vs. Cutoff

Importantly, the observed suppression could have two causes:

  1. GZK effect: Propagation losses limit the distance of observable sources
  2. Source exhaustion: Accelerators simply don't produce particles above this energy

Both effects are likely operating. The spectrum suppresses because sources run out of energy AND because high-energy particles can't travel far.

Distinguishing between them requires composition measurements and multi-messenger observations.

Beyond the Suppression

Despite the suppression, cosmic rays above 10²⁰ eV do exist — the Oh-My-God particle (3 × 10²⁰ eV) and Amaterasu (2.4 × 10²⁰ eV) prove that.

These trans-GZK events require either:

Reading the Spectrum

The cosmic ray spectrum encodes information about:

Decoding these contributions requires combining spectrum measurements with composition studies, anisotropy searches, and theoretical modeling.

The All-Particle Spectrum

The spectrum shown above is the "all-particle" spectrum — all cosmic rays regardless of composition. Experiments also measure spectra of individual elements (hydrogen, helium, carbon, iron, etc.) at lower energies where composition can be determined.

Each element has its own spectrum, knee, and maximum energy. The all-particle spectrum is the sum of these individual contributions.

At the highest energies, we can't identify individual nuclei event-by-event — only statistical composition from air shower properties.

Spectral Measurements

Different experiments measure different energy ranges:

📊 Who Measures What

  • Space-based (CALET, DAMPE, AMS-02): GeV to ~100 TeV, direct composition
  • KASCADE/KASCADE-Grande: 10¹⁵ to 10¹⁸ eV, knee region
  • IceTop/IceCube: 10¹⁵ to 10¹⁸ eV, knee region
  • Telescope Array: >10¹⁸ eV, ankle and suppression
  • Pierre Auger: >10¹⁸ eV, ankle and suppression

Connecting measurements across experiments is challenging — each has different systematic uncertainties, energy scales, and composition sensitivities.

What the Spectrum Tells Us

After a century of measurements, the cosmic ray spectrum has revealed:

Two populations: Galactic cosmic rays dominate below ~10¹⁸ eV; extragalactic above.

Fermi acceleration works: The power-law form confirms shock acceleration as the primary mechanism.

The universe is opaque: The GZK suppression confirms that ultra-high-energy particles interact with the CMB.

Extreme accelerators exist: Something produces particles beyond 10²⁰ eV — the most energetic particles in the universe.

Yet fundamental questions remain: What sources dominate at each energy? What is the composition? Why does the spectrum have the precise features it does?

The spectrum is a message from the universe's most powerful accelerators. We're still learning to read it.

Continue Reading