Magnetic Field-Line Curvature and Its Role in Particle Acceleration by Magnetically Dominated Turbulence

TL;DR

The study addresses how magnetic field-line curvature in magnetically dominated turbulence facilitates particle acceleration via curvature-drift along the motional electric field. Using fully kinetic PIC simulations across varying fluctuation-to-mean-field ratios and a guiding-center framework, the authors characterize curvature statistics, intermittency, and the energy-transfer channels, notably the $E_\perp$-driven energization. They find broad curvature PDFs with heavy tails and symmetric contraction statistics best described by a symmetric Pareto distribution, with intermittency and tail strength increasing as the mean field weakens. Curvature-drift energization emerges as a dominant mechanism for transferring energy from turbulence to nonthermal particles across a range of $\delta B_0/B_0$, offering a kinetic basis for turbulent particle acceleration in astrophysical plasmas and a pathway to connect field-line geometry with observable high-energy phenomena.

Abstract

We employ first-principles, fully kinetic particle-in-cell simulations to investigate magnetic field-line curvature in magnetically dominated turbulent plasmas and its role in particle acceleration through curvature-drift motion along the motional electric field. By varying the fluctuation-to-mean magnetic-field ratio $δB_0/B_0$, we examine curvature $κ$ statistics and their connection to particle acceleration. The curvature probability densities display broad power-law wings, scaling linearly in $κ$ below the peak and developing hard high-$κ$ tails for $δB_0/B_0 \gtrsim 1$. As the mean field strengthens, the high-$κ$ tails steepen, and large-curvature events are suppressed when $δB_0/B_0 \ll 1$. The probability density functions of magnetic field-line contraction, ${\bf v}_E \cdot {\bf κ}$, with ${\bf v}_E$ the field-line velocity, develop power-law tails well described by a symmetric Pareto distribution, characteristic of stochastic energy exchanges, with the tails becoming harder as $δB_0/B_0$ increases. Our guiding-center analysis shows that curvature-drift acceleration accounts for a substantial fraction of the energization via the motional electric field, and that it strengthens with increasing $δB_0/B_0$. For well-magnetized particles, curvature-drift acceleration typically exceeds ${\bf\nabla}B$ drift, polarization drift, and betatron contributions. These results identify curvature-drift acceleration as a principal pathway through which magnetized turbulence transfers energy to nonthermal particles in astrophysical plasmas.

Figures (10)

• Figure 1: Perpendicular and parallel spectra of (a) magnetic and (b) electric fluctuations, from the simulation with $\delta B_0/B_0=1$ at $t=2\,l_c/c$. Power-law slopes of $k_\perp^{-5/3}$ and $k_\parallel^{-2}$ are shown for reference. Here, perpendicular ($\perp$) and parallel ($\parallel$) are measured with respect to ${\bm{B}}_0$.
• Figure 2: Second-order structure function of magnetic fluctuations (a) and anisotropy of magnetic and electric fluctuations (b), from the simulation with $\delta B_0/B_0=1$ at $t=2\,l_c/c$. Here, perpendicular ($\perp$) and parallel ($\parallel$) are measured with respect to the scale-dependent local mean magnetic field, as detailed in CV2000.
• Figure 3: Probability density functions of field increments at different spatial separations, from the outer scale $r = 42d_i$ down to the ion inertial scale $r = d_i$, including intermediate scales. Results are from the reference simulation ($\delta B_0/B_0=1$) at $t=2\,l_c/c$. The increments $\Delta f(r)$ are (a) $\Delta B_x(r)$ and (b) $\Delta E_x(r)$, each normalized by their standard deviation at the corresponding $r$. At small separations, the PDFs show pronounced non-Gaussian heavy tails, while at large $r$ they approach a Gaussian distribution (dashed curve).
• Figure 4: Two-dimensional slices in the $x$–$y$ plane at fixed $z$ from the reference simulation ($\delta B_0/B_0=1$) at $t=2\,l_c/c$, showing the normalized magnetic field strength $B/B_{\rm rms}$ and the normalized magnetic field-line curvature magnitude $\kappa/\kappa_{\rm rms}$. The bottom panel shows the joint distribution of $B/B_{\rm rms}$ and $\kappa$ (measured in units of the inverse coherence length, $l_c^{-1}$), showing that regions of strong curvature tend to coincide with weaker magnetic fields. A power-law slope of $-1/2$ is provided for reference.
• Figure 5: (a) Probability density functions of the normalized magnetic field-line curvature, $\kappa l_c$, from simulations with different $\delta B_0/B_0$. Power-law slopes $\kappa^{+1}$ and $\kappa^{-2.5}$ are shown for reference. (b) Probability density functions of the normalized magnetic field-line contraction, $({\bm{v}}_E \cdot {\bm{\kappa}})l_c/c$, for the same set of simulations, with symmetric Pareto distribution fits plotted as dashed lines.