TY - JOUR A1 - Cherstvy, Andrey G. A1 - Nagel, Oliver A1 - Beta, Carsten A1 - Metzler, Ralf T1 - Non-Gaussianity, population heterogeneity, and transient superdiffusion in the spreading dynamics of amoeboid cells JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - What is the underlying diffusion process governing the spreading dynamics and search strategies employed by amoeboid cells? Based on the statistical analysis of experimental single-cell tracking data of the two-dimensional motion of the Dictyostelium discoideum amoeboid cells, we quantify their diffusive behaviour based on a number of standard and complementary statistical indicators. We compute the ensemble- and time-averaged mean-squared displacements (MSDs) of the diffusing amoebae cells and observe a pronounced spread of short-time diffusion coefficients and anomalous MSD-scaling exponents for individual cells. The distribution functions of the cell displacements, the long-tailed distribution of instantaneous speeds, and the velocity autocorrelations are also computed. In particular, we observe a systematic superdiffusive short-time behaviour for the ensemble- and time-averaged MSDs of the amoeboid cells. Also, a clear anti-correlation of scaling exponents and generalised diffusivity values for different cells is detected. Most significantly, we demonstrate that the distribution function of the cell displacements has a strongly non-Gaussian shape andusing a rescaled spatio-temporal variablethe cell-displacement data collapse onto a universal master curve. The current analysis of single-cell motions can be implemented for quantifying diffusive behaviours in other living-matter systems, in particular, when effects of active transport, non-Gaussian displacements, and heterogeneity of the population are involved in the dynamics. Y1 - 2018 U6 - https://doi.org/10.1039/c8cp04254c SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 35 SP - 23034 EP - 23054 PB - Royal Society of Chemistry CY - Cambridge 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 - TY - JOUR A1 - Sposini, Vittoria A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - First passage statistics for diffusing diffusivity JF - Journal of physics : A, Mathematical and theoretical N2 - 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. KW - diffusion KW - superstatistics KW - first passage Y1 - 2018 U6 - https://doi.org/10.1088/1751-8121/aaf6ff SN - 1751-8113 SN - 1751-8121 VL - 52 IS - 4 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Dybiec, Bartlomiej A1 - Capala, Karol A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Conservative random walks in confining potentials JF - Journal of physics : A, Mathematical and theoretical N2 - 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. KW - Levy walk KW - conservative random walks KW - Levy flight Y1 - 2018 U6 - https://doi.org/10.1088/1751-8121/aaefc2 SN - 1751-8113 SN - 1751-8121 VL - 52 IS - 1 PB - IOP Publ. Ltd. CY - Bristol 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 - 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 - Eliazar, Iddo A1 - Metzler, Ralf A1 - Reuveni, Shlomi T1 - Gumbel central limit theorem for max-min and min-max JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - The max-min and min-max of matrices arise prevalently in science and engineering. However, in many real-world situations the computation of the max-min and min-max is challenging as matrices are large and full information about their entries is lacking. Here we take a statistical-physics approach and establish limit laws—akin to the central limit theorem—for the max-min and min-max of large random matrices. The limit laws intertwine random-matrix theory and extreme-value theory, couple the matrix dimensions geometrically, and assert that Gumbel statistics emerge irrespective of the matrix entries' distribution. Due to their generality and universality, as well as their practicality, these results are expected to have a host of applications in the physical sciences and beyond. Y1 - 2019 U6 - https://doi.org/10.1103/PhysRevE.100.020104 SN - 2470-0045 SN - 2470-0053 VL - 100 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Javanainen, Matti A1 - Martinez-Seara, Hector A1 - Metzler, Ralf A1 - Vattulainen, Ilpo T1 - Diffusion of Integral Membrane Proteins in Protein-Rich Membranes JF - The journal of physical chemistry letters N2 - The lateral diffusion of embedded proteins along lipid membranes in protein-poor conditions has been successfully described in terms of the Saffman-Delbruck (SD) model, which predicts that the protein diffusion coefficient D is weakly dependent on its radius R as D proportional to ln(1/R). However, instead of being protein-poor, native cell membranes are extremely crowded with proteins. On the basis of extensive molecular simulations, we here demonstrate that protein crowding of the membrane at physiological levels leads to deviations from the SD relation and to the emergence of a stronger Stokes-like dependence D proportional to 1/R. We propose that this 1/R law mainly arises due to geometrical factors: smaller proteins are able to avoid confinement effects much better than their larger counterparts. The results highlight that the lateral dynamics in the crowded setting found in native membranes is radically different from protein-poor conditions and plays a significant role in formation of functional multiprotein complexes. Y1 - 2017 U6 - https://doi.org/10.1021/acs.jpclett.7b01758 SN - 1948-7185 VL - 8 SP - 4308 EP - 4313 PB - American Chemical Society CY - Washington ER - TY - JOUR A1 - Palyulin, Vladimir V. A1 - Mantsevich, Vladimir N. A1 - Klages, Rainer A1 - Metzler, Ralf A1 - Chechkin, Aleksei V. T1 - Comparison of pure and combined search strategies for single and multiple targets JF - The European physical journal : B, Condensed matter and complex systems N2 - 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. Y1 - 2017 U6 - https://doi.org/10.1140/epjb/e2017-80372-4 SN - 1434-6028 SN - 1434-6036 VL - 90 SP - 20 EP - 37 PB - Springer CY - New York ER - TY - JOUR A1 - Caetano, Daniel L. Z. A1 - de Carvalho, Sidney J. A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Critical adsorption of periodic and random polyampholytes onto charged surfaces JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - How different are the properties of critical adsorption of polyampholytes and polyelectrolytes onto charged surfaces? How important are the details of polyampholyte charge distribution on the onset of critical adsorption transition? What are the scaling relations governing the dependence of critical surface charge density on salt concentration in the surrounding solution? Here, we employ Metropolis Monte Carlo simulations and uncover the scaling relations for critical adsorption for quenched periodic and random charge distributions along the polyampholyte chains. We also evaluate and discuss the dependence of the adsorbed layer width on solution salinity and details of the charge distribution. We contrast our findings to the known results for polyelectrolyte adsorption onto oppositely charged surfaces, in particular, their dependence on electrolyte concentration. Y1 - 2017 U6 - https://doi.org/10.1039/c7cp04040g SN - 1463-9076 SN - 1463-9084 VL - 19 SP - 23397 EP - 23413 PB - Royal Society of Chemistry CY - Cambridge 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 - 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 - Palyulin, Vladimir V. A1 - Blackburn, George A1 - Lomholt, Michael A. A1 - Watkins, Nicholas W. A1 - Metzler, Ralf A1 - Klages, Rainer A1 - Chechkin, Aleksei V. T1 - First passage and first hitting times of Levy flights and Levy walks JF - New journal of physics : the open-access journal for physics N2 - For both Lévy flight and Lé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é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é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. KW - Levy flights KW - Levy walks KW - first-passage time KW - first-hitting time Y1 - 2019 U6 - https://doi.org/10.1088/1367-2630/ab41bb SN - 1367-2630 VL - 21 IS - 10 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Teomy, Eial A1 - Metzler, Ralf T1 - Correlations and transport in exclusion processes with general finite memory JF - Journal of statistical mechanics: theory and experiment KW - Brownian motion KW - exclusion processes Y1 - 2019 U6 - https://doi.org/10.1088/1742-5468/ab47fb SN - 1742-5468 VL - 2019 IS - 10 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Hou, Ru A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf A1 - Akimoto, Takuma T1 - Biased continuous-time random walks for ordinary and equilibrium cases BT - facilitation of diffusion, ergodicity breaking and ageing JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - We examine renewal processes with power-law waiting time distributions (WTDs) and non-zero drift via computing analytically and by computer simulations their ensemble and time averaged spreading characteristics. All possible values of the scaling exponent alpha are considered for the WTD psi(t) similar to 1/t(1+alpha). We treat continuous-time random walks (CTRWs) with 0 < alpha < 1 for which the mean waiting time diverges, and investigate the behaviour of the process for both ordinary and equilibrium CTRWs for 1 < alpha < 2 and alpha > 2. We demonstrate that in the presence of a drift CTRWs with alpha < 1 are ageing and non-ergodic in the sense of the non-equivalence of their ensemble and time averaged displacement characteristics in the limit of lag times much shorter than the trajectory length. In the sense of the equivalence of ensemble and time averages, CTRW processes with 1 < alpha < 2 are ergodic for the equilibrium and non-ergodic for the ordinary situation. Lastly, CTRW renewal processes with alpha > 2-both for the equilibrium and ordinary situation-are always ergodic. For the situations 1 < alpha < 2 and alpha > 2 the variance of the diffusion process, however, depends on the initial ensemble. For biased CTRWs with alpha > 1 we also investigate the behaviour of the ergodicity breaking parameter. In addition, we demonstrate that for biased CTRWs the Einstein relation is valid on the level of the ensemble and time averaged displacements, in the entire range of the WTD exponent alpha. Y1 - 2018 U6 - https://doi.org/10.1039/c8cp01863d SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 32 SP - 20827 EP - 20848 PB - Royal Society of Chemistry CY - Cambridge 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 - 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 - GEN A1 - Javanainen, Matti A1 - Martinez-Seara, Hector A1 - Metzler, Ralf A1 - Vattulainen, Ilpo Tapio T1 - Diffusion of Proteins and Lipids in Protein-Rich Membranesa T2 - Biophysical journal Y1 - 2018 U6 - https://doi.org/10.1016/j.bpj.2017.11.3009 SN - 0006-3495 SN - 1542-0086 VL - 114 IS - 3 SP - 551A EP - 551A PB - Cell Press CY - Cambridge 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 - GEN A1 - Gudowska-Nowak, Ewa A1 - Lindenberg, Katja A1 - Metzler, Ralf T1 - Preface: Marian Smoluchowski’s 1916 paper—a century of inspiration T2 - Journal of physics : A, Mathematical and theoretical Y1 - 2017 U6 - https://doi.org/10.1088/1751-8121/aa8529 SN - 1751-8113 SN - 1751-8121 VL - 50 IS - 38 PB - IOP Publ. Ltd. CY - Bristol ER - TY - GEN A1 - Palyulin, Vladimir V A1 - Blackburn, George A1 - Lomholt, Michael A A1 - Watkins, Nicholas W A1 - Metzler, Ralf A1 - Klages, Rainer A1 - Chechkin, Aleksei V. T1 - First passage and first hitting times of Lévy flights and Lévy walks T2 - Postprints der Universität Potsdam Mathematisch-Naturwissenschaftliche Reihe N2 - For both Lévy flight and Lé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é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é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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 785 KW - Lévy flights KW - Lévy walks KW - first-passage time KW - first-hitting time Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-439832 SN - 1866-8372 IS - 785 ER - TY - JOUR A1 - Metzler, Ralf T1 - Superstatistics and non-Gaussian diffusion JF - The European physical journal special topics N2 - Brownian motion and viscoelastic anomalous diffusion in homogeneous environments are intrinsically Gaussian processes. In a growing number of systems, however, non-Gaussian displacement distributions of these processes are being reported. The physical cause of the non-Gaussianity is typically seen in different forms of disorder. These include, for instance, imperfect "ensembles" of tracer particles, the presence of local variations of the tracer mobility in heteroegenous environments, or cases in which the speed or persistence of moving nematodes or cells are distributed. From a theoretical point of view stochastic descriptions based on distributed ("superstatistical") transport coefficients as well as time-dependent generalisations based on stochastic transport parameters with built-in finite correlation time are invoked. After a brief review of the history of Brownian motion and the famed Gaussian displacement distribution, we here provide a brief introduction to the phenomenon of non-Gaussianity and the stochastic modelling in terms of superstatistical and diffusing-diffusivity approaches. KW - Brownian diffusion KW - anomalous diffusion KW - dynamics KW - kinetic-theory KW - models KW - motion KW - nanoparticles KW - nonergodicity KW - statistics KW - subdiffusion Y1 - 2020 U6 - https://doi.org/10.1140/epjst/e2020-900210-x SN - 1951-6355 SN - 1951-6401 VL - 229 IS - 5 SP - 711 EP - 728 PB - Springer CY - Heidelberg ER - TY - GEN A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Population splitting, trapping, and non-ergodicity in heterogeneous diffusion processes N2 - We consider diffusion processes with a spatially varying diffusivity giving rise to anomalous diffusion. Such heterogeneous diffusion processes are analysed for the cases of exponential, power-law, and logarithmic dependencies of the diffusion coefficient on the particle position. Combining analytical approaches with stochastic simulations, we show that the functional form of the space-dependent diffusion coefficient and the initial conditions of the diffusing particles are vital for their statistical and ergodic properties. In all three cases a weak ergodicity breaking between the time and ensemble averaged mean squared displacements is observed. We also demonstrate a population splitting of the time averaged traces into fast and slow diffusers for the case of exponential variation of the diffusivity as well as a particle trapping in the case of the logarithmic diffusivity. Our analysis is complemented by the quantitative study of the space coverage, the diffusive spreading of the probability density, as well as the survival probability. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 236 KW - anomalous diffusion KW - disordered media KW - fractional dynamics KW - infection pathway KW - inhomogeneous-media KW - intracellular-transport KW - langevin equation KW - living cells KW - random-walks KW - single-particle tracking Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-94468 SP - 20220 EP - 20235 ER - TY - JOUR A1 - Metzler, Ralf A1 - Jeon, Jae-Hyung A1 - Cherstvy, Andrey G. A1 - Barkai, Eli T1 - Anomalous diffusion models and their properties BT - non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking JF - physical chemistry, chemical physics : PCCP N2 - Modern microscopic techniques following the stochastic motion of labelled tracer particles have uncovered significant deviations from the laws of Brownian motion in a variety of animate and inanimate systems. Such anomalous diffusion can have different physical origins, which can be identified from careful data analysis. In particular, single particle tracking provides the entire trajectory of the traced particle, which allows one to evaluate different observables to quantify the dynamics of the system under observation. We here provide an extensive overview over different popular anomalous diffusion models and their properties. We pay special attention to their ergodic properties, highlighting the fact that in several of these models the long time averaged mean squared displacement shows a distinct disparity to the regular, ensemble averaged mean squared displacement. In these cases, data obtained from time averages cannot be interpreted by the standard theoretical results for the ensemble averages. Here we therefore provide a comparison of the main properties of the time averaged mean squared displacement and its statistical behaviour in terms of the scatter of the amplitudes between the time averages obtained from different trajectories. We especially demonstrate how anomalous dynamics may be identified for systems, which, on first sight, appear to be Brownian. Moreover, we discuss the ergodicity breaking parameters for the different anomalous stochastic processes and showcase the physical origins for the various behaviours. This Perspective is intended as a guidebook for both experimentalists and theorists working on systems, which exhibit anomalous diffusion. KW - intermittent chaotic systems KW - Fokker-Planck equations KW - time random-walks KW - fluorescence photobleaching recovery KW - fluctuation-dissipation theorem KW - fractional dynamics approach KW - photon-counting statistics KW - weak ergodicity breaking KW - flight search patterns KW - levy flights Y1 - 2014 U6 - https://doi.org/10.1039/c4cp03465a SN - 1463-9076 SN - 1463-9084 VL - 2014 IS - 16 SP - 24128 EP - 24164 ER - TY - GEN A1 - Metzler, Ralf A1 - Jeon, Jae-Hyung A1 - Cherstvy, Andrey G. A1 - Barkai, Eli T1 - Anomalous diffusion models and their properties BT - non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking N2 - Modern microscopic techniques following the stochastic motion of labelled tracer particles have uncovered significant deviations from the laws of Brownian motion in a variety of animate and inanimate systems. Such anomalous diffusion can have different physical origins, which can be identified from careful data analysis. In particular, single particle tracking provides the entire trajectory of the traced particle, which allows one to evaluate different observables to quantify the dynamics of the system under observation. We here provide an extensive overview over different popular anomalous diffusion models and their properties. We pay special attention to their ergodic properties, highlighting the fact that in several of these models the long time averaged mean squared displacement shows a distinct disparity to the regular, ensemble averaged mean squared displacement. In these cases, data obtained from time averages cannot be interpreted by the standard theoretical results for the ensemble averages. Here we therefore provide a comparison of the main properties of the time averaged mean squared displacement and its statistical behaviour in terms of the scatter of the amplitudes between the time averages obtained from different trajectories. We especially demonstrate how anomalous dynamics may be identified for systems, which, on first sight, appear to be Brownian. Moreover, we discuss the ergodicity breaking parameters for the different anomalous stochastic processes and showcase the physical origins for the various behaviours. This Perspective is intended as a guidebook for both experimentalists and theorists working on systems, which exhibit anomalous diffusion. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - paper 174 KW - Fokker-Planck equations KW - flight search patterns KW - fluctuation-dissipation theorem KW - fluorescence photobleaching recovery KW - fractional dynamics approach KW - intermittent chaotic systems KW - levy flights KW - photon-counting statistics KW - time random-walks KW - weak ergodicity breaking Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-74448 SP - 24128 EP - 24164 ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity JF - Soft matter N2 - We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells. KW - anomalous diffusion KW - intracellular-transport KW - adenoassociated virus KW - infection pathway KW - escherichia-coli KW - endosomal escape KW - living cells KW - trafficking KW - cytoplasm KW - models Y1 - 2014 U6 - https://doi.org/10.1039/c3sm52846d SN - 2046-2069 VL - 2014 IS - 10 SP - 1591 EP - 1601 PB - Royal Society of Chemistry ER - TY - GEN A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity N2 - We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - paper 168 KW - adenoassociated virus KW - anomalous diffusion KW - cytoplasm KW - endosomal escape KW - escherichia-coli KW - infection pathway KW - intracellular-transport KW - living cells KW - models KW - trafficking Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-74021 IS - 168 SP - 1591 EP - 1601 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 - GEN A1 - Palyulin, Vladimir V. A1 - Ala-Nissila, Tapio A1 - Metzler, Ralf T1 - Polymer translocation: the first two decades and the recent diversification 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 179 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-76287 SP - 9016 EP - 9037 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 - 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 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf ED - Metzler, Ralf T1 - Non-universal tracer diffusion in crowded media of non-inert obstacles JF - Physical Chemistry Chemical Physics N2 - We study the diffusion of a tracer particle, which moves in continuum space between a lattice of excluded volume, immobile non-inert obstacles. In particular, we analyse how the strength of the tracer–obstacle interactions and the volume occupancy of the crowders alter the diffusive motion of the tracer. From the details of partitioning of the tracer diffusion modes between trapping states when bound to obstacles and bulk diffusion, we examine the degree of localisation of the tracer in the lattice of crowders. We study the properties of the tracer diffusion in terms of the ensemble and time averaged mean squared displacements, the trapping time distributions, the amplitude variation of the time averaged mean squared displacements, and the non-Gaussianity parameter of the diffusing tracer. We conclude that tracer–obstacle adsorption and binding triggers a transient anomalous diffusion. From a very narrow spread of recorded individual time averaged trajectories we exclude continuous type random walk processes as the underlying physical model of the tracer diffusion in our system. For moderate tracer–crowder attraction the motion is found to be fully ergodic, while at stronger attraction strength a transient disparity between ensemble and time averaged mean squared displacements occurs. We also put our results into perspective with findings from experimental single-particle tracking and simulations of the diffusion of tagged tracers in dense crowded suspensions. Our results have implications for the diffusion, transport, and spreading of chemical components in highly crowded environments inside living cells and other structured liquids. KW - fluorescence correlation spectroscopy KW - single-particle tracking KW - anomalous diffusion KW - living cells KW - physiological consequences KW - langevin equation KW - infection pathway KW - excluded volume KW - brownian-motion KW - random-walks Y1 - 2014 SN - 1463-9076 VL - 3 IS - 17 SP - 1847 EP - 1858 PB - The Royal Society of Chemistry CY - Cambridge ER - TY - GEN A1 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Non-universal tracer diffusion in crowded media of non-inert obstacles N2 - We study the diffusion of a tracer particle, which moves in continuum space between a lattice of excluded volume, immobile non-inert obstacles. In particular, we analyse how the strength of the tracer–obstacle interactions and the volume occupancy of the crowders alter the diffusive motion of the tracer. From the details of partitioning of the tracer diffusion modes between trapping states when bound to obstacles and bulk diffusion, we examine the degree of localisation of the tracer in the lattice of crowders. We study the properties of the tracer diffusion in terms of the ensemble and time averaged mean squared displacements, the trapping time distributions, the amplitude variation of the time averaged mean squared displacements, and the non-Gaussianity parameter of the diffusing tracer. We conclude that tracer–obstacle adsorption and binding triggers a transient anomalous diffusion. From a very narrow spread of recorded individual time averaged trajectories we exclude continuous type random walk processes as the underlying physical model of the tracer diffusion in our system. For moderate tracer–crowder attraction the motion is found to be fully ergodic, while at stronger attraction strength a transient disparity between ensemble and time averaged mean squared displacements occurs. We also put our results into perspective with findings from experimental single-particle tracking and simulations of the diffusion of tagged tracers in dense crowded suspensions. Our results have implications for the diffusion, transport, and spreading of chemical components in highly crowded environments inside living cells and other structured liquids. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 186 KW - escence correlation spectroscopy KW - single-particle tracking KW - anomalous diffusion KW - living cells KW - physiological consequences KW - langevin equation KW - infection pathway KW - excluded volume KW - brownian-motion KW - random-walks Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-77128 SP - 1847 EP - 1858 PB - The Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Shin, Jaeoh A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf ED - Metzler, Ralf T1 - Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size JF - Soft Matter N2 - The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping–unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions. KW - gene-regulation kinetics KW - physiological consequences KW - spatial-organization KW - anomalous diffusion KW - folding kinetics KW - living cells KW - dna coiling KW - in-vitro KW - dynamics KW - mixtures Y1 - 2014 SN - 1744-683X SP - 472 EP - 488 PB - The Royal Society of Chemistry CY - Cambridge ER - TY - GEN A1 - Shin, Jaeoh A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size N2 - The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping–unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 185 KW - gene-regulation kinetics KW - physiological consequences KW - spatial-organization KW - anomalous diffusion KW - folding kinetics KW - living cells KW - dna coiling KW - in-vitro KW - dynamics KW - mixtures Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-76961 SP - 472 EP - 488 PB - The Royal Society of Chemistry CY - Cambridge ER - TY - GEN 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 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 189 KW - gene regulatory networks KW - biological physics Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-79411 ER - TY - JOUR A1 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Petrov, Eugene P. A1 - Metzler, Ralf T1 - Interactions of rod-like particles on responsive elastic sheets JF - Soft matter N2 - What are the physical laws of the mutual interactions of objects bound to cell membranes, such as various membrane proteins or elongated virus particles? To rationalise this, we here investigate by extensive computer simulations mutual interactions of rod-like particles adsorbed on the surface of responsive elastic two-dimensional sheets. Specifically, we quantify sheet deformations as a response to adhesion of such filamentous particles. We demonstrate that tip-to-tip contacts of rods are favoured for relatively soft sheets, while side-by-side contacts are preferred for stiffer elastic substrates. These attractive orientation-dependent substrate-mediated interactions between the rod-like particles on responsive sheets can drive their aggregation and self-assembly. The optimal orientation of the membrane-bound rods is established via responding to the elastic energy profiles created around the particles. We unveil the phase diagramme of attractive–repulsive rod–rod interactions in the plane of their separation and mutual orientation. Applications of our results to other systems featuring membrane-associated particles are also discussed. Y1 - 2016 U6 - https://doi.org/10.1039/C6SM01522K SN - 1744-6848 SN - 1744-683X PB - RSC CY - London ER - TY - GEN A1 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Petrov, Eugene P. A1 - Metzler, Ralf T1 - Interactions of rod-like particles on responsive elastic sheets N2 - What are the physical laws of the mutual interactions of objects bound to cell membranes, such as various membrane proteins or elongated virus particles? To rationalise this, we here investigate by extensive computer simulations mutual interactions of rod-like particles adsorbed on the surface of responsive elastic two-dimensional sheets. Specifically, we quantify sheet deformations as a response to adhesion of such filamentous particles. We demonstrate that tip-to-tip contacts of rods are favoured for relatively soft sheets, while side-by-side contacts are preferred for stiffer elastic substrates. These attractive orientation-dependent substrate-mediated interactions between the rod-like particles on responsive sheets can drive their aggregation and self-assembly. The optimal orientation of the membrane-bound rods is established via responding to the elastic energy profiles created around the particles. We unveil the phase diagramme of attractive–repulsive rod–rod interactions in the plane of their separation and mutual orientation. Applications of our results to other systems featuring membrane-associated particles are also discussed. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 256 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-95882 ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Anomalous diffusion in time-fluctuating non-stationary diffusivity landscapes JF - Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies N2 - We investigate the ensemble and time averaged mean squared displacements for particle diffusion in a simple model for disordered media by assuming that the local diffusivity is both fluctuating in time and has a deterministic average growth or decay in time. In this study we compare computer simulations of the stochastic Langevin equation for this random diffusion process with analytical results. We explore the regimes of normal Brownian motion as well as anomalous diffusion in the sub- and superdiffusive regimes. We also consider effects of the inertial term on the particle motion. The investigation of the resulting diffusion is performed for unconfined and confined motion. Y1 - 2016 U6 - https://doi.org/10.1039/C6CP03101C SN - 1463-9084 SN - 1463-9076 VL - 18 SP - 23840 EP - 23852 PB - RSC Publ. CY - Cambridge ER - TY - GEN A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Anomalous diffusion in time-fluctuating non-stationary diffusivity landscapes N2 - We investigate the ensemble and time averaged mean squared displacements for particle diffusion in a simple model for disordered media by assuming that the local diffusivity is both fluctuating in time and has a deterministic average growth or decay in time. In this study we compare computer simulations of the stochastic Langevin equation for this random diffusion process with analytical results. We explore the regimes of normal Brownian motion as well as anomalous diffusion in the sub- and superdiffusive regimes. We also consider effects of the inertial term on the particle motion. The investigation of the resulting diffusion is performed for unconfined and confined motion. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 257 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-95901 SP - 23840 EP - 23852 ER - TY - JOUR A1 - de Carvalho, Sidney J. A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Critical adsorption of polyelectrolytes onto planar and convex highly charged surfaces BT - the nonlinear Poisson–Boltzmann approach JF - New journal of physics : the open-access journal for physics N2 - We study the adsorption–desorption transition of polyelectrolyte chains onto planar, cylindrical and spherical surfaces with arbitrarily high surface charge densities by massive Monte Carlo computer simulations. We examine in detail how the well known scaling relations for the threshold transition—demarcating the adsorbed and desorbed domains of a polyelectrolyte near weakly charged surfaces—are altered for highly charged interfaces. In virtue of high surface potentials and large surface charge densities, the Debye–Hückel approximation is often not feasible and the nonlinear Poisson–Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson–Boltzmann equation is smaller than the Debye–Hückel result, such that the required critical surface charge density for polyelectrolyte adsorption σc increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical–chemical and biophysical systems. KW - polyelectrolyte adsorption KW - electrostatic interactions KW - critical phenomena KW - Debye screening Y1 - 2016 U6 - https://doi.org/10.1088/1367-2630/18/8/083037 SN - 1367-2630 VL - 18 PB - IOP Publ. CY - London ER - TY - GEN A1 - de Carvalho, Sidney J. A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Critical adsorption of polyelectrolytes onto planar and convex highly charged surfaces BT - the nonlinear Poisson–Boltzmann approach N2 - We study the adsorption–desorption transition of polyelectrolyte chains onto planar, cylindrical and spherical surfaces with arbitrarily high surface charge densities by massive Monte Carlo computer simulations. We examine in detail how the well known scaling relations for the threshold transition—demarcating the adsorbed and desorbed domains of a polyelectrolyte near weakly charged surfaces—are altered for highly charged interfaces. In virtue of high surface potentials and large surface charge densities, the Debye–Hückel approximation is often not feasible and the nonlinear Poisson–Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson–Boltzmann equation is smaller than the Debye–Hückel result, such that the required critical surface charge density for polyelectrolyte adsorption σc increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical–chemical and biophysical systems. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 292 KW - polyelectrolyte adsorption KW - electrostatic interactions KW - critical phenomena KW - Debye screening Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-100295 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 - GEN A1 - Mardoukhi, Yousof A1 - Jeon, Jae-Hyung A1 - Metzler, Ralf T1 - Geometry controlled anomalous diffusion in random fractal geometries BT - looking beyond the infinite cluster T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - We investigate the ergodic properties of a random walker performing (anomalous) diffusion on a random fractal geometry. Extensive Monte Carlo simulations of the motion of tracer particles on an ensemble of realisations of percolation clusters are performed for a wide range of percolation densities. Single trajectories of the tracer motion are analysed to quantify the time averaged mean squared displacement (MSD) and to compare this with the ensemble averaged MSD of the particle motion. Other complementary physical observables associated with ergodicity are studied, as well. It turns out that the time averaged MSD of individual realisations exhibits non-vanishing fluctuations even in the limit of very long observation times as the percolation density approaches the critical value. This apparent non-ergodic behaviour concurs with the ergodic behaviour on the ensemble averaged level. We demonstrate how the non-vanishing fluctuations in single particle trajectories are analytically expressed in terms of the fractal dimension and the cluster size distribution of the random geometry, thus being of purely geometrical origin. Moreover, we reveal that the convergence scaling law to ergodicity, which is known to be inversely proportional to the observation time T for ergodic diffusion processes, follows a power-law ∼T−h with h < 1 due to the fractal structure of the accessible space. These results provide useful measures for differentiating the subdiffusion on random fractals from an otherwise closely related process, namely, fractional Brownian motion. Implications of our results on the analysis of single particle tracking experiments are provided. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 980 KW - plasma-membrane KW - mechanisms KW - motion KW - nonergodicity KW - models Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-474864 SN - 1866-8372 IS - 980 SP - 30134 EP - 30147 ER - TY - GEN A1 - Mardoukhi, Yousof A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Spurious ergodicity breaking in normal and fractional Ornstein–Uhlenbeck process T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 981 KW - Ornstein–Uhlenbeck process KW - stationary stochastic process KW - ensemble and time averaged mean squared displacement Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-474875 SN - 1866-8372 IS - 981 ER - TY - GEN A1 - Granado, Felipe Le Vot A1 - Abad, Enrique A1 - Metzler, Ralf A1 - Yuste, Santos B. T1 - Continuous time random walk in a velocity field BT - role of domain growth, Galilei-invariant advection-diffusion, and kinetics of particle mixing T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - We consider the emerging dynamics of a separable continuous time random walk (CTRW) in the case when the random walker is biased by a velocity field in a uniformly growing domain. Concrete examples for such domains include growing biological cells or lipid vesicles, biofilms and tissues, but also macroscopic systems such as expanding aquifers during rainy periods, or the expanding Universe. The CTRW in this study can be subdiffusive, normal diffusive or superdiffusive, including the particular case of a Lévy flight. We first consider the case when the velocity field is absent. In the subdiffusive case, we reveal an interesting time dependence of the kurtosis of the particle probability density function. In particular, for a suitable parameter choice, we find that the propagator, which is fat tailed at short times, may cross over to a Gaussian-like propagator. We subsequently incorporate the effect of the velocity field and derive a bi-fractional diffusion-advection equation encoding the time evolution of the particle distribution. We apply this equation to study the mixing kinetics of two diffusing pulses, whose peaks move towards each other under the action of velocity fields acting in opposite directions. This deterministic motion of the peaks, together with the diffusive spreading of each pulse, tends to increase particle mixing, thereby counteracting the peak separation induced by the domain growth. As a result of this competition, different regimes of mixing arise. In the case of Lévy flights, apart from the non-mixing regime, one has two different mixing regimes in the long-time limit, depending on the exact parameter choice: in one of these regimes, mixing is mainly driven by diffusive spreading, while in the other mixing is controlled by the velocity fields acting on each pulse. Possible implications for encounter–controlled reactions in real systems are discussed. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1005 KW - diffusion KW - expanding medium KW - continuous time random walk Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-479997 SN - 1866-8372 IS - 1005 SP - 28 ER - TY - JOUR A1 - Norregaard, Kamilla A1 - Metzler, Ralf A1 - Ritter, Christine M. A1 - Berg-Sorensen, Kirstine A1 - Oddershede, Lene Broeng T1 - Manipulation and Motion of Organelles and Single Molecules in Living Cells JF - Chemical reviews N2 - The biomolecule is among the most important building blocks of biological systems, and a full understanding of its function forms the scaffold for describing the mechanisms of higher order structures as organelles and cells. Force is a fundamental regulatory mechanism of biomolecular interactions driving many cellular processes. The forces on a molecular scale are exactly in the range that can be manipulated and probed with single molecule force spectroscopy. The natural environment of a biomolecule is inside a living cell, hence, this is the most relevant environment for probing their function. In vivo studies are, however, challenged by the complexity of the cell. In this review, we start with presenting relevant theoretical tools for analyzing single molecule data obtained in intracellular environments followed by a description of state-of-the art visualization techniques. The most commonly used force spectroscopy techniques, namely optical tweezers, magnetic tweezers, and atomic force microscopy, are described in detail, and their strength and limitations related to in vivo experiments are discussed. Finally, recent exciting discoveries within the field of in vivo manipulation and dynamics of single molecule and organelles are reviewed. Y1 - 2017 U6 - https://doi.org/10.1021/acs.chemrev.6b00638 SN - 0009-2665 SN - 1520-6890 VL - 117 IS - 5 SP - 4342 EP - 4375 PB - American Chemical Society CY - Washington ER - TY - GEN A1 - Metzler, Ralf T1 - Gaussianity Fair BT - the Riddle of Anomalous yet Non-Gaussian Diffusion T2 - Biophysical journal Y1 - 2017 U6 - https://doi.org/10.1016/j.bpj.2016.12.019 SN - 0006-3495 SN - 1542-0086 VL - 112 IS - 3 SP - 413 EP - 415 PB - Cell Press CY - Cambridge ER - TY - GEN A1 - Metzler, Ralf T1 - Anomalous Diffusion in Membranes and the Cytoplasm of Biological Cells T2 - Biophysical journal Y1 - 2017 U6 - https://doi.org/10.1016/j.bpj.2016.11.2577 SN - 0006-3495 SN - 1542-0086 VL - 112 IS - 3 SP - 476A EP - 476A PB - Cell Press CY - Cambridge ER - TY - JOUR A1 - Liu, Lin A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Facilitated Diffusion of Transcription Factor Proteins with Anomalous Bulk Diffusion JF - The journal of physical chemistry : B, Condensed matter, materials, surfaces, interfaces & biophysical chemistry N2 - What are the physical laws of the diffusive search of proteins for their specific binding sites on DNA in the presence of the macromolecular crowding in cells? We performed extensive computer simulations to elucidate the protein target search on DNA. The novel feature is the viscoelastic non-Brownian protein bulk diffusion recently observed experimentally. We examine the influence of the protein-DNA binding affinity and the anomalous diffusion exponent on the target search time. In all cases an optimal search time is found. The relative contribution of intermittent three-dimensional bulk diffusion and one-dimensional sliding of proteins along the DNA is quantified. Our results are discussed in the light of recent single molecule tracking experiments, aiming at a better understanding of the influence of anomalous kinetics of proteins on the facilitated diffusion mechanism. Y1 - 2017 U6 - https://doi.org/10.1021/acs.jpcb.6b12413 SN - 1520-6106 VL - 121 SP - 1284 EP - 1289 PB - American Chemical Society CY - Washington ER - TY - JOUR A1 - Godec, Aljaž A1 - Metzler, Ralf T1 - First passage time statistics for two-channel diffusion JF - Journal of physics : A, Mathematical and theoretical N2 - We present rigorous results for the mean first passage time and first passage time statistics for two-channel Markov additive diffusion in a 3-dimensional spherical domain. Inspired by biophysical examples we assume that the particle can only recognise the target in one of the modes, which is shown to effect a non-trivial first passage behaviour. We also address the scenario of intermittent immobilisation. In both cases we prove that despite the perfectly non-recurrent motion of two-channel Markov additive diffusion in 3 dimensions the first passage statistics at long times do not display Poisson-like behaviour if none of the phases has a vanishing diffusion coefficient. This stands in stark contrast to the standard (one-channel) Markov diffusion counterpart. We also discuss the relevance of our results in the context of cellular signalling. KW - first passage time KW - Markov additive processes KW - Fokker-Planck equation KW - random search processes KW - coupled initial boundary value problem KW - cellular signalling KW - asymptotic analysis Y1 - 2017 U6 - https://doi.org/10.1088/1751-8121/aa5204 SN - 1751-8113 SN - 1751-8121 VL - 50 IS - 8 PB - IOP Publ. Ltd. CY - Bristol 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 - Akimoto, Takuma A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Ergodicity, rejuvenation, enhancement, and slow relaxation of diffusion in biased continuous-time random walks JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Bias plays an important role in the enhancement of diffusion in periodic potentials. Using the continuous-time random walk in the presence of a bias, we report on an interesting phenomenon for the enhancement of diffusion by the start of the measurement in a random energy landscape. When the variance of the waiting time diverges, in contrast to the bias-free case, the dynamics with bias becomes superdiffusive. In the superdiffusive regime, we find a distinct initial ensemble dependence of the diffusivity. Moreover, the diffusivity can be increased by the aging time when the initial ensemble is not in equilibrium. We show that the time-averaged variance converges to the corresponding ensemble-averaged variance; i.e., ergodicity is preserved. However, trajectory-to-trajectory fluctuations of the time-averaged variance decay unexpectedly slowly. Our findings provide a rejuvenation phenomenon in the superdiffusive regime, that is, the diffusivity for a nonequilibrium initial ensemble gradually increases to that for an equilibrium ensemble when the start of the measurement is delayed. Y1 - 2018 U6 - https://doi.org/10.1103/PhysRevE.98.022105 SN - 2470-0045 SN - 2470-0053 VL - 98 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Mardoukhi, Yousof A1 - Jeon, Jae-Hyung A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Fluctuations of random walks in critical random environments JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - Percolation networks have been widely used in the description of porous media but are now found to be relevant to understand the motion of particles in cellular membranes or the nucleus of biological cells. Random walks on the infinite cluster at criticality of a percolation network are asymptotically ergodic. On any finite size cluster of the network stationarity is reached at finite times, depending on the cluster's size. Despite of this we here demonstrate by combination of analytical calculations and simulations that at criticality the disorder and cluster size average of the ensemble of clusters leads to a non-vanishing variance of the time averaged mean squared displacement, regardless of the measurement time. Fluctuations of this relevant experimental quantity due to the disorder average of such ensembles are thus persistent and non-negligible. The relevance of our results for single particle tracking analysis in complex and biological systems is discussed. Y1 - 2018 U6 - https://doi.org/10.1039/c8cp03212b SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 31 SP - 20427 EP - 20438 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Kosztolowicz, Tadeusz A1 - Metzler, Ralf A1 - Wąsik, Slawomir A1 - Arabski, Michal T1 - Modelling experimentally measured of ciprofloxacin antibiotic diffusion in Pseudomonas aeruginosa biofilm formed in artificial sputum medium JF - PLoS ONE N2 - We study the experimentally measured ciprofloxacin antibiotic diffusion through a gel-like artificial sputum medium (ASM) mimicking physiological conditions typical for a cystic fibrosis layer, in which regions occupied by Pseudomonas aeruginosa bacteria are present. To quantify the antibiotic diffusion dynamics we employ a phenomenological model using a subdiffusion-absorption equation with a fractional time derivative. This effective equation describes molecular diffusion in a medium structured akin Thompson’s plumpudding model; here the ‘pudding’ background represents the ASM and the ‘plums’ represent the bacterial biofilm. The pudding is a subdiffusion barrier for antibiotic molecules that can affect bacteria found in plums. For the experimental study we use an interferometric method to determine the time evolution of the amount of antibiotic that has diffused through the biofilm. The theoretical model shows that this function is qualitatively different depending on whether or not absorption of the antibiotic in the biofilm occurs. We show that the process can be divided into three successive stages: (1) only antibiotic subdiffusion with constant biofilm parameters, (2) subdiffusion and absorption of antibiotic molecules with variable biofilm transport parameters, (3) subdiffusion and absorption in the medium but the biofilm parameters are constant again. Stage 2 is interpreted as the appearance of an intensive defence build–up of bacteria against the action of the antibiotic, and in the stage 3 it is likely that the bacteria have been inactivated. Times at which stages change are determined from the experimentally obtained temporal evolution of the amount of antibiotic that has diffused through the ASM with bacteria. Our analysis shows good agreement between experimental and theoretical results and is consistent with the biologically expected biofilm response. We show that an experimental method to study the temporal evolution of the amount of a substance that has diffused through a biofilm is useful in studying the processes occurring in a biofilm. We also show that the complicated biological process of antibiotic diffusion in a biofilm can be described by a fractional subdiffusion-absorption equation with subdiffusion and absorption parameters that change over time. KW - Bacterial biofilms KW - Antibiotics KW - Biofilms KW - Cystic fibrosis KW - Absorption KW - Pseudomonas aeruginosa KW - Sputum KW - Biological defense mechanisms Y1 - 2020 U6 - https://doi.org/10.1371/journal.pone.0243003 SN - 1932-6203 VL - 15 PB - PLOS CY - San Francisco, California, US ER - TY - JOUR A1 - Aydiner, Ekrem A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Money distribution in agent-based models with position-exchange dynamics BT - the Pareto paradigm revisited JF - The European physical journal : B, Condensed matter and complex systems N2 - Wealth and income distributions are known to feature country-specific Pareto exponents for their long power-law tails. To propose a rationale for this, we introduce an agent-based dynamic model and use Monte Carlo simulations to unveil the wealth distributions in closed and open economical systems. The standard money-exchange scenario is supplemented with the position-exchange agent dynamics that vitally affects the Pareto law. Specifically, in closed systems with position-exchange dynamics the power law changes to an exponential shape, while for open systems with traps the Pareto law remains valid. KW - Statistical and Nonlinear Physics Y1 - 2019 U6 - https://doi.org/10.1140/epjb/e2019-90674-0 SN - 1434-6028 SN - 1434-6036 VL - 92 IS - 5 PB - Springer CY - New York ER - TY - JOUR A1 - Padash, Amin A1 - Chechkin, Aleksei V. A1 - Dybiec, Bartlomiej A1 - Pavlyukevich, Ilya A1 - Shokri, Babak A1 - Metzler, Ralf T1 - First-passage properties of asymmetric Levy flights JF - Journal of physics : A, Mathematical and theoretical N2 - Lévy flights are paradigmatic generalised random walk processes, in which the independent stationary increments—the 'jump lengths'—are drawn from an -stable jump length distribution with long-tailed, power-law asymptote. As a result, the variance of Lévy flights diverges and the trajectory is characterised by occasional extremely long jumps. Such long jumps significantly decrease the probability to revisit previous points of visitation, rendering Lévy flights efficient search processes in one and two dimensions. To further quantify their precise property as random search strategies we here study the first-passage time properties of Lévy flights in one-dimensional semi-infinite and bounded domains for symmetric and asymmetric jump length distributions. To obtain the full probability density function of first-passage times for these cases we employ two complementary methods. One approach is based on the space-fractional diffusion equation for the probability density function, from which the survival probability is obtained for different values of the stable index and the skewness (asymmetry) parameter . The other approach is based on the stochastic Langevin equation with -stable driving noise. Both methods have their advantages and disadvantages for explicit calculations and numerical evaluation, and the complementary approach involving both methods will be profitable for concrete applications. We also make use of the Skorokhod theorem for processes with independent increments and demonstrate that the numerical results are in good agreement with the analytical expressions for the probability density function of the first-passage times. KW - Levy flights KW - first-passage KW - search dynamics Y1 - 2019 U6 - https://doi.org/10.1088/1751-8121/ab493e SN - 1751-8113 SN - 1751-8121 VL - 52 IS - 45 PB - IOP Publ. Ltd. CY - Bristol 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 - Granado, Felipe Le Vot A1 - Abad, Enrique A1 - Metzler, Ralf A1 - Yuste, Santos B. T1 - Continuous time random walk in a velocity field BT - role of domain growth, Galilei-invariant advection-diffusion, and kinetics of particle mixing JF - New Journal of Physics N2 - We consider the emerging dynamics of a separable continuous time random walk (CTRW) in the case when the random walker is biased by a velocity field in a uniformly growing domain. Concrete examples for such domains include growing biological cells or lipid vesicles, biofilms and tissues, but also macroscopic systems such as expanding aquifers during rainy periods, or the expanding Universe. The CTRW in this study can be subdiffusive, normal diffusive or superdiffusive, including the particular case of a Lévy flight. We first consider the case when the velocity field is absent. In the subdiffusive case, we reveal an interesting time dependence of the kurtosis of the particle probability density function. In particular, for a suitable parameter choice, we find that the propagator, which is fat tailed at short times, may cross over to a Gaussian-like propagator. We subsequently incorporate the effect of the velocity field and derive a bi-fractional diffusion-advection equation encoding the time evolution of the particle distribution. We apply this equation to study the mixing kinetics of two diffusing pulses, whose peaks move towards each other under the action of velocity fields acting in opposite directions. This deterministic motion of the peaks, together with the diffusive spreading of each pulse, tends to increase particle mixing, thereby counteracting the peak separation induced by the domain growth. As a result of this competition, different regimes of mixing arise. In the case of Lévy flights, apart from the non-mixing regime, one has two different mixing regimes in the long-time limit, depending on the exact parameter choice: in one of these regimes, mixing is mainly driven by diffusive spreading, while in the other mixing is controlled by the velocity fields acting on each pulse. Possible implications for encounter–controlled reactions in real systems are discussed. KW - diffusion KW - expanding medium KW - continuous time random walk Y1 - 2020 U6 - https://doi.org/10.1088/1367-2630/ab9ae2 SN - 1367-2630 VL - 22 PB - Dt. Physikalische Ges. CY - Bad Honnef ER - TY - JOUR A1 - Metzler, Ralf T1 - Brownian motion and beyond: first-passage, power spectrum, non-Gaussianity, and anomalous diffusion JF - Journal of statistical mechanics: theory and experiment N2 - Brownian motion is a ubiquitous physical phenomenon across the sciences. After its discovery by Brown and intensive study since the first half of the 20th century, many different aspects of Brownian motion and stochastic processes in general have been addressed in Statistical Physics. In particular, there now exists a very large range of applications of stochastic processes in various disciplines. Here we provide a summary of some of the recent developments in the field of stochastic processes, highlighting both the experimental findings and theoretical frameworks. KW - 15 KW - 4 Y1 - 2019 U6 - https://doi.org/10.1088/1742-5468/ab4988 SN - 1742-5468 VL - 2019 IS - 11 PB - IOP Publ. Ltd. CY - Bristol 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 - 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 - 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 - 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 - GEN A1 - Ślęzak, Jakub A1 - Metzler, Ralf A1 - Magdziarz, Marcin T1 - Codifference can detect ergodicity breaking and non-Gaussianity T2 - Postprints der Universität Potsdam Mathematisch-Naturwissenschaftliche Reihe N2 - We show that the codifference is a useful tool in studying the ergodicity breaking and non-Gaussianity properties of stochastic time series. While the codifference is a measure of dependence that was previously studied mainly in the context of stable processes, we here extend its range of applicability to random-parameter and diffusing-diffusivity models which are important in contemporary physics, biology and financial engineering. We prove that the codifference detects forms of dependence and ergodicity breaking which are not visible from analysing the covariance and correlation functions. We also discuss a related measure of dispersion, which is a nonlinear analogue of the mean squared displacement. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 748 KW - diffusion KW - anomalous diffusion KW - stochastic time series Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-436178 IS - 748 ER - TY - GEN A1 - Sposini, Vittoria A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Single-trajectory spectral analysis of scaled Brownian motion T2 - Postprints der Universität Potsdam Mathematisch-Naturwissenschaftliche Reihe 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 753 KW - diffusion KW - anomalous diffusion KW - power spectral analysis KW - single trajectory analysis Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-436522 SN - 1866-8372 IS - 753 ER - TY - GEN A1 - Guggenberger, Tobias A1 - Pagnini, Gianni A1 - Vojta, Thomas A1 - Metzler, Ralf T1 - Fractional Brownian motion in a finite interval BT - correlations effect depletion or accretion zones of particles near boundaries T2 - Postprints der Universität Potsdam Mathematisch-Naturwissenschaftliche Reihe N2 - Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) 〈X² (t)〉 ⋍ tᵅ with 1 < α < 2, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion (0 < α < 1) this behaviour is reversed and the particle density is depleted close to the boundaries. The MSD in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 755 KW - anomalous diffusion KW - fractional Brownian motion KW - reflecting boundary conditions Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-436665 SN - 1866-8372 IS - 755 ER - TY - GEN A1 - Ślęzak, Jakub A1 - Burnecki, Krzysztof A1 - Metzler, Ralf T1 - Random coefficient autoregressive processes describe Brownian yet non-Gaussian diffusion in heterogeneous systems T2 - Postprints der Universität Potsdam Mathematisch-Naturwissenschaftliche Reihe 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 765 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 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-437923 SN - 1866-8372 IS - 765 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 - Sandev, Trifce A1 - Metzler, Ralf A1 - Chechkin, Aleksei V. T1 - From continuous time random walks to the generalized diffusion equation JF - Fractional calculus and applied analysis : an international journal for theory and applications N2 - We obtain a generalized diffusion equation in modified or Riemann-Liouville form from continuous time random walk theory. The waiting time probability density function and mean squared displacement for different forms of the equation are explicitly calculated. We show examples of generalized diffusion equations in normal or Caputo form that encode the same probability distribution functions as those obtained from the generalized diffusion equation in modified form. The obtained equations are general and many known fractional diffusion equations are included as special cases. KW - continuous time random walk (CTRW) KW - generalized diffusion equation KW - Mittag-Leffler functions KW - anomalous diffusion Y1 - 2018 U6 - https://doi.org/10.1515/fca-2018-0002 SN - 1311-0454 SN - 1314-2224 VL - 21 IS - 1 SP - 10 EP - 28 PB - De Gruyter CY - Berlin ER - TY - JOUR A1 - Wang, Wei A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Thapa, Samudrajit A1 - Seno, Flavio A1 - Liu, Xianbin A1 - Metzler, Ralf T1 - Fractional Brownian motion with random diffusivity BT - emerging residual nonergodicity below the correlation time JF - Journal of physics : A, Mathematical and theoretical N2 - Numerous examples for a priori unexpected non-Gaussian behaviour for normal and anomalous diffusion have recently been reported in single-particle tracking experiments. Here, we address the case of non-Gaussian anomalous diffusion in terms of a random-diffusivity mechanism in the presence of power-law correlated fractional Gaussian noise. We study the ergodic properties of this model via examining the ensemble- and time-averaged mean-squared displacements as well as the ergodicity breaking parameter EB quantifying the trajectory-to-trajectory fluctuations of the latter. For long measurement times, interesting crossover behaviour is found as function of the correlation time tau characterising the diffusivity dynamics. We unveil that at short lag times the EB parameter reaches a universal plateau. The corresponding residual value of EB is shown to depend only on tau and the trajectory length. The EB parameter at long lag times, however, follows the same power-law scaling as for fractional Brownian motion. We also determine a corresponding plateau at short lag times for the discrete representation of fractional Brownian motion, absent in the continuous-time formulation. These analytical predictions are in excellent agreement with results of computer simulations of the underlying stochastic processes. Our findings can help distinguishing and categorising certain nonergodic and non-Gaussian features of particle displacements, as observed in recent single-particle tracking experiments. KW - stochastic processes KW - anomalous diffusion KW - fractional Brownian motion KW - diffusing diffusivity KW - weak ergodicity breaking Y1 - 2020 U6 - https://doi.org/10.1088/1751-8121/aba467 SN - 1751-8113 SN - 1751-8121 VL - 53 IS - 47 PB - IOP Publ. Ltd. CY - Bristol ER - TY - GEN A1 - Kindler, Oliver A1 - Pulkkinen, Otto A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Burst Statistics in an Early Biofilm Quorum Sensing Mode BT - The Role of Spatial Colony-Growth Heterogeneity T2 - Postprints der Universität Potsdam Mathematisch-Naturwissenschaftliche Reihe N2 - Quorum-sensing bacteria in a growing colony of cells send out signalling molecules (so-called “autoinducers”) and themselves sense the autoinducer concentration in their vicinity. Once—due to increased local cell density inside a “cluster” of the growing colony—the concentration of autoinducers exceeds a threshold value, cells in this clusters get “induced” into a communal, multi-cell biofilm-forming mode in a cluster-wide burst event. We analyse quantitatively the influence of spatial disorder, the local heterogeneity of the spatial distribution of cells in the colony, and additional physical parameters such as the autoinducer signal range on the induction dynamics of the cell colony. Spatial inhomogeneity with higher local cell concentrations in clusters leads to earlier but more localised induction events, while homogeneous distributions lead to comparatively delayed but more concerted induction of the cell colony, and, thus, a behaviour close to the mean-field dynamics. We quantify the induction dynamics with quantifiers such as the time series of induction events and burst sizes, the grouping into induction families, and the mean autoinducer concentration levels. Consequences for different scenarios of biofilm growth are discussed, providing possible cues for biofilm control in both health care and biotechnology. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 777 Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-439099 SN - 1866-8372 IS - 777 ER - TY - JOUR A1 - Sarabadani, Jalal A1 - Metzler, Ralf A1 - Ala-Nissila, Tapio T1 - Driven polymer translocation into a channel: Isoflux tension propagation theory and Langevin dynamics simulations JF - Physical Review Research N2 - Isoflux tension propagation (IFTP) theory and Langevin dynamics (LD) simulations are employed to study the dynamics of channel-driven polymer translocation in which a polymer translocates into a narrow channel and the monomers in the channel experience a driving force fc. In the high driving force limit, regardless of the channel width, IFTP theory predicts τ ∝ f βc for the translocation time, where β = −1 is the force scaling exponent. Moreover, LD data show that for a very narrow channel fitting only a single file of monomers, the entropic force due to the subchain inside the channel does not play a significant role in the translocation dynamics and the force exponent β = −1 regardless of the force magnitude. As the channel width increases the number of possible spatial configurations of the subchain inside the channel becomes significant and the resulting entropic force causes the force exponent to drop below unity. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevResearch.4.033003 SN - 2643-1564 VL - 4 SP - 033003-1 EP - 033003-14 PB - American Physical Society CY - College Park, Maryland, USA ET - 3 ER - TY - JOUR A1 - Wang, Wei A1 - Cherstvy, Andrey G. A1 - Liu, Xianbin A1 - Metzler, Ralf T1 - Anomalous diffusion and nonergodicity for heterogeneous diffusion processes with fractional Gaussian noise JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Heterogeneous diffusion processes (HDPs) feature a space-dependent diffusivity of the form D(x) = D-0|x|(alpha). Such processes yield anomalous diffusion and weak ergodicity breaking, the asymptotic disparity between ensemble and time averaged observables, such as the mean-squared displacement. Fractional Brownian motion (FBM) with its long-range correlated yet Gaussian increments gives rise to anomalous and ergodic diffusion. Here, we study a combined model of HDPs and FBM to describe the particle dynamics in complex systems with position-dependent diffusivity driven by fractional Gaussian noise. This type of motion is, inter alia, relevant for tracer-particle diffusion in biological cells or heterogeneous complex fluids. We show that the long-time scaling behavior predicted theoretically and by simulations for the ensemble-and time-averaged mean-squared displacements couple the scaling exponents alpha of HDPs and the Hurst exponent H of FBM in a characteristic way. Our analysis of the simulated data in terms of the rescaled variable y similar to |x|(1/(2/(2-alpha)))/t(H) coupling particle position x and time t yields a simple, Gaussian probability density function (PDF), PHDP-FBM(y) = e(-y2)/root pi. Its universal shape agrees well with theoretical predictions for both uni- and bimodal PDF distributions. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevE.102.012146 SN - 2470-0045 SN - 2470-0053 SN - 1063-651X SN - 1539-3755 SN - 2470-0061 VL - 102 IS - 1 SP - 012146-1 EP - 012146-16 PB - American Physical Society CY - College Park 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 - 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 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 SP - 033055-1 EP - 033055-16 PB - American Physical Society CY - College Park, MD ER - TY - JOUR A1 - Seckler, Henrik A1 - Metzler, Ralf T1 - Bayesian deep learning for error estimation in the analysis of anomalous diffusion JF - Nature Communnications N2 - Modern single-particle-tracking techniques produce extensive time-series of diffusive motion in a wide variety of systems, from single-molecule motion in living-cells to movement ecology. The quest is to decipher the physical mechanisms encoded in the data and thus to better understand the probed systems. We here augment recently proposed machine-learning techniques for decoding anomalous-diffusion data to include an uncertainty estimate in addition to the predicted output. To avoid the Black-Box-Problem a Bayesian-Deep-Learning technique named Stochastic-Weight-Averaging-Gaussian is used to train models for both the classification of the diffusionmodel and the regression of the anomalous diffusion exponent of single-particle-trajectories. Evaluating their performance, we find that these models can achieve a wellcalibrated error estimate while maintaining high prediction accuracies. In the analysis of the output uncertainty predictions we relate these to properties of the underlying diffusion models, thus providing insights into the learning process of the machine and the relevance of the output. KW - random-walk KW - models Y1 - 2022 U6 - https://doi.org/10.1038/s41467-022-34305-6 SN - 2041-1723 VL - 13 PB - Nature Publishing Group UK CY - London ER - TY - JOUR A1 - Eliazar, Iddo A1 - Metzler, Ralf T1 - Anomalous statistics of random relaxations in random environments JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We comprehensively analyze the emergence of anomalous statistics in the context of the random relaxation ( RARE) model [Eliazar and Metzler, J. Chem. Phys. 137, 234106 ( 2012)], a recently introduced versatile model of random relaxations in random environments. The RARE model considers excitations scattered randomly across a metric space around a reaction center. The excitations react randomly with the center, the reaction rates depending on the excitations' distances from this center. Relaxation occurs upon the first reaction between an excitation and the center. Addressing both the relaxation time and the relaxation range, we explore when these random variables display anomalous statistics, namely, heavy tails at zero and at infinity that manifest, respectively, exceptionally high occurrence probabilities of very small and very large outliers. A cohesive set of closed-form analytic results is established, determining precisely when such anomalous statistics emerge. Y1 - 2013 U6 - https://doi.org/10.1103/PhysRevE.87.022141 SN - 1539-3755 SN - 1550-2376 VL - 87 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Kursawe, Jochen A1 - Schulz, Johannes H. P. A1 - Metzler, Ralf T1 - Transient aging in fractional brownian and langevin-equation motion JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Stochastic processes driven by stationary fractional Gaussian noise, that is, fractional Brownian motion and fractional Langevin-equation motion, are usually considered to be ergodic in the sense that, after an algebraic relaxation, time and ensemble averages of physical observables coincide. Recently it was demonstrated that fractional Brownian motion and fractional Langevin-equation motion under external confinement are transiently nonergodic-time and ensemble averages behave differently-from the moment when the particle starts to sense the confinement. Here we show that these processes also exhibit transient aging, that is, physical observables such as the time-averaged mean-squared displacement depend on the time lag between the initiation of the system at time t = 0 and the start of the measurement at the aging time t(a). In particular, it turns out that for fractional Langevin-equation motion the aging dependence on ta is different between the cases of free and confined motion. We obtain explicit analytical expressions for the aged moments of the particle position as well as the time-averaged mean-squared displacement and present a numerical analysis of this transient aging phenomenon. Y1 - 2013 U6 - https://doi.org/10.1103/PhysRevE.88.062124 SN - 1539-3755 SN - 1550-2376 VL - 88 IS - 6 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Bauer, Maximilian A1 - Metzler, Ralf T1 - In vivo facilitated diffusion model JF - PLoS one N2 - Under dilute in vitro conditions transcription factors rapidly locate their target sequence on DNA by using the facilitated diffusion mechanism. However, whether this strategy of alternating between three-dimensional bulk diffusion and one-dimensional sliding along the DNA contour is still beneficial in the crowded interior of cells is highly disputed. Here we use a simple model for the bacterial genome inside the cell and present a semi-analytical model for the in vivo target search of transcription factors within the facilitated diffusion framework. Without having to resort to extensive simulations we determine the mean search time of a lac repressor in a living E. coli cell by including parameters deduced from experimental measurements. The results agree very well with experimental findings, and thus the facilitated diffusion picture emerges as a quantitative approach to gene regulation in living bacteria cells. Furthermore we see that the search time is not very sensitive to the parameters characterizing the DNA configuration and that the cell seems to operate very close to optimal conditions for target localization. Local searches as implied by the colocalization mechanism are only found to mildly accelerate the mean search time within our model. Y1 - 2013 U6 - https://doi.org/10.1371/journal.pone.0053956 SN - 1932-6203 VL - 8 IS - 1 PB - PLoS CY - San Fransisco ER - TY - JOUR A1 - Godec, Aljaz A1 - Metzler, Ralf T1 - Finite-Time effects and ultraweak ergodicity breaking in superdiffusive dynamics JF - Physical review letters N2 - We study the ergodic properties of superdiffusive, spatiotemporally coupled Levy walk processes. For trajectories of finite duration, we reveal a distinct scatter of the scaling exponents of the time averaged mean squared displacement (delta x(2)) over bar around the ensemble value 3 - alpha (1 < alpha < 2) ranging from ballistic motion to subdiffusion, in strong contrast to the behavior of subdiffusive processes. In addition we find a significant dependence of the average of (delta x(2)) over bar over an ensemble of trajectories as a function of the finite measurement time. This so-called finite-time amplitude depression and the scatter of the scaling exponent is vital in the quantitative evaluation of superdiffusive processes. Comparing the long time average of the second moment with the ensemble mean squared displacement, these only differ by a constant factor, an ultraweak ergodicity breaking. Y1 - 2013 U6 - https://doi.org/10.1103/PhysRevLett.110.020603 SN - 0031-9007 VL - 110 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Schulz, Johannes H. P. A1 - Barkai, Eli A1 - Metzler, Ralf T1 - Aging effects and population splitting in single-particle trajectoryaverages JF - Physical review letters N2 - We study time averages of single particle trajectories in scale-free anomalous diffusion processes, in which the measurement starts at some time t(a) > 0 after initiation of the process at t = 0. Using aging renewal theory, we show that for such nonstationary processes a large class of observables are affected by a unique aging function, which is independent of boundary conditions or the external forces. Moreover, we discuss the implications of aging induced population splitting: with growing age ta of the process, an increasing fraction of particles remains motionless in a measurement of fixed duration. Consequences for single biomolecule tracking in live cells are discussed. Y1 - 2013 U6 - https://doi.org/10.1103/PhysRevLett.110.020602 SN - 0031-9007 SN - 1079-7114 VL - 110 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Population splitting, trapping, and non-ergodicity in heterogeneous diffusion processes JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - We consider diffusion processes with a spatially varying diffusivity giving rise to anomalous diffusion. Such heterogeneous diffusion processes are analysed for the cases of exponential, power-law, and logarithmic dependencies of the diffusion coefficient on the particle position. Combining analytical approaches with stochastic simulations, we show that the functional form of the space-dependent diffusion coefficient and the initial conditions of the diffusing particles are vital for their statistical and ergodic properties. In all three cases a weak ergodicity breaking between the time and ensemble averaged mean squared displacements is observed. We also demonstrate a population splitting of the time averaged traces into fast and slow diffusers for the case of exponential variation of the diffusivity as well as a particle trapping in the case of the logarithmic diffusivity. Our analysis is complemented by the quantitative study of the space coverage, the diffusive spreading of the probability density, as well as the survival probability. Y1 - 2013 U6 - https://doi.org/10.1039/c3cp53056f SN - 1463-9076 SN - 1463-9084 VL - 15 IS - 46 SP - 20220 EP - 20235 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Vahabi, Mahsa A1 - Schulz, Johannes H. P. A1 - Shokri, Babak A1 - Metzler, Ralf T1 - Area coverage of radial Levy flights with periodic boundary conditions JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We consider the area coverage of radial Levy flights in a finite square area with periodic boundary conditions. From simulations we show how the fractal path dimension d(f) and thus the degree of area coverage depends on the number of steps of the trajectory, the size of the area, and the resolution of the applied box counting algorithm. For sufficiently long trajectories and not too high resolution, the fractal dimension returned by the box counting method equals two, and in that sense the Levy flight fully covers the area. Otherwise, the determined fractal dimension equals the stable index of the distribution of jump lengths of the Levy flight. We provide mathematical expressions for the turnover between these two scaling regimes. As complementary methods to analyze confined Levy flights we investigate fractional order moments of the position for which we also provide scaling arguments. Finally, we study the time evolution of the probability density function and the first passage time density of Levy flights in a square area. Our findings are of interest for a general understanding of Levy flights as well as for the analysis of recorded trajectories of animals searching for food or for human motion patterns. Y1 - 2013 U6 - https://doi.org/10.1103/PhysRevE.87.042136 SN - 1539-3755 VL - 87 IS - 4 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Jeon, Jae-Hyung A1 - Leijnse, Natascha A1 - Oddershede, Lene B. A1 - Metzler, Ralf T1 - Anomalous diffusion and power-law relaxation of the time averaged mean squared displacement in worm-like micellar solutions JF - New journal of physics : the open-access journal for physics N2 - We report the results of single tracer particle tracking by optical tweezers and video microscopy in micellar solutions. From careful analysis in terms of different stochastic models, we show that the polystyrene tracer beads of size 0.52-2.5 mu m after short-time normal diffusion turn over to perform anomalous diffusion of the form < r(2)(t)> similar or equal to t(alpha) with alpha approximate to 0.3. This free anomalous diffusion is ergodic and consistent with a description in terms of the generalized Langevin equation with a power-law memory kernel. With optical tweezers tracking, we unveil a power-law relaxation over several decades in time to the thermal plateau value under the confinement of the harmonic tweezer potential, as predicted previously (Phys. Rev. E 85 021147 (2012)). After the subdiffusive motion in the millisecond range, the motion becomes faster and turns either back to normal Brownian diffusion or to even faster superdiffusion, depending on the size of the tracer beads. Y1 - 2013 U6 - https://doi.org/10.1088/1367-2630/15/4/045011 SN - 1367-2630 VL - 15 IS - 4 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Maity, Alok Kumar A1 - Bandyopadhyay, Arnab A1 - Chattopadhyay, Sudip A1 - Chaudhuri, Jyotipratim Ray A1 - Metzler, Ralf A1 - Chaudhury, Pinaki A1 - Banik, Suman K. T1 - Quantification of noise in bifunctionality-induced post-translational modification JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We present a generic analytical scheme for the quantification of fluctuations due to bifunctionality-induced signal transduction within the members of a bacterial two-component system. The proposed model takes into account post-translational modifications in terms of elementary phosphotransfer kinetics. Sources of fluctuations due to autophosphorylation, kinase, and phosphatase activity of the sensor kinase have been considered in the model via Langevin equations, which are then solved within the framework of linear noise approximation. The resultant analytical expression of phosphorylated response regulators are then used to quantify the noise profile of biologically motivated single and branched pathways. Enhancement and reduction of noise in terms of extra phosphate outflux and influx, respectively, have been analyzed for the branched system. Furthermore, the role of fluctuations of the network output in the regulation of a promoter with random activation-deactivation dynamics has been analyzed. Y1 - 2013 U6 - https://doi.org/10.1103/PhysRevE.88.032716 SN - 1539-3755 VL - 88 IS - 3 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Jeon, Jae-Hyung A1 - Barkai, Eli A1 - Metzler, Ralf T1 - Noisy continuous time random walks JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - Experimental studies of the diffusion of biomolecules within biological cells are routinely confronted with multiple sources of stochasticity, whose identification renders the detailed data analysis of single molecule trajectories quite intricate. Here, we consider subdiffusive continuous time random walks that represent a seminal model for the anomalous diffusion of tracer particles in complex environments. This motion is characterized by multiple trapping events with infinite mean sojourn time. In real physical situations, however, instead of the full immobilization predicted by the continuous time random walk model, the motion of the tracer particle shows additional jiggling, for instance, due to thermal agitation of the environment. We here present and analyze in detail an extension of the continuous time random walk model. Superimposing the multiple trapping behavior with additive Gaussian noise of variable strength, we demonstrate that the resulting process exhibits a rich variety of apparent dynamic regimes. In particular, such noisy continuous time random walks may appear ergodic, while the bare continuous time random walk exhibits weak ergodicity breaking. Detailed knowledge of this behavior will be useful for the truthful physical analysis of experimentally observed subdiffusion. Y1 - 2013 U6 - https://doi.org/10.1063/1.4816635 SN - 0021-9606 SN - 1089-7690 VL - 139 IS - 12 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Population splitting, trapping, and non-ergodicity in heterogeneous diffusion processes N2 - We consider diffusion processes with a spatially varying diffusivity giving rise to anomalous diffusion. Such heterogeneous diffusion processes are analysed for the cases of exponential, power-law, and logarithmic dependencies of the diffusion coefficient on the particle position. Combining analytical approaches with stochastic simulations, we show that the functional form of the space-dependent diffusion coefficient and the initial conditions of the diffusing particles are vital for their statistical and ergodic properties. In all three cases a weak ergodicity breaking between the time and ensemble averaged mean squared displacements is observed. We also demonstrate a population splitting of the time averaged traces into fast and slow diffusers for the case of exponential variation of the diffusivity as well as a particle trapping in the case of the logarithmic diffusivity. Our analysis is complemented by the quantitative study of the space coverage, the diffusive spreading of the probability density, as well as the survival probability. Y1 - 2008 UR - http://pubs.rsc.org/en/content/articlehtml/2013/cp/c3cp53056f U6 - https://doi.org/10.1039/C3CP53056F ER - TY - JOUR A1 - Schulz, Johannes H. P. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Correlated continuous time random walks - combining scale-invariance with long-range memory for spatial and temporal dynamics JF - Journal of physics : A, Mathematical and theoretical N2 - Standard continuous time random walk (CTRW) models are renewal processes in the sense that at each jump a new, independent pair of jump length and waiting time are chosen. Globally, anomalous diffusion emerges through scale-free forms of the jump length and/or waiting time distributions by virtue of the generalized central limit theorem. Here we present a modified version of recently proposed correlated CTRW processes, where we incorporate a power-law correlated noise on the level of both jump length and waiting time dynamics. We obtain a very general stochastic model, that encompasses key features of several paradigmatic models of anomalous diffusion: discontinuous, scale-free displacements as in Levy flights, scale-free waiting times as in subdiffusive CTRWs, and the long-range temporal correlations of fractional Brownian motion (FBM). We derive the exact solutions for the single-time probability density functions and extract the scaling behaviours. Interestingly, we find that different combinations of the model parameters lead to indistinguishable shapes of the emerging probability density functions and identical scaling laws. Our model will be useful for describing recent experimental single particle tracking data that feature a combination of CTRW and FBM properties. Y1 - 2013 U6 - https://doi.org/10.1088/1751-8113/46/47/475001 SN - 1751-8113 SN - 1751-8121 VL - 46 IS - 47 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Talukder, Srijeeta A1 - Sen, Shrabani A1 - Metzler, Ralf A1 - Banik, Suman K. A1 - Chaudhury, Pinaki T1 - Stochastic optimization-based study of dimerization kinetics JF - JOURNAL OF CHEMICAL SCIENCES N2 - We investigate the potential of numerical algorithms to decipher the kinetic parameters involved in multi-step chemical reactions. To this end, we study dimerization kinetics of protein as a model system. We follow the dimerization kinetics using a stochastic simulation algorithm and combine it with three different optimization techniques (genetic algorithm, simulated annealing and parallel tempering) to obtain the rate constants involved in each reaction step. We find good convergence of the numerical scheme to the rate constants of the process. We also perform a sensitivity test on the reaction kinetic parameters to see the relative effects of the parameters for the associated profile of the monomer/dimer distribution. KW - Stochastic optimization KW - dimerization kinetics KW - sensitivity analysis KW - stochastic simulation algorithm KW - probability distribution function Y1 - 2013 SN - 0974-3626 SN - 0973-7103 VL - 125 IS - 6 SP - 1619 EP - 1627 PB - INDIAN ACAD SCIENCES CY - BANGALORE ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Anomalous diffusion and ergodicity breaking in heterogeneous diffusion processes JF - New journal of physics : the open-access journal for physics N2 - We demonstrate the non-ergodicity of a simple Markovian stochastic process with space-dependent diffusion coefficient D(x). For power-law forms D(x) similar or equal to vertical bar x vertical bar(alpha), this process yields anomalous diffusion of the form < x(2)(t)> similar or equal to t(2/(2-alpha)). Interestingly, in both the sub- and superdiffusive regimes we observe weak ergodicity breaking: the scaling of the time-averaged mean-squared displacement <(delta(2)(Delta))over bar> remains linear in the lag time Delta and thus differs from the corresponding ensemble average < x(2)(t)>. We analyse the non-ergodic behaviour of this process in terms of the time-averaged mean- squared displacement (delta(2)) over bar and its random features, i.e. the statistical distribution of (delta(2)) over bar and the ergodicity breaking parameters. The heterogeneous diffusion model represents an alternative approach to non- ergodic, anomalous diffusion that might be particularly relevant for diffusion in heterogeneous media. Y1 - 2013 U6 - https://doi.org/10.1088/1367-2630/15/8/083039 SN - 1367-2630 VL - 15 IS - 15 PB - IOP Publ. Ltd. CY - Bristol ER - TY - GEN A1 - Barkai, Eli A1 - Garini, Yuval A1 - Metzler, Ralf T1 - Electrostatic effects in living cells Reply T2 - PHYSICS TODAY Y1 - 2013 SN - 0031-9228 SN - 1945-0699 VL - 66 IS - 7 SP - 11 EP - 11 PB - AMER INST PHYSICS CY - MELVILLE ER - TY - JOUR A1 - Lomholt, Michael A. A1 - Lizana, Ludvig A1 - Metzler, Ralf A1 - Ambjoernsson, Tobias T1 - Microscopic origin of the logarithmic time evolution of aging processes in complex systems JF - Physical review letters N2 - There exists compelling experimental evidence in numerous systems for logarithmically slow time evolution, yet its full theoretical understanding remains elusive. We here introduce and study a generic transition process in complex systems, based on nonrenewal, aging waiting times. Each state n of the system follows a local clock initiated at t = 0. The random time tau between clock ticks follows the waiting time density psi (tau). Transitions between states occur only at local clock ticks and are hence triggered by the local forward waiting time, rather than by psi (tau). For power-law forms psi (tau) similar or equal to tau(-1-alpha) (0 < alpha < 1) we obtain a logarithmic time evolution of the state number < n(t)> similar or equal to log(t/t(0)), while for alpha > 2 the process becomes normal in the sense that < n(t)> similar or equal to t. In the intermediate range 1 < alpha < 2 we find the power-law growth < n(t)> similar or equal to t(alpha-1). Our model provides a universal description for transition dynamics between aging and nonaging states. Y1 - 2013 U6 - https://doi.org/10.1103/PhysRevLett.110.208301 SN - 0031-9007 VL - 110 IS - 20 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Pulkkinen, Otto A1 - Metzler, Ralf T1 - Distance matters the impact of gene proximity in bacterial gene regulation JF - Physical review letters N2 - Following recent discoveries of colocalization of downstream-regulating genes in living cells, the impact of the spatial distance between such genes on the kinetics of gene product formation is increasingly recognized. We here show from analytical and numerical analysis that the distance between a transcription factor (TF) gene and its target gene drastically affects the speed and reliability of transcriptional regulation in bacterial cells. For an explicit model system, we develop a general theory for the interactions between a TF and a transcription unit. The observed variations in regulation efficiency are linked to the magnitude of the variation of the TF concentration peaks as a function of the binding site distance from the signal source. Our results support the role of rapid binding site search for gene colocalization and emphasize the role of local concentration differences. Y1 - 2013 U6 - https://doi.org/10.1103/PhysRevLett.110.198101 SN - 0031-9007 VL - 110 IS - 19 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Jeon, Jae-Hyung A1 - Metzler, Ralf T1 - Inequivalence of time and ensemble averages in ergodic systems: exponential versus power-law relaxation in confinement JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Single-particle tracking has become a standard tool for the investigation of diffusive properties, especially in small systems such as biological cells. Usually the resulting time series are analyzed in terms of time averages over individual trajectories. Here we study confined normal as well as anomalous diffusion, modeled by fractional Brownian motion and the fractional Langevin equation, and show that even for such ergodic systems time-averaged quantities behave differently from their ensemble-averaged counterparts, irrespective of how long the measurement time becomes. Knowledge of the exact behavior of time averages is therefore fundamental for the proper physical interpretation of measured time series, in particular, for extraction of the relaxation time scale from data. Y1 - 2012 U6 - https://doi.org/10.1103/PhysRevE.85.021147 SN - 1539-3755 VL - 85 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Palyulin, Vladimir V. A1 - Metzler, Ralf T1 - How a finite potential barrier decreases the mean first-passage time JF - Journal of statistical mechanics: theory and experiment N2 - We consider the mean first-passage time of a random walker moving in a potential landscape on a finite interval, the starting and end points being at different potentials. From analytical calculations and Monte Carlo simulations we demonstrate that the mean first-passage time for a piecewise linear curve between these two points is minimized by the introduction of a potential barrier. Due to thermal fluctuations, this barrier may be crossed. It turns out that the corresponding expense for this activation is less severe than the gain from an increased slope towards the end point. In particular, the resulting mean first-passage time is shorter than for a linear potential drop between the two points. KW - diffusion Y1 - 2012 U6 - https://doi.org/10.1088/1742-5468/2012/03/L03001 SN - 1742-5468 IS - 1 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Godec, Aljaz A1 - Bauer, Maximilian A1 - Metzler, Ralf T1 - Collective dynamics effect transient subdiffusion of inert tracers in flexible gel networks JF - New journal of physics : the open-access journal for physics N2 - Based on extensive Brownian dynamics simulations we study the thermal motion of a tracer bead in a cross-linked, flexible gel in the limit when the tracer particle size is comparable to or even larger than the equilibrium mesh size of the gel. The analysis of long individual trajectories of the tracer demonstrates the existence of pronounced transient anomalous diffusion. From the time averaged mean squared displacement and the time averaged van Hove correlation functions we elucidate the many-body origin of the non-Brownian tracer bead dynamics. Our results shed new light onto the ongoing debate over the physical origin of steric tracer interactions with structured environments. KW - anomalous diffusion KW - gel network KW - van Hove correlation Y1 - 2014 U6 - https://doi.org/10.1088/1367-2630/16/9/092002 SN - 1367-2630 VL - 16 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Bauer, Maximilian A1 - Godec, Aljaz A1 - Metzler, Ralf T1 - Diffusion of finite-size particles in two-dimensional channels with random wall configurations JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - Diffusion of chemicals or tracer molecules through complex systems containing irregularly shaped channels is important in many applications. Most theoretical studies based on the famed Fick-Jacobs equation focus on the idealised case of infinitely small particles and reflecting boundaries. In this study we use numerical simulations to consider the transport of finite-size particles through asymmetrical two-dimensional channels. Additionally, we examine transient binding of the molecules to the channel walls by applying sticky boundary conditions. We consider an ensemble of particles diffusing in independent channels, which are characterised by common structural parameters. We compare our results for the long-time effective diffusion coefficient with a recent theoretical formula obtained by Dagdug and Pineda Y1 - 2014 U6 - https://doi.org/10.1039/c3cp55160a SN - 1463-9076 SN - 1463-9084 VL - 16 IS - 13 SP - 6118 EP - 6128 PB - Royal Society of Chemistry CY - Cambridge 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 : a journal of European Chemical Societies 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 < x(2)(t) similar or equal to 2K(t)t with K(t) similar or equal to t(alpha-1) for 0 < alpha < 2. SBM may provide a seemingly adequate description in the case of unbounded diffusion, for which its probability density function coincides with that of fractional Brownian motion. Here we show that free SBM is weakly non-ergodic but does not exhibit a significant amplitude scatter of the time averaged mean squared displacement. More severely, we demonstrate that under confinement, the dynamics encoded by SBM is fundamentally different from both fractional Brownian motion and continuous time random walks. SBM is highly non-stationary and cannot provide a physical description for particles in a thermalised stationary system. Our findings have direct impact on the modelling of single particle tracking experiments, in particular, under confinement inside cellular compartments or when optical tweezers tracking methods are used. Y1 - 2014 U6 - https://doi.org/10.1039/c4cp02019g SN - 1463-9076 SN - 1463-9084 VL - 16 IS - 30 SP - 15811 EP - 15817 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - de Carvalho, Sidney J. A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Critical adsorption of polyelectrolytes onto charged Janus nanospheres JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - Based on extensive Monte Carlo simulations and analytical considerations we study the electrostatically driven adsorption of flexible polyelectrolyte chains onto charged Janus nanospheres. These net-neutral colloids are composed of two equally but oppositely charged hemispheres. The critical binding conditions for polyelectrolyte chains are analysed as function of the radius of the Janus particle and its surface charge density, as well as the salt concentration in the ambient solution. Specifically for the adsorption of finite-length polyelectrolyte chains onto Janus nanoparticles, we demonstrate that the critical adsorption conditions drastically differ when the size of the Janus particle or the screening length of the electrolyte are varied. We compare the scaling laws obtained for the adsorption-desorption threshold to the known results for uniformly charged spherical particles, observing significant disparities. We also contrast the changes to the polyelectrolyte chain conformations close to the surface of the Janus nanoparticles as compared to those for simple spherical particles. Finally, we discuss experimentally relevant physicochemical systems for which our simulations results may become important. In particular, we observe similar trends with polyelectrolyte complexation with oppositely but heterogeneously charged proteins. Y1 - 2014 U6 - https://doi.org/10.1039/c4cp02207f SN - 1463-9076 SN - 1463-9084 VL - 16 IS - 29 SP - 15539 EP - 15550 PB - Royal Society of Chemistry CY - Cambridge ER -