Fractional Brownian motion in a finite interval: correlations effect depletion or accretion zones of particles near boundaries
T Guggenberger,G Pagnini,T Vojta,R Metzler
DOI: https://doi.org/10.1088/1367-2630/ab075f
2019-02-28
New Journal of Physics
Abstract:Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) with 1 α α α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers.Diffusive transport is quite ubiquitous, ranging from quantum processes such as laser cooling over thermally activated transport in living biological cells, to the dispersal of tracer chemicals in geophysical aquifers. While theoretical works on diffusion often consider infinite or semi-infinite domains, in many cases the particle motion is restricted to a finite interval. Such a scenario is relevant, inter alia, for the dispersal of light in disordered optical cavities, for the motion of thermally driven particles in confining microfluidic chambers, gels, or for the molecular and (sub)micron tracer motion in the confines of biological cells or their organelles. For normal Brownian diffusion in an interval of length L the probability density to find the particle anywhere within this interval has the constant amplitude 1/L at sufficiently long times: in the stationary state the particle can be found everywhere equally likely. This property is so engrained in our intuition for random processes that we typically would not question its validity. Indeed, even for generalised random processes such as the continuous time random walk, such equidistributions naturally occur [1–4]. We here show that for the widely used, Gaussian process of fractional Brownian motion (FBM) describing overdamped viscoelastic diffusion, this property is strikingly violated: the inherent negative or positive correlations of FBM in a finite interval effect a pronounced depletion or accretion of the probability density in the vicinity of the boundaries, respectively4.While normal diffusion is characterised by the linear time dependence of the mean squared displacement (MSD)5, anomalous diffusion typically classifies the power law time dependenceof the MSD, where the generalised diffusion coefficient Kα has physical dimension of length2/timeα, and α is the anomalous diffusion exponent. Depending on the value of α one distinguishes between subdiffusion (0 α α > 1) [2, 4].Anomalous diffusion has been observed experimentally and numerically in a wide array of systems [2, 4]. Thus, experimentally subdiffusion is frequently observed for passive particle motion inside living biological cells [5–11] and in crowded, viscoelastic solutions [6, 12–14]. It was also reported from supercomputing studies of lipid and protein molecule motion in bilayer membranes [15–18]. Subdiffusion moreover occurs in other systems, such as for charge carrier motion in amorphous semiconductors [1] or chemical tracers in underground aquifers [19]. Superdiffusion has been observed for the motion of molecular motor-transported particles inside biological cells [20–22], the motion of tracers in two-dimensional rotating flows [23], and for bulk-mediated surface diffusion at liquid–solid interfaces [24].Normal Brownian motion is characterised by the universal Gaussian probability density [2]. Anomalous diffusion processes lose this universality, and the various possibilities to break the attraction of the basin of the central limit theorem give rise to different mathematical models [25]. For instance, the existence of power-law sojourn times between motion events with a diverging mean waiting time, as indeed measured -Abstract Truncated-
physics, multidisciplinary