@article{GrebenkovMetzlerOshanin2021,
  author    = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb},
  title     = {A molecular relay race: sequential first-passage events to the terminal reaction centre in a cascade of diffusion controlled processes},
  series = {New Journal of Physics (NJP)},
  volume    = {23},
  journal   = {New Journal of Physics (NJP)},
  publisher = {IOP - Institute of Physics Publishing},
  address   = {Bristol},
  issn      = {1367-2630},
  doi       = {10.1088/1367-2630/ac1e42},
  pages     = {18},
  year      = {2021},
  abstract  = {We consider a sequential cascade of molecular first-reaction events towards a terminal reaction centre in which each reaction step is controlled by diffusive motion of the particles. The model studied here represents a typical reaction setting encountered in diverse molecular biology systems, in which, e.g. a signal transduction proceeds via a series of consecutive 'messengers': the first messenger has to find its respective immobile target site triggering a launch of the second messenger, the second messenger seeks its own target site and provokes a launch of the third messenger and so on, resembling a relay race in human competitions. For such a molecular relay race taking place in infinite one-, two- and three-dimensional systems, we find exact expressions for the probability density function of the time instant of the terminal reaction event, conditioned on preceding successful reaction events on an ordered array of target sites. The obtained expressions pertain to the most general conditions: number of intermediate stages and the corresponding diffusion coefficients, the sizes of the target sites, the distances between them, as well as their reactivities are arbitrary.},
  language  = {en}
}
@misc{GrebenkovMetzlerOshanin2021,
  author    = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb},
  title     = {A molecular relay race: sequential first-passage events to the terminal reaction centre in a cascade of diffusion controlled processes},
  series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-52194},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-521942},
  pages     = {20},
  year      = {2021},
  abstract  = {We consider a sequential cascade of molecular first-reaction events towards a terminal reaction centre in which each reaction step is controlled by diffusive motion of the particles. The model studied here represents a typical reaction setting encountered in diverse molecular biology systems, in which, e.g. a signal transduction proceeds via a series of consecutive 'messengers': the first messenger has to find its respective immobile target site triggering a launch of the second messenger, the second messenger seeks its own target site and provokes a launch of the third messenger and so on, resembling a relay race in human competitions. For such a molecular relay race taking place in infinite one-, two- and three-dimensional systems, we find exact expressions for the probability density function of the time instant of the terminal reaction event, conditioned on preceding successful reaction events on an ordered array of target sites. The obtained expressions pertain to the most general conditions: number of intermediate stages and the corresponding diffusion coefficients, the sizes of the target sites, the distances between them, as well as their reactivities are arbitrary.},
  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{GuggenbergerChechkinMetzler2022,
  author    = {Guggenberger, Tobias and Chechkin, Aleksei and Metzler, Ralf},
  title     = {Absence of stationary states and non-Boltzmann distributions of fractional Brownian motion in shallow external potentials},
  series = {New journal of physics : the open-access journal for physics},
  volume    = {24},
  journal   = {New journal of physics : the open-access journal for physics},
  number    = {7},
  publisher = {Dt. Physikalische Ges.},
  address   = {[Bad Honnef]},
  issn      = {1367-2630},
  doi       = {10.1088/1367-2630/ac7b3c},
  pages     = {18},
  year      = {2022},
  abstract  = {We study the diffusive motion of a particle in a subharmonic potential of the form U(x) = |x|( c ) (0 < c < 2) driven by long-range correlated, stationary fractional Gaussian noise xi ( alpha )(t) with 0 < alpha <= 2. In the absence of the potential the particle exhibits free fractional Brownian motion with anomalous diffusion exponent alpha. While for an harmonic external potential the dynamics converges to a Gaussian stationary state, from extensive numerical analysis we here demonstrate that stationary states for shallower than harmonic potentials exist only as long as the relation c > 2(1 - 1/alpha) holds. We analyse the motion in terms of the mean squared displacement and (when it exists) the stationary probability density function. Moreover we discuss analogies of non-stationarity of Levy flights in shallow external potentials.},
  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}
}
@article{MutothyaXuLietal.2021,
  author    = {Mutothya, Nicholas Mwilu and Xu, Yong and Li, Yongge and Metzler, Ralf},
  title     = {Characterising stochastic motion in heterogeneous media driven by coloured non-Gaussian noise},
  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/abfba6},
  pages     = {31},
  year      = {2021},
  abstract  = {We study the stochastic motion of a test particle in a heterogeneous medium in terms of a position dependent diffusion coefficient mimicking measured deterministic diffusivity gradients in biological cells or the inherent heterogeneity of geophysical systems. Compared to previous studies we here investigate the effect of the interplay of anomalous diffusion effected by position dependent diffusion coefficients and coloured non-Gaussian noise. The latter is chosen to be distributed according to Tsallis' q-distribution, representing a popular example for a non-extensive statistic. We obtain the ensemble and time averaged mean squared displacements for this generalised process and establish its non-ergodic properties as well as analyse the non-Gaussian nature of the associated displacement distribution. We consider both non-stratified and stratified environments.},
  language  = {en}
}
@article{VilkAghionNathanetal.2022,
  author    = {Vilk, Ohad and Aghion, Erez and Nathan, Ran and Toledo, Sivan and Metzler, Ralf and Assaf, Michael},
  title     = {Classification of anomalous diffusion in animal movement data using power spectral analysis},
  series = {Journal of physics : A, Mathematical and theoretical},
  volume    = {55},
  journal   = {Journal of physics : A, Mathematical and theoretical},
  number    = {33},
  publisher = {IOP Publishing},
  address   = {Bristol},
  issn      = {1751-8113},
  doi       = {10.1088/1751-8121/ac7e8f},
  pages     = {16},
  year      = {2022},
  abstract  = {The field of movement ecology has seen a rapid increase in high-resolution data in recent years, leading to the development of numerous statistical and numerical methods to analyse relocation trajectories. Data are often collected at the level of the individual and for long periods that may encompass a range of behaviours. Here, we use the power spectral density (PSD) to characterise the random movement patterns of a black-winged kite (Elanus caeruleus) and a white stork (Ciconia ciconia). The tracks are first segmented and clustered into different behaviours (movement modes), and for each mode we measure the PSD and the ageing properties of the process. For the foraging kite we find 1/f noise, previously reported in ecological systems mainly in the context of population dynamics, but not for movement data. We further suggest plausible models for each of the behavioural modes by comparing both the measured PSD exponents and the distribution of the single-trajectory PSD to known theoretical results and simulations.},
  language  = {en}
}
@misc{ŚlęzakMetzlerMagdziarz2019,
  author    = {Ślęzak, Jakub and Metzler, Ralf and Magdziarz, Marcin},
  title     = {Codifference can detect ergodicity breaking and non-Gaussianity},
  series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  number    = {748},
  doi       = {10.25932/publishup-43617},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-436178},
  pages     = {25},
  year      = {2019},
  abstract  = {We show that the codifference is a useful tool in studying the ergodicity breaking and non-Gaussianity properties of stochastic time series. While the codifference is a measure of dependence that was previously studied mainly in the context of stable processes, we here extend its range of applicability to random-parameter and diffusing-diffusivity models which are important in contemporary physics, biology and financial engineering. We prove that the codifference detects forms of dependence and ergodicity breaking which are not visible from analysing the covariance and correlation functions. We also discuss a related measure of dispersion, which is a nonlinear analogue of the mean squared displacement.},
  language  = {en}
}
@article{ŚlęzakMetzlerMagdziarz2019,
  author    = {Ślęzak, Jakub and Metzler, Ralf and Magdziarz, Marcin},
  title     = {Codifference can detect ergodicity breaking and non-Gaussianity},
  series = {New Journal of Physics},
  volume    = {21},
  journal   = {New Journal of Physics},
  publisher = {Deutsche Physikalische Gesellschaft},
  address   = {Bad Honnef},
  issn      = {1367-2630},
  doi       = {10.1088/1367-2630/ab13f3},
  pages     = {25},
  year      = {2019},
  abstract  = {We show that the codifference is a useful tool in studying the ergodicity breaking and non-Gaussianity properties of stochastic time series. While the codifference is a measure of dependence that was previously studied mainly in the context of stable processes, we here extend its range of applicability to random-parameter and diffusing-diffusivity models which are important in contemporary physics, biology and financial engineering. We prove that the codifference detects forms of dependence and ergodicity breaking which are not visible from analysing the covariance and correlation functions. We also discuss a related measure of dispersion, which is a nonlinear analogue of the mean squared displacement.},
  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{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}
}
@article{MagdziarzMetzlerSzczotkaetal.2012,
  author    = {Magdziarz, Marcin and Metzler, Ralf and Szczotka, Wladyslaw and Zebrowski, Piotr},
  title     = {Correlated continuous-time random walks-scaling limits and Langevin picture},
  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/2012/04/P04010},
  pages     = {18},
  year      = {2012},
  abstract  = {In this paper we analyze correlated continuous-time random walks introduced recently by Tejedor and Metzler (2010 J. Phys. A: Math. Theor. 43 082002). We obtain the Langevin equations associated with this process and the corresponding scaling limits of their solutions. We prove that the limit processes are self-similar and display anomalous dynamics. Moreover, we extend the model to include external forces. Our results are confirmed by Monte Carlo simulations.},
  language  = {en}
}
@misc{GrebenkovMetzlerOshanin2021,
  author    = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb},
  title     = {Distribution of first-reaction times with target regions on boundaries of shell-like domains},
  series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  publisher = {Universit{\"a}tsverlag Potsdam},
  address   = {Potsdam},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-55754},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-557542},
  pages     = {1 -- 23},
  year      = {2021},
  abstract  = {We study the probability density function (PDF) of the first-reaction times between a diffusive ligand and a membrane-bound, immobile imperfect target region in a restricted 'onion-shell' geometry bounded by two nested membranes of arbitrary shapes. For such a setting, encountered in diverse molecular signal transduction pathways or in the narrow escape problem with additional steric constraints, we derive an exact spectral form of the PDF, as well as present its approximate form calculated by help of the so-called self-consistent approximation. For a particular case when the nested domains are concentric spheres, we get a fully explicit form of the approximated PDF, assess the accuracy of this approximation, and discuss various facets of the obtained distributions. Our results can be straightforwardly applied to describe the PDF of the terminal reaction event in multi-stage signal transduction processes.},
  language  = {en}
}
@article{GrebenkovMetzlerOshanin2021,
  author    = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb},
  title     = {Distribution of first-reaction times with target regions on boundaries of shell-like domains},
  series = {New Journal of Physics (NJP)},
  volume    = {2021},
  journal   = {New Journal of Physics (NJP)},
  edition   = {23},
  publisher = {IOP Publishing},
  address   = {London},
  issn      = {1367-2630},
  doi       = {10.1088/1367-2630/ac4282},
  pages     = {1 -- 23},
  year      = {2021},
  abstract  = {We study the probability density function (PDF) of the first-reaction times between a diffusive ligand and a membrane-bound, immobile imperfect target region in a restricted 'onion-shell' geometry bounded by two nested membranes of arbitrary shapes. For such a setting, encountered in diverse molecular signal transduction pathways or in the narrow escape problem with additional steric constraints, we derive an exact spectral form of the PDF, as well as present its approximate form calculated by help of the so-called self-consistent approximation. For a particular case when the nested domains are concentric spheres, we get a fully explicit form of the approximated PDF, assess the accuracy of this approximation, and discuss various facets of the obtained distributions. Our results can be straightforwardly applied to describe the PDF of the terminal reaction event in multi-stage signal transduction processes.},
  language  = {en}
}
@article{CherstvyMetzler2015,
  author    = {Cherstvy, Andrey G. and Metzler, Ralf},
  title     = {Ergodicity breaking, ageing, and confinement in generalized diffusion processes with position and time dependent diffusivity},
  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/2015/05/P05010},
  pages     = {20},
  year      = {2015},
  abstract  = {We study generalized anomalous diffusion processes whose diffusion coefficient D(x, t) similar to D-0x(alpha)t(beta) depends on both the position x of the test particle and the process time t. This process thus combines the features of scaled Brownian motion and heterogeneous diffusion parent processes. We compute the ensemble and time averaged mean squared displacements of this generalized diffusion process. The scaling exponent of the ensemble averaged mean squared displacement is shown to be the product of the critical exponents of the parent processes, and describes both subdiffusive and superdiffusive systems. We quantify the amplitude fluctuations of the time averaged mean squared displacement as function of the length of the time series and the lag time. In particular, we observe a weak ergodicity breaking of this generalized diffusion process: even in the long time limit the ensemble and time averaged mean squared displacements are strictly disparate. When we start to observe this process some time after its initiation we observe distinct features of ageing. We derive a universal ageing factor for the time averaged mean squared displacement containing all information on the ageing time and the measurement time. External confinement is shown to alter the magnitudes and statistics of the ensemble and time averaged mean squared displacements.},
  language  = {en}
}
@article{GrebenkovSposiniMetzleretal.2020,
  author    = {Grebenkov, Denis S. and Sposini, Vittoria and Metzler, Ralf and Oshanin, Gleb and Seno, Flavio},
  title     = {Exact distributions of the maximum and range of random diffusivity processes},
  series = {New Journal of Physics},
  volume    = {23},
  journal   = {New Journal of Physics},
  publisher = {Dt. Physikalische Ges.},
  address   = {Bad Honnef},
  issn      = {1367-2630},
  doi       = {10.1088/1367-2630/abd313},
  pages     = {23},
  year      = {2020},
  abstract  = {We study the extremal properties of a stochastic process xt defined by the Langevin equation ẋₜ =√2Dₜ ξₜ, in which ξt is a Gaussian white noise with zero mean and Dₜ is a stochastic'diffusivity', defined as a functional of independent Brownian motion Bₜ.We focus on threechoices for the random diffusivity Dₜ: cut-off Brownian motion, Dₜt ∼ Θ(Bₜ), where Θ(x) is the Heaviside step function; geometric Brownian motion, Dₜ ∼ exp(-Bₜ); and a superdiffusive process based on squared Brownian motion, Dₜ ∼ B²ₜ. For these cases we derive exact expressions for the probability density functions of the maximal positive displacement and of the range of the process xₜ on the time interval ₜ ∈ (0, T).We discuss the asymptotic behaviours of the associated probability density functions, compare these against the behaviour of the corresponding properties of standard Brownian motion with constant diffusivity (Dₜ = D0) and also analyse the typical behaviour of the probability density functions which is observed for a majority of realisations of the stochastic diffusivity process.},
  language  = {en}
}
@misc{GrebenkovSposiniMetzleretal.2020,
  author    = {Grebenkov, Denis S. and Sposini, Vittoria and Metzler, Ralf and Oshanin, Gleb and Seno, Flavio},
  title     = {Exact distributions of the maximum and range 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    = {1142},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-50397},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-503976},
  pages     = {24},
  year      = {2020},
  abstract  = {We study the extremal properties of a stochastic process xt defined by the Langevin equation ẋₜ =√2Dₜ ξₜ, in which ξt is a Gaussian white noise with zero mean and Dₜ is a stochastic'diffusivity', defined as a functional of independent Brownian motion Bₜ.We focus on threechoices for the random diffusivity Dₜ: cut-off Brownian motion, Dₜt ∼ Θ(Bₜ), where Θ(x) is the Heaviside step function; geometric Brownian motion, Dₜ ∼ exp(-Bₜ); and a superdiffusive process based on squared Brownian motion, Dₜ ∼ B²ₜ. For these cases we derive exact expressions for the probability density functions of the maximal positive displacement and of the range of the process xₜ on the time interval ₜ ∈ (0, T).We discuss the asymptotic behaviours of the associated probability density functions, compare these against the behaviour of the corresponding properties of standard Brownian motion with constant diffusivity (Dₜ = D0) and also analyse the typical behaviour of the probability density functions which is observed for a majority of realisations of the stochastic diffusivity process.},
  language  = {en}
}
@article{ShinCherstvyKimetal.2015,
  author    = {Shin, Jaeoh and Cherstvy, Andrey G. and Kim, Won Kyu and Metzler, Ralf},
  title     = {Facilitation of polymer looping and giant polymer diffusivity in crowded solutions of active particles},
  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/11/113008},
  pages     = {12},
  year      = {2015},
  abstract  = {We study the dynamics of polymer chains in a bath of self-propelled particles (SPP) by extensive Langevin dynamics simulations in a two-dimensional model system. Specifically, we analyse the polymer looping properties versus the SPP activity and investigate how the presence of the active particles alters the chain conformational statistics. We find that SPPs tend to extend flexible polymer chains, while they rather compactify stiffer semiflexible polymers, in agreement with previous results. Here we show that higher activities of SPPs yield a higher effective temperature of the bath and thus facilitate the looping kinetics of a passive polymer chain. We explicitly compute the looping probability and looping time in a wide range of the model parameters. We also analyse the motion of a monomeric tracer particle and the polymer's centre of mass in the presence of the active particles in terms of the time averaged mean squared displacement, revealing a giant diffusivity enhancement for the polymer chain via SPP pooling. Our results are applicable to rationalising the dimensions and looping kinetics of biopolymers at constantly fluctuating and often actively driven conditions inside biological cells or in suspensions of active colloidal particles or bacteria cells.},
  language  = {en}
}
@article{MutothyaXuLietal.2021,
  author    = {Mutothya, Nicholas Mwilu and Xu, Yong and Li, Yongge and Metzler, Ralf and Mutua, Nicholas Muthama},
  title     = {First passage dynamics of stochastic motion in heterogeneous media driven by correlated white Gaussian and coloured non-Gaussian noises},
  series = {Journal of physics. Complexity},
  volume    = {2},
  journal   = {Journal of physics. Complexity},
  publisher = {IOP Publishing},
  address   = {Bristol},
  issn      = {2632-072X},
  doi       = {10.1088/2632-072X/ac35b5},
  pages     = {24},
  year      = {2021},
  abstract  = {We study the first passage dynamics for a diffusing particle experiencing a spatially varying diffusion coefficient while driven by correlated additive Gaussian white noise and multiplicative coloured non-Gaussian noise. We consider three functional forms for position dependence of the diffusion coefficient: power-law, exponential, and logarithmic. The coloured non-Gaussian noise is distributed according to Tsallis' q-distribution. Tracks of the non-Markovian systems are numerically simulated by using the fourth-order Runge-Kutta algorithm and the first passage times (FPTs) are recorded. The FPT density is determined along with the mean FPT (MFPT). Effects of the noise intensity and self-correlation of the multiplicative noise, the intensity of the additive noise, the cross-correlation strength, and the non-extensivity parameter on the MFPT are discussed.},
  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}
}
@misc{GrebenkovMetzlerOshanin2020,
  author    = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb},
  title     = {From single-particle stochastic kinetics to macroscopic reaction rates},
  series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  number    = {1018},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-48405},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-484059},
  pages     = {29},
  year      = {2020},
  abstract  = {We consider the first-passage problem for N identical independent particles that are initially released uniformly in a finite domain Ω and then diffuse toward a reactive area Γ, which can be part of the outer boundary of Ω or a reaction centre in the interior of Ω. For both cases of perfect and partial reactions, we obtain the explicit formulas for the first two moments of the fastest first-passage time (fFPT), i.e., the time when the first out of the N particles reacts with Γ. Moreover, we investigate the full probability density of the fFPT. We discuss a significant role of the initial condition in the scaling of the average fFPT with the particle number N, namely, a much stronger dependence (1/N and 1/N² for partially and perfectly reactive targets, respectively), in contrast to the well known inverse-logarithmic behaviour found when all particles are released from the same fixed point. We combine analytic solutions with scaling arguments and stochastic simulations to rationalise our results, which open new perspectives for studying the relevance of multiple searchers in various situations of molecular reactions, in particular, in living cells.},
  language  = {en}
}
@article{GrebenkovMetzlerOshanin2020,
  author    = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb},
  title     = {From single-particle stochastic kinetics to macroscopic reaction rates},
  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/abb1de},
  pages     = {28},
  year      = {2020},
  abstract  = {We consider the first-passage problem for N identical independent particles that are initially released uniformly in a finite domain Ω and then diffuse toward a reactive area Γ, which can be part of the outer boundary of Ω or a reaction centre in the interior of Ω. For both cases of perfect and partial reactions, we obtain the explicit formulas for the first two moments of the fastest first-passage time (fFPT), i.e., the time when the first out of the N particles reacts with Γ. Moreover, we investigate the full probability density of the fFPT. We discuss a significant role of the initial condition in the scaling of the average fFPT with the particle number N, namely, a much stronger dependence (1/N and 1/N² for partially and perfectly reactive targets, respectively), in contrast to the well known inverse-logarithmic behaviour found when all particles are released from the same fixed point. We combine analytic solutions with scaling arguments and stochastic simulations to rationalise our results, which open new perspectives for studying the relevance of multiple searchers in various situations of molecular reactions, in particular, in living cells.},
  language  = {en}
}
@article{PalyulinMetzler2012,
  author    = {Palyulin, Vladimir V. and Metzler, Ralf},
  title     = {How a finite potential barrier decreases the mean first-passage time},
  series = {Journal of statistical mechanics: theory and experiment},
  journal   = {Journal of statistical mechanics: theory and experiment},
  number    = {1},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {1742-5468},
  doi       = {10.1088/1742-5468/2012/03/L03001},
  pages     = {10},
  year      = {2012},
  abstract  = {We consider the mean first-passage time of a random walker moving in a potential landscape on a finite interval, the starting and end points being at different potentials. From analytical calculations and Monte Carlo simulations we demonstrate that the mean first-passage time for a piecewise linear curve between these two points is minimized by the introduction of a potential barrier. Due to thermal fluctuations, this barrier may be crossed. It turns out that the corresponding expense for this activation is less severe than the gain from an increased slope towards the end point. In particular, the resulting mean first-passage time is shorter than for a linear potential drop between the two points.},
  language  = {en}
}
@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}
}
@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}
}
@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{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{Ś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}
}
@article{SinghMetzlerSandev2020,
  author    = {Singh, Rishu Kumar and Metzler, Ralf and Sandev, Trifce},
  title     = {Resetting dynamics in a confining potential},
  series = {Journal of physics : A, Mathematical and theoretical},
  volume    = {53},
  journal   = {Journal of physics : A, Mathematical and theoretical},
  number    = {50},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {1751-8113},
  doi       = {10.1088/1751-8121/abc83a},
  pages     = {28},
  year      = {2020},
  abstract  = {We study Brownian motion in a confining potential under a constant-rate resetting to a reset position x(0). The relaxation of this system to the steady-state exhibits a dynamic phase transition, and is achieved in a light cone region which grows linearly with time. When an absorbing boundary is introduced, effecting a symmetry breaking of the system, we find that resetting aids the barrier escape only when the particle starts on the same side as the barrier with respect to the origin. We find that the optimal resetting rate exhibits a continuous phase transition with critical exponent of unity. Exact expressions are derived for the mean escape time, the second moment, and the coefficient of variation (CV).},
  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}
}
@article{SposiniMetzlerOshanin2019,
  author    = {Sposini, Vittoria and Metzler, Ralf and Oshanin, Gleb},
  title     = {Single-trajectory spectral analysis of scaled Brownian motion},
  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/ab2f52},
  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}
}
@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{XuZhouMetzleretal.2022,
  author    = {Xu, Pengbo and Zhou, Tian and Metzler, Ralf and Deng, Weihua},
  title     = {Stochastic harmonic trapping of a L{\´e}vy walk},
  series = {New journal of physics : the open-access journal for physics / Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics},
  volume    = {24},
  journal   = {New journal of physics : the open-access journal for physics / Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics},
  number    = {3},
  publisher = {Deutsche Physikalische Gesellschaft},
  address   = {Bad Honnef},
  issn      = {1367-2630},
  doi       = {10.1088/1367-2630/ac5282},
  pages     = {1 -- 28},
  year      = {2022},
  abstract  = {We introduce and study a L{\´e}vy walk (LW) model of particle spreading with a finite propagation speed combined with soft resets, stochastically occurring periods in which an harmonic external potential is switched on and forces the particle towards a specific position. Soft resets avoid instantaneous relocation of particles that in certain physical settings may be considered unphysical. Moreover, soft resets do not have a specific resetting point but lead the particle towards a resetting point by a restoring Hookean force. Depending on the exact choice for the LW waiting time density and the probability density of the periods when the harmonic potential is switched on, we demonstrate a rich emerging response behaviour including ballistic motion and superdiffusion. When the confinement periods of the soft-reset events are dominant, we observe a particle localisation with an associated non-equilibrium steady state. In this case the stationary particle probability density function turns out to acquire multimodal states. Our derivations are based on Markov chain ideas and LWs with multiple internal states, an approach that may be useful and flexible for the investigation of other generalised random walks with soft and hard resets. The spreading efficiency of soft-rest LWs is characterised by the first-passage time statistic.},
  language  = {en}
}
@misc{XuZhouMetzleretal.2022,
  author    = {Xu, Pengbo and Zhou, Tian and Metzler, Ralf and Deng, Weihua},
  title     = {Stochastic harmonic trapping of a L{\´e}vy walk: transport and first-passage dynamics under soft resetting strategies},
  series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  publisher = {Universit{\"a}tsverlag Potsdam},
  address   = {Potsdam},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-56040},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-560402},
  pages     = {1 -- 28},
  year      = {2022},
  abstract  = {We introduce and study a L{\´e}vy walk (LW) model of particle spreading with a finite propagation speed combined with soft resets, stochastically occurring periods in which an harmonic external potential is switched on and forces the particle towards a specific position. Soft resets avoid instantaneous relocation of particles that in certain physical settings may be considered unphysical. Moreover, soft resets do not have a specific resetting point but lead the particle towards a resetting point by a restoring Hookean force. Depending on the exact choice for the LW waiting time density and the probability density of the periods when the harmonic potential is switched on, we demonstrate a rich emerging response behaviour including ballistic motion and superdiffusion. When the confinement periods of the soft-reset events are dominant, we observe a particle localisation with an associated non-equilibrium steady state. In this case the stationary particle probability density function turns out to acquire multimodal states. Our derivations are based on Markov chain ideas and LWs with multiple internal states, an approach that may be useful and flexible for the investigation of other generalised random walks with soft and hard resets. The spreading efficiency of soft-rest LWs is characterised by the first-passage time statistic.},
  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{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{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{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}
}