TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - A molecular relay race: sequential first-passage events to the terminal reaction centre in a cascade of diffusion controlled processes JF - New Journal of Physics (NJP) N2 - We consider a sequential cascade of molecular first-reaction events towards a terminal reaction centre in which each reaction step is controlled by diffusive motion of the particles. The model studied here represents a typical reaction setting encountered in diverse molecular biology systems, in which, e.g. a signal transduction proceeds via a series of consecutive 'messengers': the first messenger has to find its respective immobile target site triggering a launch of the second messenger, the second messenger seeks its own target site and provokes a launch of the third messenger and so on, resembling a relay race in human competitions. For such a molecular relay race taking place in infinite one-, two- and three-dimensional systems, we find exact expressions for the probability density function of the time instant of the terminal reaction event, conditioned on preceding successful reaction events on an ordered array of target sites. The obtained expressions pertain to the most general conditions: number of intermediate stages and the corresponding diffusion coefficients, the sizes of the target sites, the distances between them, as well as their reactivities are arbitrary. KW - diffusion KW - reaction cascade KW - first passage time Y1 - 2021 U6 - https://doi.org/10.1088/1367-2630/ac1e42 SN - 1367-2630 VL - 23 PB - IOP - Institute of Physics Publishing CY - Bristol ER - TY - JOUR A1 - Herrmann, Carl J. J. A1 - Metzler, Ralf A1 - Engbert, Ralf T1 - A self-avoiding walk with neural delays as a model of fixational eye movements JF - Scientific reports N2 - Fixational eye movements show scaling behaviour of the positional mean-squared displacement with a characteristic transition from persistence to antipersistence for increasing time-lag. These statistical patterns were found to be mainly shaped by microsaccades (fast, small-amplitude movements). However, our re-analysis of fixational eye-movement data provides evidence that the slow component (physiological drift) of the eyes exhibits scaling behaviour of the mean-squared displacement that varies across human participants. These results suggest that drift is a correlated movement that interacts with microsaccades. Moreover, on the long time scale, the mean-squared displacement of the drift shows oscillations, which is also present in the displacement auto-correlation function. This finding lends support to the presence of time-delayed feedback in the control of drift movements. Based on an earlier non-linear delayed feedback model of fixational eye movements, we propose and discuss different versions of a new model that combines a self-avoiding walk with time delay. As a result, we identify a model that reproduces oscillatory correlation functions, the transition from persistence to antipersistence, and microsaccades. Y1 - 2017 U6 - https://doi.org/10.1038/s41598-017-13489-8 SN - 2045-2322 VL - 7 PB - Nature Publ. Group CY - London ER - TY - JOUR A1 - Herrmann, Carl J. J. A1 - Metzler, Ralf T1 - A self-avoiding walk with neural delays as a model of fixational eye movements JF - Scientific reports N2 - Fixational eye movements show scaling behaviour of the positional mean-squared displacement with a characteristic transition from persistence to antipersistence for increasing time-lag. These statistical patterns were found to be mainly shaped by microsaccades (fast, small-amplitude movements). However, our re-analysis of fixational eye-movement data provides evidence that the slow component (physiological drift) of the eyes exhibits scaling behaviour of the mean-squared displacement that varies across human participants. These results suggest that drift is a correlated movement that interacts with microsaccades. Moreover, on the long time scale, the mean-squared displacement of the drift shows oscillations, which is also present in the displacement auto-correlation function. This finding lends support to the presence of time-delayed feedback in the control of drift movements. Based on an earlier non-linear delayed feedback model of fixational eye movements, we propose and discuss different versions of a new model that combines a self-avoiding walk with time delay. As a result, we identify a model that reproduces oscillatory correlation functions, the transition from persistence to antipersistence, and microsaccades. Y1 - 2017 U6 - https://doi.org/10.1038/s41598-017-13489-8 SN - 2045-2322 VL - 7 SP - 1 EP - 17 PB - Springer Nature CY - London ER - TY - JOUR A1 - Schwarzl, Maria A1 - Godec, Aljaz A1 - Oshanin, Gleb A1 - Metzler, Ralf T1 - A single predator charging a herd of prey: effects of self volume and predator-prey decision-making JF - Journal of physics : A, Mathematical and theoretical N2 - We study the degree of success of a single predator hunting a herd of prey on a two-dimensional square lattice landscape. We explicitly consider the self volume of the prey restraining their dynamics on the lattice. The movement of both predator and prey is chosen to include an intelligent, decision making step based on their respective sighting ranges, the radius in which they can detect the other species (prey cannot recognise each other besides the self volume interaction): after spotting each other the motion of prey and predator turns from a nearest neighbour random walk into directed escape or chase, respectively. We consider a large range of prey densities and sighting ranges and compute the mean first passage time for a predator to catch a prey as well as characterise the effective dynamics of the hunted prey. We find that the prey's sighting range dominates their life expectancy and the predator profits more from a bad eyesight of the prey than from his own good eye sight. We characterise the dynamics in terms of the mean distance between the predator and the nearest prey. It turns out that effectively the dynamics of this distance coordinate can be captured in terms of a simple Ornstein–Uhlenbeck picture. Reducing the many-body problem to a simple two-body problem by imagining predator and nearest prey to be connected by an effective Hookean bond, all features of the model such as prey density and sighting ranges merge into the effective binding constant. KW - first passage process KW - diffusion KW - predator-prey model Y1 - 2016 U6 - https://doi.org/10.1088/1751-8113/49/22/225601 SN - 1751-8113 SN - 1751-8121 VL - 49 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Sereshki, L. E. A1 - Lomholt, M. A. A1 - Metzler, Ralf T1 - A solution to the subdiffusion-efficiency paradox inactive states enhance reaction efficiency at subdiffusion conditions in living cells JF - epl : a letters journal exploring the frontiers of physics N2 - Macromolecular crowding in living biological cells effects subdiffusion of larger biomolecules such as proteins and enzymes. Mimicking this subdiffusion in terms of random walks on a critical percolation cluster, we here present a case study of EcoRV restriction enzymes involved in vital cellular defence. We show that due to its so far elusive propensity to an inactive state the enzyme avoids non-specific binding and remains well-distributed in the bulk cytoplasm of the cell. Despite the reduced volume exploration capability of subdiffusion processes, this mechanism guarantees a high efficiency of the enzyme. By variation of the non-specific binding constant and the bond occupation probability on the percolation network, we demonstrate that reduced nonspecific binding are beneficial for efficient subdiffusive enzyme activity even in relatively small bacteria cells. Our results corroborate a more local picture of cellular regulation. Y1 - 2012 U6 - https://doi.org/10.1209/0295-5075/97/20008 SN - 0295-5075 VL - 97 IS - 2 PB - EDP Sciences CY - Mulhouse ER - TY - JOUR A1 - Guggenberger, Tobias A1 - Chechkin, Aleksei A1 - Metzler, Ralf T1 - Absence of stationary states and non-Boltzmann distributions of fractional Brownian motion in shallow external potentials JF - New journal of physics : the open-access journal for physics N2 - We study the diffusive motion of a particle in a subharmonic potential of the form U(x) = |x|( c ) (0 < c < 2) driven by long-range correlated, stationary fractional Gaussian noise xi ( alpha )(t) with 0 < alpha <= 2. In the absence of the potential the particle exhibits free fractional Brownian motion with anomalous diffusion exponent alpha. While for an harmonic external potential the dynamics converges to a Gaussian stationary state, from extensive numerical analysis we here demonstrate that stationary states for shallower than harmonic potentials exist only as long as the relation c > 2(1 - 1/alpha) holds. We analyse the motion in terms of the mean squared displacement and (when it exists) the stationary probability density function. Moreover we discuss analogies of non-stationarity of Levy flights in shallow external potentials. KW - diffusion KW - Boltzmann distribution KW - fractional Brownian motion Y1 - 2022 U6 - https://doi.org/10.1088/1367-2630/ac7b3c SN - 1367-2630 VL - 24 IS - 7 PB - Dt. Physikalische Ges. CY - [Bad Honnef] ER - TY - JOUR A1 - Kar, Prathitha A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Acceleration of bursty multiprotein target search kinetics on DNA by colocalisation JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - Proteins are capable of locating specific targets on DNA by employing a facilitated diffusion process with intermittent 1D and 3D search steps. Gene colocalisation and coregulation-i.e. the spatial proximity of two communicating genes-is one factor capable of accelerating the target search process along the DNA. We perform Monte Carlo computer simulations and demonstrate the benefits of gene colocalisation for minimising the search time in a model DNA-protein system. We use a simple diffusion model to mimic the search for targets by proteins, produced initially in bursts of multiple proteins and performing the first-passage search on the DNA chain. The behaviour of the mean first-passage times to the target is studied as a function of distance between the initial position of proteins and the DNA target position, as well as versus the concentration of proteins. We also examine the properties of bursty target search kinetics for varying physical-chemical protein-DNA binding affinity. Our findings underline the relevance of colocalisation of production and binding sites for protein search inside biological cells. Y1 - 2017 U6 - https://doi.org/10.1039/c7cp06922g SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 12 SP - 7931 EP - 7946 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Godec, Aljaz A1 - Metzler, Ralf T1 - Active transport improves the precision of linear long distance molecular signalling JF - Journal of physics : A, Mathematical and theoretical N2 - Molecular signalling in living cells occurs at low copy numbers and is thereby inherently limited by the noise imposed by thermal diffusion. The precision at which biochemical receptors can count signalling molecules is intimately related to the noise correlation time. In addition to passive thermal diffusion, messenger RNA and vesicle-engulfed signalling molecules can transiently bind to molecular motors and are actively transported across biological cells. Active transport is most beneficial when trafficking occurs over large distances, for instance up to the order of 1 metre in neurons. Here we explain how intermittent active transport allows for faster equilibration upon a change in concentration triggered by biochemical stimuli. Moreover, we show how intermittent active excursions induce qualitative changes in the noise in effectively one-dimensional systems such as dendrites. Thereby they allow for significantly improved signalling precision in the sense of a smaller relative deviation in the concentration read-out by the receptor. On the basis of linear response theory we derive the exact mean field precision limit for counting actively transported molecules. We explain how intermittent active excursions disrupt the recurrence in the molecular motion, thereby facilitating improved signalling accuracy. Our results provide a deeper understanding of how recurrence affects molecular signalling precision in biological cells and novel medical-diagnostic devices. KW - noise in biochemical signalling KW - Brownian motion KW - active transport KW - linear response theory KW - fluctuation-dissipation theorem KW - generalised Langevin equation KW - recurrence Y1 - 2016 U6 - https://doi.org/10.1088/1751-8113/49/36/364001 SN - 1751-8113 SN - 1751-8121 VL - 49 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Ageing and confinement in non-ergodic heterogeneous diffusion processes JF - Journal of physics : A, Mathematical and theoretical N2 - We study the effects of ageing-the time delay between initiation of the physical process at t = 0 and start of observation at some time t(a) > 0-and spatial confinement on the properties of heterogeneous diffusion processes (HDPs) with deterministic power-law space-dependent diffusivities, D(x) = D-0 vertical bar x vertical bar(alpha). From analysis of the ensemble and time averaged mean squared displacements and the ergodicity breaking parameter quantifying the inherent degree of irreproducibility of individual realizations of the HDP we obtain striking similarities to ageing subdiffusive continuous time random walks with scale-free waiting time distributions. We also explore how both processes can be distinguished. For confined HDPs we study the long-time saturation of the ensemble and time averaged particle displacements as well as the magnitude of the inherent scatter of time averaged displacements and contrast the outcomes to the results known for other anomalous diffusion processes under confinement. KW - stochastic processes KW - anomalous diffusion KW - ageing KW - weak ergodicity breaking Y1 - 2014 U6 - https://doi.org/10.1088/1751-8113/47/48/485002 SN - 1751-8113 SN - 1751-8121 VL - 47 IS - 48 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Chechkin, Aleksei V. A1 - Kantz, Holger A1 - Metzler, Ralf T1 - Ageing effects in ultraslow continuous time random walks JF - The European physical journal : B, Condensed matter and complex systems N2 - In ageing systems physical observables explicitly depend on the time span elapsing between the original initiation of the system and the actual start of the recording of the particle motion. We here study the signatures of ageing in the framework of ultraslow continuous time random walk processes with super-heavy tailed waiting time densities. We derive the density for the forward or recurrent waiting time of the motion as function of the ageing time, generalise the Montroll-Weiss equation for this process, and analyse the ageing behaviour of the ensemble and time averaged mean squared displacements. Y1 - 2017 U6 - https://doi.org/10.1140/epjb/e2017-80270-9 SN - 1434-6028 SN - 1434-6036 VL - 90 PB - Springer CY - New York ER - TY - JOUR A1 - Krüsemann, Henning A1 - Godec, Aljaz A1 - Metzler, Ralf T1 - Ageing first passage time density in continuous time random walks and quenched energy landscapes JF - Journal of physics : A, Mathematical and theoretical N2 - We study the first passage dynamics of an ageing stochastic process in the continuous time random walk (CTRW) framework. In such CTRW processes the test particle performs a random walk, in which successive steps are separated by random waiting times distributed in terms of the waiting time probability density function Psi (t) similar or equal to t(-1-alpha) (0 <= alpha <= 2). An ageing stochastic process is defined by the explicit dependence of its dynamic quantities on the ageing time t(a), the time elapsed between its preparation and the start of the observation. Subdiffusive ageing CTRWs with 0 < alpha < 1 describe systems such as charge carriers in amorphous semiconducters, tracer dispersion in geological and biological systems, or the dynamics of blinking quantum dots. We derive the exact forms of the first passage time density for an ageing subdiffusive CTRW in the semi-infinite, confined, and biased case, finding different scaling regimes for weakly, intermediately, and strongly aged systems: these regimes, with different scaling laws, are also found when the scaling exponent is in the range 1 < alpha < 2, for sufficiently long ta. We compare our results with the ageing motion of a test particle in a quenched energy landscape. We test our theoretical results in the quenched landscape against simulations: only when the bias is strong enough, the correlations from returning to previously visited sites become insignificant and the results approach the ageing CTRW results. With small bias or without bias, the ageing effects disappear and a change in the exponent compared to the case of a completely annealed landscape can be found, reflecting the build-up of correlations in the quenched landscape. KW - first passage KW - random walks KW - anomalous diffusion Y1 - 2015 U6 - https://doi.org/10.1088/1751-8113/48/28/285001 SN - 1751-8113 SN - 1751-8121 VL - 48 IS - 28 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Kruesemann, Henning A1 - Schwarzl, Richard A1 - Metzler, Ralf T1 - Ageing Scher-Montroll Transport JF - Transport in Porous Media N2 - We study the properties of ageing Scher-Montroll transport in terms of a biased subdiffusive continuous time random walk in which the waiting times between consecutive jumps of the charge carriers are distributed according to the power law probability with . As we show, the dynamical properties of the Scher-Montroll transport depend on the ageing time span between the initial preparation of the system and the start of the observation. The Scher-Montroll transport theory was originally shown to describe the photocurrent in amorphous solids in the presence of an external electric field, but it has since been used in many other fields of physical sciences, in particular also in the geophysical context for the description of the transport of tracer particles in subsurface aquifers. In the absence of ageing () the photocurrent of the classical Scher-Montroll model or the breakthrough curves in the groundwater context exhibit a crossover between two power law regimes in time with the scaling exponents and . In the presence of ageing a new power law regime and an initial plateau regime of the current emerge. We derive the different power law regimes and crossover times of the ageing Scher-Montroll transport and show excellent agreement with simulations of the process. Experimental data of ageing Scher-Montroll transport in polymeric semiconductors are shown to agree well with the predictions of our theory. KW - Anomalous diffusion KW - Ageing KW - Scher-Montroll transport Y1 - 2016 U6 - https://doi.org/10.1007/s11242-016-0686-y SN - 0169-3913 SN - 1573-1634 VL - 115 SP - 327 EP - 344 PB - Springer CY - New York ER - TY - JOUR A1 - Metzler, Ralf A1 - Sanders, L. A1 - Lomholt, M. A. A1 - Lizana, L. A1 - Fogelmark, K. A1 - Ambjoernsson, Tobias T1 - Ageing single file motion JF - The European physical journal Y1 - 2014 U6 - https://doi.org/10.1140/epjst/e2014-02333-5 SN - 1951-6355 SN - 1951-6401 VL - 223 IS - 14 SP - 3287 EP - 3293 PB - Springer CY - Heidelberg 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 - Schulz, Johannes H. P. A1 - Barkai, Eli A1 - Metzler, Ralf T1 - Aging renewal theory and application to random walks JF - Physical review : X, Expanding access N2 - We discuss a renewal process in which successive events are separated by scale-free waiting time periods. Among other ubiquitous long-time properties, this process exhibits aging: events counted initially in a time interval [0, t] statistically strongly differ from those observed at later times [t(a,) t(a) + t]. The versatility of renewal theory is owed to its abstract formulation. Renewals can be interpreted as steps of a random walk, switching events in two-state models, domain crossings of a random motion, etc. In complex, disordered media, processes with scale-free waiting times play a particularly prominent role. We set up a unified analytical foundation for such anomalous dynamics by discussing in detail the distribution of the aging renewal process. We analyze its half-discrete, half-continuous nature and study its aging time evolution. These results are readily used to discuss a scale-free anomalous diffusion process, the continuous-time random walk. By this, we not only shed light on the profound origins of its characteristic features, such as weak ergodicity breaking, along the way, we also add an extended discussion on aging effects. In particular, we find that the aging behavior of time and ensemble averages is conceptually very distinct, but their time scaling is identical at high ages. Finally, we show how more complex motion models are readily constructed on the basis of aging renewal dynamics. Y1 - 2014 U6 - https://doi.org/10.1103/PhysRevX.4.011028 SN - 2160-3308 VL - 4 IS - 1 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Safdari, Hadiseh A1 - Chechkin, Aleksei V. A1 - Jafari, Gholamreza R. A1 - Metzler, Ralf T1 - Aging scaled Brownian motion JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Scaled Brownian motion (SBM) is widely used to model anomalous diffusion of passive tracers in complex and biological systems. It is a highly nonstationary process governed by the Langevin equation for Brownian motion, however, with a power-law time dependence of the noise strength. Here we study the aging properties of SBM for both unconfined and confined motion. Specifically, we derive the ensemble and time averaged mean squared displacements and analyze their behavior in the regimes of weak, intermediate, and strong aging. A very rich behavior is revealed for confined aging SBM depending on different aging times and whether the process is sub- or superdiffusive. We demonstrate that the information on the aging factorizes with respect to the lag time and exhibits a functional form that is identical to the aging behavior of scale-free continuous time random walk processes. While SBM exhibits a disparity between ensemble and time averaged observables and is thus weakly nonergodic, strong aging is shown to effect a convergence of the ensemble and time averaged mean squared displacement. Finally, we derive the density of first passage times in the semi-infinite domain that features a crossover defined by the aging time. Y1 - 2015 U6 - https://doi.org/10.1103/PhysRevE.91.042107 SN - 1539-3755 SN - 1550-2376 VL - 91 IS - 4 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Safdari, Hadiseh A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Bodrova, Anna A1 - Metzler, Ralf T1 - Aging underdamped scaled Brownian motion BT - Ensemble- and time-averaged particle displacements, nonergodicity, and the failure of the overdamping approximation JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We investigate both analytically and by computer simulations the ensemble- and time-averaged, nonergodic, and aging properties of massive particles diffusing in a medium with a time dependent diffusivity. We call this stochastic diffusion process the (aging) underdamped scaled Brownian motion (UDSBM). We demonstrate how the mean squared displacement (MSD) and the time-averaged MSD of UDSBM are affected by the inertial term in the Langevin equation, both at short, intermediate, and even long diffusion times. In particular, we quantify the ballistic regime for the MSD and the time-averaged MSD as well as the spread of individual time-averaged MSD trajectories. One of the main effects we observe is that, both for the MSD and the time-averaged MSD, for superdiffusive UDSBM the ballistic regime is much shorter than for ordinary Brownian motion. In contrast, for subdiffusive UDSBM, the ballistic region extends to much longer diffusion times. Therefore, particular care needs to be taken under what conditions the overdamped limit indeed provides a correct description, even in the long time limit. We also analyze to what extent ergodicity in the Boltzmann-Khinchin sense in this nonstationary system is broken, both for subdiffusive and superdiffusive UDSBM. Finally, the limiting case of ultraslow UDSBM is considered, with a mixed logarithmic and power-law dependence of the ensemble-and time-averaged MSDs of the particles. In the limit of strong aging, remarkably, the ordinary UDSBM and the ultraslow UDSBM behave similarly in the short time ballistic limit. The approaches developed here open ways for considering other stochastic processes under physically important conditions when a finite particle mass and aging in the system cannot be neglected. Y1 - 2017 U6 - https://doi.org/10.1103/PhysRevE.95.012120 SN - 2470-0045 SN - 2470-0053 VL - 95 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Javanainen, Matti A1 - Hammaren, Henrik A1 - Monticelli, Luca A1 - Jeon, Jae-Hyung A1 - Miettinen, Markus S. A1 - Martinez-Seara, Hector A1 - Metzler, Ralf A1 - Vattulainen, Ilpo T1 - Anomalous and normal diffusion of proteins and lipids in crowded lipid membranes JF - Faraday discussions N2 - Lateral diffusion plays a crucial role in numerous processes that take place in cell membranes, yet it is quite poorly understood in native membranes characterized by, e.g., domain formation and large concentration of proteins. In this article, we use atomistic and coarse-grained simulations to consider how packing of membranes and crowding with proteins affect the lateral dynamics of lipids and membrane proteins. We find that both packing and protein crowding have a profound effect on lateral diffusion, slowing it down. Anomalous diffusion is observed to be an inherent property in both protein-free and protein-rich membranes, and the time scales of anomalous diffusion and the exponent associated with anomalous diffusion are found to strongly depend on packing and crowding. Crowding with proteins also has a striking effect on the decay rate of dynamical correlations associated with lateral single-particle motion, as the transition from anomalous to normal diffusion is found to take place at macroscopic time scales: while in protein-poor conditions normal diffusion is typically observed in hundreds of nanoseconds, in protein-rich conditions the onset of normal diffusion is tens of microseconds, and in the most crowded systems as large as milliseconds. The computational challenge which results from these time scales is not easy to deal with, not even in coarse-grained simulations. We also briefly discuss the physical limits of protein motion. Our results suggest that protein concentration is anything but constant in the plane of cell membranes. Instead, it is strongly dependent on proteins' preference for aggregation. Y1 - 2013 U6 - https://doi.org/10.1039/c2fd20085f SN - 1359-6640 VL - 161 IS - 1 SP - 397 EP - 417 PB - Royal Society of Chemistry CY - Cambridge 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 - 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 - 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 - 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 - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Anomalous diffusion in time-fluctuating non-stationary diffusivity landscapes JF - Physical chemistry, chemical physics : 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-9076 SN - 1463-9084 VL - 18 SP - 23840 EP - 23852 PB - Royal Society of Chemistry CY - Cambridge 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 - JOUR A1 - Metzler, Ralf A1 - Jeon, Jae-Hyung A1 - Cherstvy, Andrey G. A1 - Barkai, Eli T1 - Anomalous diffusion models and their properties: non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking JF - Physical chemistry, chemical physics : a journal of European Chemical Societies 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. Y1 - 2014 U6 - https://doi.org/10.1039/c4cp03465a SN - 1463-9076 SN - 1463-9084 VL - 16 IS - 44 SP - 24128 EP - 24164 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Jeon, Jae-Hyung A1 - Monne, Hector Martinez-Seara A1 - Javanainen, Matti A1 - Metzler, Ralf T1 - Anomalous diffusion of phospholipids and cholesterols in a lipid bilayer and its origins JF - Physical review letters N2 - Combining extensive molecular dynamics simulations of lipid bilayer systems of varying chemical compositions with single-trajectory analyses, we systematically elucidate the stochastic nature of the lipid motion. We observe subdiffusion over more than 4 orders of magnitude in time, clearly stretching into the submicrosecond domain. The lipid motion depends on the lipid chemistry, the lipid phase, and especially the presence of cholesterol. We demonstrate that fractional Langevin equation motion universally describes the lipid motion in all phases, including the gel phase, and in the presence of cholesterol. The results underline the relevance of anomalous diffusion in lipid bilayers and the strong effects of the membrane composition. Y1 - 2012 U6 - https://doi.org/10.1103/PhysRevLett.109.188103 SN - 0031-9007 VL - 109 IS - 18 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Wang, Wei A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Anomalous diffusion, aging, and nonergodicity of scaled Brownian motion with fractional Gaussian noise: overview of related experimental observations and models JF - Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies N2 - How does a systematic time-dependence of the diffusion coefficient D(t) affect the ergodic and statistical characteristics of fractional Brownian motion (FBM)? Here, we answer this question via studying the characteristics of a set of standard statistical quantifiers relevant to single-particle-tracking (SPT) experiments. We examine, for instance, how the behavior of the ensemble- and time-averaged mean-squared displacements-denoted as the standard MSD < x(2)(Delta)> and TAMSD <<(delta(2)(Delta))over bar>> quantifiers-of FBM featuring < x(2) (Delta >> = <<(delta(2)(Delta >)over bar>> proportional to Delta(2H) (where H is the Hurst exponent and Delta is the [lag] time) changes in the presence of a power-law deterministically varying diffusivity D-proportional to(t) proportional to t(alpha-1) -germane to the process of scaled Brownian motion (SBM)-determining the strength of fractional Gaussian noise. The resulting compound "scaled-fractional" Brownian motion or FBM-SBM is found to be nonergodic, with < x(2)(Delta >> proportional to Delta(alpha+)(2H)(-1) and <(delta 2(Delta >) over bar > proportional to Delta(2H). We also detect a stalling behavior of the MSDs for very subdiffusive SBM and FBM, when alpha + 2H - 1 < 0. The distribution of particle displacements for FBM-SBM remains Gaussian, as that for the parent processes of FBM and SBM, in the entire region of scaling exponents (0 < alpha < 2 and 0 < H < 1). The FBM-SBM process is aging in a manner similar to SBM. The velocity autocorrelation function (ACF) of particle increments of FBM-SBM exhibits a dip when the parent FBM process is subdiffusive. Both for sub- and superdiffusive FBM contributions to the FBM-SBM process, the SBM exponent affects the long-time decay exponent of the ACF. Applications of the FBM-SBM-amalgamated process to the analysis of SPT data are discussed. A comparative tabulated overview of recent experimental (mainly SPT) and computational datasets amenable for interpretation in terms of FBM-, SBM-, and FBM-SBM-like models of diffusion culminates the presentation. The statistical aspects of the dynamics of a wide range of biological systems is compared in the table, from nanosized beads in living cells, to chromosomal loci, to water diffusion in the brain, and, finally, to patterns of animal movements. Y1 - 2022 U6 - https://doi.org/10.1039/d2cp01741e SN - 1463-9076 SN - 1463-9084 VL - 24 IS - 31 SP - 18482 EP - 18504 PB - RSC Publ. CY - Cambridge ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Safdari, Hadiseh A1 - Metzler, Ralf T1 - Anomalous diffusion, nonergodicity, and ageing for exponentially and logarithmically time-dependent diffusivity BT - striking differences for massive versus massless particles JF - Journal of physics. D, Applied physics N2 - We investigate a diffusion process with a time-dependent diffusion coefficient, both exponentially increasing and decreasing in time, D(t)=D-0(e +/- 2 alpha t). For this (hypothetical) nonstationary diffusion process we compute-both analytically and from extensive stochastic simulations-the behavior of the ensemble- and time-averaged mean-squared displacements (MSDs) of the particles, both in the over- and underdamped limits. Simple asymptotic relations derived for the short- and long-time behaviors are shown to be in excellent agreement with the results of simulations. The diffusive characteristics in the presence of ageing are also considered, with dramatic differences of the over- versus underdamped regime. Our results for D(t)=D-0(e +/- 2 alpha t) extend and generalize the class of diffusive systems obeying scaled Brownian motion featuring a power-law-like variation of the diffusivity with time, D(t) similar to t(alpha-1). We also examine the logarithmically increasing diffusivity, D(t)=D(0)log[t/tau(0)], as another fundamental functional dependence (in addition to the power-law and exponential) and as an example of diffusivity slowly varying in time. One of the main conclusions is that the behavior of the massive particles is predominantly ergodic, while weak ergodicity breaking is repeatedly found for the time-dependent diffusion of the massless particles at short times. The latter manifests itself in the nonequivalence of the (both nonaged and aged) MSD and the mean time-averaged MSD. The current findings are potentially applicable to a class of physical systems out of thermal equilibrium where a rapid increase or decrease of the particles' diffusivity is inherently realized. One biological system potentially featuring all three types of time-dependent diffusion (power-law-like, exponential, and logarithmic) is water diffusion in the brain tissues, as we thoroughly discuss in the end. KW - anomalous diffusion KW - scaled Brownian motion KW - stochastic processes KW - nonstationary diffusivity KW - water diffusion in the brain KW - nonergodicity Y1 - 2021 U6 - https://doi.org/10.1088/1361-6463/abdff0 SN - 0022-3727 SN - 1361-6463 VL - 54 IS - 19 PB - IOP Publ. Ltd. CY - Bristol 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 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Grebenkov, Denis S. A1 - Metzler, Ralf T1 - Anomalous, non-Gaussian tracer diffusion in crowded two-dimensional environments JF - NEW JOURNAL OF PHYSICS N2 - A topic of intense current investigation pursues the question of how the highly crowded environment of biological cells affects the dynamic properties of passively diffusing particles. Motivated by recent experiments we report results of extensive simulations of the motion of a finite sized tracer particle in a heterogeneously crowded environment made up of quenched distributions of monodisperse crowders of varying sizes in finite circular two-dimensional domains. For given spatial distributions of monodisperse crowders we demonstrate how anomalous diffusion with strongly non-Gaussian features arises in this model system. We investigate both biologically relevant situations of particles released either at the surface of an inner domain or at the outer boundary, exhibiting distinctly different features of the observed anomalous diffusion for heterogeneous distributions of crowders. Specifically we reveal an asymmetric spreading of tracers even at moderate crowding. In addition to the mean squared displacement (MSD) and local diffusion exponent we investigate the magnitude and the amplitude scatter of the time averaged MSD of individual tracer trajectories, the non-Gaussianity parameter, and the van Hove correlation function. We also quantify how the average tracer diffusivity varies with the position in the domain with a heterogeneous radial distribution of crowders and examine the behaviour of the survival probability and the dynamics of the tracer survival probability. Inter alia, the systems we investigate are related to the passive transport of lipid molecules and proteins in two-dimensional crowded membranes or the motion in colloidal solutions or emulsions in effectively two-dimensional geometries, as well as inside supercrowded, surface adhered cells. KW - anomalous diffusion KW - crowded fluids KW - stochastic processes Y1 - 2016 U6 - https://doi.org/10.1088/1367-2630/18/1/013027 SN - 1367-2630 VL - 18 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Doerries, Timo J. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Apparent anomalous diffusion and non-Gaussian distributions in a simple mobile-immobile transport model with Poissonian switching JF - Interface : journal of the Royal Society N2 - We analyse mobile-immobile transport of particles that switch between the mobile and immobile phases with finite rates. Despite this seemingly simple assumption of Poissonian switching, we unveil a rich transport dynamics including significant transient anomalous diffusion and non-Gaussian displacement distributions. Our discussion is based on experimental parameters for tau proteins in neuronal cells, but the results obtained here are expected to be of relevance for a broad class of processes in complex systems. Specifically, we obtain that, when the mean binding time is significantly longer than the mean mobile time, transient anomalous diffusion is observed at short and intermediate time scales, with a strong dependence on the fraction of initially mobile and immobile particles. We unveil a Laplace distribution of particle displacements at relevant intermediate time scales. For any initial fraction of mobile particles, the respective mean squared displacement (MSD) displays a plateau. Moreover, we demonstrate a short-time cubic time dependence of the MSD for immobile tracers when initially all particles are immobile. KW - diffusion KW - mobile-immobile model KW - tau proteins Y1 - 2022 U6 - https://doi.org/10.1098/rsif.2022.0233 SN - 1742-5689 SN - 1742-5662 VL - 19 IS - 192 PB - Royal Society CY - London 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 - Padash, Amin A1 - Sandev, Trifce A1 - Kantz, Holger A1 - Metzler, Ralf A1 - Chechkin, Aleksei T1 - Asymmetric Levy flights are more efficient in random search JF - Fractal and fractional N2 - We study the first-arrival (first-hitting) dynamics and efficiency of a one-dimensional random search model performing asymmetric Levy flights by leveraging the Fokker-Planck equation with a delta-sink and an asymmetric space-fractional derivative operator with stable index alpha and asymmetry (skewness) parameter beta. We find exact analytical results for the probability density of first-arrival times and the search efficiency, and we analyse their behaviour within the limits of short and long times. We find that when the starting point of the searcher is to the right of the target, random search by Brownian motion is more efficient than Levy flights with beta <= 0 (with a rightward bias) for short initial distances, while for beta>0 (with a leftward bias) Levy flights with alpha -> 1 are more efficient. When increasing the initial distance of the searcher to the target, Levy flight search (except for alpha=1 with beta=0) is more efficient than the Brownian search. Moreover, the asymmetry in jumps leads to essentially higher efficiency of the Levy search compared to symmetric Levy flights at both short and long distances, and the effect is more pronounced for stable indices alpha close to unity. KW - asymmetric Levy flights KW - first-arrival density KW - search efficiency Y1 - 2022 U6 - https://doi.org/10.3390/fractalfract6050260 SN - 2504-3110 VL - 6 IS - 5 PB - MDPI CY - Basel ER - TY - JOUR A1 - Thapa, Samudrajit A1 - Lomholt, Michael Andersen A1 - Krog, Jens A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Bayesian analysis of single-particle tracking data using the nested-sampling algorithm: maximum-likelihood model selection applied to stochastic-diffusivity data JF - Physical chemistry, chemical physics : PCCP ; a journal of European Chemical Societies N2 - We employ Bayesian statistics using the nested-sampling algorithm to compare and rank multiple models of ergodic diffusion (including anomalous diffusion) as well as to assess their optimal parameters for in silico-generated and real time-series. We focus on the recently-introduced model of Brownian motion with "diffusing diffusivity'-giving rise to widely-observed non-Gaussian displacement statistics-and its comparison to Brownian and fractional Brownian motion, also for the time-series with some measurement noise. We conduct this model-assessment analysis using Bayesian statistics and the nested-sampling algorithm on the level of individual particle trajectories. We evaluate relative model probabilities and compute best-parameter sets for each diffusion model, comparing the estimated parameters to the true ones. We test the performance of the nested-sampling algorithm and its predictive power both for computer-generated (idealised) trajectories as well as for real single-particle-tracking trajectories. Our approach delivers new important insight into the objective selection of the most suitable stochastic model for a given time-series. We also present first model-ranking results in application to experimental data of tracer diffusion in polymer-based hydrogels. Y1 - 2018 U6 - https://doi.org/10.1039/c8cp04043e SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 46 SP - 29018 EP - 29037 PB - Royal Society of Chemistry CY - Cambridge 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 - Seckler, Henrik A1 - Metzler, Ralf T1 - Bayesian deep learning for error estimation in the analysis of anomalous diffusion JF - Nature Communications 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 diffusion model and the regression of the anomalous diffusion exponent of single-particle-trajectories. Evaluating their performance, we find that these models can achieve a well-calibrated 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.
Diffusive motions in complex environments such as living biological cells or soft matter systems can be analyzed with single-particle-tracking approaches, where accuracy of output may vary. The authors involve a machine-learning technique for decoding anomalous-diffusion data and provide an uncertainty estimate together with predicted output. Y1 - 2022 U6 - https://doi.org/10.1038/s41467-022-34305-6 SN - 2041-1723 VL - 13 IS - 1 PB - Nature portfolio CY - Berlin ER - TY - JOUR A1 - Thapa, Samudrajit A1 - Park, Seongyu A1 - Kim, Yeongjin A1 - Jeon, Jae-Hyung A1 - Metzler, Ralf A1 - Lomholt, Michael A. T1 - Bayesian inference of scaled versus fractional Brownian motion JF - Journal of physics : A, mathematical and theoretical N2 - We present a Bayesian inference scheme for scaled Brownian motion, and investigate its performance on synthetic data for parameter estimation and model selection in a combined inference with fractional Brownian motion. We include the possibility of measurement noise in both models. We find that for trajectories of a few hundred time points the procedure is able to resolve well the true model and parameters. Using the prior of the synthetic data generation process also for the inference, the approach is optimal based on decision theory. We include a comparison with inference using a prior different from the data generating one. KW - Bayesian inference KW - scaled Brownian motion KW - single particle tracking Y1 - 2022 U6 - https://doi.org/10.1088/1751-8121/ac60e7 SN - 1751-8113 SN - 1751-8121 VL - 55 IS - 19 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Sandev, Trifce A1 - Sokolov, Igor M. A1 - Metzler, Ralf A1 - Chechkin, Aleksei V. T1 - Beyond monofractional kinetics JF - Chaos, solitons & fractals : applications in science and engineering ; an interdisciplinary journal of nonlinear science N2 - We discuss generalized integro-differential diffusion equations whose integral kernels are not of a simple power law form, and thus these equations themselves do not belong to the family of fractional diffusion equations exhibiting a monoscaling behavior. They instead generate a broad class of anomalous nonscaling patterns, which correspond either to crossovers between different power laws, or to a non-power-law behavior as exemplified by the logarithmic growth of the width of the distribution. We consider normal and modified forms of these generalized diffusion equations and provide a brief discussion of three generic types of integral kernels for each form, namely, distributed order, truncated power law and truncated distributed order kernels. For each of the cases considered we prove the non-negativity of the solution of the corresponding generalized diffusion equation and calculate the mean squared displacement. (C) 2017 Elsevier Ltd. All rights reserved. KW - Distributed order diffusion-wave equations KW - Complete Bernstein function KW - Completely monotone function Y1 - 2017 U6 - https://doi.org/10.1016/j.chaos.2017.05.001 SN - 0960-0779 SN - 1873-2887 VL - 102 SP - 210 EP - 217 PB - Elsevier CY - Oxford 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 - Talukder, Srijeeta A1 - Sen, Shrabani A1 - Chakraborti, Prantik A1 - Metzler, Ralf A1 - Banik, Suman K. A1 - Chaudhury, Pinaki T1 - Breathing dynamics based parameter sensitivity analysis of hetero-polymeric DNA JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - We study the parameter sensitivity of hetero-polymeric DNA within the purview of DNA breathing dynamics. The degree of correlation between the mean bubble size and the model parameters is estimated for this purpose for three different DNA sequences. The analysis leads us to a better understanding of the sequence dependent nature of the breathing dynamics of hetero-polymeric DNA. Out of the 14 model parameters for DNA stability in the statistical Poland-Scheraga approach, the hydrogen bond interaction epsilon(hb)(AT) for an AT base pair and the ring factor. turn out to be the most sensitive parameters. In addition, the stacking interaction epsilon(st)(TA-TA) for an TA-TA nearest neighbor pair of base-pairs is found to be the most sensitive one among all stacking interactions. Moreover, we also establish that the nature of stacking interaction has a deciding effect on the DNA breathing dynamics, not the number of times a particular stacking interaction appears in a sequence. We show that the sensitivity analysis can be used as an effective measure to guide a stochastic optimization technique to find the kinetic rate constants related to the dynamics as opposed to the case where the rate constants are measured using the conventional unbiased way of optimization. (c) 2014 AIP Publishing LLC. Y1 - 2014 U6 - https://doi.org/10.1063/1.4869112 SN - 0021-9606 SN - 1089-7690 VL - 140 IS - 12 PB - American Institute of Physics CY - Melville 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 - Chechkin, Aleksei V. A1 - Seno, Flavio A1 - Metzler, Ralf A1 - Sokolov, Igor M. T1 - Brownian yet Non-Gaussian Diffusion: From Superstatistics to Subordination of Diffusing Diffusivities JF - Physical review : X, Expanding access N2 - A growing number of biological, soft, and active matter systems are observed to exhibit normal diffusive dynamics with a linear growth of the mean-squared displacement, yet with a non-Gaussian distribution of increments. Based on the Chubinsky-Slater idea of a diffusing diffusivity, we here establish and analyze a minimal model framework of diffusion processes with fluctuating diffusivity. In particular, we demonstrate the equivalence of the diffusing diffusivity process with a superstatistical approach with a distribution of diffusivities, at times shorter than the diffusivity correlation time. At longer times, a crossover to a Gaussian distribution with an effective diffusivity emerges. Specifically, we establish a subordination picture of Brownian but non-Gaussian diffusion processes, which can be used for a wide class of diffusivity fluctuation statistics. Our results are shown to be in excellent agreement with simulations and numerical evaluations. Y1 - 2017 U6 - https://doi.org/10.1103/PhysRevX.7.021002 SN - 2160-3308 VL - 7 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Emanuel, Marc D. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf A1 - Gompper, Gerhard T1 - Buckling transitions and soft-phase invasion of two-component icosahedral shells JF - Physical review / publ. by The American Physical Society. E, Statistical, nonlinear, and soft matter physics N2 - What is the optimal distribution of two types of crystalline phases on the surface of icosahedral shells, such as of many viral capsids? We here investigate the distribution of a thin layer of soft material on a crystalline convex icosahedral shell. We demonstrate how the shapes of spherical viruses can be understood from the perspective of elasticity theory of thin two-component shells. We develop a theory of shape transformations of an icosahedral shell upon addition of a softer, but still crystalline, material onto its surface. We show how the soft component "invades" the regions with the highest elastic energy and stress imposed by the 12 topological defects on the surface. We explore the phase diagram as a function of the surface fraction of the soft material, the shell size, and the incommensurability of the elastic moduli of the rigid and soft phases. We find that, as expected, progressive filling of the rigid shell by the soft phase starts from the most deformed regions of the icosahedron. With a progressively increasing soft-phase coverage, the spherical segments of domes are filled first (12 vertices of the shell), then the cylindrical segments connecting the domes (30 edges) are invaded, and, ultimately, the 20 flat faces of the icosahedral shell tend to be occupied by the soft material. We present a detailed theoretical investigation of the first two stages of this invasion process and develop a model of morphological changes of the cone structure that permits noncircular cross sections. In conclusion, we discuss the biological relevance of some structures predicted from our calculations, in particular for the shape of viral capsids. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevE.102.062104 SN - 2470-0045 SN - 2470-0053 SN - 2470-0061 SN - 1538-4519 VL - 102 IS - 6 PB - Woodbury CY - New York ER - TY - JOUR A1 - Chechkin, Aleksei V. A1 - Zaid, Irwin M. A1 - Lomholt, Michael A. A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Bulk-mediated diffusion on a planar surface full solution JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We consider the effective surface motion of a particle that intermittently unbinds from a planar surface and performs bulk excursions. Based on a random-walk approach, we derive the diffusion equations for surface and bulk diffusion including the surface-bulk coupling. From these exact dynamic equations, we analytically obtain the propagator of the effective surface motion. This approach allows us to deduce a superdiffusive, Cauchy-type behavior on the surface, together with exact cutoffs limiting the Cauchy form. Moreover, we study the long-time dynamics for the surface motion. Y1 - 2012 U6 - https://doi.org/10.1103/PhysRevE.86.041101 SN - 1539-3755 VL - 86 IS - 4 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Chechkin, Aleksei V. A1 - Zaid, I. M. A1 - Lomholt, M. A. A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Bulk-mediated surface diffusion on a cylinder in the fast exchange limit JF - Mathematical modelling of natural phenomena N2 - In various biological systems and small scale technological applications particles transiently bind to a cylindrical surface. Upon unbinding the particles diffuse in the vicinal bulk before rebinding to the surface. Such bulk-mediated excursions give rise to an effective surface translation, for which we here derive and discuss the dynamic equations, including additional surface diffusion. We discuss the time evolution of the number of surface-bound particles, the effective surface mean squared displacement, and the surface propagator. In particular, we observe sub- and superdiffusive regimes. A plateau of the surface mean-squared displacement reflects a stalling of the surface diffusion at longer times. Finally, the corresponding first passage problem for the cylindrical geometry is analysed. KW - Bulk-mediated diffusion KW - anomalous diffusion KW - Levy flights KW - stochastic processes Y1 - 2013 U6 - https://doi.org/10.1051/mmnp/20138208 SN - 0973-5348 VL - 8 IS - 2 SP - 114 EP - 126 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR 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 JF - Scientific Reports 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. Y1 - 2019 U6 - https://doi.org/10.1038/s41598-019-48525-2 SN - 2045-2322 VL - 9 PB - Macmillan Publishers Limited part of Springer Nature CY - London ER - TY - JOUR A1 - Varghese, Alan J. A1 - Chechkin, Aleksei A1 - Metzler, Ralf A1 - Sujith, Raman I. T1 - Capturing multifractality of pressure fluctuations in thermoacoustic systems using fractional-order derivatives JF - Chaos : an interdisciplinary journal of nonlinear science N2 - The stable operation of a turbulent combustor is not completely silent; instead, there is a background of small amplitude aperiodic acoustic fluctuations known as combustion noise. Pressure fluctuations during this state of combustion noise are multifractal due to the presence of multiple temporal scales that contribute to its dynamics. However, existing models are unable to capture the multifractality in the pressure fluctuations. We conjecture an underlying fractional dynamics for the thermoacoustic system and obtain a fractional-order model for pressure fluctuations. The data from this model has remarkable visual similarity to the experimental data and also has a wide multifractal spectrum during the state of combustion noise. Quantitative similarity with the experimental data in terms of the Hurst exponent and the multifractal spectrum is observed during the state of combustion noise. This model is also able to produce pressure fluctuations that are qualitatively similar to the experimental data acquired during intermittency and thermoacoustic instability. Furthermore, we argue that the fractional dynamics vanish as we approach the state of thermoacoustic instability. Y1 - 2021 U6 - https://doi.org/10.1063/5.0032585 SN - 1054-1500 SN - 1089-7682 VL - 31 IS - 3 PB - American Institute of Physics, AIP CY - Melville ER - TY - JOUR A1 - Mutothya, Nicholas Mwilu A1 - Xu, Yong A1 - Li, Yongge A1 - Metzler, Ralf T1 - Characterising stochastic motion in heterogeneous media driven by coloured non-Gaussian noise JF - Journal of physics : A, Mathematical and theoretical N2 - We study the stochastic motion of a test particle in a heterogeneous medium in terms of a position dependent diffusion coefficient mimicking measured deterministic diffusivity gradients in biological cells or the inherent heterogeneity of geophysical systems. Compared to previous studies we here investigate the effect of the interplay of anomalous diffusion effected by position dependent diffusion coefficients and coloured non-Gaussian noise. The latter is chosen to be distributed according to Tsallis' q-distribution, representing a popular example for a non-extensive statistic. We obtain the ensemble and time averaged mean squared displacements for this generalised process and establish its non-ergodic properties as well as analyse the non-Gaussian nature of the associated displacement distribution. We consider both non-stratified and stratified environments. KW - diffusion KW - anomalous diffusion KW - non-extensive statistics KW - coloured KW - noise KW - heterogeneous diffusion process Y1 - 2021 U6 - https://doi.org/10.1088/1751-8121/abfba6 SN - 1751-8113 SN - 1751-8121 VL - 54 IS - 29 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Vilk, Ohad A1 - Aghion, Erez A1 - Nathan, Ran A1 - Toledo, Sivan A1 - Metzler, Ralf A1 - Assaf, Michael T1 - Classification of anomalous diffusion in animal movement data using power spectral analysis JF - Journal of physics : A, Mathematical and theoretical N2 - The field of movement ecology has seen a rapid increase in high-resolution data in recent years, leading to the development of numerous statistical and numerical methods to analyse relocation trajectories. Data are often collected at the level of the individual and for long periods that may encompass a range of behaviours. Here, we use the power spectral density (PSD) to characterise the random movement patterns of a black-winged kite (Elanus caeruleus) and a white stork (Ciconia ciconia). The tracks are first segmented and clustered into different behaviours (movement modes), and for each mode we measure the PSD and the ageing properties of the process. For the foraging kite we find 1/f noise, previously reported in ecological systems mainly in the context of population dynamics, but not for movement data. We further suggest plausible models for each of the behavioural modes by comparing both the measured PSD exponents and the distribution of the single-trajectory PSD to known theoretical results and simulations. KW - diffusion KW - anomalous diffusion KW - power spectral analysis KW - ecological KW - movement data Y1 - 2022 U6 - https://doi.org/10.1088/1751-8121/ac7e8f SN - 1751-8113 SN - 1751-8121 VL - 55 IS - 33 PB - IOP Publishing CY - Bristol ER - TY - JOUR A1 - Awad, Emad A1 - Metzler, Ralf T1 - Closed-form multi-dimensional solutions and asymptotic behaviours for subdiffusive processes with crossovers: II. Accelerating case JF - Journal of physics : A, Mathematical and theoretical N2 - Anomalous diffusion with a power-law time dependence vertical bar R vertical bar(2)(t) similar or equal to t(alpha i) of the mean squared displacement occurs quite ubiquitously in numerous complex systems. Often, this anomalous diffusion is characterised by crossovers between regimes with different anomalous diffusion exponents alpha(i). Here we consider the case when such a crossover occurs from a first regime with alpha(1) to a second regime with alpha(2) such that alpha(2) > alpha(1), i.e., accelerating anomalous diffusion. A widely used framework to describe such crossovers in a one-dimensional setting is the bi-fractional diffusion equation of the so-called modified type, involving two time-fractional derivatives defined in the Riemann-Liouville sense. We here generalise this bi-fractional diffusion equation to higher dimensions and derive its multidimensional propagator (Green's function) for the general case when also a space fractional derivative is present, taking into consideration long-ranged jumps (Levy flights). We derive the asymptotic behaviours for this propagator in both the short- and long-time as well the short- and long-distance regimes. Finally, we also calculate the mean squared displacement, skewness and kurtosis in all dimensions, demonstrating that in the general case the non-Gaussian shape of the probability density function changes. KW - multidimensional fractional diffusion equation KW - continuous time random KW - walks KW - crossover anomalous diffusion dynamics KW - non-Gaussian probability KW - density Y1 - 2022 U6 - https://doi.org/10.1088/1751-8121/ac5a90 SN - 1751-8113 SN - 1751-8121 VL - 55 IS - 20 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Ślęzak, Jakub A1 - Metzler, Ralf A1 - Magdziarz, Marcin T1 - Codifference can detect ergodicity breaking and non-Gaussianity JF - New Journal of Physics 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. KW - diffusion KW - stochastic time series KW - anomalous diffusion Y1 - 2019 U6 - https://doi.org/10.1088/1367-2630/ab13f3 SN - 1367-2630 VL - 21 PB - Deutsche Physikalische Gesellschaft CY - Bad Honnef 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 - Sandev, Trifce A1 - Iomin, Alexander A1 - Kantz, Holger A1 - Metzler, Ralf A1 - Chechkin, Aleksei V. T1 - Comb Model with Slow and Ultraslow Diffusion JF - Mathematical modelling of natural phenomena N2 - We consider a generalized diffusion equation in two dimensions for modeling diffusion on a comb-like structures. We analyze the probability distribution functions and we derive the mean squared displacement in x and y directions. Different forms of the memory kernels (Dirac delta, power-law, and distributed order) are considered. It is shown that anomalous diffusion may occur along both x and y directions. Ultraslow diffusion and some more general diffusive processes are observed as well. We give the corresponding continuous time random walk model for the considered two dimensional diffusion-like equation on a comb, and we derive the probability distribution functions which subordinate the process governed by this equation to the Wiener process. KW - comb-like model KW - anomalous diffusion KW - mean squared displacement KW - probability density function Y1 - 2016 U6 - https://doi.org/10.1051/mmnp/201611302 SN - 0973-5348 SN - 1760-6101 VL - 11 SP - 18 EP - 33 PB - EDP Sciences CY - Les Ulis 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 - Blavatska, Viktoria A1 - Metzler, Ralf T1 - Conformational properties of complex polymers: rosette versus star-like structures JF - Journal of physics : A, Mathematical and theoretical N2 - Multiple loop formation in polymer macromolecules is an important feature of the chromatin organization and DNA compactification in the nuclei. We analyse the size and shape characteristics of complex polymer structures, containing in general f(1) loops (petals) and f(2) linear chains (branches). Within the frames of continuous model of Gaussian macromolecule, we apply the path integration method and obtain the estimates for gyration radius R-g and asphericity (A) over cap of typical conformation as functions of parameters f(1), f(2). In particular, our results qualitatively reveal the extent of anisotropy of star-like topologies as compared to the rosette structures of the same total molecular weight. KW - polymers KW - path integration KW - conformational properties Y1 - 2015 U6 - https://doi.org/10.1088/1751-8113/48/13/135001 SN - 1751-8113 SN - 1751-8121 VL - 48 IS - 13 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 - 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 - 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 - Magdziarz, Marcin A1 - Metzler, Ralf A1 - Szczotka, Wladyslaw A1 - Zebrowski, Piotr T1 - Correlated continuous-time random walks in external force fields JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We study the anomalous diffusion of a particle in an external force field whose motion is governed by nonrenewal continuous time random walks with correlated waiting times. In this model the current waiting time T-i is equal to the previous waiting time Ti-1 plus a small increment. Based on the associated coupled Langevin equations the force field is systematically introduced. We show that in a confining potential the relaxation dynamics follows power-law or stretched exponential pattern, depending on the model parameters. The process obeys a generalized Einstein-Stokes-Smoluchowski relation and observes the second Einstein relation. The stationary solution is of Boltzmann-Gibbs form. The case of an harmonic potential is discussed in some detail. We also show that the process exhibits aging and ergodicity breaking. Y1 - 2012 U6 - https://doi.org/10.1103/PhysRevE.85.051103 SN - 1539-3755 SN - 1550-2376 VL - 85 IS - 5 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Magdziarz, Marcin A1 - Metzler, Ralf A1 - Szczotka, Wladyslaw A1 - Zebrowski, Piotr T1 - Correlated continuous-time random walks-scaling limits and Langevin picture JF - Journal of statistical mechanics: theory and experiment N2 - In this paper we analyze correlated continuous-time random walks introduced recently by Tejedor and Metzler (2010 J. Phys. A: Math. Theor. 43 082002). We obtain the Langevin equations associated with this process and the corresponding scaling limits of their solutions. We prove that the limit processes are self-similar and display anomalous dynamics. Moreover, we extend the model to include external forces. Our results are confirmed by Monte Carlo simulations. KW - stochastic processes (theory) KW - diffusion Y1 - 2012 U6 - https://doi.org/10.1088/1742-5468/2012/04/P04010 SN - 1742-5468 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Sandev, Trifce A1 - Metzler, Ralf A1 - Tomovski, Zivorad T1 - Correlation functions for the fractional generalized Langevin equation in the presence of internal and external noise JF - Journal of mathematical physics N2 - We study generalized fractional Langevin equations in the presence of a harmonic potential. General expressions for the mean velocity and particle displacement, the mean squared displacement, position and velocity correlation functions, as well as normalized displacement correlation function are derived. We report exact results for the cases of internal and external friction, that is, when the driving noise is either internal and thus the fluctuation-dissipation relation is fulfilled or when the noise is external. The asymptotic behavior of the generalized stochastic oscillator is investigated, and the case of high viscous damping (overdamped limit) is considered. Additional behaviors of the normalized displacement correlation functions different from those for the regular damped harmonic oscillator are observed. In addition, the cases of a constant external force and the force free case are obtained. The validity of the generalized Einstein relation for this process is discussed. The considered fractional generalized Langevin equation may be used to model anomalous diffusive processes including single file-type diffusion. Y1 - 2014 U6 - https://doi.org/10.1063/1.4863478 SN - 0022-2488 SN - 1089-7658 VL - 55 IS - 2 PB - American Institute of Physics CY - Melville 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 - Caetano, Daniel L. Z. A1 - Carvalho, Sidney Jurado de A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Critical adsorption of multiple polyelectrolytes onto a nanosphere BT - splitting the adsorption-desorption transition boundary JF - Interface : journal of the Royal Society N2 - Employing extensive Monte Carlo computer simulations, we investigate in detail the properties of multichain adsorption of charged flexible polyelectrolytes (PEs) onto oppositely charged spherical nanoparticles (SNPs). We quantify the conditions of critical adsorption-the phase-separation curve between the adsorbed and desorbed states of the PEs-as a function of the SNP surface-charge density and the concentration of added salt. We study the degree of fluctuations of the PE-SNP electrostatic binding energy, which we use to quantify the emergence of the phase subtransitions, including a series of partially adsorbed PE configurations. We demonstrate how the phase-separation adsorption-desorption boundary shifts and splits into multiple subtransitions at low-salt conditions, thereby generalizing and extending the results for critical adsorption of a single PE onto the SNP. The current findings are relevant for finite concentrations of PEs around the attracting SNP, such as the conditions for PE adsorption onto globular proteins carrying opposite electric charges. KW - nanoparticles KW - polyelectrolytes KW - electrostatics KW - critical adsorption KW - phase-transition boundary Y1 - 2020 U6 - https://doi.org/10.1098/rsif.2020.0199 SN - 1742-5689 SN - 1742-5662 VL - 17 IS - 167 PB - Royal Society CY - London 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 - 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 - 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 - 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: the nonlinear Poisson-Boltzmann approach JF - NEW JOURNAL OF 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-Huckel approximation is often not feasible and the nonlinear Poisson-Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson-Boltzmann equation is smaller than the Debye-Huckel result, such that the required critical surface charge density for polyelectrolyte adsorption sigma(c) increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical-chemical and biophysical systems. 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. Ltd. CY - Bristol ER - TY - JOUR A1 - Awad, Emad A1 - Metzler, Ralf T1 - Crossover dynamics from superdiffusion to subdiffusion BT - models and solutions JF - Fractional calculus and applied analysis : an international journal for theory and applications N2 - The Cattaneo or telegrapher's equation describes the crossover from initial ballistic to normal diffusion. Here we study and survey time-fractional generalisations of this equation that are shown to produce the crossover of the mean squared displacement from superdiffusion to subdiffusion. Conditional solutions are derived in terms of Fox H-functions and the dth-order moments as well as the diffusive flux of the different models are derived. Moreover, the concept of the distribution-like is proposed as an alternative to the probability density function. KW - Cattaneo equation KW - telegrapher's equation KW - crossover dynamics KW - fractional dynamic equations KW - anomalous diffusion KW - superdiffusion and KW - subdiffusion KW - Fox H-functions Y1 - 2020 U6 - https://doi.org/10.1515/fca-2020-0003 SN - 1311-0454 SN - 1314-2224 VL - 23 IS - 1 SP - 55 EP - 102 PB - De Gruyter CY - Berlin ; Boston ER - TY - JOUR A1 - Molina-Garcia, Daniel A1 - Sandev, Trifce A1 - Safdari, Hadiseh A1 - Pagnini, Gianni A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Crossover from anomalous to normal diffusion BT - truncated power-law noise correlations and applications to dynamics in lipid bilayers JF - New Journal of Physics N2 - Abstract The emerging diffusive dynamics in many complex systems show a characteristic crossover behaviour from anomalous to normal diffusion which is otherwise fitted by two independent power-laws. A prominent example for a subdiffusive–diffusive crossover are viscoelastic systems such as lipid bilayer membranes, while superdiffusive–diffusive crossovers occur in systems of actively moving biological cells. We here consider the general dynamics of a stochastic particle driven by so-called tempered fractional Gaussian noise, that is noise with Gaussian amplitude and power-law correlations, which are cut off at some mesoscopic time scale. Concretely we consider such noise with built-in exponential or power-law tempering, driving an overdamped Langevin equation (fractional Brownian motion) and fractional Langevin equation motion. We derive explicit expressions for the mean squared displacement and correlation functions, including different shapes of the crossover behaviour depending on the concrete tempering, and discuss the physical meaning of the tempering. In the case of power-law tempering we also find a crossover behaviour from faster to slower superdiffusion and slower to faster subdiffusion. As a direct application of our model we demonstrate that the obtained dynamics quantitatively describes the subdiffusion–diffusion and subdiffusion–subdiffusion crossover in lipid bilayer systems. We also show that a model of tempered fractional Brownian motion recently proposed by Sabzikar and Meerschaert leads to physically very different behaviour with a seemingly paradoxical ballistic long time scaling. KW - anomalous diffusion KW - truncated power-law correlated noise KW - lipid bilayer membrane dynamics Y1 - 2018 U6 - https://doi.org/10.1088/1367-2630/aae4b2 SN - 1367-2630 VL - 20 PB - IOP Publishing Ltd CY - London und Bad Honnef ER - TY - JOUR A1 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Deformation propagation in responsive polymer network films JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - We study the elastic deformations in a cross-linked polymer network film triggered by the binding of submicron particles with a sticky surface, mimicking the interactions of viral pathogens with thin films of stimulus-responsive polymeric materials such as hydrogels. From extensive Langevin Dynamics simulations we quantify how far the network deformations propagate depending on the elasticity parameters of the network and the adhesion strength of the particles. We examine the dynamics of the collective area shrinkage of the network and obtain some simple relations for the associated characteristic decay lengths. A detailed analysis elucidates how the elastic energy of the network is distributed between stretching and compression modes in response to the particle binding. We also examine the force-distance curves of the repulsion or attraction interactions for a pair of sticky particles in the polymer network film as a function of the particle-particle separation. The results of this computational study provide new insight into collective phenomena in soft polymer network films and may, in particular, be applied to applications for visual detection of pathogens such as viruses via a macroscopic response of thin films of cross-linked hydrogels. (C) 2014 AIP Publishing LLC. Y1 - 2014 U6 - https://doi.org/10.1063/1.4893056 SN - 0021-9606 SN - 1089-7690 VL - 141 IS - 7 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Sandev, Trifce A1 - Chechkin, Aleksei V. A1 - Kantz, Holger A1 - Metzler, Ralf T1 - Diffusion and fokker-planck-smoluchowski equations with generalized memory kernel JF - Fractional calculus and applied analysis : an international journal for theory and applications N2 - We consider anomalous stochastic processes based on the renewal continuous time random walk model with different forms for the probability density of waiting times between individual jumps. In the corresponding continuum limit we derive the generalized diffusion and Fokker-Planck-Smoluchowski equations with the corresponding memory kernels. We calculate the qth order moments in the unbiased and biased cases, and demonstrate that the generalized Einstein relation for the considered dynamics remains valid. The relaxation of modes in the case of an external harmonic potential and the convergence of the mean squared displacement to the thermal plateau are analyzed. KW - continuous time random walk (CTRW) KW - Fokker-Planck-Smoluchowski equation KW - Mittag-Leffler functions KW - anomalous diffusion KW - multi-scaling Y1 - 2015 U6 - https://doi.org/10.1515/fca-2015-0059 SN - 1311-0454 SN - 1314-2224 VL - 18 IS - 4 SP - 1006 EP - 1038 PB - De Gruyter CY - Berlin ER - TY - JOUR A1 - Leijnse, N. A1 - Jeon, J. -H. A1 - Loft, S. A1 - Metzler, Ralf A1 - Oddershede, L. B. T1 - Diffusion inside living human cells JF - European physical journal special topics N2 - Naturally occurring lipid granules diffuse in the cytoplasm and can be used as tracers to map out the viscoelastic landscape inside living cells. Using optical trapping and single particle tracking we found that lipid granules exhibit anomalous diffusion inside human umbilical vein endothelial cells. For these cells the exact diffusional pattern of a particular granule depends on the physiological state of the cell and on the localization of the granule within the cytoplasm. Granules located close to the actin rich periphery of the cell move less than those located towards to the center of the cell or within the nucleus. Also, granules in cells which are stressed by intense laser illumination or which have attached to a surface for a long period of time move in a more restricted fashion than those within healthy cells. For granules diffusing in healthy cells, in regions away from the cell periphery, occurrences of weak ergodicity breaking are observed, similar to the recent observations inside living fission yeast cells [1]. Y1 - 2012 U6 - https://doi.org/10.1140/epjst/e2012-01553-y SN - 1951-6355 VL - 204 IS - 1 SP - 75 EP - 84 PB - Springer CY - Heidelberg ER - TY - JOUR A1 - Kosztolowicz, Tadeusz A1 - Metzler, Ralf T1 - Diffusion of antibiotics through a biofilm in the presence of diffusion and absorption barriers JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We propose a model of antibiotic diffusion through a bacterial biofilm when diffusion and/or absorption barriers develop in the biofilm. The idea of this model is: We deduce details of the diffusion process in a medium in which direct experimental study is difficult, based on probing diffusion in external regions. Since a biofilm has a gel-like consistency, we suppose that subdiffusion of particles in the biofilm may occur. To describe this process we use a fractional subdiffusion-absorption equation with an adjustable anomalous diffusion exponent. The boundary conditions at the boundaries of the biofilm are derived by means of a particle random walk model on a discrete lattice leading to an expression involving a fractional time derivative. We show that the temporal evolution of the total amount of substance that has diffused through the biofilm explicitly depends on whether there is antibiotic absorption in the biofilm. This fact is used to experimentally check for antibiotic absorption in the biofilm and if subdiffusion and absorption parameters of the biofilm change over time. We propose a four-stage model of antibiotic diffusion in biofilm based on the following physical characteristics: whether there is absorption of the antibiotic in the biofilm and whether all biofilm parameters remain unchanged over time. The biological interpretation of the stages, in particular their relation with the bacterial defense mechanisms, is discussed. Theoretical results are compared with empirical results of ciprofloxacin diffusion through Pseudomonas aeruginosa biofilm, and ciprofloxacin and gentamicin diffusion through Proteus mirabilis biofilm. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevE.102.032408 SN - 2470-0045 SN - 2470-0053 VL - 102 IS - 3 PB - American Physical Society CY - Melville, NY 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 - Bauer, Maximilian A1 - Godec, Aljaž A1 - Metzler, Ralf T1 - Diffusion of finite-size particles in two-dimensional channels with random wall configurations JF - Physical chemistry, chemical physics : PCCP ; 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 [J. Chem. Phys., 2012, 137, 024107]. KW - anomalous diffusion KW - fractional dynamics KW - transport KW - nonergodicity KW - coefficient Y1 - 2014 U6 - https://doi.org/10.1039/C3CP55160A SN - 1463-9084 SN - 1463-9076 VL - 16 IS - 13 SP - 6118 EP - 6128 PB - RSC Publications CY - Cambridge 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 - 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 - Sandev, Trifce A1 - Chechkin, Aleksei V. A1 - Korabel, Nickolay A1 - Kantz, Holger A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Distributed-order diffusion equations and multifractality: Models and solutions JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We study distributed-order time fractional diffusion equations characterized by multifractal memory kernels, in contrast to the simple power-law kernel of common time fractional diffusion equations. Based on the physical approach to anomalous diffusion provided by the seminal Scher-Montroll-Weiss continuous time random walk, we analyze both natural and modified-form distributed-order time fractional diffusion equations and compare the two approaches. The mean squared displacement is obtained and its limiting behavior analyzed. We derive the connection between the Wiener process, described by the conventional Langevin equation and the dynamics encoded by the distributed-order time fractional diffusion equation in terms of a generalized subordination of time. A detailed analysis of the multifractal properties of distributed-order diffusion equations is provided. Y1 - 2015 U6 - https://doi.org/10.1103/PhysRevE.92.042117 SN - 1539-3755 SN - 1550-2376 VL - 92 IS - 4 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Distribution of first-reaction times with target regions on boundaries of shell-like domains JF - New Journal of Physics (NJP) N2 - We study the probability density function (PDF) of the first-reaction times between a diffusive ligand and a membrane-bound, immobile imperfect target region in a restricted 'onion-shell' geometry bounded by two nested membranes of arbitrary shapes. For such a setting, encountered in diverse molecular signal transduction pathways or in the narrow escape problem with additional steric constraints, we derive an exact spectral form of the PDF, as well as present its approximate form calculated by help of the so-called self-consistent approximation. For a particular case when the nested domains are concentric spheres, we get a fully explicit form of the approximated PDF, assess the accuracy of this approximation, and discuss various facets of the obtained distributions. Our results can be straightforwardly applied to describe the PDF of the terminal reaction event in multi-stage signal transduction processes. KW - diffusion KW - first-passage time KW - first-reaction time KW - shell-like geometries KW - approximate methods Y1 - 2021 U6 - https://doi.org/10.1088/1367-2630/ac4282 SN - 1367-2630 VL - 2021 SP - 1 EP - 23 PB - IOP Publishing CY - London ET - 23 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 - Mejia-Monasterio, Carlos A1 - Metzler, Ralf A1 - Vollmer, Jürgen T1 - Editorial: anomalous transport BT - applications, mathematical perspectives, and big data JF - Frontiers in Physics KW - anomalous (or non-Fickian) diffusion KW - anomalous heat conduction KW - stochastic dynamics KW - molecular overcrowding KW - dynamical systems Y1 - 2020 U6 - https://doi.org/10.3389/fphy.2020.622417 SN - 2296-424X VL - 8 PB - Frontiers Media CY - Lausanne ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains JF - New journal of physics : the open-access journal for physics N2 - We study the mean first passage time (MFPT) to a reaction event on a specific site in a cylindrical geometry-characteristic, for instance, for bacterial cells, with a concentric inner cylinder representing the nuclear region of the bacterial cell. A similar problem emerges in the description of a diffusive search by a transcription factor protein for a specific binding region on a single strand of DNA. We develop a unified theoretical approach to study the underlying boundary value problem which is based on a self-consistent approximation of the mixed boundary condition. Our approach permits us to derive explicit, novel, closed-form expressions for the MFPT valid for a generic setting with an arbitrary relation between the system parameters. We analyse this general result in the asymptotic limits appropriate for the above-mentioned biophysical problems. Our investigation reveals the crucial role of the target aspect ratio and of the intrinsic reactivity of the binding region, which were disregarded in previous studies. Theoretical predictions are confirmed by numerical simulations. KW - first passage time KW - cylindrical geometry KW - aspect ratio KW - protein search Y1 - 2017 U6 - https://doi.org/10.1088/1367-2630/aa8ed9 SN - 1367-2630 VL - 19 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains JF - New journal of physics N2 - Westudy the mean first passage time (MFPT) to a reaction event on a specific site in a cylindrical geometry—characteristic, for instance, for bacterial cells, with a concentric inner cylinder representing the nuclear region of the bacterial cell. Asimilar problem emerges in the description of a diffusive search by a transcription factor protein for a specific binding region on a single strand of DNA.We develop a unified theoretical approach to study the underlying boundary value problem which is based on a self-consistent approximation of the mixed boundary condition. Our approach permits us to derive explicit, novel, closed-form expressions for the MFPT valid for a generic setting with an arbitrary relation between the system parameters.Weanalyse this general result in the asymptotic limits appropriate for the above-mentioned biophysical problems. Our investigation reveals the crucial role of the target aspect ratio and of the intrinsic reactivity of the binding region, which were disregarded in previous studies. Theoretical predictions are confirmed by numerical simulations. KW - first passage time KW - cylindrical geometry KW - aspect ratio KW - protein search Y1 - 2017 U6 - https://doi.org/10.1088/1367-2630/aa8ed9 SN - 1367-2630 VL - 19 SP - 1 EP - 11 PB - IOP CY - London ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Ergodicity breaking and particle spreading in noisy heterogeneous diffusion processes JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - We study noisy heterogeneous diffusion processes with a position dependent diffusivity of the form D(x) similar to D-0 vertical bar x vertical bar (alpha 0) in the presence of annealed and quenched disorder of the environment, corresponding to an effective variation of the exponent a in time and space. In the case of annealed disorder, for which effectively alpha(0) = alpha(0)(t), we show how the long time scaling of the ensemble mean squared displacement (MSD) and the amplitude variation of individual realizations of the time averaged MSD are affected by the disorder strength. For the case of quenched disorder, the long time behavior becomes effectively Brownian after a number of jumps between the domains of a stratified medium. In the latter situation, the averages are taken over both an ensemble of particles and different realizations of the disorder. As physical observables, we analyze in detail the ensemble and time averaged MSDs, the ergodicity breaking parameter, and higher order moments of the time averages. (C) 2015 AIP Publishing LLC. Y1 - 2015 U6 - https://doi.org/10.1063/1.4917077 SN - 0021-9606 SN - 1089-7690 VL - 142 IS - 14 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Ergodicity breaking, ageing, and confinement in generalized diffusion processes with position and time dependent diffusivity JF - Journal of statistical mechanics: theory and experiment N2 - We study generalized anomalous diffusion processes whose diffusion coefficient D(x, t) similar to D-0x(alpha)t(beta) depends on both the position x of the test particle and the process time t. This process thus combines the features of scaled Brownian motion and heterogeneous diffusion parent processes. We compute the ensemble and time averaged mean squared displacements of this generalized diffusion process. The scaling exponent of the ensemble averaged mean squared displacement is shown to be the product of the critical exponents of the parent processes, and describes both subdiffusive and superdiffusive systems. We quantify the amplitude fluctuations of the time averaged mean squared displacement as function of the length of the time series and the lag time. In particular, we observe a weak ergodicity breaking of this generalized diffusion process: even in the long time limit the ensemble and time averaged mean squared displacements are strictly disparate. When we start to observe this process some time after its initiation we observe distinct features of ageing. We derive a universal ageing factor for the time averaged mean squared displacement containing all information on the ageing time and the measurement time. External confinement is shown to alter the magnitudes and statistics of the ensemble and time averaged mean squared displacements. KW - diffusion Y1 - 2015 U6 - https://doi.org/10.1088/1742-5468/2015/05/P05010 SN - 1742-5468 PB - IOP Publ. Ltd. CY - Bristol 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 - Grebenkov, Denis S. A1 - Sposini, Vittoria A1 - Metzler, Ralf A1 - Oshanin, Gleb A1 - Seno, Flavio T1 - Exact distributions of the maximum and range of random diffusivity processes JF - New Journal of Physics N2 - We study the extremal properties of a stochastic process xt defined by the Langevin equation ẋₜ =√2Dₜ ξₜ, in which ξt is a Gaussian white noise with zero mean and Dₜ is a stochastic‘diffusivity’, defined as a functional of independent Brownian motion Bₜ.We focus on threechoices for the random diffusivity Dₜ: cut-off Brownian motion, Dₜt ∼ Θ(Bₜ), where Θ(x) is the Heaviside step function; geometric Brownian motion, Dₜ ∼ exp(−Bₜ); and a superdiffusive process based on squared Brownian motion, Dₜ ∼ B²ₜ. For these cases we derive exact expressions for the probability density functions of the maximal positive displacement and of the range of the process xₜ on the time interval ₜ ∈ (0, T).We discuss the asymptotic behaviours of the associated probability density functions, compare these against the behaviour of the corresponding properties of standard Brownian motion with constant diffusivity (Dₜ = D0) and also analyse the typical behaviour of the probability density functions which is observed for a majority of realisations of the stochastic diffusivity process. KW - random diffusivity KW - extremal values KW - maximum and range KW - diffusion KW - Brownian motion Y1 - 2021 U6 - https://doi.org/10.1088/1367-2630/abd313 SN - 1367-2630 VL - 23 PB - Dt. Physikalische Ges. CY - Bad Honnef ER - TY - JOUR A1 - Scott, Shane A1 - Weiss, Matthias A1 - Selhuber-Unkel, Christine A1 - Barooji, Younes F. A1 - Sabri, Adal A1 - Erler, Janine T. A1 - Metzler, Ralf A1 - Oddershede, Lene B. T1 - Extracting, quantifying, and comparing dynamical and biomechanical properties of living matter through single particle tracking JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - A panoply of new tools for tracking single particles and molecules has led to an explosion of experimental data, leading to novel insights into physical properties of living matter governing cellular development and function, health and disease. In this Perspective, we present tools to investigate the dynamics and mechanics of living systems from the molecular to cellular scale via single-particle techniques. In particular, we focus on methods to measure, interpret, and analyse complex data sets that are associated with forces, materials properties, transport, and emergent organisation phenomena within biological and soft-matter systems. Current approaches, challenges, and existing solutions in the associated fields are outlined in order to support the growing community of researchers at the interface of physics and the life sciences. Each section focuses not only on the general physical principles and the potential for understanding living matter, but also on details of practical data extraction and analysis, discussing limitations, interpretation, and comparison across different experimental realisations and theoretical frameworks. Particularly relevant results are introduced as examples. While this Perspective describes living matter from a physical perspective, highlighting experimental and theoretical physics techniques relevant for such systems, it is also meant to serve as a solid starting point for researchers in the life sciences interested in the implementation of biophysical methods. Y1 - 2022 U6 - https://doi.org/10.1039/d2cp01384c SN - 1463-9076 SN - 1463-9084 VL - 25 IS - 3 SP - 1513 EP - 1537 PB - RSC Publ. 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 - Shin, Jaeoh A1 - Cherstvy, Andrey G. A1 - Kim, Won Kyu A1 - Metzler, Ralf T1 - Facilitation of polymer looping and giant polymer diffusivity in crowded solutions of active particles JF - New journal of physics : the open-access journal for physics N2 - We study the dynamics of polymer chains in a bath of self-propelled particles (SPP) by extensive Langevin dynamics simulations in a two-dimensional model system. Specifically, we analyse the polymer looping properties versus the SPP activity and investigate how the presence of the active particles alters the chain conformational statistics. We find that SPPs tend to extend flexible polymer chains, while they rather compactify stiffer semiflexible polymers, in agreement with previous results. Here we show that higher activities of SPPs yield a higher effective temperature of the bath and thus facilitate the looping kinetics of a passive polymer chain. We explicitly compute the looping probability and looping time in a wide range of the model parameters. We also analyse the motion of a monomeric tracer particle and the polymer's centre of mass in the presence of the active particles in terms of the time averaged mean squared displacement, revealing a giant diffusivity enhancement for the polymer chain via SPP pooling. Our results are applicable to rationalising the dimensions and looping kinetics of biopolymers at constantly fluctuating and often actively driven conditions inside biological cells or in suspensions of active colloidal particles or bacteria cells. KW - diffusion KW - active transport KW - polymers Y1 - 2015 U6 - https://doi.org/10.1088/1367-2630/17/11/113008 SN - 1367-2630 VL - 17 PB - IOP Publ. Ltd. CY - Bristol 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 - 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 - 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 JF - New Journal of 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 - Lévy flights KW - Lévy 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 PB - Dt. Physikalische Ges. CY - Bad Honnef ER - TY - JOUR A1 - Mutothya, Nicholas Mwilu A1 - Xu, Yong A1 - Li, Yongge A1 - Metzler, Ralf A1 - Mutua, Nicholas Muthama T1 - First passage dynamics of stochastic motion in heterogeneous media driven by correlated white Gaussian and coloured non-Gaussian noises JF - Journal of physics. Complexity N2 - We study the first passage dynamics for a diffusing particle experiencing a spatially varying diffusion coefficient while driven by correlated additive Gaussian white noise and multiplicative coloured non-Gaussian noise. We consider three functional forms for position dependence of the diffusion coefficient: power-law, exponential, and logarithmic. The coloured non-Gaussian noise is distributed according to Tsallis' q-distribution. Tracks of the non-Markovian systems are numerically simulated by using the fourth-order Runge-Kutta algorithm and the first passage times (FPTs) are recorded. The FPT density is determined along with the mean FPT (MFPT). Effects of the noise intensity and self-correlation of the multiplicative noise, the intensity of the additive noise, the cross-correlation strength, and the non-extensivity parameter on the MFPT are discussed. KW - first passage KW - diffusion KW - non-Gaussian KW - correlated noise Y1 - 2021 U6 - https://doi.org/10.1088/2632-072X/ac35b5 SN - 2632-072X VL - 2 PB - IOP Publishing CY - Bristol 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 - Godec, Aljaz A1 - Metzler, Ralf T1 - First passage time distribution in heterogeneity controlled kinetics: going beyond the mean first passage time JF - Scientific reports N2 - The first passage is a generic concept for quantifying when a random quantity such as the position of a diffusing molecule or the value of a stock crosses a preset threshold (target) for the first time. The last decade saw an enlightening series of new results focusing mostly on the so-called mean and global first passage time (MFPT and GFPT, respectively) of such processes. Here we push the understanding of first passage processes one step further. For a simple heterogeneous system we derive rigorously the complete distribution of first passage times (FPTs). Our results demonstrate that the typical FPT significantly differs from the MFPT, which corresponds to the long time behaviour of the FPT distribution. Conversely, the short time behaviour is shown to correspond to trajectories connecting directly from the initial value to the target. Remarkably, we reveal a previously overlooked third characteristic time scale of the first passage dynamics mirroring brief excursion away from the target. Y1 - 2016 U6 - https://doi.org/10.1038/srep20349 SN - 2045-2322 VL - 6 PB - Nature Publ. Group CY - London ER - TY - JOUR A1 - Abdolvahab, Rouhollah Haji A1 - Metzler, Ralf A1 - Ejtehadi, Mohammad Reza T1 - First passage time distribution of chaperone driven polymer translocation through a nanopore homopolymer and heteropolymer cases JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - Combining the advection-diffusion equation approach with Monte Carlo simulations we study chaperone driven polymer translocation of a stiff polymer through a nanopore. We demonstrate that the probability density function of first passage times across the pore depends solely on the Peclet number, a dimensionless parameter comparing drift strength and diffusivity. Moreover it is shown that the characteristic exponent in the power-law dependence of the translocation time on the chain length, a function of the chaperone-polymer binding energy, the chaperone concentration, and the chain length, is also effectively determined by the Peclet number. We investigate the effect of the chaperone size on the translocation process. In particular, for large chaperone size, the translocation progress and the mean waiting time as function of the reaction coordinate exhibit pronounced sawtooth-shapes. The effects of a heterogeneous polymer sequence on the translocation dynamics is studied in terms of the translocation velocity, the probability distribution for the translocation progress, and the monomer waiting times. (C) 2011 American Institute of Physics. Y1 - 2011 U6 - https://doi.org/10.1063/1.3669427 SN - 0021-9606 VL - 135 IS - 24 PB - American Institute of Physics CY - Melville 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 - Mattos, Thiago G. A1 - Mejia-Monasterio, Carlos A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - First passages in bounded domains When is the mean first passage time meaningful? JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We study the first passage statistics to adsorbing boundaries of a Brownian motion in bounded two-dimensional domains of different shapes and configurations of the adsorbing and reflecting boundaries. From extensive numerical analysis we obtain the probability P(omega) distribution of the random variable omega = tau(1)/(tau(1) + tau(2)), which is a measure for how similar the first passage times tau(1) and tau(2) are of two independent realizations of a Brownian walk starting at the same location. We construct a chart for each domain, determining whether P(omega) represents a unimodal, bell-shaped form, or a bimodal, M-shaped behavior. While in the former case the mean first passage time (MFPT) is a valid characteristic of the first passage behavior, in the latter case it is an insufficient measure for the process. Strikingly we find a distinct turnover between the two modes of P(omega), characteristic for the domain shape and the respective location of absorbing and reflective boundaries. Our results demonstrate that large fluctuations of the first passage times may occur frequently in two-dimensional domains, rendering quite vague the general use of the MFPT as a robust measure of the actual behavior even in bounded domains, in which all moments of the first passage distribution exist. Y1 - 2012 U6 - https://doi.org/10.1103/PhysRevE.86.031143 SN - 1539-3755 VL - 86 IS - 3 PB - American Physical Society CY - College Park 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 -