An official website of the United States government
Stochastic processes exhibiting power-law slopes in the frequency domain are frequently well modeled by fractional Brownian motion (fBm), with the spectral slope at high frequencies being associated with the degree of small-scale roughness or fractal dimension. However, a broad class of real-world signals have a high-frequency slope, like fBm, but a plateau in the vicinity of zero frequency. This low-frequency plateau, it is shown, implies that the temporal integral of the process exhibits diffusive behavior, dispersing from its initial location at a constant rate. Such processes are not well modeled by fBm, which has a singularity at zero frequency corresponding to an unbounded rate of dispersion. A more appropriate stochastic model is a much lesser-known random process called the Matérn process, which is shown herein to be a damped version of fractional Brownian motion. This article first provides a thorough introduction to fractional Brownian motion, then examines the details of the Matérn process and its relationship to fBm. An algorithm for the simulation of the Matérn process in O(NlogN) operations is given. Unlike fBm, the Matérn process is found to provide an excellent match to modeling velocities from particle trajectories in an application to two-dimensional fluid turbulence.
more »
« less
Evaluating Gaussianity of heterogeneous fractional Brownian motion
Abstract Heterogeneous diffusion processes are prevalent in various fields, including the motion of proteins in living cells, the migratory movement of birds and mammals, and finance. These processes are often characterized by time-varying dynamics, where interactions with the environment evolve, and the system undergoes fluctuations in diffusivity. Moreover, in many complex systems anomalous diffusion is observed, where the mean square displacement (MSD) exhibits non-linear scaling with time. Among the models used to describe this phenomenon, fractional Brownian motion (FBM) is a widely applied stochastic process, particularly for systems exhibiting long-range temporal correlations. Although FBM is characterized by Gaussian increments, heterogeneous processes with FBM-like characteristics may deviate from Gaussianity. In this article, we study the non-Gaussian behavior of switching fractional Brownian motion (SFBM), a model in which the diffusivity of the FBM process varies while temporal correlations are maintained. To characterize non-Gaussianity, we evaluate the kurtosis, a common tool used to quantify deviations from the normal distribution. We derive exact expressions for the kurtosis of the considered heterogeneous anomalous diffusion process and investigate how it can identify non-Gaussian behavior. We also compare the kurtosis results with those obtained using the Hellinger distance, a classical measure of divergence between probability density functions. Through both analytical and numerical methods, we demonstrate the potential of kurtosis as a metric for detecting non-Gaussianity in heterogeneous anomalous diffusion processes.
more »
« less
- Award ID(s):
- 2102832
- PAR ID:
- 10581654
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Journal of Physics A: Mathematical and Theoretical
- ISSN:
- 1751-8113
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We introduce the stochastic process of incremental multifractional Brownian motion (IMFBM), which locally behaves like fractional Brownian motion with a given local Hurst exponent and diffusivity. When these parameters change as function of time the process responds to the evolution gradually: only new increments are governed by the new parameters, while still retaining a power-law dependence on the past of the process. We obtain the mean squared displacement and correlations of IMFBM which are given by elementary formulas. We also provide a comparison with simulations and introduce estimation methods for IMFBM. This mathematically simple process is useful in the description of anomalous diffusion dynamics in changing environments, e.g. in viscoelastic systems, or when an actively moving particle changes its degree of persistence or its mobility.more » « less
-
Fractional Brownian Motion (FBM) is a stochastic process with long-time correlations which has been used to model anomalous diffusion in numerous biological systems. Recently, it has been used to study the distribution of serotonergic fibers in the brain [1,2]. To better represent the biological process we are trying to simulate, we introduce the concept of branching FBM (bFBM). In this stochastic process, individual particles perform FBM but may randomly split into two. Here, we study bFBM in one space dimension in the subdiffusive and superdiffusive regimes, both in free space and on finite intervals with reflecting boundaries. We examine three possible types of behavior of the correlations (memory) at a branching event: both particles keep the memory of the previous steps, only one particle keeps the memory, and no particles keep the memory. We calculate the mean-square particle displacement, the corresponding probability distribution, and displacement correlation function. We find that the qualitative features of the bFBM process strongly depend on the type of branching event. We also discuss implications of our results for the distribution of serotonergic fibers, and we discuss possible future refinements of the model, including interactions between different fibers. [1] T. Vojta, S. Halladay, S. Skinner, S. Janusonis, T. Guggenberger, R. Metzler, Phys. Rev. E 102, 032108 (2020). [2] S. Janusonis, N. Detering, R. Metzler, T. Vojta, Front. Comput. Neurosci. 14, 56 (2020). This work was supported in part by the NSF under grant no. IIS-2112862 and by a Cottrell SEED award from Research Corporation.more » « less
-
Abstract The stochastic trajectories of molecules in living cells, as well as the dynamics in many other complex systems, often exhibit memory in their path over long periods of time. In addition, these systems can show dynamic heterogeneities due to which the motion changes along the trajectories. Such effects manifest themselves as spatiotemporal correlations. Despite the broad occurrence of heterogeneous complex systems in nature, their analysis is still quite poorly understood and tools to model them are largely missing. We contribute to tackling this problem by employing an integral representation of Mandelbrot’s fractional Brownian motion that is compliant with varying motion parameters while maintaining long memory. Two types of switching fractional Brownian motion are analysed, with transitions arising from a Markovian stochastic process and scale-free intermittent processes. We obtain simple formulas for classical statistics of the processes, namely the mean squared displacement and the power spectral density. Further, a method to identify switching fractional Brownian motion based on the distribution of displacements is described. A validation of the model is given for experimental measurements of the motion of quantum dots in the cytoplasm of live mammalian cells that were obtained by single-particle tracking.more » « less
-
In vertebrate brains, virtually all neural circuits operate inside a dense matrix of axons (fibers) that have a strongly stochastic character. These fibers originate in the brainstem raphe region, produce highly tortuous trajectories, and release serotonin (5-hydroxytryptamine, 5-HT), with other neurotransmitters. They can robustly regenerate in the adult mammalian brain and appear to support neuroplasticity [1], with implications for mental disorders [2] and artificial neural networks [3]. The self-organization of this “serotonergic” matrix remains poorly understood. In our previous study, we have shown that serotonergic fibers can be modeled as paths of fractional Brownian motion (FBM), a continuous-time stochastic process. FBM is parametrized by the Hurst index, which defines three distinctive regimes: subdiffusion (H < 0.5), normal diffusion (H = 0.5), and superdiffusion (H > 0.5). In two-dimensional (2D) shapes based on the adult mouse brain, simulated FBM-fibers (with H = 0.8) have produced regional distributions similar to those of the actual serotonergic fibers [4]. However, increments of superdiffusive FBM trajectories have long-range positive correlations, which implies that a fiber path in one 2D-section depends on its history in other sections. In a major extension of this study, we used a supercomputing simulation to generate 960 fibers in a complex, three-dimensional shape based on the late-embryonic mouse brain (at embryonic day 17.5). The fibers were modeled as paths of reflected FBM with H = 0.8. The reflection was caused by natural neuroanatomical borders such as the pia and ventricles. The resultant regional densities were compared to the actual fiber densities in the corresponding neuroanatomically-defined regions, based on a recently published comprehensive map [5]. The relative simulated densities showed strong similarities to the actual densities in the telencephalon, diencephalon, and mesencephalon. The current simulation does not include tissue heterogeneities, but it can be further improved with novel models of multifractional FBM, such as the one introduced by our group [6]. The study demonstrates that serotonergic fiber densities can be strongly influenced by the geometry of the brain, with implications for neurodevelopment, neuroplasticity, and brain evolution. Acknowledgements: This research was funded by an NSF-BMBF CRCNS grant (NSF #2112862 to SJ & TV; BMBF #STAXS to RM). References: 1. Lesch KP, Waider J. Serotonin in the modulation of neural plasticity and networks: implications for neurodevelopmental disorders. Neuron. 2012, 76, 175-191. 2. Daws RE, Timmermann C, Giribaldi B, et al. Increased global integration in the brain after psilocybin therapy for depression. Nat. Med. 2022, 28, 844-851. 3. Lee C, Zhang Z, Janušonis S. Brain serotonergic fibers suggest anomalous diffusion-based dropout in artificial neural networks. Front. Neurosci. 2022, 16, 949934. 4. Janušonis S, Detering N, Metzler R, Vojta T. Serotonergic axons as fractional Brownian motion paths: Insights Into the self-organization of regional densities. Front. Comput. Neurosci. 2020, 14, 56. 5. Awasthi JR, Tamada K, Overton ETN, Takumi T. Comprehensive topographical map of the serotonergic fibers in the male mouse brain. J. Comp. Neurol. 2021, 529, 1391-1429. 6. Wang W, Balcerek M, Burnecki K, et al. Memory-multi-fractional Brownian motion with continuous correlation. arXiv. 2023, 2303.01551.more » « less
