Self-mediation of runaway electrons via self-excited wave-wave and wave-particle interactions

Nonlinear dynamics of runaway electron induced wave instabilities can significantly modify the runaway distribution critical to tokamak operations. Here we present the first-ever fully kinetic simulations of runaway-driven instabilities towards nonlinear saturation in a warm plasma where collisional damping is subdominant. It is found that the slow-X modes grow an order of magnitude faster than the whistler modes, and they parametrically decay to produce whistlers much faster than those directly driven by runaways. These parent-daughter waves, as well as secondary and tertiary wave instabilities, initiate a chain of wave-particle resonances that strongly diffuse runaways to the backward direction. This reduces almost half of the current carried by high-energy runaways, over a time scale orders of magnitude faster than experimental shot duration. These results beyond quasilinear analysis may impact anisotropic energetic electrons broadly in laboratory, space and astrophysics.

💡 Research Summary

This paper presents the first fully kinetic particle‑in‑cell (PIC) investigation of runaway‑electron (RE) driven electromagnetic instabilities and their nonlinear saturation in a warm tokamak‑relevant plasma where collisional damping is negligible. Using a realistic set of parameters (J_RE = 2 MA m⁻², n_e = 0.6×10¹⁹ m⁻³, T_e = 320 eV, B = 1.45 T, |ω_ce|/ω_pe ≈ 1.84) the authors first compute a self‑consistent RE distribution with a drift‑kinetic Fokker‑Planck solver under a strong electric field (E = 65 E_c). This distribution is then fed into the VPIC code, which resolves three‑dimensional velocity space in a one‑dimensional spatial domain with periodic boundaries, realistic mass ratio, and a fine grid (Δx = 0.0125 d_e). Macro‑particle weighting is adjusted to give ten‑fold statistical resolution of the RE tail.

Linear dispersion analysis predicts that the slow‑X (extraordinary) branch at an angle θ = 40° has a growth rate ≈4×10⁻³ ω_pe, roughly an order of magnitude larger than the fastest whistler branch (≈10⁻⁴ ω_pe). The simulations confirm this: the slow‑X mode rapidly reaches large amplitude (t ω_pe ≈ 5 000) and, because of its high amplitude, undergoes parametric decay into two pairs of daughter waves. One pair consists of a low‑frequency whistler (ω ≈ 0.13 ω_pe, k d_e ≈ 0.35) and a high‑frequency slow‑X (ω ≈ 1.47 |ω_ce|, k d_e ≈ 1.88); the other pair produces higher‑k whistlers. These daughter whistlers appear much earlier and with larger amplitudes than the primary whistlers that would be directly driven by the REs.

The authors analyse wave‑particle interaction using the resonance condition ω − k_∥ v_ξ = n ω_ce/γ and a gradient operator ĤL acting on the distribution f(p,ξ). The sign of ĤL f determines whether a wave extracts energy from or gives energy to the particles, and the associated diffusion direction is given by a unit vector ĝ. Early in the simulation, the slow‑X mode and its daughters strongly diffuse the high‑energy tail (p > 10 m_ec) toward larger pitch angles (ξ ≈ −1), creating a “finger” in f(p,ξ). This finger both damps the slow‑X waves and seeds secondary whistler modes through n = −1 resonances at moderate momenta (p ≈ 5 m_ec). These secondary whistlers, in turn, undergo their own parametric decay, generating further whistlers that resonate via n = −2, −3, etc., propagating the diffusion cascade to even higher energies.

At a medium time scale (t ω_pe ≈ 3×10⁴) the cascade of secondary whistlers diffuses moderate‑energy REs backward, reinforcing the finger structure. At long times (t ω_pe ≈ 3×10⁵) a broad spectrum of whistlers produced by successive decays continuously scatters the high‑energy REs (p ≈ 20 m_ec, ξ ≈ −0.35) through a series of resonances (n = 2 → −2). The combined effect of all these resonances drives a substantial portion of the RE population toward negative parallel velocity, effectively reducing the net RE current by roughly 50 % within a time that is orders of magnitude shorter than the collisional current‑dissipation time (τ_coll ≈ 10⁷–10⁸ ω_pe⁻¹).

Because collisions are deliberately omitted (their timescales far exceed the simulated interval), the results isolate the intrinsic wave‑driven dynamics. The use of a long periodic domain (L_x = 1344 d_e) ensures a dense k‑space, allowing realistic representation of the continuous spectrum needed for parametric processes. The weighted macro‑particle technique guarantees sufficient statistics for the sparse RE tail without prohibitive computational cost.

The study yields several key insights: (1) Slow‑X modes dominate the early linear phase and act as a powerful energy reservoir that feeds whistler turbulence via rapid parametric decay. (2) Nonlinear wave‑wave coupling and the ensuing chain of wave‑particle resonances provide a highly efficient mechanism for RE pitch‑angle scattering and backward momentum transport, far exceeding predictions from quasilinear theory. (3) The resulting self‑mediated RE current mitigation occurs on sub‑microsecond scales in tokamak parameters, suggesting a potentially exploitable pathway for disruption mitigation or RE control. (4) The identified physics is generic to any plasma where relativistic electrons coexist with low‑frequency electromagnetic modes, implying relevance to solar flares, Earth’s magnetosphere, and intracluster media.

In conclusion, the authors demonstrate that fully kinetic simulations reveal a rich hierarchy of wave‑wave and wave‑particle interactions that can self‑regulate runaway electron populations far more rapidly than collisional processes. This advances our understanding of RE dynamics beyond quasilinear approximations and opens new avenues for both theoretical modeling and experimental diagnostics in fusion and space plasma environments.