Nonthermal electron acceleration in turbulent post-flare coronal loops

Technical Summary: Non-thermal electron acceleration in turbulent, post-flare coronal loops

Problem Statement
The generation of energetic, non-thermal electrons during solar flares is a critical mechanism for energy transport from the corona to the chromosphere, manifesting as intense X-ray emission. While turbulence in post-flare loops, particularly driven by Kelvin–Helmholtz instabilities (KHI), has been proposed as a mechanism for trapping and accelerating electrons, the precise energization processes remain incompletely understood. Previous work (e.g., Bacchini et al. 2024) demonstrated that KHI-driven turbulence can confine electrons but did not fully document the specific mechanisms of particle energization. Furthermore, standard numerical approaches using the Guiding-Centre Approximation (GCA) often fail to strictly conserve energy in turbulent environments, leading to inconsistencies in quantifying acceleration rates and the relative roles of magnetic effects.

Methodology
The authors employed a test-particle approach coupled with a 2.5D (two spatial dimensions, three vector components) ideal Magnetohydrodynamic (MHD) model of a time-evolving post-flare coronal looptop. The MHD background, simulated using the MPI-AMRVAC 3.0 code, models chromospheric evaporation flows meeting at the looptop to generate turbulence via KHI, without explicitly modeling the primary reconnection current sheet.

A key methodological contribution is the development of an improved GCA formalism. The authors identified that standard GCA implementations, which enforce magnetic moment conservation ($\mu = \text{const}$), fail to conserve total particle energy in the perpendicular direction when integrated numerically. To resolve this, they reformulated the equations of motion to explicitly solve for perpendicular energy evolution (derived from the magnetic moment conservation condition) rather than enforcing $\mu$ conservation directly. This approach sacrifices exact magnetic moment invariance (allowing $\mu$ to oscillate slightly) to ensure strict energy conservation, which is essential for accurately analyzing energization in turbulent plasmas.

Approximately $10^6$ test electrons, initialized with a Maxwellian distribution at 20 MK, were evolved within the turbulent MHD background starting from a developed turbulence state ($t = 2.5t_0$). The study focused on the energy dynamics up to $t = 3.2t_0$, analyzing the contributions of parallel electric fields (absent in ideal MHD), curvature effects, and perpendicular gradient effects.

Key Results

1. Energy Conservation Correction: In a test case involving Alfvén wave propagation, the standard GCA approach was shown to fail in conserving perpendicular energy, whereas the new energy-conserving formulation successfully maintained energy balance, with magnetic moment $\mu$ varying only mildly ($\sim 10^{-4}$).
2. Particle Energization: Electrons in the turbulent looptop developed supra-thermal energy distributions with tails compatible with hard X-ray (HXR) emission. The bulk of the distribution remained within soft X-ray energies ($\gamma \lesssim 1.02$), while non-thermal tails reached $\gamma \sim [1.02 - 2]$. The maximum energies observed ($\gamma \sim 10^2$) were noted as likely numerical artifacts of the 2D geometry and long confinement times, but the moderate-energy range is consistent with observations.
3. Acceleration Mechanisms: The dominant energization channel was identified as perpendicular gradient effects ($W_{\nabla B}$), which drive a second-order, Fermi-like stochastic acceleration. Curvature effects ($W_{\text{curv}}$) played a secondary role, primarily influencing particles on long trajectories or those exploring large domains. Magnetic mirror effects ($W_{\text{mirr}}$) were found to be strong in magnitude but acted only to exchange energy between parallel and perpendicular degrees of freedom, producing no net energization.
4. Trajectory Correlation: Statistical analysis revealed that the strongest acceleration occurs for electrons trapped in bouncing motions within turbulent magnetic structures. This is confirmed by a correlation between high energy gain ($\Delta \gamma$) and large pitch angles, indicating that particles trapped in magnetic mirrors experience the most efficient stochastic acceleration.
5. Temporal Evolution: The energization process was characterized by three phases: (I) an injection phase where thermal electrons adjust to the turbulent environment; (II) a sustained acceleration phase with efficient energy gain; and (III) a saturation phase where energy gain halts.

Significance and Claims
The paper claims to provide a coherent numerical framework for describing particle energetics in turbulent plasmas by correcting the energy conservation deficiencies in standard GCA implementations. By applying this improved formalism, the authors dissect the relative roles of magnetic effects, clarifying that perpendicular-gradient effects are the primary driver of stochastic, Fermi-like energization in post-flare loops.

The results support the interpretation that KHI-induced turbulence in coronal looptops is an efficient mechanism for producing and sustaining populations of trapped non-thermal electrons, thereby contributing to observed X-ray sources. The authors emphasize that while their 2.5D setup may overestimate confinement times and maximum energies compared to 3D reality, the findings offer a validated, reproducible test for future benchmarking of GCA implementations and a foundation for understanding magnetic-to-particle energy conversion in flare modeling. The study does not claim to resolve all kinetic feedback effects but establishes a robust method for quantifying acceleration in pre-existing turbulent states.