Research Article | Open Access

Long Shi, Aiguo Xiao, "Modeling Anomalous Diffusion by a Subordinated Integrated Brownian Motion", *Advances in Mathematical Physics*, vol. 2017, Article ID 7246865, 7 pages, 2017. https://doi.org/10.1155/2017/7246865

# Modeling Anomalous Diffusion by a Subordinated Integrated Brownian Motion

**Academic Editor:**Giorgio Kaniadakis

#### Abstract

We consider a particular type of continuous time random walk where the jump lengths between subsequent waiting times are correlated. In a continuum limit, the process can be defined by an integrated Brownian motion subordinated by an inverse -stable subordinator. We compute the mean square displacement of the proposed process and show that the process exhibits subdiffusion when , normal diffusion when , and superdiffusion when . The time-averaged mean square displacement is also employed to show weak ergodicity breaking occurring in the proposed process. An extension to the fractional case is also considered.

#### 1. Introduction

Anomalous diffusion is found in a wide diversity of systems (see review articles [1–4] and references therein). In one dimension, it is characterized by a mean square displacement (MSD) of the formwith , which deviates from the linear dependence on time found in normal diffusion. The coefficient is generalized diffusion constant. It is called subdiffusion for and superdiffusion for [2].

A fundamental account to anomalous diffusion is provided by a stochastic process called continuous time random walk (CTRW), which was originally introduced by Montroll and Weiss in 1965 [5]. In a continuum limit, the process has been considered by Fogedby [6] via coupled Langevin equationswhere is a white Gaussian noise with , , and is a white -stable Lévy noise, taking positive values only and independent of .

In (2), the random walk is parametrized in terms of a continuous variable , which is subjected to a random time change. This random time change to the physical time is described by the second equation. The combined process in the physical time is then given by , where is the inverse process to defined as

Mathematically, the fundamental approach to describe the combined process is based on subordination technique, which was first introduced by Bochner [7]. Using the notation of subordination, the process , , and are named parent process, subordinator, and inverse subordinator, respectively.

In recent years, (2) consisting of Brownian motions with or without external field and inverse -stable subordinator are becoming a hot topic [8–17]. There are also several other processes considered as parent processes within the subordination framework, for example, Lévy-stable process [18, 19], arithmetic Brownian motion [20], geometric Brownian motion [21, 22], Ornstein-Uhlenbeck process [23, 24], tempered stable process [25], fractional Brownian motion [26, 27], and fractional Lévy-stable process [28]. Here, we note that, apart from inverse -stable subordinator, inverse tempered -stable subordinator and infinitely divisible subordinators are also considered in the literatures [16, 20, 25–31].

In the simplest CTRW process, after each jumps, a new pair of waiting time and jump length is drawn from the associated distributions, independent of the previous values. This independence giving rise to a renewal process is not always justified, for instance, by observations of human motion patterns [32] and active biological movements [33] or in financial market dynamics [34]. Recently, three correlated CTRW models are introduced to model the random walks with some forms of memory [35–37]. Some advances in the field of CTRWs with correlated temporal or/and spatial structure can be also found in [38–45].

In this work, we consider a jump-correlated CTRW model which has the subordination form . Here, the parent process is an integrated Brownian motion, defined byand inverse subordinator is the inverse of one-side -stable Lévy process , defined by

The integrated Brownian motion is called the random acceleration process in the physical literature and has been studied by many authors. For instance, it appears in the continuum description of the equilibrium Boltzmann weight of a semiflexible polymer chain [46]. It also appears in the description of statistical properties of the Burgers equation with Brownian initial velocity [47]. Some further results of the integrated Brownian motion can be found in the paper [48] reviewing this subject.

The structure of the paper is as follows. In Section 2, we introduce the jump-correlated CTRW model. In Section 3, we compute MSD of the proposed process and observe the corresponding anomalous diffusive behaviors. The time-averaged MSD is also employed to show weak ergodicity breaking occurring in the proposed process. In Section 4, we generalize the integrated Brownian motion to the fractional integral of Brownian motion and compute the corresponding MSD. The conclusions are given in Section 5.

#### 2. Model

We begin by recalling the general framework for CTRW theory. Let be the sequence of nonnegative independent identically distributed (IID) random variable representing waiting times between jumps of a particle. We set and , that is, the time of the th jump. Let be the sequence of IID jump lengths of the particle, which are assumed to be independent of waiting times. We set and , that is, the position of the particle after the th jump. Then, the position of the particle at time is given bywhere is the number of jumps up to time . The process is called CTRW.

In what follows, we analyze a particular type of CTRW where the jump lengths are correlated. Assume that each jump is equal towhere are IID random variables with finite second moment (for simplicity we assume that their second moment is equal to 1). Moreover, we assume that each waiting time is nonnegative IID random variable, whose characteristic function is given by

In the continuous limit, we get the following set of coupled Langevin equations for the position and time of the CTRWwhere and are the same as those in (2) and is the standard Brownian motion with , .

An equivalent representation of (9) in the form of subordination isHere the parent process has the formand the inverse subordinator is defined bywhere is an -stable totally skewed Lévy motion with characteristic function

#### 3. Discussions

At first, let us compute the MSD of subordinated process .

Assume that , , and are PDFs of subordinated process , parent process , and inverse subordinator , respectively. In terms of subordination, we have

Since the first moment of parent process and the second momentwe obtain

Thus, the MSD of the subordinated process is

Let us turn to the inverse subordinator . Observing the equivalence from (12) we obtain the relationwhich gives the formula for the PDF in terms of the PDF :Taking the Laplace transform for (22) about variable , we get

Thus, the MSD of the subordinated process in Laplace space iswhich implies that the MSD of is

It is easy to observe from (25) that the process is subdiffusive when , normally diffusive when , and superdiffusive when .

It is well-known that the MSD of the process given by (2) is of the formComparing (25) with (26), we see that Fogedby’s model can only represent anomalous subdiffusion, but our model can represent subdiffusion, normal diffusion, and superdiffusion.

Next, we study weak ergodicity breaking of the subordinated process .

In an ergodic system, one can find the equivalenceHere, is the time-averaged MSD of the process , defined aswhere is the lag time and is the overall measure time.

For anomalous diffusion, the behavior of the ensemble MSD and the time-averaged MSD (28) may be fundamentally different. The disparity is usually called weak ergodicity breaking (or weak nonergodicity) [49]. In recent years, weak nonergodicity of anomalous diffusion process attracts more and more attentions [49–55].

Since, for any , parent process satisfieswhere means an equality in distribution, we have

Thus,where .

In the limit ,

Hence,

Sincecomparing (33) with (34), we see that the linear lag time dependence of is different from the power-law form of , which implies that subordinated process is weakly nonergodic.

At last, we consider the propagator associated with the subordinated process . By the total probability formula, we obtain an integral representation of :

For fixed , the random variable is normally distributed. From (15) and (16), we have

It follows from and the Laplace transform for that we obtainAfter taking the inverse Laplace transform for , we getwhere is the Mittag-Leffler function with parameter [56].

#### 4. An Extension to the Fractional Case

In this section, we introduce the dependent sequence of jump lengths in the following manner:where is a memory function. The continuous limit is of the form

Integrating (42) we get

After taking () and using the integration by parts, (43) can be written aswhere is the Riemann-Liouville fractional integration operator of order , defined by [56]

As a result, the jump-correlated CTRW has the subordination form , where parent process is of the form (44), and inverse subordinator is defined by (12).

Here, we are interested in the competition between the memory parameter and stability index . In what follows, we will not discuss any properties of motion other than the MSD.

In terms of (44), we getBy denoting and exchanging the order of quadratic integral , we obtainThus,where and is Beta function.

We observe from (48) that, in the limiting case , memory function , the parent process defined by (43) reduces to the form defined by (11), and the second moment of computed by (48) reduces to , where , the same form as (16). We are also interested in the limiting case . At the moment, the memory function is a Dirac -function; defined by (43) reduces to the standard Brownian motion.

Let us turn to the MSD of the subordinated process . In terms of (18) and (48), we obtainIn the Laplace space, the MSD is of the formTaking the inverse Laplace transform for , we havewhere . In the limiting case , the parameter reduces to Thus, (51) reduces to (25).

It is easy to observe from (51) that there exists a competition between the memory parameter and stability index . For the case , the subordinated process exhibits subdiffusive behaviors. For the case , the process is subdiffusive when , normal diffusive when , and superdiffusive when .

#### 5. Conclusions

We introduce an integrated Brownian motion subordinated by inverse -stable one-sided Lévy motion, which is a continuous limit of a jump-correlated CTRW. In terms of the ensemble MSD of the proposed process, we conclude that the process is subdiffusive when , normal diffusive when , and superdiffusive when . The time-averaged MSD is also employed to show weak ergodicity breaking occurring in the proposed process.

We also generalize the process to the case, where the parent process is fractional integral of Brownian motion. In terms of the MSD, we observe a competition between the memory parameter and stability index . Other types of inverse subordinators may be also candidates.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant no. 11671343) and the Scientific Research Project of Hunan Provincial Education Department (no. 17B258).

#### References

- J.-P. Bouchaud and A. Georges, “Anomalous diffusion in disordered media: statistical mechanisms, models and physical applications,”
*Physics Reports*, vol. 195, no. 4-5, pp. 127–293, 1990. View at: Publisher Site | Google Scholar | MathSciNet - R. Metzler and J. Klafter, “The random walk's guide to anomalous diffusion: a fractional dynamics approach,”
*Physics Reports*, vol. 339, no. 1, pp. 1–77, 2000. View at: Publisher Site | Google Scholar | MathSciNet - R. Metzler and J. Klafter, “The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics,”
*Journal of Physics. A. Mathematical and General*, vol. 37, no. 31, pp. R161–R208, 2004. View at: Publisher Site | Google Scholar | MathSciNet - I. Eliazar and J. Klafter, “Anomalous is ubiquitous,”
*Annals of Physics*, vol. 326, no. 9, pp. 2517–2531, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH - E. W. Montroll and G. H. Weiss, “Random walks on lattices. II,”
*Journal of Mathematical Physics*, vol. 6, pp. 167–181, 1965. View at: Publisher Site | Google Scholar | MathSciNet - H. C. Fogedby, “Langevin equations for continuous time Lévy flights,”
*Physical Review E*, vol. 50, no. 2, pp. 1657–1660, 1994. View at: Publisher Site | Google Scholar - S. Bochner, “Diffusion equation and stochastic processes,”
*Proceedings of the National Academy of Sciences*, vol. 35, no. 7, pp. 368–370, 1949. View at: Publisher Site | Google Scholar - A. Baule and R. Friedrich, “Joint probability distributions for a class of non-Markovian processes,”
*Physical Review E*, vol. 71, no. 2, Article ID 026101, 2005. View at: Publisher Site | Google Scholar - A. Piryatinska, A. I. Saichev, and W. A. Woyczynski, “Models of anomalous diffusion: the subdiffusive case,”
*Physica A: Statistical Mechanics and Its Applications*, vol. 349, no. 3-4, pp. 375–420, 2005. View at: Publisher Site | Google Scholar - M. Magdziarz, A. Weron, and K. Weron, “Fractional Fokker-Planck dynamics: stochastic representation and computer simulation,”
*Physical Review E—Statistical, Nonlinear, and Soft Matter Physics*, vol. 75, no. 1, Article ID 016708, 2007. View at: Publisher Site | Google Scholar - D. Kleinhans and R. Friedrich, “Continuous-time random walks: simulation of continuous trajectories,”
*Physical Review E*, vol. 76, no. 6, Article ID 061102, 2007. View at: Publisher Site | Google Scholar - A. Weron, M. Magdziarz, and K. Weron, “Modeling of subdiffusion in space-time-dependent force fields beyond the fractional Fokker-Planck equation,”
*Physical Review E*, vol. 77, no. 3, Article ID 036704, 2008. View at: Publisher Site | Google Scholar - M. Magdziarz, A. Weron, and J. Klafter, “Equivalence of the fractional fokker-planck and subordinated langevin equations: the case of a time-dependent force,”
*Physical Review Letters*, vol. 101, no. 21, Article ID 210601, 2008. View at: Publisher Site | Google Scholar - A. Weron and S. Orzel, “Itô formula for subordinated Langevin equation. A case of time dependent force,”
*Acta Physica Polonica B*, vol. 40, no. 5, pp. 1271–1277, 2009. View at: Google Scholar - S. Eule and R. Friedrich, “Subordinated Langevin equations for anomalous diffusion in external potentials—biasing and decoupled external forces,”
*EPL*, vol. 86, no. 3, Article ID 30008, 2009. View at: Publisher Site | Google Scholar - M. Magdziarz, “Langevin picture of subdiffusion with infinitely divisible waiting times,”
*Journal of Statistical Physics*, vol. 135, no. 4, pp. 763–772, 2009. View at: Publisher Site | Google Scholar - M. Magdziarz, “Stochastic representation of subdiffusion processes with time-dependent drift,”
*Stochastic Processes and Their Applications*, vol. 119, no. 10, pp. 3238–3252, 2009. View at: Publisher Site | Google Scholar | MathSciNet - M. Magdziarz and A. Weron, “Competition between subdiffusion and Lévy flights: a Monte Carlo approach,”
*Physical Review E*, vol. 75, Article ID 056702, 2007. View at: Publisher Site | Google Scholar - B. o. Dybiec and E. Gudowska-Nowak, “Subordinated diffusion and continuous time random walk asymptotics,”
*Chaos*, vol. 20, no. 4, Article ID 043129, 2010. View at: Publisher Site | Google Scholar | MathSciNet - A. Wyłomańska, “Arithmetic Brownian motion subordinated by tempered stable and inverse tempered stable processes,”
*Physica A: Statistical Mechanics and Its Applications*, vol. 391, no. 22, pp. 5685–5696, 2012. View at: Publisher Site | Google Scholar - H. Gu, J.-R. Liang, and Y.-X. Zhang, “On a time-changed geometric Brownian motion and its application in financial market,”
*Acta Physica Polonica B*, vol. 43, no. 8, pp. 1667–1681, 2012. View at: Publisher Site | Google Scholar | MathSciNet - J. Gajda and A. Wyłomańska, “Geometric Brownian motion with tempered stable waiting times,”
*Journal of Statistical Physics*, vol. 148, no. 2, pp. 296–305, 2012. View at: Publisher Site | Google Scholar - J. Janczura, S. Orzeł, and A. Wyłomańska, “Subordinated
*αα*-stable Ornstein-Uhlenbeck process as a tool for financial data description,”*Physica A: Statistical Mechanics and its Applications*, vol. 390, no. 23-24, pp. 4379–4387, 2011. View at: Publisher Site | Google Scholar - J. Gajda and A. Wyłomańska, “Time-changed Ornstein–Uhlenbeck process,”
*Journal of Physics A: Mathematical and Theoretical*, vol. 48, no. 13, Article ID 135004, 2015. View at: Publisher Site | Google Scholar - A. Wyłomańska, “The tempered stable process with infinitely divisible inverse subordinators,”
*Journal of Statistical Mechanics: Theory and Experiment*, vol. 2013, no. 10, Article ID P10011, 2013. View at: Publisher Site | Google Scholar - Y.-X. Zhang, H. Gu, and J.-R. Liang, “Fokker-planck type equations associated with subordinated processes controlled by tempered
*α*-stable processes,”*Journal of Statistical Physics*, vol. 152, no. 4, pp. 742–752, 2013. View at: Publisher Site | Google Scholar | MathSciNet - J. Gajda and A. Wyłomańska, “Fokker–Planck type equations associated with fractional Brownian motion controlled by infinitely divisible processes,”
*Physica A: Statistical Mechanics and Its Applications*, vol. 405, pp. 104–113, 2014. View at: Publisher Site | Google Scholar - M. Teuerle, A. Wyłomańska, and G. Sikora, “Modeling anomalous diffusion by a subordinated fractional Lévy-stable process,”
*Journal of Statistical Mechanics: Theory and Experiment*, vol. 2013, no. 5, Article ID P05016, 2013. View at: Publisher Site | Google Scholar - J. Gajda and M. Magdziarz, “Fractional Fokker-Planck equation with tempered
*α*-stable waiting times: Langevin picture and computer simulation,”*Physical Review E. Statistical, Nonlinear, and Soft Matter Physics*, vol. 82, no. 1, Article ID 011117, 2010. View at: Publisher Site | Google Scholar | MathSciNet - J. Janczura and A. Wyłomańska, “Anomalous diffusion models: different types of subordinator distribution,”
*Acta Physica Polonica B*, vol. 43, no. 5, pp. 1001–1016, 2012. View at: Publisher Site | Google Scholar - J. Gajda, “Fractional Fokker–Planck equation with space dependent drift and diffusion: the case of tempered
*α*-stable waiting-times,”*Jagellonian University. Institute of Physics. Acta Physica Polonica B*, vol. 44, no. 5, pp. 1149–1161, 2013. View at: Publisher Site | Google Scholar | MathSciNet - C. Song, T. Koren, P. Wang, and A. Barabási, “Modelling the scaling properties of human mobility,”
*Nature Physics*, vol. 6, no. 10, pp. 818–823, 2010. View at: Publisher Site | Google Scholar - A. Maye, C.-H. Hsieh, G. Sugihara, and B. Brembs, “Order in spontaneous behavior,”
*PLoS ONE*, vol. 2, no. 5, article e443, 2007. View at: Publisher Site | Google Scholar - E. Scalas, “The application of continuous-time random walks in finance and economics,”
*Physica A: Statistical Mechanics and its Applications*, vol. 362, no. 2, pp. 225–239, 2006. View at: Publisher Site | Google Scholar - A. Chechkin, M. Hofmann, and I. M. Sokolov, “Continuous-time random walk with correlated waiting times,”
*Physical Review E. Statistical, Nonlinear, and Soft Matter Physics*, vol. 80, no. 3, Article ID 031112, 2009. View at: Publisher Site | Google Scholar | MathSciNet - M. M. Meerschaert, E. Nane, and Y. Xiao, “Correlated continuous time random walks,”
*Statistics & Probability Letters*, vol. 79, no. 9, pp. 1194–1202, 2009. View at: Publisher Site | Google Scholar | MathSciNet - V. Tejedor and R. Metzler, “Anomalous diffusion in correlated continuous time random walks,”
*Journal of Physics. A. Mathematical and Theoretical*, vol. 43, no. 8, Article ID 082002, 2010. View at: Publisher Site | Google Scholar | MathSciNet - M. Magdziarz, R. Metzler, W. Szczotka, and P. Zebrowski, “Correlated continuous-time random walks in external force fields,”
*Physical Review E*, vol. 85, no. 5, Article ID 051103, 2012. View at: Publisher Site | Google Scholar - M. Magdziarz, R. Metzler, W. Szczotka, and P. Zebrowski, “Correlated continuous-time random walks—scaling limits and Langevin picture,”
*Journal of Statistical Mechanics: Theory and Experiment*, vol. 2012, no. 4, Article ID P04010, 2012. View at: Publisher Site | Google Scholar - M. Magdziarz, W. Szczotka, and P. Zebrowski, “Asymptotic behaviour of random walks with correlated temporal structure,”
*Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences*, vol. 469, no. 2159, Article ID 20130419, 2013. View at: Publisher Site | Google Scholar - J. H. Schulz, A. V. Chechkin, and R. Metzler, “Correlated continuous time random walks: combining scale-invariance with long-range memory for spatial and temporal dynamics,”
*Journal of Physics. A. Mathematical and Theoretical*, vol. 46, no. 47, Article ID 475001, 2013. View at: Publisher Site | Google Scholar | MathSciNet - L. Shi, Z. Yu, H. Huang, Z. Mao, and A. Xiao, “The subordinated processes controlled by a family of subordinators and corresponding Fokker–Planck type equations,”
*Journal of Statistical Mechanics: Theory and Experiment*, vol. 2014, no. 12, Article ID P12002, 2014. View at: Publisher Site | Google Scholar - J. Wang, J. Zhou, L. Lv, W. Qiu, and F. Ren, “Heterogeneous memorized continuous time random walks in an external force fields,”
*Journal of Statistical Physics*, vol. 156, no. 6, pp. 1111–1124, 2014. View at: Publisher Site | Google Scholar - F.-Y. Ren, J. Wang, L.-J. Lv, H. Pan, and W.-Y. Qiu, “Effect of different waiting time processes with memory to anomalous diffusion dynamics in an external force fields,”
*Physica A: Statistical Mechanics and Its Applications*, vol. 417, pp. 202–214, 2015. View at: Publisher Site | Google Scholar - L. Lv, F.-Y. Ren, J. Wang, and J. Xiao, “Correlated continuous time random walk with time averaged waiting time,”
*Physica A. Statistical Mechanics and Its Applications*, vol. 422, pp. 101–106, 2015. View at: Publisher Site | Google Scholar | MathSciNet - T. W. Burkhardt, “Semiflexible polymer in the half plane and statistics of the integral of a Brownian curve,”
*Journal of Physics A: Mathematical and General*, vol. 26, no. 22, pp. L1157–L1162, 1993. View at: Publisher Site | Google Scholar - P. Valageas, “Statistical properties of the burgers equation with brownian initial velocity,”
*Journal of Statistical Physics*, vol. 134, no. 3, pp. 589–640, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH - T. W. Burkhardt, “The random acceleration process in bounded geometries,”
*Journal of Statistical Mechanics: Theory and Experiment*, vol. 2007, no. 7, Article ID P07004, 2007. View at: Publisher Site | Google Scholar - J. Jeon, A. V. Chechkin, and R. Metzler, “Scaled Brownian motion: a paradoxical process with a time dependent diffusivity for the description of anomalous diffusion,”
*Physical Chemistry Chemical Physics*, vol. 16, no. 30, pp. 15811–15817, 2014. View at: Publisher Site | Google Scholar - Y. He, S. Burov, R. Metzler, and E. Barkai, “Random time-scale invariant diffusion and transport coefficients,”
*Physical Review Letters*, vol. 101, no. 5, Article ID 058101, 2008. View at: Publisher Site | Google Scholar - W. Deng and E. Barkai, “Ergodic properties of fractional Brownian-Langevin motion,”
*Physical Review E*, vol. 79, no. 1, Article ID 011112, 2009. View at: Publisher Site | Google Scholar - J.-H. Jeon and R. Metzler, “Fractional Brownian motion and motion governed by the fractional Langevin equation in confined geometries,”
*Physical Review E*, vol. 81, no. 2, Article ID 021103, 2010. View at: Publisher Site | Google Scholar | MathSciNet - E. Barkai, Y. Garini, and R. Metzler, “Strange kinetics of single molecules in living cells,”
*Physics Today*, vol. 65, no. 8, pp. 29–35, 2012. View at: Publisher Site | Google Scholar - J. Kursawe, J. Schulz, and R. Metzler, “Transient aging in fractional Brownian and Langevin-equation motion,”
*Physical Review E*, vol. 88, no. 6, Article ID 062124, 2013. View at: Publisher Site | Google Scholar - R. Metzler, J.-H. Jeon, A. G. Cherstvy, and E. Barkai, “Anomalous diffusion models and their properties: non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking,”
*Physical Chemistry Chemical Physics*, vol. 16, no. 44, pp. 24128–24164, 2014. View at: Publisher Site | Google Scholar - I. Podlubny,
*Fractional Differential Equations*, Academic Press, San Diego, Calif, USA, 1999.

#### Copyright

Copyright © 2017 Long Shi and Aiguo Xiao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.