Refine
Has Fulltext
- no (38)
Year of publication
Document Type
- Article (38) (remove)
Language
- English (38)
Is part of the Bibliography
- yes (38)
Keywords
- anomalous diffusion (38) (remove)
Institute
Recent advances in single particle tracking and supercomputing techniques demonstrate the emergence of normal or anomalous, viscoelastic diffusion in conjunction with non-Gaussian distributions in soft, biological, and active matter systems. We here formulate a stochastic model based on a generalised Langevin equation in which non-Gaussian shapes of the probability density function and normal or anomalous diffusion have a common origin, namely a random parametrisation of the stochastic force. We perform a detailed analysis demonstrating how various types of parameter distributions for the memory kernel result in exponential, power law, or power-log law tails of the memory functions. The studied system is also shown to exhibit a further unusual property: the velocity has a Gaussian one point probability density but non-Gaussian joint distributions. This behaviour is reflected in the relaxation from a Gaussian to a non-Gaussian distribution observed for the position variable. We show that our theoretical results are in excellent agreement with stochastic simulations.
We show that the codifference is a useful tool in studying the ergodicity breaking and non-Gaussianity properties of stochastic time series. While the codifference is a measure of dependence that was previously studied mainly in the context of stable processes, we here extend its range of applicability to random-parameter and diffusing-diffusivity models which are important in contemporary physics, biology and financial engineering. We prove that the codifference detects forms of dependence and ergodicity breaking which are not visible from analysing the covariance and correlation functions. We also discuss a related measure of dispersion, which is a nonlinear analogue of the mean squared displacement.
We introduce and study a Lévy walk (LW) model of particle spreading with a finite propagation speed combined with soft resets, stochastically occurring periods in which an harmonic external potential is switched on and forces the particle towards a specific position. Soft resets avoid instantaneous relocation of particles that in certain physical settings may be considered unphysical. Moreover, soft resets do not have a specific resetting point but lead the particle towards a resetting point by a restoring Hookean force. Depending on the exact choice for the LW waiting time density and the probability density of the periods when the harmonic potential is switched on, we demonstrate a rich emerging response behaviour including ballistic motion and superdiffusion. When the confinement periods of the soft-reset events are dominant, we observe a particle localisation with an associated non-equilibrium steady state. In this case the stationary particle probability density function turns out to acquire multimodal states. Our derivations are based on Markov chain ideas and LWs with multiple internal states, an approach that may be useful and flexible for the investigation of other generalised random walks with soft and hard resets. The spreading efficiency of soft-rest LWs is characterised by the first-passage time statistic.
To analyze stochastic processes, one often uses integral transform (Fourier and Laplace) methods. However, for the time-space coupled cases, e.g. the Levy walk, sometimes the integral transform method may fail. Here we provide a Hermite polynomial expansion approach, being complementary to the integral transform method, to the Levy walk. Two approaches are compared for some already known results. We also consider the generalized Levy walk with parameter dependent velocity. Namely, we consider the Levy walk with velocity which depends on the walking length or on the duration of each step. Some interesting features of the generalized Levy walk are observed, including the special shapes of the probability density function, the first passage time distributions, and various diffusive behaviors of the mean squared displacement.
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.
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.
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.
Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories
(2021)
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.
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.
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.
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.
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.
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.
Velocity and displacement correlation functions for fractional generalized Langevin equations
(2012)
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.
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.
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.
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.
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.
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.
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.