Jae-Hyung
Jeon
a,
Aleksei V.
Chechkin
bcd and
Ralf
Metzler
*ad
aDepartment of Physics, Tampere University of Technology, FI-33101 Tampere, Finland
bInstitute for Theoretical Physics, Kharkov Institute of Physics and Technology, Kharkov 61108, Ukraine
cMax-Planck Institute for the Physics of Complex Systems, Nöthnitzer Straße 38, 01187 Dresden, Germany
dInstitute of Physics & Astronomy, University of Potsdam, 14776 Potsdam-Golm, Germany. E-mail: rmetzler@uni-potsdam.de
First published on 18th June 2014
Anomalous diffusion is frequently described by scaled Brownian motion (SBM), a Gaussian process with a power-law time dependent diffusion coefficient. Its mean squared displacement is 〈x2(t)〉 ≃ 2(t)t with (t) ≃ tα−1 for 0 < α < 2. SBM may provide a seemingly adequate description in the case of unbounded diffusion, for which its probability density function coincides with that of fractional Brownian motion. Here we show that free SBM is weakly non-ergodic but does not exhibit a significant amplitude scatter of the time averaged mean squared displacement. More severely, we demonstrate that under confinement, the dynamics encoded by SBM is fundamentally different from both fractional Brownian motion and continuous time random walks. SBM is highly non-stationary and cannot provide a physical description for particles in a thermalised stationary system. Our findings have direct impact on the modelling of single particle tracking experiments, in particular, under confinement inside cellular compartments or when optical tweezers tracking methods are used.
〈x2(t)〉 ≃ 2Kαtα | (1) |
In biological contexts one often invokes the concept of apparent subdiffusion, which is in fact a transient crossover between free normal diffusion with α = 1 and the thermal plateau of the MSD in confinement.6 However, there is clear experimental evidence of long-ranged subdiffusion (1) in living biological cells4,7–12 and in artificially crowded or structured liquids.13–17 The physical mechanisms for such anomalous diffusion are very varied, even considering the most common approaches: we mention the continuous time random walk approach,3,5,18 in which subdiffusion is effected by multiple trapping or sticking events with an emerging power-law distribution of waiting times, as seen in the data of ref. 9, 10 and 12. Alternatively, subdiffusion arises in fractal geometries with their dead ends and bottlenecks across spatial scales,19 an additional mechanism pinpointed in ref. 12. We also mention diffusion processes in which the anomaly arises from the spatial dependence of the diffusion coefficients,20–22 or from interaction of the tracer particle with a structured environment.23–26
A special role play the Gaussian models for anomalous diffusion which belong to two families, that can best be distinguished on the Langevin equation level. The Mandelbrot–van Ness fractional Brownian motion (FBM) follows the stochastic equation ẋ(t) = ζfGn(t), where ζfGn(t) represents fractional Gaussian noise characterised by the covariance 〈ζfGn(t1)ζfGn(t2)〉 ≃ α(α − 1)|t1 − t2|α−2 which has long-ranged negative or positive correlations, respectively, for subdiffusion (0 < α < 1) or superdiffusion (1 < α < 2).27,28 FBM and the closely related fractional Langevin equation motion are physically related to the motion in a viscoelastic environment.29 In various experiments this type of motion was identified as single or partial component.9–11,15,30
Here we focus on the second family of Gaussian anomalous diffusion models, namely, scaled Brownian motion (SBM).31 SBM is governed by the Langevin equation
(2) |
(t) = αKαtα−1, | (3) |
For unconfined SBM the probability density function is the Gaussian31
(4) |
(5) |
For SBM the mean of the time averaged MSD over multiple trajectories i, can be derived from the Langevin eqn (2), yielding48
(6) |
(7) |
(8) |
In Fig. 1 we show results from simulations of the SBM process for both sub- and superdiffusion, observing excellent agreement with the analytical findings. We also see that the scatter between individual trajectories is very small. Such scatter is a characteristic for anomalous diffusion models and can be used, e.g., to reliably distinguish FBM from subdiffusive CTRW processes.40,46 For CTRW subdiffusion we find pronounced scatter between from individual trajectories even in the limit of extremely long trajectories t → ∞,5,40,41 while for FBM the scatter vanishes for longer t, a signature of the ergodic nature of FBM.29,40,46,50 The amplitude scatter of individual can be characterised in terms of the dimensionless variable . If ξ has a narrow distribution around ξ = 1 and the width decreases with increasing t, the process is usually considered ergodic. This width is characterised in terms of the ergodicity breaking parameter EB = 〈ξ2〉 − 〈ξ〉2. For SBM it was found49
(9) |
Fig. 1 MSD and time averaged MSD for SBM with α = 1/2 (top) and α = 3/2 (bottom). The simulations of the SBM-Langevin eqn (2) show excellent agreement with the analytical results, eqn (1) and (6). We also show results for the time averaged MSD for 20 individual trajectories. Apart from the region where the lag time Δ approaches the length t of the time series and statistics worsen, there is hardly any amplitude scatter between individual . |
From the analysis so far for free motion we can see that the SBM process has a truly split personality. Thus, on the one hand SBM's PDF is identical to that for free FBM. On the other hand, however, in contrast to the ergodic behaviour of FBM for sufficiently long t,40,50 SBM exhibits weak ergodicity breaking , as demonstrated by comparison of eqn (7) with eqn (1). This scaling form exactly matches the result for CTRW subdiffusion5,40,41 or diffusion processes with space-dependent diffusivity.20,21 Unlike the weakly non-ergodic dynamics of the latter two, for SBM the fluctuations around the mean measured by the distribution ϕ(ξ) are narrow and decrease with longer t. We now show that the behaviour of confined SBM is also unconventional.
(10) |
〈x2(t)〉 = 2Kαtαe−2ktM(α,1 + α,2kt), | (11) |
〈x2(t)〉 ∼ αKαk−1tα−1 | (12) |
The result for the time averaged MSD is similarly remarkable. As shown in Fig. 2 for SBM simulations based on the Langevin eqn (2) with the Hookean forcing, it exhibits a pronounced apparent plateau for lag times Δ ≫ 1/k for both sub- and superdiffusive SBM. This behaviour is in excellent agreement with the full analytical solution (16) provided in the Appendix in terms of Kummer functions. Taking the limit Δ ≪ t we obtain
(13) |
Fig. 3 analyses this behaviour in the harmonic potential further by showing the time series x(t) for subdiffusive SBM, Brownian motion, and superdiffusive SBM. Indeed we see that for subdiffusive and superdiffusive confined SBM the fluctuations continue to decrease and increase with time, while those in the Brownian limit become stationary. This is the direct effect of the time dependent temperature or viscosity encoded in the SBM diffusivity (t). If the system is stationary or thermalised this behaviour clearly underlines the unsuitability of SBM for the description of anomalous diffusion.
The behaviour of SBM dynamics under confinement is the central result of our study. The continued temporal decay or increase of the MSD that we obtained for SBM is in stark contrast to the behaviour of confined CTRW subdiffusion, for which the MSD 〈x2〉 saturates to the thermal plateau, while the time averaged MSD continues to grow in the power-law form up until the lag time Δ approaches t.40,43 It strongly differs from FBM, which relaxes to a plateau for both the MSD and the time averaged MSD, and for which only a transient disparity between 〈x2〉 and exists.15 Finally, SBM is also at variance with heterogeneous diffusion processes with a space dependent diffusion coefficient that relax to a plateau for both 〈x2〉 and .53
For free SBM the probability density function P(x,t) equals that of FBM, despite the fact that both processes are governed by different stochastic (Langevin) equations. SBM's time averaged MSD scales equally to those of the weakly non-ergodic CTRW subdiffusion and diffusion processes with space-dependent diffusivity. Despite this weakly non-ergodic character of the mean time averaged MSD , the amplitude scatter between the time averaged MSD of individual realisations is small and the distribution is bell shaped, as otherwise observed for the ergodic FBM or for Brownian motion of finite time t. In that sense SBM represents a new class of non-ergodic processes. The most striking behaviour of SBM is, however, its strongly non-stationary behaviour under confinement. Instead of relaxing to a plateau the MSD acquires a power-law decay or growth mirroring a continuously decreasing or increasing temperature encoded in SBM's time dependent diffusion coefficient (t).
SBM is thus at variance with the currently available experimental observations in complex liquids using single particle tracking by video microscopy or by optical tweezers tracking of single submicron particles. Thus the free anomalous diffusion data garnered so far was classified into FBM-like and CTRW-like behaviour, or combinations thereof.8,10,12 Note that also for fluorescence correlation spectroscopy experiments recent analysis tools corroborated an FBM nature of the data.13 For optical tweezers tracking of lipid granules in different complex liquids the subdiffusive time averaged MSD either continues to grow under confinement, reflecting the non-ergodic features of CTRWs,5,9 or it relaxes to a plateau value mirroring an ergodic dynamics.15,54
How can FBM and SBM have the same distribution (4) in free space? Simply put for FBM, the viscoelastic properties of the environment effecting the long-range correlations of FBM lead to a frequency dependent response of the environment to a disturbance, while the materials properties remain unchanged in time. For free FBM this gives rise to the subdiffusive MSD (1), in which the time dependent diffusivity effectively encodes the frequency dependent response of the viscoelastic environment. The distribution (4) for FBM and its description in terms of a Fokker–Planck equation of the form (10) is treacherous, however. This can already be seen when we use the PDF (4) or the dynamics eqn (10) to calculate the first passage behaviour. This procedure leads to the wrong result,55 and the full analytical description of FBM in the presence of boundaries remains elusive, a difficulty imposed by the highly correlated fractional Gaussian noise driving its Langevin equation. In contrast, SBM, according to its Langevin description (2), is driven by uncorrelated noise but the environment itself is changing as function of time, effecting an extremely non-stationary process. The equivalence of the PDF (4) of both processes is thus simply due to the fact that a Gaussian PDF is completely defined by its second moment, the MSD (1).
With its interesting behaviour that is so different from the other conventional anomalous diffusion models, SBM may indeed have relevant applications to weakly coupled or fully adiabatic systems, as well as for active systems, in which the existence of a temperature is not meaningful. In particular, in the superdiffusive case SBM or analogous dynamics with other increasing effective diffusivities may represent an alternative approach to active Brownian motion.56
The difference of SBM to other processes can also be seen in Fig. 4. For the subdiffusive case a sample trajectory of SBM is compared to that of FBM and the noisy CTRW,47 in which the pure subdiffusive CTRW is superimposed with Ornstein–Uhlenbeck noise to accommodate the thermal noise of the environment observed in experiments.25 We see that SBM with its Gaussian probability density function and uncorrelated driving noise appears more similar to the noisy CTRW motion, albeit it has a more pronounced tendency to reach larger amplitudes than the CTRW with its waiting time periods. The fluctuations of SBM are dramatically less than those of the anticorrelated FBM, which also shows the largest amplitudes in x(t). For superdiffusion, we compare SBM with a noisy version of the Lévy walk process57 and again with FBM. Note the vastly different size of the window shown on the ordinate. This time, the persistent FBM shows pronouncedly larger excursions and a distinct reduction of the noise compared to SBM. The shape of x(t) of the latter does not appear to be qualitatively different from the subdiffusive case. The noisy CTRW is fundamentally different from both SBM and FBM. Despite their disparate physical nature and their dissimilarity when setting one against the other in a direct comparison, we note that generally it is difficult to identify a stochastic process solely from the appearance of the recorded trajectory.
Fig. 4 Individual trajectories of anomalous stochastic processes: FBM, SBM, and noisy CTRW (nCTRW), for subdiffusion with α = 1/2 (top) and superdiffusion with α = 3/2 (right). |
From this discussion it is obvious that one thing is to have at our disposal a handy and easy-to-use description for anomalous diffusion processes: SBM with its Gaussian and uncorrelated nature appears deceivingly simple and is therefore easy to implement in numerical analyses and descriptions such as diffusion limited reactions. However, it is yet another thing to consider the physical relevance of a model process. Using SBM with its time-dependent diffusion coefficient violates the physical setting in typical experiments, in which the system is held at approximately constant temperature and its predictions are at odds with actual observations. The unphysical nature for the kind of processes we have in mind is most obvious for confined motion. In contrast, FBM and fractional Langevin equation motion are ergodic processes with the physical background of an effective particle motion in a viscoelastic multi-body environment. FBM and fractional Langevin equation motion exhibit a transient disparity between the MSD and the time averaged MSD under confinement.15 Weakly non-ergodic CTRW processes emerge due to immobilisation periods that are imprinted on the dynamics by the structure of the environment or binding events to the environment. Finally, diffusion processes with space dependent diffusivities arise naturally in non-homogeneous systems such as biological cells or subsurface aquifers. To extract physically meaningful information from anomalous diffusion data one needs to have some physical insight into the observed process and properly analyse the data using complementary tools5,10–12,15,47,58–60 before settling for the appropriate physical model.
We finally note that it will be interesting to compare the predictions of the SBM model in the superdiffusive range 1 < α < 2 with that of active processes in viscoelastic environments, which in the relevant overdamped limit are of FBM type.61,62
〈x(t1)x(t2)〉 = 2Kαtα1e−k(t1+t2)M(α,1 + α,2kt1) | (14) |
(15) |
(16) |
(17) |
In Fig. 5 we illustrate the convergence of the time averaged MSD to the value of the MSD 〈x2(t)〉 in the limit Δ → t.
Fig. 5 Convergence of the time averaged MSD to the MSD 〈x2(Δ)〉 when the lag time Δ approaches the measurement time t. The remarkable behaviour is due to the pole in the defining expression (5). We show the result for α = 1/2 and α = 3/2. The measurement time t = 105. Top: free SBM. For the superdiffusive case the dashed black line of unity slope aids in demonstrating the non-linear behaviour of . Bottom: SBM in an harmonic potential with k = 0.1. Note the different ranges of the abscissa in the two cases. For the extremely small window for the subdiffusive case needed to illustrate the cusp at Δ → t and the relatively large variation of the MSD 〈x2(Δ)〉 appears almost constant. |
In Fig. 6 we demonstrate the bell shaped scatter around the ergodic value ξ = 1 of the relative amplitude . For smaller values of Δ/t, the scatter clearly decreases and is almost Gaussian. At the larger Δ value the distribution is somewhat more asymmetric. The dashed lines in both panels of Fig. 6 indicate the Gaussian shape calculated for the ergodic FBM process in ref. 46. Clearly, the distribution ϕ(ξ) for SBM is broader.
Fig. 6 Scatter PDF ϕ(ξ) of the relative amplitude around the ergodic value ξ = 1, for Δ = 1 (top) and Δ = 10 (bottom), as well as α = 1/2. |
We acknowledge funding from the Academy of Finland within the FiDiPro scheme.
This journal is © the Owner Societies 2014 |