Particle Acceleration in Magnetized Shear-Driven Turbulence

TL;DR

The paper addresses how particles gain energy in subsonic, magnetized shear-driven turbulence excited by the Kelvin-Helmholtz instability while including full particle backreaction. It uses 2D MHD-PIC simulations with a flow-aligned magnetic field and continuous external driving to generate a quasi-steady turbulent state, modeling the energetics of a large particle population. The main finding is that sustained particle acceleration requires driving, and the mechanism is a second-order Fermi process arising from systematic distortion of gyro-orbits by turbulent motional electric fields, yielding $\ig<\Delta \epsilon\big> \propto U_{\text{shear}}^2$ and $D_{\gamma} \propto \gamma^2$, along with non-thermal tails and log-normal energy distributions for crossers; high-energy particles show pitch-angle anisotropy with $\mu \to 0$ as $r_g$ grows beyond the turbulence scale. The study provides a general, backreaction-inclusive model for cosmic-ray energization in shear flows applicable to AGN jets, ICM turbulence, and accretion-flow boundaries, and suggests extensions to 3D and relativistic regimes.

Abstract

Shear flows, ubiquitous in space and astrophysical plasmas, can accelerate particles through turbulence excited by the Kelvin-Helmholtz instability. We present the first numerical study of particle acceleration in non-relativistic, magnetized, and purely shear-driven turbulence that includes full particle backreaction. Using two-dimensional MHD-PIC simulations with an initially uniform flow-aligned magnetic field and external stirring force, we demonstrate that sustained particle acceleration requires continuously driven turbulence, whereas freely decaying turbulence rapidly depletes its energy reservoirs and halts the acceleration. The acceleration mechanism operates through the systematic distortion of gyro-orbits by turbulent electric fields: acceleration phases extend the particle trajectory along the electric force, increasing the energy gain, while deceleration phases shorten the trajectory, reducing the energy loss. This asymmetry produces net energy gain despite stochastic fluctuations, with the mean energy change scaling quadratically with shear velocity, characteristic of second-order Fermi acceleration. Initially monoenergetic particles develop substantial non-thermal tails after the turbulence onset. For particles repeatedly crossing shear layers, their energization follows geometric Brownian motion with weak systematic drift, yielding a log-normal distribution. High-energy particles exhibit pitch-angle anisotropy, becoming preferentially perpendicular to the flow-aligned magnetic field as their gyroradii exceed the turbulent layer width. These results establish shear-driven turbulence as a viable particle acceleration mechanism, providing a general model for particle energization in shear flows.

Figures (17)

• Figure 1: (a) Initial shear profile of background plasma's velocity field, as a function of y positions. (b) $X$-averaged mean flow profile at different times, ranging from $\mathrm{a}/\mathrm{U_0}$ = 0 to 100. Both profiles are normalized by the characteristic flow velocity $U_0$. Fluctuations near the velocity transition region are caused by turbulent eddies.
• Figure 2: Evolution of the Kelvin-Helmholtz instability showing snapshots at $a/U_0$ = 40 (top) and 80 (bottom). Left column: background fluid density $\rho$. Middle column: velocity field $\boldsymbol{u}$ with streamlines indicating flow direction. Right column: magnetic field $\boldsymbol{B}$ with field lines shown. Colors represent their respective normalized magnitudes.
• Figure 3: Top: Energy power spectrum of shear-driven turbulence in steady state as a function of wavenumber $k$. Red and blue curves show kinetic ($E_k$) and magnetic ($E_m$) energy densities, respectively; black curve combines both components ($E_k + E_m$). Dashed line indicates $k^{-2}$ scaling in the inertial range. Bottom: Magnetic-to-kinetic energy ratio $E_m/E_k$. Gray vertical line marks the energy injection scale $\lambda_{\rm inj}$ in both panels.
• Figure 4: Time evolution of volume-averaged energy densities in driven (solid) and freely decaying (dashed) turbulence. Top: Background fluid kinetic ($\varepsilon_k$, red) and magnetic ($\varepsilon_m$, blue) energy densities vs. time $a/U_0$. Bottom: Particle energy density $\varepsilon_p$ vs. time.
• Figure 5: Evolution of the particle energy distribution $dN/d\gamma$ from $a/U_0 = 0$ to $2400$ (color-coded by time). The initially monoenergetic distribution develops a substantial non-thermal tail, with the high-energy cutoff advancing continuously through sustained acceleration.