TSWJ The Scientific World Journal 1537-744X Hindawi Publishing Corporation 182508 10.1155/2014/182508 182508 Research Article A Directed Continuous Time Random Walk Model with Jump Length Depending on Waiting Time Shi Long 1, 2 Yu Zuguo 1, 3 Mao Zhi 1 Xiao Aiguo 1 Li C. Liu F. Yuste S. B. 1 Hunan Key Laboratory for Computation and Simulation in Science and Engineering and Key Laboratory of Intelligent Computing and Information Processing of Ministry of Education Xiangtan University Xiangtan, Hunan 411105 China xtu.edu.cn 2 Institute of Mathematics and Physics Central South University of Forest and Technology, Changsha Hunan 410004 China cfstujwc.com 3 School of Mathematical Sciences Queensland University of Technology G.P.O. Box 2434, Brisbane, QLD 4001 Australia qut.edu.au 2014 13 03 2014 2014 03 01 2014 10 02 2014 13 3 2014 2014 Copyright © 2014 Long Shi et al. 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.

In continuum one-dimensional space, a coupled directed continuous time random walk model is proposed, where the random walker jumps toward one direction and the waiting time between jumps affects the subsequent jump. In the proposed model, the Laplace-Laplace transform of the probability density function P ( x , t ) of finding the walker at position x at time t is completely determined by the Laplace transform of the probability density function φ ( t ) of the waiting time. In terms of the probability density function of the waiting time in the Laplace domain, the limit distribution of the random process and the corresponding evolving equations are derived.

1. Introduction

The continuous time random walk (CTRW) theory, which was introduced by Montroll and Weiss  to study random walks on a lattice, has been applied successfully in many fields (see, e.g., the reviews  and references therein).

In continuum one-dimensional space, a CTRW process is generated by a sequence of independent identically distributed (IID) positive waiting times T 1 , T 2 , T 3 , , and a sequence of IID random jump lengths X 1 , X 2 , X 3 , . Each waiting time has the same probability density function (PDF) φ ( t ) , t 0 , and each jump length has the same PDF λ ( x ) (usually chosen to be symmetric λ ( x ) = λ ( - x ) ). Setting t 0 = 0 , t n = T 1 + T 2 + + T n for n N and x 0 = 0 , x n = X 1 + X 2 + + X n , x ( t ) = x n for t n t < t n + 1 , we get a microscopic description of the diffusion process . If { X n } and { T n } are independent, the CTRW is called decoupled. Otherwise it is called coupled CTRW . The decoupled CTRW, which is completely determined by mutually independent random jump length and random waiting time, has been widely studied in recent years .

In some applications it becomes important to consider coupled CTRW [7, 8]. The coupled CTRW can be described by the joint PDF ϕ ( x , t ) of jump length and waiting time. Because ϕ ( x , t ) d x d t is the probability of a jump to be in the interval ( x , x + d x ) in the time interval ( t , t + d t ) , the waiting time PDF φ ( t ) = - + ϕ ( x , t ) d x and the jump length PDF λ ( x ) = 0 + ϕ ( x , t ) d t can be deduced. Some kinds of couplings and correlations were proposed in , where the symmetric jump length PDF is chosen. For the coupled CTRW, there exist two coupled forms: ϕ ( x , t ) = λ ( x ) φ ( t x ) and ϕ ( x , t ) = φ ( t ) λ ( x t ) . The first coupled form has been studied sufficiently in many literatures [8, 2123]. The famous model is Lévy walk. Recently, we considered the second coupled form, discussed the asymptotic behaviors of the coupled jump probability density function in the Fourier-Laplace domain, and derived the corresponding fractional diffusion equations from the given asymptotic behaviors .

In this work, we introduce a directed CTRW model with jump length depending on waiting time (i.e., ϕ ( x , t ) = φ ( t ) λ ( x t ) , x > 0 , t > 0 ). In our model, the Laplace-Laplace transform  of    P ( x , t )    of finding the walker at position x at time t is completely determined by the Laplace transform of φ ( t ) . Generally, CTRW processes can be categorised by the mean waiting time T = 0 + t φ ( t ) d t being finite or infinite. Here we find that the long-time limit distributions of the PDF P ( x , t ) are a Dirac delta function for finite T and a beta-like density for infinite T , the corresponding evolving equations are a standard advection equation for finite T and a pseudodifferential equation with fractional power of coupled space and time derivative for infinite T .

This paper is organized as follows. In Section 2, we introduce the basic concepts of the coupled CTRW. In Section 3, a coupled directed CTRW model is introduced. In Section 4, the limit distributions and the corresponding evolving equations of the coupled directed CTRW model are derived. The conclusions are given in Section 5.

2. The Coupled Continuous Time Random Walk

Now we recall briefly the general theory of CTRW . Let η ( x , t ) be the PDF of just having arrived at position x at time t . It can be expressed by η ( x , t ) (the PDF of just having arrived at position x at time t < t ) as (1) η ( x , t ) = - + d x 0 + d t η ( x , t ) ϕ ( x - x , t - t ) + δ ( x ) δ ( t ) . Then, the PDF P ( x , t ) with the initial condition P ( x , 0 ) = δ ( x ) can be described by the following integral equation : (2) P ( x , t ) = 0 t η ( x , t ) ω ( t - t ) d t , where ω ( t ) = 1 - 0 t φ ( τ ) d τ is the probability of not having made a jump until time t .

Let f ^ ( k ) and g ~ ( s ) be the transforms of Fourier and Laplace of sufficiently well-behaved (generalized) functions f ( x ) and g ( t ) , respectively, defined by (3) f ^ ( k ) = { f ( x ) ; k } = - + f ( x ) e i k x d x , k R , g ~ ( s ) = { g ( t ) ; s } = 0 + g ( t ) e - s t d t , s > s 0 .

After using the Fourier-Laplace transforms and the convolution theorems for integral equation (2), one can obtain the following famous algebraic relation : (4) P ~ ^ ( k , s ) = 1 - φ ~ ( s ) s · 1 1 - ϕ ~ ^ ( k , s ) .

3. A Coupled Directed CTRW Model

In , the author considered a CTRW model with waiting time depending on the preceding jump length, where the author supposed that the PDF of the waiting time is a function of a preceding jump length. In that model, the author introduced a natural “physiological” analogy: after making a jump one needs time to rest and recover. The longer the jump distance is, the longer the recovery and the waiting time needed are. This is an interesting hypothetical physiological example. Motivated by this, we consider a directed CTRW model with jump length depending on the waiting time and give an analogue physiological explanation.

A directed CTRW model with jump length depending on the waiting time can be generated by a sequence of IID positive waiting times T 1 , T 2 , T 3 , , and a sequence of jumps X 1 , X 2 , X 3 , ; each waiting time has the same PDF φ ( t ) , t 0 . Every time jump has the same direction and each jump length has the same conditional PDF λ ( x t ) , x 0 , which is the PDF of the random walker making a jump of length x following a waiting time t .

A natural assumption is that the jump length is proportional to the waiting time. So we can take the simplest jump length PDF as λ ( x t ) = δ ( x - v t ) , v > 0 . Without loss of generality, we take v = 1 in the following discussion. Setting t 0 = 0 , t n = T 1 + T 2 + + T n for n N and x 0 = 0 , x n = X 1 + X 2 + + X n , x ( t ) = x n for t n t < t n + 1 , we get a directed CTRW process, where the joint PDF ϕ ( x , t ) can be expressed by ϕ ( x , t ) = φ ( t ) δ ( x - t ) . A physiological explanation can be made as follows: the walker has a random time for a rest to supplement energy and then makes a jump. The longer the rest time is, the longer the jump length can be.

Since the variable x takes positive values in proposed directed CTRW model, it is convenient to replace the Fourier transform for variable x in formula (4) with the Laplace transform (i.e.,   f ~ ( k ) = { f ( x ) ; k } = 0 + f ( x ) e - k x d x ) to obtain the following Laplace-Laplace relation : (5) P ~ ~ ( k , s ) = 1 - φ ~ ( s ) s · 1 1 - ϕ ~ ~ ( k , s ) . Since (6) ϕ ~ ~ ( k , s ) = 0 + d t 0 + ϕ ( x , t ) e - k x - s t d x = 0 + d t 0 + φ ( t ) δ ( x - t ) e - k x - s t d x = 0 + φ ( t ) e - ( s + k ) t d t = φ ~ ( s + k ) , (5) is recast into (7) P ~ ~ ( k , s ) = 1 - φ ~ ( s ) s · 1 1 - φ ~ ( s + k ) .

The    n th ( n = 1,2 ) moment of P ( x , t ) is given by (8) x n ( t ) = 0 + x n ( t ) P ( x , t ) d x = ( - 1 ) n n k n P ~ ( k , t ) | k = 0 = - 1 { 1 - φ ~ ( s ) s · ( - 1 ) n n k n 1 1 - φ ~ ( s + k ) | k = 0 } .

In the following section, we will study the possible behaviors of P ( x , t ) and its n th ( n = 1,2 ) moment.

4. The Limit Distributions of the Coupled Directed CTRW Model

From (7), we can see that the Laplace-Laplace transform of PDF P ( x , t ) is completely determined by the Laplace transform of the waiting time PDF φ ( t ) . Usually, the random waiting time is characterized by its mean value T . It may be finite or infinite.

For finite mean waiting time T , the Laplace transform of φ ( t ) is of the form (9) φ ~ ( s ) = 1 - s T + o ( s ) , s 0 .

Substituting (9) into (7), in the limit ( k , s ) ( 0,0 ) , we get the asymptotic relation (10) P ~ ~ ( k , s ) ~ 1 - ( 1 - s T ) s · 1 1 - ( 1 - ( s + k ) T ) = 1 s + k .

After taking the inverse Laplace transforms for (10) about k and s , we have (11) P ( x , t ) = δ ( x - t ) .

For long times (12) x ( t ) = t , x 2 ( t ) = t 2 .

From (10), we get (13) s P ~ ~ ( k , s ) - 1 + k P ~ ~ ( k , s ) = 0 .

Using { P ( x , t ) / t ; s } = s P ~ ( x , s ) - P ( x , 0 ) , { P ( x , t ) / x ; k } = k P ~ ( k , t ) - P ( 0 , t ) , initial condition P ( x , 0 ) = δ ( x ) , and natural boundary conditions, we obtain the partial differential equation (14) P ( x , t ) t + P ( x , t ) x = 0 , which is the standard advection equation.

In many applications, one needs to consider a long waiting time (i.e., T is infinite); it is natural to generalize (9) to the following form: (15) φ ~ ( s ) = 1 - s β + o ( s β ) , s 0 , 0 < β 1 .

Inserting (15) into (7), in the limit ( k , s ) ( 0,0 ) , we get the asymptotic relation (16) P ~ ~ ( k , s ) ~ 1 - ( 1 - s β ) s · 1 1 - ( 1 - ( s + k ) β ) = s β - 1 ( s + k ) β .

After taking the Laplace inverse transform for (16) about s , one has (17) P ~ ( k , t ) = t - β Γ ( 1 - β ) * [ e - k t t β - 1 Γ ( β ) ] = 0 t e - k τ τ β - 1 ( t - τ ) - β Γ ( β ) Γ ( 1 - β ) d τ , where we use the formulas { t β - 1 ; s } = Γ ( β ) / s β for β > 0 , { e - a t g ( t ) ; s } = g ~ ( s + a ) , and { ( f * g ) ( t ) ; s } = f ~ ( s ) g ~ ( s ) .

According to formula (8) and (17), for long times, one gets (18) x ( t ) = β t , x 2 ( t ) = β ( β + 1 ) 2 t 2 .

Then taking the Laplace inverse transform for (17) about k , the following form is obtained: (19) P ( x , t ) = 0 t δ ( x - τ ) τ β - 1 ( t - τ ) - β Γ ( β ) Γ ( 1 - β ) d τ = x β - 1 ( t - x ) - β Γ ( β ) Γ ( 1 - β ) , which is the density of a random variable t B , where B has a Beta distribution with parameters β and 1 - β .

From (17), we can also obtain (20) ( s + k ) β P ~ ~ ( k , s ) = s β - 1 , which leads to the pseudodifferential equation [27, 28] (21) ( t + x ) β P ( x , t ) = δ ( x ) t - β Γ ( 1 - β ) with a coupled space-time fractional derivative operator on the left-hand side.

Equation (21) is useful to model flow in porous media and other physical systems characterized by a link between the waiting time and the jump length.

5. Conclusions

In this work, we introduce a directed CTRW model with jump lengths depending on waiting times. By the Laplace-Laplace transform technique, we find that the PDF P ( x , t ) is determined only by the waiting times PDF φ ( t ) . For finite and infinite mean waiting time, we deduce the limit distributions of P ( x , t ) from the asymptotic behaviors of φ ( t ) in the Laplace domain, respectively. The corresponding evolving equations are also derived. For finite mean waiting time, the limit behavior of the PDF P ( x , t ) is governed by a standard advection equation. For infinite mean waiting time, the limit behavior of the PDF P ( x , t ) is governed by a pseudodifferential equation with coupled space-time fractional derivative. We also calculate the first-order moment x ( t ) and the second-order moment x 2 ( t ) of P ( x , t ) . An interesting phenomenon is obtained: there exist the relations x ( t ) ~ t , x 2 ( t ) ~ t 2 , whether the mean waiting time is finite or not.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This project was supported by the Natural Science Foundation of China (Grant nos. 11371016 and 11271311), the Chinese Program for Changjiang Scholars and Innovative Research Team in University (PCSIRT) (Grant no. IRT1179), the Research Foundation of Education Commission of Hunan Province of China (Grant no. 11A122), and the Lotus Scholars Program of Hunan province of China.

Montroll E. W. Weiss G. H. Random walks on lattices. II Journal of Mathematical Physics 1965 6 2 167 181 2-s2.0-0002567920 Bouchaud J.-P. Georges A. Anomalous diffusion in disordered media: statistical mechanisms, models and physical applications Physics Report 1990 195 4-5 127 293 2-s2.0-0040307478 Metzler R. Klafter J. The random walk's guide to anomalous diffusion: a fractional dynamics approach Physics Report 2000 339 1 1 77 2-s2.0-0002641421 Metzler R. Klafter J. The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics Journal of Physics A 2004 37 31 R161 R208 2-s2.0-4043151477 10.1088/0305-4470/37/31/R01 Gorenflo R. Vivoli A. Mainardi F. Discrete and continuous random walk models for space-time fractional diffusion Nonlinear Dynamics 2004 38 1–4 101 116 2-s2.0-15544368331 10.1007/s11071-004-3749-5 Meerschaert M. M. Nane E. Xiao Y. Correlated continuous time random walks Statistics and Probability Letters 2009 79 9 1194 1202 2-s2.0-63549130217 10.1016/j.spl.2009.01.007 Shlesinger M. F. Klafter J. Wong Y. M. Random walks with infinite spatial and temporal moments Journal of Statistical Physics 1982 27 3 499 512 2-s2.0-0001479044 10.1007/BF01011089 Klafter J. Blumen A. Shlesinger M. F. Stochastic pathway to anomalous diffusion Physical Review A 1987 35 7 3081 3085 2-s2.0-0000550285 10.1103/PhysRevA.35.3081 Weiss G. H. Aspects and Applications of the Random Walk. 1994 Amsterdam, The Netherlands North Holland Roman H. E. Alemany P. A. Continuous-time random walks and the fractional diffusion equation Journal of Physics A 1994 27 10, article 017 3407 3410 2-s2.0-21344498117 10.1088/0305-4470/27/10/017 Hilfer R. Anton L. Fractional master equations and fractal time random walks Physical Review E 1995 51 2 R848 R851 2-s2.0-0000437484 10.1103/PhysRevE.51.R848 Scalas E. Gorenflo R. Mainardi F. Fractional calculus and continuous-time finance Physica A 2000 284 1 376 384 2-s2.0-0034275979 10.1016/S0378-4371(00)00255-7 Barkai E. CTRW pathways to the fractional diffusion equation Chemical Physics 2002 284 1-2 13 27 2-s2.0-0036828743 10.1016/S0301-0104(02)00533-5 Hilfer R. On fractional diffusion and continuous time random walks Physica A 2003 329 1-2 35 40 2-s2.0-0141955949 10.1016/S0378-4371(03)00583-1 Scalas E. Gorenflo R. Mainardi F. Uncoupled continuous-time random walks: solution and limiting behavior of the master equation Physical Review E 2004 69 1 2-s2.0-1842531858 011107 Scalas E. The application of continuous-time random walks in finance and economics Physica A 2006 362 2 225 239 2-s2.0-32544456755 10.1016/j.physa.2005.11.024 Gorenflo R. Mainardi F. Vivoli A. Continuous-time random walk and parametric subordination in fractional diffusion Chaos, Solitons and Fractals 2007 34 1 87 103 2-s2.0-34047155650 10.1016/j.chaos.2007.01.052 Chechkin A. V. Hofmann M. Sokolov I. M. Continuous-time random walk with correlated waiting times Physical Review E 2009 80 3 2-s2.0-70349768154 10.1103/PhysRevE.80.031112 031112 Tejedor V. Metzler R. Anomalous diffusion in correlated continuous time random walks Journal of Physics A 2010 43 8 2-s2.0-76349090018 10.1088/1751-8113/43/8/082002 082002 Fa K. S. Uncoupled continuous-time random walk: finite jump length probability density function Journal of Physics A 2012 45 19 2-s2.0-84860452557 10.1088/1751-8113/45/19/195002 195002 Blumen A. Zumofen G. Klafter J. Transport aspects in anomalous diffusion: Lévy walks Physical Review A 1989 40 7 3964 3973 2-s2.0-0000914304 10.1103/PhysRevA.40.3964 Zumofen G. Klafter J. Scale-invariant motion in intermittent chaotic systems Physical Review E 1993 47 2 851 863 2-s2.0-33750026674 10.1103/PhysRevE.47.851 Zaburdaev V. Y. Random walk model with waiting times depending on the preceding jump length Journal of Statistical Physics 2006 123 4 871 881 2-s2.0-33747654448 10.1007/s10955-006-9104-0 Liu J. Bao J. D. Continuous time random walk with jump length correlated with waiting time Physica A 2013 392 612 617 Shi L. Yu Z. Mao Z. Xiao A. Huang H. Space-time fractional diffusion equations and asymptotic behaviors of a coupled continuous time random walk model Physica A 2013 392 5801 5807 Goreno R. Mainardi F. Rogosin S. Laplace-Laplace analysis of the fractional poisson process Analytical Methods of Analysis and Differential Equations 2012 43 58 Becher-Kern P. Meerschaert M. M. Scheffer H. P. Limit theorems for coupled continuous time random walks The Annals of Probability 2004 32 730 756 Jurlewicz A. Kern P. Meerschaert M. M. Scheffer H. P. Fractional governing equations for coupled random walks Computers & Mathematics with Applications 2012 64 3021 3036