TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 351, Number 10, Pages 3915–3934 S 0002-9947(99)02421-6 Article electronically published on May 21, 1999
RATES OF CONVERGENCE OF DIFFUSIONS WITH DRIFTED BROWNIAN POTENTIALS
YUEYUN HU, ZHAN SHI, AND MARC YOR
Abstract. We are interested in the asymptotic behaviour of a diffusion pro- cess with drifted Brownian potential. The model is a continuous time analogue to the random walk in random environment studied in the classical paper of Kesten, Kozlov, and Spitzer. We not only recover the convergence of the dif- fusion process which was previously established by Kawazu and Tanaka, but also obtain all the possible convergence rates. An interesting feature of our ap- proach is that it shows a clear relationship between drifted Brownian potentials and Bessel processes.
1. Introduction Modelling a class of random phenomena in physics and biology (cf. Bouchaud et al. [4]), the simple Random Walk in Random Environment (RWRE) is defined as def follows: Let Ξ = ξi i Z be a sequence of independent and identically distributed random variables{ taking} ∈ values in (0,1), and any realization of Ξ is labelled as “environment”. Define the RWRE Sn n 0 by S0 =0and { } ≥ ξi, if j = i +1, P Sn+1 = j Sn = i;Ξ = 1 ξ, if j = i 1, − i − 0, otherwise, for any n 1andi Z. Among the best known works on RWRE, Solomon [30] treats the≥ basic recurrence/transience∈ problem, Sinai [29] and Kesten et al. [21] es- tablish limit theorems respectively in recurrent and transient situations, Deheuvels and R´ev´esz [9] investigate the almost sure behaviour in Sinai’s setting, and Cs¨org˝o et al. [8] study the local time. We refer to the books of R´ev´esz [25] and Hughes [16] for a detailed account. There is also a renewed interest concerning various aspects of RWRE. Let us mention recent papers by Alili [1], Comets et al. [7], Dembo et al. [10], Gantert and Zeitouni [12], Greven and den Hollander [13], Hu and Shi [14]–[15], and Shi [28]. In the recurrent case treated in the celebrated paper of Sinai [29], it is shown that RWRE moves far slower than the usual random walk (in non-random environment). On the other hand, the convergence and rates of RWRE in the transient case are completely characterized by Kesten et al. [21]. In this paper, we are interested in the continuous time analogue of RWRE, namely, a diffusion process with random potential, introduced by Schumacher [27]
Received by the editors November 17, 1997 and, in revised form, July 3, 1998. 1991 Mathematics Subject Classification. Primary 60J60, 60F05. Key words and phrases. Diffusion with random potential, Bessel process, rate of convergence.
c 1999 American Mathematical Society
3915
License or copyright restrictions may apply to redistribution; see https://www.ams.org/journal-terms-of-use 3916 YUEYUN HU, ZHAN SHI, AND MARC YOR
and Brox [5]. In what follows, W (x); x R denotes a two-sided Brownian motion with W (0) = 0, and { ∈ } κ (1.1) V (x) def= W (x) x, x R, − 2 ∈ is a Brownian motion with drift ( κ/2). Consider the formal solution of the fol- lowing stochastic differential equation:− 1 dX(t)=dβ(t) V0(X(t)) d t, (1.2) −2 X(0) = 0, where β denotes a Brownian motion starting from 0, independent of W (hence of the random potential V ). We note that changing κ into ( κ) leads to the same formal equation for ( X), since β(t); t 0 and W ( −x); x R are independent Brownian motions.− Hence,{− we may,≥ without} loss{ of− generality,∈ } restrict ourselves to κ 0. In fact, we shall only discuss the case κ>0, which characterizes the transience≥ of X (see Brox [5] for the recurrence with κ = 0, and Kawazu and Tanaka [19] for the transience with κ>0). Rigorously speaking, instead of writing the formal derivative of V in (1.2), we should consider X as a diffusion process with generator 1 d d (1.3) eV (x) e V (x) . 2 d x − d x A representation of X(t); t 0 using time change is given in Section 3. We refer to Kawazu and Tanaka{ [19],≥ Mathieu} [23], and Carmona [6] for related studies and further extensions. This paper aims at completing previous work of Kawazu and Tanaka [19], [20], which respectively obtain the convergence of X, and its convergence rate in the particular case κ>2. Our main result is a complete characterization of all possible convergence rates, therefore providing the continuous time analogue of what Kesten et al. [21] established for transient RWRE. Our approach is new, mainly based on Bessel process techniques by showing how Bessel and Jacobi processes can be used in the study of random potentials. To state our result, we define
(1.4) H(r) def=inf t>0: X(t)>r ,r>0, { } the first hitting time process associated with X. Throughout the paper, stands ca N for a Gaussian (0, 1) variable, α (for 0 <α<2) a completely asymmetric N caS stable variable of index α,and λ (for λ>0) a (positive) completely asymmetric Cauchy variable of parameter Cλ (the superscript “ca” standing for “completely ca asymmetric”). When 0 <α<1, α is a positive variable. The distributions of ca and ca are characterized by theirS respective characteristic functions: for t R, Sα Cλ ∈ πα E exp(it ca)=exp t α 1 i sgn(t)tan , Sα −| | − 2 h 2 i E exp(it ca)=exp λ t +it log t . Cλ − || π | | h i The symbols “ law ”, “ P. ”and“a.s.” denote convergences, respectively, in distri- bution, in probability−→ −→ and with probability−→ one.
License or copyright restrictions may apply to redistribution; see https://www.ams.org/journal-terms-of-use DIFFUSIONS WITH DRIFTED BROWNIAN POTENTIALS 3917
Throughout the paper, we write π (1.5) c def=2γ+log , 0 4 πx 1/x (1.6) ψ(x) def= , 0
License or copyright restrictions may apply to redistribution; see https://www.ams.org/journal-terms-of-use 3918 YUEYUN HU, ZHAN SHI, AND MARC YOR
Notation. The usual symbol a(x) b(x)(x x0) means limx x0 a(x)/b(x)= 1. We use BES(d) (resp. BESQ(d))∼ to denote→ a Bessel process→ (resp. squared
Bessel process) of dimension d. Furthermore, log2 stands for the iterated logarithm function. Unless stated otherwise, ε( ) denotes an asymptotically negligible random · term in the sense that ε(t) P. 0(ast ), whose value may change from one −→ →∞ line to another. Finally, we shall also use, whenever possible, the notation d0 for d/2, which will often be convenient in our formulae involving dimension parameters d1 and d2. Acknowledgment. We are grateful to the referee for valuable comments which helped us to improve the presentation of our paper and give more accurate refer- ences.
2. Bessel and Jacobi processes
Throughout this section, R1(t); t 0 and R2(t); t 0 denote independent Bessel processes, of dimensions{ d and≥d }respectively.{ We≥ assume} R (0) = r 0 1 2 1 1 ≥ and R2(0) = r2 > 0. Our aim is to study the asymptotic behaviour of t R2(s) (2.1) Ξ(t) def= 1 d s, R4(s) Z0 2 as t tends to infinity. Observe that, in order to ensure the finiteness of Ξ(t), we have to take d2 2, a hypothesis which, unless stated otherwise, we assume throughout this section.≥ Here is the main estimate of this section, which plays an important role in the def def proof of Theorem 1.1 in Section 3. Recall that d10 = d1/2andd20 = d2/2.
Proposition 2.1. Let d1 > 0 and d2 > 2, and let Ξ(t) be as in (2.1).Whent tends to infinity (writing h(t) def= d (log t)/(d 2)(d 4) for brevity), 1 2 − 2 − Ξ(t) law ca (2.2) c 2 , 2
1/(d0 1) def 1 8 2− c = ψ(d0 1) , 2 2 2 − (d 2)2 B(d ,d 1) 2 − 10 20 − 2 1 h i def d1 d0 1 d1(c0 log(d1 +2)) c = y 1− log(1 y)dy+ − , 3 8 − 4 Z0 1/2 def d1 1 1 2 d2 2 c4 = + + + − , (d20 1)(d2 4) d2 2 d2 4 d2 6 d1(d2 6) − − h − − − − i and B( , ) denotes the beta function. · ·
License or copyright restrictions may apply to redistribution; see https://www.ams.org/journal-terms-of-use DIFFUSIONS WITH DRIFTED BROWNIAN POTENTIALS 3919
Remark 2.2. We shall not use (2.6) in the paper. It is presented here only for the sake of completeness. It is also possible to treat the case d1 > 0, d2 = 2, for which one can prove that log Ξ(t) 1 law ,t , log t −→ →∞ E where denotes an exponential variable with mean 1. E The proof of Proposition 2.1 relies on some preliminaries, stated below as Lem- mas 2.3 and 2.8. Lemma 2.3. Let β be a standard Brownian motion, with local time Lx; t 0,x { t ≥ ∈ def 0 R .Letτ(r)=inf s>0: Ls >r (for r>0) denote its inverse local time at 0. Fix}a>0,b>0,µ{1and λ>0.Wedefinethefunction} ≥ y d x S(y)def= ,y(0, 1), xa(1 x)b+1 ∈ Z1/2 − and the process 1 def a+µ 1 b µ 1 S(y)/t J(t) = y − (1 y) − − L d y, t > 0. − τ(λ) Z0 Then we have λB(a+µ, b µ)+δ(t),b>µ, − 1 Lx λ Lx 1τ(λ)− ∞ τ(λ) J(t)=c5λ+bλlog t + 0 x d x + 1 x d x + ε(t),b=µ, µ/b 2 ∞ µ/b 2 x µ/b 1 b − 0 x − RLτ (λ) d x + ε(t) tR − ,b<µ, a.s. where δ(t) 0(as t R ),B( , )is the usual beta function, and −→ →∞ · · 1 def log b a+b 2 c = +(a+b 1) y − log(1 y)dy. 5 b − − Z0 Proof of Lemma 2.3. Looking at the integral expression in J(t), it is seen that the only possible explosion comes from the neighbourhood of 1. However, if b>µ,the 1 a+µ 1 b µ 1 integral 0 y − (1 y) − − d y being finite, Lemma 2.3 follows trivially from the dominated convergence− theorem. It remainsR to treat the case b µ.Write ≤ 1/2 1 a+µ 1 b µ 1 S(y)/t J(t)= + y − (1 y) − − Lτ(λ) d y (2.7) 0 1/2 − Z Z def = J1(t)+J2(t).
Again, we can apply the dominated convergence theorem to J1(t) to arrive at 1/2 a+µ 1 b µ 1 (2.8) J (t)=λ y − (1 y) − − d y + ε(t). 1 − Z0 1 We have to deal with J2(t). Writing S− for the inverse function of S,andbya change of variables x = S(y)/t,
∞ 1 2a+µ 1 1 2b µ x (2.9) J (t)=t (S− (tx)) − (1 S− (tx)) − L d x. 2 − τ (λ) Z0
License or copyright restrictions may apply to redistribution; see https://www.ams.org/journal-terms-of-use 3920 YUEYUN HU, ZHAN SHI, AND MARC YOR
Observe that 1 1/b (2.10) 1 S− (v) (bv)− ,v . − ∼ →∞ Therefore, in case b<µ, we have, by the dominated convergence theorem, 1 2b µ µ/b 1 µ/b 2 ∞ 1 2a+µ 1 (1 S− (tx)) − x J (t)=t − b − (S− (tx)) − − L d x 2 (bt)µ/b 2 τ (λ) Z0 − Lx = tµ/b 1bµ/b 2 ∞ τ (λ) d x + ε(t) , − − 2 µ/b 0 x − Z which, in view of (2.7)–(2.8), yields Lemma 2.3 (the constant term in J1(t)being absorbed by ε(t)). Finally, let us assume b = µ. By (2.9),
∞ 1 2a+b 1 1 b x J2(t)=t (S− (tx)) − (1 S− (tx)) Lτ (λ) d x 1 − Z 1 1 2a+b 1 1 b x + t (S− (tx)) − (1 S− (tx)) (Lτ (λ) λ)dx 0 − − Z 1 1 2a+b 1 1 b +λt (S− (tx)) − (1 S− (tx)) d x. − Z0 Now, on the right hand side, we apply the dominated convergence theorem (via observation (2.10), of course) to the first two terms, and use a change of variables x = S(y)/t for the third term. Accordingly, x 1 x 1 ∞ Lτ(λ) 1 Lτ(λ) λ J2(t)= dx+ − dx+ε(t) b 1 x b 0 x (2.11) Z 1 Z S− (t) ya+b 1 +λ − dy. 1 y Z1/2 − To treat the last integral term on the right hand side, let us recall (2.10) and use integration by parts, which leads to: 1 S− (t) ya+b 1 λ log t λ − d y = + o(1) 1 y b (2.12) Z1/2 − log b log 2 1 + λ +(a+b 1) ya+b 2 log(1 y)dy . a+b 1 − b − 2 − − 1/2 − Z Combining (2.7)–(2.8) and (2.11)–(2.12) completes the proof of Lemma 2.3. To get our second preliminary estimate, we have to recall the definition of the Jacobi processes (cf., for example, Karlin and Taylor [18]): d Y (t)=2 Y(t)(1 Y (t)) d B(t)+(d (d +d )Y(t)) d t, − 1 − 1 2 Y (0) = a (0, 1), p∈ where B is a standard Brownian motion. In the literature, Y (t); t 0 is some- { ≥ } times referred to as a Jacobi process of dimensions (d1,d2), starting from a. We need the following known results for Jacobi processes and Brownian local times. Fact 2.4 is due to Warren and Yor [33]. Fact 2.5 is a consequence of the ergodic theorem on path space, cf. for example Revuz and Yor [26, p. 411]. Fact 2.7 is due to McKean [22] and Ray [24]; see, e.g. Itˆo and McKean [17, p. 65 et seq.]; this estimate has been extended by Barlow and Marcus–Rosen to a large class of
License or copyright restrictions may apply to redistribution; see https://www.ams.org/journal-terms-of-use DIFFUSIONS WITH DRIFTED BROWNIAN POTENTIALS 3921
L´evy processes. Fact 2.6, which clearly shows where the stable laws in Theorem 1.1 come from, can be found in Biane and Yor [3], except for (2.18), whose proof is postponed until Section 4. We point out that the value of c0 in (2.16) is given in (1.5); this corrects a misprint made in Biane and Yor [3, p. 66], in passing from the formula on their (correct) 4th line to their (wrong) 8th line.
Fact 2.4 (Warren and Yor). Let R1(t); t 0 and R2(t); t 0 be independent Bessel processes of dimensions d and{ d respectively,≥ } with{ d +d≥ }2, R (0) = r 1 2 1 2 ≥ 1 1 ≥ 0andR2(0) = r2 > 0. There exists a Jacobi process Y (t); t 0 of dimensions 2 2 { 2 ≥2 } (d1,d2), starting from r1/ r1 + r2, independent of R1(t)+R2(t); t 0 ,such that for all t 0, { ≥ } ≥ p 2 t R1(t) ds (2.13) 2 2 =Y 2 2 . R1(t)+R2(t) 0 R1(s)+R2(s) Z Fact 2.5 (Revuz and Yor). If R(t); t 0 is a Bessel process of dimension d>2, starting from r>0, then { ≥ } 1 t d u 1 (2.14) a.s. ,t . log t R2(u) −→ d 2 →∞ Z0 − Moreover, 1 t d u log t 2 (2.15) law ,t . 1/2 2 3/2 (log t) 0 R (u) − d 2 −→ ( d 2) N →∞ Z − − Fact 2.6 (Biane and Yor). In the notation of Lemma 2.3, for any λ>0, 0
0 d β(s) = α − (α 1) ψ(α) λ α , { } − S Z0 where “law= ” denotes identity in law. Fact 2.7 (McKean, Ray). In the notation of Lemma 2.3, for any T>0, there exists a (finite) random variable Θ(T ), depending only on T , such that
x y 1 sup Ls Ls Θ(T ) x y log , 0 s T | − |≤ s | − | x y ≤ ≤ | − | for all (x, y) R2 with x y 1/2. ∈ | − |≤ Our second preliminary estimate concerns a class of processes which includes Ξ(t) (cf. (2.1)) as a special case.
Lemma 2.8. Let R1 and R2 be independent Bessel processes of dimensions d1 > 0 and d > 2 respectively, with R (0) = r 0 and R (0) = r > 0.Fixµ 1and 2 1 1 ≥ 2 2 ≥ t 2µ def R1 (s) Ξµ(t) = d s, t > 0. R2µ+2(s) Z0 2
License or copyright restrictions may apply to redistribution; see https://www.ams.org/journal-terms-of-use 3922 YUEYUN HU, ZHAN SHI, AND MARC YOR
We have, as t goes to infinity,
Ξµ(t) a.s. B ( d 10 + µ, d20 µ 1) (2.19) − − ,d0>µ+1, log t −→ (d + d 2)B(d ,d ) 2 1 2 − 10 20 Ξµ(t) c6 (log t)(log2 t) law ca (2.20) − c 7 + c 8 ,d20=µ+1, log t −→ C 4/B(d10 ,d20 )
Ξµ(t) law ca (2.21) c 9 , 1