@article{ThapaWyłomańskaSikoraetal.2021, author = {Thapa, Samudrajit and Wyłomańska, Agnieszka and Sikora, Grzegorz and Wagner, Caroline E. and Krapf, Diego and Kantz, Holger and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories}, series = {New Journal of Physics}, volume = {23}, journal = {New Journal of Physics}, publisher = {Dt. Physikalische Ges. ; IOP}, address = {Bad Honnef ; London}, issn = {1367-2630}, doi = {10.1088/1367-2630/abd50e}, pages = {22}, year = {2021}, abstract = {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.}, language = {en} } @article{SposiniChechkinMetzler2018, author = {Sposini, Vittoria and Chechkin, Aleksei V. and Metzler, Ralf}, title = {First passage statistics for diffusing diffusivity}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {52}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {4}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aaf6ff}, pages = {11}, year = {2018}, abstract = {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.}, language = {en} } @article{DybiecCapalaChechkinetal.2018, author = {Dybiec, Bartlomiej and Capala, Karol and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Conservative random walks in confining potentials}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {52}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aaefc2}, pages = {25}, year = {2018}, abstract = {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.}, language = {en} } @article{MardoukhiChechkinMetzler2020, author = {Mardoukhi, Yousof and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Spurious ergodicity breaking in normal and fractional Ornstein-Uhlenbeck process}, series = {New Journal of Physics}, volume = {22}, journal = {New Journal of Physics}, publisher = {IOP}, address = {London}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab950b}, pages = {18}, year = {2020}, abstract = {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.}, language = {en} } @article{PalyulinMantsevichKlagesetal.2017, author = {Palyulin, Vladimir V. and Mantsevich, Vladimir N. and Klages, Rainer and Metzler, Ralf and Chechkin, Aleksei V.}, title = {Comparison of pure and combined search strategies for single and multiple targets}, series = {The European physical journal : B, Condensed matter and complex systems}, volume = {90}, journal = {The European physical journal : B, Condensed matter and complex systems}, publisher = {Springer}, address = {New York}, issn = {1434-6028}, doi = {10.1140/epjb/e2017-80372-4}, pages = {20 -- 37}, year = {2017}, abstract = {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.}, language = {en} } @article{PalyulinBlackburnLomholtetal.2019, author = {Palyulin, Vladimir V. and Blackburn, George and Lomholt, Michael A. and Watkins, Nicholas W. and Metzler, Ralf and Klages, Rainer and Chechkin, Aleksei V.}, title = {First passage and first hitting times of Levy flights and Levy walks}, series = {New journal of physics : the open-access journal for physics}, volume = {21}, journal = {New journal of physics : the open-access journal for physics}, number = {10}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab41bb}, pages = {23}, year = {2019}, abstract = {For both L{\´e}vy flight and L{\´e}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{\´e}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{\´e}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.}, language = {en} } @article{CherstvyThapaMardoukhietal.2018, author = {Cherstvy, Andrey G. and Thapa, Samudrajit and Mardoukhi, Yousof and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Time averages and their statistical variation for the Ornstein-Uhlenbeck process}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {98}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {2}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.98.022134}, pages = {15}, year = {2018}, abstract = {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.}, language = {en} } @misc{PalyulinBlackburnLomholtetal.2019, author = {Palyulin, Vladimir V and Blackburn, George and Lomholt, Michael A and Watkins, Nicholas W and Metzler, Ralf and Klages, Rainer and Chechkin, Aleksei V.}, title = {First passage and first hitting times of L{\´e}vy flights and L{\´e}vy walks}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {785}, issn = {1866-8372}, doi = {10.25932/publishup-43983}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-439832}, pages = {25}, year = {2019}, abstract = {For both L{\´e}vy flight and L{\´e}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{\´e}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{\´e}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.}, language = {en} } @article{CherstvyChechkinMetzler2014, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity}, series = {Soft matter}, volume = {2014}, journal = {Soft matter}, number = {10}, publisher = {Royal Society of Chemistry}, issn = {2046-2069}, doi = {10.1039/c3sm52846d}, pages = {1591 -- 1601}, year = {2014}, abstract = {We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells.}, language = {en} } @misc{CherstvyChechkinMetzler2014, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity}, number = {168}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-74021}, pages = {1591 -- 1601}, year = {2014}, abstract = {We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells.}, language = {en} } @article{JeonChechkinMetzler2014, author = {Jeon, Jae-Hyung and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Scaled Brownian motion: a paradoxical process with a time dependent diffusivity for the description of anomalous diffusion}, series = {Physical chemistry, chemical physics : PCCP}, volume = {30}, journal = {Physical chemistry, chemical physics : PCCP}, number = {16}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, doi = {10.1039/C4CP02019G}, pages = {15811 -- 15817}, year = {2014}, abstract = {Anomalous diffusion is frequently described by scaled Brownian motion (SBM){,} a Gaussian process with a power-law time dependent diffusion coefficient. Its mean squared displacement is ?x2(t)? [similar{,} equals] 2K(t)t with K(t) [similar{,} equals] t[small alpha]-1 for 0 < [small alpha] < 2. SBM may provide a seemingly adequate description in the case of unbounded diffusion{,} for which its probability density function coincides with that of fractional Brownian motion. Here we show that free SBM is weakly non-ergodic but does not exhibit a significant amplitude scatter of the time averaged mean squared displacement. More severely{,} we demonstrate that under confinement{,} the dynamics encoded by SBM is fundamentally different from both fractional Brownian motion and continuous time random walks. SBM is highly non-stationary and cannot provide a physical description for particles in a thermalised stationary system. Our findings have direct impact on the modelling of single particle tracking experiments{,} in particular{,} under confinement inside cellular compartments or when optical tweezers tracking methods are used.}, language = {en} } @misc{MardoukhiChechkinMetzler2020, author = {Mardoukhi, Yousof and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Spurious ergodicity breaking in normal and fractional Ornstein-Uhlenbeck process}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {981}, issn = {1866-8372}, doi = {10.25932/publishup-47487}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-474875}, pages = {20}, year = {2020}, abstract = {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.}, language = {en} } @article{MardoukhiJeonChechkinetal.2018, author = {Mardoukhi, Yousof and Jeon, Jae-Hyung and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Fluctuations of random walks in critical random environments}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {20}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {31}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c8cp03212b}, pages = {20427 -- 20438}, year = {2018}, abstract = {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.}, language = {en} } @article{PadashChechkinDybiecetal.2019, author = {Padash, Amin and Chechkin, Aleksei V. and Dybiec, Bartlomiej and Pavlyukevich, Ilya and Shokri, Babak and Metzler, Ralf}, title = {First-passage properties of asymmetric Levy flights}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {52}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {45}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/ab493e}, pages = {48}, year = {2019}, abstract = {L{\´e}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{\´e}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{\´e}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{\´e}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.}, language = {en} } @article{SandevMetzlerChechkin2018, author = {Sandev, Trifce and Metzler, Ralf and Chechkin, Aleksei V.}, title = {From continuous time random walks to the generalized diffusion equation}, series = {Fractional calculus and applied analysis : an international journal for theory and applications}, volume = {21}, journal = {Fractional calculus and applied analysis : an international journal for theory and applications}, number = {1}, publisher = {De Gruyter}, address = {Berlin}, issn = {1311-0454}, doi = {10.1515/fca-2018-0002}, pages = {10 -- 28}, year = {2018}, abstract = {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.}, language = {en} } @article{WangCherstvyChechkinetal.2020, author = {Wang, Wei and Cherstvy, Andrey G. and Chechkin, Aleksei V. and Thapa, Samudrajit and Seno, Flavio and Liu, Xianbin and Metzler, Ralf}, title = {Fractional Brownian motion with random diffusivity}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {53}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {47}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aba467}, pages = {34}, year = {2020}, abstract = {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.}, language = {en} } @article{SchulzChechkinMetzler2013, author = {Schulz, Johannes H. P. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Correlated continuous time random walks - combining scale-invariance with long-range memory for spatial and temporal dynamics}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {46}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {47}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/46/47/475001}, pages = {22}, year = {2013}, abstract = {Standard continuous time random walk (CTRW) models are renewal processes in the sense that at each jump a new, independent pair of jump length and waiting time are chosen. Globally, anomalous diffusion emerges through scale-free forms of the jump length and/or waiting time distributions by virtue of the generalized central limit theorem. Here we present a modified version of recently proposed correlated CTRW processes, where we incorporate a power-law correlated noise on the level of both jump length and waiting time dynamics. We obtain a very general stochastic model, that encompasses key features of several paradigmatic models of anomalous diffusion: discontinuous, scale-free displacements as in Levy flights, scale-free waiting times as in subdiffusive CTRWs, and the long-range temporal correlations of fractional Brownian motion (FBM). We derive the exact solutions for the single-time probability density functions and extract the scaling behaviours. Interestingly, we find that different combinations of the model parameters lead to indistinguishable shapes of the emerging probability density functions and identical scaling laws. Our model will be useful for describing recent experimental single particle tracking data that feature a combination of CTRW and FBM properties.}, language = {en} } @article{CherstvyChechkinMetzler2013, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Anomalous diffusion and ergodicity breaking in heterogeneous diffusion processes}, series = {New journal of physics : the open-access journal for physics}, volume = {15}, journal = {New journal of physics : the open-access journal for physics}, number = {15}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/15/8/083039}, pages = {13}, year = {2013}, abstract = {We demonstrate the non-ergodicity of a simple Markovian stochastic process with space-dependent diffusion coefficient D(x). For power-law forms D(x) similar or equal to vertical bar x vertical bar(alpha), this process yields anomalous diffusion of the form < x(2)(t)> similar or equal to t(2/(2-alpha)). Interestingly, in both the sub- and superdiffusive regimes we observe weak ergodicity breaking: the scaling of the time-averaged mean-squared displacement <(delta(2)(Delta))over bar> remains linear in the lag time Delta and thus differs from the corresponding ensemble average < x(2)(t)>. We analyse the non-ergodic behaviour of this process in terms of the time-averaged mean- squared displacement (delta(2)) over bar and its random features, i.e. the statistical distribution of (delta(2)) over bar and the ergodicity breaking parameters. The heterogeneous diffusion model represents an alternative approach to non- ergodic, anomalous diffusion that might be particularly relevant for diffusion in heterogeneous media.}, language = {en} } @article{JeonChechkinMetzler2014, author = {Jeon, Jae-Hyung and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Scaled Brownian motion: a paradoxical process with a time dependent diffusivity for the description of anomalous diffusion}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {16}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {30}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c4cp02019g}, pages = {15811 -- 15817}, year = {2014}, abstract = {Anomalous diffusion is frequently described by scaled Brownian motion (SBM), a Gaussian process with a power-law time dependent diffusion coefficient. Its mean squared displacement is < x(2)(t) similar or equal to 2K(t)t with K(t) similar or equal to t(alpha-1) for 0 < alpha < 2. SBM may provide a seemingly adequate description in the case of unbounded diffusion, for which its probability density function coincides with that of fractional Brownian motion. Here we show that free SBM is weakly non-ergodic but does not exhibit a significant amplitude scatter of the time averaged mean squared displacement. More severely, we demonstrate that under confinement, the dynamics encoded by SBM is fundamentally different from both fractional Brownian motion and continuous time random walks. SBM is highly non-stationary and cannot provide a physical description for particles in a thermalised stationary system. Our findings have direct impact on the modelling of single particle tracking experiments, in particular, under confinement inside cellular compartments or when optical tweezers tracking methods are used.}, language = {en} } @article{PalyulinChechkinMetzler2014, author = {Palyulin, Vladimir V. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Space-fractional Fokker-Planck equation and optimization of random search processes in the presence of an external bias}, series = {Journal of statistical mechanics: theory and experiment}, journal = {Journal of statistical mechanics: theory and experiment}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1742-5468}, doi = {10.1088/1742-5468/2014/11/P11031}, pages = {32}, year = {2014}, abstract = {Based on the space-fractional Fokker-Planck equation with a delta-sink term, we study the efficiency of random search processes based on Levy flights with power-law distributed jump lengths in the presence of an external drift, for instance, an underwater current, an airflow, or simply the preference of the searcher based on prior experience. While Levy flights turn out to be efficient search processes when the target is upstream relative to the starting point, in the downstream scenario, regular Brownian motion turns out to be advantageous. This is caused by the occurrence of leapovers of Levy flights, due to which Levy flights typically overshoot a point or small interval. Studying the solution of the fractional Fokker-Planck equation, we establish criteria when the combination of the external stream and the initial distance between the starting point and the target favours Levy flights over the regular Brownian search. Contrary to the common belief that Levy flights with a Levy index alpha = 1 (i.e. Cauchy flights) are optimal for sparse targets, we find that the optimal value for alpha may range in the entire interval (1, 2) and explicitly include Brownian motion as the most efficient search strategy overall.}, language = {en} } @article{NezhadhaghighiChechkinMetzler2014, author = {Nezhadhaghighi, M. Ghasemi and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Numerical approach to unbiased and driven generalized elastic model}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {140}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {2}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.4858425}, pages = {9}, year = {2014}, abstract = {From scaling arguments and numerical simulations, we investigate the properties of the generalized elastic model (GEM) that is used to describe various physical systems such as polymers, membranes, single-file systems, or rough interfaces. We compare analytical and numerical results for the subdiffusion exponent beta characterizing the growth of the mean squared displacement <(delta h)(2)> of the field h described by the GEM dynamic equation. We study the scaling properties of the qth order moments with time, finding that the interface fluctuations show no intermittent behavior. We also investigate the ergodic properties of the process h in terms of the ergodicity breaking parameter and the distribution of the time averaged mean squared displacement. Finally, we study numerically the driven GEM with a constant, localized perturbation and extract the characteristics of the average drift for a tagged probe.}, language = {en} } @article{BodrovaChechkinCherstvyetal.2015, author = {Bodrova, Anna S. and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Ultraslow scaled Brownian motion}, series = {New journal of physics : the open-access journal for physics}, volume = {17}, journal = {New journal of physics : the open-access journal for physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/17/6/063038}, pages = {16}, year = {2015}, abstract = {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.}, language = {en} } @article{PalyulinChechkinMetzler2014, author = {Palyulin, Vladimir V. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Levy flights do not always optimize random blind search for sparse targets}, series = {Proceedings of the National Academy of Sciences of the United States of America}, volume = {111}, journal = {Proceedings of the National Academy of Sciences of the United States of America}, number = {8}, publisher = {National Acad. of Sciences}, address = {Washington}, issn = {0027-8424}, doi = {10.1073/pnas.1320424111}, pages = {2931 -- 2936}, year = {2014}, abstract = {It is generally believed that random search processes based on scale-free, Levy stable jump length distributions (Levy flights) optimize the search for sparse targets. Here we show that this popular search advantage is less universal than commonly assumed. We study the efficiency of a minimalist search model based on Levy flights in the absence and presence of an external drift (underwater current, atmospheric wind, a preference of the walker owing to prior experience, or a general bias in an abstract search space) based on two different optimization criteria with respect to minimal search time and search reliability (cumulative arrival probability). Although Levy flights turn out to be efficient search processes when the target is far from the starting point, or when relative to the starting point the target is upstream, we show that for close targets and for downstream target positioning regular Brownian motion turns out to be the advantageous search strategy. Contrary to claims that Levy flights with a critical exponent alpha = 1 are optimal for the search of sparse targets in different settings, based on our optimization parameters the optimal a may range in the entire interval (1, 2) and especially include Brownian motion as the overall most efficient search strategy.}, language = {en} } @article{SandevChechkinKantzetal.2015, author = {Sandev, Trifce and Chechkin, Aleksei V. and Kantz, Holger and Metzler, Ralf}, title = {Diffusion and fokker-planck-smoluchowski equations with generalized memory kernel}, series = {Fractional calculus and applied analysis : an international journal for theory and applications}, volume = {18}, journal = {Fractional calculus and applied analysis : an international journal for theory and applications}, number = {4}, publisher = {De Gruyter}, address = {Berlin}, issn = {1311-0454}, doi = {10.1515/fca-2015-0059}, pages = {1006 -- 1038}, year = {2015}, abstract = {We consider anomalous stochastic processes based on the renewal continuous time random walk model with different forms for the probability density of waiting times between individual jumps. In the corresponding continuum limit we derive the generalized diffusion and Fokker-Planck-Smoluchowski equations with the corresponding memory kernels. We calculate the qth order moments in the unbiased and biased cases, and demonstrate that the generalized Einstein relation for the considered dynamics remains valid. The relaxation of modes in the case of an external harmonic potential and the convergence of the mean squared displacement to the thermal plateau are analyzed.}, language = {en} } @article{GodecChechkinBarkaietal.2014, author = {Godec, Aljaz and Chechkin, Aleksei V. and Barkai, Eli and Kantz, Holger and Metzler, Ralf}, title = {Localisation and universal fluctuations in ultraslow diffusion processes}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {47}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {49}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/47/49/492002}, pages = {10}, year = {2014}, abstract = {We study ultraslow diffusion processes with logarithmic mean squared displacement (MSD) < x(2)(t)> similar or equal to log(gamma)t. Comparison of annealed (renewal) continuous time random walks (CTRWs) with logarithmic waiting time distribution psi(tau) similar or equal to 1/(tau log(1+gamma)tau) and Sinai diffusion in quenched random landscapes reveals striking similarities, despite the great differences in their physical nature. In particular, they exhibit a weakly non-ergodic disparity of the time-averaged and ensemble-averaged MSDs. Remarkably, for the CTRW we observe that the fluctuations of time averages become universal, with an exponential suppression of mobile trajectories. We discuss the fundamental connection between the Golosov localization effect and non-ergodicity in the sense of the disparity between ensemble-averaged MSD and time-averaged MSD.}, language = {en} } @article{CherstvyChechkinMetzler2014, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Ageing and confinement in non-ergodic heterogeneous diffusion processes}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {47}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {48}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/47/48/485002}, pages = {18}, year = {2014}, abstract = {We study the effects of ageing-the time delay between initiation of the physical process at t = 0 and start of observation at some time t(a) > 0-and spatial confinement on the properties of heterogeneous diffusion processes (HDPs) with deterministic power-law space-dependent diffusivities, D(x) = D-0 vertical bar x vertical bar(alpha). From analysis of the ensemble and time averaged mean squared displacements and the ergodicity breaking parameter quantifying the inherent degree of irreproducibility of individual realizations of the HDP we obtain striking similarities to ageing subdiffusive continuous time random walks with scale-free waiting time distributions. We also explore how both processes can be distinguished. For confined HDPs we study the long-time saturation of the ensemble and time averaged particle displacements as well as the magnitude of the inherent scatter of time averaged displacements and contrast the outcomes to the results known for other anomalous diffusion processes under confinement.}, language = {en} } @article{CherstvyChechkinMetzler2014, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity}, series = {Soft matter}, volume = {10}, journal = {Soft matter}, number = {10}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, doi = {10.1039/c3sm52846d}, pages = {1591 -- 1601}, year = {2014}, abstract = {We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells.}, language = {en} } @article{SposiniChechkinSenoetal.2018, author = {Sposini, Vittoria and Chechkin, Aleksei V. and Seno, Flavio and Pagnini, Gianni and Metzler, Ralf}, title = {Random diffusivity from stochastic equations}, series = {New Journal of Physics}, journal = {New Journal of Physics}, publisher = {Deutsche Physikalische Gesellschaft / Institute of Physics}, address = {Bad Honnef und London}, issn = {1367-2630}, doi = {10.1088/1367-2630/aab696}, pages = {1 -- 33}, year = {2018}, abstract = {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.}, language = {en} } @misc{SposiniChechkinFlavioetal.2018, author = {Sposini, Vittoria and Chechkin, Aleksei V. and Flavio, Seno and Pagnini, Gianni and Metzler, Ralf}, title = {Random diffusivity from stochastic equations}, series = {New Journal of Physics}, journal = {New Journal of Physics}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-409743}, pages = {33}, year = {2018}, abstract = {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.}, language = {en} } @article{PalyulinBlackburnLomholtetal.2019, author = {Palyulin, Vladimir V and Blackburn, George and Lomholt, Michael A and Watkins, Nicholas W and Metzler, Ralf and Klages, Rainer and Chechkin, Aleksei V.}, title = {First passage and first hitting times of L{\´e}vy flights and L{\´e}vy walks}, series = {New Journal of Physics}, volume = {21}, journal = {New Journal of Physics}, publisher = {Dt. Physikalische Ges.}, address = {Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab41bb}, pages = {24}, year = {2019}, abstract = {For both L{\´e}vy flight and L{\´e}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{\´e}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{\´e}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.}, language = {en} } @article{PalyulinChechkinKlagesetal.2016, author = {Palyulin, Vladimir V. and Chechkin, Aleksei V. and Klages, Rainer and Metzler, Ralf}, title = {Search reliability and search efficiency of combined Levy-Brownian motion: long relocations mingled with thorough local exploration}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {49}, journal = {Journal of physics : A, Mathematical and theoretical}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/49/39/394002}, pages = {2189 -- 2193}, year = {2016}, abstract = {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.}, language = {en} } @misc{JeonChechkinMetzler2014, author = {Jeon, Jae-Hyung and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Scaled Brownian motion: a paradoxical process with a time dependent diffusivity for the description of anomalous diffusion}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-76302}, pages = {15811 -- 15817}, year = {2014}, abstract = {Anomalous diffusion is frequently described by scaled Brownian motion (SBM){,} a Gaussian process with a power-law time dependent diffusion coefficient. Its mean squared displacement is ?x2(t)? [similar{,} equals] 2K(t)t with K(t) [similar{,} equals] t[small alpha]-1 for 0 < [small alpha] < 2. SBM may provide a seemingly adequate description in the case of unbounded diffusion{,} for which its probability density function coincides with that of fractional Brownian motion. Here we show that free SBM is weakly non-ergodic but does not exhibit a significant amplitude scatter of the time averaged mean squared displacement. More severely{,} we demonstrate that under confinement{,} the dynamics encoded by SBM is fundamentally different from both fractional Brownian motion and continuous time random walks. SBM is highly non-stationary and cannot provide a physical description for particles in a thermalised stationary system. Our findings have direct impact on the modelling of single particle tracking experiments{,} in particular{,} under confinement inside cellular compartments or when optical tweezers tracking methods are used.}, language = {en} } @article{BodrovaChechkinCherstvyetal.2015, author = {Bodrova, Anna and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Quantifying non-ergodic dynamics of force-free granular gases}, series = {Physical chemistry, chemical physics : PCCP ; a journal of European Chemical Societies}, journal = {Physical chemistry, chemical physics : PCCP ; a journal of European Chemical Societies}, number = {17}, issn = {1463-9084}, doi = {10.1039/C5CP02824H}, pages = {21791 -- 21798}, year = {2015}, abstract = {Brownianmotion 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 e. 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 e on the impact velocity of particles.}, language = {en} } @misc{BodrovaChechkinCherstvyetal.2015, author = {Bodrova, Anna and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Quantifying non-ergodic dynamics of force-free granular gases}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-85200}, year = {2015}, abstract = {Brownianmotion 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 e. 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 e on the impact velocity of particles.}, language = {en} } @article{MetzlerCherstvyChechkinetal.2015, author = {Metzler, Ralf and Cherstvy, Andrey G. and Chechkin, Aleksei V. and Bodrova, Anna S.}, title = {Ultraslow scaled Brownian motion}, series = {New journal of physics : the open-access journal for physics}, volume = {17}, journal = {New journal of physics : the open-access journal for physics}, number = {063038}, publisher = {Dt. Physikalische Ges., IOP}, address = {Bad Honnef, London}, issn = {1367-2630}, doi = {10.1088/1367-2630/17/6/063038}, year = {2015}, abstract = {We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form . 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.}, language = {en} } @misc{MetzlerCherstvyChechkinetal.2015, author = {Metzler, Ralf and Cherstvy, Andrey G. and Chechkin, Aleksei V. and Bodrova, Anna S.}, title = {Ultraslow scaled Brownian motion}, series = {New journal of physics : the open-access journal for physics}, journal = {New journal of physics : the open-access journal for physics}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-78618}, year = {2015}, abstract = {We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form . 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.}, language = {en} } @article{CherstvyVinodAghionetal.2017, author = {Cherstvy, Andrey G. and Vinod, Deepak and Aghion, Erez and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Time averaging, ageing and delay analysis of financial time series}, series = {New journal of physics}, volume = {19}, journal = {New journal of physics}, publisher = {IOP}, address = {London}, issn = {1367-2630}, doi = {10.1088/1367-2630/aa7199}, pages = {1 -- 11}, year = {2017}, abstract = {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.}, language = {en} } @misc{CherstvyVinodAghionetal.2017, author = {Cherstvy, Andrey G. and Vinod, Deepak and Aghion, Erez and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Time averaging, ageing and delay analysis of financial time series}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-400541}, pages = {11}, year = {2017}, abstract = {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.}, language = {en} } @article{BodrovaChechkinCherstvyetal.2015, author = {Bodrova, Anna and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Quantifying non-ergodic dynamics of force-free granular gases}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {17}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {34}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c5cp02824h}, pages = {21791 -- 21798}, year = {2015}, abstract = {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.}, language = {en} } @article{ChechkinKantzMetzler2017, author = {Chechkin, Aleksei V. and Kantz, Holger and Metzler, Ralf}, title = {Ageing effects in ultraslow continuous time random walks}, series = {The European physical journal : B, Condensed matter and complex systems}, volume = {90}, journal = {The European physical journal : B, Condensed matter and complex systems}, publisher = {Springer}, address = {New York}, issn = {1434-6028}, doi = {10.1140/epjb/e2017-80270-9}, pages = {12}, year = {2017}, abstract = {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.}, language = {en} } @article{CherstvyVinodAghionetal.2017, author = {Cherstvy, Andrey G. and Vinod, Deepak and Aghion, Erez and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Time averaging, ageing and delay analysis of financial time series}, series = {New journal of physics : the open-access journal for physics}, volume = {19}, journal = {New journal of physics : the open-access journal for physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/aa7199}, pages = {135 -- 147}, year = {2017}, abstract = {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.}, language = {en} } @misc{ThapaWyłomańskaSikoraetal.2021, author = {Thapa, Samudrajit and Wyłomańska, Agnieszka and Sikora, Grzegorz and Wagner, Caroline E. and Krapf, Diego and Kantz, Holger and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1118}, issn = {1866-8372}, doi = {10.25932/publishup-49349}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-493494}, pages = {24}, year = {2021}, abstract = {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.}, language = {en} } @article{GuggenbergerChechkinMetzler2021, author = {Guggenberger, Tobias and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Fractional Brownian motion in superharmonic potentials and non-Boltzmann stationary distributions}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {54}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {29}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/ac019b}, pages = {17}, year = {2021}, abstract = {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.}, language = {en} } @article{SandevIominKantzetal.2016, author = {Sandev, T. and Iomin, Alexander and Kantz, Holger and Metzler, Ralf and Chechkin, Aleksei V.}, title = {Comb Model with Slow and Ultraslow Diffusion}, series = {Mathematical modelling of natural phenomena}, volume = {11}, journal = {Mathematical modelling of natural phenomena}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0973-5348}, doi = {10.1051/mmnp/201611302}, pages = {18 -- 33}, year = {2016}, abstract = {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.}, language = {en} } @article{CapałaPadashChechkinetal.2020, author = {Capała, Karol and Padash, Amin and Chechkin, Aleksei V. and Shokri, Babak and Metzler, Ralf and Dybiec, Bartłomiej}, title = {Levy noise-driven escape from arctangent potential wells}, series = {Chaos : an interdisciplinary journal of nonlinear science}, volume = {30}, journal = {Chaos : an interdisciplinary journal of nonlinear science}, number = {12}, publisher = {American Institute of Physics}, address = {Woodbury, NY}, issn = {1054-1500}, doi = {10.1063/5.0021795}, pages = {15}, year = {2020}, abstract = {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.}, language = {en} } @article{ChechkinZaidLomholtetal.2013, author = {Chechkin, Aleksei V. and Zaid, I. M. and Lomholt, M. A. and Sokolov, Igor M. and Metzler, Ralf}, title = {Bulk-mediated surface diffusion on a cylinder in the fast exchange limit}, series = {Mathematical modelling of natural phenomena}, volume = {8}, journal = {Mathematical modelling of natural phenomena}, number = {2}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0973-5348}, doi = {10.1051/mmnp/20138208}, pages = {114 -- 126}, year = {2013}, abstract = {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.}, language = {en} } @article{SandevChechkinKorabeletal.2015, author = {Sandev, Trifce and Chechkin, Aleksei V. and Korabel, Nickolay and Kantz, Holger and Sokolov, Igor M. and Metzler, Ralf}, title = {Distributed-order diffusion equations and multifractality: Models and solutions}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {92}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {4}, publisher = {American Physical Society}, address = {College Park}, issn = {1539-3755}, doi = {10.1103/PhysRevE.92.042117}, pages = {19}, year = {2015}, abstract = {We study distributed-order time fractional diffusion equations characterized by multifractal memory kernels, in contrast to the simple power-law kernel of common time fractional diffusion equations. Based on the physical approach to anomalous diffusion provided by the seminal Scher-Montroll-Weiss continuous time random walk, we analyze both natural and modified-form distributed-order time fractional diffusion equations and compare the two approaches. The mean squared displacement is obtained and its limiting behavior analyzed. We derive the connection between the Wiener process, described by the conventional Langevin equation and the dynamics encoded by the distributed-order time fractional diffusion equation in terms of a generalized subordination of time. A detailed analysis of the multifractal properties of distributed-order diffusion equations is provided.}, language = {en} } @article{ChechkinZaidLomholtetal.2012, author = {Chechkin, Aleksei V. and Zaid, Irwin M. and Lomholt, Michael A. and Sokolov, Igor M. and Metzler, Ralf}, title = {Bulk-mediated diffusion on a planar surface full solution}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {86}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {4}, publisher = {American Physical Society}, address = {College Park}, issn = {1539-3755}, doi = {10.1103/PhysRevE.86.041101}, pages = {11}, year = {2012}, abstract = {We consider the effective surface motion of a particle that intermittently unbinds from a planar surface and performs bulk excursions. Based on a random-walk approach, we derive the diffusion equations for surface and bulk diffusion including the surface-bulk coupling. From these exact dynamic equations, we analytically obtain the propagator of the effective surface motion. This approach allows us to deduce a superdiffusive, Cauchy-type behavior on the surface, together with exact cutoffs limiting the Cauchy form. Moreover, we study the long-time dynamics for the surface motion.}, language = {en} } @article{SandevSokolovMetzleretal.2017, author = {Sandev, Trifce and Sokolov, Igor M. and Metzler, Ralf and Chechkin, Aleksei V.}, title = {Beyond monofractional kinetics}, series = {Chaos, solitons \& fractals : applications in science and engineering ; an interdisciplinary journal of nonlinear science}, volume = {102}, journal = {Chaos, solitons \& fractals : applications in science and engineering ; an interdisciplinary journal of nonlinear science}, publisher = {Elsevier}, address = {Oxford}, issn = {0960-0779}, doi = {10.1016/j.chaos.2017.05.001}, pages = {210 -- 217}, year = {2017}, abstract = {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.}, language = {en} } @article{WangSenoSokolovetal.2020, author = {Wang, Wei and Seno, Flavio and Sokolov, Igor M. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Unexpected crossovers in correlated random-diffusivity processes}, series = {New Journal of Physics}, volume = {22}, journal = {New Journal of Physics}, publisher = {Dt. Physikalische Ges.}, address = {Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/aba390}, pages = {17}, year = {2020}, abstract = {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.}, language = {en} } @article{ChechkinSenoMetzleretal.2017, author = {Chechkin, Aleksei V. and Seno, Flavio and Metzler, Ralf and Sokolov, Igor M.}, title = {Brownian yet Non-Gaussian Diffusion: From Superstatistics to Subordination of Diffusing Diffusivities}, series = {Physical review : X, Expanding access}, volume = {7}, journal = {Physical review : X, Expanding access}, publisher = {American Physical Society}, address = {College Park}, issn = {2160-3308}, doi = {10.1103/PhysRevX.7.021002}, pages = {20}, year = {2017}, abstract = {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.}, language = {en} } @misc{WangSenoSokolovetal.2020, author = {Wang, Wei and Seno, Flavio and Sokolov, Igor M. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Unexpected crossovers in correlated random-diffusivity processes}, number = {1006}, issn = {1866-8372}, doi = {10.25932/publishup-48004}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-480049}, pages = {18}, year = {2020}, abstract = {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.}, language = {en} } @misc{ChechkinZaidLomholtetal.2013, author = {Chechkin, Aleksei V. and Zaid, Irwin M. and Lomholt, Michael A. and Sokolov, Igor M. and Metzler, Ralf}, title = {Bulk-mediated surface diffusion on a cylinder in the fast exchange limit}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, number = {593}, issn = {1866-8372}, doi = {10.25932/publishup-41548}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-415480}, pages = {114 -- 126}, year = {2013}, abstract = {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.}, language = {en} } @article{DoerriesChechkinMetzler2022, author = {Doerries, Timo J. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Apparent anomalous diffusion and non-Gaussian distributions in a simple mobile-immobile transport model with Poissonian switching}, series = {Interface : journal of the Royal Society}, volume = {19}, journal = {Interface : journal of the Royal Society}, number = {192}, publisher = {Royal Society}, address = {London}, issn = {1742-5689}, doi = {10.1098/rsif.2022.0233}, pages = {14}, year = {2022}, abstract = {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.}, language = {en} } @misc{BodrovaChechkinCherstvyetal.2016, author = {Bodrova, Anna S. and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Safdari, Hadiseh and Sokolov, Igor M. and Metzler, Ralf}, title = {Underdamped scaled Brownian motion}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-97158}, pages = {16}, year = {2016}, abstract = {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.}, language = {en} } @article{BodrovaChechkinCherstvyetal.2016, author = {Bodrova, Anna S. and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Safdari, Hadiseh and Sokolov, Igor M. and Metzler, Ralf}, title = {Underdamped scaled Brownian motion}, series = {Scientific reports}, volume = {6}, journal = {Scientific reports}, publisher = {Nature Publishing Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/srep30520}, year = {2016}, abstract = {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.}, language = {en} } @article{SafdariCherstvyChechkinetal.2015, author = {Safdari, Hadiseh and Cherstvy, Andrey G. and Chechkin, Aleksei V. and Thiel, Felix and Sokolov, Igor M. and Metzler, Ralf}, title = {Quantifying the non-ergodicity of scaled Brownian motion}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {48}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {37}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/48/37/375002}, pages = {18}, year = {2015}, abstract = {We examine the non-ergodic properties of scaled Brownian motion (SBM), a non-stationary stochastic process with a time dependent diffusivity of the form D(t) similar or equal to t(alpha-1). We compute the ergodicity breaking parameter EB in the entire range of scaling exponents a, both analytically and via extensive computer simulations of the stochastic Langevin equation. We demonstrate that in the limit of long trajectory lengths T and short lag times Delta the EB parameter as function of the scaling exponent a has no divergence at alpha - 1/2 and present the asymptotes for EB in different limits. We generalize the analytical and simulations results for the time averaged and ergodic properties of SBM in the presence of ageing, that is, when the observation of the system starts only a finite time span after its initiation. The approach developed here for the calculation of the higher time averaged moments of the particle displacement can be applied to derive the ergodic properties of other stochastic processes such as fractional Brownian motion.}, language = {en} } @article{SafdariChechkinJafarietal.2015, author = {Safdari, Hadiseh and Chechkin, Aleksei V. and Jafari, Gholamreza R. and Metzler, Ralf}, title = {Aging scaled Brownian motion}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {91}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {4}, publisher = {American Physical Society}, address = {College Park}, issn = {1539-3755}, doi = {10.1103/PhysRevE.91.042107}, pages = {9}, year = {2015}, abstract = {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.}, language = {en} } @misc{MolinaGarciaSandevSafdarietal.2019, author = {Molina-Garcia, Daniel and Sandev, Trifce and Safdari, Hadiseh and Pagnini, Gianni and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Crossover from anomalous to normal diffusion}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {507}, issn = {1866-8372}, doi = {10.25932/publishup-42259}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-422590}, pages = {28}, year = {2019}, abstract = {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.}, language = {en} } @article{MolinaGarciaSandevSafdarietal.2018, author = {Molina-Garcia, Daniel and Sandev, Trifce and Safdari, Hadiseh and Pagnini, Gianni and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Crossover from anomalous to normal diffusion}, series = {New Journal of Physics}, volume = {20}, journal = {New Journal of Physics}, publisher = {IOP Publishing Ltd}, address = {London und Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/aae4b2}, pages = {28}, year = {2018}, abstract = {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.}, language = {en} } @article{BodrovaChechkinCherstvyetal.2016, author = {Bodrova, Anna S. and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Safdari, Hadiseh and Sokolov, Igor M. and Metzler, Ralf}, title = {Underdamped scaled Brownian motion: (non-)existence of the overdamped limit in anomalous diffusion}, series = {Scientific reports}, volume = {6}, journal = {Scientific reports}, publisher = {Nature Publ. Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/srep30520}, pages = {16}, year = {2016}, abstract = {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.}, language = {en} } @article{SafdariCherstvyChechkinetal.2017, author = {Safdari, Hadiseh and Cherstvy, Andrey G. and Chechkin, Aleksei V. and Bodrova, Anna and Metzler, Ralf}, title = {Aging underdamped scaled Brownian motion}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {95}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.95.012120}, pages = {15}, year = {2017}, abstract = {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.}, language = {en} }