1/f noise and anomalous scaling in L´evy on-off intermittency 3
S(f) = ReiftC(t)dt, according to the Wiener-Khinchin theorem [40, 41]. This has been
exploited to explain the S(f)∝f−1/2range observed at low frequencies for small µ > 0,
cf. [42, 43]. Such behaviour, namely the existence of a wide range in log(f), at small f,
for which the PSD S(f) is of power-law form with exponent smaller than 0 and greater
than −2, is generically called 1/f noise, also known as Flicker noise, or pink noise.
It has been observed in a wide variety of systems, ranging from voltage and current
fluctuations in vacuum tubes and transistors, where this behaviour was first recognised
[44, 45, 46], to astrophysical magnetic fields [47] and biological systems [48], climate
[49], turbulent flows [50, 51, 52], reversing flows [53, 54, 55, 56], traffic [57], as well as
music and speech [58, 59], to name a few, and is also found in fractional renewal models
[60]. In addition, 1/f noise has also been observed for L´evy flights in inhomogeneous
environments [61, 62], but these studies did not consider any bifurcation points.
While the above-described case of Gaussian noise has been studied in depth, non-
Gaussian fluctuations arise in many systems. For example, out-of-equilibrium dynamics,
such as turbulent fluid flows, typically exhibit non-Gaussian statistics, see e.g. [63],
implying that instabilities developing on a turbulent background generally exhibit non-
Gaussian growth rates, cf. [19, 64]. Power-law-distributed fluctuations in particular
are found in a variety of systems, including the human brain [65], climate [66], finance
[67] and beyond. An important example of random motion resulting from additive
non-Gaussian noise is given by L´evy flights (a term coined by Mandelbrot [68]), which
are driven by L´evy noise. L´evy noise follows a heavy-tailed α-stable distribution that
depends on a stability parameter α∈(0,2] and a skewness parameter β∈[−1,1].
Stable distributions come in different forms: the case α= 2 corresponds to the Gaussian
distribution, while at α < 2 the distribution has power-law tails with exponent −1−α.
The main interest lies in the parameter regime 1 < α ≤2, where there is a finite mean,
but an infinite variance. While the parameter regime 0 < α ≤1 is formally admissible,
it is of little practical interest, since the noise distribution has a diverging mean in
this case. The reason why Gaussian random variables are common in physics is their
stability: by the central limit theorem [69], the Gaussian distribution constitutes an
attractor in the space of PDFs with finite variance. Similarly, by the generalised central
limit theorem [70, 71], non-Gaussian α-stable distributions constitute an attractor in
the space of PDFs whose variance does not exist. Stable distributions can be symmetric
(β= 0) or asymmetric (β6= 0), giving rise to symmetric and asymmetric L´evy flights.
L´evy flights have since found numerous applications in many areas both in physics
[72, 73, 74, 75, 76, 77] and beyond, including climatology [78], finance [79], ecology [80]
and human travel [81].
L´evy statistics and on-off intermittency can be present in the same system.
Examples include experiments of human balancing motion [12, 82, 83], blinking quantum
dots [84, 85] and the intermittent growth of three-dimensional instabilities in quasi-
two-dimensional turbulence [19]. In a recent study [86], the problem of L´evy on-off
intermittency was formally introduced as the case where f(t) in equation (1) is L´evy
noise with 1 < α < 2. In this case, if X(t) solves (1), then log(X(t)) performs a L´evy