@article{ThapaWyłomańskaSikoraetal.2021, author = {Thapa, Samudrajit and Wyłomańska, Agnieszka and Sikora, Grzegorz and Wagner, Caroline E. and Krapf, Diego and Kantz, Holger and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories}, series = {New Journal of Physics}, volume = {23}, journal = {New Journal of Physics}, publisher = {Dt. Physikalische Ges. ; IOP}, address = {Bad Honnef ; London}, issn = {1367-2630}, doi = {10.1088/1367-2630/abd50e}, pages = {22}, year = {2021}, abstract = {Extensive time-series encoding the position of particles such as viruses, vesicles, or individualproteins are routinely garnered insingle-particle tracking experiments or supercomputing studies.They contain vital clues on how viruses spread or drugs may be delivered in biological cells.Similar time-series are being recorded of stock values in financial markets and of climate data.Such time-series are most typically evaluated in terms of time-averaged mean-squareddisplacements (TAMSDs), which remain random variables for finite measurement times. Theirstatistical properties are different for differentphysical stochastic processes, thus allowing us toextract valuable information on the stochastic process itself. To exploit the full potential of thestatistical information encoded in measured time-series we here propose an easy-to-implementand computationally inexpensive new methodology, based on deviations of the TAMSD from itsensemble average counterpart. Specifically, we use the upper bound of these deviations forBrownian motion (BM) to check the applicability of this approach to simulated and real data sets.By comparing the probability of deviations fordifferent data sets, we demonstrate how thetheoretical bound for BM reveals additional information about observed stochastic processes. Weapply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracerbeads measured in mucin hydrogels, and of geographic surface temperature anomalies. Ouranalysis shows how the large-deviation properties can be efficiently used as a simple yet effectiveroutine test to reject the BM hypothesis and unveil relevant information on statistical propertiessuch as ergodicity breaking and short-time correlations.}, language = {en} } @article{SposiniChechkinMetzler2018, author = {Sposini, Vittoria and Chechkin, Aleksei V. and Metzler, Ralf}, title = {First passage statistics for diffusing diffusivity}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {52}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {4}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aaf6ff}, pages = {11}, year = {2018}, abstract = {A rapidly increasing number of systems is identified in which the stochastic motion of tracer particles follows the Brownian law < r(2)(t)> similar or equal to Dt yet the distribution of particle displacements is strongly non-Gaussian. A central approach to describe this effect is the diffusing diffusivity (DD) model in which the diffusion coefficient itself is a stochastic quantity, mimicking heterogeneities of the environment encountered by the tracer particle on its path. We here quantify in terms of analytical and numerical approaches the first passage behaviour of the DD model. We observe significant modifications compared to Brownian-Gaussian diffusion, in particular that the DD model may have a faster first passage dynamics. Moreover we find a universal crossover point of the survival probability independent of the initial condition.}, language = {en} } @article{DybiecCapalaChechkinetal.2018, author = {Dybiec, Bartlomiej and Capala, Karol and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Conservative random walks in confining potentials}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {52}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aaefc2}, pages = {25}, year = {2018}, abstract = {Levy walks are continuous time random walks with spatio-temporal coupling of jump lengths and waiting times, often used to model superdiffusive spreading processes such as animals searching for food, tracer motion in weakly chaotic systems, or even the dynamics in quantum systems such as cold atoms. In the simplest version Levy walks move with a finite speed. Here, we present an extension of the Levy walk scenario for the case when external force fields influence the motion. The resulting motion is a combination of the response to the deterministic force acting on the particle, changing its velocity according to the principle of total energy conservation, and random velocity reversals governed by the distribution of waiting times. For the fact that the motion stays conservative, that is, on a constant energy surface, our scenario is fundamentally different from thermal motion in the same external potentials. In particular, we present results for the velocity and position distributions for single well potentials of different steepness. The observed dynamics with its continuous velocity changes enriches the theory of Levy walk processes and will be of use in a variety of systems, for which the particles are externally confined.}, language = {en} } @article{SandevTomovskiDubbeldametal.2018, author = {Sandev, Trifce and Tomovski, Zivorad and Dubbeldam, Johan L. A. and Chechkin, Aleksei V.}, title = {Generalized diffusion-wave equation with memory kernel}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {52}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aaefa3}, pages = {22}, year = {2018}, abstract = {We study generalized diffusion-wave equation in which the second order time derivative is replaced by an integro-differential operator. It yields time fractional and distributed order time fractional diffusion-wave equations as particular cases. We consider different memory kernels of the integro-differential operator, derive corresponding fundamental solutions, specify the conditions of their non-negativity and calculate the mean squared displacement for all cases. In particular, we introduce and study generalized diffusion-wave equations with a regularized Prabhakar derivative of single and distributed orders. The equations considered can be used for modeling the broad spectrum of anomalous diffusion processes and various transitions between different diffusion regimes.}, language = {en} } @article{GajdaWylomanskaKantzetal.2018, author = {Gajda, J. and Wylomanska, Agnieszka and Kantz, Holger and Chechkin, Aleksei V. and Sikora, Grzegorz}, title = {Large deviations of time-averaged statistics for Gaussian processes}, series = {Statistics \& Probability Letters}, volume = {143}, journal = {Statistics \& Probability Letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0167-7152}, doi = {10.1016/j.spl.2018.07.013}, pages = {47 -- 55}, year = {2018}, abstract = {In this paper we study the large deviations of time averaged mean square displacement (TAMSD) for Gaussian processes. The theory of large deviations is related to the exponential decay of probabilities of large fluctuations in random systems. From the mathematical point of view a given statistics satisfies the large deviation principle, if the probability that it belongs to a certain range decreases exponentially. The TAMSD is one of the main statistics used in the problem of anomalous diffusion detection. Applying the theory of generalized chi-squared distribution and sub-gamma random variables we prove the upper bound for large deviations of TAMSD for Gaussian processes. As a special case we consider fractional Brownian motion, one of the most popular models of anomalous diffusion. Moreover, we derive the upper bound for large deviations of the estimator for the anomalous diffusion exponent. (C) 2018 Elsevier B.V. All rights reserved.}, language = {en} } @article{MardoukhiChechkinMetzler2020, author = {Mardoukhi, Yousof and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Spurious ergodicity breaking in normal and fractional Ornstein-Uhlenbeck process}, series = {New Journal of Physics}, volume = {22}, journal = {New Journal of Physics}, publisher = {IOP}, address = {London}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab950b}, pages = {18}, year = {2020}, abstract = {The Ornstein-Uhlenbeck process is a stationary and ergodic Gaussian process, that is fully determined by its covariance function and mean. We show here that the generic definitions of the ensemble- and time-averaged mean squared displacements fail to capture these properties consistently, leading to a spurious ergodicity breaking. We propose to remedy this failure by redefining the mean squared displacements such that they reflect unambiguously the statistical properties of any stochastic process. In particular we study the effect of the initial condition in the Ornstein-Uhlenbeck process and its fractional extension. For the fractional Ornstein-Uhlenbeck process representing typical experimental situations in crowded environments such as living biological cells, we show that the stationarity of the process delicately depends on the initial condition.}, language = {en} } @article{DelleSideNassisiPennettaetal.2017, author = {Delle Side, Domenico and Nassisi, Vincenzo and Pennetta, Cecilia and Alifano, Pietro and Di Salvo, Marco and Tala, Adelfia and Chechkin, Aleksei V. and Seno, Flavio and Trovato, Antonio}, title = {Bacterial bioluminescence onset and quenching: a dynamical model for a quorum sensing-mediated property}, series = {Royal Society Open Science}, volume = {4}, journal = {Royal Society Open Science}, publisher = {Royal Society}, address = {London}, issn = {2054-5703}, doi = {10.1098/rsos.171586}, pages = {12}, year = {2017}, abstract = {We present an effective dynamical model for the onset of bacterial bioluminescence, one of the most studied quorum sensing-mediated traits. Our model is built upon simple equations that describe the growth of the bacterial colony, the production and accumulation of autoinducer signal molecules, their sensing within bacterial cells, and the ensuing quorum activation mechanism that triggers bioluminescent emission. The model is directly tested to quantitatively reproduce the experimental distributions of photon emission times, previously measured for bacterial colonies of Vibrio jasicida, a luminescent bacterium belonging to the Harveyi clade, growing in a highly drying environment. A distinctive and novel feature of the proposed model is bioluminescence 'quenching' after a given time elapsed from activation. Using an advanced fitting procedure based on the simulated annealing algorithm, we are able to infer from the experimental observations the biochemical parameters used in the model. Such parameters are in good agreement with the literature data. As a further result, we find that, at least in our experimental conditions, light emission in bioluminescent bacteria appears to originate from a subtle balance between colony growth and quorum activation due to autoinducers diffusion, with the two phenomena occurring on the same time scale. This finding is consistent with a negative feedback mechanism previously reported for Vibrio harveyi.}, language = {en} } @article{PalyulinMantsevichKlagesetal.2017, author = {Palyulin, Vladimir V. and Mantsevich, Vladimir N. and Klages, Rainer and Metzler, Ralf and Chechkin, Aleksei V.}, title = {Comparison of pure and combined search strategies for single and multiple targets}, series = {The European physical journal : B, Condensed matter and complex systems}, volume = {90}, journal = {The European physical journal : B, Condensed matter and complex systems}, publisher = {Springer}, address = {New York}, issn = {1434-6028}, doi = {10.1140/epjb/e2017-80372-4}, pages = {20 -- 37}, year = {2017}, abstract = {We address the generic problem of random search for a point-like target on a line. Using the measures of search reliability and efficiency to quantify the random search quality, we compare Brownian search with Levy search based on long-tailed jump length distributions. We then compare these results with a search process combined of two different long-tailed jump length distributions. Moreover, we study the case of multiple targets located by a Levy searcher.}, language = {en} } @article{PalyulinBlackburnLomholtetal.2019, author = {Palyulin, Vladimir V. and Blackburn, George and Lomholt, Michael A. and Watkins, Nicholas W. and Metzler, Ralf and Klages, Rainer and Chechkin, Aleksei V.}, title = {First passage and first hitting times of Levy flights and Levy walks}, series = {New journal of physics : the open-access journal for physics}, volume = {21}, journal = {New journal of physics : the open-access journal for physics}, number = {10}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab41bb}, pages = {23}, year = {2019}, abstract = {For both L{\´e}vy flight and L{\´e}vy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For L{\´e}vy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the L{\´e}vy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.}, language = {en} } @article{CherstvyThapaMardoukhietal.2018, author = {Cherstvy, Andrey G. and Thapa, Samudrajit and Mardoukhi, Yousof and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Time averages and their statistical variation for the Ornstein-Uhlenbeck process}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {98}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {2}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.98.022134}, pages = {15}, year = {2018}, abstract = {How ergodic is diffusion under harmonic confinements? How strongly do ensemble- and time-averaged displacements differ for a thermally-agitated particle performing confined motion for different initial conditions? We here study these questions for the generic Ornstein-Uhlenbeck (OU) process and derive the analytical expressions for the second and fourth moment. These quantifiers are particularly relevant for the increasing number of single-particle tracking experiments using optical traps. For a fixed starting position, we discuss the definitions underlying the ensemble averages. We also quantify effects of equilibrium and nonequilibrium initial particle distributions onto the relaxation properties and emerging nonequivalence of the ensemble- and time-averaged displacements (even in the limit of long trajectories). We derive analytical expressions for the ergodicity breaking parameter quantifying the amplitude scatter of individual time-averaged trajectories, both for equilibrium and outof-equilibrium initial particle positions, in the entire range of lag times. Our analytical predictions are in excellent agreement with results of computer simulations of the Langevin equation in a parabolic potential. We also examine the validity of the Einstein relation for the ensemble- and time-averaged moments of the OU-particle. Some physical systems, in which the relaxation and nonergodic features we unveiled may be observable, are discussed.}, language = {en} }