Reinterpretation of the Fermi acceleration of cosmic rays in terms of the ballistic surfing acceleration in supernova shocks

The applicability of first-order Fermi acceleration in explaining the cosmic ray spectrum has been reexamined using recent results on shock acceleration mechanisms from the Multiscale Magnetospheric mission in Earth’s bow shock. It is demonstrated that the Fermi mechanism is a crude approximation of the ballistic surfing acceleration (BSA) mechanism. While both mechanisms yield similar expressions for the energy gain of a particle after encountering a shock once, leading to similar power-law distributions of the cosmic ray energy spectrum, the Fermi mechanism is found to be inconsistent with fundamental equations of electrodynamics. It is shown that the spectral index of cosmic rays is determined by the average magnetic field compression rather than the density compression, as in the Fermi model. It is shown that the knee observed in the spectrum at an energy of 5x10^{15} eV could correspond to ions with a gyroradius comparable to the size of shocks in supernova remnants. The BSA mechanism can accurately reproduce the observed spectral index s = -2.5 below the knee energy, as well as a steeper spectrum, s = -3, above the knee. The acceleration time up to the knee, as implied by BSA, is on the order of 300 years. First-order Fermi acceleration does not represent a physically valid mechanism and should be replaced by ballistic surfing acceleration in applications or models related to quasi-perpendicular shocks in space. It is noted that BSA, which operates outside of shocks, was previously misattributed to shock drift acceleration (SDA), which operates within shocks.

💡 Research Summary

The paper revisits the long‑standing problem of how Galactic cosmic rays acquire their characteristic power‑law energy spectrum, especially the change of slope at the “knee” around 5 × 10¹⁵ eV. Using high‑resolution measurements from the Magnetospheric Multiscale (MMS) mission of Earth’s bow shock, the author identifies four fundamental processes that operate in collisionless shocks: stochastic wave energization (SWE), transit‑time thermalization (TTT), quasi‑adiabatic heating (QAH), and ballistic surfing acceleration (BSA). The focus of the study is on BSA, a mechanism that becomes effective when a particle’s gyroradius r_c greatly exceeds the shock width D (r_c ≫ D). In this regime the particle does not feel the magnetic gradient inside the shock ramp; instead it drifts in the large‑scale convection electric field E_y that exists upstream and downstream of the shock.

Starting from the relativistic momentum equation dp/dt = q(E + v × B) in a shock‑frame where x is the shock normal and y points along E_y, the author derives the energy change ΔK = q E·v dt. By separating the integration domain into three parts—outside the ramp (|x| ≥ D), inside the ramp (|x| < D) for shock‑drift acceleration (SDA), and the contribution of the cross‑shock electric field E_S—the total energy gain per gyration is expressed as

ΔK ≈ 2 q E_y (r_c^up − r_c^down).

For relativistic particles K ≈ pc and r_c = p_⊥/(qB), this leads to a fractional energy increase

ΔK/K = (c_B − 1) V_u/c,

where c_B = B_d/B_u is the magnetic‑field compression ratio and V_u = E_y/B_u is the upstream convection speed. This expression (Eq. 4 in the paper) is mathematically identical in form to the classic first‑order Fermi acceleration formula, but crucially depends on magnetic compression rather than the density compression c_N used in the conventional DSA derivation. The author argues that the standard DSA derivation mixes inertial frames and therefore computes a scalar Lorentz‑transformation energy difference, not the work done by the electric field in a single frame. Consequently, the physically correct acceleration law is the BSA expression that involves only c_B.

Assuming that after each interaction a particle remains in the acceleration region with probability P ≈ 1 − V_u/c, and that the energy multiplication factor per interaction is h = 1 + (1 − 1/c_B) V_u/c, the steady‑state particle distribution follows

N(K) ∝ K^s, s = ln P/ln h − 1 = −(2c_B − 1)/(c_B − 1).

For a typical magnetic compression c_B ≈ 4, the model yields s ≈ −2.5, matching the observed spectrum below the knee. Above the knee the gyroradius becomes comparable to the size of a supernova‑remnant shock, the acceleration efficiency drops, and the effective index steepens to s ≈ −3, again in agreement with observations.

The paper also estimates the acceleration time. Using typical supernova‑remnant parameters (E_y ≈ V_u B_u, shock size L ≈ 1 pc), the fractional energy gain per gyration is of order (c_B − 1) V_u/c ≈ 10⁻³. Reaching 5 × 10¹⁵ eV therefore requires a few hundred gyrations, corresponding to an acceleration timescale of roughly 300 years (≈10⁸ s). This is far shorter than the ∼10⁴‑year timescales often quoted for diffusive shock acceleration, and it aligns with the age constraints inferred from cosmic‑ray composition studies.

Importantly, the author distinguishes BSA from shock‑drift acceleration (SDA). SDA relies on ∇B drifts within the shock ramp and is only significant for particles with r_c ≪ D. In contrast, BSA operates entirely outside the ramp, where the magnetic field is essentially uniform and the only driver is the large‑scale convection electric field. The MMS observations show that high‑energy ions gain energy continuously while moving upstream and downstream, with negligible energy change at the actual shock crossing itself—direct evidence that the dominant energization mechanism is BSA, not the reflection‑based processes assumed in the classic Fermi picture.

In conclusion, the study provides a self‑consistent, electrodynamics‑based acceleration mechanism that reproduces the observed cosmic‑ray spectrum, explains the knee, and yields realistic acceleration timescales. It challenges the physical validity of first‑order Fermi (DSA) acceleration for quasi‑perpendicular shocks and argues that ballistic surfing acceleration should replace it in theoretical models and numerical simulations of high‑energy astrophysical shocks.