Electron Injection by Whistler Waves in Non-relativistic Shocks

Electron acceleration to non-thermal, ultra-relativistic energies (~ 10-100 TeV) is revealed by radio and X-ray observations of shocks in young supernova remnants (SNRs). The diffusive shock acceleration (DSA) mechanism is usually invoked to explain this acceleration, but the way in which electrons are initially energized or ‘injected’ into this acceleration process starting from thermal energies is an unresolved problem. In this paper we study the initial acceleration of electrons in non-relativistic shocks from first principles, using two- and three-dimensional particle-in-cell (PIC) plasma simulations. We systematically explore the space of shock parameters (the Alfv'enic Mach number, M_A, the shock velocity, v_{sh}, the angle between the upstream magnetic field and the shock normal, theta_{Bn}, and the ion to electron mass ratio, m_i/m_e). We find that significant non-thermal acceleration occurs due to the growth of oblique whistler waves in the foot of quasi-perpendicular shocks. The obtained electron energy distributions show power law tails with spectral indices up to alpha ~ 3-4. The maximum energies of the accelerated particles are consistent with the electron Larmor radii being comparable to that of the ions, indicating potential injection into the subsequent DSA process. This injection mechanism, however, requires the shock waves to have fairly low Alf'enic Mach numbers, M_A <~ 20, which is consistent with the theoretical conditions for the growth of whistler waves in the shock foot (M_A <~ (m_i/m_e)^{1/2}). Thus, if the whistler mechanism is the only robust electron injection process at work in SNR shocks, then SNRs that display non-thermal emission must have significantly amplified upstream magnetic fields. Such field amplification is likely achieved by the escaping cosmic rays, so electron and proton acceleration in SNR shocks must be interconnected.

💡 Research Summary

This paper tackles the long‑standing problem of how thermal electrons are injected into the diffusive shock acceleration (DSA) process at non‑relativistic astrophysical shocks, such as those found in young supernova remnants (SNRs). Using state‑of‑the‑art two‑ and three‑dimensional particle‑in‑cell (PIC) simulations, the authors systematically vary four key shock parameters: the Alfvénic Mach number (M_A), the shock speed (v_sh), the obliquity angle between the upstream magnetic field and the shock normal (θ_Bn), and the ion‑to‑electron mass ratio (m_i/m_e). Their primary discovery is that, in quasi‑perpendicular shocks (θ_Bn ≈ 70°–80°), oblique whistler waves grow in the foot region and provide the necessary electric fields to reflect and repeatedly scatter electrons. This interaction produces a non‑thermal tail in the electron energy distribution that follows a power‑law with spectral indices α ranging from about 3 to 4, depending on the exact parameters. The maximum electron energies attained correspond to Larmor radii comparable to those of the ions, indicating that the electrons have been pre‑accelerated to energies sufficient for subsequent DSA. Crucially, the growth of whistler waves—and therefore efficient electron injection—requires relatively low Alfvénic Mach numbers, specifically M_A ≲ 20, which matches the theoretical condition M_A ≲ (m_i/m_e)^{1/2}. If this whistler‑driven mechanism is the dominant injection channel in SNR shocks, then the upstream magnetic field must be substantially amplified to keep M_A in the required range. The authors argue that such amplification is naturally supplied by the streaming of escaping cosmic‑ray ions, linking electron and ion acceleration processes. The paper concludes that whistler‑mediated electron injection offers a self‑consistent explanation for the observed non‑thermal X‑ray and radio emission from SNRs, but it also implies that efficient electron acceleration is intimately tied to cosmic‑ray‑driven magnetic field amplification. Future work is suggested to explore higher mass‑ratio simulations, larger spatial domains, and direct comparison with multi‑wavelength observations to further validate the model.