arxiv: v2 [math.ds] 13 May 2012

Transcription

1 ESCAPE PROBABILITY FOR STOCHASTIC DYNAMICAL SYSTEMS WITH JUMPS* HUIJIE QIAO 1, XINGYE KAN AND JINQIAO DUAN arxiv: v [math.ds] 13 May 1 1. Department of Mathematics, Southeast University Nanjing, Jiangsu 11189, China hjqiaogean@yahoo.com.cn. Institute for Pure and Applied Mathematics, University of California Los Angeles, CA 995, USA & Department of Applied Mathematics, Illinois Institute of Technology Chicago, IL 6616, USA kan@hawk.iit.edu, duan@iit.edu Abstract. The escape probability is a deterministic concept that quantifies some aspects of stochastic dynamics. This issue has been investigated previously for dynamical systems driven by Gaussian Brownian motions. The present work considers escape probabilities for dynamical systems driven by non-gaussian Lévy motions, especially symmetric α-stable Lévy motions. The escape probabilities are characterized as solutions of the Balayage-Dirichlet problems of certain partial differential-integral equations. Differences between escape probabilities for dynamical systems driven by Gaussian and non-gaussian noises are highlighted. In certain special cases, analytic results for escape probabilities are given. Dedicated to Professor David Nualart on the occasion of his 6th birthday 1. Introduction Stochastic dynamical systems arise as mathematical models for comple phenomena in biological, geophysical, physical and chemical sciences, under random fluctuations. A specific orbit (or trajectory) for such a system could vary wildly from one realization to another, unlike the situation for deterministic dynamical systems. It is desirable to have different concepts for quantifying stochastic dynamical behaviors. The escape probability is such a concept. Brownian motions are Gaussian stochastic processes and thus are appropriate for modeling Gaussian random fluctuations. Almost all sample paths of Brownian motions are Date: February 7, 18(Revised version). AMS Subject Classification(1): 6H1, 6J75, 35S15; 31C5. Keywords: Escape probability, Balayage-Dirichlet problem, Lévy processes, discontinuous stochastic dynamical systems, nonlocal differential equation, non-gaussian noise. *This work is partially supported by the NSF of China (No and No. 1181) and the NSF grant DMS

2 continuous in time. For a dynamical system driven by Brownian motions, almost all orbits (or paths or trajectories) are thus continuous in time. The escape probability is the likelihood that an orbit, starting inside an open domain D, eits this domain first through a specific part Γ of the boundary D. This concept helps understand various phenomena in sciences. One eample is in molecular genetics [6]. The frequency of collisions of two single strands of long helical DNA molecules that leads to a double-stranded molecule is of interest and can be computed by virtu of solving an escape probability problem. It turns out that the escape probability satisfies an elliptic partial differential equation with properly chosen boundary conditions [6, 18, 5, 5]. Non-Gaussian random fluctuations are widely observed in various areas such as physics, biology, seismology, electrical engineering and finance [9, 16, ]. Lévy motions are a large class of non-gaussian stochastic processes whose sample paths are discontinuous in time. For a dynamical system driven by Lévy motions, almost all the orbits X t are discontinuous in time. In fact, these orbits are càdlàg (right-continuous with left limit at each time instant), i.e., each of these orbits has countable jumps in time. Due to these jumps, an orbit could escape an open domain without passing through its boundary. In this case, the escape probability is the likelihood that an orbit, starting inside an open domain D, eits this domain first by landing in a target domain U in D c (the complement of domain D). As we see, the escape probability is defined slightly differently for dynamical systems driven by Gaussian or non-gaussian processes. Although the escape probability for the former has been investigated etensively, the characterization for the escape probability for the latter has not been well documented as a dynamical systems analysis tool for applied mathematics and science communities. See our recent works [6, 11] for numerical analysis of escape probability and mean eit time for dynamical systems driven by symmetric α-stable Lévy motions. In this paper, we carefully derive a partial differential-integral equation to be satisfied by the escape probability for a class of dynamical systems driven by Lévy motions, especially symmetric α-stable Lévy motions. Namely the escape probability is a solution of a nonlocal differential equation. We highlight the differences between escape probabilities for dynamical systems driven by Gaussian and non-gaussian processes. These are illustrated in a few eamples. More precisely, let {X t,t } be a R d -valued Markov process defined on a complete filtered probability space (Ω,F,{F t } t,p). Let D be an open domain in R d. Define the eit time τ D c := inf{t > : X t D c }, where D c is the complement of D in R d. Namely, τ D c is the first time when X t hits D c. When X t has almost surely continuous paths, i.e., X t is either a Brownian motion or a solution process for a dynamical system driven by Brownian motions, a path starting at D will hit D c by hitting D first (assume for the moment that D is smooth). Thus τ D c = τ D. Let Γ be a subset of the boundary D. The likelihood that X t, starting at

3 , eits from D first through Γ is called the escape probability from D to Γ, denoted as p(). That is, p() = P{X τ D Γ}. We will verify that (Section 3.) the escape probability p() solves the following Dirichlet boundary value problem: { Lp =, D, (1) p D = ψ, wherelistheinfinitesimal generatoroftheprocessx t andtheboundarydataψ isdefined as follows { 1, Γ, ψ() =, D \Γ. When X t has càdlàg paths which have countable jumps in time, i.e., X t could be either a Lévy motion or a solution process of a dynamical system driven by Lévy motions, the first hitting of D c may occur somewhere in D c. For this reason, we take a subset U of D c, and define the likelihood that X t eits firstly from D by landing in the target set U as the escape probability from D to U, also denoted by p(). That is, p() = P{X τd c U}. We will demonstrate that (Section 3.4) the escape probability p() solves the following Balayage-Dirichlet boundary value problem: { Ap =, D, () p D c = ϕ, where A is the characteristic operator of X t and ϕ is defined as follows { 1, U, ϕ() =, D c \U. Therefore by solving a deterministic boundary value problem (1) or (), we obtain the escape probability p(). This paper is arranged as follows. In Section, we introduce Balayage-Dirichlet problem for discontinuous Markov processes, and also define Lévy motions. The main result is stated and proved in Section 3. In Section 4, we present analytic solutions for escape probabilities in a few special cases.. Preliminaries In this section, we recall basic concepts and results that will be needed throughout the paper. 3

4 .1. Balayage-Dirichlet problem for discontinuous Markov processes. The following materials are from [3, 17, 7, 7, 19, 1]. Let G be a locally compact space with a countable base and G be the Borel σ-field of G. Also, ς is adjoined to G as the point at infinity if G is noncompact, and as an isolated point if G is compact. Furthermore, let G ς be the σ-field of Borel sets of G ς = G {ς}. Definition.1. A Markov process Y with state space (G, G) is called a Hunt process provided: (i) The paths functions t Y t are right continuous on [, ) and have left-hand limits on [,ζ) almost surely, where ζ := inf{t : Y t = ς}. (ii) Y is strong Markov. (iii) Y is quasi-left-continuous: whenever {τ n } is an increasing sequence of F t -stopping times with limit τ, then almost surely Y τn Y τ on {τ < }. Definition.. Let G be an open subset of G and Y t () be a Hunt process starting at G. A nonnegative function h defined on G is said to be harmonic with respect to Y t in G if for every compact set K G, E[h(Y τk c())] = h(), G. Definition.3. Let f be nonnegative on G c. We say h defined on G solves the Balayage- Dirichlet problem for G with boundary value f, denoted by (G,f), if h = f on G c, h is harmonic with respect to Y t in G and further satisfies the following boundary condition: z G, h(y) f(z), as y z from inside G. A point z G is called regular for G c with respect to Y t (z) if P{τ G c = } = 1. Here G is said to be regular if any z G is regular for G c. Let ρ be a metric on G compatible with the given topology. Let I G be the family of functions g bounded on G and lower semicontinuous in G such that G, there is a number Ag() satisfying E[g(Y τε ())] g() E[τ ε ] Ag(), as ε, where τ ε := inf{t > : ρ(y t (),) > ε}. We call A with domain I G the characteristic operator of Y t relative to G. If L with domain D G is the infinitesimal generator of Y t relative to G, D G I G and Af = Lf, f D G. (c.f. [1]) We quote the following result about the eistence and regularity of the solution for the Balayage-Dirichlet problem. Theorem.4. ([17]) Suppose that G is relatively compact and regular, and f is nonnegative and bounded on 4

5 G c. If f is continuous at any z G, then h() = E[f(Y τg c())] is the unique solution to the Balayage-Dirichlet problem (G,f), and Ah() = for h I G... Lévy motions. Definition.5. A process L t, with L = a.s. is a d-dimensional Lévy process or Lévy motion if (i) L t has independent increments; that is, L t L s is independent of L v L u if (u,v) (s,t) = ; (ii) L t has stationary increments; that is, L t L s has the same distribution as L v L u if t s = v u > ; (iii) L t is stochastically continuous; (iv) L t is right continuous with left limit. The characteristic function for L t is given by E(ep{i z,l t }) = ep{tψ(z)}, z R d, where, is the scalar product in R d. The function Ψ : R d C is called the characteristic eponent of the Lévy process L t. By the Lévy-Khintchine formula, there eist a nonnegative-definite d d matri Q, a measure ν on R d satisfying ν({}) = and ( u 1)ν(du) <, and γ R d such that R d \{} Ψ(z) = i z,γ 1 z,qz ( ) + e i z,u 1 i z,u 1 u 1 ν(du). (3) R d \{} The measure ν is called the Lévy measure of L t, Q is the diffusion matri, and γ is the drift vector. We now introduce a special class of Lévy motions, i.e., the symmetric α-stable Lévy motions L α t. Definition.6. For α (,). A d-dimensional symmetric α-stable Lévy motion L α t is a Lévy process with characteristic eponent where Ψ(z) = C z α, z R d, (4) C = π 1/Γ((1+α)/)Γ(d/). Γ((d+α)/) (c.f. [4, Page 115] for the above formula of C.) 5

6 Thus, for a d-dimensional symmetric α-stable Lévy motion L α t, the diffusion matri Q =, the drift vector γ =, and the Lévy measure ν is given by where ν(du) = C d,α u d+αdu, C d,α = αγ((d+α)/) 1 α π d/ Γ(1 α/). (c.f. [8, Page 131] for the above formula of C d,α.) Moreover, comparing (4) with (3), we obtain ( ) C z α = e i z,u C d,α 1 i z,u 1 u 1 u d+αdu. R d \{} Let C (R d ) be the space of continuous functions f on R d satisfying lim f() = with norm f C (R d ) = sup f(). Let C (Rd ) be the set of f C (R d ) such that f is two times R d differentiable and the first and second order partial derivatives of f belong to C (R d ). Let Cc (R d ) stand for the space of all infinitely differentiable functions on R d with compact supports. Define ( ) C d,α (L α f)() := f(+u) f() f(),u 1 u 1 u d+αdu R d \{} on C (R d ). And then for ξ R d (L α e i,ξ )() = e i,ξ R d \{} ( e i u,ξ 1 i ξ,u 1 u 1 ) C d,α u d+αdu. By Courrège s second theorem ([1, Theorem 3.5.5, p.183]), for every f Cc (R d ) 1 (L α f)() = e i z, [ e i,z (L (π) d/ α e i,z )() ] ˆf(z)dz R d [ ] 1 ( ) = e i z, e i u,z C d,α 1 i z,u 1 (π) d/ u 1 R d R d \{} u d+αdu ˆf(z)dz = C e i z, z α ˆf(z)dz (π) d/ R d = C [ ( ) α/ f](). Set p t := L t L t. Then p t defines a stationary (F t )-adapted Poisson point process with values in R d \{} ([13]). And the characteristic measure of p is the Lévy measure ν. Let N p ((,t],du) be the counting measure of p t, i.e., for B B(R d \{}) N p ((,t],b) := #{ < s t : p s B}, where # denotes the cardinality of a set. The compensator measure of N p is given by Ñ p ((,t],du) := N p ((,t],du) tν(du). 6

7 The Lévy-Itô theorem states that for a symmetric α-stable process L t, (i) for 1 α <, t t L t = uñp(ds,du)+ un p (ds,du), (ii) for < α < 1, u 1 L t = t R d \{} u >1 un p (ds,du). 3. Boundary value problems for escape probability In this section, we formulate boundary value problems for the escape probability associated with Brownian motions, SDEs driven by Brownian motions, Lévy motions and SDEs driven by Lévy motions. For Lévy motions, in particular, we consider symmetric α-stable Lévy motions. We will see that the escape probability can be found by solving deterministic partial differential equations or partial differential-integral equations, with properly chosen boundary conditions. a D δ +δ b Figure 1. A particle eecuting unbiased random walk in a bounded interval 3.1. Boundary value problem for escape probability of Brownian motions. Suppose that a particle eecutes an unbiased random walk on a straight line. Let D = (a,b). Figure 1 shows the random walk scenario. That is, a particle moves according to the following rules ([18]): (i) During the passage of a certain fied time interval, a particle takes 1 step of a certain fied length δ along the ais. (ii) It is equally probable that the step is to the right or to the left. If the particle starting from D eventually escapes D by crossing the boundary b, then it must have moved to one of the two points adjacent to first and then crossed the boundary. Thus p() = 1 [p( δ)+p(+δ)], for D. By Taylor epansion on the right hand side to the second order, we have 1 p () =. The boundary conditions are lim b p() = 1, limp() =, a 7

8 since the nearer the particle starts to b, the more likely it will first cross the boundary through b. Note that the limit of the random walk is a standard Brownian motion W t, that is, (i) W has independent increments; (ii) for < s < t, W t W s is a Gaussian random variable with mean zero and variance (t s). Thus, the escape probability p() of a standard Brownian motion from D through the boundary b satisfies 1 p() =, p(b) = 1, p(a) =, where 1 = 1 is the infinitesimal generator for a scalar standard Brownian motion W t. 3.. Boundary value problem for escape probability of SDEs driven by Brownian motions. Some results in this subsection can be found in [1, Chapter 9]. Let {W(t)} t be an m-dimensional standard F t -adapted Brownian motion. Consider the following stochastic differential equation (SDE) in R d : X t () = + t b(x s ())ds+ t σ(x s ())dw s. (5) We make the following assumptions about the drift b : R d R d and the diffusion coefficient σ : R d R d R m. (H 1 b,σ ) b() b(y) λ( y ), σ() σ(y) γ( y ). Here λ and γ are increasing concave functions with the properties λ() = γ() =, 1 du = 1 du =. + λ(u) + γ (u) Under (H 1 b,σ ), it is well known that there eists a unique strong solution to Eq.(5)([3]). This solution is denoted by X t (). We also make the following assumption. (H σ ) There eists a ξ > such that for any,y D y,σσ ()y ξ y. This condition guarantees that the infinitesimal generator d L := b i () d + a ij () i i j i=1 i,j=1 for Eq.(5) is uniformly elliptic in D, since then the eigenvalues of σσ are away from in D. Here the matri [a ij ] := 1 σ()σ (). 8

9 D Γ Figure. Escape probability for SDEs driven by Brownian motions: An annular open domain D with a subset Γ of its boundary D Let D be an open annular domain as in Figure. In one dimensional case, it is just an open interval. Let Γ be its inner (or outer) boundary. Taking ψ() = { 1, Γ,, D \Γ, (6) we have E[ψ(X τ D ())] = ψ(x τ D ())dp(ω) {ω:x τ D () Γ} + {ω:x τ D () D\Γ} = P{ω : X τ D () Γ} = p(). ψ(x τ D ())dp(ω) This means that, for this specific ψ, E[ψ(X τ D ())] is the escape probability p(), which we are looking for. We need to use [1, Theorem 9..14] or [9] in order to see that the escape probability p() is closely related to a harmonic function with respect to X t. This requires that the boundary data ψ to be bounded and continuous on D. For the domain D taken as in Figure, with Γ the inner boundary (or outer) boundary, the above chosen ψ in (6) is indeed bounded and continuous on D. Thus, we have the following result by [1, Theorem 9..14]. Theorem 3.1. The escape probability p() from an open annular domain D to its inner (or outer) boundary Γ, for the dynamical system driven by Brownian motions (5), is the 9

10 solution to the following Dirichlet boundary value problem Lp =, p Γ = 1, p D\Γ =. a D +y b Figure 3. A particle eecuting Lévy motion in a bounded interval 3.3. Boundary value problem for escape probability of symmetric α-stable Lévy motions. Assume a particle is taking a one-dimensional Lévy flight, where the distribution of step sizes is a symmetric α-stable distribution (Figure 3). Let p() denote the escape probability of the particle starting at in D = (a,b) and then first escapes D over the right boundary b. It could first move to somewhere inside D, say +y D, and then achieve its goal by jumping over the right boundary b from the new starting point +y. More precisely, p() = P{the first step length is y}p(+y)dy. (7) R\{} According to [], the symmetric α-stable probability density function is the following { 1 ( 1) k Γ(αk+1) π k=1 sin[k( απ αargy)], < α < 1, k! y y f α, (y) = αk 1 π k= ( 1)kΓ(k+1 α ) k!α yk cos[k( π )], 1 < α <, where argy = π when y <. For < α <, the asymptotic epansion has also been given by [] as follows f α, (y) = = { 1 n ( 1) k π k=1 k! 1 π Γ(αk+1) y y αk sin [ k( απ αargy)] +o( y α(n+1) 1 ), y, n k= ( 1)kΓ(k+1 α ) k!α yk cos [ k( π)] +o( y n+1 ), y, { C 1 (α)/ y 1+α +o( y α 1 ), y, C (α)+o( y ), y, where C 1 (α) = 1 π sin(πα)γ(1+α) and C (α) = 1 Γ(1/α). Take N > large enough and fi π α it. Thus, = f α, (y)[p(+y) p()]dy R\{} = f α, (y)[p(+y) p()]dy ( N,N)\{} 1

11 + f α, (y)[p(+y) p()]dy R\( N,N) =: I 1 +I. For I 1, by self-affine property in [8], we obtain I 1 = f α, (N y ( N,N)\{} N )[p(+y) p()]dy f α, (N) = ( N,N)\{} ( y /N) 1+α[p(+y) p()]dy C 1 (α) = ( N,N)\{} N 1+α ( y /N) 1+α[p(+y) p()]dy C 1 (α) = y 1+α[p(+y) p()]dy. (8) For I, we calculate R\{} I = = = = N N N ( N,N)\{} f α, (y)[1 p()]dy + [ ( C1 (α) 1 y +o 1+α N C 1 (α) dy y1+α R\[ N,N] y 1+α N )] [1 p()]dy [ C1 (α) +o ( y) 1+α R\[ N,N] C 1 (α)p(+y) dy y 1+α f α, (y)[ p()]dy ( )] 1 p()dy ( y) 1+α C 1 (α) y 1+αp()dy R\[ N,N] C 1 (α) y 1+αp()dy (9) Notethat for < α < 1, by the fact that the integral of an oddfunction on a symmetric interval is zero, it holds that p C 1 (α) ()yi { y 1} y 1+αdy = y p ()C 1 (α) y 1+αdy =. (1) { y 1}\{} Thus, putting (8), (9) and (1) together, we have for < α < 1 [ ] p(+y) p() p C 1 (α) ()yi { y 1} y 1+αdy =. R\{} Moreover, C 1 (α) = C 1,α. For α [1,), we only divide I 1 into two parts I 11 and I 1, where I 11 := f α, (y)[p(+y) p()]dy, { y }\{} 11

12 I 1 := ( N,N)\(, ) f α, (y)[p(+y) p()]dy, and > is a small enough constant. For I 11, by Taylor epansion and self-affine property in [6], we get f α, (y)[p(+y) p()]dy { y }\{} = f α, ( 1 y)p ()ydy = = = { y }\{} { y }\{} { y }\{} R\{} f α, ( 1 ) ( y ) 1+αp ()ydy C 1 (α) ( 1 )1+α ( y ) 1+αp ()ydy p ()yi { y } C 1 (α) y 1+αdy. For I 1, we apply the same technique as that in dealing with I 1 for α (,1). Net, by the similar calculation to that for α (,1), we obtain for α [1,) [ ] p(+y) p() p C 1 (α)dy ()yi { y 1} =. y 1+α R\{} Since the limit of the Lévy flight is a symmetric α-stable Lévy motion L α t, the escape probability p() of a symmetric α-stable Lévy motion, from D to [b, ) satisfies ( ) α p() =, p() [b, ) = 1, p() (,a] =. Note that ( ) α is the infinitesimal generator for a scalar symmetric α-stable Lévy motion L α t Boundary value problem for escape probability of SDEs driven by general Lévy motions. Let L t bealévy process independent ofw t. Consider thefollowing SDE in R d X t () = + t b(x s ())ds+ t σ(x s ())dw s +L t. (11) Assume that the drift b and the diffusion σ satisfy the following conditions: (H b ) there eists a constant C b > such that for,y R d b() b(y) C b y log( y 1 +e); (H σ ) there eists a constant C σ > such that for,y R d σ() σ(y) C σ y log( y 1 +e). 1

13 Under (H b ) and (H σ ), it is well known that there eists a unique strong solution to Eq.(11)(see [3]). This solution will be denoted by X t (). Moreover, X t () is continuous in. Lemma 3.. The solution process X t () of the SDE (11) is a strong Markov process. Proof. Let η be a (F t ) t -stopping time. Set G t := σ{w η+t W η,l η+t L η } N, t, where N is of all P-zero sets. That is, G t is a completed σ-algebragenerated by W η+t W η and L η+t L η. Besides, G t is independent of F t. Let X(,η,η + t) denote the unique solution of the the following SDE X(,η,η+t) = + η+t η b(x(,η,s))ds+ η+t η σ(x(,η,s))dw s +L η+t L η. (1) Moreover, X(,η,η + t) is G t -measurable and X(,,t) = X t (). By the uniqueness of the solution to (1), we have X(,,η+t) = X(X(,,η),η,η+t), a.s.. For any bounded measurable function g, E[g(X η+t ()) F η ] = E[g(X(,,η+t)) F η ] = E[g(X(X(,,η),η,η+t)) F η ] = E[g(X(y,η,η+t))] y=x(,,η) = E[g(X(y,,t))] y=x(,,η). (13) Here the last equality holds because the distribution of X(y,η,η+t) is the same to that of X(y,,t). The proof is completed since (13) implies that E[g(X η+t ()) F η ] = E[g(X η+t ()) X η ()]. Because L t has càdlàg and quasi-left-continuous paths([4]), X t () also has càdlàg and quasi-left-continuous paths. Thus by Lemma 3. and Definition.1 above, we see that X t () is a Hunt process. Let D be a relatively compact and regular open domain (Figure 4 or Figure 5). Theorem.4 implies that E[ϕ(X τd c())] is the unique solution to the Balayage-Dirichlet problem (D, ϕ), under the condition that ϕ is nonnegative and bounded on D c. Set { 1, U, ϕ() =, D c \U. Then ϕ is nonnegative and bounded on D c. We observe that E[ϕ(X τd c())] = ϕ(x τd c())dp(ω) {ω:x τd c() U} 13

14 D U Figure 4. Escape probability for SDEs driven by Lévy motions: an open annular domain D, withits inner part U (which isind c ) asatargetdomain D U Figure 5. Escape probability for SDEs driven by Lévy motions: a general open domain D, with a target domain U in D c + ϕ(x τd c())dp(ω) {ω:x τd c() D c \U} = P{ω : X τd c() U} 14

15 = p(). This means that, for this specific ϕ, E[ϕ(X τd c())] is the escape probability p() that we are looking for. By the definition of the characteristic operator, p I D and by Theorem.4, Ap() =. Thus we obtain the following theorem. Theorem 3.3. Let D be a relatively compact and regular open domain, and let U be a set in D c. Then the escape probability p(), for the dynamical system driven by Lévy motions (11), from D to U, is the solution of the following Balayage-Dirichlet problem Ap =, p U = 1, p D c \U =, where A is the characteristic operator for this system. Remark 3.4. Unlike the SDEs driven by Brownian motions, a typical open domain D here could be a quite general open domain (Figure 5), as well as an annular domain (Figure 4). This is due to the jumping properties of the solution paths. It is also due to the fact that, in Theorem.4, the function f is only required to be continuous on the boundary D (not on the domain D c ). Finally we consider the representation of the characteristic operator A, for a SDE driven by a symmetric α-stable Lévy process L α t, with α (,): X t () = + t b(x s ())ds+ t σ(x s ())dw s +L α t. (14) Let us first consider the case of 1 α <. For f C (Rd ), applying the Itô formula to f(x τε ()), we obtain f(x τε ()) f() = τε + + τε τε y f(x s ),b(x s ) ds+ y f(x s ),σ(x s )dw s (f(x s +u) f(x s ))Ñp(ds,du) τε u 1 u >1 τε τε (f(x s +u) f(x s ))N p (ds,du) ( ) f(x s ) σ ik (X s )σ kj (X s )ds y i y j ( f(xs +u) f(x s ) u 1 y f(x s ),u ) C d,α u d+αduds. 15

16 Here and hereafter, we use the convention that repeated indices imply summation from 1 to d. Taking epectation on both sides, we get So, E[f(X τε ())] f() τε = E y f(x s ),b(x s ) ds+ 1 τε ( ) E f(x s ) σ ik (X s )σ kj (X s )ds y i y j τε ( +E f(xs +u) f(x s ) y f(x s ),u ) C d,α u d+αduds. R d \{} The infinitesimal generator L of Eq.(11) is as follows [1]: (Lf)() = f(),b() + 1 ( ) f() σ ik ()σ kj () i j + (f(+u) f() f(),u ) C d,α u d+αdu. E[f(X τε ())] f() E[τ ε ] R d \{} (Lf)() = Because (Lf)() is continuous in, E τ ε (Lf)(X s)ds E τ ε (Lf)()ds E[τ ε ] E[τ ε ] E τ ε (Lf)(X s) (Lf)() ds E[τ ε ] sup (Lf)(y) (Lf)(). y <ε E[f(X τε ())] f() Af() = lim = (Lf)(). ε E[τ ε ] Similarly, we also have A = L for < α < 1. Remark 3.5. The above deduction tells us Af = Lf for f C (R d ). If the considered driving process is not a symmetric α-stable Lévy motion, the domain of L is unclear and thus A = L may not be true. The corresponding escape probability p() is the solution of the following Balayage-Dirichlet problem (in terms of operator L, instead of A): Lp =, p U = 1, p D c \U =. 4. Eamples In this section we consider a few eamples. 16

17 p() Figure 6. Escape probability for 1-dimensional Brownian motion in Eample 4.1, r = Eample 4.1. In 1-dimensional case, take D = ( r,r) and Γ = {r}. For each D, the escape probability p() of X t = +W t from D to Γ satisfies the following differential equation 1 p () =, ( r,r), p(r) = 1, p( r) =. We obtain that p() = +r r for [ r,r]. It is a straight line (See Figure 6). In -dimensional case, take D = { R ;r < < R} and Γ = { R ; = r}. For every D, the escape probability p() of X t = + W t from D to Γ satisfies the following elliptic partial differential equation 1 p() =, D, p() =r = 1, p() =R =. By solving this equation, we obtain that p() = logr log. It is plotted in Figure 7. logr logr Eample 4.. Consider the following SDE driven by Brownian motions: dx t = b(x t )dt+σ(x t )dw t, where b and (non-zero) σ are real functions. When b and σ satisfy (H 1 b,σ ), the equation has a unique solution which is denoted as X t. We take D = ( r,r) and Γ = {r}. For 17

18 1.8 p( 1, ) Figure 7. Escape probability for -dimensional Brownian motion in Eample 4.1: r =,R = 4 each D, under the condition (H σ ), the escape probability p() satisfies The solution is for [ r,r]. See Figure 8. 1 σ ()p ()+b()p () =, ( r,r), p(r) = 1, p( r) =. p() = y b(z) r e r σ dz (z) dy y b(z) r r e r σ dz (z) dy Eample 4.3. In 1-dimensional case, take D = ( r,r) and U = [r, ). For each D and a symmetric α-stable Lévy process L α t, the escape probability p() of X t = + L α t from D to U satisfies the following differential-integral equation ( ) α p() =, ( r,r), p() [r, ) = 1, p() (, r] =. It is difficult to deal with this equation because of the fractional Laplacian operator. But we can solve it via Poisson kernel. From [14], for ( r,r), p() = sin πα π r (r ) α/ (y r ) α/ 1 (y ) dy. 18

19 p() Figure 8. Escape probability in Eample 4.: b() =,σ() = 1,r = Obviously, p( r) =. To justify p(r) = 1, we apply the substitution y = (r v)( v) 1 to obtain p(r) = sin πα π = sin πα π r r 1 (r v) α 1 (r v ) α dv (1 v) α 1 v 1 α 1 dv = sin πα π B(α,1 α ) = sin πα π Γ(α )Γ(1 α ) = 1, where the Beta and Gamma functions and their properties are used in the last two steps. The escape probability p() is plotted in Figure 9 for various α values. Acknowledgement. We have benefited from our previous collaboration with Ting Gao, Xiaofan Li and Renming Song [11]. We thank Ming Liao, Renming Song and Zhen-Qing Chen for helpful discussions. This work was done while Huijie Qiao was visiting the Institute for Pure and Applied Mathematics (IPAM), Los Angeles. 19

20 p p Α Α p p Α Α Figure 9. Escape probability in Eample 4.3: r = References [1] D. Applebaum: Lévy Processes and Stochastic Calculus, Second Edition, Cambridge Univ. Press, Cambridge, 9. [] H. Bergström: On some epansions of stable distribution functions, Ark. Math., (195) [3] R. M. Blumenthal and R. K. Getoor: Markov Processes and Potential Theory, Academic Press, New York, [4] R. M. Blumenthal, R. K. Getoor and D. B. Ray: On the distribution of first hits for the symmetric stable processes, Transactions of the American Mathematical Society, 99(1961) [5] J. Brannan, J. Duan and V. Ervin: Escape Probability, Mean Residence Time and Geophysical Fluid Particle Dynamics, Physica D 133 (1999), [6] H. Chen, J. Duan, X. Li and C. Zhang: A computational analysis for mean eit time under non- Gaussian Levy noises. Appl. Math. Comput., volume 18, issue 5, 11, pp [7] Z. Q. Chen: On notions of harmonicity, Proceedings of the American Mathematical Society, 137(9) [8] Z. Chen, P. Kim and R. Song: Heat kernel estimates for the Dirichlet fractional Laplacian, J. Eur. Math. Soc., 1(1) [9] Z. Q. Chen and Z. Zhao: Diffusion processes and second order elliptic operators with singular coefficients for lower order terms. Math. Ann. 3, (1995). [1] E. B. Dynkin: Markov Processes, Vol. I. Springer-Verlag, [11] T. Gao, J. Duan, X. Li and R. Song: Mean eit time and escape probability for dynamical systems driven by Lévy noise. Preprint, 11. [1] Q. Guan and Z. Ma: Boundary problems for fractional Laplacians, Stochastics and Dynamics, 5(5) [13] N. Ikeda and S. Watanabe: Stochastic differential equations and diffusion processes, nd ed., North- Holland/Kodanska, Amsterdam/Tokyo, [14] N. Jacob: Pseudo Differential Operators and Markov Processes, Volume III, Markov Processes and Applications, Imperial College Press, 5. [15] F. Klebaner: Introduction to Stochastic Calculus with Application. Imperial College Press, Second Edition, 5. [16] T. Koren, A. V. Chechkin and J. Klafter: On the first passage time and leapover properties of Lévy motions, Physica A 379, (7)1-.

21 [17] M. Liao: The Dirichlet problem of a discontinuous Markov process, Acta Mathematica Sinica, 5(1989)9-15. [18] C. C. Lin and L. A. Segel: Mathematics Applied to Deterministic Problems in the Natural Sciences, SIAM, [19] Z. Ma, R. Zhu and X. Zhu: On notions of harmonicity for non-symmetric Dirichlet form, Science in China (Mathematics) 53(1), [] B. Matkowsky, Z. Schuss, and C. Tier, Asymptotic Methods for Markov Jump Processes, A. M. S. Lectures in Applied Mathematics 7, 1991, p [1] B. Oksendal: Stochastic Differential Equations, An Introduction with Applications, Sith ed., Springer-Verlag, Berlin, 3. [] P. Protter: Stochastic Integration and Differential Equations, nd ed., Springer-Verlag, Berlin, 4. [3] H. Qiao and X. Zhang: Homeomorphism flows for non-lipschitz stochastic differential equations with jumps, Stochastic Processes and their Applications, 118(8) [4] K. Sato: Lévy processes and infinitely divisible distributions, Cambridge University Press, [5] Z. Schuss, Theory and Applications of Stochastic Differential Equations, Wiley & Sons, New York, 198. [6] W. Smith and G. S. Watson: Diffusion out of a triangle, J. Appl. Prob., 4(1967) [7] R. Song, Probabilistic approach to the Dirichlet problem of perturbed stable processes. Probab. Theory Related Fields, 95(3) (1993), [8] G. M. Viswanathan, V. Afanasyev, S. V. Buldyrev, S. Havlin, M. Daluz, E. Raposo and H. Stanley, Lévy flights in random searches, Physica A 8 (), 1-1. [9] W. A. Woyczynski, Lévy processes in the physical sciences. In Lévy Processes: Theory and Applications, O. E. Barndorff-Nielsen, T. Mikosch and S. I. Resnick (Eds.), 41-66, Birkhäuser, Boston, 1. [3] T. Yamada and S. Watanabe: On the uniqueness of solutions of stochastic differential equations, J. Math. Kyoto Univ., 11(1971)