TY - JOUR A1 - Vitali, Silvia A1 - Sposini, Vittoria A1 - Sliusarenko, Oleksii A1 - Paradisi, Paolo A1 - Castellani, Gastone A1 - Pagnini, Gianni T1 - Langevin equation in complex media and anomalous diffusion JF - Interface : journal of the Royal Society N2 - 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. KW - anomalous diffusion KW - heterogeneous media KW - biological transport KW - Gaussian processes KW - space-time fractional diffusion equation KW - fractional Brownian motion Y1 - 2018 U6 - https://doi.org/10.1098/rsif.2018.0282 SN - 1742-5689 SN - 1742-5662 VL - 15 IS - 145 PB - Royal Society CY - London ER - TY - JOUR A1 - Vilk, Ohad A1 - Aghion, Erez A1 - Nathan, Ran A1 - Toledo, Sivan A1 - Metzler, Ralf A1 - Assaf, Michael T1 - Classification of anomalous diffusion in animal movement data using power spectral analysis JF - Journal of physics : A, Mathematical and theoretical N2 - The field of movement ecology has seen a rapid increase in high-resolution data in recent years, leading to the development of numerous statistical and numerical methods to analyse relocation trajectories. Data are often collected at the level of the individual and for long periods that may encompass a range of behaviours. Here, we use the power spectral density (PSD) to characterise the random movement patterns of a black-winged kite (Elanus caeruleus) and a white stork (Ciconia ciconia). The tracks are first segmented and clustered into different behaviours (movement modes), and for each mode we measure the PSD and the ageing properties of the process. For the foraging kite we find 1/f noise, previously reported in ecological systems mainly in the context of population dynamics, but not for movement data. We further suggest plausible models for each of the behavioural modes by comparing both the measured PSD exponents and the distribution of the single-trajectory PSD to known theoretical results and simulations. KW - diffusion KW - anomalous diffusion KW - power spectral analysis KW - ecological KW - movement data Y1 - 2022 U6 - https://doi.org/10.1088/1751-8121/ac7e8f SN - 1751-8113 SN - 1751-8121 VL - 55 IS - 33 PB - IOP Publishing CY - Bristol ER - TY - JOUR A1 - Thapa, Samudrajit A1 - Wyłomańska, Agnieszka A1 - Sikora, Grzegorz A1 - Wagner, Caroline E. A1 - Krapf, Diego A1 - Kantz, Holger A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories JF - New Journal of Physics N2 - Extensive time-series encoding the position of particles such as viruses, vesicles, or individualproteins are routinely garnered insingle-particle tracking experiments or supercomputing studies.They contain vital clues on how viruses spread or drugs may be delivered in biological cells.Similar time-series are being recorded of stock values in financial markets and of climate data.Such time-series are most typically evaluated in terms of time-averaged mean-squareddisplacements (TAMSDs), which remain random variables for finite measurement times. Theirstatistical properties are different for differentphysical stochastic processes, thus allowing us toextract valuable information on the stochastic process itself. To exploit the full potential of thestatistical information encoded in measured time-series we here propose an easy-to-implementand computationally inexpensive new methodology, based on deviations of the TAMSD from itsensemble average counterpart. Specifically, we use the upper bound of these deviations forBrownian motion (BM) to check the applicability of this approach to simulated and real data sets.By comparing the probability of deviations fordifferent data sets, we demonstrate how thetheoretical bound for BM reveals additional information about observed stochastic processes. Weapply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracerbeads measured in mucin hydrogels, and of geographic surface temperature anomalies. Ouranalysis shows how the large-deviation properties can be efficiently used as a simple yet effectiveroutine test to reject the BM hypothesis and unveil relevant information on statistical propertiessuch as ergodicity breaking and short-time correlations. KW - diffusion KW - anomalous diffusion KW - large-deviation statistic KW - time-averaged mean squared displacement KW - Chebyshev inequality Y1 - 2020 U6 - https://doi.org/10.1088/1367-2630/abd50e SN - 1367-2630 VL - 23 PB - Dt. Physikalische Ges. ; IOP CY - Bad Honnef ; London ER - TY - GEN A1 - Thapa, Samudrajit A1 - Wyłomańska, Agnieszka A1 - Sikora, Grzegorz A1 - Wagner, Caroline E. A1 - Krapf, Diego A1 - Kantz, Holger A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - Extensive time-series encoding the position of particles such as viruses, vesicles, or individualproteins are routinely garnered insingle-particle tracking experiments or supercomputing studies.They contain vital clues on how viruses spread or drugs may be delivered in biological cells.Similar time-series are being recorded of stock values in financial markets and of climate data.Such time-series are most typically evaluated in terms of time-averaged mean-squareddisplacements (TAMSDs), which remain random variables for finite measurement times. Theirstatistical properties are different for differentphysical stochastic processes, thus allowing us toextract valuable information on the stochastic process itself. To exploit the full potential of thestatistical information encoded in measured time-series we here propose an easy-to-implementand computationally inexpensive new methodology, based on deviations of the TAMSD from itsensemble average counterpart. Specifically, we use the upper bound of these deviations forBrownian motion (BM) to check the applicability of this approach to simulated and real data sets.By comparing the probability of deviations fordifferent data sets, we demonstrate how thetheoretical bound for BM reveals additional information about observed stochastic processes. Weapply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracerbeads measured in mucin hydrogels, and of geographic surface temperature anomalies. Ouranalysis shows how the large-deviation properties can be efficiently used as a simple yet effectiveroutine test to reject the BM hypothesis and unveil relevant information on statistical propertiessuch as ergodicity breaking and short-time correlations. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1118 KW - diffusion KW - anomalous diffusion KW - large-deviation statistic KW - time-averaged mean squared displacement KW - Chebyshev inequality Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-493494 SN - 1866-8372 IS - 1118 ER - TY - THES A1 - Thapa, Samudrajit T1 - Deciphering anomalous diffusion in complex systems using Bayesian inference and large deviation theory N2 - 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. N2 - Die Entwicklung von Methoden wie der superauflösenden Mikroskopie (Nobelpreis für Chemie, 2014) und der Multiskalen-Computermodellierung (Nobelpreis für Chemie, 2013) hat Wis- senschaftlern mächtige Werkzeuge zur Untersuchung mikroskopischer Systeme zur Verfügung gestellt. Submikrometer Partikel und sogar einzelne fluoreszent markierte Moleküle können heute über lange Beobachtungszeiten in einer Vielzahl von Systemen, wie z.B. lebenden Zellen, biologischen Membranen und kolloidalen Suspensionen, mit bisher unerreichter räumlicher und zeitlicher Auflösung verfolgt werden. Neben solchen Einzelpartikelverfolgungsexperi- menten ermöglichen Supercomputer die Simulation großer atomistischer oder coarse-grained Systeme, wie z..B. biologisch relevante Membranen oder Proteine, über wenige Picosekunden bis hin zu einigen Sekunden, wobei große Datenmengen produziert werden. Diese haben zu einem beispiellosen Anstieg in der Zahl berichteter Fälle von anomaler Diffusion geführt, bei welcher die charakteristischen Eigenschaften der Brownschen Diffusion—nämlich das lineare Wachstum der mittleren quadratischen Verschiebung mit der Zeit und die Gaußsche Form der Wahrscheinlichkeitsdichtefunktion ein Partikel an einem gegebenen Ort und zu gegebener Zeit zu finden—verletzt sind. Dies stellt eine große Herausforderung bei der Identifizierung des zugrundeliegenden stochastischen Prozesses und der Schätzung der zugehörigen Prozess- parameter dar, was zur vollständigen Beschreibung des beobachteten Verhaltens nötig ist. Das Auffinden des korrekten physikalischen Mechanismus, welcher zum beobachteten Verhal- ten führt, ist von überragender Bedeutung, z.B. beim Verständnis der, die Genregulation steuernden, first-arrival time von Transkriptionsfaktoren oder der Überlebensfunktion eines Pathogens in einer biologischen Zelle nach Medikamentengabe. Die statistische Physik stellt nützliche Methoden bereit, die angewendet werden können, um solch wichtige Informationen zu erhalten. Der Schwerpunkt der vorliegenden, auf fünf Publikationen basierenden, kumulativen Dissertation liegt auf der Entwicklung, Implementierung und Anwendung solcher Methoden, mit einem besonderen Schwerpunkt auf der Bayesschen Inferenz und der Theorie der großen Abweichungen. Zusammen mit der Implementierung eines Bayesschen Modellvergleichs und Methoden zur Parameterschätzung für Diffusionsmodelle werden ergänzende Methoden en- twickelt, welche auf unterschiedliche Observablen und der Theorie der großen Abweichungen basieren, um stochastische Prozesse zu klassifizieren und wichtige Informationen zu erhalten. Die Bayessche Analyse der Bewegungsdaten von Mikrometerpartikeln, welche in Mucin- hydrogelen mit verschiedenen pH-Werten verfolgt wurden, enthüllte mehrere interessante Eigenschaften und wir haben z.B. Einsichten darüber gewonnen, wie der Übergang von basischen zu sauren pH-Werten die Heterogenität des Hydrogels erhöht und Phasentrennung einsetzen kann, was zur beobachteten Nicht-Ergodizität (Inäquivalenz von Zeit- und En- semblemittelwert) und nicht-Gaußscher Wahrscheinlichkeitsdichte führt. Mit einer auf der Theorie der großen Abweichungen basierenden Analyse konnten wir, z.B., nicht-Gaußsches Verhalten bei der scheinbaren Brownschen Diffusion von Partikeln in wässriger Lösung, anisotrope Bewegung von Partikeln in Mucin bei neutralem pH-Wert und Kurzzeitkorrela- tionen in Klimadaten detektieren. Folglich werden durch die Anwendung der entwickelten Methoden auf biologische und meteorologische Daten entscheidende Informationen über die zugrundeliegenden stochastischen Prozesse gesammelt und bedeutende Erkenntnisse für das Verständnis der Eigenschaften dieser Systeme erhalten. KW - anomalous diffusion KW - Bayesian inference KW - large deviation theory KW - statistical physics Y1 - 2020 ER - TY - GEN A1 - Sposini, Vittoria A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Single-trajectory spectral analysis of scaled Brownian motion T2 - Postprints der Universität Potsdam Mathematisch-Naturwissenschaftliche Reihe N2 - Astandard approach to study time-dependent stochastic processes is the power spectral density (PSD), an ensemble-averaged property defined as the Fourier transform of the autocorrelation function of the process in the asymptotic limit of long observation times, T → ∞. In many experimental situations one is able to garner only relatively few stochastic time series of finite T, such that practically neither an ensemble average nor the asymptotic limit T → ∞ can be achieved. To accommodate for a meaningful analysis of such finite-length data we here develop the framework of single-trajectory spectral analysis for one of the standard models of anomalous diffusion, scaled Brownian motion.Wedemonstrate that the frequency dependence of the single-trajectory PSD is exactly the same as for standard Brownian motion, which may lead one to the erroneous conclusion that the observed motion is normal-diffusive. However, a distinctive feature is shown to be provided by the explicit dependence on the measurement time T, and this ageing phenomenon can be used to deduce the anomalous diffusion exponent.Wealso compare our results to the single-trajectory PSD behaviour of another standard anomalous diffusion process, fractional Brownian motion, and work out the commonalities and differences. Our results represent an important step in establishing singletrajectory PSDs as an alternative (or complement) to analyses based on the time-averaged mean squared displacement. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 753 KW - diffusion KW - anomalous diffusion KW - power spectral analysis KW - single trajectory analysis Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-436522 SN - 1866-8372 IS - 753 ER - TY - JOUR A1 - Sposini, Vittoria A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Single-trajectory spectral analysis of scaled Brownian motion JF - New Journal of Physics N2 - Astandard approach to study time-dependent stochastic processes is the power spectral density (PSD), an ensemble-averaged property defined as the Fourier transform of the autocorrelation function of the process in the asymptotic limit of long observation times, T → ∞. In many experimental situations one is able to garner only relatively few stochastic time series of finite T, such that practically neither an ensemble average nor the asymptotic limit T → ∞ can be achieved. To accommodate for a meaningful analysis of such finite-length data we here develop the framework of single-trajectory spectral analysis for one of the standard models of anomalous diffusion, scaled Brownian motion.Wedemonstrate that the frequency dependence of the single-trajectory PSD is exactly the same as for standard Brownian motion, which may lead one to the erroneous conclusion that the observed motion is normal-diffusive. However, a distinctive feature is shown to be provided by the explicit dependence on the measurement time T, and this ageing phenomenon can be used to deduce the anomalous diffusion exponent.Wealso compare our results to the single-trajectory PSD behaviour of another standard anomalous diffusion process, fractional Brownian motion, and work out the commonalities and differences. Our results represent an important step in establishing singletrajectory PSDs as an alternative (or complement) to analyses based on the time-averaged mean squared displacement. KW - diffusion KW - anomalous diffusion KW - power spectral analysis KW - single trajectory analysis Y1 - 2019 U6 - https://doi.org/10.1088/1367-2630/ab2f52 SN - 1367-2630 VL - 21 PB - Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics CY - Bad Honnef und London ER - TY - JOUR A1 - Sliusarenko, Oleksii Yu A1 - Vitali, Silvia A1 - Sposini, Vittoria A1 - Paradisi, Paolo A1 - Chechkin, Aleksei V. A1 - Castellani, Gastone A1 - Pagnini, Gianni T1 - Finite-energy Levy-type motion through heterogeneous ensemble of Brownian particles JF - Journal of physics : A, Mathematical and theoretical N2 - 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. KW - anomalous diffusion KW - heterogeneous ensemble of Brownian particles KW - biological transport KW - Langevin equation KW - Gaussian processes KW - fractional diffusion KW - Levy walk Y1 - 2019 U6 - https://doi.org/10.1088/1751-8121/aafe90 SN - 1751-8113 SN - 1751-8121 VL - 52 IS - 9 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Shin, Jaeoh A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf ED - Metzler, Ralf T1 - Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size JF - Soft Matter N2 - The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping–unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions. KW - gene-regulation kinetics KW - physiological consequences KW - spatial-organization KW - anomalous diffusion KW - folding kinetics KW - living cells KW - dna coiling KW - in-vitro KW - dynamics KW - mixtures Y1 - 2014 SN - 1744-683X SP - 472 EP - 488 PB - The Royal Society of Chemistry CY - Cambridge ER - TY - GEN A1 - Shin, Jaeoh A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size N2 - The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping–unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 185 KW - gene-regulation kinetics KW - physiological consequences KW - spatial-organization KW - anomalous diffusion KW - folding kinetics KW - living cells KW - dna coiling KW - in-vitro KW - dynamics KW - mixtures Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-76961 SP - 472 EP - 488 PB - The Royal Society of Chemistry CY - Cambridge ER -