@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} } @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 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} } @article{GodecChechkinBarkaietal.2014, author = {Godec, Aljaz and Chechkin, Aleksei V. and Barkai, Eli and Kantz, Holger and Metzler, Ralf}, title = {Localisation and universal fluctuations in ultraslow diffusion processes}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {47}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {49}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/47/49/492002}, pages = {10}, year = {2014}, abstract = {We study ultraslow diffusion processes with logarithmic mean squared displacement (MSD) < x(2)(t)> 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{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{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{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} } @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{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{GhoshCherstvyMetzler2014, author = {Ghosh, Surya K. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Non-universal tracer diffusion in crowded media of non-inert obstacles}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-77128}, pages = {1847 -- 1858}, year = {2014}, abstract = {We study the diffusion of a tracer particle, which moves in continuum space between a lattice of excluded volume, immobile non-inert obstacles. In particular, we analyse how the strength of the tracer-obstacle interactions and the volume occupancy of the crowders alter the diffusive motion of the tracer. From the details of partitioning of the tracer diffusion modes between trapping states when bound to obstacles and bulk diffusion, we examine the degree of localisation of the tracer in the lattice of crowders. We study the properties of the tracer diffusion in terms of the ensemble and time averaged mean squared displacements, the trapping time distributions, the amplitude variation of the time averaged mean squared displacements, and the non-Gaussianity parameter of the diffusing tracer. We conclude that tracer-obstacle adsorption and binding triggers a transient anomalous diffusion. From a very narrow spread of recorded individual time averaged trajectories we exclude continuous type random walk processes as the underlying physical model of the tracer diffusion in our system. For moderate tracer-crowder attraction the motion is found to be fully ergodic, while at stronger attraction strength a transient disparity between ensemble and time averaged mean squared displacements occurs. We also put our results into perspective with findings from experimental single-particle tracking and simulations of the diffusion of tagged tracers in dense crowded suspensions. Our results have implications for the diffusion, transport, and spreading of chemical components in highly crowded environments inside living cells and other structured liquids.}, language = {en} } @article{GhoshCherstvyMetzler2014, author = {Ghosh, Surya K. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Non-universal tracer diffusion in crowded media of non-inert obstacles}, series = {Physical Chemistry Chemical Physics}, volume = {3}, journal = {Physical Chemistry Chemical Physics}, number = {17}, editor = {Metzler, Ralf}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, pages = {1847 -- 1858}, year = {2014}, abstract = {We study the diffusion of a tracer particle, which moves in continuum space between a lattice of excluded volume, immobile non-inert obstacles. In particular, we analyse how the strength of the tracer-obstacle interactions and the volume occupancy of the crowders alter the diffusive motion of the tracer. From the details of partitioning of the tracer diffusion modes between trapping states when bound to obstacles and bulk diffusion, we examine the degree of localisation of the tracer in the lattice of crowders. We study the properties of the tracer diffusion in terms of the ensemble and time averaged mean squared displacements, the trapping time distributions, the amplitude variation of the time averaged mean squared displacements, and the non-Gaussianity parameter of the diffusing tracer. We conclude that tracer-obstacle adsorption and binding triggers a transient anomalous diffusion. From a very narrow spread of recorded individual time averaged trajectories we exclude continuous type random walk processes as the underlying physical model of the tracer diffusion in our system. For moderate tracer-crowder attraction the motion is found to be fully ergodic, while at stronger attraction strength a transient disparity between ensemble and time averaged mean squared displacements occurs. We also put our results into perspective with findings from experimental single-particle tracking and simulations of the diffusion of tagged tracers in dense crowded suspensions. Our results have implications for the diffusion, transport, and spreading of chemical components in highly crowded environments inside living cells and other structured liquids.}, language = {en} }