@article{AwadMetzler2020, author = {Awad, Emad and Metzler, Ralf}, title = {Crossover dynamics from superdiffusion to subdiffusion}, series = {Fractional calculus and applied analysis : an international journal for theory and applications}, volume = {23}, journal = {Fractional calculus and applied analysis : an international journal for theory and applications}, number = {1}, publisher = {De Gruyter}, address = {Berlin ; Boston}, issn = {1311-0454}, doi = {10.1515/fca-2020-0003}, pages = {55 -- 102}, year = {2020}, abstract = {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.}, language = {en} } @article{BauerGodecMetzler2014, author = {Bauer, Maximilian and Godec, Aljaž and Metzler, Ralf}, title = {Diffusion of finite-size particles in two-dimensional channels with random wall configurations}, series = {Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies}, volume = {16}, journal = {Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies}, number = {13}, publisher = {RSC Publications}, address = {Cambridge}, issn = {1463-9084}, doi = {10.1039/C3CP55160A}, pages = {6118 -- 6128}, year = {2014}, abstract = {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].}, language = {en} } @misc{BauerGodecMetzler2014, author = {Bauer, Maximilian and Godec, Aljaž and Metzler, Ralf}, title = {Diffusion of finite-size particles in two-dimensional channels with random wall configurations}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-76199}, year = {2014}, abstract = {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].}, language = {en} } @article{BodrovaChechkinCherstvyetal.2015, author = {Bodrova, Anna S. and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Ultraslow scaled Brownian motion}, series = {New journal of physics : the open-access journal for physics}, volume = {17}, journal = {New journal of physics : the open-access journal for physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/17/6/063038}, pages = {16}, year = {2015}, abstract = {We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form D(t) similar or equal to 1/t. For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In a harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations.}, language = {en} } @article{ChechkinZaidLomholtetal.2013, author = {Chechkin, Aleksei V. and Zaid, I. M. and Lomholt, M. A. and Sokolov, Igor M. and Metzler, Ralf}, title = {Bulk-mediated surface diffusion on a cylinder in the fast exchange limit}, series = {Mathematical modelling of natural phenomena}, volume = {8}, journal = {Mathematical modelling of natural phenomena}, number = {2}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0973-5348}, doi = {10.1051/mmnp/20138208}, pages = {114 -- 126}, year = {2013}, abstract = {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.}, language = {en} } @misc{ChechkinZaidLomholtetal.2013, author = {Chechkin, Aleksei V. and Zaid, Irwin M. and Lomholt, Michael A. and Sokolov, Igor M. and Metzler, Ralf}, title = {Bulk-mediated surface diffusion on a cylinder in the fast exchange limit}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, number = {593}, issn = {1866-8372}, doi = {10.25932/publishup-41548}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-415480}, pages = {114 -- 126}, year = {2013}, abstract = {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.}, language = {en} } @article{CherstvyChechkinMetzler2014, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity}, series = {Soft matter}, volume = {2014}, journal = {Soft matter}, number = {10}, publisher = {Royal Society of Chemistry}, issn = {2046-2069}, doi = {10.1039/c3sm52846d}, pages = {1591 -- 1601}, year = {2014}, abstract = {We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells.}, language = {en} } @misc{CherstvyChechkinMetzler2014, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity}, number = {168}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-74021}, pages = {1591 -- 1601}, year = {2014}, abstract = {We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells.}, language = {en} } @article{CherstvyChechkinMetzler2014, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Ageing and confinement in non-ergodic heterogeneous diffusion processes}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {47}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {48}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/47/48/485002}, pages = {18}, year = {2014}, abstract = {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.}, language = {en} } @misc{CherstvyMetzler2013, author = {Cherstvy, Andrey G. and Metzler, Ralf}, title = {Population splitting, trapping, and non-ergodicity in heterogeneous diffusion processes}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-94468}, pages = {20220 -- 20235}, year = {2013}, abstract = {We consider diffusion processes with a spatially varying diffusivity giving rise to anomalous diffusion. Such heterogeneous diffusion processes are analysed for the cases of exponential, power-law, and logarithmic dependencies of the diffusion coefficient on the particle position. Combining analytical approaches with stochastic simulations, we show that the functional form of the space-dependent diffusion coefficient and the initial conditions of the diffusing particles are vital for their statistical and ergodic properties. In all three cases a weak ergodicity breaking between the time and ensemble averaged mean squared displacements is observed. We also demonstrate a population splitting of the time averaged traces into fast and slow diffusers for the case of exponential variation of the diffusivity as well as a particle trapping in the case of the logarithmic diffusivity. Our analysis is complemented by the quantitative study of the space coverage, the diffusive spreading of the probability density, as well as the survival probability.}, language = {en} } @article{CherstvySafdariMetzler2021, author = {Cherstvy, Andrey G. and Safdari, Hadiseh and Metzler, Ralf}, title = {Anomalous diffusion, nonergodicity, and ageing for exponentially and logarithmically time-dependent diffusivity}, series = {Journal of physics. D, Applied physics}, volume = {54}, journal = {Journal of physics. D, Applied physics}, number = {19}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0022-3727}, doi = {10.1088/1361-6463/abdff0}, pages = {18}, year = {2021}, abstract = {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.}, language = {en} } @article{DieterichKlagesChechkin2015, author = {Dieterich, Peter and Klages, Rainer and Chechkin, Aleksei V.}, title = {Fluctuation relations for anomalous dynamics generated by time-fractional Fokker-Planck equations}, series = {New journal of physics : the open-access journal for physics}, volume = {17}, journal = {New journal of physics : the open-access journal for physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/17/7/075004}, pages = {14}, year = {2015}, abstract = {Anomalous dynamics characterized by non-Gaussian probability distributions (PDFs) and/or temporal long-range correlations can cause subtle modifications of conventional fluctuation relations (FRs). As prototypes we study three variants of a generic time-fractional Fokker-Planck equation with constant force. Type A generates superdiffusion, type B subdiffusion and type C both super-and subdiffusion depending on parameter variation. Furthermore type C obeys a fluctuation-dissipation relation whereas A and B do not. We calculate analytically the position PDFs for all three cases and explore numerically their strongly non-Gaussian shapes. While for type C we obtain the conventional transient work FR, type A and type B both yield deviations by featuring a coefficient that depends on time and by a nonlinear dependence on the work. We discuss possible applications of these types of dynamics and FRs to experiments.}, language = {en} } @phdthesis{Doerries2024, author = {D{\"o}rries, Timo Julian}, title = {Anomalous transport and non-Gaussian dynamics in mobile-immobile models}, doi = {10.25932/publishup-63495}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-634959}, school = {Universit{\"a}t Potsdam}, pages = {ii, 177}, year = {2024}, abstract = {The mobile-immobile model (MIM) has been established in geoscience in the context of contaminant transport in groundwater. Here the tracer particles effectively immobilise, e.g., due to diffusion into dead-end pores or sorption. The main idea of the MIM is to split the total particle density into a mobile and an immobile density. Individual tracers switch between the mobile and immobile state following a two-state telegraph process, i.e., the residence times in each state are distributed exponentially. In geoscience the focus lies on the breakthrough curve (BTC), which is the concentration at a fixed location over time. We apply the MIM to biological experiments with a special focus on anomalous scaling regimes of the mean squared displacement (MSD) and non-Gaussian displacement distributions. As an exemplary system, we have analysed the motion of tau proteins, that diffuse freely inside axons of neurons. Their free diffusion thereby corresponds to the mobile state of the MIM. Tau proteins stochastically bind to microtubules, which effectively immobilises the tau proteins until they unbind and continue diffusing. Long immobilisation durations compared to the mobile durations give rise to distinct non-Gaussian Laplace shaped distributions. It is accompanied by a plateau in the MSD for initially mobile tracer particles at relevant intermediate timescales. An equilibrium fraction of initially mobile tracers gives rise to non-Gaussian displacements at intermediate timescales, while the MSD remains linear at all times. In another setting bio molecules diffuse in a biosensor and transiently bind to specific receptors, where advection becomes relevant in the mobile state. The plateau in the MSD observed for the advection-free setting and long immobilisation durations persists also for the case with advection. We find a new clear regime of anomalous diffusion with non-Gaussian distributions and a cubic scaling of the MSD. This regime emerges for initially mobile and for initially immobile tracers. For an equilibrium fraction of initially mobile tracers we observe an intermittent ballistic scaling of the MSD. The long-time effective diffusion coefficient is enhanced by advection, which we physically explain with the variance of mobile durations. Finally, we generalize the MIM to incorporate arbitrary immobilisation time distributions and focus on a Mittag-Leffler immobilisation time distribution with power-law tail ~ t^(-1-mu) with 0 similar or equal to log(gamma)t. Comparison of annealed (renewal) continuous time random walks (CTRWs) with logarithmic waiting time distribution psi(tau) similar or equal to 1/(tau log(1+gamma)tau) and Sinai diffusion in quenched random landscapes reveals striking similarities, despite the great differences in their physical nature. In particular, they exhibit a weakly non-ergodic disparity of the time-averaged and ensemble-averaged MSDs. Remarkably, for the CTRW we observe that the fluctuations of time averages become universal, with an exponential suppression of mobile trajectories. We discuss the fundamental connection between the Golosov localization effect and non-ergodicity in the sense of the disparity between ensemble-averaged MSD and time-averaged MSD.}, language = {en} } @article{GuggenbergerChechkinMetzler2021, author = {Guggenberger, Tobias and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Fractional Brownian motion in superharmonic potentials and non-Boltzmann stationary distributions}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {54}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {29}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/ac019b}, pages = {17}, year = {2021}, abstract = {We study the stochastic motion of particles driven by long-range correlated fractional Gaussian noise (FGN) in a superharmonic external potential of the form U(x) proportional to x(2n) (n is an element of N). When the noise is considered to be external, the resulting overdamped motion is described by the non-Markovian Langevin equation for fractional Brownian motion. For this case we show the existence of long time, stationary probability density functions (PDFs) the shape of which strongly deviates from the naively expected Boltzmann PDF in the confining potential U(x). We analyse in detail the temporal approach to stationarity as well as the shape of the non-Boltzmann stationary PDF. A typical characteristic is that subdiffusive, antipersistent (with negative autocorrelation) motion tends to effect an accumulation of probability close to the origin as compared to the corresponding Boltzmann distribution while the opposite trend occurs for superdiffusive (persistent) motion. For this latter case this leads to distinct bimodal shapes of the PDF. This property is compared to a similar phenomenon observed for Markovian Levy flights in superharmonic potentials. We also demonstrate that the motion encoded in the fractional Langevin equation driven by FGN always relaxes to the Boltzmann distribution, as in this case the fluctuation-dissipation theorem is fulfilled.}, language = {en} } @misc{GuggenbergerPagniniVojtaetal.2019, author = {Guggenberger, Tobias and Pagnini, Gianni and Vojta, Thomas and Metzler, Ralf}, title = {Fractional Brownian motion in a finite interval}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {755}, issn = {1866-8372}, doi = {10.25932/publishup-43666}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-436665}, pages = {13}, year = {2019}, abstract = {Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) 〈X² (t)〉 ⋍ tᵅ with 1 < α < 2, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion (0 < α < 1) this behaviour is reversed and the particle density is depleted close to the boundaries. The MSD in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers.}, language = {en} } @article{GuggenbergerPagniniVojtaetal.2019, author = {Guggenberger, Tobias and Pagnini, Gianni and Vojta, Thomas and Metzler, Ralf}, title = {Fractional Brownian motion in a finite interval}, series = {New Journal of Physics}, volume = {21}, journal = {New Journal of Physics}, publisher = {Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics}, address = {Bad Honnef und London}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab075f}, pages = {13}, year = {2019}, abstract = {Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) 〈X² (t)〉 ⋍ tᵅ with 1 < α < 2, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion (0 < α < 1) this behaviour is reversed and the particle density is depleted close to the boundaries. The MSD in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers.}, language = {en} } @article{KruesemannGodecMetzler2015, author = {Kr{\"u}semann, Henning and Godec, Aljaz and Metzler, Ralf}, title = {Ageing first passage time density in continuous time random walks and quenched energy landscapes}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {48}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {28}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/48/28/285001}, pages = {20}, year = {2015}, abstract = {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.}, language = {en} } @misc{Metzler2020, author = {Metzler, Ralf}, title = {Superstatistics and non-Gaussian diffusion}, series = {The European physical journal special topics}, volume = {229}, journal = {The European physical journal special topics}, number = {5}, publisher = {Springer}, address = {Heidelberg}, issn = {1951-6355}, doi = {10.1140/epjst/e2020-900210-x}, pages = {711 -- 728}, year = {2020}, abstract = {Brownian motion and viscoelastic anomalous diffusion in homogeneous environments are intrinsically Gaussian processes. In a growing number of systems, however, non-Gaussian displacement distributions of these processes are being reported. The physical cause of the non-Gaussianity is typically seen in different forms of disorder. These include, for instance, imperfect "ensembles" of tracer particles, the presence of local variations of the tracer mobility in heteroegenous environments, or cases in which the speed or persistence of moving nematodes or cells are distributed. From a theoretical point of view stochastic descriptions based on distributed ("superstatistical") transport coefficients as well as time-dependent generalisations based on stochastic transport parameters with built-in finite correlation time are invoked. After a brief review of the history of Brownian motion and the famed Gaussian displacement distribution, we here provide a brief introduction to the phenomenon of non-Gaussianity and the stochastic modelling in terms of superstatistical and diffusing-diffusivity approaches.}, language = {en} } @article{MetzlerCherstvyChechkinetal.2015, author = {Metzler, Ralf and Cherstvy, Andrey G. and Chechkin, Aleksei V. and Bodrova, Anna S.}, title = {Ultraslow scaled Brownian motion}, series = {New journal of physics : the open-access journal for physics}, volume = {17}, journal = {New journal of physics : the open-access journal for physics}, number = {063038}, publisher = {Dt. Physikalische Ges., IOP}, address = {Bad Honnef, London}, issn = {1367-2630}, doi = {10.1088/1367-2630/17/6/063038}, year = {2015}, abstract = {We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form . For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In a harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations.}, language = {en} } @misc{MetzlerCherstvyChechkinetal.2015, author = {Metzler, Ralf and Cherstvy, Andrey G. and Chechkin, Aleksei V. and Bodrova, Anna S.}, title = {Ultraslow scaled Brownian motion}, series = {New journal of physics : the open-access journal for physics}, journal = {New journal of physics : the open-access journal for physics}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-78618}, year = {2015}, abstract = {We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form . For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In a harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations.}, language = {en} } @article{MolinaGarciaSandevSafdarietal.2018, author = {Molina-Garcia, Daniel and Sandev, Trifce and Safdari, Hadiseh and Pagnini, Gianni and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Crossover from anomalous to normal diffusion}, series = {New Journal of Physics}, volume = {20}, journal = {New Journal of Physics}, publisher = {IOP Publishing Ltd}, address = {London und Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/aae4b2}, pages = {28}, year = {2018}, abstract = {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.}, language = {en} } @misc{MolinaGarciaSandevSafdarietal.2019, author = {Molina-Garcia, Daniel and Sandev, Trifce and Safdari, Hadiseh and Pagnini, Gianni and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Crossover from anomalous to normal diffusion}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {507}, issn = {1866-8372}, doi = {10.25932/publishup-42259}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-422590}, pages = {28}, year = {2019}, abstract = {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.}, language = {en} } @article{MutothyaXuLietal.2021, author = {Mutothya, Nicholas Mwilu and Xu, Yong and Li, Yongge and Metzler, Ralf}, title = {Characterising stochastic motion in heterogeneous media driven by coloured non-Gaussian noise}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {54}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {29}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/abfba6}, pages = {31}, year = {2021}, abstract = {We study the stochastic motion of a test particle in a heterogeneous medium in terms of a position dependent diffusion coefficient mimicking measured deterministic diffusivity gradients in biological cells or the inherent heterogeneity of geophysical systems. Compared to previous studies we here investigate the effect of the interplay of anomalous diffusion effected by position dependent diffusion coefficients and coloured non-Gaussian noise. The latter is chosen to be distributed according to Tsallis' q-distribution, representing a popular example for a non-extensive statistic. We obtain the ensemble and time averaged mean squared displacements for this generalised process and establish its non-ergodic properties as well as analyse the non-Gaussian nature of the associated displacement distribution. We consider both non-stratified and stratified environments.}, language = {en} } @article{PalyulinAlaNissilaMetzler2014, author = {Palyulin, Vladimir V. and Ala-Nissila, Tapio and Metzler, Ralf}, title = {Polymer translocation: the first two decades and the recent diversification}, series = {Soft matter}, volume = {45}, journal = {Soft matter}, number = {10}, editor = {Metzler, Ralf}, publisher = {the Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-76266}, pages = {9016 -- 9037}, year = {2014}, abstract = {Probably no other field of statistical physics at the borderline of soft matter and biological physics has caused such a flurry of papers as polymer translocation since the 1994 landmark paper by Bezrukov, Vodyanoy, and Parsegian and the study of Kasianowicz in 1996. Experiments, simulations, and theoretical approaches are still contributing novel insights to date, while no universal consensus on the statistical understanding of polymer translocation has been reached. We here collect the published results, in particular, the famous-infamous debate on the scaling exponents governing the translocation process. We put these results into perspective and discuss where the field is going. In particular, we argue that the phenomenon of polymer translocation is non-universal and highly sensitive to the exact specifications of the models and experiments used towards its analysis.}, language = {en} } @misc{PalyulinAlaNissilaMetzler2014, author = {Palyulin, Vladimir V. and Ala-Nissila, Tapio and Metzler, Ralf}, title = {Polymer translocation: the first two decades and the recent diversification}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-76287}, pages = {9016 -- 9037}, year = {2014}, abstract = {Probably no other field of statistical physics at the borderline of soft matter and biological physics has caused such a flurry of papers as polymer translocation since the 1994 landmark paper by Bezrukov, Vodyanoy, and Parsegian and the study of Kasianowicz in 1996. Experiments, simulations, and theoretical approaches are still contributing novel insights to date, while no universal consensus on the statistical understanding of polymer translocation has been reached. We here collect the published results, in particular, the famous-infamous debate on the scaling exponents governing the translocation process. We put these results into perspective and discuss where the field is going. In particular, we argue that the phenomenon of polymer translocation is non-universal and highly sensitive to the exact specifications of the models and experiments used towards its analysis.}, language = {en} } @article{PalyulinMetzler2014, author = {Palyulin, Vladimir V. and Metzler, Ralf}, title = {Speeding up the first-passage for subdiffusion by introducing a finite potential barrier}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {47}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {3}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/47/3/032002}, pages = {13}, year = {2014}, abstract = {We show that for a subdiffusive continuous time random walk with scale-free waiting time distribution the first-passage dynamics on a finite interval can be optimized by introduction of a piecewise linear potential barrier. Analytical results for the survival probability and first-passage density based on the fractional Fokker-Planck equation are shown to agree well with Monte Carlo simulations results. As an application we discuss an improved design for efficient translocation of gradient copolymers compared to homopolymer translocation in a quasi-equilibrium approximation.}, language = {en} } @article{SafdariCherstvyChechkinetal.2015, author = {Safdari, Hadiseh and Cherstvy, Andrey G. and Chechkin, Aleksei V. and Thiel, Felix and Sokolov, Igor M. and Metzler, Ralf}, title = {Quantifying the non-ergodicity of scaled Brownian motion}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {48}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {37}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/48/37/375002}, pages = {18}, year = {2015}, abstract = {We examine the non-ergodic properties of scaled Brownian motion (SBM), a non-stationary stochastic process with a time dependent diffusivity of the form D(t) similar or equal to t(alpha-1). We compute the ergodicity breaking parameter EB in the entire range of scaling exponents a, both analytically and via extensive computer simulations of the stochastic Langevin equation. We demonstrate that in the limit of long trajectory lengths T and short lag times Delta the EB parameter as function of the scaling exponent a has no divergence at alpha - 1/2 and present the asymptotes for EB in different limits. We generalize the analytical and simulations results for the time averaged and ergodic properties of SBM in the presence of ageing, that is, when the observation of the system starts only a finite time span after its initiation. The approach developed here for the calculation of the higher time averaged moments of the particle displacement can be applied to derive the ergodic properties of other stochastic processes such as fractional Brownian motion.}, language = {en} } @article{SandevChechkinKantzetal.2015, author = {Sandev, Trifce and Chechkin, Aleksei V. and Kantz, Holger and Metzler, Ralf}, title = {Diffusion and fokker-planck-smoluchowski equations with generalized memory kernel}, series = {Fractional calculus and applied analysis : an international journal for theory and applications}, volume = {18}, journal = {Fractional calculus and applied analysis : an international journal for theory and applications}, number = {4}, publisher = {De Gruyter}, address = {Berlin}, issn = {1311-0454}, doi = {10.1515/fca-2015-0059}, pages = {1006 -- 1038}, year = {2015}, abstract = {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.}, language = {en} } @article{SandevIominKantzetal.2016, author = {Sandev, Trifce and Iomin, Alexander and Kantz, Holger and Metzler, Ralf and Chechkin, Aleksei V.}, title = {Comb Model with Slow and Ultraslow Diffusion}, series = {Mathematical modelling of natural phenomena}, volume = {11}, journal = {Mathematical modelling of natural phenomena}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0973-5348}, doi = {10.1051/mmnp/201611302}, pages = {18 -- 33}, year = {2016}, abstract = {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.}, language = {en} } @article{SandevMetzlerChechkin2018, author = {Sandev, Trifce and Metzler, Ralf and Chechkin, Aleksei V.}, title = {From continuous time random walks to the generalized diffusion equation}, series = {Fractional calculus and applied analysis : an international journal for theory and applications}, volume = {21}, journal = {Fractional calculus and applied analysis : an international journal for theory and applications}, number = {1}, publisher = {De Gruyter}, address = {Berlin}, issn = {1311-0454}, doi = {10.1515/fca-2018-0002}, pages = {10 -- 28}, year = {2018}, abstract = {We obtain a generalized diffusion equation in modified or Riemann-Liouville form from continuous time random walk theory. The waiting time probability density function and mean squared displacement for different forms of the equation are explicitly calculated. We show examples of generalized diffusion equations in normal or Caputo form that encode the same probability distribution functions as those obtained from the generalized diffusion equation in modified form. The obtained equations are general and many known fractional diffusion equations are included as special cases.}, language = {en} } @article{SandevMetzlerTomovski2012, author = {Sandev, Trifce and Metzler, Ralf and Tomovski, Zivorad}, title = {Velocity and displacement correlation functions for fractional generalized Langevin equations}, series = {Fractional calculus and applied analysis : an international journal for theory and applications}, volume = {15}, journal = {Fractional calculus and applied analysis : an international journal for theory and applications}, number = {3}, publisher = {Versita}, address = {Warsaw}, issn = {1311-0454}, doi = {10.2478/s13540-012-0031-2}, pages = {426 -- 450}, year = {2012}, abstract = {We study analytically a generalized fractional Langevin equation. General formulas for calculation of variances and the mean square displacement are derived. Cases with a three parameter Mittag-Leffler frictional memory kernel are considered. Exact results in terms of the Mittag-Leffler type functions for the relaxation functions, average velocity and average particle displacement are obtained. The mean square displacement and variances are investigated analytically. Asymptotic behaviors of the particle in the short and long time limit are found. The model considered in this paper may be used for modeling anomalous diffusive processes in complex media including phenomena similar to single file diffusion or possible generalizations thereof. We show the importance of the initial conditions on the anomalous diffusive behavior of the particle.}, language = {en} } @article{SandevTomovskiDubbeldametal.2018, author = {Sandev, Trifce and Tomovski, Zivorad and Dubbeldam, Johan L. A. and Chechkin, Aleksei V.}, title = {Generalized diffusion-wave equation with memory kernel}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {52}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aaefa3}, pages = {22}, year = {2018}, abstract = {We study generalized diffusion-wave equation in which the second order time derivative is replaced by an integro-differential operator. It yields time fractional and distributed order time fractional diffusion-wave equations as particular cases. We consider different memory kernels of the integro-differential operator, derive corresponding fundamental solutions, specify the conditions of their non-negativity and calculate the mean squared displacement for all cases. In particular, we introduce and study generalized diffusion-wave equations with a regularized Prabhakar derivative of single and distributed orders. The equations considered can be used for modeling the broad spectrum of anomalous diffusion processes and various transitions between different diffusion regimes.}, language = {en} } @article{ShinCherstvyMetzler2014, author = {Shin, Jaeoh and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size}, series = {Soft Matter}, journal = {Soft Matter}, editor = {Metzler, Ralf}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, pages = {472 -- 488}, year = {2014}, abstract = {The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.}, language = {en} } @misc{ShinCherstvyMetzler2014, author = {Shin, Jaeoh and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-76961}, pages = {472 -- 488}, year = {2014}, abstract = {The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.}, language = {en} } @article{ShinCherstvyMetzler2015, author = {Shin, Jaeoh and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Self-subdiffusion in solutions of star-shaped crowders: non-monotonic effects of inter-particle interactions}, series = {New journal of physics : the open-access journal for physics}, volume = {17}, journal = {New journal of physics : the open-access journal for physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/17/11/113028}, pages = {12}, year = {2015}, abstract = {We examine by extensive computer simulations the self-diffusion of anisotropic star-like particles in crowded two-dimensional solutions. We investigate the implications of the area coverage fraction phi of the crowders and the crowder-crowder adhesion properties on the regime of transient anomalous diffusion. We systematically compute the mean squared displacement (MSD) of the particles, their time averaged MSD, and the effective diffusion coefficient. The diffusion is ergodic in the limit of long traces, such that the mean time averaged MSD converges towards the ensemble averaged MSD, and features a small residual amplitude spread of the time averaged MSD from individual trajectories. At intermediate time scales, we quantify the anomalous diffusion in the system. Also, we show that the translational-but not rotational-diffusivity of the particles Dis a nonmonotonic function of the attraction strength between them. Both diffusion coefficients decrease as the power law D(phi) similar to (1 - phi/phi*)(2 ... 2.4) with the area fraction phi occupied by the crowders and the critical value phi*. Our results might be applicable to rationalising the experimental observations of non-Brownian diffusion for a number of standard macromolecular crowders used in vitro to mimic the cytoplasmic conditions of living cells.}, language = {en} } @article{SliusarenkoVitaliSposinietal.2019, author = {Sliusarenko, Oleksii Yu and Vitali, Silvia and Sposini, Vittoria and Paradisi, Paolo and Chechkin, Aleksei V. and Castellani, Gastone and Pagnini, Gianni}, title = {Finite-energy Levy-type motion through heterogeneous ensemble of Brownian particles}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {52}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {9}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aafe90}, pages = {27}, year = {2019}, abstract = {Complex systems are known to display anomalous diffusion, whose signature is a space/time scaling x similar to t(delta) with delta not equal 1/2 in the probability density function (PDF). Anomalous diffusion can emerge jointly with both Gaussian, e.g. fractional Brownian motion, and power-law decaying distributions, e.g. Levy Flights or Levy Walks (LWs). Levy flights get anomalous scaling, but, being jumps of any size allowed even at short times, have infinite position variance, infinite energy and discontinuous paths. LWs, which are based on random trapping events, overcome these limitations: they resemble a Levy-type power-law distribution that is truncated in the large displacement range and have finite moments, finite energy and, even with discontinuous velocity, they are continuous. However, LWs do not take into account the role of strong heterogeneity in many complex systems, such as biological transport in the crowded cell environment. In this work we propose and discuss a model describing a heterogeneous ensemble of Brownian particles (HEBP). Velocity of each single particle obeys a standard underdamped Langevin equation for the velocity, with linear friction term and additive Gaussian noise. Each particle is characterized by its own relaxation time and velocity diffusivity. We show that, for proper distributions of relaxation time and velocity diffusivity, the HEBP resembles some LW statistical features, in particular power-law decaying PDF, long-range correlations and anomalous diffusion, at the same time keeping finite position moments and finite energy. The main differences between the HEBP model and two different LWs are investigated, finding that, even when both velocity and position PDFs are similar, they differ in four main aspects: (i) LWs are biscaling, while HEBP is monoscaling; (ii) a transition from anomalous (delta = 1/2) to normal (delta = 1/2) diffusion in the long-time regime is seen in the HEBP and not in LWs; (iii) the power-law index of the position PDF and the space/time diffusion scaling are independent in the HEBP, while they both depend on the scaling of the interevent time PDF in LWs; (iv) at variance with LWs, our HEBP model obeys a fluctuation-dissipation theorem.}, language = {en} } @misc{SposiniMetzlerOshanin2019, author = {Sposini, Vittoria and Metzler, Ralf and Oshanin, Gleb}, title = {Single-trajectory spectral analysis of scaled Brownian motion}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {753}, issn = {1866-8372}, doi = {10.25932/publishup-43652}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-436522}, pages = {16}, year = {2019}, abstract = {Astandard approach to study time-dependent stochastic processes is the power spectral density (PSD), an ensemble-averaged property defined as the Fourier transform of the autocorrelation function of the process in the asymptotic limit of long observation times, T → ∞. In many experimental situations one is able to garner only relatively few stochastic time series of finite T, such that practically neither an ensemble average nor the asymptotic limit T → ∞ can be achieved. To accommodate for a meaningful analysis of such finite-length data we here develop the framework of single-trajectory spectral analysis for one of the standard models of anomalous diffusion, scaled Brownian motion.Wedemonstrate that the frequency dependence of the single-trajectory PSD is exactly the same as for standard Brownian motion, which may lead one to the erroneous conclusion that the observed motion is normal-diffusive. However, a distinctive feature is shown to be provided by the explicit dependence on the measurement time T, and this ageing phenomenon can be used to deduce the anomalous diffusion exponent.Wealso compare our results to the single-trajectory PSD behaviour of another standard anomalous diffusion process, fractional Brownian motion, and work out the commonalities and differences. Our results represent an important step in establishing singletrajectory PSDs as an alternative (or complement) to analyses based on the time-averaged mean squared displacement.}, language = {en} } @article{SposiniMetzlerOshanin2019, author = {Sposini, Vittoria and Metzler, Ralf and Oshanin, Gleb}, title = {Single-trajectory spectral analysis of scaled Brownian motion}, series = {New Journal of Physics}, volume = {21}, journal = {New Journal of Physics}, publisher = {Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics}, address = {Bad Honnef und London}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab2f52}, pages = {16}, year = {2019}, abstract = {Astandard approach to study time-dependent stochastic processes is the power spectral density (PSD), an ensemble-averaged property defined as the Fourier transform of the autocorrelation function of the process in the asymptotic limit of long observation times, T → ∞. In many experimental situations one is able to garner only relatively few stochastic time series of finite T, such that practically neither an ensemble average nor the asymptotic limit T → ∞ can be achieved. To accommodate for a meaningful analysis of such finite-length data we here develop the framework of single-trajectory spectral analysis for one of the standard models of anomalous diffusion, scaled Brownian motion.Wedemonstrate that the frequency dependence of the single-trajectory PSD is exactly the same as for standard Brownian motion, which may lead one to the erroneous conclusion that the observed motion is normal-diffusive. However, a distinctive feature is shown to be provided by the explicit dependence on the measurement time T, and this ageing phenomenon can be used to deduce the anomalous diffusion exponent.Wealso compare our results to the single-trajectory PSD behaviour of another standard anomalous diffusion process, fractional Brownian motion, and work out the commonalities and differences. Our results represent an important step in establishing singletrajectory PSDs as an alternative (or complement) to analyses based on the time-averaged mean squared displacement.}, language = {en} } @phdthesis{Thapa2020, author = {Thapa, Samudrajit}, title = {Deciphering anomalous diffusion in complex systems using Bayesian inference and large deviation theory}, pages = {xx, 186}, year = {2020}, abstract = {The development of methods such as super-resolution microscopy (Nobel prize in Chemistry, 2014) and multi-scale computer modelling (Nobel prize in Chemistry, 2013) have provided scientists with powerful tools to study microscopic systems. Sub-micron particles or even fluorescently labelled single molecules can now be tracked for long times in a variety of systems such as living cells, biological membranes, colloidal solutions etc. at spatial and temporal resolutions previously inaccessible. Parallel to such single-particle tracking experiments, super-computing techniques enable simulations of large atomistic or coarse-grained systems such as biologically relevant membranes or proteins from picoseconds to seconds, generating large volume of data. These have led to an unprecedented rise in the number of reported cases of anomalous diffusion wherein the characteristic features of Brownian motion—namely linear growth of the mean squared displacement with time and the Gaussian form of the probability density function (PDF) to find a particle at a given position at some fixed time—are routinely violated. This presents a big challenge in identifying the underlying stochastic process and also estimating the corresponding parameters of the process to completely describe the observed behaviour. Finding the correct physical mechanism which leads to the observed dynamics is of paramount importance, for example, to understand the first-arrival time of transcription factors which govern gene regulation, or the survival probability of a pathogen in a biological cell post drug administration. Statistical Physics provides useful methods that can be applied to extract such vital information. This cumulative dissertation, based on five publications, focuses on the development, implementation and application of such tools with special emphasis on Bayesian inference and large deviation theory. Together with the implementation of Bayesian model comparison and parameter estimation methods for models of diffusion, complementary tools are developed based on different observables and large deviation theory to classify stochastic processes and gather pivotal information. Bayesian analysis of the data of micron-sized particles traced in mucin hydrogels at different pH conditions unveiled several interesting features and we gained insights into, for example, how in going from basic to acidic pH, the hydrogel becomes more heterogeneous and phase separation can set in, leading to observed non-ergodicity (non-equivalence of time and ensemble averages) and non-Gaussian PDF. With large deviation theory based analysis we could detect, for instance, non-Gaussianity in seeming Brownian diffusion of beads in aqueous solution, anisotropic motion of the beads in mucin at neutral pH conditions, and short-time correlations in climate data. Thus through the application of the developed methods to biological and meteorological datasets crucial information is garnered about the underlying stochastic processes and significant insights are obtained in understanding the physical nature of these systems.}, language = {en} } @article{ThapaWyłomańskaSikoraetal.2021, author = {Thapa, Samudrajit and Wyłomańska, Agnieszka and Sikora, Grzegorz and Wagner, Caroline E. and Krapf, Diego and Kantz, Holger and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories}, series = {New Journal of Physics}, volume = {23}, journal = {New Journal of Physics}, publisher = {Dt. Physikalische Ges. ; IOP}, address = {Bad Honnef ; London}, issn = {1367-2630}, doi = {10.1088/1367-2630/abd50e}, pages = {22}, year = {2021}, abstract = {Extensive time-series encoding the position of particles such as viruses, vesicles, or individualproteins are routinely garnered insingle-particle tracking experiments or supercomputing studies.They contain vital clues on how viruses spread or drugs may be delivered in biological cells.Similar time-series are being recorded of stock values in financial markets and of climate data.Such time-series are most typically evaluated in terms of time-averaged mean-squareddisplacements (TAMSDs), which remain random variables for finite measurement times. Theirstatistical properties are different for differentphysical stochastic processes, thus allowing us toextract valuable information on the stochastic process itself. To exploit the full potential of thestatistical information encoded in measured time-series we here propose an easy-to-implementand computationally inexpensive new methodology, based on deviations of the TAMSD from itsensemble average counterpart. Specifically, we use the upper bound of these deviations forBrownian motion (BM) to check the applicability of this approach to simulated and real data sets.By comparing the probability of deviations fordifferent data sets, we demonstrate how thetheoretical bound for BM reveals additional information about observed stochastic processes. Weapply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracerbeads measured in mucin hydrogels, and of geographic surface temperature anomalies. Ouranalysis shows how the large-deviation properties can be efficiently used as a simple yet effectiveroutine test to reject the BM hypothesis and unveil relevant information on statistical propertiessuch as ergodicity breaking and short-time correlations.}, language = {en} } @misc{ThapaWyłomańskaSikoraetal.2021, author = {Thapa, Samudrajit and Wyłomańska, Agnieszka and Sikora, Grzegorz and Wagner, Caroline E. and Krapf, Diego and Kantz, Holger and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1118}, issn = {1866-8372}, doi = {10.25932/publishup-49349}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-493494}, pages = {24}, year = {2021}, abstract = {Extensive time-series encoding the position of particles such as viruses, vesicles, or individualproteins are routinely garnered insingle-particle tracking experiments or supercomputing studies.They contain vital clues on how viruses spread or drugs may be delivered in biological cells.Similar time-series are being recorded of stock values in financial markets and of climate data.Such time-series are most typically evaluated in terms of time-averaged mean-squareddisplacements (TAMSDs), which remain random variables for finite measurement times. Theirstatistical properties are different for differentphysical stochastic processes, thus allowing us toextract valuable information on the stochastic process itself. To exploit the full potential of thestatistical information encoded in measured time-series we here propose an easy-to-implementand computationally inexpensive new methodology, based on deviations of the TAMSD from itsensemble average counterpart. Specifically, we use the upper bound of these deviations forBrownian motion (BM) to check the applicability of this approach to simulated and real data sets.By comparing the probability of deviations fordifferent data sets, we demonstrate how thetheoretical bound for BM reveals additional information about observed stochastic processes. Weapply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracerbeads measured in mucin hydrogels, and of geographic surface temperature anomalies. Ouranalysis shows how the large-deviation properties can be efficiently used as a simple yet effectiveroutine test to reject the BM hypothesis and unveil relevant information on statistical propertiessuch as ergodicity breaking and short-time correlations.}, language = {en} } @article{VilkAghionNathanetal.2022, author = {Vilk, Ohad and Aghion, Erez and Nathan, Ran and Toledo, Sivan and Metzler, Ralf and Assaf, Michael}, title = {Classification of anomalous diffusion in animal movement data using power spectral analysis}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {55}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {33}, publisher = {IOP Publishing}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/ac7e8f}, pages = {16}, year = {2022}, abstract = {The field of movement ecology has seen a rapid increase in high-resolution data in recent years, leading to the development of numerous statistical and numerical methods to analyse relocation trajectories. Data are often collected at the level of the individual and for long periods that may encompass a range of behaviours. Here, we use the power spectral density (PSD) to characterise the random movement patterns of a black-winged kite (Elanus caeruleus) and a white stork (Ciconia ciconia). The tracks are first segmented and clustered into different behaviours (movement modes), and for each mode we measure the PSD and the ageing properties of the process. For the foraging kite we find 1/f noise, previously reported in ecological systems mainly in the context of population dynamics, but not for movement data. We further suggest plausible models for each of the behavioural modes by comparing both the measured PSD exponents and the distribution of the single-trajectory PSD to known theoretical results and simulations.}, language = {en} } @article{VitaliSposiniSliusarenkoetal.2018, author = {Vitali, Silvia and Sposini, Vittoria and Sliusarenko, Oleksii and Paradisi, Paolo and Castellani, Gastone and Pagnini, Gianni}, title = {Langevin equation in complex media and anomalous diffusion}, series = {Interface : journal of the Royal Society}, volume = {15}, journal = {Interface : journal of the Royal Society}, number = {145}, publisher = {Royal Society}, address = {London}, issn = {1742-5689}, doi = {10.1098/rsif.2018.0282}, pages = {10}, year = {2018}, abstract = {The problem of biological motion is a very intriguing and topical issue. Many efforts are being focused on the development of novel modelling approaches for the description of anomalous diffusion in biological systems, such as the very complex and heterogeneous cell environment. Nevertheless, many questions are still open, such as the joint manifestation of statistical features in agreement with different models that can also be somewhat alternative to each other, e.g. continuous time random walk and fractional Brownian motion. To overcome these limitations, we propose a stochastic diffusion model with additive noise and linear friction force (linear Langevin equation), thus involving the explicit modelling of velocity dynamics. The complexity of the medium is parametrized via a population of intensity parameters (relaxation time and diffusivity of velocity), thus introducing an additional randomness, in addition to white noise, in the particle's dynamics. We prove that, for proper distributions of these parameters, we can get both Gaussian anomalous diffusion, fractional diffusion and its generalizations.}, language = {en} } @article{WangCherstvyChechkinetal.2020, author = {Wang, Wei and Cherstvy, Andrey G. and Chechkin, Aleksei V. and Thapa, Samudrajit and Seno, Flavio and Liu, Xianbin and Metzler, Ralf}, title = {Fractional Brownian motion with random diffusivity}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {53}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {47}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aba467}, pages = {34}, year = {2020}, abstract = {Numerous examples for a priori unexpected non-Gaussian behaviour for normal and anomalous diffusion have recently been reported in single-particle tracking experiments. Here, we address the case of non-Gaussian anomalous diffusion in terms of a random-diffusivity mechanism in the presence of power-law correlated fractional Gaussian noise. We study the ergodic properties of this model via examining the ensemble- and time-averaged mean-squared displacements as well as the ergodicity breaking parameter EB quantifying the trajectory-to-trajectory fluctuations of the latter. For long measurement times, interesting crossover behaviour is found as function of the correlation time tau characterising the diffusivity dynamics. We unveil that at short lag times the EB parameter reaches a universal plateau. The corresponding residual value of EB is shown to depend only on tau and the trajectory length. The EB parameter at long lag times, however, follows the same power-law scaling as for fractional Brownian motion. We also determine a corresponding plateau at short lag times for the discrete representation of fractional Brownian motion, absent in the continuous-time formulation. These analytical predictions are in excellent agreement with results of computer simulations of the underlying stochastic processes. Our findings can help distinguishing and categorising certain nonergodic and non-Gaussian features of particle displacements, as observed in recent single-particle tracking experiments.}, language = {en} } @article{WangSenoSokolovetal.2020, author = {Wang, Wei and Seno, Flavio and Sokolov, Igor M. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Unexpected crossovers in correlated random-diffusivity processes}, series = {New Journal of Physics}, volume = {22}, journal = {New Journal of Physics}, publisher = {Dt. Physikalische Ges.}, address = {Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/aba390}, pages = {17}, year = {2020}, abstract = {The passive and active motion of micron-sized tracer particles in crowded liquids and inside living biological cells is ubiquitously characterised by 'viscoelastic' anomalous diffusion, in which the increments of the motion feature long-ranged negative and positive correlations. While viscoelastic anomalous diffusion is typically modelled by a Gaussian process with correlated increments, so-called fractional Gaussian noise, an increasing number of systems are reported, in which viscoelastic anomalous diffusion is paired with non-Gaussian displacement distributions. Following recent advances in Brownian yet non-Gaussian diffusion we here introduce and discuss several possible versions of random-diffusivity models with long-ranged correlations. While all these models show a crossover from non-Gaussian to Gaussian distributions beyond some correlation time, their mean squared displacements exhibit strikingly different behaviours: depending on the model crossovers from anomalous to normal diffusion are observed, as well as a priori unexpected dependencies of the effective diffusion coefficient on the correlation exponent. Our observations of the non-universality of random-diffusivity viscoelastic anomalous diffusion are important for the analysis of experiments and a better understanding of the physical origins of 'viscoelastic yet non-Gaussian' diffusion.}, language = {en} }