# A Directed Continuous Time Random Walk Model with Jump Length Depending on Waiting Time

Authors: Long Shi; Zuguo Yu; Zhi Mao; Aiguo Xiao
Journal: The Scientific World Journal
 (2014)

---

## Abstract

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 functionP
(
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.


---

## Body

## 1. Introduction

The continuous time random walk (CTRW) theory, which was introduced by Montroll and Weiss [1] to study random walks on a lattice, has been applied successfully in many fields (see, e.g., the reviews [2–4] and references therein).In continuum one-dimensional space, a CTRW process is generated by a sequence of independent identically distributed (IID) positive waiting timesT
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 [5]. If {
X
n
} and {
T
n
} are independent, the CTRW is called decoupled. Otherwise it is called coupled CTRW [6]. The decoupled CTRW, which is completely determined by mutually independent random jump length and random waiting time, has been widely studied in recent years [3–20].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 [21–25], 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, 21–23]. 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 [25].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 [26] 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 Section2, 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 [3]. 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 [3]:
(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.Letf
^
(
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 [3]:
(4)
P
~
^
(
k
,
s
)
=
1
-
φ
~
(
s
)
s
·
1
1
-
ϕ
~
^
(
k
,
s
)
.


## 3. A Coupled Directed CTRW Model

In [23], 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 timesT
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 variablex 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 [26]:
(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
)
.Thenth (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 ofP
(
x
,
t
) and its nth (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 timeT, 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 PDFP
(
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.



---

*Source: 182508-2014-03-13.xml*