Distinguishing between fractional Brownian motion with random and constant Hurst exponent using sample autocovariance-based statistics.

Fractional Brownian motion (FBM) is a canonical model for describing dynamics in various complex systems. It is characterized by the Hurst exponent, which is responsible for the correlation between FBM increments, its self-similarity property, and anomalous diffusion behavior. However, recent research indicates that the classical model may be insufficient in describing experimental observations when the anomalous diffusion exponent varies from trajectory to trajectory. As a result, modifications of the classical FBM have been considered in the literature, with a natural extension being the FBM with a random Hurst exponent. In this paper, we discuss the problem of distinguishing between two models: (i) FBM with the constant Hurst exponent and (ii) FBM with random Hurst exponent, by analyzing the probabilistic properties of statistics represented by the quadratic forms. These statistics have recently found application in Gaussian processes and have proven to serve as efficient tools for hypothesis testing. Here, we examine two statistics-the sample autocovariance function and the empirical anomaly measure-utilizing the correlation properties of the considered models. Based on these statistics, we introduce a testing procedure to differentiate between the two models. We present analytical and simulation results considering the two-point and beta distributions as exemplary distributions of the random Hurst exponent. Finally, to demonstrate the utility of the presented methodology, we analyze real-world datasets from the financial market and single particle tracking experiment in biological gels.

Erratum: "Fractional Brownian motion with random Hurst exponent: Accelerating diffusion and persistence transitions" [Chaos 32, 093114 (2022)].

We point out a minor mistake in Fig. 10 in the published version of our paper [M. Balcerek et al., Chaos 32, 093114 (2022)]. The conclusions drawn from the illustration remain the same.

Fractional Brownian motion with random Hurst exponent: Accelerating diffusion and persistence transitions.

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.

Objective comparison of methods to decode anomalous diffusion.

Deviations from Brownian motion leading to anomalous diffusion are found in transport dynamics from quantum physics to life sciences. The characterization of anomalous diffusion from the measurement of an individual trajectory is a challenging task, which traditionally relies on calculating the trajectory mean squared displacement. However, this approach breaks down for cases of practical interest, e.g., short or noisy trajectories, heterogeneous behaviour, or non-ergodic processes. Recently, several new approaches have been proposed, mostly building on the ongoing machine-learning revolution. To perform an objective comparison of methods, we gathered the community and organized an open competition, the Anomalous Diffusion challenge (AnDi). Participating teams applied their algorithms to a commonly-defined dataset including diverse conditions. Although no single method performed best across all scenarios, machine-learning-based approaches achieved superior performance for all tasks. The discussion of the challenge results provides practical advice for users and a benchmark for developers.

Transient superdiffusion of polydisperse vacuoles in highly motile amoeboid cells.

We perform a detailed statistical analysis of diffusive trajectories of membrane-enclosed vesicles (vacuoles) in the supercrowded cytoplasm of living Acanthamoeba castellanii cells. From the vacuole traces recorded in the center-of-area frame of moving amoebae, we examine the statistics of the time-averaged mean-squared displacements of vacuoles, their generalized diffusion coefficients and anomalous scaling exponents, the ergodicity breaking parameter, the non-Gaussian features of displacement distributions of vacuoles, the displacement autocorrelation function, as well as the distributions of speeds and positions of vacuoles inside the amoeba cells. Our findings deliver novel insights into the internal dynamics of cellular structures in these infectious pathogens.

Non-Gaussian, non-ergodic, and non-Fickian diffusion of tracers in mucin hydrogels.

Native mucus is polymer-based soft-matter material of paramount biological importance. How non-Gaussian and non-ergodic is the diffusive spreading of pathogens in mucus? We study the passive, thermally driven motion of micron-sized tracers in hydrogels of mucins, the main polymeric component of mucus. We report the results of the Bayesian analysis for ranking several diffusion models for a set of tracer trajectories [C. E. Wagner et al., Biomacromolecules, 2017, 18, 3654]. The models with "diffusing diffusivity", fractional and standard Brownian motion are used. The likelihood functions and evidences of each model are computed, ranking the significance of each model for individual traces. We find that viscoelastic anomalous diffusion is often most probable, followed by Brownian motion, while the model with a diffusing diffusion coefficient is only realised rarely. Our analysis also clarifies the distribution of time-averaged displacements, correlations of scaling exponents and diffusion coefficients, and the degree of non-Gaussianity of displacements at varying pH levels. Weak ergodicity breaking is also quantified. We conclude that-consistent with the original study-diffusion of tracers in the mucin gels is most non-Gaussian and non-ergodic at low pH that corresponds to the most heterogeneous networks. Using the Bayesian approach with the nested-sampling algorithm, together with the quantitative analysis of multiple statistical measures, we report new insights into possible physical mechanisms of diffusion in mucin gels.

Bayesian analysis of single-particle tracking data using the nested-sampling algorithm: maximum-likelihood model selection applied to stochastic-diffusivity data.

We employ Bayesian statistics using the nested-sampling algorithm to compare and rank multiple models of ergodic diffusion (including anomalous diffusion) as well as to assess their optimal parameters for in silico-generated and real time-series. We focus on the recently-introduced model of Brownian motion with "diffusing diffusivity"-giving rise to widely-observed non-Gaussian displacement statistics-and its comparison to Brownian and fractional Brownian motion, also for the time-series with some measurement noise. We conduct this model-assessment analysis using Bayesian statistics and the nested-sampling algorithm on the level of individual particle trajectories. We evaluate relative model probabilities and compute best-parameter sets for each diffusion model, comparing the estimated parameters to the true ones. We test the performance of the nested-sampling algorithm and its predictive power both for computer-generated (idealised) trajectories as well as for real single-particle-tracking trajectories. Our approach delivers new important insight into the objective selection of the most suitable stochastic model for a given time-series. We also present first model-ranking results in application to experimental data of tracer diffusion in polymer-based hydrogels.

Time averages and their statistical variation for the Ornstein-Uhlenbeck process: Role of initial particle distributions and relaxation to stationarity.

How ergodic is diffusion under harmonic confinements? How strongly do ensemble- and time-averaged displacements differ for a thermally-agitated particle performing confined motion for different initial conditions? We here study these questions for the generic Ornstein-Uhlenbeck (OU) process and derive the analytical expressions for the second and fourth moment. These quantifiers are particularly relevant for the increasing number of single-particle tracking experiments using optical traps. For a fixed starting position, we discuss the definitions underlying the ensemble averages. We also quantify effects of equilibrium and nonequilibrium initial particle distributions onto the relaxation properties and emerging nonequivalence of the ensemble- and time-averaged displacements (even in the limit of long trajectories). We derive analytical expressions for the ergodicity breaking parameter quantifying the amplitude scatter of individual time-averaged trajectories, both for equilibrium and out-of-equilibrium initial particle positions, in the entire range of lag times. Our analytical predictions are in excellent agreement with results of computer simulations of the Langevin equation in a parabolic potential. We also examine the validity of the Einstein relation for the ensemble- and time-averaged moments of the OU-particle. Some physical systems, in which the relaxation and nonergodic features we unveiled may be observable, are discussed.