Diffusive Shock Acceleration in Test-Particle Regime

We examine the test-particle solution for diffusive shock acceleration, based on simple models for thermal leakage injection and Alfv’enic drift. The critical injection rate, \xi_c, above which the cosmic ray (CR) pressure becomes dynamically significant, depends mainly on the sonic shock Mach number, M, and preshock gas temperature, T_1. In the hot-phase interstellar medium (ISM) and intracluster medium, \xi_c < 10^{-3} for shocks with M < 5, while \xi_c ~ 10^{-4}(T_1/10^6 K)^{1/2} for shocks with M > 10. For T_1=10^6 K, for example, the test-particle solution would be valid if the injection momentum, p_{inj} > 3.8 p_{th}. This leads to the postshock CR pressure less than 10% of the shock ram pressure. If the Alfv’en speed is comparable to the sound speed in the preshock flow, as in the hot-phase ISM, the power-law slope of CR spectrum can be significantly softer than the canonical test-particle slope. Then the CR spectrum at the shock can be approximated by the revised test-particle power-law with an exponential cutoff at the highest accelerated momentum, p_{max}(t). An analytic form of the exponential cutoff is also suggested.

💡 Research Summary

The paper revisits the classic test‑particle solution of diffusive shock acceleration (DSA) by incorporating two physically motivated refinements: a thermal‑leakage injection prescription and the effect of Alfvénic drift in the upstream plasma. The central question is under what conditions the test‑particle approximation remains valid, i.e., when the pressure contributed by the accelerated cosmic‑ray (CR) population stays dynamically insignificant compared with the bulk shock ram pressure.

To answer this, the authors define an injection efficiency ξ as the fraction of thermal particles that are injected into the acceleration process at a momentum p_inj. Using a simple model in which particles from the high‑energy tail of a Maxwellian distribution leak across the shock when their momentum exceeds p_inj, they derive a critical injection efficiency ξ_c. If ξ < ξ_c, the CR pressure P_CR ≪ ρ₁ u₁² (where ρ₁ and u₁ are the upstream density and velocity), and the shock structure is essentially unchanged – the test‑particle regime. Conversely, ξ > ξ_c would lead to a non‑linear modification of the shock.

The analysis shows that ξ_c depends primarily on the sonic Mach number M of the shock and the upstream temperature T₁. For weak to moderate shocks (M < 5) typical of the hot‑phase interstellar medium (ISM) and intracluster medium (ICM), ξ_c is below 10⁻³ regardless of temperature. For strong shocks (M > 10) the dependence on temperature becomes explicit:

\