Investigating the interplay between mechanisms of anomalous diffusion via fractional Brownian walks on a comb-like structure
The comb model is a simplified description for anomalous diffusion under geometric constraints. It represents particles spreading out in a two-dimensional space where the motions in the -direction are allowed only when the coordinate of the particle is zero. Here, we propose an extension for the comb model via Langevin-like equations driven by fractional Gaussian noises (long-range correlated). By carrying out computer simulations, we show that the correlations in the direction affect the diffusive behavior in the direction in a non-trivial fashion, resulting in a quite rich diffusive scenario characterized by usual, superdiffusive or subdiffusive scaling of second moment in the direction. We further show that the long-range correlations affect the probability distribution of the particle positions in the direction, making their tails longer when noise in the direction is persistent and shorter for anti-persistent noise. Our model thus combines and allows the study/analysis of the interplay between different mechanisms of anomalous diffusion (geometric constraints and long-range correlations) and may find direct applications for describing diffusion in complex systems such as living cells.
pacs:05.40.-a, 87.10.Mn, 89.75.-k, 05.40.Fb
Keywords: anomalous diffusion, comb model, fractional Brownian motion, long-range correlations
Understanding diffusive motions is a long-standing problem in the physicists’ agenda. Since the works of Einstein, Smoluchowski and Langevin, we know that one of the most striking patterns of classical (usual) free diffusion is the linear growth with time of the second moment of the particle positions. While this feature is ubiquitous in nature, many other patterns have been also observed for the second moment, both from analytical calculations and experimental data . These deviants behaviors are often called anomalous diffusion.
Nowadays, efforts have been mainly focused on understanding the underlying physical mechanisms that lead to such deviations, since only the classification of the diffusive behaviors of a given system is not enough to fully characterize it. In this context, two common pathways for anomalous diffusion are long-range correlations in the particle positions and geometric constraints related with the structural complexity of the environment where the particles are moving. The former may be related to memory affects and an example is the fractional Brownian motion [2, 3], where the second moment has a power-law dependence with time, that is, . For this case, the Hurst exponent classifies the anomalous diffusion in subdiffusion () or superdiffusion () [4, 5] and also recovers usual () and ballistic () diffusions as limiting cases. Anomalous diffusion due to geometric constraints is well exemplified by the comb model. In this model, diffusive particles walk on a two-dimensional space; however, motions in the -direction are only allowed when the coordinate of the particle positions is zero. The result of this simple constraint is a subdiffusive motion (precisely, with ) and the appearance of a backbone (at ) and teeth (along the axis direction), which were originally proposed to mimic the quasilinear structure and dangling ends of percolation clusters [6, 7, 8].
Mainly due to its simplicity and ability of mimicking diffusive aspects of highly disordered systems, the comb model has been extensively studied and extended by means of analytical calculations [9, 10, 11, 12, 13, 14, 15] and also employed as simplified description of natural phenomena such as cancer proliferation [16, 17], transport of spiny dendrites  and diffusion of ultra cold atoms [19, 20]. Moreover, as crowded environments tend to slow down the diffusive motion of particles by obstruction and/or trapping , the comb model could also be used as a toy model for describing diffusion in intracellular processes, where several researchers have reported the existence of subdiffusion and ergodicity breaking in the random walk of different biochemical compounds [22, 23, 24, 25, 27, 26, 28, 21, 29, 30, 31, 32] by using single particle tracking . In the case of intracellular processes, the power-law exponents describing the second moment ranges from to (evaluated both from ensemble and time average) and, in the case of ultra cold atoms, appears as a function of lattice depth and assumes values larger than one (, superdiffusion).
In addition to the anomalous scaling of the second moment, the work of Weiss  has also shown that the diffusive processes of intracellular fluids may present long-range corrections as in the fractional Brownian motion, an empirical finding that have not been explored within the comb model framework. Also, recently, Yamamoto et al.  have provided evidence (from molecular dynamics simulation) that the anomalous diffusion of water molecules on cell membrane surface arises from both divergent mean trapping time and long-range correlated noises. Here, we extend the comb model by considering a Langevin-like equation where the noises are long-range correlated. Our model thus combines two of the main mechanisms of anomalous diffusion: geometric constraints and long-range correlations, allowing a more complete description of complex diffusive motions and also the study of the interplay between these two mechanisms. By carrying out computer simulations, we show that the correlations in the direction affect the diffusive behavior in the direction in a non-trivial fashion and that a quite rich diffusive scenario emerges from these long-range correlations. We report that diffusive process in the direction of this generalized version of the comb model can be usual, superdiffusive or subdiffusive, depending on both corrections in the and directions. We further show that the long-range correlations affect the probability distribution of the particle positions in the direction: the distribution becomes more leptokurtic when the noise in the direction is persistent whereas, for anti-persistent noise, it becomes less leptokurtic and approaches the Gaussian distribution as a limiting case.
2 Generalized comb and long-range correlations
We start by writing the generalized diffusion equation describing the motion under the comb structure :
here is the joint probability of the particle positions ( and coordinates) as a function of time , and are the diffusive coefficients in the and directions. Note the presence of the Dirac delta multiplying the spatial derivative with respect to and, consequently, limiting the diffusion in the direction only for . In order to investigate the role of long-range correlations in this diffusive process, we propose to numerically solve the following coupled Langevin-like equations
Notice that the noises and are power-law correlated with scaling exponents and , where and are the so-called Hurst exponents. For , the correlation functions have negative signs and the noises are anti-persistent, meaning that positive values are followed by negative values (or vice-versa) more frequently than by chance; while for the correlations are positive, meaning that positive values are followed by positive values and negative values are followed by negative values more frequently than by chance. Fully persistent noises occur as a limiting case when , whereas normal Brownian (white) noises correspond to the limit of .
When running out the simulations of the Langevin equations, we have considered a narrow band of thickness along the axis, inside which the diffusion in the direction occurs. This strip mimics the effect of the Dirac delta (that appears multiplying the noise in Eq. 2) and we have verified that the value has no influence in the diffusive process, as long as and noises amplitudes and are of the same order of magnitude. This restriction ensures that the variable presents no dynamics inside the strip. For our proposes, we have fixed . The fractional Gaussian noises and were generated by following the Hosking method . Figure 1 shows examples of simulated trajectories for three different values of and keeping . By visual inspection, we note that depending on whether the motion in is anti-persistent [Fig 1(a)] or persistent [Fig 1(c)], the diffusive behavior in the direction drastically changes. The anti-persistence in causes the particle to return more often to the strip where the diffusive motion in the direction is allowed; consequently, covers a larger interval along the axis when compared with the usual comb [, Fig 1(b)]. On the other hand, persistence in causes the particle to return less often to the strip and, consequently, covers a small interval along the axis.
Before quantifing the effects of the long-range correlations on the variances of and , we have first verified whether the Langevin equations [Eq. (2)] are actually equivalent to the generalized diffusion equation [Eq. (1)] when the noises and are uncorrelated (). We thus simulate an ensemble of particle positions considering and . By using these data, we evaluate the temporal evolution of the variances
where denotes ensemble average. Figures 2(a) and (b) show the behavior of and versus time for several values of the maximum integration time in the Langevin equations, ranging from to as indicated in these plots. In both cases, we note a finite-size-like effect characterized by a crossover time where and change their behaviors. Despite that and also discarding an initial transient regime (), we observe that the variance profiles of and are in well agreement with the power laws predicted by the generalized diffusion equation of the usual comb [Eq. (1)], that is, and with and [36, 38, 39]. We have also analyzed the behavior of the crossover time in function of the maximum integration time . Figures 2(c) and (d) show that grows linearly with for and , confirming that the simulated power-law regimes prevail when tends to infinity.
In order to further strengthen the equivalence between our Langevin equations (2) and the diffusion equation (1), we have also evaluated the marginal probability functions and from the simulated data. Figures 3(a) and 3(b) show these distributions for several values of where we note the spreading of these distributions. For a better comparison with the analytical predictions of Eq. (1), we have also calculated the marginal distributions considering the normalized positions and . Figures 3(c) and 3(d) reveal a good collapse of the empirical distributions and also a good agreement with the analytical distributions (dashed lines) predicted by Eq. (1), which is a standard Gaussian for and can be expressed in terms of the Fox H function for [36, 38, 39] or by the following summation [5, 30]
Thus, for the case where no correlations are present in the noises and , the results shown in Fig. 3 corroborate with the equivalence between our simulations of the Langevin equations and the diffusion equation for the comb model.
Convinced of the equivalence in the usual case, we now address the role of the long-range correlations on the diffusive properties of a random walker described by our Langevin equations (2). We first investigate how the memory effects in the direction modify the diffusion in the and directions. In order to do so, we have built an ensemble of particle positions considering and for each one of nineteen values of homogeneously distributed between and . For each value of , we evaluate the temporal evolution of the variances and . Figures 4(a) and 4(b) show the variance profiles in the and directions, respectively. In these plots, the dashed lines represent the behavior for the usual comb, that is, when no correlations are present. In the case of , we note that the more the noise in is anti-persistent (), the larger is the power-law exponent describing the behavior of ; in contrast, persistent noise in slows down the diffusion and decreases the power-law exponent . The extreme cases where and are characterized by and , respectively. Thus, as we have first noticed in Fig. 1, anti-persistent noise in forces the walker to return more frequently to the strip where it can move in the direction and, consequently, the diffusion is enhanced when compared with the case usual comb (). On the other hand, persistent noise in drives the walker away from strip and, consequently, the diffusion slows down. In the case of , the behavior is straightforward: anti-persistent noise in slows down the diffusion while persistent noise enhances it [see Fig. 4(b)]. In order to quantify the role of long-range correlations in on and , we have evaluated the dependence of the power-law exponents and on . In Fig. 4(c), we observe a monotonic but nonlinear decay of as increases, while in Fig. 4(d), we confirm the straightforward behavior of as function of ().
Another interesting question is whether the combinations of long-range correlated noises in the and directions affect the diffusive behaviors in a nontrivial way. To answer this question, we have built an ensemble of particle positions for each possible pair of and homogeneously distributed between and with step size of . For each combination of and , we calculate the profile of the variances and and also the power-law exponents and describing the main tendency of these curves (by using the same procedure of Fig. 4). Figure 5(a) shows a contour plot of as a function of and and Fig. 5(b) represents the same for . For the Fig. 5(a), we note a quite rich diffusive scenery where, depending on the values of and , we may have subdiffusion slower () or faster () than the usual comb, usual diffusion () and even superdiffusion (). We further observe from this figure that cannot be written as a linear combination of and (which would be represented by straight lines in this plot) and thus, the correlations in and present a nontrivial interplay in the behavior of . On the other hand, the case of Fig 5(b) is rather simple because the correlations in do not affect the diffusive behavior in the direction, this becomes evident by noting that the level curves of are horizontal straight lines.
It is also important to investigate whether the long-range correlations affect the distribution profiles of the particle positions. In order to do so, we first calculate the marginal distributions and for there different values of (, and ) and keeping . Figures 6(a) and 6(b) show these distributions. We note that the distribution displays tails longer than those of usual comb (dark shaded region) when the noise in the direction is persistent () whereas, for anti-persistent noise (), the tails are shorter than those of usual comb and close to the Gaussian distribution (light shaded region). The distributions of the normalized positions are not affected by long-range correlated noises in the direction. We thus confirm that the long-range correlations in have influence on the profiles of . For a better and quantitative characterization of the role of the correlated noise on the shape of , we evaluate the coefficient of kurtosis in function of and . For the distribution is leptokurtic (peaked) while for the distribution is platykurtic (flat) also, for a Gaussian distribution and for the usual comb (distribution of the Eq. 5). Figure 6(c) shows a contour plot of as a function of and , where we observe that the kurtosis ranges from very close to ( for ) up to almost ( for ). It is worth noting that, similar to the variance case [Fig. 5(a)], the kurtosis is not a trivial combination of and .
3 Summary and Conclusions
We have proposed an extension for the comb model via Langevin equations driven by long-range correlated noises (fractional Gaussian noises). We initially showed that our Langevin equations are equivalent to the generalized diffusion equation (describing the usual comb) through the comparison of the numerical-obtained variances and marginal distributions with the analytical results predicted by the diffusion equation. Next, we have presented an extensive characterization of the diffusive properties of the particle positions in function of the Hurst exponents and defining the correlated noises. Our results show that the noise in the direction affects the variances and marginal distributions along the direction. Specifically, long-range persistence in slows down the diffusion in by reducing the power-law exponent and makes the tails of the marginal distribution longer as increases; whereas long-range anti-persistence in enhances the diffusion in by increasing the value of and makes the tails of shorter as decreases. We have also investigated the interplay between the long-range correlations in the and directions, where we found that both and the kurtosis coefficient of are not a trivial (such as a linear combination) function of and , but instead our results suggest a more complex coupling between the Hurst exponents and the diffusion in the direction. In summary, a quite rich diffusive scenery emerges from our model which may be directly applied for describing (at least in a first approximation) some of the recent empirical findings related to diffusion in complex systems such as living cells. Furthermore, our model combines and allows the study/analysis of the interplay between different mechanisms of anomalous diffusion, geometric constraints and long-range correlations, an important fact that had not been explored within the comb model framework and that may provide a better understanding of these processes as a combination of different mechanisms.
We are all grateful to Capes, CNPq and Fundação Araucária for financial support.
- : New J. Phys.
- Klafter J and Sokolov I M 2005 Phys. World 18 29
- Mandelbrot B B and van Ness J W 1968 SIAM Rev. 10 422
- Gripenberg G and Norros I 1996 J. Appl. Probab. 33 400
- Bouchad J-P and Georges A 1990 Phys. Rep. 195 127
- Metzler R and Klafter J 2000 Phys. Rep. 339 1
- White S R and Barma M 1984 J. Phys. A 17 2995
- Ben-Avraham D and Havlin S 2000 Diffusion and Reactions in Fractals and Disordered Systems (Cambridge University Press, Cambridge)
- Weiss G H and Havlin S 1986 Physica A 134 474
- Baskin E and Iomin A 2004 Phys. Rev. Lett. 93 120603
- Arkhincheev V E 2007 Chaos 17 043102
- Iomin A 2011 Phys. Rev. E 83 052106
- Tateishi A A, Lenzi E K, Ribeiro H V, Evangelista L R, Mendes R S and da Silva LR 2011 J. Stat. Mech P02022
- Villamaina D, Sarracino A, Gradenigo G, Puglisi A and Vulpiani A 2011 J. Stat. Mech L01002
- Tateishi A A, Michels F S, dos Santos M A F, Lenzi E K and Ribeiro H V 2013 J. Stat. Mech P09017
- Lenzi E K, da Silva L R, Tateishi A A, Lenzi M K and Ribeiro H V 2013 Phys. Rev. E 87 012121
- Iomin A 2006 Phys. Rev. E 73 061918
- Iomin A 2005 J. Phys.: Conf. Ser. 7 57
- Mendez V and Iomin A 2013 Chaos Solitons Fract. 53 46
- Iomin A 2012 Phys. Rev. E 86 032101
- Sagi Y, Brook M, Almog I and Davidson N (2012) Phys. Rev. Lett. 108 093002
- Höfling F and Franosch T 2013 Rep. Prog. Phys. 76 046602
- Condamin S, Tejedor V, Voituriez R, Bénichou O and Klafter J 2008 Proc. Natl. Acad. Sci. USA 105 5675
- Selhuber-Unkel C, Yde P, Berg-Sørensen K and Oddershede L B 2009 Phys. Biol. 6 025015
- Jeon J-H, Tejedor V, Burov S, Barkai E, Selhuber-Unkel C, Berg-Sørensen K, Oddershede L and Metzler R 2011 Phys. Rev. Lett. 106 048103
- Burov S, Jeon J-H, Metzler R, and Barkai E 2011 Phys. Chem. Chem. Phys. 13 1800
- Sokolov I M 2012 Soft Matter 8 9043
- Leijnse N, Jeon J-H, Loft S, Metzler R and Oddershede L B 2012 Eur. Phys. J. Special Topics 204 75
- Cherstvy A G, Chechkin A V and Metzler R 2013 New J. Phys. 15 083039
- Brackley C A, Cates M E and Marenduzzo D 2013 Phys. Rev. Lett. 111 108101
- Girst S, Hable V, Drexler G A, Greubel C, Siebenwirth C, Haum M, Friedl A A and Dollinger G 2013 Sci. Rep. 3 2511
- Tabei S M A, Burov S, Kim H Y, Kuznetsov A, Huynh T, Jureller J, Philipson L H, Dinner A R and Scherer N F 2013 Proc. Natl. Acad. Sci. USA 110 4911
- Cherstvy A G, Chechkin A V and Metzler R 2014 Soft Matter 10 1591
- Braüchle C, Lamb C D and Michaelis J (eds.) 2010 Single Particle Tracking and Single Molecule Energy Transfer (Wiley-VCH, Weinheim)
- Weiss M 2013 Phys. Rev. E 88 010101(R)
- Yamamoto E, Akimoto T, Yasui M and Yasuoka K (2014) Sci. Rep. 4 4720
- Arkhincheev V E and Baskin E M 1991 Zh. Eksp. Teor. Fiz. 100 292 [1991 Sov. Phys. JETP 73 161].
- Hosking J R M 1984 Water Resour. Res. 20 1898
- Arkhincheev V E 1999 JETP 88 710
- da Silva L R, Tateishi A A, Lenzi M K, Lenzi E K and da Silva P C 2009 Braz. J. Phys. 39 483