  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.
