Charged particle dynamics in the presence of non-Gaussian L\'evy electrostatic fluctuations

Full orbit dynamics of charged particles in a $3$-dimensional helical magnetic field in the presence of $\alpha$-stable L\'evy electrostatic fluctuations and linear friction modeling collisional Coulomb drag is studied via Monte Carlo numerical simulations. The L\'evy fluctuations are introduced to model the effect of non-local transport due to fractional diffusion in velocity space resulting from intermittent electrostatic turbulence. The probability distribution functions of energy, particle displacements, and Larmor radii are computed and showed to exhibit a transition from exponential decay, in the case of Gaussian fluctuations, to power law decay in the case of L\'evy fluctuations. The absolute value of the power law decay exponents are linearly proportional to the L\'evy index $\alpha$. The observed anomalous non-Gaussian statistics of the particles' Larmor radii (resulting from outlier transport events) indicate that, when electrostatic turbulent fluctuations exhibit non-Gaussian L\'evy statistics, gyro-averaging and guiding centre approximations might face limitations and full particle orbit effects should be taken into account.