Particle injection in three-dimensional relativistic magnetic reconnection

Relativistic magnetic reconnection has been proposed as an important nonthermal particle acceleration (NTPA) mechanism that generates power-law spectra and high-energy emissions. Power-law particle spectra are in general characterized by three parameters: the power-law index, the high-energy cutoff, and the low-energy cutoff (i.e., the injection energy). Particle injection into the nonthermal power law, despite also being a critical step in the NTPA chain, has received considerably less attention than the subsequent acceleration to high energies. Open questions on particle injection that are important for both physical understanding and astronomical observations include how the upstream magnetization~$σ$ influences the injection energy and the contributions of the known injection mechanisms (i.e., direct acceleration by the reconnection electric field, Fermi kicks, and pickup acceleration) to the injected particle population. Using fully kinetic particle-in-cell simulations, we uncover these relationships by systematically measuring the injection energy and calculating the contributions of each acceleration mechanism to the total injected particle population. We also present a theoretical model to explain these results. Additionally, we compare two- and three-dimensional simulations to assess the impact of the flux-rope kink and drift-kink instability on particle injection. We conclude with comparisons with previous work and outlook for future work.

💡 Research Summary

This paper presents a comprehensive investigation of particle injection—the first step of non‑thermal particle acceleration (NTPA)—in relativistic magnetic reconnection, using fully kinetic particle‑in‑cell (PIC) simulations in both two and three dimensions. The authors begin by defining an injection criterion: a particle is considered injected when its gyroradius exceeds half the thickness of the reconnection layer at an X‑point (r_g ≥ δ/2). By expressing the layer thickness δ in terms of the relativistic electron skin depth (δ ≈ d_e,rel ≈ c/ω_pe √Γ) and introducing the mean Lorentz factor inside the layer Γ = Γ_th + k σ, they derive analytic expressions for the injection energy γ_inj. Two regimes emerge: a thermally‑dominated regime (1 ≲ σ_h ≲ k⁻¹) where γ_inj ≈ σ (1 + b_g²)⁻¹ᐟ² Γ_th, and a magnetically‑dominated regime (σ_h ≫ k⁻¹) where γ_inj ≈ σ (1 + b_g²)⁻¹ᐟ² k. Both scalings predict that stronger guide fields (larger b_g) lower the required injection energy because the effective magnetic field strength in the gyroradius formula is reduced.

The simulation campaign explores pair plasmas (electron‑positron) with upstream magnetizations σ = 10, 30, 100, both with zero guide field and with a modest guide field (b_g = 0.3). Harris current sheets are initialized, and particle trajectories are recorded to separate three injection mechanisms: (i) direct acceleration by the reconnection electric field E_rec near X‑points, (ii) Fermi “kicks” arising from the curvature drift aligned with the motional electric field (E_m = −u × B) as newly‑reconnected field lines relax, and (iii) pickup acceleration that occurs when particles cross the separatrix and experience a sudden increase in perpendicular momentum p′_⊥. Analytic estimates give W_direct ≈ η_rec β_A ω_ce Δt, W_Fermi ≈ 2σ/(1 + σ b_g²), and a comparable order‑of‑magnitude for pickup work.

Results confirm the theoretical scalings: γ_inj grows roughly linearly with σ and follows the (1 + b_g²)⁻¹ᐟ² dependence. For σ = 100 the measured injection energy is ≈ 30 m_ec², marking the transition from the thermally‑dominated to the magnetically‑dominated regime. Regarding mechanism contributions, direct acceleration accounts for ~20 % of injected particles at low σ, while Fermi kicks dominate with 50‑60 % across all σ values, and pickup acceleration contributes a steady ~10‑15 %. As σ increases, the relative importance of the Fermi process rises, consistent with the larger curvature drift speeds in higher‑σ outflows.

Three‑dimensional runs reveal the development of flux‑rope kink and drift‑kink instabilities, which distort the X‑point geometry and reduce the volume of direct‑acceleration regions. Consequently, the overall injection fraction drops by ≈ 12 % relative to 2‑D, but the γ_inj scaling with σ and b_g remains essentially unchanged. When a modest guide field is present, the kink instability is suppressed, and the 3‑D injection efficiency matches the 2‑D case.

The authors compare their findings with prior work that focused mainly on high‑energy acceleration, highlighting that the injection stage controls the low‑energy cutoff of the power‑law spectrum and therefore the dynamic range of observable emission. They argue that a proper description of NTPA must include the injection physics elucidated here. The paper concludes that (1) the injection energy follows γ_inj ∝ σ (1 + b_g²)⁻¹ᐟ² q(σ), (2) Fermi kicks are the primary injector across a broad σ range, and (3) three‑dimensional effects modestly diminish injection efficiency without altering the fundamental scaling. Future directions include extending the study to electron‑ion plasmas, incorporating radiative cooling, and linking the injection‑controlled spectra to astrophysical observations such as blazar flares and pulsar wind nebulae emission.