Revisiting the role of the streaming instability for the cosmic-ray spectrum in the GeV to TeV range

A complete understanding of the cosmic-ray energy spectrum remains a challenge to theory that must be met by comprehensive modeling efforts. One of these is the subject of the present study, namely, an explanation of the recently discovered spectral hardening at $\sim 300$ GeV with self-consistently treated cosmic-ray diffusion, where self-generated waves resulting from the streaming instability impact the diffusion of high-energy particles. We revisit the corresponding model by Blasi et al. (2012), perform an extensive parameter study, and determine an optimal range of parameters that best fit the cosmic-ray data. We conclude that self-consistently treated cosmic-ray transport remains a competitive alternative to explain the spectral hardening of the cosmic-ray energy spectrum at a few hundred GeV.

💡 Research Summary

This paper presents a comprehensive re-examination of the self-consistent cosmic-ray (CR) diffusion model proposed by Blasi et al. (2012) to explain the observed spectral hardening in the CR energy spectrum at around 300 GeV. The central hypothesis is that the interaction between CRs and self-generated Alfvén waves via the streaming instability can modify the diffusion coefficient in a way that naturally produces a spectral break in the GeV to TeV range, offering an alternative to explanations based on local sources or source spectral features.

The authors construct a coupled, nonlinear model consisting of two core equations: 1) A steady-state transport equation for CR protons (Eq. 1), describing their diffusion, advection with Alfvén waves, and adiabatic effects, with a power-law injection at the Galactic disk. 2) An evolution equation for the spectrum of Alfvén waves (Eq. 8), which includes wave cascading (either Kolmogorov- or Kraichnan-type phenomenology) and wave growth driven by the CR gradient via the streaming instability (Eq. 10). The critical link is that the CR diffusion coefficient (Eq. 6) depends inversely on the wave spectral power, while the wave growth rate depends on the gradient of the CR distribution. This creates a feedback loop where CRs affect their own diffusion.

To solve this system, the authors derive implicit integral solutions for both the CR distribution at the disk (Eq. 13) and the wave spectrum (Eq. 15). They implement a numerical iterative scheme using trapezoidal integration over logarithmically spaced grids in momentum and wavenumber, continuing until convergence is achieved.

The study first provides an analytical estimate for the break energy (Eq. 16), showing its dependence on key parameters: injection normalization (A), magnetic field strength (B0), turbulence level (η), Kolmogorov constant (c_K), and injection slope (α). A quantitative parameter study follows, visually demonstrating how varying each parameter shifts the resulting CR spectrum and wave spectrum (Figures 3 & 4). A key finding is that the Kraichnan-type cascade typically predicts a break at higher energies than the Kolmogorov-type.

The core result of the paper is the determination, through an extensive parameter study (hinting at methods like MCMC fitting), of an optimal parameter range that best reproduces the observed CR proton spectrum. The authors conclude that despite the model’s simplifications—such as the one-dimensional flux-tube approximation and the exclusive focus on protons—the self-consistent treatment of CR transport, where waves are generated by the CRs themselves, remains a viable and competitive framework for explaining the spectral hardening feature. It successfully demonstrates that the interplay between wave growth and cascade can imprint a observable signature on the CR spectrum without invoking new source classes or exotic injection spectra. Future work is suggested to include heavier nuclei and more complex galactic geometry.