@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} } @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} } @article{GranadoAbadMetzleretal.2020, author = {Granado, Felipe Le Vot and Abad, Enrique and Metzler, Ralf and Yuste, Santos B.}, title = {Continuous time random walk in a velocity field}, 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/ab9ae2}, pages = {27}, year = {2020}, abstract = {We consider the emerging dynamics of a separable continuous time random walk (CTRW) in the case when the random walker is biased by a velocity field in a uniformly growing domain. Concrete examples for such domains include growing biological cells or lipid vesicles, biofilms and tissues, but also macroscopic systems such as expanding aquifers during rainy periods, or the expanding Universe. The CTRW in this study can be subdiffusive, normal diffusive or superdiffusive, including the particular case of a L{\´e}vy flight. We first consider the case when the velocity field is absent. In the subdiffusive case, we reveal an interesting time dependence of the kurtosis of the particle probability density function. In particular, for a suitable parameter choice, we find that the propagator, which is fat tailed at short times, may cross over to a Gaussian-like propagator. We subsequently incorporate the effect of the velocity field and derive a bi-fractional diffusion-advection equation encoding the time evolution of the particle distribution. We apply this equation to study the mixing kinetics of two diffusing pulses, whose peaks move towards each other under the action of velocity fields acting in opposite directions. This deterministic motion of the peaks, together with the diffusive spreading of each pulse, tends to increase particle mixing, thereby counteracting the peak separation induced by the domain growth. As a result of this competition, different regimes of mixing arise. In the case of L{\´e}vy flights, apart from the non-mixing regime, one has two different mixing regimes in the long-time limit, depending on the exact parameter choice: in one of these regimes, mixing is mainly driven by diffusive spreading, while in the other mixing is controlled by the velocity fields acting on each pulse. Possible implications for encounter-controlled reactions in real systems are discussed.}, language = {en} } @misc{GranadoAbadMetzleretal.2020, author = {Granado, Felipe Le Vot and Abad, Enrique and Metzler, Ralf and Yuste, Santos B.}, title = {Continuous time random walk in a velocity field}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1005}, issn = {1866-8372}, doi = {10.25932/publishup-47999}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-479997}, pages = {28}, year = {2020}, abstract = {We consider the emerging dynamics of a separable continuous time random walk (CTRW) in the case when the random walker is biased by a velocity field in a uniformly growing domain. Concrete examples for such domains include growing biological cells or lipid vesicles, biofilms and tissues, but also macroscopic systems such as expanding aquifers during rainy periods, or the expanding Universe. The CTRW in this study can be subdiffusive, normal diffusive or superdiffusive, including the particular case of a L{\´e}vy flight. We first consider the case when the velocity field is absent. In the subdiffusive case, we reveal an interesting time dependence of the kurtosis of the particle probability density function. In particular, for a suitable parameter choice, we find that the propagator, which is fat tailed at short times, may cross over to a Gaussian-like propagator. We subsequently incorporate the effect of the velocity field and derive a bi-fractional diffusion-advection equation encoding the time evolution of the particle distribution. We apply this equation to study the mixing kinetics of two diffusing pulses, whose peaks move towards each other under the action of velocity fields acting in opposite directions. This deterministic motion of the peaks, together with the diffusive spreading of each pulse, tends to increase particle mixing, thereby counteracting the peak separation induced by the domain growth. As a result of this competition, different regimes of mixing arise. In the case of L{\´e}vy flights, apart from the non-mixing regime, one has two different mixing regimes in the long-time limit, depending on the exact parameter choice: in one of these regimes, mixing is mainly driven by diffusive spreading, while in the other mixing is controlled by the velocity fields acting on each pulse. Possible implications for encounter-controlled reactions in real systems are discussed.}, language = {en} } @article{VojtaSkinnerMetzler2019, author = {Vojta, Thomas and Skinner, Sarah and Metzler, Ralf}, title = {Probability density of the fractional Langevin equation with reflecting walls}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {100}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {4}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.100.042142}, pages = {11}, year = {2019}, abstract = {We investigate anomalous diffusion processes governed by the fractional Langevin equation and confined to a finite or semi-infinite interval by reflecting potential barriers. As the random and damping forces in the fractional Langevin equation fulfill the appropriate fluctuation-dissipation relation, the probability density on a finite interval converges for long times towards the expected uniform distribution prescribed by thermal equilibrium. In contrast, on a semi-infinite interval with a reflecting wall at the origin, the probability density shows pronounced deviations from the Gaussian behavior observed for normal diffusion. If the correlations of the random force are persistent (positive), particles accumulate at the reflecting wall while antipersistent (negative) correlations lead to a depletion of particles near the wall. We compare and contrast these results with the strong accumulation and depletion effects recently observed for nonthermal fractional Brownian motion with reflecting walls, and we discuss broader implications.}, language = {en} } @article{Metzler2019, author = {Metzler, Ralf}, title = {Brownian motion and beyond: first-passage, power spectrum, non-Gaussianity, and anomalous diffusion}, series = {Journal of statistical mechanics: theory and experiment}, volume = {2019}, journal = {Journal of statistical mechanics: theory and experiment}, number = {11}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1742-5468}, doi = {10.1088/1742-5468/ab4988}, pages = {18}, year = {2019}, abstract = {Brownian motion is a ubiquitous physical phenomenon across the sciences. After its discovery by Brown and intensive study since the first half of the 20th century, many different aspects of Brownian motion and stochastic processes in general have been addressed in Statistical Physics. In particular, there now exists a very large range of applications of stochastic processes in various disciplines. Here we provide a summary of some of the recent developments in the field of stochastic processes, highlighting both the experimental findings and theoretical frameworks.}, 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{SposiniGrebenkovMetzleretal.2020, author = {Sposini, Vittoria and Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb and Seno, Flavio}, title = {Universal spectral features of different classes of random-diffusivity processes}, series = {New Journal of Physics}, volume = {22}, journal = {New Journal of Physics}, number = {6}, publisher = {Dt. Physikalische Ges.}, address = {Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab9200}, pages = {26}, year = {2020}, abstract = {Stochastic models based on random diffusivities, such as the diffusing-diffusivity approach, are popular concepts for the description of non-Gaussian diffusion in heterogeneous media. Studies of these models typically focus on the moments and the displacement probability density function. Here we develop the complementary power spectral description for a broad class of random-diffusivity processes. In our approach we cater for typical single particle tracking data in which a small number of trajectories with finite duration are garnered. Apart from the diffusing-diffusivity model we study a range of previously unconsidered random-diffusivity processes, for which we obtain exact forms of the probability density function. These new processes are different versions of jump processes as well as functionals of Brownian motion. The resulting behaviour subtly depends on the specific model details. Thus, the central part of the probability density function may be Gaussian or non-Gaussian, and the tails may assume Gaussian, exponential, log-normal, or even power-law forms. For all these models we derive analytically the moment-generating function for the single-trajectory power spectral density. We establish the generic 1/f²-scaling of the power spectral density as function of frequency in all cases. Moreover, we establish the probability density for the amplitudes of the random power spectral density of individual trajectories. The latter functions reflect the very specific properties of the different random-diffusivity models considered here. Our exact results are in excellent agreement with extensive numerical simulations.}, language = {en} } @misc{SposiniGrebenkovMetzleretal.2020, author = {Sposini, Vittoria and Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb and Seno, Flavio}, title = {Universal spectral features of different classes of random-diffusivity processes}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {999}, issn = {1866-8372}, doi = {10.25932/publishup-47696}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-476960}, pages = {27}, year = {2020}, abstract = {Stochastic models based on random diffusivities, such as the diffusing-diffusivity approach, are popular concepts for the description of non-Gaussian diffusion in heterogeneous media. Studies of these models typically focus on the moments and the displacement probability density function. Here we develop the complementary power spectral description for a broad class of random-diffusivity processes. In our approach we cater for typical single particle tracking data in which a small number of trajectories with finite duration are garnered. Apart from the diffusing-diffusivity model we study a range of previously unconsidered random-diffusivity processes, for which we obtain exact forms of the probability density function. These new processes are different versions of jump processes as well as functionals of Brownian motion. The resulting behaviour subtly depends on the specific model details. Thus, the central part of the probability density function may be Gaussian or non-Gaussian, and the tails may assume Gaussian, exponential, log-normal, or even power-law forms. For all these models we derive analytically the moment-generating function for the single-trajectory power spectral density. We establish the generic 1/f²-scaling of the power spectral density as function of frequency in all cases. Moreover, we establish the probability density for the amplitudes of the random power spectral density of individual trajectories. The latter functions reflect the very specific properties of the different random-diffusivity models considered here. Our exact results are in excellent agreement with extensive numerical simulations.}, 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} } @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{LiMeiXuetal.2020, author = {Li, Yongge and Mei, Ruoxing and Xu, Yong and Kurths, J{\"u}rgen and Duan, Jinqiao and Metzler, Ralf}, title = {Particle dynamics and transport enhancement in a confined channel with position-dependent diffusivity}, 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/ab81b9}, pages = {27}, year = {2020}, abstract = {This work focuses on the dynamics of particles in a confined geometry with position-dependent diffusivity, where the confinement is modelled by a periodic channel consisting of unit cells connected by narrow passage ways. We consider three functional forms for the diffusivity, corresponding to the scenarios of a constant (D ₀), as well as a low (D ₘ) and a high (D d) mobility diffusion in cell centre of the longitudinally symmetric cells. Due to the interaction among the diffusivity, channel shape and external force, the system exhibits complex and interesting phenomena. By calculating the probability density function, mean velocity and mean first exit time with the It{\^o} calculus form, we find that in the absence of external forces the diffusivity D d will redistribute particles near the channel wall, while the diffusivity D ₘ will trap them near the cell centre. The superposition of external forces will break their static distributions. Besides, our results demonstrate that for the diffusivity D d, a high dependence on the x coordinate (parallel with the central channel line) will improve the mean velocity of the particles. In contrast, for the diffusivity D ₘ, a weak dependence on the x coordinate will dramatically accelerate the moving speed. In addition, it shows that a large external force can weaken the influences of different diffusivities; inversely, for a small external force, the types of diffusivity affect significantly the particle dynamics. In practice, one can apply these results to achieve a prominent enhancement of the particle transport in two- or three-dimensional channels by modulating the local tracer diffusivity via an engineered gel of varying porosity or by adding a cold tube to cool down the diffusivity along the central line, which may be a relevant effect in engineering applications. Effects of different stochastic calculi in the evaluation of the underlying multiplicative stochastic equation for different physical scenarios are discussed.}, language = {en} } @misc{LiMeiXuetal.2020, author = {Li, Yongge and Mei, Ruoxing and Xu, Yong and Kurths, J{\"u}rgen and Duan, Jinqiao and Metzler, Ralf}, title = {Particle dynamics and transport enhancement in a confined channel with position-dependent diffusivity}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {974}, issn = {1866-8372}, doi = {10.25932/publishup-47454}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-474542}, pages = {28}, year = {2020}, abstract = {This work focuses on the dynamics of particles in a confined geometry with position-dependent diffusivity, where the confinement is modelled by a periodic channel consisting of unit cells connected by narrow passage ways. We consider three functional forms for the diffusivity, corresponding to the scenarios of a constant (D ₀), as well as a low (D ₘ) and a high (D d) mobility diffusion in cell centre of the longitudinally symmetric cells. Due to the interaction among the diffusivity, channel shape and external force, the system exhibits complex and interesting phenomena. By calculating the probability density function, mean velocity and mean first exit time with the It{\^o} calculus form, we find that in the absence of external forces the diffusivity D d will redistribute particles near the channel wall, while the diffusivity D ₘ will trap them near the cell centre. The superposition of external forces will break their static distributions. Besides, our results demonstrate that for the diffusivity D d, a high dependence on the x coordinate (parallel with the central channel line) will improve the mean velocity of the particles. In contrast, for the diffusivity D ₘ, a weak dependence on the x coordinate will dramatically accelerate the moving speed. In addition, it shows that a large external force can weaken the influences of different diffusivities; inversely, for a small external force, the types of diffusivity affect significantly the particle dynamics. In practice, one can apply these results to achieve a prominent enhancement of the particle transport in two- or three-dimensional channels by modulating the local tracer diffusivity via an engineered gel of varying porosity or by adding a cold tube to cool down the diffusivity along the central line, which may be a relevant effect in engineering applications. Effects of different stochastic calculi in the evaluation of the underlying multiplicative stochastic equation for different physical scenarios are discussed.}, language = {en} } @article{ThapaLomholtKrogetal.2018, author = {Thapa, Samudrajit and Lomholt, Michael Andersen and Krog, Jens and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Bayesian analysis of single-particle tracking data using the nested-sampling algorithm: maximum-likelihood model selection applied to stochastic-diffusivity data}, series = {Physical chemistry, chemical physics : PCCP ; a journal of European Chemical Societies}, volume = {20}, journal = {Physical chemistry, chemical physics : PCCP ; a journal of European Chemical Societies}, number = {46}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c8cp04043e}, pages = {29018 -- 29037}, year = {2018}, abstract = {We employ Bayesian statistics using the nested-sampling algorithm to compare and rank multiple models of ergodic diffusion (including anomalous diffusion) as well as to assess their optimal parameters for in silico-generated and real time-series. We focus on the recently-introduced model of Brownian motion with "diffusing diffusivity'-giving rise to widely-observed non-Gaussian displacement statistics-and its comparison to Brownian and fractional Brownian motion, also for the time-series with some measurement noise. We conduct this model-assessment analysis using Bayesian statistics and the nested-sampling algorithm on the level of individual particle trajectories. We evaluate relative model probabilities and compute best-parameter sets for each diffusion model, comparing the estimated parameters to the true ones. We test the performance of the nested-sampling algorithm and its predictive power both for computer-generated (idealised) trajectories as well as for real single-particle-tracking trajectories. Our approach delivers new important insight into the objective selection of the most suitable stochastic model for a given time-series. We also present first model-ranking results in application to experimental data of tracer diffusion in polymer-based hydrogels.}, 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} } @article{SchwarzlGodecMetzler2017, author = {Schwarzl, Maria and Godec, Aljaz and Metzler, Ralf}, title = {Quantifying non-ergodicity of anomalous diffusion with higher order moments}, series = {Scientific reports}, volume = {7}, journal = {Scientific reports}, publisher = {Nature Publ. Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-017-03712-x}, pages = {18}, year = {2017}, abstract = {Anomalous diffusion is being discovered in a fast growing number of systems. The exact nature of this anomalous diffusion provides important information on the physical laws governing the studied system. One of the central properties analysed for finite particle motion time series is the intrinsic variability of the apparent diffusivity, typically quantified by the ergodicity breaking parameter EB. Here we demonstrate that frequently EB is insufficient to provide a meaningful measure for the observed variability of the data. Instead, important additional information is provided by the higher order moments entering by the skewness and kurtosis. We analyse these quantities for three popular anomalous diffusion models. In particular, we find that even for the Gaussian fractional Brownian motion a significant skewness in the results of physical measurements occurs and needs to be taken into account. Interestingly, the kurtosis and skewness may also provide sensitive estimates of the anomalous diffusion exponent underlying the data. We also derive a new result for the EB parameter of fractional Brownian motion valid for the whole range of the anomalous diffusion parameter. Our results are important for the analysis of anomalous diffusion but also provide new insights into the theory of anomalous stochastic processes.}, 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} } @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{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{JavanainenMartinezSearaMetzleretal.2017, author = {Javanainen, Matti and Martinez-Seara, Hector and Metzler, Ralf and Vattulainen, Ilpo}, title = {Diffusion of Integral Membrane Proteins in Protein-Rich Membranes}, series = {The journal of physical chemistry letters}, volume = {8}, journal = {The journal of physical chemistry letters}, publisher = {American Chemical Society}, address = {Washington}, issn = {1948-7185}, doi = {10.1021/acs.jpclett.7b01758}, pages = {4308 -- 4313}, year = {2017}, abstract = {The lateral diffusion of embedded proteins along lipid membranes in protein-poor conditions has been successfully described in terms of the Saffman-Delbruck (SD) model, which predicts that the protein diffusion coefficient D is weakly dependent on its radius R as D proportional to ln(1/R). However, instead of being protein-poor, native cell membranes are extremely crowded with proteins. On the basis of extensive molecular simulations, we here demonstrate that protein crowding of the membrane at physiological levels leads to deviations from the SD relation and to the emergence of a stronger Stokes-like dependence D proportional to 1/R. We propose that this 1/R law mainly arises due to geometrical factors: smaller proteins are able to avoid confinement effects much better than their larger counterparts. The results highlight that the lateral dynamics in the crowded setting found in native membranes is radically different from protein-poor conditions and plays a significant role in formation of functional multiprotein complexes.}, language = {en} } @article{CaetanodeCarvalhoMetzleretal.2017, author = {Caetano, Daniel L. Z. and de Carvalho, Sidney J. and Metzler, Ralf and Cherstvy, Andrey G.}, title = {Critical adsorption of periodic and random polyampholytes onto charged surfaces}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {19}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c7cp04040g}, pages = {23397 -- 23413}, year = {2017}, abstract = {How different are the properties of critical adsorption of polyampholytes and polyelectrolytes onto charged surfaces? How important are the details of polyampholyte charge distribution on the onset of critical adsorption transition? What are the scaling relations governing the dependence of critical surface charge density on salt concentration in the surrounding solution? Here, we employ Metropolis Monte Carlo simulations and uncover the scaling relations for critical adsorption for quenched periodic and random charge distributions along the polyampholyte chains. We also evaluate and discuss the dependence of the adsorbed layer width on solution salinity and details of the charge distribution. We contrast our findings to the known results for polyelectrolyte adsorption onto oppositely charged surfaces, in particular, their dependence on electrolyte concentration.}, language = {en} } @article{HerrmannMetzlerEngbert2017, author = {Herrmann, Carl J. J. and Metzler, Ralf and Engbert, Ralf}, title = {A self-avoiding walk with neural delays as a model of fixational eye movements}, series = {Scientific reports}, volume = {7}, journal = {Scientific reports}, publisher = {Nature Publ. Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-017-13489-8}, pages = {17}, year = {2017}, abstract = {Fixational eye movements show scaling behaviour of the positional mean-squared displacement with a characteristic transition from persistence to antipersistence for increasing time-lag. These statistical patterns were found to be mainly shaped by microsaccades (fast, small-amplitude movements). However, our re-analysis of fixational eye-movement data provides evidence that the slow component (physiological drift) of the eyes exhibits scaling behaviour of the mean-squared displacement that varies across human participants. These results suggest that drift is a correlated movement that interacts with microsaccades. Moreover, on the long time scale, the mean-squared displacement of the drift shows oscillations, which is also present in the displacement auto-correlation function. This finding lends support to the presence of time-delayed feedback in the control of drift movements. Based on an earlier non-linear delayed feedback model of fixational eye movements, we propose and discuss different versions of a new model that combines a self-avoiding walk with time delay. As a result, we identify a model that reproduces oscillatory correlation functions, the transition from persistence to antipersistence, and microsaccades.}, language = {en} } @article{GrebenkovMetzlerOshanin2017, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains}, 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/aa8ed9}, pages = {11}, year = {2017}, abstract = {We study the mean first passage time (MFPT) to a reaction event on a specific site in a cylindrical geometry-characteristic, for instance, for bacterial cells, with a concentric inner cylinder representing the nuclear region of the bacterial cell. A similar problem emerges in the description of a diffusive search by a transcription factor protein for a specific binding region on a single strand of DNA. We develop a unified theoretical approach to study the underlying boundary value problem which is based on a self-consistent approximation of the mixed boundary condition. Our approach permits us to derive explicit, novel, closed-form expressions for the MFPT valid for a generic setting with an arbitrary relation between the system parameters. We analyse this general result in the asymptotic limits appropriate for the above-mentioned biophysical problems. Our investigation reveals the crucial role of the target aspect ratio and of the intrinsic reactivity of the binding region, which were disregarded in previous studies. Theoretical predictions are confirmed by numerical simulations.}, 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{SandevIominKantzetal.2016, author = {Sandev, Trifce 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{GhoshCherstvyPetrovetal.2016, author = {Ghosh, Surya K. and Cherstvy, Andrey G. and Petrov, Eugene P. and Metzler, Ralf}, title = {Interactions of rod-like particles on responsive elastic sheets}, series = {Soft matter}, volume = {12}, journal = {Soft matter}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, doi = {10.1039/c6sm01522k}, pages = {7908 -- 7919}, year = {2016}, abstract = {What are the physical laws of the mutual interactions of objects bound to cell membranes, such as various membrane proteins or elongated virus particles? To rationalise this, we here investigate by extensive computer simulations mutual interactions of rod-like particles adsorbed on the surface of responsive elastic two-dimensional sheets. Specifically, we quantify sheet deformations as a response to adhesion of such filamentous particles. We demonstrate that tip-to-tip contacts of rods are favoured for relatively soft sheets, while side-by-side contacts are preferred for stiffer elastic substrates. These attractive orientation-dependent substrate-mediated interactions between the rod-like particles on responsive sheets can drive their aggregation and self-assembly. The optimal orientation of the membrane-bound rods is established via responding to the elastic energy profiles created around the particles. We unveil the phase diagramme of attractive-repulsive rod-rod interactions in the plane of their separation and mutual orientation. Applications of our results to other systems featuring membrane-associated particles are also discussed.}, language = {en} } @article{GhoshCherstvyGrebenkovetal.2016, author = {Ghosh, Surya K. and Cherstvy, Andrey G. and Grebenkov, Denis S. and Metzler, Ralf}, title = {Anomalous, non-Gaussian tracer diffusion in crowded two-dimensional environments}, series = {NEW JOURNAL OF PHYSICS}, volume = {18}, journal = {NEW JOURNAL OF PHYSICS}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/18/1/013027}, pages = {16}, year = {2016}, abstract = {A topic of intense current investigation pursues the question of how the highly crowded environment of biological cells affects the dynamic properties of passively diffusing particles. Motivated by recent experiments we report results of extensive simulations of the motion of a finite sized tracer particle in a heterogeneously crowded environment made up of quenched distributions of monodisperse crowders of varying sizes in finite circular two-dimensional domains. For given spatial distributions of monodisperse crowders we demonstrate how anomalous diffusion with strongly non-Gaussian features arises in this model system. We investigate both biologically relevant situations of particles released either at the surface of an inner domain or at the outer boundary, exhibiting distinctly different features of the observed anomalous diffusion for heterogeneous distributions of crowders. Specifically we reveal an asymmetric spreading of tracers even at moderate crowding. In addition to the mean squared displacement (MSD) and local diffusion exponent we investigate the magnitude and the amplitude scatter of the time averaged MSD of individual tracer trajectories, the non-Gaussianity parameter, and the van Hove correlation function. We also quantify how the average tracer diffusivity varies with the position in the domain with a heterogeneous radial distribution of crowders and examine the behaviour of the survival probability and the dynamics of the tracer survival probability. Inter alia, the systems we investigate are related to the passive transport of lipid molecules and proteins in two-dimensional crowded membranes or the motion in colloidal solutions or emulsions in effectively two-dimensional geometries, as well as inside supercrowded, surface adhered cells.}, language = {en} } @misc{Metzler2016, author = {Metzler, Ralf}, title = {PROTEIN PHYSICS Forever ageing}, series = {Nature physics}, volume = {12}, journal = {Nature physics}, publisher = {Nature Publ. Group}, address = {London}, issn = {1745-2473}, doi = {10.1038/nphys3585}, pages = {113 -- 114}, year = {2016}, abstract = {Single-molecule techniques have long given us insight into the motion and interactions of individual molecules. But simulations now show that the dynamics inside single proteins is not as simple as we thought — and that proteins are forever changing.}, language = {en} } @article{GodecMetzler2016, author = {Godec, Aljaz and Metzler, Ralf}, title = {First passage time distribution in heterogeneity controlled kinetics: going beyond the mean first passage time}, series = {Scientific reports}, volume = {6}, journal = {Scientific reports}, publisher = {Nature Publ. Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/srep20349}, pages = {11}, year = {2016}, abstract = {The first passage is a generic concept for quantifying when a random quantity such as the position of a diffusing molecule or the value of a stock crosses a preset threshold (target) for the first time. The last decade saw an enlightening series of new results focusing mostly on the so-called mean and global first passage time (MFPT and GFPT, respectively) of such processes. Here we push the understanding of first passage processes one step further. For a simple heterogeneous system we derive rigorously the complete distribution of first passage times (FPTs). Our results demonstrate that the typical FPT significantly differs from the MFPT, which corresponds to the long time behaviour of the FPT distribution. Conversely, the short time behaviour is shown to correspond to trajectories connecting directly from the initial value to the target. Remarkably, we reveal a previously overlooked third characteristic time scale of the first passage dynamics mirroring brief excursion away from the target.}, language = {en} } @article{JeonJavanainenMartinezSearaetal.2016, author = {Jeon, Jae-Hyung and Javanainen, Matti and Martinez-Seara, Hector and Metzler, Ralf and Vattulainen, Ilpo}, title = {Protein Crowding in Lipid Bilayers Gives Rise to Non-Gaussian Anomalous Lateral Diffusion of Phospholipids and Proteins}, series = {Physical review : X, Expanding access}, volume = {6}, journal = {Physical review : X, Expanding access}, publisher = {American Physical Society}, address = {College Park}, issn = {2160-3308}, doi = {10.1103/PhysRevX.6.021006}, pages = {17}, year = {2016}, abstract = {Biomembranes are exceptionally crowded with proteins with typical protein-to-lipid ratios being around 1:50 - 1:100. Protein crowding has a decisive role in lateral membrane dynamics as shown by recent experimental and computational studies that have reported anomalous lateral diffusion of phospholipids and membrane proteins in crowded lipid membranes. Based on extensive simulations and stochastic modeling of the simulated trajectories, we here investigate in detail how increasing crowding by membrane proteins reshapes the stochastic characteristics of the anomalous lateral diffusion in lipid membranes. We observe that correlated Gaussian processes of the fractional Langevin equation type, identified as the stochastic mechanism behind lipid motion in noncrowded bilayer, no longer adequately describe the lipid and protein motion in crowded but otherwise identical membranes. It turns out that protein crowding gives rise to a multifractal, non-Gaussian, and spatiotemporally heterogeneous anomalous lateral diffusion on time scales from nanoseconds to, at least, tens of microseconds. Our investigation strongly suggests that the macromolecular complexity and spatiotemporal membrane heterogeneity in cellular membranes play critical roles in determining the stochastic nature of the lateral diffusion and, consequently, the associated dynamic phenomena within membranes. Clarifying the exact stochastic mechanism for various kinds of biological membranes is an important step towards a quantitative understanding of numerous intramembrane dynamic phenomena.}, language = {en} } @article{SchwarzlGodecOshaninetal.2016, author = {Schwarzl, Maria and Godec, Aljaz and Oshanin, Gleb and Metzler, Ralf}, title = {A single predator charging a herd of prey: effects of self volume and predator-prey decision-making}, 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/22/225601}, pages = {19}, year = {2016}, abstract = {We study the degree of success of a single predator hunting a herd of prey on a two-dimensional square lattice landscape. We explicitly consider the self volume of the prey restraining their dynamics on the lattice. The movement of both predator and prey is chosen to include an intelligent, decision making step based on their respective sighting ranges, the radius in which they can detect the other species (prey cannot recognise each other besides the self volume interaction): after spotting each other the motion of prey and predator turns from a nearest neighbour random walk into directed escape or chase, respectively. We consider a large range of prey densities and sighting ranges and compute the mean first passage time for a predator to catch a prey as well as characterise the effective dynamics of the hunted prey. We find that the prey's sighting range dominates their life expectancy and the predator profits more from a bad eyesight of the prey than from his own good eye sight. We characterise the dynamics in terms of the mean distance between the predator and the nearest prey. It turns out that effectively the dynamics of this distance coordinate can be captured in terms of a simple Ornstein-Uhlenbeck picture. Reducing the many-body problem to a simple two-body problem by imagining predator and nearest prey to be connected by an effective Hookean bond, all features of the model such as prey density and sighting ranges merge into the effective binding constant.}, 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{deCarvalhoMetzlerCherstvy2016, author = {de Carvalho, Sidney J. and Metzler, Ralf and Cherstvy, Andrey G.}, title = {Critical adsorption of polyelectrolytes onto planar and convex highly charged surfaces: the nonlinear Poisson-Boltzmann approach}, series = {NEW JOURNAL OF PHYSICS}, volume = {18}, journal = {NEW JOURNAL OF PHYSICS}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/18/8/083037}, pages = {17}, year = {2016}, abstract = {We study the adsorption-desorption transition of polyelectrolyte chains onto planar, cylindrical and spherical surfaces with arbitrarily high surface charge densities by massive Monte Carlo computer simulations. We examine in detail how the well known scaling relations for the threshold transition demarcating the adsorbed and desorbed domains of a polyelectrolyte near weakly charged surfaces-are altered for highly charged interfaces. In virtue of high surface potentials and large surface charge densities, the Debye-Huckel approximation is often not feasible and the nonlinear Poisson-Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson-Boltzmann equation is smaller than the Debye-Huckel result, such that the required critical surface charge density for polyelectrolyte adsorption sigma(c) increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical-chemical and biophysical systems.}, language = {en} } @article{GodecMetzler2016, author = {Godec, Aljaz and Metzler, Ralf}, title = {Active transport improves the precision of linear long distance molecular signalling}, 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/36/364001}, pages = {11}, year = {2016}, abstract = {Molecular signalling in living cells occurs at low copy numbers and is thereby inherently limited by the noise imposed by thermal diffusion. The precision at which biochemical receptors can count signalling molecules is intimately related to the noise correlation time. In addition to passive thermal diffusion, messenger RNA and vesicle-engulfed signalling molecules can transiently bind to molecular motors and are actively transported across biological cells. Active transport is most beneficial when trafficking occurs over large distances, for instance up to the order of 1 metre in neurons. Here we explain how intermittent active transport allows for faster equilibration upon a change in concentration triggered by biochemical stimuli. Moreover, we show how intermittent active excursions induce qualitative changes in the noise in effectively one-dimensional systems such as dendrites. Thereby they allow for significantly improved signalling precision in the sense of a smaller relative deviation in the concentration read-out by the receptor. On the basis of linear response theory we derive the exact mean field precision limit for counting actively transported molecules. We explain how intermittent active excursions disrupt the recurrence in the molecular motion, thereby facilitating improved signalling accuracy. Our results provide a deeper understanding of how recurrence affects molecular signalling precision in biological cells and novel medical-diagnostic devices.}, language = {en} } @article{CherstvyMetzler2016, author = {Cherstvy, Andrey G. and Metzler, Ralf}, title = {Anomalous diffusion in time-fluctuating non-stationary diffusivity landscapes}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {18}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c6cp03101c}, pages = {23840 -- 23852}, year = {2016}, abstract = {We investigate the ensemble and time averaged mean squared displacements for particle diffusion in a simple model for disordered media by assuming that the local diffusivity is both fluctuating in time and has a deterministic average growth or decay in time. In this study we compare computer simulations of the stochastic Langevin equation for this random diffusion process with analytical results. We explore the regimes of normal Brownian motion as well as anomalous diffusion in the sub- and superdiffusive regimes. We also consider effects of the inertial term on the particle motion. The investigation of the resulting diffusion is performed for unconfined and confined motion.}, 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{MetzlerJeonCherstvy2016, author = {Metzler, Ralf and Jeon, J. -H. and Cherstvy, Andrey G.}, title = {Non-Brownian diffusion in lipid membranes: Experiments and simulations}, series = {Biochimica et biophysica acta : Biomembranes}, volume = {1858}, journal = {Biochimica et biophysica acta : Biomembranes}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0005-2736}, doi = {10.1016/j.bbamem.2016.01.022}, pages = {2451 -- 2467}, year = {2016}, abstract = {The dynamics of constituents and the surface response of cellular membranes also in connection to the binding of various particles and macromolecules to the membrane are still a matter of controversy in the membrane biophysics community, particularly with respect to crowded membranes of living biological cells. We here put into perspective recent single particle tracking experiments in the plasma membranes of living cells and supercomputing studies of lipid bilayer model membranes with and without protein crowding. Special emphasis is put on the observation of anomalous, non-Brownian diffusion of both lipid molecules and proteins embedded in the lipid bilayer. While single component, pure lipid bilayers in simulations exhibit only transient anomalous diffusion of lipid molecules on nanosecond time scales, the persistence of anomalous diffusion becomes significantly longer ranged on the addition of disorder through the addition of cholesterol or proteins and on passing of the membrane lipids to the gel phase. Concurrently, experiments demonstrate the anomalous diffusion of membrane embedded proteins up to macroscopic time scales in the minute time range. Particular emphasis will be put on the physical character of the anomalous diffusion, in particular, the occurrence of ageing observed in the experiments the effective diffusivity of the measured particles is a decreasing function of time. Moreover, we present results for the time dependent local scaling exponent of the mean squared displacement of the monitored particles. Recent results finding deviations from the commonly assumed Gaussian diffusion patterns in protein crowded membranes are reported. The properties of the displacement autocorrelation function of the lipid molecules are discussed in the light of their appropriate physical anomalous diffusion models, both for non-crowded and crowded membranes. In the last part of this review we address the upcoming field of membrane distortion by elongated membrane-binding particles. We discuss how membrane compartmentalisation and the particle-membrane binding energy may impact the dynamics and response of lipid membranes. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Rog. (C) 2016 The Authors. Published by Elsevier B.V.}, language = {en} } @article{KruesemannSchwarzlMetzler2016, author = {Kruesemann, Henning and Schwarzl, Richard and Metzler, Ralf}, title = {Ageing Scher-Montroll Transport}, series = {Transport in Porous Media}, volume = {115}, journal = {Transport in Porous Media}, publisher = {Springer}, address = {New York}, issn = {0169-3913}, doi = {10.1007/s11242-016-0686-y}, pages = {327 -- 344}, year = {2016}, abstract = {We study the properties of ageing Scher-Montroll transport in terms of a biased subdiffusive continuous time random walk in which the waiting times between consecutive jumps of the charge carriers are distributed according to the power law probability with . As we show, the dynamical properties of the Scher-Montroll transport depend on the ageing time span between the initial preparation of the system and the start of the observation. The Scher-Montroll transport theory was originally shown to describe the photocurrent in amorphous solids in the presence of an external electric field, but it has since been used in many other fields of physical sciences, in particular also in the geophysical context for the description of the transport of tracer particles in subsurface aquifers. In the absence of ageing () the photocurrent of the classical Scher-Montroll model or the breakthrough curves in the groundwater context exhibit a crossover between two power law regimes in time with the scaling exponents and . In the presence of ageing a new power law regime and an initial plateau regime of the current emerge. We derive the different power law regimes and crossover times of the ageing Scher-Montroll transport and show excellent agreement with simulations of the process. Experimental data of ageing Scher-Montroll transport in polymeric semiconductors are shown to agree well with the predictions of our theory.}, language = {en} } @article{GodecMetzler2016, author = {Godec, Aljaz and Metzler, Ralf}, title = {Universal Proximity Effect in Target Search Kinetics in the Few-Encounter Limit}, series = {Physical review : X, Expanding access}, volume = {6}, journal = {Physical review : X, Expanding access}, publisher = {American Physical Society}, address = {College Park}, issn = {2160-3308}, doi = {10.1103/PhysRevX.6.041037}, pages = {11}, year = {2016}, abstract = {When does a diffusing particle reach its target for the first time? This first-passage time (FPT) problem is central to the kinetics of molecular reactions in chemistry and molecular biology. Here, we explain the behavior of smooth FPT densities, for which all moments are finite, and demonstrate universal yet generally non-Poissonian long-time asymptotics for a broad variety of transport processes. While Poisson-like asymptotics arise generically in the presence of an effective repulsion in the immediate vicinity of the target, a time-scale separation between direct and reflected indirect trajectories gives rise to a universal proximity effect: Direct paths, heading more or less straight from the point of release to the target, become typical and focused, with a narrow spread of the corresponding first-passage times. Conversely, statistically dominant indirect paths exploring the entire system tend to be massively dissimilar. The initial distance to the target particularly impacts gene regulatory or competitive stochastic processes, for which few binding events often determine the regulatory outcome. The proximity effect is independent of details of the transport, highlighting the robust character of the FPT features uncovered here.}, language = {en} } @article{GrebenkovMetzlerOshanin2019, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Full distribution of first exit times in the narrow escape problem}, 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/ab5de4}, pages = {23}, year = {2019}, abstract = {In the scenario of the narrow escape problem (NEP) a particle diffuses in a finite container and eventually leaves it through a small 'escape window' in the otherwise impermeable boundary, once it arrives to this window and crosses an entropic barrier at the entrance to it. This generic problem is mathematically identical to that of a diffusion-mediated reaction with a partially-reactive site on the container's boundary. Considerable knowledge is available on the dependence of the mean first-reaction time (FRT) on the pertinent parameters. We here go a distinct step further and derive the full FRT distribution for the NEP. We demonstrate that typical FRTs may be orders of magnitude shorter than the mean one, thus resulting in a strong defocusing of characteristic temporal scales. We unveil the geometry-control of the typical times, emphasising the role of the initial distance to the target as a decisive parameter. A crucial finding is the further FRT defocusing due to the barrier, necessitating repeated escape or reaction attempts interspersed with bulk excursions. These results add new perspectives and offer a broad comprehension of various features of the by-now classical NEP that are relevant for numerous biological and technological systems.}, language = {en} } @misc{GrebenkovMetzlerOshanin2019, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Full distribution of first exit times in the narrow escape problem}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {810}, issn = {1866-8372}, doi = {10.25932/publishup-44288}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-442883}, pages = {24}, year = {2019}, abstract = {In the scenario of the narrow escape problem (NEP) a particle diffuses in a finite container and eventually leaves it through a small 'escape window' in the otherwise impermeable boundary, once it arrives to this window and crosses an entropic barrier at the entrance to it. This generic problem is mathematically identical to that of a diffusion-mediated reaction with a partially-reactive site on the container's boundary. Considerable knowledge is available on the dependence of the mean first-reaction time (FRT) on the pertinent parameters. We here go a distinct step further and derive the full FRT distribution for the NEP. We demonstrate that typical FRTs may be orders of magnitude shorter than the mean one, thus resulting in a strong defocusing of characteristic temporal scales. We unveil the geometry-control of the typical times, emphasising the role of the initial distance to the target as a decisive parameter. A crucial finding is the further FRT defocusing due to the barrier, necessitating repeated escape or reaction attempts interspersed with bulk excursions. These results add new perspectives and offer a broad comprehension of various features of the by-now classical NEP that are relevant for numerous biological and technological systems.}, 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} } @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} } @misc{KindlerPulkkinenCherstvyetal.2019, author = {Kindler, Oliver and Pulkkinen, Otto and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Burst Statistics in an Early Biofilm Quorum Sensing Mode}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {777}, issn = {1866-8372}, doi = {10.25932/publishup-43909}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-439099}, pages = {21}, year = {2019}, abstract = {Quorum-sensing bacteria in a growing colony of cells send out signalling molecules (so-called "autoinducers") and themselves sense the autoinducer concentration in their vicinity. Once—due to increased local cell density inside a "cluster" of the growing colony—the concentration of autoinducers exceeds a threshold value, cells in this clusters get "induced" into a communal, multi-cell biofilm-forming mode in a cluster-wide burst event. We analyse quantitatively the influence of spatial disorder, the local heterogeneity of the spatial distribution of cells in the colony, and additional physical parameters such as the autoinducer signal range on the induction dynamics of the cell colony. Spatial inhomogeneity with higher local cell concentrations in clusters leads to earlier but more localised induction events, while homogeneous distributions lead to comparatively delayed but more concerted induction of the cell colony, and, thus, a behaviour close to the mean-field dynamics. We quantify the induction dynamics with quantifiers such as the time series of induction events and burst sizes, the grouping into induction families, and the mean autoinducer concentration levels. Consequences for different scenarios of biofilm growth are discussed, providing possible cues for biofilm control in both health care and biotechnology.}, language = {en} } @article{KindlerPulkkinenCherstvyetal.2019, author = {Kindler, Oliver and Pulkkinen, Otto and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Burst Statistics in an Early Biofilm Quorum Sensing Mode}, series = {Scientific Reports}, volume = {9}, journal = {Scientific Reports}, publisher = {Macmillan Publishers Limited part of Springer Nature}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-019-48525-2}, pages = {19}, year = {2019}, abstract = {Quorum-sensing bacteria in a growing colony of cells send out signalling molecules (so-called "autoinducers") and themselves sense the autoinducer concentration in their vicinity. Once—due to increased local cell density inside a "cluster" of the growing colony—the concentration of autoinducers exceeds a threshold value, cells in this clusters get "induced" into a communal, multi-cell biofilm-forming mode in a cluster-wide burst event. We analyse quantitatively the influence of spatial disorder, the local heterogeneity of the spatial distribution of cells in the colony, and additional physical parameters such as the autoinducer signal range on the induction dynamics of the cell colony. Spatial inhomogeneity with higher local cell concentrations in clusters leads to earlier but more localised induction events, while homogeneous distributions lead to comparatively delayed but more concerted induction of the cell colony, and, thus, a behaviour close to the mean-field dynamics. We quantify the induction dynamics with quantifiers such as the time series of induction events and burst sizes, the grouping into induction families, and the mean autoinducer concentration levels. Consequences for different scenarios of biofilm growth are discussed, providing possible cues for biofilm control in both health care and biotechnology.}, language = {en} } @misc{ŚlęzakBurneckiMetzler2019, author = {Ślęzak, Jakub and Burnecki, Krzysztof and Metzler, Ralf}, title = {Random coefficient autoregressive processes describe Brownian yet non-Gaussian diffusion in heterogeneous systems}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {765}, issn = {1866-8372}, doi = {10.25932/publishup-43792}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-437923}, pages = {18}, year = {2019}, abstract = {Many studies on biological and soft matter systems report the joint presence of a linear mean-squared displacement and a non-Gaussian probability density exhibiting, for instance, exponential or stretched-Gaussian tails. This phenomenon is ascribed to the heterogeneity of the medium and is captured by random parameter models such as 'superstatistics' or 'diffusing diffusivity'. Independently, scientists working in the area of time series analysis and statistics have studied a class of discrete-time processes with similar properties, namely, random coefficient autoregressive models. In this work we try to reconcile these two approaches and thus provide a bridge between physical stochastic processes and autoregressive models.Westart from the basic Langevin equation of motion with time-varying damping or diffusion coefficients and establish the link to random coefficient autoregressive processes. By exploring that link we gain access to efficient statistical methods which can help to identify data exhibiting Brownian yet non-Gaussian diffusion.}, language = {en} } @article{ŚlęzakBurneckiMetzler2019, author = {Ślęzak, Jakub and Burnecki, Krzysztof and Metzler, Ralf}, title = {Random coefficient autoregressive processes describe Brownian yet non-Gaussian diffusion in heterogeneous systems}, series = {New Journal of Physics}, volume = {21}, journal = {New Journal of Physics}, publisher = {Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics}, address = {Bad Honnef und London}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab3366}, pages = {18}, year = {2019}, abstract = {Many studies on biological and soft matter systems report the joint presence of a linear mean-squared displacement and a non-Gaussian probability density exhibiting, for instance, exponential or stretched-Gaussian tails. This phenomenon is ascribed to the heterogeneity of the medium and is captured by random parameter models such as 'superstatistics' or 'diffusing diffusivity'. Independently, scientists working in the area of time series analysis and statistics have studied a class of discrete-time processes with similar properties, namely, random coefficient autoregressive models. In this work we try to reconcile these two approaches and thus provide a bridge between physical stochastic processes and autoregressive models.Westart from the basic Langevin equation of motion with time-varying damping or diffusion coefficients and establish the link to random coefficient autoregressive processes. By exploring that link we gain access to efficient statistical methods which can help to identify data exhibiting Brownian yet non-Gaussian diffusion.}, language = {en} } @misc{GuggenbergerPagniniVojtaetal.2019, author = {Guggenberger, Tobias and Pagnini, Gianni and Vojta, Thomas and Metzler, Ralf}, title = {Fractional Brownian motion in a finite interval}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {755}, issn = {1866-8372}, doi = {10.25932/publishup-43666}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-436665}, pages = {13}, year = {2019}, abstract = {Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) 〈X² (t)〉 ⋍ tᵅ with 1 < α < 2, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion (0 < α < 1) this behaviour is reversed and the particle density is depleted close to the boundaries. The MSD in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers.}, language = {en} } @article{GuggenbergerPagniniVojtaetal.2019, author = {Guggenberger, Tobias and Pagnini, Gianni and Vojta, Thomas and Metzler, Ralf}, title = {Fractional Brownian motion in a finite interval}, series = {New Journal of Physics}, volume = {21}, journal = {New Journal of Physics}, publisher = {Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics}, address = {Bad Honnef und London}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab075f}, pages = {13}, year = {2019}, abstract = {Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) 〈X² (t)〉 ⋍ tᵅ with 1 < α < 2, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion (0 < α < 1) this behaviour is reversed and the particle density is depleted close to the boundaries. The MSD in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers.}, language = {en} } @misc{SposiniMetzlerOshanin2019, author = {Sposini, Vittoria and Metzler, Ralf and Oshanin, Gleb}, title = {Single-trajectory spectral analysis of scaled Brownian motion}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {753}, issn = {1866-8372}, doi = {10.25932/publishup-43652}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-436522}, pages = {16}, year = {2019}, abstract = {Astandard approach to study time-dependent stochastic processes is the power spectral density (PSD), an ensemble-averaged property defined as the Fourier transform of the autocorrelation function of the process in the asymptotic limit of long observation times, T → ∞. In many experimental situations one is able to garner only relatively few stochastic time series of finite T, such that practically neither an ensemble average nor the asymptotic limit T → ∞ can be achieved. To accommodate for a meaningful analysis of such finite-length data we here develop the framework of single-trajectory spectral analysis for one of the standard models of anomalous diffusion, scaled Brownian motion.Wedemonstrate that the frequency dependence of the single-trajectory PSD is exactly the same as for standard Brownian motion, which may lead one to the erroneous conclusion that the observed motion is normal-diffusive. However, a distinctive feature is shown to be provided by the explicit dependence on the measurement time T, and this ageing phenomenon can be used to deduce the anomalous diffusion exponent.Wealso compare our results to the single-trajectory PSD behaviour of another standard anomalous diffusion process, fractional Brownian motion, and work out the commonalities and differences. Our results represent an important step in establishing singletrajectory PSDs as an alternative (or complement) to analyses based on the time-averaged mean squared displacement.}, language = {en} }