Relativistic Electron Shock Drift Acceleration in Low Mach Number Galaxy Cluster Shocks

An extreme case of electron shock drift acceleration in low Mach number collisionless shocks is investigated as a plausible mechanism of initial acceleration of relativistic electrons in large-scale shocks in galaxy clusters where upstream plasma temperature is of the order of 10 keV and a degree of magnetization is not too small. One-dimensional electromagnetic full particle simulations reveal that, even though a shock is rather moderate, a part of thermal incoming electrons are accelerated and reflected through relativistic shock drift acceleration and form a local nonthermal population just upstream of the shock. The accelerated electrons can self-generate local coherent waves and further be back-scattered toward the shock by those waves. This may be a scenario for the first stage of the electron shock acceleration occurring at the large-scale shocks in galaxy clusters such as CIZA J2242.8+5301 which has well defined radio relics.

💡 Research Summary

**
This paper investigates a limiting case of electron shock‑drift acceleration (SDA) in low‑Mach number (M ≈ 2–5) collisionless shocks that are typical of the outskirts of galaxy clusters. The authors focus on the specific plasma conditions found in these environments: upstream electron temperatures of order 10 keV, plasma beta values β ≈ 0.1–10, and magnetization parameters σ ranging from 10⁻³ to 10⁻¹. Under such conditions the electron thermal speed becomes a sizable fraction of the speed of light, which dramatically changes the reflection and acceleration physics compared with the more familiar high‑Mach number supernova‑remnant shocks.

The theoretical part revisits the classic SDA mechanism. In the normal‑incidence frame (NIF) an electron drifts along the shock surface because of the magnetic‑field gradient, gaining energy from the motional electric field E₀. Transforming to the de Hoffmann‑Teller (HT) frame eliminates E₀, making the total energy W′ a constant of motion. The electron can be reflected by the magnetic mirror if its pitch angle α′₁ exceeds a loss‑cone angle α′_c, which depends on the upstream/downstream magnetic field ratio, the shock potential, and the Mach number. By incorporating relativistic corrections the authors derive a compact condition (Eq. 4–5) that relates the required velocity shell factor η_c to the upstream flow speed U₁, the shock obliquity Θ_Bn, the magnetic compression B₁/B₂, the plasma beta β_e, and the magnetization σ. When the HT speed V_HT ≈ U₁ tan Θ_Bn approaches the speed of light, η_c becomes comparable to the ratio of the electron thermal speed to c, allowing a non‑negligible fraction of the thermal population to be reflected and accelerated to relativistic energies.

To test these analytic expectations, the authors perform one‑dimensional electromagnetic particle‑in‑cell (PIC) simulations using the “injection‑wall” method to generate a steady shock. The upstream plasma is continuously injected from the left boundary, reflected at the right boundary, and the resulting interaction creates a shock that propagates leftward in the downstream rest frame. Simulation parameters are chosen to mimic observed cluster shocks: B₁ ≈ 3 µG, n₁ ≈ 10⁻⁴ cm⁻³, T_e = T_i ≈ 10 keV, Mach number M ≈ 3–4, and an obliquity angle Θ_Bn ≈ 45°–60°. The spatial grid resolves the electron Debye length (Δx ≈ 0.24 λ_De) and the time step satisfies the Courant condition.

The PIC results confirm the analytic picture. A small but significant fraction of incoming electrons (∼10⁻³ of the total) are reflected at the shock transition layer. Their trajectories show a drift along the shock surface, gaining energy over a timescale of order one ion gyroperiod (∼Ω_i⁻¹). The reflected electrons acquire Lorentz factors γ ≈ 2–5, i.e., relativistic energies, and form a distinct non‑thermal “upstream tail” in the electron distribution function. This tail drives the growth of coherent electromagnetic waves, identified as whistler‑mode fluctuations with wavelengths comparable to the electron inertial length. The waves scatter the reflected electrons back toward the shock, establishing a feedback loop: SDA → non‑thermal tail → wave growth → back‑scattering → further SDA. The wave amplitudes reach a few percent of the background magnetic field, sufficient to isotropize the reflected electrons and to trap them near the shock front.

The authors argue that this self‑generated turbulence provides the missing “injection” mechanism required for diffusive shock acceleration (DSA). In standard DSA theory, a pre‑existing suprathermal population must already be present to cross the shock repeatedly; the SDA‑induced tail naturally supplies such particles. Consequently, the low‑Mach number cluster shocks can seed DSA and produce the extended power‑law electron spectra inferred from radio relic observations.

The paper applies these findings to the well‑studied radio relic in the galaxy cluster CIZA J2242.8+5301. Observations of this relic show a spectral index consistent with a power‑law electron distribution of index ≈ 2.2–2.5, and a sharp spectral steepening downstream of the shock, indicative of radiative losses after acceleration. The simulated electron spectra and the characteristics of the self‑generated whistler waves match these observational constraints, supporting the scenario that relativistic SDA followed by wave‑mediated scattering initiates the DSA process in such environments.

In summary, the study demonstrates that even modest‑Mach number shocks in hot, weakly magnetized intracluster plasma can efficiently reflect and accelerate a fraction of thermal electrons to relativistic energies via shock‑drift acceleration. The reflected electrons self‑generate coherent electromagnetic fluctuations that scatter them back to the shock, thereby providing a natural injection channel for diffusive shock acceleration. This mechanism offers a plausible explanation for the origin of relativistic electrons responsible for the bright radio relics observed in merging galaxy clusters. Future work should extend the simulations to two and three dimensions, explore the dependence on plasma parameters more systematically, and compare directly with high‑resolution radio and X‑ray observations to further validate the model.