@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{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{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}
}
@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}
}
@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}
}
@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}
}
@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}
}