The Heating of Test Particles in Numerical Simulations of Alfvenic Turbulence

We study the heating of charged test particles in three-dimensional numerical simulations of weakly compressible magnetohydrodynamic (MHD) turbulence (Alfvenic turbulence''); these results are relevant to particle heating and acceleration in the solar wind, solar flares, accretion disks onto black holes, and other astrophysics and heliospheric environments. The physics of particle heating depends on whether the gyrofrequency of a particle is comparable to the frequency of a turbulent fluctuation that is resolved on the computational domain. Particles with these frequencies nearly equal undergo strong perpendicular heating (relative to the local magnetic field) and pitch angle scattering. By contrast, particles with large gyrofrequency undergo strong parallel heating. Simulations with a finite resistivity produce additional parallel heating due to parallel electric fields in small-scale current sheets. Many of our results are consistent with linear theory predictions for the particle heating produced by the Alfven and slow magnetosonic waves that make up Alfvenic turbulence. However, in contrast to linear theory predictions, energy exchange is not dominated by discrete resonances between particles and waves; instead, the resonances are substantially broadened.’’ We discuss the implications of our results for solar and astrophysics problems, in particular the thermodynamics of the near-Earth solar wind. We conclude that Alfvenic turbulence produces significant parallel heating via the interaction between particles and magnetic field compressions (``slow waves’’). However, on scales above the proton Larmor radius, Alfvenic turbulence does not produce significant perpendicular heating of protons or minor ions.

💡 Research Summary

This paper investigates how charged test particles are heated in three‑dimensional, weakly compressible magnetohydrodynamic (MHD) simulations that emulate Alfvénic turbulence. The authors generate a turbulent cascade containing both Alfvén and slow magnetosonic modes, then follow a large ensemble of non‑interacting test particles as they move under the Lorentz force produced by the simulated electric and magnetic fields. The central parameter is the particle gyrofrequency (Ωc) relative to the characteristic frequencies (ω) of the turbulent fluctuations that are resolved on the computational grid.

When Ωc is comparable to ω, particles experience strong perpendicular heating (relative to the local magnetic field) and rapid pitch‑angle scattering. This behavior matches the linear resonance picture in which particles exchange energy with Alfvénic or slow‑mode fluctuations whose phase speed along the field matches the particle’s parallel velocity. However, the simulations reveal that the resonances are not delta‑function‑like; instead they are substantially broadened. The broadening arises from nonlinear wave‑wave coupling, the stochastic nature of the turbulent field, and the presence of thin current sheets that locally modify the wave spectrum. Consequently, energy transfer occurs over a wide range of particle velocities rather than being confined to a narrow resonant shell.

In contrast, particles whose gyrofrequency far exceeds the turbulent frequencies (e.g., electrons or high‑energy ions) do not resonate with the Alfvénic fluctuations. Their heating is dominated by parallel electric fields associated with compressive slow‑mode structures. When a finite resistivity is introduced, the simulation develops small‑scale current sheets where strong parallel electric fields (E∥) appear. These fields provide an additional channel for parallel heating, especially for the high‑Ωc species. The result is a pronounced increase in the parallel temperature (T∥) relative to the perpendicular temperature (T⊥).

The authors compare their numerical findings with linear theory predictions for Alfvén and slow‑mode damping/heating rates. While the overall partition of heating between perpendicular (Alfvén‑driven) and parallel (slow‑mode‑driven) channels agrees with theory, the detailed spectral distribution of energy exchange does not. The simulations show that, on scales larger than the proton Larmor radius, Alfvénic turbulence does not produce significant perpendicular heating of protons or minor ions; instead, the dominant heating channel is parallel heating mediated by magnetic‑field compressions (slow waves).

Implications for heliophysics and astrophysics are discussed. In the near‑Earth solar wind, the observed T∥ > T⊥ anisotropy of ions can be explained by the parallel heating identified here, without invoking kinetic Alfvén wave (KAW) dissipation at sub‑ion scales. In solar flares, the presence of resistive current sheets may provide the parallel electric fields needed for rapid electron energization. For accretion disks around black holes, where turbulence is expected to be strong and compressive, the same mechanisms could contribute to the heating of electrons that ultimately radiate the observed X‑ray emission.

In summary, the study demonstrates that Alfvénic turbulence heats particles through two distinct pathways that depend on the ratio Ωc/ω: (1) near‑resonant particles receive perpendicular heating and pitch‑angle scattering, with broadened resonances due to turbulence; (2) high‑frequency particles receive parallel heating from slow‑mode compressions and, when resistivity is present, from parallel electric fields in current sheets. These results refine our understanding of plasma heating in space and astrophysical environments and highlight the importance of including both wave‑particle resonance broadening and small‑scale dissipative structures in kinetic models of turbulent plasmas.