TY - JOUR A1 - Aydiner, Ekrem A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Wealth distribution, Pareto law, and stretched exponential decay of money BT - Computer simulations analysis of agent-based models JF - Physica : europhysics journal ; A, Statistical mechanics and its applications N2 - We study by Monte Carlo simulations a kinetic exchange trading model for both fixed and distributed saving propensities of the agents and rationalize the person and wealth distributions. We show that the newly introduced wealth distribution – that may be more amenable in certain situations – features a different power-law exponent, particularly for distributed saving propensities of the agents. For open agent-based systems, we analyze the person and wealth distributions and find that the presence of trap agents alters their amplitude, leaving however the scaling exponents nearly unaffected. For an open system, we show that the total wealth – for different trap agent densities and saving propensities of the agents – decreases in time according to the classical Kohlrausch–Williams–Watts stretched exponential law. Interestingly, this decay does not depend on the trap agent density, but rather on saving propensities. The system relaxation for fixed and distributed saving schemes are found to be different. KW - Econophysics KW - Wealth and income distribution KW - Pareto law KW - Scaling exponents Y1 - 2017 U6 - https://doi.org/10.1016/j.physa.2017.08.017 SN - 0378-4371 SN - 1873-2119 VL - 490 SP - 278 EP - 288 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Pulkkinen, Otto A1 - Metzler, Ralf T1 - Variance-corrected Michaelis-Menten equation predicts transient rates of single-enzyme reactions and response times in bacterial gene-regulation JF - Scientific reports N2 - Many chemical reactions in biological cells occur at very low concentrations of constituent molecules. Thus, transcriptional gene-regulation is often controlled by poorly expressed transcription-factors, such as E.coli lac repressor with few tens of copies. Here we study the effects of inherent concentration fluctuations of substrate-molecules on the seminal Michaelis-Menten scheme of biochemical reactions. We present a universal correction to the Michaelis-Menten equation for the reaction-rates. The relevance and validity of this correction for enzymatic reactions and intracellular gene-regulation is demonstrated. Our analytical theory and simulation results confirm that the proposed variance-corrected Michaelis-Menten equation predicts the rate of reactions with remarkable accuracy even in the presence of large non-equilibrium concentration fluctuations. The major advantage of our approach is that it involves only the mean and variance of the substrate-molecule concentration. Our theory is therefore accessible to experiments and not specific to the exact source of the concentration fluctuations. Y1 - 2015 U6 - https://doi.org/10.1038/srep17820 SN - 2045-2322 IS - 5 PB - Nature Publishing Group CY - London ER - TY - JOUR A1 - Vilk, Ohad A1 - Aghion, Erez A1 - Avgar, Tal A1 - Beta, Carsten A1 - Nagel, Oliver A1 - Sabri, Adal A1 - Sarfati, Raphael A1 - Schwartz, Daniel K. A1 - Weiß, Matthias A1 - Krapf, Diego A1 - Nathan, Ran A1 - Metzler, Ralf A1 - Assaf, Michael T1 - Unravelling the origins of anomalous diffusion BT - from molecules to migrating storks JF - Physical review research / American Physical Society N2 - Anomalous diffusion or, more generally, anomalous transport, with nonlinear dependence of the mean-squared displacement on the measurement time, is ubiquitous in nature. It has been observed in processes ranging from microscopic movement of molecules to macroscopic, large-scale paths of migrating birds. Using data from multiple empirical systems, spanning 12 orders of magnitude in length and 8 orders of magnitude in time, we employ a method to detect the individual underlying origins of anomalous diffusion and transport in the data. This method decomposes anomalous transport into three primary effects: long-range correlations (“Joseph effect”), fat-tailed probability density of increments (“Noah effect”), and nonstationarity (“Moses effect”). We show that such a decomposition of real-life data allows us to infer nontrivial behavioral predictions and to resolve open questions in the fields of single-particle tracking in living cells and movement ecology. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevResearch.4.033055 SN - 2643-1564 VL - 4 IS - 3 PB - American Physical Society CY - College Park, MD ER - TY - JOUR A1 - Ritschel, Stefan A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Universality of delay-time averages for financial time series BT - analytical results, computer simulations, and analysis of historical stock-market prices JF - Journal of physics. Complexity N2 - We analyze historical data of stock-market prices for multiple financial indices using the concept of delay-time averaging for the financial time series (FTS). The region of validity of our recent theoretical predictions [Cherstvy A G et al 2017 New J. Phys. 19 063045] for the standard and delayed time-averaged mean-squared 'displacements' (TAMSDs) of the historical FTS is extended to all lag times. As the first novel element, we perform extensive computer simulations of the stochastic differential equation describing geometric Brownian motion (GBM) which demonstrate a quantitative agreement with the analytical long-term price-evolution predictions in terms of the delayed TAMSD (for all stock-market indices in crisis-free times). Secondly, we present a robust procedure of determination of the model parameters of GBM via fitting the features of the price-evolution dynamics in the FTS for stocks and cryptocurrencies. The employed concept of single-trajectory-based time averaging can serve as a predictive tool (proxy) for a mathematically based assessment and rationalization of probabilistic trends in the evolution of stock-market prices. KW - econophysics KW - geometric Brownian motion KW - time-series analysis Y1 - 2021 U6 - https://doi.org/10.1088/2632-072X/ac2220 SN - 2632-072X VL - 2 IS - 4 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Sposini, Vittoria A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb A1 - Seno, Flavio T1 - Universal spectral features of different classes of random-diffusivity processes JF - New Journal of Physics N2 - 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. KW - diffusion KW - power spectrum KW - random diffusivity KW - single trajectories Y1 - 2020 U6 - https://doi.org/10.1088/1367-2630/ab9200 SN - 1367-2630 VL - 22 IS - 6 PB - Dt. Physikalische Ges. CY - Bad Honnef ER - TY - JOUR A1 - Wang, Wei A1 - Seno, Flavio A1 - Sokolov, Igor M. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Unexpected crossovers in correlated random-diffusivity processes JF - New Journal of Physics N2 - 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. KW - diffusion KW - anomalous diffusion KW - non-Gaussianity KW - fractional Brownian motion Y1 - 2020 U6 - https://doi.org/10.1088/1367-2630/aba390 SN - 1367-2630 VL - 22 PB - Dt. Physikalische Ges. CY - Bad Honnef ER - TY - JOUR A1 - Bodrova, Anna S. A1 - Chechkin, Aleksei V. A1 - Cherstvy, Andrey G. A1 - Safdari, Hadiseh A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Underdamped scaled Brownian motion BT - (non-)existence of the overdamped limit in anomalous diffusion JF - Scientific reports N2 - It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases. Y1 - 2016 U6 - https://doi.org/10.1038/srep30520 SN - 2045-2322 VL - 6 PB - Nature Publishing Group CY - London ER - TY - JOUR A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Bodrova, Anna S. T1 - Ultraslow scaled Brownian motion JF - New journal of physics : the open-access journal for physics N2 - We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form . For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In a harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations. KW - anomalous diffusion KW - stochastic processes KW - ageing Y1 - 2015 U6 - https://doi.org/10.1088/1367-2630/17/6/063038 SN - 1367-2630 VL - 17 IS - 063038 PB - Dt. Physikalische Ges., IOP CY - Bad Honnef, London ER - TY - JOUR A1 - Petreska, Irina A1 - Pejov, Ljupco A1 - Sandev, Trifce A1 - Kocarev, Ljupčo A1 - Metzler, Ralf T1 - Tuning of the dielectric relaxation and complex susceptibility in a system of polar molecules: a generalised model based on rotational diffusion with resetting JF - Fractal and fractional N2 - The application of the fractional calculus in the mathematical modelling of relaxation processes in complex heterogeneous media has attracted a considerable amount of interest lately. The reason for this is the successful implementation of fractional stochastic and kinetic equations in the studies of non-Debye relaxation. In this work, we consider the rotational diffusion equation with a generalised memory kernel in the context of dielectric relaxation processes in a medium composed of polar molecules. We give an overview of existing models on non-exponential relaxation and introduce an exponential resetting dynamic in the corresponding process. The autocorrelation function and complex susceptibility are analysed in detail. We show that stochastic resetting leads to a saturation of the autocorrelation function to a constant value, in contrast to the case without resetting, for which it decays to zero. The behaviour of the autocorrelation function, as well as the complex susceptibility in the presence of resetting, confirms that the dielectric relaxation dynamics can be tuned by an appropriate choice of the resetting rate. The presented results are general and flexible, and they will be of interest for the theoretical description of non-trivial relaxation dynamics in heterogeneous systems composed of polar molecules. KW - rotational diffusion KW - memory kernel KW - Fokker-Planck equation KW - non-exponential relaxation KW - autocorrelation function KW - complex KW - susceptibility Y1 - 2022 U6 - https://doi.org/10.3390/fractalfract6020088 SN - 2504-3110 VL - 6 IS - 2 PB - MDPI AG, Fractal Fract Editorial Office CY - Basel ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb A1 - Dagdug, Leonardo A1 - Berezhkovskii, Alexander M. A1 - Skvortsov, Alexei T. T1 - Trapping of diffusing particles by periodic absorbing rings on a cylindrical tube JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr Y1 - 2019 U6 - https://doi.org/10.1063/1.5098390 SN - 0021-9606 SN - 1089-7690 VL - 150 IS - 20 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Teomy, Eial A1 - Metzler, Ralf T1 - Transport in exclusion processes with one-step memory: density dependence and optimal acceleration JF - Journal of physics : A, Mathematical and theoretical N2 - We study a lattice gas of persistent walkers, in which each site is occupied by at most one particle and the direction each particle attempts to move to depends on its last step. We analyse the mean squared displacement (MSD) of the particles as a function of the particle density and their persistence (the tendency to continue moving in the same direction). For positive persistence the MSD behaves as expected: it increases with the persistence and decreases with the density. However, for strong anti-persistence we find two different regimes, in which the dependence of the MSD on the density is non-monotonic. For very strong anti-persistence there is an optimal density at which the MSD reaches a maximum. In an intermediate regime, the MSD as a function of the density exhibits both a minimum and a maximum, a phenomenon which has not been observed before. We derive a mean-field theory which qualitatively explains this behaviour. KW - exclusion process KW - persistence KW - lattice gas Y1 - 2019 U6 - https://doi.org/10.1088/1751-8121/ab37e4 SN - 1751-8113 SN - 1751-8121 VL - 52 IS - 38 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Li, Hua A1 - Xu, Yong A1 - Li, Yongge A1 - Metzler, Ralf T1 - Transition path dynamics across rough inverted parabolic potential barrier JF - The European physical journal : Plus N2 - Transition path dynamics have been widely studied in chemical, physical, and technological systems. Mostly, the transition path dynamics is obtained for smooth barrier potentials, for instance, generic inverse-parabolic shapes. We here present analytical results for the mean transition path time, the distribution of transition path times, the mean transition path velocity, and the mean transition path shape in a rough inverted parabolic potential function under the driving of Gaussian white noise. These are validated against extensive simulations using the forward flux sampling scheme in parallel computations. We observe how precisely the potential roughness, the barrier height, and the noise intensity contribute to the particle transition in the rough inverted barrier potential. Y1 - 2020 U6 - https://doi.org/10.1140/epjp/s13360-020-00752-7 SN - 2190-5444 VL - 135 IS - 9 PB - Springer CY - Berlin ; Heidelberg ER - TY - JOUR A1 - Thapa, Samudrajit A1 - Lukat, Nils A1 - Selhuber-Unkel, Christine A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Transient superdiffusion of polydisperse vacuoles in highly motile amoeboid cells JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - We perform a detailed statistical analysis of diffusive trajectories of membrane-enclosed vesicles (vacuoles) in the supercrowded cytoplasm of living Acanthamoeba castellanii cells. From the vacuole traces recorded in the center-of-area frame of moving amoebae, we examine the statistics of the time-averaged mean-squared displacements of vacuoles, their generalized diffusion coefficients and anomalous scaling exponents, the ergodicity breaking parameter, the non-Gaussian features of displacement distributions of vacuoles, the displacement autocorrelation function, as well as the distributions of speeds and positions of vacuoles inside the amoeba cells. Our findings deliver novel insights into the internal dynamics of cellular structures in these infectious pathogens. Published under license by AIP Publishing. Y1 - 2019 U6 - https://doi.org/10.1063/1.5086269 SN - 0021-9606 SN - 1089-7690 VL - 150 IS - 14 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Sposini, Vittoria A1 - Krapf, Diego A1 - Marinari, Enzo A1 - Sunyer, Raimon A1 - Ritort, Felix A1 - Taheri, Fereydoon A1 - Selhuber-Unkel, Christine A1 - Benelli, Rebecca A1 - Weiss, Matthias A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Towards a robust criterion of anomalous diffusion JF - Communications Physics N2 - Anomalous-diffusion, the departure of the spreading dynamics of diffusing particles from the traditional law of Brownian-motion, is a signature feature of a large number of complex soft-matter and biological systems. Anomalous-diffusion emerges due to a variety of physical mechanisms, e.g., trapping interactions or the viscoelasticity of the environment. However, sometimes systems dynamics are erroneously claimed to be anomalous, despite the fact that the true motion is Brownian—or vice versa. This ambiguity in establishing whether the dynamics as normal or anomalous can have far-reaching consequences, e.g., in predictions for reaction- or relaxation-laws. Demonstrating that a system exhibits normal- or anomalous-diffusion is highly desirable for a vast host of applications. Here, we present a criterion for anomalous-diffusion based on the method of power-spectral analysis of single trajectories. The robustness of this criterion is studied for trajectories of fractional-Brownian-motion, a ubiquitous stochastic process for the description of anomalous-diffusion, in the presence of two types of measurement errors. In particular, we find that our criterion is very robust for subdiffusion. Various tests on surrogate data in absence or presence of additional positional noise demonstrate the efficacy of this method in practical contexts. Finally, we provide a proof-of-concept based on diverse experiments exhibiting both normal and anomalous-diffusion. Y1 - 2022 U6 - https://doi.org/10.1038/s42005-022-01079-8 SN - 2399-3650 VL - 5 PB - Springer Nature CY - London ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Towards a full quantitative description of single-molecule reaction kinetics in biological cells JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - The first-passage time (FPT), i.e., the moment when a stochastic process reaches a given threshold value for the first time, is a fundamental mathematical concept with immediate applications. In particular, it quantifies the statistics of instances when biomolecules in a biological cell reach their specific binding sites and trigger cellular regulation. Typically, the first-passage properties are given in terms of mean first-passage times. However, modern experiments now monitor single-molecular binding-processes in living cells and thus provide access to the full statistics of the underlying first-passage events, in particular, inherent cell-to-cell fluctuations. We here present a robust explicit approach for obtaining the distribution of FPTs to a small partially reactive target in cylindrical-annulus domains, which represent typical bacterial and neuronal cell shapes. We investigate various asymptotic behaviours of this FPT distribution and show that it is typically very broad in many biological situations, thus, the mean FPT can differ from the most probable FPT by orders of magnitude. The most probable FPT is shown to strongly depend only on the starting position within the geometry and to be almost independent of the target size and reactivity. These findings demonstrate the dramatic relevance of knowing the full distribution of FPTs and thus open new perspectives for a more reliable description of many intracellular processes initiated by the arrival of one or few biomolecules to a small, spatially localised region inside the cell. Y1 - 2018 U6 - https://doi.org/10.1039/c8cp02043d SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 24 SP - 16393 EP - 16401 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Vinod, Deepak A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf A1 - Sokolov, Igor M. T1 - Time-averaging and nonergodicity of reset geometric Brownian motion with drift JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - How do near-bankruptcy events in the past affect the dynamics of stock-market prices in the future? Specifically, what are the long-time properties of a time-local exponential growth of stock-market prices under the influence of stochastically occurring economic crashes? Here, we derive the ensemble- and time-averaged properties of the respective "economic" or geometric Brownian motion (GBM) with a nonzero drift exposed to a Poissonian constant-rate price-restarting process of "resetting." We examine-based both on thorough analytical calculations and on findings from systematic stochastic computer simulations-the general situation of reset GBM with a nonzero [positive] drift and for all special cases emerging for varying parameters of drift, volatility, and reset rate in the model. We derive and summarize all short- and long-time dependencies for the mean-squared displacement (MSD), the variance, and the mean time-averaged MSD (TAMSD) of the process of Poisson-reset GBM under the conditions of both rare and frequent resetting. We consider three main regions of model parameters and categorize the crossovers between different functional behaviors of the statistical quantifiers of this process. The analytical relations are fully supported by the results of computer simulations. In particular, we obtain that Poisson-reset GBM is a nonergodic stochastic process, with generally MSD(Delta) not equal TAMSD(Delta) and Variance(Delta) not equal TAMSD(Delta) at short lag times Delta and for long trajectory lengths T. We investigate the behavior of the ergodicity-breaking parameter in each of the three regions of parameters and examine its dependence on the rate of reset at Delta/T << 1. Applications of these theoretical results to the analysis of prices of reset-containing options are pertinent. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevE.106.034137 SN - 2470-0045 SN - 2470-0053 VL - 106 IS - 3 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Vinod, Deepak A1 - Aghion, Erez A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Time averaging, ageing and delay analysis of financial time series JF - New journal of physics N2 - 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. KW - time averaging KW - diffusion KW - geometric Brownian motion KW - financial time series Y1 - 2017 U6 - https://doi.org/10.1088/1367-2630/aa7199 SN - 1367-2630 VL - 19 SP - 1 EP - 11 PB - IOP CY - London ER - TY - JOUR A1 - Wang, Wei A1 - Cherstvy, Andrey G. A1 - Kantz, Holger A1 - Metzler, Ralf A1 - Sokolov, Igor M. T1 - Time averaging and emerging nonergodicity upon resetting of fractional Brownian motion and heterogeneous diffusion processes JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - How different are the results of constant-rate resetting of anomalous-diffusion processes in terms of their ensemble-averaged versus time-averaged mean-squared displacements (MSDs versus TAMSDs) and how does stochastic resetting impact nonergodicity? We examine, both analytically and by simulations, the implications of resetting on the MSD- and TAMSD-based spreading dynamics of particles executing fractional Brownian motion (FBM) with a long-time memory, heterogeneous diffusion processes (HDPs) with a power-law space-dependent diffusivity D(x) = D0|x|gamma and their "combined" process of HDP-FBM. We find, inter alia, that the resetting dynamics of originally ergodic FBM for superdiffusive Hurst exponents develops disparities in scaling and magnitudes of the MSDs and mean TAMSDs indicating weak ergodicity breaking. For subdiffusive HDPs we also quantify the nonequivalence of the MSD and TAMSD and observe a new trimodal form of the probability density function. For reset FBM, HDPs and HDP-FBM we compute analytically and verify by simulations the short-time MSD and TAMSD asymptotes and long-time plateaus reminiscent of those for processes under confinement. We show that certain characteristics of these reset processes are functionally similar despite a different stochastic nature of their nonreset variants. Importantly, we discover nonmonotonicity of the ergodicitybreaking parameter EB as a function of the resetting rate r. For all reset processes studied we unveil a pronounced resetting-induced nonergodicity with a maximum of EB at intermediate r and EB similar to(1/r )-decay at large r. Alongside the emerging MSD-versus-TAMSD disparity, this r-dependence of EB can be an experimentally testable prediction. We conclude by discussing some implications to experimental systems featuring resetting dynamics. Y1 - 2021 U6 - https://doi.org/10.1103/PhysRevE.104.024105 SN - 2470-0045 SN - 2470-0053 VL - 104 IS - 2 PB - American Institute of Physics CY - Woodbury, NY ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Thapa, Samudrajit A1 - Mardoukhi, Yousof A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Time averages and their statistical variation for the Ornstein-Uhlenbeck process BT - Role of initial particle distributions and relaxation to stationarity JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - 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. Y1 - 2018 U6 - https://doi.org/10.1103/PhysRevE.98.022134 SN - 2470-0045 SN - 2470-0053 VL - 98 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Fernandez, Amanda Diez A1 - Charchar, Patrick A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf A1 - Finnis, Michael W. T1 - The diffusion of doxorubicin drug molecules in silica nanoslits is non-Gaussian, intermittent and anticorrelated JF - Physical chemistry, chemical physics N2 - In this study we investigate, using all-atom molecular-dynamics computer simulations, the in-plane diffusion of a doxorubicin drug molecule in a thin film of water confined between two silica surfaces. We find that the molecule diffuses along the channel in the manner of a Gaussian diffusion process, but with parameters that vary according to its varying transversal position. Our analysis identifies that four Gaussians, each describing particle motion in a given transversal region, are needed to adequately describe the data. Each of these processes by itself evolves with time at a rate slower than that associated with classical Brownian motion due to a predominance of anticorrelated displacements. Long adsorption events lead to ageing, a property observed when the diffusion is intermittently hindered for periods of time with an average duration which is theoretically infinite. This study presents a simple system in which many interesting features of anomalous diffusion can be explored. It exposes the complexity of diffusion in nanoconfinement and highlights the need to develop new understanding. Y1 - 2020 U6 - https://doi.org/10.1039/d0cp03849k SN - 1463-9076 SN - 1463-9084 VL - 22 IS - 48 SP - 27955 EP - 27965 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Ślęzak, Jakub A1 - Metzler, Ralf A1 - Magdziarz, Marcin T1 - Superstatistical generalised Langevin equation BT - non-Gaussian viscoelastic anomalous diffusion JF - New Journal of Physics N2 - Recent advances in single particle tracking and supercomputing techniques demonstrate the emergence of normal or anomalous, viscoelastic diffusion in conjunction with non-Gaussian distributions in soft, biological, and active matter systems. We here formulate a stochastic model based on a generalised Langevin equation in which non-Gaussian shapes of the probability density function and normal or anomalous diffusion have a common origin, namely a random parametrisation of the stochastic force. We perform a detailed analysis demonstrating how various types of parameter distributions for the memory kernel result in exponential, power law, or power-log law tails of the memory functions. The studied system is also shown to exhibit a further unusual property: the velocity has a Gaussian one point probability density but non-Gaussian joint distributions. This behaviour is reflected in the relaxation from a Gaussian to a non-Gaussian distribution observed for the position variable. We show that our theoretical results are in excellent agreement with stochastic simulations. KW - anomalous diffusion KW - generalised langevin equation KW - superstatistics KW - non-Gaussian diffusion Y1 - 2018 U6 - https://doi.org/10.1088/1367-2630/aaa3d4 SN - 1367-2630 VL - 20 IS - 023026 SP - 1 EP - 25 PB - Deutsche Physikalische Gesellschaft / Institute of Physics CY - Bad Honnef und London ER - TY - JOUR A1 - Krapf, Diego A1 - Metzler, Ralf T1 - Strange interfacial molecular dynamics JF - Physics today Y1 - 2019 U6 - https://doi.org/10.1063/PT.3.4294 SN - 0031-9228 SN - 1945-0699 VL - 72 IS - 9 SP - 48 EP - 54 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Xu, Pengbo A1 - Zhou, Tian A1 - Metzler, Ralf A1 - Deng, Weihua T1 - Stochastic harmonic trapping of a Lévy walk BT - transport and first-passage dynamics under soft resetting strategies JF - New journal of physics : the open-access journal for physics / Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics N2 - We introduce and study a Lé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. KW - diffusion KW - anomalous diffusion KW - stochastic resetting KW - Levy walks Y1 - 2022 U6 - https://doi.org/10.1088/1367-2630/ac5282 SN - 1367-2630 VL - 24 IS - 3 SP - 1 EP - 28 PB - Deutsche Physikalische Gesellschaft CY - Bad Honnef ER - TY - JOUR A1 - Mardoukhi, Yousof A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Spurious ergodicity breaking in normal and fractional Ornstein–Uhlenbeck process JF - New Journal of Physics N2 - 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. KW - Ornstein–Uhlenbeck process KW - stationary stochastic process KW - ensemble and time averaged mean squared displacement Y1 - 2020 U6 - https://doi.org/10.1088/1367-2630/ab950b SN - 1367-2630 VL - 22 PB - IOP CY - London ER - TY - JOUR A1 - Krapf, Diego A1 - Lukat, Nils A1 - Marinari, Enzo A1 - Metzler, Ralf A1 - Oshanin, Gleb A1 - Selhuber-Unkel, Christine A1 - Squarcini, Alessio A1 - Stadler, Lorenz A1 - Weiss, Matthias A1 - Xu, Xinran T1 - Spectral Content of a Single Non-Brownian Trajectory JF - Physical review : X, Expanding access N2 - Time-dependent processes are often analyzed using the power spectral density (PSD) calculated by taking an appropriate Fourier transform of individual trajectories and finding the associated ensemble average. Frequently, the available experimental datasets are too small for such ensemble averages, and hence, it is of a great conceptual and practical importance to understand to which extent relevant information can be gained from S(f, T), the PSD of a single trajectory. Here we focus on the behavior of this random, realization-dependent variable parametrized by frequency f and observation time T, for a broad family of anomalous diffusions-fractional Brownian motion with Hurst index H-and derive exactly its probability density function. We show that S(f, T) is proportional-up to a random numerical factor whose universal distribution we determine-to the ensemble-averaged PSD. For subdiffusion (H < 1/2), we find that S(f, T) similar to A/f(2H+1) with random amplitude A. In sharp contrast, for superdiffusion (H > 1/2) S(f, T) similar to BT2H-1/f(2) with random amplitude B. Remarkably, for H > 1/2 the PSD exhibits the same frequency dependence as Brownian motion, a deceptive property that may lead to false conclusions when interpreting experimental data. Notably, for H > 1/2 the PSD is ageing and is dependent on T. Our predictions for both sub-and superdiffusion are confirmed by experiments in live cells and in agarose hydrogels and by extensive simulations. KW - Biological Physics KW - Interdisciplinary Physics KW - Statistical Physics Y1 - 2019 U6 - https://doi.org/10.1103/PhysRevX.9.011019 SN - 2160-3308 VL - 9 IS - 1 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Sposini, Vittoria A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Single-trajectory spectral analysis of scaled Brownian motion JF - New Journal of Physics N2 - 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. KW - diffusion KW - anomalous diffusion KW - power spectral analysis KW - single trajectory analysis Y1 - 2019 U6 - https://doi.org/10.1088/1367-2630/ab2f52 SN - 1367-2630 VL - 21 PB - Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics CY - Bad Honnef und London ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Search efficiency in the Adam-Delbruck reduction-of-dimensionality scenario versus direct diffusive search JF - New journal of physics : the open-access journal for physics N2 - The time instant-the first-passage time (FPT)-when a diffusive particle (e.g., a ligand such as oxygen or a signalling protein) for the first time reaches an immobile target located on the surface of a bounded three-dimensional domain (e.g., a hemoglobin molecule or the cellular nucleus) is a decisive characteristic time-scale in diverse biophysical and biochemical processes, as well as in intermediate stages of various inter- and intra-cellular signal transduction pathways. Adam and Delbruck put forth the reduction-of-dimensionality concept, according to which a ligand first binds non-specifically to any point of the surface on which the target is placed and then diffuses along this surface until it locates the target. In this work, we analyse the efficiency of such a scenario and confront it with the efficiency of a direct search process, in which the target is approached directly from the bulk and not aided by surface diffusion. We consider two situations: (i) a single ligand is launched from a fixed or a random position and searches for the target, and (ii) the case of 'amplified' signals when N ligands start either from the same point or from random positions, and the search terminates when the fastest of them arrives to the target. For such settings, we go beyond the conventional analyses, which compare only the mean values of the corresponding FPTs. Instead, we calculate the full probability density function of FPTs for both scenarios and study its integral characteristic-the 'survival' probability of a target up to time t. On this basis, we examine how the efficiencies of both scenarios are controlled by a variety of parameters and single out realistic conditions in which the reduction-of-dimensionality scenario outperforms the direct search. KW - first-passage times KW - Adam-Delbruck scenario KW - dimensional reduction KW - bulk KW - and surface diffusion Y1 - 2022 U6 - https://doi.org/10.1088/1367-2630/ac8824 SN - 1367-2630 VL - 24 IS - 8 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Vinod, Deepak A1 - Aghion, Erez A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Scaled geometric Brownian motion features sub- or superexponential ensemble-averaged, but linear time-averaged mean-squared displacements JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Various mathematical Black-Scholes-Merton-like models of option pricing employ the paradigmatic stochastic process of geometric Brownian motion (GBM). The innate property of such models and of real stock-market prices is the roughly exponential growth of prices with time [on average, in crisis-free times]. We here explore the ensemble- and time averages of a multiplicative-noise stochastic process with power-law-like time-dependent volatility, sigma(t) similar to t(alpha), named scaled GBM (SGBM). For SGBM, the mean-squared displacement (MSD) computed for an ensemble of statistically equivalent trajectories can grow faster than exponentially in time, while the time-averaged MSD (TAMSD)-based on a sliding-window averaging along a single trajectory-is always linear at short lag times Delta. The proportionality factor between these the two averages of the time series is Delta/T at short lag times, where T is the trajectory length, similarly to GBM. This discrepancy of the scaling relations and pronounced nonequivalence of the MSD and TAMSD at Delta/T << 1 is a manifestation of weak ergodicity breaking for standard GBM and for SGBM with s (t)-modulation, the main focus of our analysis. The analytical predictions for the MSD and mean TAMSD for SGBM are in quantitative agreement with the results of stochastic computer simulations. Y1 - 2021 U6 - https://doi.org/10.1103/PhysRevE.103.062127 SN - 2470-0045 SN - 2470-0053 VL - 103 IS - 6 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Jeon, Jae-Hyung A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Scaled Brownian motion: a paradoxical process with a time dependent diffusivity for the description of anomalous diffusion JF - Physical chemistry, chemical physics : PCCP N2 - Anomalous diffusion is frequently described by scaled Brownian motion (SBM){,} a Gaussian process with a power-law time dependent diffusion coefficient. Its mean squared displacement is ?x2(t)? [similar{,} equals] 2K(t)t with K(t) [similar{,} equals] t[small alpha]-1 for 0 < [small alpha] < 2. SBM may provide a seemingly adequate description in the case of unbounded diffusion{,} for which its probability density function coincides with that of fractional Brownian motion. Here we show that free SBM is weakly non-ergodic but does not exhibit a significant amplitude scatter of the time averaged mean squared displacement. More severely{,} we demonstrate that under confinement{,} the dynamics encoded by SBM is fundamentally different from both fractional Brownian motion and continuous time random walks. SBM is highly non-stationary and cannot provide a physical description for particles in a thermalised stationary system. Our findings have direct impact on the modelling of single particle tracking experiments{,} in particular{,} under confinement inside cellular compartments or when optical tweezers tracking methods are used. KW - single-particle tracking KW - living cells KW - random-walks KW - subdiffusion KW - dynamics KW - nonergodicity KW - coefficients KW - transport KW - membrane KW - behavior Y1 - 2014 U6 - https://doi.org/10.1039/C4CP02019G VL - 30 IS - 16 SP - 15811 EP - 15817 PB - The Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Wang, Wei A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf A1 - Sokolov, Igor M. T1 - Restoring ergodicity of stochastically reset anomalous-diffusion processes JF - Physical Review Research N2 - How do different reset protocols affect ergodicity of a diffusion process in single-particle-tracking experiments? We here address the problem of resetting of an arbitrary stochastic anomalous-diffusion process (ADP) from the general mathematical points of view and assess ergodicity of such reset ADPs for an arbitrary resetting protocol. The process of stochastic resetting describes the events of the instantaneous restart of a particle’s motion via randomly distributed returns to a preset initial position (or a set of those). The waiting times of such resetting events obey the Poissonian, Gamma, or more generic distributions with specified conditions regarding the existence of moments. Within these general approaches, we derive general analytical results and support them by computer simulations for the behavior of the reset mean-squared displacement (MSD), the new reset increment-MSD (iMSD), and the mean reset time-averaged MSD (TAMSD). For parental nonreset ADPs with the MSD(t)∝ tμ we find a generic behavior and a switch of the short-time growth of the reset iMSD and mean reset TAMSDs from ∝ _μ for subdiffusive to ∝ _1 for superdiffusive reset ADPs. The critical condition for a reset ADP that recovers its ergodicity is found to be more general than that for the nonequilibrium stationary state, where obviously the iMSD and the mean TAMSD are equal. The consideration of the new statistical quantifier, the iMSD—as compared to the standard MSD—restores the ergodicity of an arbitrary reset ADP in all situations when the μth moment of the waiting-time distribution of resetting events is finite. Potential applications of these new resetting results are, inter alia, in the area of biophysical and soft-matter systems. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevResearch.4.013161 SN - 2643-1564 VL - 4 SP - 013161-1 EP - 013161-13 PB - American Physical Society CY - College Park, Maryland, United States ET - 1 ER - TY - JOUR A1 - Singh, Rishu Kumar A1 - Metzler, Ralf A1 - Sandev, Trifce T1 - Resetting dynamics in a confining potential JF - Journal of physics : A, Mathematical and theoretical N2 - 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). KW - diffusion KW - resetting KW - barrier escape KW - first-passage Y1 - 2020 U6 - https://doi.org/10.1088/1751-8121/abc83a SN - 1751-8113 SN - 1751-8121 VL - 53 IS - 50 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Metzler, Ralf A1 - Bauer, Maximilian A1 - Rasmussen, Emil S. A1 - Lomholt, Michael A. T1 - Real sequence effects on the search dynamics of transcription factors on DNA JF - Scientific Reports N2 - Recent experiments show that transcription factors (TFs) indeed use the facilitated diffusion mechanism to locate their target sequences on DNA in living bacteria cells: TFs alternate between sliding motion along DNA and relocation events through the cytoplasm. From simulations and theoretical analysis we study the TF-sliding motion for a large section of the DNA-sequence of a common E. coli strain, based on the two-state TF-model with a fast-sliding search state and a recognition state enabling target detection. For the probability to detect the target before dissociating from DNA the TF-search times self-consistently depend heavily on whether or not an auxiliary operator (an accessible sequence similar to the main operator) is present in the genome section. Importantly, within our model the extent to which the interconversion rates between search and recognition states depend on the underlying nucleotide sequence is varied. A moderate dependence maximises the capability to distinguish between the main operator and similar sequences. Moreover, these auxiliary operators serve as starting points for DNA looping with the main operator, yielding a spectrum of target detection times spanning several orders of magnitude. Auxiliary operators are shown to act as funnels facilitating target detection by TFs. KW - gene regulatory networks KW - biological physics Y1 - 2015 U6 - https://doi.org/10.1038/srep10072 SN - 2045-2322 VL - 5 IS - 10072 PB - Nature Publishing Group CY - London ER - TY - JOUR A1 - Doerries, Timo J. A1 - Chechkin, Aleksei A1 - Schumer, Rina A1 - Metzler, Ralf T1 - Rate equations, spatial moments, and concentration profiles for mobile-immobile models with power-law and mixed waiting time distributions JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We present a framework for systems in which diffusion-advection transport of a tracer substance in a mobile zone is interrupted by trapping in an immobile zone. Our model unifies different model approaches based on distributed-order diffusion equations, exciton diffusion rate models, and random-walk models for multirate mobile-immobile mass transport. We study various forms for the trapping time dynamics and their effects on the tracer mass in the mobile zone. Moreover, we find the associated breakthrough curves, the tracer density at a fixed point in space as a function of time, and the mobile and immobile concentration profiles and the respective moments of the transport. Specifically, we derive explicit forms for the anomalous transport dynamics and an asymptotic power-law decay of the mobile mass for a Mittag-Leffler trapping time distribution. In our analysis we point out that even for exponential trapping time densities, transient anomalous transport is observed. Our results have direct applications in geophysical contexts, but also in biological, soft matter, and solid state systems. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevE.105.014105 SN - 2470-0045 SN - 2470-0053 VL - 105 IS - 1 PB - The American Institute of Physics CY - Woodbury, NY ER - TY - JOUR A1 - Estrada, Ernesto A1 - Delvenne, Jean-Charles A1 - Hatano, Naomichi A1 - Mateos, Jose L. A1 - Metzler, Ralf A1 - Riascos, Alejandro P. A1 - Schaub, Michael T. T1 - Random multi-hopper model BT - super-fast random walks on graphs JF - Journal of Complex Networks N2 - We develop a mathematical model considering a random walker with long-range hops on arbitrary graphs. The random multi-hopper can jump to any node of the graph from an initial position, with a probability that decays as a function of the shortest-path distance between the two nodes in the graph. We consider here two decaying functions in the form of Laplace and Mellin transforms of the shortest-path distances. We prove that when the parameters of these transforms approach zero asymptotically, the hitting time in the multi-hopper approaches the minimum possible value for a normal random walker. We show by computational experiments that the multi-hopper explores a graph with clusters or skewed degree distributions more efficiently than a normal random walker. We provide computational evidences of the advantages of the random multi-hopper model with respect to the normal random walk by studying deterministic, random and real-world networks. Y1 - 2018 U6 - https://doi.org/10.1093/comnet/cnx043 SN - 2051-1310 SN - 2051-1329 VL - 6 IS - 3 SP - 382 EP - 403 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Sposini, Vittoria A1 - Chechkin, Aleksei V. A1 - Seno, Flavio A1 - Pagnini, Gianni A1 - Metzler, Ralf T1 - Random diffusivity from stochastic equations BT - comparison of two models for Brownian yet non-Gaussian diffusion JF - New Journal of Physics N2 - A considerable number of systems have recently been reported in which Brownian yet non-Gaussian dynamics was observed. These are processes characterised by a linear growth in time of the mean squared displacement, yet the probability density function of the particle displacement is distinctly non-Gaussian, and often of exponential(Laplace) shape. This apparently ubiquitous behaviour observed in very different physical systems has been interpreted as resulting from diffusion in inhomogeneous environments and mathematically represented through a variable, stochastic diffusion coefficient. Indeed different models describing a fluctuating diffusivity have been studied. Here we present a new view of the stochastic basis describing time dependent random diffusivities within a broad spectrum of distributions. Concretely, our study is based on the very generic class of the generalised Gamma distribution. Two models for the particle spreading in such random diffusivity settings are studied. The first belongs to the class of generalised grey Brownian motion while the second follows from the idea of diffusing diffusivities. The two processes exhibit significant characteristics which reproduce experimental results from different biological and physical systems. We promote these two physical models for the description of stochastic particle motion in complex environments. Y1 - 2018 U6 - https://doi.org/10.1088/1367-2630/aab696 SN - 1367-2630 SP - 1 EP - 33 PB - Deutsche Physikalische Gesellschaft / Institute of Physics CY - Bad Honnef und London ER - TY - JOUR A1 - Ślęzak, Jakub A1 - Burnecki, Krzysztof A1 - Metzler, Ralf T1 - Random coefficient autoregressive processes describe Brownian yet non-Gaussian diffusion in heterogeneous systems JF - New Journal of Physics N2 - 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. KW - diffusion KW - Langevin equation KW - Brownian yet non-Gaussian diffusion KW - diffusing diffusivity KW - superstatistics KW - autoregressive models KW - time series analysis KW - codifference Y1 - 2019 U6 - https://doi.org/10.1088/1367-2630/ab3366 SN - 1367-2630 VL - 21 PB - Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics CY - Bad Honnef und London ER - TY - JOUR A1 - Schwarzl, Maria A1 - Godec, Aljaž A1 - Metzler, Ralf T1 - Quantifying non-ergodicity of anomalous diffusion with higher order moments JF - Scientific reports N2 - Anomalous diffusion is being discovered in a fast growing number of systems. The exact nature of this anomalous diffusion provides important information on the physical laws governing the studied system. One of the central properties analysed for finite particle motion time series is the intrinsic variability of the apparent diffusivity, typically quantified by the ergodicity breaking parameter EB. Here we demonstrate that frequently EB is insufficient to provide a meaningful measure for the observed variability of the data. Instead, important additional information is provided by the higher order moments entering by the skewness and kurtosis. We analyse these quantities for three popular anomalous diffusion models. In particular, we find that even for the Gaussian fractional Brownian motion a significant skewness in the results of physical measurements occurs and needs to be taken into account. Interestingly, the kurtosis and skewness may also provide sensitive estimates of the anomalous diffusion exponent underlying the data. We also derive a new result for the EB parameter of fractional Brownian motion valid for the whole range of the anomalous diffusion parameter. Our results are important for the analysis of anomalous diffusion but also provide new insights into the theory of anomalous stochastic processes. Y1 - 2017 U6 - https://doi.org/10.1038/s41598-017-03712-x VL - 7 PB - Macmillan Publishers Limited CY - London ER - TY - JOUR A1 - Vojta, Thomas A1 - Skinner, Sarah A1 - Metzler, Ralf T1 - Probability density of the fractional Langevin equation with reflecting walls JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We investigate anomalous diffusion processes governed by the fractional Langevin equation and confined to a finite or semi-infinite interval by reflecting potential barriers. As the random and damping forces in the fractional Langevin equation fulfill the appropriate fluctuation-dissipation relation, the probability density on a finite interval converges for long times towards the expected uniform distribution prescribed by thermal equilibrium. In contrast, on a semi-infinite interval with a reflecting wall at the origin, the probability density shows pronounced deviations from the Gaussian behavior observed for normal diffusion. If the correlations of the random force are persistent (positive), particles accumulate at the reflecting wall while antipersistent (negative) correlations lead to a depletion of particles near the wall. We compare and contrast these results with the strong accumulation and depletion effects recently observed for nonthermal fractional Brownian motion with reflecting walls, and we discuss broader implications. Y1 - 2019 U6 - https://doi.org/10.1103/PhysRevE.100.042142 SN - 2470-0045 SN - 2470-0053 VL - 100 IS - 4 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Krapf, Diego A1 - Marinari, Enzo A1 - Metzler, Ralf A1 - Oshanin, Gleb A1 - Xu, Xinran A1 - Squarcini, Alessio T1 - Power spectral density of a single Brownian trajectory BT - what one can and cannot learn from it JF - New journal of physics : the open-access journal for physics N2 - The power spectral density (PSD) of any time-dependent stochastic processX (t) is ameaningful feature of its spectral content. In its text-book definition, the PSD is the Fourier transform of the covariance function of X-t over an infinitely large observation timeT, that is, it is defined as an ensemble-averaged property taken in the limitT -> infinity. Alegitimate question is what information on the PSD can be reliably obtained from single-trajectory experiments, if one goes beyond the standard definition and analyzes the PSD of a single trajectory recorded for a finite observation timeT. In quest for this answer, for a d-dimensional Brownian motion (BM) we calculate the probability density function of a single-trajectory PSD for arbitrary frequency f, finite observation time T and arbitrary number k of projections of the trajectory on different axes. We show analytically that the scaling exponent for the frequency-dependence of the PSD specific to an ensemble of BM trajectories can be already obtained from a single trajectory, while the numerical amplitude in the relation between the ensemble-averaged and single-trajectory PSDs is afluctuating property which varies from realization to realization. The distribution of this amplitude is calculated exactly and is discussed in detail. Our results are confirmed by numerical simulations and single-particle tracking experiments, with remarkably good agreement. In addition we consider a truncated Wiener representation of BM, and the case of a discrete-time lattice random walk. We highlight some differences in the behavior of a single-trajectory PSD for BM and for the two latter situations. The framework developed herein will allow for meaningful physical analysis of experimental stochastic trajectories. KW - power spectral density KW - single-trajectory analysis KW - probability density function KW - exact results Y1 - 2018 U6 - https://doi.org/10.1088/1367-2630/aaa67c SN - 1367-2630 VL - 20 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Palyulin, Vladimir V. A1 - Ala-Nissila, Tapio A1 - Metzler, Ralf ED - Metzler, Ralf T1 - Polymer translocation: the first two decades and the recent diversification JF - Soft matter N2 - Probably no other field of statistical physics at the borderline of soft matter and biological physics has caused such a flurry of papers as polymer translocation since the 1994 landmark paper by Bezrukov, Vodyanoy, and Parsegian and the study of Kasianowicz in 1996. Experiments, simulations, and theoretical approaches are still contributing novel insights to date, while no universal consensus on the statistical understanding of polymer translocation has been reached. We here collect the published results, in particular, the famous–infamous debate on the scaling exponents governing the translocation process. We put these results into perspective and discuss where the field is going. In particular, we argue that the phenomenon of polymer translocation is non-universal and highly sensitive to the exact specifications of the models and experiments used towards its analysis. KW - solid-state nanopores KW - single-stranded-dna KW - posttranslational protein translocation KW - anomalous diffusion KW - monte-carlo KW - structured polynucleotides KW - dynamics simulation KW - equation approach KW - osmotic-pressure KW - membrane channel Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-76266 SN - 1744-683X VL - 45 IS - 10 SP - 9016 EP - 9037 PB - the Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Eliazar, Iddo A1 - Metzler, Ralf A1 - Reuveni, Shlomi T1 - Poisson-process limit laws yield Gumbel max-min and min-max JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - “A chain is only as strong as its weakest link” says the proverb. But what about a collection of statistically identical chains: How long till all chains fail? The answer to this question is given by the max-min of a matrix whose (i,j)entry is the failure time of link j of chain i: take the minimum of each row, and then the maximum of the rows' minima. The corresponding min-max is obtained by taking the maximum of each column, and then the minimum of the columns' maxima. The min-max applies to the storage of critical data. Indeed, consider multiple backup copies of a set of critical data items, and consider the (i,j) matrix entry to be the time at which item j on copy i is lost; then, the min-max is the time at which the first critical data item is lost. In this paper we address random matrices whose entries are independent and identically distributed random variables. We establish Poisson-process limit laws for the row's minima and for the columns' maxima. Then, we further establish Gumbel limit laws for the max-min and for the min-max. The limit laws hold whenever the entries' distribution has a density, and yield highly applicable approximation tools and design tools for the max-min and min-max of large random matrices. A brief of the results presented herein is given in: Gumbel central limit theorem for max-min and min-max Y1 - 2019 U6 - https://doi.org/10.1103/PhysRevE.100.022129 SN - 2470-0045 SN - 2470-0053 VL - 100 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Li, Yongge A1 - Mei, Ruoxing A1 - Xu, Yong A1 - Kurths, Jürgen A1 - Duan, Jinqiao A1 - Metzler, Ralf T1 - Particle dynamics and transport enhancement in a confined channel with position-dependent diffusivity JF - New Journal of Physics N2 - 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ô 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. KW - diffusion KW - channel KW - space-dependent diffusivity Y1 - 2020 U6 - https://doi.org/10.1088/1367-2630/ab81b9 SN - 1367-2630 VL - 22 PB - Dt. Physikalische Ges. CY - Bad Honnef ER - TY - JOUR A1 - Vinod, Deepak A1 - Cherstvy, Andrey G. A1 - Wang, Wei A1 - Metzler, Ralf A1 - Sokolov, Igor M. T1 - Nonergodicity of reset geometric Brownian motion JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We derive. the ensemble-and time-averaged mean-squared displacements (MSD, TAMSD) for Poisson-reset geometric Brownian motion (GBM), in agreement with simulations. We find MSD and TAMSD saturation for frequent resetting, quantify the spread of TAMSDs via the ergodicity-breaking parameter and compute distributions of prices. General MSD-TAMSD nonequivalence proves reset GBM nonergodic. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevE.105.L012106 SN - 2470-0045 SN - 2470-0053 VL - 105 IS - 1 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Klett, Kolja A1 - Cherstvy, Andrey G. A1 - Shin, Jaeoh A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Non-Gaussian, transiently anomalous, and ergodic self-diffusion of flexible dumbbells in crowded two-dimensional environments BT - coupled translational and rotational motions JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We employ Langevin-dynamics simulations to unveil non-Brownian and non-Gaussian center-of-mass self-diffusion of massive flexible dumbbell-shaped particles in crowded two-dimensional solutions. We study the intradumbbell dynamics of the relative motion of the two constituent elastically coupled disks. Our main focus is on effects of the crowding fraction phi and of the particle structure on the diffusion characteristics. We evaluate the time-averaged mean-squared displacement (TAMSD), the displacement probability-density function (PDF), and the displacement autocorrelation function (ACF) of the dimers. For the TAMSD at highly crowded conditions of dumbbells, e.g., we observe a transition from the short-time ballistic behavior, via an intermediate subdiffusive regime, to long-time Brownian-like spreading dynamics. The crowded system of dimers exhibits two distinct diffusion regimes distinguished by the scaling exponent of the TAMSD, the dependence of the diffusivity on phi, and the features of the displacement-ACF. We attribute these regimes to a crowding-induced transition from viscous to viscoelastic diffusion upon growing phi. We also analyze the relative motion in the dimers, finding that larger phi suppress their vibrations and yield strongly non-Gaussian PDFs of rotational displacements. For the diffusion coefficients D(phi) of translational and rotational motion of the dumbbells an exponential decay with phi for weak and a power-law variation D(phi) proportional to (phi - phi(star))(2.4) for strong crowding is found. A comparison of simulation results with theoretical predictions for D(phi) is discussed and some relevant experimental systems are overviewed. Y1 - 2021 U6 - https://doi.org/10.1103/PhysRevE.104.064603 SN - 2470-0045 SN - 2470-0053 VL - 104 IS - 6 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Thapa, Samudrajit A1 - Wagner, Caroline E. A1 - Metzler, Ralf T1 - Non-Gaussian, non-ergodic, and non-Fickian diffusion of tracers in mucin hydrogels JF - Soft matter N2 - Native mucus is polymer-based soft-matter material of paramount biological importance. How non-Gaussian and non-ergodic is the diffusive spreading of pathogens in mucus? We study the passive, thermally driven motion of micron-sized tracers in hydrogels of mucins, the main polymeric component of mucus. We report the results of the Bayesian analysis for ranking several diffusion models for a set of tracer trajectories [C. E. Wagner et al., Biomacromolecules, 2017, 18, 3654]. The models with "diffusing diffusivity', fractional and standard Brownian motion are used. The likelihood functions and evidences of each model are computed, ranking the significance of each model for individual traces. We find that viscoelastic anomalous diffusion is often most probable, followed by Brownian motion, while the model with a diffusing diffusion coefficient is only realised rarely. Our analysis also clarifies the distribution of time-averaged displacements, correlations of scaling exponents and diffusion coefficients, and the degree of non-Gaussianity of displacements at varying pH levels. Weak ergodicity breaking is also quantified. We conclude that-consistent with the original study-diffusion of tracers in the mucin gels is most non-Gaussian and non-ergodic at low pH that corresponds to the most heterogeneous networks. Using the Bayesian approach with the nested-sampling algorithm, together with the quantitative analysis of multiple statistical measures, we report new insights into possible physical mechanisms of diffusion in mucin gels. Y1 - 2019 U6 - https://doi.org/10.1039/c8sm02096e SN - 1744-683X SN - 1744-6848 VL - 15 IS - 12 SP - 2526 EP - 2551 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Goychuk, Igor A. A1 - Kharchenko, Vasyl O. A1 - Metzler, Ralf T1 - Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion JF - Physical Chemistry Chemical Physics N2 - The discovery of anomalous diffusion of larger biopolymers and submicron tracers such as endogenous granules, organelles, or virus capsids in living cells, attributed to the viscoelastic nature of the cytoplasm, provokes the question whether this complex environment equally impacts the active intracellular transport of submicron cargos by molecular motors such as kinesins: does the passive anomalous diffusion of free cargo always imply its anomalously slow active transport by motors, the mean transport distance along microtubule growing sublinearly rather than linearly in time? Here we analyze this question within the widely used two-state Brownian ratchet model of kinesin motors based on the continuous-state diffusion along microtubules driven by a flashing binding potential, where the cargo particle is elastically attached to the motor. Depending on the cargo size, the loading force, the amplitude of the binding potential, the turnover frequency of the molecular motor enzyme, and the linker stiffness we demonstrate that the motor transport may turn out either normal or anomalous, as indeed measured experimentally. We show how a highly efficient normal active transport mediated by motors may emerge despite the passive anomalous diffusion of the cargo, and study the intricate effects of the elastic linker. Under different, well specified conditions the microtubule-based motor transport becomes anomalously slow and thus significantly less efficient. KW - royal soc chemistry KW - thomas graham house KW - science park KW - milton rd KW - cambridge cb4 0wf KW - cambs KW - england Y1 - 2014 SN - 1463-9076 IS - 16 SP - 16524 EP - 16535 PB - the Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Xu, Pengbo A1 - Zhou, Tian A1 - Metzler, Ralf A1 - Deng, Weihua T1 - Lévy walk dynamics in an external harmonic potential JF - Physical review : E, Statistical, nonlinear, and soft matter physics N2 - Levy walks (LWs) are spatiotemporally coupled random-walk processes describing superdiffusive heat conduction in solids, propagation of light in disordered optical materials, motion of molecular motors in living cells, or motion of animals, humans, robots, and viruses. We here investigate a key feature of LWs-their response to an external harmonic potential. In this generic setting for confined motion we demonstrate that LWs equilibrate exponentially and may assume a bimodal stationary distribution. We also show that the stationary distribution has a horizontal slope next to a reflecting boundary placed at the origin, in contrast to correlated superdiffusive processes. Our results generalize LWs to confining forces and settle some longstanding puzzles around LWs. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevE.101.062127 SN - 2470-0045 SN - 2470-0053 SN - 1550-2376 SN - 1063-651X SN - 1539-3755 VL - 101 IS - 6 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Padash, Amin A1 - Aghion, Erez A1 - Schulz, Alexander A1 - Barkai, Eli A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf A1 - Kantz, Holger T1 - Local equilibrium properties of ultraslow diffusion in the Sinai model JF - New journal of physics N2 - We perform numerical studies of a thermally driven, overdamped particle in a random quenched force field, known as the Sinai model. We compare the unbounded motion on an infinite 1-dimensional domain to the motion in bounded domains with reflecting boundaries and show that the unbounded motion is at every time close to the equilibrium state of a finite system of growing size. This is due to time scale separation: inside wells of the random potential, there is relatively fast equilibration, while the motion across major potential barriers is ultraslow. Quantities studied by us are the time dependent mean squared displacement, the time dependent mean energy of an ensemble of particles, and the time dependent entropy of the probability distribution. Using a very fast numerical algorithm, we can explore times up top 10(17) steps and thereby also study finite-time crossover phenomena. KW - Sinai diffusion KW - clustering KW - local equilibrium Y1 - 2022 U6 - https://doi.org/10.1088/1367-2630/ac7df8 SN - 1367-2630 VL - 24 IS - 7 PB - IOP Publishing CY - Bristol ER - TY - JOUR A1 - Capała, Karol A1 - Padash, Amin A1 - Chechkin, Aleksei V. A1 - Shokri, Babak A1 - Metzler, Ralf A1 - Dybiec, Bartłomiej T1 - Levy noise-driven escape from arctangent potential wells JF - Chaos : an interdisciplinary journal of nonlinear science N2 - The escape from a potential well is an archetypal problem in the study of stochastic dynamical systems, representing real-world situations from chemical reactions to leaving an established home range in movement ecology. Concurrently, Levy noise is a well-established approach to model systems characterized by statistical outliers and diverging higher order moments, ranging from gene expression control to the movement patterns of animals and humans. Here, we study the problem of Levy noise-driven escape from an almost rectangular, arctangent potential well restricted by two absorbing boundaries, mostly under the action of the Cauchy noise. We unveil analogies of the observed transient dynamics to the general properties of stationary states of Levy processes in single-well potentials. The first-escape dynamics is shown to exhibit exponential tails. We examine the dependence of the escape on the shape parameters, steepness, and height of the arctangent potential. Finally, we explore in detail the behavior of the probability densities of the first-escape time and the last-hitting point. Y1 - 2020 U6 - https://doi.org/10.1063/5.0021795 SN - 1054-1500 SN - 1089-7682 VL - 30 IS - 12 PB - American Institute of Physics CY - Woodbury, NY ER - TY - JOUR A1 - Thapa, Samudrajit A1 - Wyłomańska, Agnieszka A1 - Sikora, Grzegorz A1 - Wagner, Caroline E. A1 - Krapf, Diego A1 - Kantz, Holger A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories JF - New Journal of Physics N2 - 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. KW - diffusion KW - anomalous diffusion KW - large-deviation statistic KW - time-averaged mean squared displacement KW - Chebyshev inequality Y1 - 2020 U6 - https://doi.org/10.1088/1367-2630/abd50e SN - 1367-2630 VL - 23 PB - Dt. Physikalische Ges. ; IOP CY - Bad Honnef ; London ER -