Lévy Flights and Leaky Boxes: Anomalous Diffusion of Cosmic Rays

Abstract

In classical diffusion, particle step-sizes have a Gaussian distribution. However, in superdiffusion, they have power-law tails, with transport dominated by rare, long Lévy flights. Similarly, if the time interval between scattering events has power-law tails, subdiffusion occurs. Both forms of anomalous diffusion are seen in cosmic ray (CR) particle tracking simulations in turbulent magnetic fields. They also likely occur if CRs are scattered by discrete intermittent structures. Anomalous diffusion mimics a scale-dependent diffusion coefficient, with potentially wide-ranging consequences. However, the finite size of galaxies implies an upper bound on step-sizes before CRs escape. This truncation results in eventual convergence to Gaussian statistics by the central limit theorem. Using Monte-Carlo simulations, we show that this occurs in both standard finite-thickness halo models, or when CR diffusion transitions to advection or streaming-dominated regimes. While optically thick intermittent structures produce power-law trapping times and thus subdiffusion, gaussianization also eventually occurs on timescales longer than the maximum trapping time. Anomalous diffusion is a transient, and converges to standard diffusion on the (usually short) timescale of particle escape, either from confining structures (subdiffusion), or the system as a whole (superdiffusion). Thus, standard assumptions of classical diffusion are physically justified in most applications, despite growing simulation evidence for anomalous diffusion. However, if escape times are long, this is no longer true. For instance, anomalous diffusion in the CGM or ICM would change CR pressure profiles. Finally, we show the standard diagnostic for anomalous diffusion, $⟨d^2⟩\propto t^{\alpha }$ with $\alpha \ne 1$, is not justified for truncated Lévy flights, and propose an alternative robust measure.

References in corpus (13)

• Lévy walks
• Extended gamma-ray sources around pulsars constrain the origin of the positron flux at Earth
• Matrix approach to discrete fractional calculus II: partial fractional differential equations
• Cosmic Ray Propagation: Nonlinear Diffusion Parallel and Perpendicular to Mean Magnetic Field
• Kinetic turbulence in relativistic plasma: from thermal bath to non-thermal continuum
• Cosmic ray-driven galactic winds: streaming or diffusion?
• First-passage and first-hitting times of Levy flights and Levy walks
• The optical depth of the Universe to ultrahigh energy cosmic ray scattering in the magnetized large scale structure
• A Random Walk to a Non-Ergodic Equilibrium Concept
• Fractional Transport in Strongly Turbulent Plasmas
• Particle Acceleration and Fractional Transport in Turbulent Reconnection
• Test of the superdiffusion model in the interstellar medium around the Geminga pulsar
• Cumulant Approach of Arbitrary Truncated Levy Flight

Cited by in corpus (2)

• Cosmic-ray transport in inhomogeneous media
• Modelling cosmic-ray transport: magnetised versus unmagnetised motion in astrophysical magnetic turbulence