Asymptotic densities of ballistic Lévy walks

Biased Lévy Walk Enables Light Gradient Sensing in Euglena gracilis

We examine the navigation behavior of the photosensitive alga Euglena gracilis in confined environments. Under uniform lighting conditions, E. gracilis exhibits stochastic movements with nearly straight runs interrupted by abrupt directional changes. The lengths of these runs follow a long-tailed distribution typical of a Lévy walk, with scaling exponents that vary with light intensity. In gradient lighting conditions, the cells modulate their run durations-extending them upon detecting an increase in light intensity and shortening them when a decrease is detected. This adjustment effectively biases the Lévy walk, enabling the cells to ascend the spatial light gradient. This behavior mirrors well-known prokaryotic stochastic navigation strategies, such as bacterial chemotaxis, offering a eukaryotic parallel. The experimental observations under varied lighting conditions are consistently replicated through an agent-based model.

Generalized two-state random walk model: Nontrivial anomalous diffusion, aging, and ergodicity breaking

The intermittent stochastic motion is a dichotomous process that alternates between two distinct states. This phenomenon, observed across various physical and biological systems, is attracting increasing interest and highlighting the need for comprehensive theories to describe it. In this paper, we introduce a generalized intermittent random walk model based on a renewal process that alternates between the continuous time random walk (CTRW) state and the generalized Lévy walk (gLW) state. Notably, the nonlinear space-time coupling inherent in the gLW state allows this generalized model to encompass a variety of random walk models and makes it applicable to diverse systems. By deriving the velocity correlation function and utilizing the scaling Green-Kubo relation, the ensemble-averaged and time-averaged mean-squared displacement (MSD) is calculated, and the anomalous diffusive behavior, aging effect, and ergodic property of the model are further analyzed and discussed. The results reveal that, due to the intermittent nature, there are two diffusive terms in the expression of the MSD, and the diffusion can be intermediately characterized by the diffusive term with the largest diffusion coefficient instead of the diffusive term with the largest diffusion exponent, which is significantly different from single-state stochastic process. We demonstrate that, due to the power-law distribution of sojourn times, nonlinear space-time coupling, and intermittent characteristics, both ergodicity and nonergodicity can coexist in intermittent stochastic processes.

Langevin picture of anomalous diffusion processes in expanding medium

The expanding medium is very common in many different fields, such as biology and cosmology. It brings a nonnegligible influence on particle's diffusion, which is quite different from the effect of an external force field. The dynamic mechanism of a particle's motion in an expanding medium has only been investigated in the framework of a continuous-time random walk. To focus on more diffusion processes and physical observables, we build the Langevin picture of anomalous diffusion in an expanding medium, and conduct detailed analyses in the framework of the Langevin equation. With the help of a subordinator, both subdiffusion process and superdiffusion process in the expanding medium are discussed. We find that the expanding medium with different changing rate (exponential form and power-law form) leads to quite different diffusion phenomena. The particle's intrinsic diffusion behavior also plays an important role. Our detailed theoretical analyses and simulations present a panoramic view of investigating anomalous diffusion in an expanding medium under the framework of the Langevin equation.

Strong anomalous diffusive behaviors of the two-state random walk process

The phenomenon of the two-state process is observed in various systems and is increasingly attracting attention, such that there is a need for a theoretical model of the process. In this paper, we present a prototypal two-state random walk (TSRW) model of a renewal process alternating between the continuous-time random walk (CTRW) state and Lévy walk (LW) state. The jump length distribution of the CTRW state is assumed to be Gaussian whereas the time distributions of the two states are both considered to follow a power law. The diffusive behavior is analyzed and discussed by calculating the mean squared displacement (MSD) analytically and numerically. The results reveal that it displays strong anomalous diffusive behaviors caused by random motions of both states, i.e., two anomalous diffusion terms coexist in the expression of the MSD, and the time distribution which has the heavier tail determines their forms. Moreover, because the two diffusion terms originate from different mechanisms, we find that the diffusion can be characterized by either the term with the largest diffusion exponent or the term with the largest diffusion coefficient at long timescales, which shows very different properties from the single-state process. In addition, the two-state nature of the process of the particle moving in a velocity field makes the TSRW model applicable to describe it. Results obtained from the two-state model reveal that the diffusion can even exhibit subdiffusive behavior, which is significantly different from known results obtained using the single-state model.

Nonergodicity of d-dimensional generalized Lévy walks and their relation to other space-time coupled models

We investigate the nonergodicity of the generalized Lévy walk introduced by Shlesinger et al. [Phys. Rev. Lett. 58, 1100 (1987)PRLTAO0031-900710.1103/PhysRevLett.58.1100] with respect to the squared displacements. We present detailed analytical derivations of our previous findings outlined in a recent letter [Phys. Rev. Lett. 120, 104501 (2018)PRLTAO0031-900710.1103/PhysRevLett.120.104501], give detailed interpretations, and in particular emphasize three surprising results. First, we find that the mean-squared displacements can diverge for a certain range of parameter values. Second, we show that an ensemble of trajectories can spread subdiffusively, whereas individual time-averaged squared displacements show superdiffusion. Third, we recognize that the fluctuations of the time-averaged squared displacements can become so large that the ergodicity breaking parameter diverges, what we call infinitely strong ergodicity breaking. This phenomenon can also occur for paramter values where the lag-time dependence of the mean-squared displacements is linear indicating normal diffusion. In order to numerically determine the full distribution of time-averaged squared displacements, we use importance sampling. For an embedding of our findings into existing results in the literature, we define a more general model which we call variable speed generalized Lévy walk and which includes well-known models from the literature as special cases such as the space-time coupled Lévy flight or the anomalous Drude model. We discuss and interpret our findings regarding the generalized Lévy walk in detail and compare them with the nonergodicity of the other space-time coupled models following from the more general model.

Lévy walks with rests: Long-time analysis

In this paper we analyze the asymptotic behavior of Lévy walks with rests. Applying recent results in the field of functional convergence of continuous-time random walks we find the corresponding limiting processes. Depending on the parameters of the model, we show that in the limit we can obtain standard Lévy walk or the process describing competition between subdiffusion and Lévy flights. Some other more complicated limit forms are also possible to obtain. Finally we present some numerical results, which confirm our findings.

Gaussian process and Lévy walk under stochastic noninstantaneous resetting and stochastic rest

A stochastic process with movement, return, and rest phases is considered in this paper. For the movement phase, the particles move following the dynamics of the Gaussian process or the ballistic type of Lévy walk, and the time of each movement is random. For the return phase, the particles will move back to the origin with a constant velocity or acceleration or under the action of a harmonic force after each movement, so that this phase can also be treated as a noninstantaneous resetting. After each return, a rest with a random time at the origin follows. The asymptotic behaviors of the mean-squared displacements with different kinds of movement dynamics, returning, and random resting time are discussed. The stationary distributions are also considered when the process is localized. In addition, the mean first passage time is considered when the dynamic of the movement phase is Brownian motion.

Numerical Solution to Anomalous Diffusion Equations for Levy Walks

The process of Levy random walks is considered in view of the constant velocity of a particle. A kinetic equation is obtained that describes the process of walks, and fractional differential equations are obtained that describe the asymptotic behavior of the process. It is shown that, in the case of finite and infinite mathematical expectation of paths, these equations have a completely different form. To solve the obtained equations, the method of local estimation of the Monte Carlo method is described. The solution algorithm is described and the advantages and disadvantages of the considered method are indicated.

Asymptotic densities of planar Lévy walks: A nonisotropic case

Lévy walks are a particular type of continuous-time random walks which results in a super-diffusive spreading of an initially localized packet. The original one-dimensional model has a simple schematization that is based on starting a new unidirectional motion event either in the positive or in the negative direction. We consider two-dimensional generalization of Lévy walks in the form of the so-called XY model. It describes a particle moving with a constant velocity along one of the four basic directions and randomly switching between them when starting a new motion event. We address the ballistic regime and derive solutions for the asymptotic density profiles. The solutions have a form of first-order integrals which can be evaluated numerically. For specific values of parameters we derive an exact expression. The analytic results are in agreement with the results of finite-time numerical samplings.

Arcsine law and multistable Brownian dynamics in a tilted periodic potential

Multistability is one of the most important phenomena in dynamical systems, e.g., bistability enables the implementation of logic gates and therefore computation. Recently multistability has attracted a greatly renewed interest related to memristors and graphene structures, to name only a few. We investigate tristability in velocity dynamics of a Brownian particle subjected to a tilted periodic potential. It is demonstrated that the origin of this effect is attributed to the arcsine law for the velocity dynamics at the zero temperature limit. We analyze the impact of thermal fluctuations and construct the phase diagram for the stability of the velocity dynamics. It suggests an efficient strategy to control the multistability by changing solely the force acting on the particle or temperature of the system. Our findings for the paradigmatic model of nonequilibrium statistical physics apply to, inter alia, Brownian motors, Josephson junctions, cold atoms dwelling in optical lattices, and colloidal systems.

Self-reinforcing directionality generates truncated Lévy walks without the power-law assumption

We introduce a persistent random walk model with finite velocity and self-reinforcing directionality, which explains how exponentially distributed runs self-organize into truncated Lévy walks observed in active intracellular transport by Chen et al. [Nature Mater., 14, 589 (2015)10.1038/nmat4239]. We derive the nonhomogeneous in space and time, hyperbolic partial differential equation for the probability density function (PDF) of particle position. This PDF exhibits a bimodal density (aggregation phenomena) in the superdiffusive regime, which is not observed in classical linear hyperbolic and Lévy walk models. We find the exact solutions for the first and second moments and criteria for the transition to superdiffusion.

Large deviations of the ballistic Lévy walk model

We study the ballistic Lévy walk stemming from an infinite mean traveling time between collision events. Our study focuses on the density of spreading particles all starting from a common origin, which is limited by a "light" cone -v_{0}t<x<v_{0}t. In particular we study this density close to its maximum in the vicinity of the light cone. The spreading density follows the Lamperti-arcsine law describing typical fluctuations. However, this law blows up in the vicinity of the spreading horizon, which is nonphysical in the sense that any finite-time observation will never diverge. We claim that one can find two laws for the spatial density: The first one is the mentioned Lamperti-arcsine law describing the central part of the distribution, and the second is an infinite density illustrating the dynamics for x≃v_{0}t. We identify the relationship between a large position and the longest traveling time describing the single big jump principle. From the renewal theory we find that the distribution of rare events of the position is related to the derivative of the average of the number of renewals at a short "time" using a rate formalism.

Lévy walks in nonhomogeneous environments

The Lévy walk process with rests is discussed. The jumping time is governed by an α-stable distribution with α>1 while a waiting time distribution is Poissonian and involves a position-dependent rate which reflects a nonhomogeneous trap distribution. The master equation is derived and solved in the asymptotic limit for a power-law form of the jumping rate. The relative density of resting and flying particles appears time-dependent, and the asymptotic form of both distributions obeys a stretched-exponential shape at large time. The diffusion properties are discussed, and it is demonstrated that, due to the heterogeneous trap structure, the enhanced diffusion, observed for the homogeneous case, may turn to a subdiffusion. The density distributions and mean squared displacements are also evaluated from Monte Carlo simulations of individual trajectories.

Aging in mortal superdiffusive Lévy walkers

A growing body of literature examines the effects of superdiffusive subballistic movement premeasurement (aging or time lag) on observations arising from single-particle tracking. A neglected aspect is the finite lifetime of these Lévy walkers, be they proteins, cells, or larger structures. We examine the effects of aging on the motility of mortal walkers, and discuss the means by which permanent stopping of walkers may be categorized as arising from "natural" death or experimental artifacts such as low photostability or radiation damage. This is done by comparison of the walkers' mean squared displacement (MSD) with the front velocity of propagation of a group of walkers, which is found to be invariant under time lags. For any running time distribution of a mortal random walker, the MSD is tempered by the stopping rate θ. This provides a physical interpretation for truncated heavy-tailed diffusion processes and serves as a tool by which to better classify the underlying running time distributions of random walkers. Tempering of aged MSDs raises the issue of misinterpreting superdiffusive motion which appears Brownian or subdiffusive over certain time scales.

Using Hilbert transform and classical chains to simulate quantum walks

We propose a simulation strategy which uses a classical device of linearly coupled chain of springs to simulate quantum dynamics, in particular quantum walks. Through this strategy, we obtain the quantum wave function from the classical evolution. Specially, this goal is achieved with the classical momenta of the particles on the chain and their Hilbert transform, from which we construct the many-body momentum and Hilbert transformed momentum pair correlation functions yielding the real and imaginary parts of the wave function, respectively. With such a wave function, we show that the classical chain's energy and heat spreading densities can be related to the wave function's modulus square. This relation provides a new perspective to understand ballistic heat transport. The results here may give a definite answer to Feynman's idea of using a classical device to simulate quantum physics.

Lévy flights versus Lévy walks in bounded domains

Lévy flights and Lévy walks serve as two paradigms of random walks resembling common features but also bearing fundamental differences. One of the main dissimilarities is the discontinuity versus continuity of their trajectories and infinite versus finite propagation velocity. As a consequence, a well-developed theory of Lévy flights is associated with their pathological physical properties, which in turn are resolved by the concept of Lévy walks. Here, we explore Lévy flight and Lévy walk models on bounded domains, examining their differences and analogies. We investigate analytically and numerically whether and under which conditions both approaches yield similar results in terms of selected statistical observables characterizing the motion: the survival probability, mean first passage time, and stationary probability density functions. It is demonstrated that the similarity of the models is affected by the type of boundary conditions and the value of the stability index defining the asymptotics of the jump length distribution.

Aging ballistic Lévy walks

Aging can be observed for numerous physical systems. In such systems statistical properties [like probability distribution, mean square displacement (MSD), first-passage time] depend on a time span t_{a} between the initialization and the beginning of observations. In this paper we study aging properties of ballistic Lévy walks and two closely related jump models: wait-first and jump-first. We calculate explicitly their probability distributions and MSDs. It turns out that despite similarities these models react very differently to the delay t_{a}. Aging weakly affects the shape of probability density function and MSD of standard Lévy walks. For the jump models the shape of the probability density function is changed drastically. Moreover for the wait-first jump model we observe a different behavior of MSD when t_{a}≪t and t_{a}≫t.

Superdiffusive Dispersals Impart the Geometry of Underlying Random Walks

It is recognized now that a variety of real-life phenomena ranging from diffusion of cold atoms to the motion of humans exhibit dispersal faster than normal diffusion. Lévy walks is a model that excelled in describing such superdiffusive behaviors albeit in one dimension. Here we show that, in contrast to standard random walks, the microscopic geometry of planar superdiffusive Lévy walks is imprinted in the asymptotic distribution of the walkers. The geometry of the underlying walk can be inferred from trajectories of the walkers by calculating the analogue of the Pearson coefficient.

Explicit densities of multidimensional ballistic Lévy walks

Lévy walks have proved to be useful models of stochastic dynamics with a number of applications in the modeling of real-life phenomena. In this paper we derive explicit formulas for densities of the two- (2D) and three-dimensional (3D) ballistic Lévy walks, which are most important in applications. It turns out that in the 3D case the densities are given by elementary functions. The densities of the 2D Lévy walks are expressed in terms of hypergeometric functions and the right-side Riemann-Liouville fractional derivative, which allows us to efficiently evaluate them numerically. The theoretical results agree perfectly with Monte Carlo simulations.

Single integrodifferential wave equation for a Lévy walk

We derive the single integrodifferential wave equation for the probability density function of the position of a classical one-dimensional Lévy walk with continuous sample paths. This equation involves a classical wave operator together with memory integrals describing the spatiotemporal coupling of the Lévy walk. It is valid at all times, not only in the long time limit, and it does not involve any large-scale approximations. It generalizes the well-known telegraph or Cattaneo equation for the persistent random walk with the exponential switching time distribution. Several non-Markovian cases are considered when the particle's velocity alternates at the gamma and power-law distributed random times. In the strong anomalous case we obtain the asymptotic solution to the integrodifferential wave equation. We implement the nonlinear reaction term of Kolmogorov-Petrovsky-Piskounov type into our equation and develop the theory of wave propagation in reaction-transport systems involving Lévy diffusion.

Scaling forms of particle densities for Lévy walks and strong anomalous diffusion

We study the scaling behavior of particle densities for Lévy walks whose transition length r is coupled with the transition time t as |r|∝t^{α} with an exponent α>0. The transition-time distribution behaves as ψ(t)∝t^{-1-β} with β>0. For 1<β<2α and α≥1, particle displacements are characterized by a finite transition time and confinement to |r|<t^{α} while the marginal distribution of transition lengths is heavy tailed. These characteristics give rise to the existence of two scaling forms for the particle density, one that is valid at particle displacements |r|≪t^{α} and one at |r|≲t^{α}. As a consequence, the Lévy walk displays strong anomalous diffusion in the sense that the average absolute moments 〈|r|^{q}〉 scale as t^{qν(q)} with ν(q) piecewise linear above and below a critical value q_{c}. We derive explicit expressions for the scaling forms of the particle densities and determine the scaling of the average absolute moments. We find that 〈|r|^{q}〉∝t^{qα/β} for q<q_{c}=β/α and 〈|r|^{q}〉∝t^{1+αq-β} for q>q_{c}. These results give insight into the possible origins of strong anomalous diffusion and anomalous behaviors in disordered systems in general.