Refine
Year of publication
Language
- English (78)
Is part of the Bibliography
- yes (78)
Keywords
- anomalous diffusion (24)
- diffusion (11)
- stochastic processes (10)
- Levy flights (6)
- ageing (5)
- living cells (4)
- financial time series (3)
- first-hitting time (3)
- first-passage time (3)
- fractional Brownian motion (3)
Institute
Neutrophil granulocytes are essential for the first host defense. After leaving the blood circulation they migrate efficiently towards sites of inflammation. They are guided by chemoattractants released from cells within the inflammatory foci. On a cellular level, directional migration is a consequence of cellular front-rear asymmetry which is induced by the concentration gradient of the chemoattractants. The generation and maintenance of this asymmetry, however, is not yet fully understood. Here we analyzed the paths of chemotacting neutrophils with different stochastic models to gain further insight into the underlying mechanisms. Wildtype chemotacting neutrophils show an anomalous superdiffusive behavior. CXCR2 blockade and TRPC6-knockout cause the tempering of temporal correlations and a reduction of chemotaxis. Importantly, such tempering is found both in vitro and in vivo. These findings indicate that the maintenance of anomalous dynamics is crucial for chemotactic behavior and the search efficiency of neutrophils.
The motility of neutrophils and their ability to sense and to react to chemoattractants in their environment are of central importance for the innate immunity. Neutrophils are guided towards sites of inflammation following the activation of G-protein coupled chemoattractant receptors such as CXCR2 whose signaling strongly depends on the activity of Ca2+ permeable TRPC6 channels. It is the aim of this study to analyze data sets obtained in vitro (murine neutrophils) and in vivo (zebrafish neutrophils) with a stochastic mathematical model to gain deeper insight into the underlying mechanisms. The model is based on the analysis of trajectories of individual neutrophils. Bayesian data analysis, including the covariances of positions for fractional Brownian motion as well as for exponentially and power-law tempered model variants, allows the estimation of parameters and model selection. Our model-based analysis reveals that wildtype neutrophils show pure superdiffusive fractional Brownian motion. This so-called anomalous dynamics is characterized by temporal long-range correlations for the movement into the direction of the chemotactic CXCL1 gradient. Pure superdiffusion is absent vertically to this gradient. This points to an asymmetric 'memory' of the migratory machinery, which is found both in vitro and in vivo. CXCR2 blockade and TRPC6-knockout cause tempering of temporal correlations in the chemotactic gradient. This can be interpreted as a progressive loss of memory, which leads to a marked reduction of chemotaxis and search efficiency of neutrophils. In summary, our findings indicate that spatially differential regulation of anomalous dynamics appears to play a central role in guiding efficient chemotactic behavior.
We analyse mobile-immobile transport of particles that switch between the mobile and immobile phases with finite rates. Despite this seemingly simple assumption of Poissonian switching, we unveil a rich transport dynamics including significant transient anomalous diffusion and non-Gaussian displacement distributions. Our discussion is based on experimental parameters for tau proteins in neuronal cells, but the results obtained here are expected to be of relevance for a broad class of processes in complex systems. Specifically, we obtain that, when the mean binding time is significantly longer than the mean mobile time, transient anomalous diffusion is observed at short and intermediate time scales, with a strong dependence on the fraction of initially mobile and immobile particles. We unveil a Laplace distribution of particle displacements at relevant intermediate time scales. For any initial fraction of mobile particles, the respective mean squared displacement (MSD) displays a plateau. Moreover, we demonstrate a short-time cubic time dependence of the MSD for immobile tracers when initially all particles are immobile.
Generalized (non-Markovian) diffusion equations with different memory kernels and subordination schemes based on random time change in the Brownian diffusion process are popular mathematical tools for description of a variety of non-Fickian diffusion processes in physics, biology, and earth sciences. Some of such processes (notably, the fluid limits of continuous time random walks) allow for either kind of description, but other ones do not. In the present work we discuss the conditions under which a generalized diffusion equation does correspond to a subordination scheme, and the conditions under which a subordination scheme does possess the corresponding generalized diffusion equation. Moreover, we discuss examples of random processes for which only one, or both kinds of description are applicable.
The escape from a potential well is an archetypal problem in the study of stochastic dynamical systems, representing real-world situations from chemical reactions to leaving an established home range in movement ecology. Concurrently, Levy noise is a well-established approach to model systems characterized by statistical outliers and diverging higher order moments, ranging from gene expression control to the movement patterns of animals and humans. Here, we study the problem of Levy noise-driven escape from an almost rectangular, arctangent potential well restricted by two absorbing boundaries, mostly under the action of the Cauchy noise. We unveil analogies of the observed transient dynamics to the general properties of stationary states of Levy processes in single-well potentials. The first-escape dynamics is shown to exhibit exponential tails. We examine the dependence of the escape on the shape parameters, steepness, and height of the arctangent potential. Finally, we explore in detail the behavior of the probability densities of the first-escape time and the last-hitting point.
Fractional Brownian motion in superharmonic potentials and non-Boltzmann stationary distributions
(2021)
We study the stochastic motion of particles driven by long-range correlated fractional Gaussian noise (FGN) in a superharmonic external potential of the form U(x) proportional to x(2n) (n is an element of N). When the noise is considered to be external, the resulting overdamped motion is described by the non-Markovian Langevin equation for fractional Brownian motion. For this case we show the existence of long time, stationary probability density functions (PDFs) the shape of which strongly deviates from the naively expected Boltzmann PDF in the confining potential U(x). We analyse in detail the temporal approach to stationarity as well as the shape of the non-Boltzmann stationary PDF. A typical characteristic is that subdiffusive, antipersistent (with negative autocorrelation) motion tends to effect an accumulation of probability close to the origin as compared to the corresponding Boltzmann distribution while the opposite trend occurs for superdiffusive (persistent) motion. For this latter case this leads to distinct bimodal shapes of the PDF. This property is compared to a similar phenomenon observed for Markovian Levy flights in superharmonic potentials. We also demonstrate that the motion encoded in the fractional Langevin equation driven by FGN always relaxes to the Boltzmann distribution, as in this case the fluctuation-dissipation theorem is fulfilled.
Numerous examples for a priori unexpected non-Gaussian behaviour for normal and anomalous diffusion have recently been reported in single-particle tracking experiments. Here, we address the case of non-Gaussian anomalous diffusion in terms of a random-diffusivity mechanism in the presence of power-law correlated fractional Gaussian noise. We study the ergodic properties of this model via examining the ensemble- and time-averaged mean-squared displacements as well as the ergodicity breaking parameter EB quantifying the trajectory-to-trajectory fluctuations of the latter. For long measurement times, interesting crossover behaviour is found as function of the correlation time tau characterising the diffusivity dynamics. We unveil that at short lag times the EB parameter reaches a universal plateau. The corresponding residual value of EB is shown to depend only on tau and the trajectory length. The EB parameter at long lag times, however, follows the same power-law scaling as for fractional Brownian motion. We also determine a corresponding plateau at short lag times for the discrete representation of fractional Brownian motion, absent in the continuous-time formulation. These analytical predictions are in excellent agreement with results of computer simulations of the underlying stochastic processes. Our findings can help distinguishing and categorising certain nonergodic and non-Gaussian features of particle displacements, as observed in recent single-particle tracking experiments.
We present a diffusion-based simulation and theoretical models for explanation of the photoluminescence (PL) emission intensity in semiconductor nanoplatelets. It is shown that the shape of the PL intensity curves can be reproduced by the interplay of recombination, diffusion and trapping of excitons. The emission intensity at short times is purely exponential and is defined by recombination. At long times, it is governed by the release of excitons from surface traps and is characterized by a power-law tail. We show that the crossover from one limit to another is controlled by diffusion properties. This intermediate region exhibits a rich behaviour depending on the value of diffusivity. The proposed approach reproduces all the features of experimental curves measured for different nanoplatelet systems.
We investigate both analytically and by computer simulations the ensemble- and time-averaged, nonergodic, and aging properties of massive particles diffusing in a medium with a time dependent diffusivity. We call this stochastic diffusion process the (aging) underdamped scaled Brownian motion (UDSBM). We demonstrate how the mean squared displacement (MSD) and the time-averaged MSD of UDSBM are affected by the inertial term in the Langevin equation, both at short, intermediate, and even long diffusion times. In particular, we quantify the ballistic regime for the MSD and the time-averaged MSD as well as the spread of individual time-averaged MSD trajectories. One of the main effects we observe is that, both for the MSD and the time-averaged MSD, for superdiffusive UDSBM the ballistic regime is much shorter than for ordinary Brownian motion. In contrast, for subdiffusive UDSBM, the ballistic region extends to much longer diffusion times. Therefore, particular care needs to be taken under what conditions the overdamped limit indeed provides a correct description, even in the long time limit. We also analyze to what extent ergodicity in the Boltzmann-Khinchin sense in this nonstationary system is broken, both for subdiffusive and superdiffusive UDSBM. Finally, the limiting case of ultraslow UDSBM is considered, with a mixed logarithmic and power-law dependence of the ensemble-and time-averaged MSDs of the particles. In the limit of strong aging, remarkably, the ordinary UDSBM and the ultraslow UDSBM behave similarly in the short time ballistic limit. The approaches developed here open ways for considering other stochastic processes under physically important conditions when a finite particle mass and aging in the system cannot be neglected.
We obtain a generalized diffusion equation in modified or Riemann-Liouville form from continuous time random walk theory. The waiting time probability density function and mean squared displacement for different forms of the equation are explicitly calculated. We show examples of generalized diffusion equations in normal or Caputo form that encode the same probability distribution functions as those obtained from the generalized diffusion equation in modified form. The obtained equations are general and many known fractional diffusion equations are included as special cases.
Random search with resetting
(2018)
We provide a unified renewal approach to the problem of random search for several targets under resetting. This framework does not rely on specific properties of the search process and resetting procedure, allows for simpler derivation of known results, and leads to new ones. Concentrating on minimizing the mean hitting time, we show that resetting at a constant pace is the best possible option if resetting helps at all, and derive the equation for the optimal resetting pace. No resetting may be a better strategy if without resetting the probability of not finding a target decays with time to zero exponentially or faster. We also calculate splitting probabilities between the targets, and define the limits in which these can be manipulated by changing the resetting procedure. We moreover show that the number of moments of the hitting time distribution under resetting is not less than the sum of the numbers of moments of the resetting time distribution and the hitting time distribution without resetting.
Percolation networks have been widely used in the description of porous media but are now found to be relevant to understand the motion of particles in cellular membranes or the nucleus of biological cells. Random walks on the infinite cluster at criticality of a percolation network are asymptotically ergodic. On any finite size cluster of the network stationarity is reached at finite times, depending on the cluster's size. Despite of this we here demonstrate by combination of analytical calculations and simulations that at criticality the disorder and cluster size average of the ensemble of clusters leads to a non-vanishing variance of the time averaged mean squared displacement, regardless of the measurement time. Fluctuations of this relevant experimental quantity due to the disorder average of such ensembles are thus persistent and non-negligible. The relevance of our results for single particle tracking analysis in complex and biological systems is discussed.
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 outof-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.
In this paper we study the large deviations of time averaged mean square displacement (TAMSD) for Gaussian processes. The theory of large deviations is related to the exponential decay of probabilities of large fluctuations in random systems. From the mathematical point of view a given statistics satisfies the large deviation principle, if the probability that it belongs to a certain range decreases exponentially. The TAMSD is one of the main statistics used in the problem of anomalous diffusion detection. Applying the theory of generalized chi-squared distribution and sub-gamma random variables we prove the upper bound for large deviations of TAMSD for Gaussian processes. As a special case we consider fractional Brownian motion, one of the most popular models of anomalous diffusion. Moreover, we derive the upper bound for large deviations of the estimator for the anomalous diffusion exponent. (C) 2018 Elsevier B.V. All rights reserved.
Levy walks are continuous time random walks with spatio-temporal coupling of jump lengths and waiting times, often used to model superdiffusive spreading processes such as animals searching for food, tracer motion in weakly chaotic systems, or even the dynamics in quantum systems such as cold atoms. In the simplest version Levy walks move with a finite speed. Here, we present an extension of the Levy walk scenario for the case when external force fields influence the motion. The resulting motion is a combination of the response to the deterministic force acting on the particle, changing its velocity according to the principle of total energy conservation, and random velocity reversals governed by the distribution of waiting times. For the fact that the motion stays conservative, that is, on a constant energy surface, our scenario is fundamentally different from thermal motion in the same external potentials. In particular, we present results for the velocity and position distributions for single well potentials of different steepness. The observed dynamics with its continuous velocity changes enriches the theory of Levy walk processes and will be of use in a variety of systems, for which the particles are externally confined.
We study generalized diffusion-wave equation in which the second order time derivative is replaced by an integro-differential operator. It yields time fractional and distributed order time fractional diffusion-wave equations as particular cases. We consider different memory kernels of the integro-differential operator, derive corresponding fundamental solutions, specify the conditions of their non-negativity and calculate the mean squared displacement for all cases. In particular, we introduce and study generalized diffusion-wave equations with a regularized Prabhakar derivative of single and distributed orders. The equations considered can be used for modeling the broad spectrum of anomalous diffusion processes and various transitions between different diffusion regimes.
A rapidly increasing number of systems is identified in which the stochastic motion of tracer particles follows the Brownian law < r(2)(t)> similar or equal to Dt yet the distribution of particle displacements is strongly non-Gaussian. A central approach to describe this effect is the diffusing diffusivity (DD) model in which the diffusion coefficient itself is a stochastic quantity, mimicking heterogeneities of the environment encountered by the tracer particle on its path. We here quantify in terms of analytical and numerical approaches the first passage behaviour of the DD model. We observe significant modifications compared to Brownian-Gaussian diffusion, in particular that the DD model may have a faster first passage dynamics. Moreover we find a universal crossover point of the survival probability independent of the initial condition.
Complex systems are known to display anomalous diffusion, whose signature is a space/time scaling x similar to t(delta) with delta not equal 1/2 in the probability density function (PDF). Anomalous diffusion can emerge jointly with both Gaussian, e.g. fractional Brownian motion, and power-law decaying distributions, e.g. Levy Flights or Levy Walks (LWs). Levy flights get anomalous scaling, but, being jumps of any size allowed even at short times, have infinite position variance, infinite energy and discontinuous paths. LWs, which are based on random trapping events, overcome these limitations: they resemble a Levy-type power-law distribution that is truncated in the large displacement range and have finite moments, finite energy and, even with discontinuous velocity, they are continuous. However, LWs do not take into account the role of strong heterogeneity in many complex systems, such as biological transport in the crowded cell environment. In this work we propose and discuss a model describing a heterogeneous ensemble of Brownian particles (HEBP). Velocity of each single particle obeys a standard underdamped Langevin equation for the velocity, with linear friction term and additive Gaussian noise. Each particle is characterized by its own relaxation time and velocity diffusivity. We show that, for proper distributions of relaxation time and velocity diffusivity, the HEBP resembles some LW statistical features, in particular power-law decaying PDF, long-range correlations and anomalous diffusion, at the same time keeping finite position moments and finite energy. The main differences between the HEBP model and two different LWs are investigated, finding that, even when both velocity and position PDFs are similar, they differ in four main aspects: (i) LWs are biscaling, while HEBP is monoscaling; (ii) a transition from anomalous (delta = 1/2) to normal (delta = 1/2) diffusion in the long-time regime is seen in the HEBP and not in LWs; (iii) the power-law index of the position PDF and the space/time diffusion scaling are independent in the HEBP, while they both depend on the scaling of the interevent time PDF in LWs; (iv) at variance with LWs, our HEBP model obeys a fluctuation-dissipation theorem.
Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories
(2021)
Extensive time-series encoding the position of particles such as viruses, vesicles, or individualproteins are routinely garnered insingle-particle tracking experiments or supercomputing studies.They contain vital clues on how viruses spread or drugs may be delivered in biological cells.Similar time-series are being recorded of stock values in financial markets and of climate data.Such time-series are most typically evaluated in terms of time-averaged mean-squareddisplacements (TAMSDs), which remain random variables for finite measurement times. Theirstatistical properties are different for differentphysical stochastic processes, thus allowing us toextract valuable information on the stochastic process itself. To exploit the full potential of thestatistical information encoded in measured time-series we here propose an easy-to-implementand computationally inexpensive new methodology, based on deviations of the TAMSD from itsensemble average counterpart. Specifically, we use the upper bound of these deviations forBrownian motion (BM) to check the applicability of this approach to simulated and real data sets.By comparing the probability of deviations fordifferent data sets, we demonstrate how thetheoretical bound for BM reveals additional information about observed stochastic processes. Weapply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracerbeads measured in mucin hydrogels, and of geographic surface temperature anomalies. Ouranalysis shows how the large-deviation properties can be efficiently used as a simple yet effectiveroutine test to reject the BM hypothesis and unveil relevant information on statistical propertiessuch as ergodicity breaking and short-time correlations.
Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories
(2021)
Extensive time-series encoding the position of particles such as viruses, vesicles, or individualproteins are routinely garnered insingle-particle tracking experiments or supercomputing studies.They contain vital clues on how viruses spread or drugs may be delivered in biological cells.Similar time-series are being recorded of stock values in financial markets and of climate data.Such time-series are most typically evaluated in terms of time-averaged mean-squareddisplacements (TAMSDs), which remain random variables for finite measurement times. Theirstatistical properties are different for differentphysical stochastic processes, thus allowing us toextract valuable information on the stochastic process itself. To exploit the full potential of thestatistical information encoded in measured time-series we here propose an easy-to-implementand computationally inexpensive new methodology, based on deviations of the TAMSD from itsensemble average counterpart. Specifically, we use the upper bound of these deviations forBrownian motion (BM) to check the applicability of this approach to simulated and real data sets.By comparing the probability of deviations fordifferent data sets, we demonstrate how thetheoretical bound for BM reveals additional information about observed stochastic processes. Weapply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracerbeads measured in mucin hydrogels, and of geographic surface temperature anomalies. Ouranalysis shows how the large-deviation properties can be efficiently used as a simple yet effectiveroutine test to reject the BM hypothesis and unveil relevant information on statistical propertiessuch as ergodicity breaking and short-time correlations.
We investigate an intermittent stochastic process in which the diffusive motion with time-dependent diffusion coefficient D(t)∼tα−1 with α>0 (scaled Brownian motion) is stochastically reset to its initial position, and starts anew. In the present work we discuss the situation in which the memory on the value of the diffusion coefficient at a resetting time is erased, so that the whole process is a fully renewal one. The situation when the resetting of the coordinate does not affect the diffusion coefficient's time dependence is considered in the other work of this series [A. S. Bodrova et al., Phys. Rev. E 100, 012119 (2019)]. We show that the properties of the probability densities in such processes (erasing or retaining the memory on the diffusion coefficient) are vastly different. In addition we discuss the first-passage properties of the scaled Brownian motion with renewal resetting and consider the dependence of the efficiency of search on the parameters of the process.
We investigate an intermittent stochastic process in which diffusive motion with a time-dependent diffusion coefficient, D(t)∼tα−1, α>0 (scaled Brownian motion), is stochastically reset to its initial position and starts anew. The resetting follows a renewal process with either an exponential or a power-law distribution of the waiting times between successive renewals. The resetting events, however, do not affect the time dependence of the diffusion coefficient, so that the whole process appears to be a nonrenewal one. We discuss the mean squared displacement of a particle and the probability density function of its positions in this process. We show that scaled Brownian motion with resetting demonstrates rich behavior whose properties essentially depend on the interplay of the parameters of the resetting process and the particle's displacement infree motion. The motion of particles can remain almost unaffected by resetting but can also get slowed down or even be completely suppressed. Especially interesting are the nonstationary situations in which the mean squared displacement stagnates but the distribution of positions does not tend to any steady state. This behavior is compared to the situation [discussed in the companion paper; A. S. Bodrova et al., Phys. Rev. E 100, 012120 (2019)] in which the memory of the value of the diffusion coefficient at a resetting time is erased, so that the whole process is a fully renewal one. We show that the properties of the probability densities in such processes (erasing or retaining the memory on the diffusion coefficient) are vastly different.
For both Lévy flight and Lévy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For Lévy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the Lévy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.
The passive and active motion of micron-sized tracer particles in crowded liquids and inside living biological cells is ubiquitously characterised by 'viscoelastic' anomalous diffusion, in which the increments of the motion feature long-ranged negative and positive correlations. While viscoelastic anomalous diffusion is typically modelled by a Gaussian process with correlated increments, so-called fractional Gaussian noise, an increasing number of systems are reported, in which viscoelastic anomalous diffusion is paired with non-Gaussian displacement distributions. Following recent advances in Brownian yet non-Gaussian diffusion we here introduce and discuss several possible versions of random-diffusivity models with long-ranged correlations. While all these models show a crossover from non-Gaussian to Gaussian distributions beyond some correlation time, their mean squared displacements exhibit strikingly different behaviours: depending on the model crossovers from anomalous to normal diffusion are observed, as well as a priori unexpected dependencies of the effective diffusion coefficient on the correlation exponent. Our observations of the non-universality of random-diffusivity viscoelastic anomalous diffusion are important for the analysis of experiments and a better understanding of the physical origins of 'viscoelastic yet non-Gaussian' diffusion.
The passive and active motion of micron-sized tracer particles in crowded liquids and inside living biological cells is ubiquitously characterised by 'viscoelastic' anomalous diffusion, in which the increments of the motion feature long-ranged negative and positive correlations. While viscoelastic anomalous diffusion is typically modelled by a Gaussian process with correlated increments, so-called fractional Gaussian noise, an increasing number of systems are reported, in which viscoelastic anomalous diffusion is paired with non-Gaussian displacement distributions. Following recent advances in Brownian yet non-Gaussian diffusion we here introduce and discuss several possible versions of random-diffusivity models with long-ranged correlations. While all these models show a crossover from non-Gaussian to Gaussian distributions beyond some correlation time, their mean squared displacements exhibit strikingly different behaviours: depending on the model crossovers from anomalous to normal diffusion are observed, as well as a priori unexpected dependencies of the effective diffusion coefficient on the correlation exponent. Our observations of the non-universality of random-diffusivity viscoelastic anomalous diffusion are important for the analysis of experiments and a better understanding of the physical origins of 'viscoelastic yet non-Gaussian' diffusion.
Lévy flights are paradigmatic generalised random walk processes, in which the independent stationary increments—the 'jump lengths'—are drawn from an -stable jump length distribution with long-tailed, power-law asymptote. As a result, the variance of Lévy flights diverges and the trajectory is characterised by occasional extremely long jumps. Such long jumps significantly decrease the probability to revisit previous points of visitation, rendering Lévy flights efficient search processes in one and two dimensions. To further quantify their precise property as random search strategies we here study the first-passage time properties of Lévy flights in one-dimensional semi-infinite and bounded domains for symmetric and asymmetric jump length distributions. To obtain the full probability density function of first-passage times for these cases we employ two complementary methods. One approach is based on the space-fractional diffusion equation for the probability density function, from which the survival probability is obtained for different values of the stable index and the skewness (asymmetry) parameter . The other approach is based on the stochastic Langevin equation with -stable driving noise. Both methods have their advantages and disadvantages for explicit calculations and numerical evaluation, and the complementary approach involving both methods will be profitable for concrete applications. We also make use of the Skorokhod theorem for processes with independent increments and demonstrate that the numerical results are in good agreement with the analytical expressions for the probability density function of the first-passage times.
The Ornstein–Uhlenbeck process is a stationary and ergodic Gaussian process, that is fully determined by its covariance function and mean. We show here that the generic definitions of the ensemble- and time-averaged mean squared displacements fail to capture these properties consistently, leading to a spurious ergodicity breaking. We propose to remedy this failure by redefining the mean squared displacements such that they reflect unambiguously the statistical properties of any stochastic process. In particular we study the effect of the initial condition in the Ornstein–Uhlenbeck process and its fractional extension. For the fractional Ornstein–Uhlenbeck process representing typical experimental situations in crowded environments such as living biological cells, we show that the stationarity of the process delicately depends on the initial condition.
The Ornstein–Uhlenbeck process is a stationary and ergodic Gaussian process, that is fully determined by its covariance function and mean. We show here that the generic definitions of the ensemble- and time-averaged mean squared displacements fail to capture these properties consistently, leading to a spurious ergodicity breaking. We propose to remedy this failure by redefining the mean squared displacements such that they reflect unambiguously the statistical properties of any stochastic process. In particular we study the effect of the initial condition in the Ornstein–Uhlenbeck process and its fractional extension. For the fractional Ornstein–Uhlenbeck process representing typical experimental situations in crowded environments such as living biological cells, we show that the stationarity of the process delicately depends on the initial condition.
A growing number of biological, soft, and active matter systems are observed to exhibit normal diffusive dynamics with a linear growth of the mean-squared displacement, yet with a non-Gaussian distribution of increments. Based on the Chubinsky-Slater idea of a diffusing diffusivity, we here establish and analyze a minimal model framework of diffusion processes with fluctuating diffusivity. In particular, we demonstrate the equivalence of the diffusing diffusivity process with a superstatistical approach with a distribution of diffusivities, at times shorter than the diffusivity correlation time. At longer times, a crossover to a Gaussian distribution with an effective diffusivity emerges. Specifically, we establish a subordination picture of Brownian but non-Gaussian diffusion processes, which can be used for a wide class of diffusivity fluctuation statistics. Our results are shown to be in excellent agreement with simulations and numerical evaluations.
We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black-Scholes-Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics.
We discuss generalized integro-differential diffusion equations whose integral kernels are not of a simple power law form, and thus these equations themselves do not belong to the family of fractional diffusion equations exhibiting a monoscaling behavior. They instead generate a broad class of anomalous nonscaling patterns, which correspond either to crossovers between different power laws, or to a non-power-law behavior as exemplified by the logarithmic growth of the width of the distribution. We consider normal and modified forms of these generalized diffusion equations and provide a brief discussion of three generic types of integral kernels for each form, namely, distributed order, truncated power law and truncated distributed order kernels. For each of the cases considered we prove the non-negativity of the solution of the corresponding generalized diffusion equation and calculate the mean squared displacement. (C) 2017 Elsevier Ltd. All rights reserved.
We address the generic problem of random search for a point-like target on a line. Using the measures of search reliability and efficiency to quantify the random search quality, we compare Brownian search with Levy search based on long-tailed jump length distributions. We then compare these results with a search process combined of two different long-tailed jump length distributions. Moreover, we study the case of multiple targets located by a Levy searcher.
We present an effective dynamical model for the onset of bacterial bioluminescence, one of the most studied quorum sensing-mediated traits. Our model is built upon simple equations that describe the growth of the bacterial colony, the production and accumulation of autoinducer signal molecules, their sensing within bacterial cells, and the ensuing quorum activation mechanism that triggers bioluminescent emission. The model is directly tested to quantitatively reproduce the experimental distributions of photon emission times, previously measured for bacterial colonies of Vibrio jasicida, a luminescent bacterium belonging to the Harveyi clade, growing in a highly drying environment. A distinctive and novel feature of the proposed model is bioluminescence ‘quenching’ after a given time elapsed from activation. Using an advanced fitting procedure based on the simulated annealing algorithm, we are able to infer from the experimental observations the biochemical parameters used in the model. Such parameters are in good agreement with the literature data. As a further result, we find that, at least in our experimental conditions, light emission in bioluminescent bacteria appears to originate from a subtle balance between colony growth and quorum activation due to autoinducers diffusion, with the two phenomena occurring on the same time scale. This finding is consistent with a negative feedback mechanism previously reported for Vibrio harveyi.
In ageing systems physical observables explicitly depend on the time span elapsing between the original initiation of the system and the actual start of the recording of the particle motion. We here study the signatures of ageing in the framework of ultraslow continuous time random walk processes with super-heavy tailed waiting time densities. We derive the density for the forward or recurrent waiting time of the motion as function of the ageing time, generalise the Montroll-Weiss equation for this process, and analyse the ageing behaviour of the ensemble and time averaged mean squared displacements.
We consider a generalized diffusion equation in two dimensions for modeling diffusion on a comb-like structures. We analyze the probability distribution functions and we derive the mean squared displacement in x and y directions. Different forms of the memory kernels (Dirac delta, power-law, and distributed order) are considered. It is shown that anomalous diffusion may occur along both x and y directions. Ultraslow diffusion and some more general diffusive processes are observed as well. We give the corresponding continuous time random walk model for the considered two dimensional diffusion-like equation on a comb, and we derive the probability distribution functions which subordinate the process governed by this equation to the Wiener process.
Underdamped scaled Brownian motion: (non-)existence of the overdamped limit in anomalous diffusion
(2016)
It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases.
A combined dynamics consisting of Brownian motion and Levy flights is exhibited by a variety of biological systems performing search processes. Assessing the search reliability of ever locating the target and the search efficiency of doing so economically of such dynamics thus poses an important problem. Here we model this dynamics by a one-dimensional fractional Fokker-Planck equation combining unbiased Brownian motion and Levy flights. By solving this equation both analytically and numerically we show that the superposition of recurrent Brownian motion and Levy flights with stable exponent alpha < 1, by itself implying zero probability of hitting a point on a line, leads to transient motion with finite probability of hitting any point on the line. We present results for the exact dependence of the values of both the search reliability and the search efficiency on the distance between the starting and target positions as well as the choice of the scaling exponent a of the Levy flight component.
For both Lévy flight and Lévy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For Lévy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the Lévy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.
For both Lévy flight and Lévy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For Lévy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the Lévy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.
In various biological systems and small scale technological applications particles transiently bind to a cylindrical surface. Upon unbinding the particles diffuse in the vicinal bulk before rebinding to the surface. Such bulk-mediated excursions give rise to an effective surface translation, for which we here derive and discuss the dynamic equations, including additional surface diffusion. We discuss the time evolution of the number of surface-bound particles, the effective surface mean squared displacement, and the surface propagator. In particular, we observe sub- and superdiffusive regimes. A plateau of the surface mean-squared displacement reflects a stalling of the surface diffusion at longer times. Finally, the corresponding first passage problem for the cylindrical geometry is analysed.
In this paper we propose an algorithm to distinguish between light- and heavy-tailed probability laws underlying random datasets. The idea of the algorithm, which is visual and easy to implement, is to check whether the underlying law belongs to the domain of attraction of the Gaussian or non-Gaussian stable distribution by examining its rate of convergence. The method allows to discriminate between stable and various non-stable distributions. The test allows to differentiate between distributions, which appear the same according to standard Kolmogorov-Smirnov test. In particular, it helps to distinguish between stable and Student's t probability laws as well as between the stable and tempered stable, the cases which are considered in the literature as very cumbersome. Finally, we illustrate the procedure on plasma data to identify cases with so-called L-H transition.
Abstract
The emerging diffusive dynamics in many complex systems show a characteristic crossover behaviour from anomalous to normal diffusion which is otherwise fitted by two independent power-laws. A prominent example for a subdiffusive–diffusive crossover are viscoelastic systems such as lipid bilayer membranes, while superdiffusive–diffusive crossovers occur in systems of actively moving biological cells. We here consider the general dynamics of a stochastic particle driven by so-called tempered fractional Gaussian noise, that is noise with Gaussian amplitude and power-law correlations, which are cut off at some mesoscopic time scale. Concretely we consider such noise with built-in exponential or power-law tempering, driving an overdamped Langevin equation (fractional Brownian motion) and fractional Langevin equation motion. We derive explicit expressions for the mean squared displacement and correlation functions, including different shapes of the crossover behaviour depending on the concrete tempering, and discuss the physical meaning of the tempering. In the case of power-law tempering we also find a crossover behaviour from faster to slower superdiffusion and slower to faster subdiffusion. As a direct application of our model we demonstrate that the obtained dynamics quantitatively describes the subdiffusion–diffusion and subdiffusion–subdiffusion crossover in lipid bilayer systems. We also show that a model of tempered fractional Brownian motion recently proposed by Sabzikar and Meerschaert leads to physically very different behaviour with a seemingly paradoxical ballistic long time scaling.
Abstract
The emerging diffusive dynamics in many complex systems show a characteristic crossover behaviour from anomalous to normal diffusion which is otherwise fitted by two independent power-laws. A prominent example for a subdiffusive–diffusive crossover are viscoelastic systems such as lipid bilayer membranes, while superdiffusive–diffusive crossovers occur in systems of actively moving biological cells. We here consider the general dynamics of a stochastic particle driven by so-called tempered fractional Gaussian noise, that is noise with Gaussian amplitude and power-law correlations, which are cut off at some mesoscopic time scale. Concretely we consider such noise with built-in exponential or power-law tempering, driving an overdamped Langevin equation (fractional Brownian motion) and fractional Langevin equation motion. We derive explicit expressions for the mean squared displacement and correlation functions, including different shapes of the crossover behaviour depending on the concrete tempering, and discuss the physical meaning of the tempering. In the case of power-law tempering we also find a crossover behaviour from faster to slower superdiffusion and slower to faster subdiffusion. As a direct application of our model we demonstrate that the obtained dynamics quantitatively describes the subdiffusion–diffusion and subdiffusion–subdiffusion crossover in lipid bilayer systems. We also show that a model of tempered fractional Brownian motion recently proposed by Sabzikar and Meerschaert leads to physically very different behaviour with a seemingly paradoxical ballistic long time scaling.
A considerable number of systems have recently been reported in which
Brownian yet non-Gaussian dynamics was observed. These are processes characterised by a linear growth in time of the mean squared displacement, yet the probability density function of the particle displacement is distinctly non-Gaussian, and often of exponential(Laplace) shape. This apparently ubiquitous behaviour observed in very different physical systems has been interpreted as resulting from diffusion in inhomogeneous environments and mathematically represented through a variable, stochastic diffusion coefficient. Indeed different models describing a fluctuating diffusivity have been studied. Here we present a new view of the stochastic basis describing time dependent random diffusivities within a broad spectrum of distributions. Concretely, our study is based on the very generic class of the generalised Gamma distribution. Two models for the particle spreading in such random diffusivity settings are studied. The first belongs to the class of generalised grey Brownian motion while the second follows from the idea of diffusing diffusivities. The two processes exhibit significant characteristics which reproduce experimental results from different biological and physical systems. We promote these two physical models for the description of stochastic particle motion in complex environments.
Brownian yet non-Gaussian dynamics was observed. These are processes characterised by a linear growth in time of the mean squared displacement, yet the probability density function of the particle displacement is distinctly non-Gaussian, and often of exponential(Laplace) shape. This apparently ubiquitous behaviour observed in very different physical systems has been interpreted as resulting from diffusion in inhomogeneous environments and mathematically represented through a variable, stochastic diffusion coefficient. Indeed different models describing a fluctuating diffusivity have been studied. Here we present a new view of the stochastic basis describing time dependent random diffusivities within a broad spectrum of distributions. Concretely, our study is based on the very generic class of the generalised Gamma distribution. Two models for the particle spreading in such random diffusivity settings are studied. The first belongs to the class of generalised grey Brownian motion while the second follows from the idea of diffusing diffusivities. The two processes exhibit significant characteristics which reproduce experimental results from different biological and physical systems. We promote these two physical models for the description of stochastic particle motion in complex environments.
We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black–Scholes–Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics.
We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black–Scholes–Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics.
Brownian motion is ergodic in the Boltzmann-Khinchin sense that long time averages of physical observables such as the mean squared displacement provide the same information as the corresponding ensemble average, even at out-of-equilibrium conditions. This property is the fundamental prerequisite for single particle tracking and its analysis in simple liquids. We study analytically and by event-driven molecular dynamics simulations the dynamics of force-free cooling granular gases and reveal a violation of ergodicity in this Boltzmann-Khinchin sense as well as distinct ageing of the system. Such granular gases comprise materials such as dilute gases of stones, sand, various types of powders, or large molecules, and their mixtures are ubiquitous in Nature and technology, in particular in Space. We treat-depending on the physical-chemical properties of the inter-particle interaction upon their pair collisions-both a constant and a velocity-dependent (viscoelastic) restitution coefficient epsilon. Moreover we compare the granular gas dynamics with an effective single particle stochastic model based on an underdamped Langevin equation with time dependent diffusivity. We find that both models share the same behaviour of the ensemble mean squared displacement (MSD) and the velocity correlations in the limit of weak dissipation. Qualitatively, the reported non-ergodic behaviour is generic for granular gases with any realistic dependence of epsilon on the impact velocity of particles.
Aging scaled Brownian motion
(2015)
Scaled Brownian motion (SBM) is widely used to model anomalous diffusion of passive tracers in complex and biological systems. It is a highly nonstationary process governed by the Langevin equation for Brownian motion, however, with a power-law time dependence of the noise strength. Here we study the aging properties of SBM for both unconfined and confined motion. Specifically, we derive the ensemble and time averaged mean squared displacements and analyze their behavior in the regimes of weak, intermediate, and strong aging. A very rich behavior is revealed for confined aging SBM depending on different aging times and whether the process is sub- or superdiffusive. We demonstrate that the information on the aging factorizes with respect to the lag time and exhibits a functional form that is identical to the aging behavior of scale-free continuous time random walk processes. While SBM exhibits a disparity between ensemble and time averaged observables and is thus weakly nonergodic, strong aging is shown to effect a convergence of the ensemble and time averaged mean squared displacement. Finally, we derive the density of first passage times in the semi-infinite domain that features a crossover defined by the aging time.
We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form D(t) similar or equal to 1/t. For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In a harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations.
Fluctuation relations for anomalous dynamics generated by time-fractional Fokker-Planck equations
(2015)
Anomalous dynamics characterized by non-Gaussian probability distributions (PDFs) and/or temporal long-range correlations can cause subtle modifications of conventional fluctuation relations (FRs). As prototypes we study three variants of a generic time-fractional Fokker-Planck equation with constant force. Type A generates superdiffusion, type B subdiffusion and type C both super-and subdiffusion depending on parameter variation. Furthermore type C obeys a fluctuation-dissipation relation whereas A and B do not. We calculate analytically the position PDFs for all three cases and explore numerically their strongly non-Gaussian shapes. While for type C we obtain the conventional transient work FR, type A and type B both yield deviations by featuring a coefficient that depends on time and by a nonlinear dependence on the work. We discuss possible applications of these types of dynamics and FRs to experiments.