Early-time velocity autocorrelation for charged particles diffusion and drift in static magnetic turbulence

Using test-particle simulations, we investigate the temporal dependence of the two-point velocity correlation function for charged particles scattering in a time-independent spatially fluctuating magnetic field derived from a three-dimensional isotropic turbulence power spectrum. Such a correlation function allowed us to compute the spatial coefficients of diffusion both parallel and perpendicular to the average magnetic field. Our simulations confirm the dependence of the perpendicular diffusion coefficient on turbulence energy density and particle energy predicted previously by a model for early-time charged particle transport. Using the computed diffusion coefficients, we exploit the particle velocity autocorrelation to investigate the time-scale over which the particles “decorrelate” from the solution to the unperturbed equation of motion. Decorrelation time-scales are evaluated for parallel and perpendicular motions, including the drift of the particles from the local magnetic field line. The regimes of strong and weak magnetic turbulence are compared for various values of the ratio of the particle gyroradius to the correlation length of the magnetic turbulence. Our simulation parameters can be applied to energetic particles in the interplanetary space, cosmic rays at the supernova shocks, and cosmic-rays transport in the intergalactic medium.

💡 Research Summary

The authors investigate the early‑time dynamics of charged particles moving in a static, three‑dimensional magnetic turbulence field. Using test‑particle simulations, they generate a spatially fluctuating magnetic field from an isotropic power‑law spectrum (Kolmogorov‑type, (P(k)\propto k^{-5/3})) and follow ensembles of particles with a range of gyroradii (r_g) relative to the turbulence correlation length (L_c). The key diagnostic is the two‑point velocity autocorrelation function (R_{ij}(t)=\langle v_i(0)v_j(t)\rangle), which they compute directly for parallel ((i=j=\parallel)) and perpendicular ((i=j=\perp)) components. By integrating (R_{ij}(t)) over time they obtain the parallel and perpendicular diffusion coefficients (\kappa_{\parallel}) and (\kappa_{\perp}) via the Taylor‑Green‑Kubo relation.

The simulations cover a broad parameter space: turbulence strength (\delta B/B_0) ranging from weak ((\ll1)) to strong ((\gtrsim1)) and the dimensionless rigidity (\rho=r_g/L_c) from 0.01 to 10. This allows the authors to emulate conditions relevant to solar‑wind energetic particles, cosmic rays accelerated at supernova remnant shocks, and cosmic‑ray transport in the intergalactic medium.

Results show that the velocity autocorrelation decays rapidly within a few gyro‑periods. In the weak‑turbulence regime the decay is well described by exponential laws, (R_{\parallel}(t)\sim \exp(-\nu t)) and (R_{\perp}(t)\sim \exp(-\nu_{\perp} t)\cos(\Omega t)), with decay rates (\nu,\nu_{\perp}) proportional to (\delta B/B_0) and to the rigidity (\rho). In strong turbulence the decay becomes non‑exponential, reflecting the breakdown of the Markovian assumption underlying many analytic diffusion models.

Integrating the autocorrelation yields diffusion coefficients that reproduce previously predicted scalings: (\kappa_{\perp}\propto (\delta B/B_0)^2 (r_g/L_c)^2). Importantly, the authors quantify the “drift” of particles away from their local magnetic field line—a process that dominates (\kappa_{\perp}) when (\rho\ll1). This drift manifests as a gradual decorrelation from the unperturbed helical orbit, leading to a finite perpendicular step size even before the particle’s motion becomes fully stochastic.

To characterize the decorrelation process, the authors define parallel and perpendicular decorrelation times, (\tau^{\parallel}{\rm dec}) and (\tau^{\perp}{\rm dec}). They find that (\tau^{\parallel}{\rm dec}) decreases roughly linearly with increasing turbulence amplitude, whereas (\tau^{\perp}{\rm dec}) exhibits a strong, non‑linear dependence on rigidity, dropping sharply near (\rho\approx1). These timescales provide a physical measure of how long particles retain memory of the solution to the unperturbed equation of motion before the turbulent field dominates their dynamics.

Comparing weak and strong turbulence regimes, the study reveals that in weak turbulence particles remain tightly bound to the mean field line, and perpendicular transport is governed by small‑amplitude, cumulative drifts. In strong turbulence the magnetic geometry is heavily distorted, leading to rapid loss of correlation, large (\kappa_{\perp}), and short decorrelation times—effectively a transition to a quasi‑ballistic perpendicular transport regime.

The authors discuss the astrophysical implications of their findings. For solar‑wind energetic particles, the derived (\kappa_{\perp}) and (\tau^{\perp}{\rm dec}) can improve models of longitudinal spread and cross‑field diffusion observed in SEP events. For supernova‑remnant shocks, the rigidity dependence of (\kappa{\perp}) informs the escape of high‑energy cosmic rays from the acceleration region, influencing the observed spectrum. In the intergalactic medium, where turbulence is expected to be weak and correlation lengths large, the drift‑dominated perpendicular diffusion may dominate cosmic‑ray propagation over cosmological distances.

In summary, by directly measuring the early‑time velocity autocorrelation in static magnetic turbulence, the paper provides a robust, simulation‑based validation of theoretical diffusion scalings, introduces quantitative decorrelation times for both parallel and perpendicular motions, and highlights the role of drift from local field lines in shaping perpendicular transport. These results bridge the gap between idealized diffusion theory and realistic particle transport in a variety of astrophysical plasma environments.