We propose a generalization of the widely used fractional Brownian motion (FBM), memory-multi-FBM (MMFBM), to describe viscoelastic or persistent anomalous diffusion with time-dependent memory exponent α(t ) in a changing environment. In MMFBM the built-in, long-range memory is continuously modulated by α(t ). We derive the essential statistical properties of MMFBM such as its response function, mean-squared displacement (MSD), autocovariance function, and Gaussian distribution. In contrast to existing forms of FBM with time-varying memory exponents but a reset memory structure, the instantaneous dynamic of MMFBM is influenced by the process history, e.g., we show that after a steplike change of α(t ) the scaling exponent of the MSD after the α step may be determined by the value of α(t ) before the change. MMFBM is a versatile and useful process for correlated physical systems with nonequilibrium initial conditions in a changing environment.
In this paper we propose an algorithm to distinguish between light- and heavy-tailed probability laws underlying random datasets. The idea of the algorithm, which is visual and easy to implement, is to check whether the underlying law belongs to the domain of attraction of the Gaussian or non-Gaussian stable distribution by examining its rate of convergence. The method allows to discriminate between stable and various non-stable distributions. The test allows to differentiate between distributions, which appear the same according to standard Kolmogorov-Smirnov test. In particular, it helps to distinguish between stable and Student's t probability laws as well as between the stable and tempered stable, the cases which are considered in the literature as very cumbersome. Finally, we illustrate the procedure on plasma data to identify cases with so-called L-H transition.
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.
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.
We address the problem of recognizing alpha-stable Levy distribution with Levy index close to 2 from experimental data. We are interested in the case when the sample size of available data is not large, thus the power law asymptotics of the distribution is not clearly detectable, and the shape of the empirical probability density function is close to a Gaussian. We propose a testing procedure combining a simple visual test based on empirical fourth moment with the Anderson-Darling and Jarque-Bera statistical tests and we check the efficiency of the method on simulated data. Furthermore, we apply our method to the analysis of turbulent plasma density and potential fluctuations measured in the stellarator-type fusion device and demonstrate that the phenomenon of the L-H transition from low confinement, L mode, to a high confinement, H mode, which occurs in this device is accompanied by the transition from Levy to Gaussian fluctuation statistics.
In this paper we study the large deviations of time averaged mean square displacement (TAMSD) for Gaussian processes. The theory of large deviations is related to the exponential decay of probabilities of large fluctuations in random systems. From the mathematical point of view a given statistics satisfies the large deviation principle, if the probability that it belongs to a certain range decreases exponentially. The TAMSD is one of the main statistics used in the problem of anomalous diffusion detection. Applying the theory of generalized chi-squared distribution and sub-gamma random variables we prove the upper bound for large deviations of TAMSD for Gaussian processes. As a special case we consider fractional Brownian motion, one of the most popular models of anomalous diffusion. Moreover, we derive the upper bound for large deviations of the estimator for the anomalous diffusion exponent. (C) 2018 Elsevier B.V. All rights reserved.
Fractional Brownian motion, a Gaussian non-Markovian self-similar process with stationary long-correlated increments, has been identified to give rise to the anomalous diffusion behavior in a great variety of physical systems. The correlation and diffusion properties of this random motion are fully characterized by its index of self-similarity or the Hurst exponent.
However, recent single-particle tracking experiments in biological cells revealed highly complicated anomalous diffusion phenomena that cannot be attributed to a class of self-similar random processes.
Inspired by these observations, we here study the process that preserves the properties of the fractional Brownian motion at a single trajectory level; however, the Hurst index randomly changes from trajectory to trajectory.
We provide a general mathematical framework for analytical, numerical, and statistical analysis of the fractional Brownian motion with the random Hurst exponent.
The explicit formulas for probability density function, mean-squared displacement, and autocovariance function of the increments are presented for three generic distributions of the Hurst exponent, namely, two-point, uniform, and beta distributions.
The important features of the process studied here are accelerating diffusion and persistence transition, which we demonstrate analytically and numerically.