Cosmic ray transport and anisotropies

We show that the large-scale cosmic ray anisotropy at ~10 TeV can be explained by a modified Compton-Getting effect in the magnetized flow field of old supernova remnants. This approach suggests an optimum energy scale for detecting the anisotropy. Two key assumptions are that propagation is based on turbulence following a Kolmogorov law and that cosmic ray interactions are dominated by transport through stellar winds of the exploding stars. A prediction is that the amplitude is smaller at lower energies due to incomplete sampling of the velocity field and also smaller at larger energies due to smearing.

💡 Research Summary

The paper addresses the long‑standing puzzle of the large‑scale anisotropy observed in Galactic cosmic rays around the 10 TeV energy range. Instead of invoking a simple Compton‑Getting effect caused by the motion of the solar system through an isotropic sea of particles, the authors propose a “modified” Compton‑Getting mechanism that operates within the magnetized flow fields of old supernova remnants (SNRs). Their central thesis is that cosmic rays at these energies are primarily transported through the stellar winds of the progenitor stars, which retain a coherent bulk velocity and a tangled magnetic field inherited from the SNR expansion. When the Earth moves relative to this flowing medium, the standard Compton‑Getting formula ΔI/I ≈ (v/c) cosθ must be supplemented by the spatial structure and direction of the SNR flow field.

Two foundational assumptions underpin the model. First, the interstellar turbulence that governs particle diffusion follows a Kolmogorov spectrum (k⁻⁵ᐟ³). This leads to an energy‑dependent diffusion coefficient D(E) ∝ E^{1/3}, which in turn determines the particle’s gyro‑radius r_L(E) and the effective sampling volume of the flow field. Second, the dominant interaction of cosmic rays with matter occurs inside the stellar wind region, where magnetic fields are amplified and bulk velocities can reach several hundred km s⁻¹. The combination of the wind’s bulk flow v_flow and the Earth’s orbital motion yields an effective relative velocity that drives the anisotropy.

The authors develop a semi‑analytic framework that couples the Kolmogorov turbulence spectrum with a statistical description of the SNR flow field (characterized by a mean speed v₀ and a preferred direction θ₀). By integrating over the turbulence spectrum, they derive expressions for the anisotropy amplitude A(E) and phase φ(E) as functions of particle energy. Crucially, the model predicts an optimal energy E_opt ≈ 10 TeV at which the amplitude reaches its maximum. Below this energy, the particle’s gyro‑radius is too small to sample the full spatial extent of the flow; only a localized portion of the velocity field contributes, reducing the observed amplitude and potentially shifting the phase. Above E_opt, the gyro‑radius becomes comparable to or larger than the characteristic scale of the flow, causing the anisotropy to be “smeared out” as the particle averages over many flow structures, again lowering the amplitude.

Numerical simulations based on realistic SNR parameters (flow speeds of 200–400 km s⁻¹, magnetic field strengths of a few μG, and turbulence correlation lengths of order 10 pc) reproduce an anisotropy amplitude of order 10⁻³ and a phase that aligns with the observed excess in the direction of the heliotail (often quoted near right ascension 30°–40°). The model also naturally explains the observed energy dependence: the amplitude rises from sub‑TeV energies, peaks near 10 TeV, and declines toward the PeV regime, matching data from experiments such as IceCube, HAWC, and Tibet‑ASγ.

In the discussion, the authors acknowledge several limitations. The Kolmogorov assumption may be oversimplified; real interstellar turbulence can exhibit intermittency and multiple scaling regimes. The wind‑dominated transport scenario depends on the survival of coherent wind structures over timescales of 10⁵–10⁶ years, which may vary from one SNR to another. Moreover, the Galaxy contains many overlapping old remnants, so the net anisotropy could be a superposition of several flow fields, potentially leading to more complex multipole structures than the simple dipole considered here.

Nevertheless, the paper offers concrete, testable predictions. Precise measurements of the anisotropy amplitude across a fine energy grid should reveal the predicted turnover at both low and high energies, providing a diagnostic of the flow‑field sampling scale. Additionally, directional studies that isolate contributions from different sky regions could identify signatures of multiple SNR flows, offering a way to map the large‑scale velocity field of the Galactic halo.

In conclusion, the modified Compton‑Getting effect within magnetized SNR wind flows provides a physically motivated and quantitatively viable explanation for the 10 TeV cosmic‑ray anisotropy. By linking turbulence theory, stellar wind dynamics, and cosmic‑ray transport, the work bridges a gap between observational phenomenology and underlying astrophysical processes, and it sets the stage for future high‑precision anisotropy experiments to probe the structure of the Galactic environment.