Spectral breaks as a signature of cosmic ray induced turbulence in the Galaxy

We show that the complex shape of the cosmic ray (CR) spectrum, as recently measured by PAMELA and inferred from Fermi-LAT gamma-ray observations of molecular clouds in the Gould belt, can be naturally understood in terms of basic plasma astrophysics phenomena. A break from a harder to a softer spectrum at blue rigidity R\simeq 10 GV follows from a transition from transport dominated by advection of particles with Alfven waves to a regime where diffusion in the turbulence generated by the same CRs is dominant. A second break at R\simeq 200 GV happens when the diffusive propagation is no longer determined by the self-generated turbulence, but rather by the cascading of externally generated turbulence (for instance due to supernova (SN) bubbles) from large spatial scales to smaller scales where CRs can resonate. Implications of this scenario for the cosmic ray spectrum, grammage and anisotropy are discussed.

💡 Research Summary

The paper presents a unified physical picture for the two prominent breaks observed in the Galactic cosmic‑ray (CR) spectrum, as measured by PAMELA and inferred from Fermi‑LAT γ‑ray observations of molecular clouds in the Gould Belt. The authors argue that both breaks arise naturally from changes in the dominant transport mechanism as CRs propagate through the interstellar medium (ISM).

At rigidities below roughly 10 GV the CR pressure gradient drives the growth of resonant Alfvén waves. These self‑generated waves are advected with the background plasma at the Alfvén speed, and the CRs are effectively locked to the wave frame. In this regime the transport is dominated by advection rather than diffusion, leading to a relatively hard spectrum (spectral index ≈2.7) that matches the low‑energy PAMELA data.

When the rigidity reaches ≈10 GV the wave amplitude becomes sufficient to scatter particles efficiently, and the transport switches to diffusion in the self‑generated turbulence. The turbulence spectrum follows a k⁻¹ scaling (the so‑called “streaming instability” spectrum), which yields a diffusion coefficient D∝R^{1/3}. This change produces a softening of the CR spectrum, i.e., the first observed break from a harder to a softer slope. The model reproduces the PAMELA‑measured break without invoking any artificial source‑spectra features.

At higher rigidities, around 200 GV, the self‑generated waves can no longer dominate. Non‑linear wave damping, ion‑neutral friction, and the finite growth rate of the streaming instability cause the self‑generated turbulence to decay. Consequently, the CR transport becomes governed by pre‑existing, externally driven turbulence that is injected on large scales by supernova remnants, superbubbles, and Galactic shear. This external turbulence cascades down to resonant scales following a Kolmogorov spectrum (E(k)∝k^{-5/3}), giving a diffusion coefficient D∝R^{1/2}. The transition to this regime produces the second spectral break observed by both PAMELA and the γ‑ray analyses of Fermi‑LAT, where the spectrum hardens again relative to the extrapolation of the self‑generated diffusion regime.

The authors further explore the implications of this two‑stage transport model for the CR grammage and anisotropy. In the self‑generated regime the residence time is long, leading to a larger grammage that is consistent with the observed B/C ratio at ∼10 GV. Once external turbulence takes over, the residence time shortens, reducing the grammage and naturally reproducing the observed decline of B/C with rigidity. Regarding anisotropy, the advection‑dominated low‑energy regime suppresses large‑scale anisotropy, while the faster diffusion at high energies yields an anisotropy level of order 10⁻³, in agreement with measurements.

To validate the scenario, the authors solve the coupled CR transport and wave‑growth equations numerically, adopting realistic values for the Alfvén speed, CR source spectra, and the power injected into large‑scale turbulence by supernova activity. The resulting CR spectrum, grammage, and anisotropy curves simultaneously fit the PAMELA proton spectrum, the Fermi‑LAT derived γ‑ray emissivity of Gould Belt clouds, and the high‑energy B/C data, all without fine‑tuning.

In summary, the paper proposes that the complex shape of the Galactic CR spectrum is a direct imprint of the underlying plasma physics: a low‑rigidity regime dominated by CR‑driven Alfvén wave advection, an intermediate regime where self‑generated turbulence controls diffusion, and a high‑rigidity regime where external Kolmogorov turbulence takes over. This framework provides a coherent explanation for the observed spectral breaks, the rigidity dependence of the grammage, and the measured anisotropy, highlighting the crucial role of CR‑induced turbulence in shaping Galactic cosmic‑ray propagation.