Memory-induced diffusive-superdiffusive transition: Ensemble and time-averaged observables

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.

Superdiffusion in self-reinforcing run-and-tumble model with rests

This paper introduces a run-and-tumble model with self-reinforcing directionality and rests. We derive a single governing hyperbolic partial differential equation for the probability density of random-walk position, from which we obtain the second moment in the long-time limit. We find the criteria for the transition between superdiffusion and diffusion caused by the addition of a rest state. The emergence of superdiffusion depends on both the parameter representing the strength of self-reinforcement and the ratio between mean running and resting times. The mean running time must be at least 2/3 of the mean resting time for superdiffusion to be possible. Monte Carlo simulations validate this theoretical result. This work demonstrates the possibility of extending the telegrapher's (or Cattaneo) equation by adding self-reinforcing directionality so that superdiffusion occurs even when rests are introduced.

Reconstruction of Diffusion Coefficients and Power Exponents from Single Lagrangian Trajectories

Using discrete wavelets, a novel technique is developed to estimate turbulent diffusion coefficients and power exponents from single Lagrangian particle trajectories. The technique differs from the classical approach (Davis (1991)’s technique) because averaging over a statistical ensemble of the mean square displacement (<X2>) is replaced by averaging along a single Lagrangian trajectory X(t) = {X(t), Y(t)}. Metzler et al. (2014) have demonstrated that for an ergodic (for example, normal diffusion) flow, the mean square displacement is <X2> = limT→∞τX2(T,s), where τX2 (T, s) = 1/(T − s) ∫0T−s(X(t+Δt) − X(t))2 dt, T and s are observational and lag times but for weak non-ergodic (such as super-diffusion and sub-diffusion) flows <X2> = limT→∞≪τX2(T,s)≫, where ≪…≫ is some additional averaging. Numerical calculations for surface drifters in the Black Sea and isobaric RAFOS floats deployed at mid depths in the California Current system demonstrated that the reconstructed diffusion coefficients were smaller than those calculated by Davis (1991)’s technique. This difference is caused by the choice of the Lagrangian mean. The technique proposed here is applied to the analysis of Lagrangian motions in the Black Sea (horizontal diffusion coefficients varied from 105 to 106 cm2/s) and for the sub-diffusion of two RAFOS floats in the California Current system where power exponents varied from 0.65 to 0.72. RAFOS float motions were found to be strongly non-ergodic and non-Gaussian.

Giant leaps and long excursions: Fluctuation mechanisms in systems with long-range memory

We analyze large deviations of time-averaged quantities in stochastic processes with long-range memory, where the dynamics at time t depends itself on the value q_{t} of the time-averaged quantity. First we consider the elephant random walk and a Gaussian variant of this model, identifying two mechanisms for unusual fluctuation behavior, which differ from the Markovian case. In particular, the memory can lead to large-deviation principles with reduced speeds and to nonanalytic rate functions. We then explain how the mechanisms operating in these two models are generic for memory-dependent dynamics and show other examples including a non-Markovian simple exclusion process.

Infinite invariant density in a semi-Markov process with continuous state variables

We report on a fundamental role of a non-normalized formal steady state, i.e., an infinite invariant density, in a semi-Markov process where the state is determined by the interevent time of successive renewals. The state describes certain observables found in models of anomalous diffusion, e.g., the velocity in the generalized Lévy walk model and the energy of a particle in the trap model. In our model, the interevent-time distribution follows a fat-tailed distribution, which makes the state value more likely to be zero because long interevent times imply small state values. We find two scaling laws describing the density for the state value, which accumulates in the vicinity of zero in the long-time limit. These laws provide universal behaviors in the accumulation process and give the exact expression of the infinite invariant density. Moreover, we provide two distributional limit theorems for time-averaged observables in these nonstationary processes. We show that the infinite invariant density plays an important role in determining the distribution of time averages.

Time-Dependent Diffusion Coefficients for Chaotic Advection due to Fluctuations of Convective Rolls

The properties of chaotic advection arising from defect turbulence, that is, weak turbulence in the electroconvection of nematic liquid crystals, were experimentally investigated. Defect turbulence is a phenomenon in which fluctuations of convective rolls arise and are globally disturbed while maintaining convective rolls locally. The time-dependent diffusion coefficient, as measured from the motion of a tagged particle driven by the turbulence, was used to clarify the dependence of the type of diffusion on coarse-graining time. The results showed that, as coarse-graining time increases, the type of diffusion changes from superdiffusion → subdiffusion → normal diffusion. The change in diffusive properties over the observed timescale reflects the coexistence of local order and global disorder in the defect turbulence.

Biased continuous-time random walks for ordinary and equilibrium cases: facilitation of diffusion, ergodicity breaking and ageing

We examine renewal processes with power-law waiting time distributions (WTDs) and non-zero drift via computing analytically and by computer simulations their ensemble and time averaged spreading characteristics. All possible values of the scaling exponent α are considered for the WTD ψ(t) ∼ 1/t1+α. We treat continuous-time random walks (CTRWs) with 0 < α < 1 for which the mean waiting time diverges, and investigate the behaviour of the process for both ordinary and equilibrium CTRWs for 1 < α < 2 and α > 2. We demonstrate that in the presence of a drift CTRWs with α < 1 are ageing and non-ergodic in the sense of the non-equivalence of their ensemble and time averaged displacement characteristics in the limit of lag times much shorter than the trajectory length. In the sense of the equivalence of ensemble and time averages, CTRW processes with 1 < α < 2 are ergodic for the equilibrium and non-ergodic for the ordinary situation. Lastly, CTRW renewal processes with α > 2-both for the equilibrium and ordinary situation-are always ergodic. For the situations 1 < α < 2 and α > 2 the variance of the diffusion process, however, depends on the initial ensemble. For biased CTRWs with α > 1 we also investigate the behaviour of the ergodicity breaking parameter. In addition, we demonstrate that for biased CTRWs the Einstein relation is valid on the level of the ensemble and time averaged displacements, in the entire range of the WTD exponent α.

Time averages and their statistical variation for the Ornstein-Uhlenbeck process: Role of initial particle distributions and relaxation to stationarity

How ergodic is diffusion under harmonic confinements? How strongly do ensemble- and time-averaged displacements differ for a thermally-agitated particle performing confined motion for different initial conditions? We here study these questions for the generic Ornstein-Uhlenbeck (OU) process and derive the analytical expressions for the second and fourth moment. These quantifiers are particularly relevant for the increasing number of single-particle tracking experiments using optical traps. For a fixed starting position, we discuss the definitions underlying the ensemble averages. We also quantify effects of equilibrium and nonequilibrium initial particle distributions onto the relaxation properties and emerging nonequivalence of the ensemble- and time-averaged displacements (even in the limit of long trajectories). We derive analytical expressions for the ergodicity breaking parameter quantifying the amplitude scatter of individual time-averaged trajectories, both for equilibrium and out-of-equilibrium initial particle positions, in the entire range of lag times. Our analytical predictions are in excellent agreement with results of computer simulations of the Langevin equation in a parabolic potential. We also examine the validity of the Einstein relation for the ensemble- and time-averaged moments of the OU-particle. Some physical systems, in which the relaxation and nonergodic features we unveiled may be observable, are discussed.

Exact Results for the Nonergodicity of d-Dimensional Generalized Lévy Walks

We provide analytical results for the ensemble-averaged and time-averaged squared displacement, and the randomness of the latter, in the full two-dimensional parameter space of the d-dimensional generalized Lévy walk introduced by Shlesinger et al. [Phys. Rev. Lett. 58, 1100 (1987)PRLTAO0031-900710.1103/PhysRevLett.58.1100]. In certain regions of the parameter plane, we obtain surprising results such as the divergence of the mean-squared displacements, the divergence of the ergodicity breaking parameter despite a finite mean-squared displacement, and subdiffusion which appears superdiffusive when one only considers time averages.

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.

Nonstationary Markovian replication process causing diverse diffusions

We introduce a single generative mechanism that can be used to describe diverse nonstationary diffusions. A nonstationary Markovian replication process for steps is considered for which we derive analytically the time evolution of the probability distribution of the walker's displacement and the generalized telegrapher equation with time-varying coefficients, and we find that diffusivity can be determined by temporal changes of replication of an immediate step. By controlling the replications, we realize diverse diffusions such as alternating diffusion, superdiffusion, subdiffusion, and marginal diffusion, which originate from oscillating, increasing, decreasing, and slowly increasing or decreasing replications with time, respectively.