Lévy flight

Lévy flights are, by construction, Markov processes. For general distributions of the step-size, satisfying the power-like condition, the distance from the origin of the random walk tends, after a large number of steps, to a stable distribution due to the generalized central limit theorem, enabling many processes to be modeled using Lévy flights.

The probability densities for particles undergoing a Levy flight can be modeled using a generalized version of the Fokker–Planck equation, which is usually used to model Brownian motion. The equation requires the use of fractional derivatives. For jump lengths which have a symmetric probability distribution, the equation takes a simple form in terms of the Riesz fractional derivative. In one dimension, the equation reads as

${\displaystyle {\frac {\partial \varphi (x,t)}{\partial t}}=-{\frac {\partial }{\partial x}}f(x,t)\varphi (x,t)+\gamma {\frac {\partial ^{\alpha }\varphi (x,t)}{\partial |x|^{\alpha }}}}$

where γ is a constant akin to the diffusion constant, α is the stability parameter and f(x,t) is the potential. The Riesz derivative can be understood in terms of its Fourier Transform.

${\displaystyle F_{k}\left[{\frac {\partial ^{\alpha }\varphi (x,t)}{\partial |x|^{\alpha }}}\right]=-|k|^{\alpha }F_{k}[\varphi (x,t)]}$

This can be easily extended to multiple dimensions.

Another important property of the Lévy flight is that of diverging variances in all cases except that of α = 2, i.e. Brownian motion. In general, the θ fractional moment of the distribution diverges if α ≤ θ. Also,

${\displaystyle \left\langle |x|^{\theta }\right\rangle \propto t^{\theta /\alpha }\quad {\text{if }}\theta <\alpha .}$

The exponential scaling of the step lengths gives Lévy flights a scale invariant property, and they are used to model data that exhibits clustering.

Figure 1. An example of 1000 steps of a Lévy flight in two dimensions. The origin of the motion is at [0,0], the angular direction is uniformly distributed and the step size is distributed according to a Lévy (i.e. stable) distribution with α = 1 and β = 0 which is a Cauchy distribution. Note the presence of large jumps in location compared to the Brownian motion illustrated in Figure 2.

Figure 2. An example of 1000 steps of an approximation to a Brownian motion type of Lévy flight in two dimensions. The origin of the motion is at [0, 0], the angular direction is uniformly distributed and the step size is distributed according to a Lévy (i.e. stable) distribution with α = 2 and β = 0 (i.e., a normal distribution).

The definition of a Lévy flight stems from the mathematics related to chaos theory and is useful in stochastic measurement and simulations for random or pseudo-random natural phenomena. Examples include earthquake data analysis, financial mathematics, cryptography, signals analysis as well as many applications in astronomy, biology, and physics.

It has been found that jumping between climate states observed in the paleoclimatic record can be described as a Lévy flight or an alpha-stable process[7] Another application is the Lévy flight foraging hypothesis. When sharks and other ocean predators cannot find food, they abandon Brownian motion, the random motion seen in swirling gas molecules, for Lévy flight — a mix of long trajectories and short, random movements found in turbulent fluids. Researchers analyzed over 12 million movements recorded over 5,700 days in 55 data-logger-tagged animals from 14 ocean predator species in the Atlantic and Pacific Oceans, including silky sharks, yellowfin tuna, blue marlin and swordfish. The data showed that Lévy flights interspersed with Brownian motion can describe the animals' hunting patterns.[8][9][10][11] Birds and other animals (including humans)[12] follow paths that have been modeled using Lévy flight (e.g. when searching for food).[13] Biological flight data can also apparently be mimicked by other models such as composite correlated random walks, which grow across scales to converge on optimal Lévy walks.[13] Composite Brownian walks can be finely tuned to theoretically optimal Lévy walks but they are not as efficient as Lévy search across most landscapes types, suggesting selection pressure for Lévy walk characteristics is more likely than multi-scaled normal diffusive patterns.[14]

Efficient routing in a network can be performed by links having a Levy flight length distribution with specific values of alpha.[2]

• Anomalous diffusion
• Fat-tailed distribution
• Heavy-tailed distribution
• Lévy process
• Lévy alpha-stable distribution
• Lévy flight foraging hypothesis

• Mandelbrot, Benoit B. (1982). The Fractal Geometry of Nature (Updated and augm. ed.). New York: W. H. Freeman. ISBN 0-7167-1186-9. OCLC 7876824.

• Cheng, Z.; Savit, R. (1987). "Fractal and nonfractal behavior in Levy flights" (PDF). Journal of Mathematical Physics. 28 (3): 592. Bibcode:1987JMP....28..592C. doi:10.1063/1.527644. hdl:2027.42/70735.
• Shlesinger, Michael F.; Klafter, Joseph; Zumofen, Gert (December 1999). "Above, below and beyond Brownian motion" (PDF). American Journal of Physics. 67 (12): 1253–1259. Bibcode:1999AmJPh..67.1253S. doi:10.1119/1.19112. Archived from the original (PDF) on 2012-03-28.

• A comparison of the paintings of Jackson Pollock to a Lévy flight model