Magnetised turbulent plasmas as high-energy particle accelerators

This proceedings paper reports on the theoretical modelling of particle acceleration in magnetised turbulent plasmas. It briefly reviews some recent findings obtained from fully kinetic numerical simulations of large-amplitude, semi to fully relativistic turbulence. The paper then argues that these findings can be understood within the framework of a ``generalised Fermi’’ picture of stochastic acceleration, which it summarises. The dominant contributions to acceleration appear to arise from particle interactions with sharp, dynamic bends of the magnetic field lines and regions of velocity compression. Interestingly, the acceleration rate is spatially inhomogeneous and its probability distribution follows a broken power law extending up to large values. This makes relativistic, large-amplitude turbulence an extreme particle accelerator. Some implications for particle transport and the shape of the particle energy spectrum in the presence of radiative losses and over long timescales are also discussed.

💡 Research Summary

This proceedings paper presents a comprehensive study of particle acceleration in magnetised turbulent plasmas, focusing on the regime of large‑amplitude (δB_rms ≈ B₀) and relativistic eddy velocities (δv_rms ∼ c). Using state‑of‑the‑art fully kinetic particle‑in‑cell (PIC) simulations, the author explores a parameter space characterized by a magnetisation σ ≈ 10, which corresponds to relativistic Alfvénic turbulence. The simulations inject test particles into a three‑dimensional turbulent box and follow their evolution without radiative or escape losses. The resulting particle energy distribution exhibits a clear thermal peak at ε_th ≈ σ m c², followed by a non‑thermal power‑law tail dN/dε ∝ ε⁻ˢ that extends over several decades in energy up to a cut‑off ε_c. The spectral index s varies from ≈4 in sub‑relativistic turbulence to ≈2 in the highly relativistic limit, in agreement with earlier kinetic studies but in stark contrast with the log‑normal spectra predicted by a purely diffusive Fokker‑Planck description using the measured diffusion coefficient D_εε ≈ 0.1 σ ε² c/ℓ_c.

A key finding is that the acceleration rate is not spatially uniform. By measuring the local acceleration rate α = dε/dt for individual particles, the probability density function f(α) is shown to follow a broken power law: it peaks near the mean value but possesses a high‑α tail that extends over orders of magnitude. This inhomogeneity explains why a single mean diffusion coefficient fails to reproduce the simulated spectra; a subset of particles experiences transient, very rapid energisation while most particles evolve at the average rate.

To interpret these results, the paper introduces a “generalised Fermi” framework. Instead of treating turbulence as a collection of discrete scattering centres (the classic Fermi pinball), the author follows particles through a continuous sequence of instantaneous local rest frames R/E defined by the local plasma velocity v_E. In each R/E frame the electric field vanishes, and particle energy changes arise solely from inertial forces associated with spatial and temporal gradients of v_E. This formalism naturally incorporates acceleration mechanisms that have no analogue in wave‑particle resonance theory, notably curvature‑drift‑like acceleration when particles encounter sharp, dynamic bends of magnetic field lines, and first‑order energisation in compressive regions. The resulting transport equation contains both an energy‑diffusion term and an energy‑advection term; the latter can become negative over certain energy ranges, reflecting the observed “super‑diffusive” bursts of acceleration.

The paper further discusses astrophysical implications. The broken‑power‑law distribution of acceleration rates implies that, even in the presence of strong radiative cooling (synchrotron, inverse‑Compton) or particle escape, the high‑energy tail can remain hard, because a fraction of particles continuously samples the high‑α tail. The cut‑off energy ε_c is set by a Hillas‑type condition r_g(ε_c) ≈ ℓ_c, linking the maximum attainable energy to the outer scale of the turbulence rather than microscopic plasma parameters. Consequently, relativistic, large‑amplitude turbulence emerges as an “extreme accelerator” capable of producing ultra‑high‑energy cosmic rays, PeV neutrinos, and TeV–PeV photons in environments such as black‑hole coronae, relativistic jets, and galaxy‑cluster shocks.

In summary, the work bridges kinetic simulation results with an analytically tractable, covariant stochastic acceleration theory, highlighting the dominant role of magnetic‑field curvature and compressive motions, the crucial importance of spatially intermittent acceleration, and the resulting hard, extended particle spectra that are relevant for multi‑messenger astrophysics.