Proton-Helium Spectral Anomaly as a Signature of Cosmic Ray Accelerator

The much-anticipated proof of cosmic ray (CR) acceleration in supernova remnants (SNR) must hinge on full consistency of acceleration theory with the observations; direct proof is impossible because of the orbit scrambling of CR particles. The PAMELA orbital telescope revealed deviation between helium and proton CR spectra deemed inconsistent with the theory, since the latter does not differentiate between elements of ultrarelativistic rigidity. By considering an initial (injection-) phase of the diffusive shock acceleration (DSA), where elemental similarity does not apply, we demonstrate that the spectral difference is, in fact, a unique signature of the DSA. Collisionless plasma SNR shocks inject more He2+ relative to protons when they are stronger and so produce harder helium spectra. The injection bias is due to Alfven waves driven by the more abundant protons, so the He2+ ions are harder to trap by these waves because of the larger gyroradii. By fitting the p/He ratio to the PAMELA data, we bolster the DSA-case for resolving the century-old mystery of CR origin.

💡 Research Summary

The paper addresses a long‑standing puzzle in cosmic‑ray (CR) physics: the observed spectral difference between protons and helium nuclei reported by the PAMELA, ATIC and CREAM experiments. Standard diffusive shock acceleration (DSA) theory predicts that, at ultrarelativistic rigidities, particles of different charge‑to‑mass ratios should acquire identical power‑law spectra because the acceleration depends only on rigidity (R = pc/eZ). Consequently, the modest but statistically significant hardening of the helium spectrum (Δq ≈ 0.1) seemed to contradict the DSA paradigm and cast doubt on supernova remnants (SNRs) as the dominant Galactic CR sources.

Malkov, Diamond and Sagdeev propose that the discrepancy originates not in the steady‑state DSA phase but in the preceding injection stage, where only a tiny fraction of the thermal upstream plasma is able to re‑cross the shock and become eligible for further acceleration. This injection efficiency is highly sensitive to the shock Mach number (M) and to the geometry of the magnetic field relative to the shock normal. In a collisionless plasma, protons generate upstream Alfvénic turbulence that acts as a moving magnetic mirror. Because the trapping efficiency of this turbulence scales with particle gyroradius, He²⁺ ions—having twice the charge and twice the mass of protons at the same velocity—possess a gyroradius twice as large. They are therefore less efficiently trapped by the proton‑driven waves and have a higher probability of escaping upstream, especially when the shock is strong (large M).

The authors formalize this idea by introducing injection efficiencies ηₚ∝M^{-σₚ} and η_{He}∝M^{-σ_{He}}. Empirical and simulation‑based considerations suggest σₚ≈0.6–0.9 and σ_{He}≈0.15–0.3. The spectral index of the accelerated particles is the usual DSA result q(M)=4/(1−M^{-2}), which approaches q≈4 for strong shocks. The total number of CRs of species α (p or He) contributed by a single SNR over its Sedov‑Taylor phase is then

N_α(p)=∫{M{min}}^{M_{max}} η_α(M) (R_{inj}/R)^{q(M)} dM,

where R is the particle rigidity and R_{inj} is a reference injection rigidity (set to 1 GV). By changing variables to x=4 ln(R/R_{inj}) and integrating over t=M^{-2}, the authors obtain an analytic expression for the proton‑to‑helium ratio:

Nₚ/N_{He}=C · F(σₚ,x)/F(σ_{He},x),

with C fixed by the ambient p/He abundance. For sufficiently large x (i.e., R≫R_{inj}), the ratio simplifies to a power law in ln R:

Nₚ/N_{He}∝