Subdiffusion in an external potential: Anomalous effects hiding behind normal behavior

Anomalous diffusion in external-force-affected deterministic systems

This study investigates the impact of external forces on the movement of particles, specifically focusing on a type of box piecewise linear map that generates normal diffusion akin to Brownian motion. Through numerical methods, the research delves into the effects of two distinct external forces: linear forces linked to the particle's current position and periodic sinusoidal forces related to time. The results uncover anomalous dynamical behavior characterized by nonlinear growth in the ensemble-averaged mean-squared displacement (EAMSD), aging, and ergodicity breaking. Notably, the diffusion pattern of particles under linear external forces resembles an Ehrenfest double urn model, with its asymptotic EAMSD coinciding with the Langevin equation under linear potential. Meanwhile, particle movement influenced by periodic sinusoidal forces corresponds to an inhomogeneous Markov chain, with its external force amplitude and diffusion coefficient function exhibiting a "multipeak" fractal structure. The study also provides insights into the formation of this structure through the turnstiles dynamics.

Ergodic criterion of a random diffusivity model

The random diffusivity, initially proposed to explain Brownian yet non-Gaussian diffusion, has garnered significant attention due to its capacity not only for elucidating the internal physical mechanism of non-Gaussian diffusion, but also for establishing an analytical framework to characterize particle motion in complex environments. In this paper, based on the correlation function C(t_{1},t_{2})=〈D(t_{1})D(t_{2})〉 of random diffusivity D(t), we quantitatively propose a general criterion of determining the ergodic property of the Langevin equation with the arbitrary random diffusivity D(t). Due to the critical role of correlation function C(t_{1},t_{2}), we derive the criterion for the two cases with stationary diffusivity or nonstationary diffusivity, respectively. By utilizing the quantitative criterion, we can directly judge the ergodic properties of the random diffusivity model based on the correlation function C(t_{1},t_{2}) of random diffusivity D(t). Several typical diffusivities, including the common square of the Brownian motion and of the (fractional) Ornstein-Uhlenbeck process, are found to contribute to different ergodic properties, which validates our proposed criterion built on the correlation function C(t_{1},t_{2}).

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.

Random diffusivity processes in an external force field

Brownian yet non-Gaussian processes have recently been observed in numerous biological systems, and corresponding theories have been constructed based on random diffusivity models. Considering the particularity of random diffusivity, this paper studies the effect of an external force acting on two kinds of random diffusivity models whose difference is embodied in whether the fluctuation-dissipation theorem is valid. Based on the two random diffusivity models, we derive the Fokker-Planck equations with an arbitrary external force, and we analyze various observables in the case with a constant force, including the Einstein relation, the moments, the kurtosis, and the asymptotic behaviors of the probability density function of particle displacement at different timescales. Both the theoretical results and numerical simulations of these observables show a significant difference between the two kinds of random diffusivity models, which implies the important role of the fluctuation-dissipation theorem in random diffusivity systems.

Langevin picture of Lévy walk in a constant force field

Lévy walk is a practical model and has wide applications in various fields. Here we focus on the effect of an external constant force on the Lévy walk with the exponent of the power-law-distributed flight time α∈(0,2). We add the term Fη(s) [η(s) is the Lévy noise] on a subordinated Langevin system to characterize such a constant force, as it is effective on the velocity process for all physical time after the subordination. We clearly show the effect of the constant force F on this Langevin system and find this system is like the continuous limit of the collision model. The first moments of velocity processes for these two models are consistent. In particular, based on the velocity correlation function derived from our subordinated Langevin equation, we investigate more interesting statistical quantities, such as the ensemble- and time-averaged mean-squared displacements. Under the influence of constant force, the diffusion of particles becomes faster. Finally, the superballistic diffusion and the nonergodic behavior are verified by the simulations with different α.

Subdiffusive continuous-time random walks with stochastic resetting

We analyze two models of subdiffusion with stochastic resetting. Each of them consists of two parts: subdiffusion based on the continuous-time random walk scheme and independent resetting events generated uniformly in time according to the Poisson point process. In the first model the whole process is reset to the initial state, whereas in the second model only the position is subject to resets. The distinction between these two models arises from the non-Markovian character of the subdiffusive process. We derive exact expressions for the two lowest moments of the full propagator, stationary distributions, and first hitting time statistics. We also show, with an example of a constant drift, how these models can be generalized to include external forces. Possible applications to data analysis and modeling of biological systems are also discussed.

Subdiffusion in an external force field

The phenomena of subdiffusion are widely observed in physical and biological systems. To investigate the effects of external potentials, say, harmonic potential, linear potential, and time-dependent force, we study the subdiffusion described by the subordinated Langevin equation with white Gaussian noise or, equivalently, by the single Langevin equation with compound noise. If the force acts on the subordinated process, it keeps working all the time; otherwise, the force just exerts an influence on the system at the moments of jump. Some common statistical quantities, such as the ensemble- and time-averaged mean squared displacements, position autocorrelation function, correlation coefficient, and generalized Einstein relation, are discussed to distinguish the effects of various forces and different patterns of acting. The corresponding Fokker-Planck equations are also presented. All the stochastic processes discussed here are nonstationary, nonergodic, and aging.

Non-Markovian intracellular transport with sub-diffusion and run-length dependent detachment rate

Intracellular transport of organelles is fundamental to cell function and health. The mounting evidence suggests that this transport is in fact anomalous. However, the reasons for the anomaly is still under debate. We examined experimental trajectories of organelles inside a living cell and propose a mathematical model that describes the previously reported transition from sub-diffusive to super-diffusive motion. In order to explain super-diffusive behaviour at long times, we introduce non-Markovian detachment kinetics of the cargo: the rate of detachment is inversely proportional to the time since the last attachment. Recently, we observed the non-Markovian detachment rate experimentally in eukaryotic cells. Here we further discuss different scenarios of how this effective non-Markovian detachment rate could arise. The non-Markovian model is successful in simultaneously describing the time averaged variance (the time averaged mean squared displacement corrected for directed motion), the mean first passage time of trajectories and the multiple peaks observed in the distributions of cargo velocities. We argue that non-Markovian kinetics could be biologically beneficial compared to the Markovian kinetics commonly used for modelling, by increasing the average distance the cargoes travel when a microtubule is blocked by other filaments. In turn, sub-diffusion allows cargoes to reach neighbouring filaments with higher probability, which promotes active motion along the microtubules.

Emergence of Lévy walks in systems of interacting individuals

We propose a model of superdiffusive Lévy walk as an emergent nonlinear phenomenon in systems of interacting individuals. The aim is to provide a qualitative explanation of recent experiments [G. Ariel et al., Nat. Commun. 6, 8396 (2015)2041-172310.1038/ncomms9396] revealing an intriguing behavior: swarming bacteria fundamentally change their collective motion from simple diffusion into a superdiffusive Lévy walk dynamics. We introduce microscopic mean-field kinetic equations in which we combine two key ingredients: (1) alignment interactions between individuals and (2) non-Markovian effects. Our interacting run-and-tumble model leads to the superdiffusive growth of the mean-squared displacement and the power-law distribution of run length with infinite variance. The main result is that the superdiffusive behavior emerges as a cooperative effect without using the standard assumption of the power-law distribution of run distances from the inception. At the same time, we find that the collision and repulsion interactions lead to the density-dependent exponential tempering of power-law distributions. This qualitatively explains the experimentally observed transition from superdiffusion to the diffusion of mussels as their density increases [M. de Jager et al., Proc. R. Soc. B 281, 20132605 (2014)PRLBA40962-845210.1098/rspb.2013.2605].

Langevin formulation of a subdiffusive continuous-time random walk in physical time

Systems living in complex nonequilibrated environments often exhibit subdiffusion characterized by a sublinear power-law scaling of the mean square displacement. One of the most common models to describe such subdiffusive dynamics is the continuous-time random walk (CTRW). Stochastic trajectories of a CTRW can be described in terms of the subordination of a normal diffusive process by an inverse Lévy-stable process. Here, we propose an equivalent Langevin formulation of a force-free CTRW without subordination. By introducing a different type of non-Gaussian noise, we are able to express the CTRW dynamics in terms of a single Langevin equation in physical time with additive noise. We derive the full multipoint statistics of this noise and compare it with the scaled Brownian motion (SBM), an alternative stochastic model describing subdiffusive dynamics. Interestingly, these two noises are identical up to the second order correlation functions, but different in the higher order statistics. We extend our formalism to general waiting time distributions and force fields and compare our results with those of the SBM. In the presence of external forces, our proposed noise generates a different class of stochastic processes, resembling a CTRW but with forces acting at all times.