Microphysics of diffusive shock acceleration: impact on the spectrum of accelerated particles

Diffusive shock acceleration at collisionless shocks remains the most likely process for accelerating particles in a variety of astrophysical sources. While the standard prediction for strong shocks is that the spectrum of accelerated particles is universal, $f(p)\propto p^{-4}$, numerous phenomena affect this simple conclusion. In general, the non-linear dynamical reaction of accelerated particles leads to a concave spectrum, steeper than $p^{-4}$ at momenta below a few tens of GeV/c and harder than the standard prediction at high energies. However, the non-linear effects become important in the presence of magnetic field amplification, which in turn leads to higher values of the maximum momentum $p_{max}$. It was recently discovered that the self-generated perturbations that enhance particle scattering, when advected downstream, move in the same direction as the background plasma, so that the effective compression factor at the shock decreases and the spectrum becomes steeper. We investigate the implications of the excitation of the non-resonant streaming instability on these spectral deformations, the dependence of the spectral steepening on the shock velocity and the role played by the injection momentum.

💡 Research Summary

This paper investigates a critical microphysical effect in Diffusive Shock Acceleration (DSA) that can lead to steeper energy spectra of accelerated particles than the canonical p^-4 prediction. While non-linear DSA theories including the dynamical back-reaction of cosmic rays (CRs) typically predict concave spectra harder than p^-4 at high energies, observations of supernova remnants (SNRs) often require spectra steeper than p^-4. This work focuses on the recently identified role of self-generated magnetic fields and the resulting “postcursor” effect.

The standard test-particle DSA theory predicts a universal p^-4 spectrum for strong shocks. The study highlights that the excitation of the non-resonant streaming (Bell) instability upstream is crucial for magnetic field amplification, enabling particles to reach high energies (high p_max). However, hybrid simulations have revealed a key detail: the magnetic perturbations generated by this instability, once advected downstream of the shock, continue to drift in the same direction as the plasma flow at a speed close to the local Alfvén speed (v_A) calculated in the amplified field. This drift motion of the scattering centers creates a post-shock “postcursor” region.

This postcursor effect fundamentally alters the kinematics of the acceleration cycle. In the standard picture, the relevant velocity for particles diffusing back to the shock from downstream is the fluid speed u2. With the postcursor, the effective speed these particles must overcome becomes u2 + v_A (since the scattering centers themselves are moving away from the shock). This reduces the effective compression ratio experienced by the particles from R = u1/u2 to R_eff ≈ u1/(u2 + v_A). A lower compression ratio directly translates into a steeper power-law spectrum (f(p) ∝ p^(-α) with α > 4).

The authors develop a semi-analytical model to quantify this spectral steepening. The model self-consistently calculates the amplified magnetic field strength (both upstream and downstream), the downstream Alfvén speed v_A, the shock modification due to CR pressure (total compression ratio R_tot), and the resulting spectral index α. They explore the dependence of α on key parameters: shock velocity (v_sh), acceleration efficiency (ξ_CR, the fraction of shock ram pressure converted to CRs), and the injection momentum (p_inj). Their analysis shows that the steepening effect is more prominent for faster shocks and higher acceleration efficiencies, as both conditions lead to stronger magnetic field amplification and thus a larger v_A. They also find that when spectra become very steep (α approaching 5), the energy content becomes dominated by low-energy particles, making the precise value of p_inj a significant parameter in the problem.

In summary, the paper provides a compelling physical mechanism—rooted in the microphysics of self-generated magnetic turbulence and the consequent postcursor drift—that can naturally produce steeper-than-p^-4 particle spectra in efficient non-linear DSA scenarios. This offers a potential solution to the long-standing discrepancy between the standard DSA prediction and observations of steep spectra in various astrophysical shocks.