Stochastic solutions for fractional wave equations
|
|
- Chloe Rice
- 5 years ago
- Views:
Transcription
1 Nonlinear Dyn DOI.7/s z ORIGINAL PAPER Stochastic solutions for fractional wave equations Mark M. Meerschaert René L. Schilling Alla Sikorskii Received: 6 November 3 / Accepted: 5 February 4 Springer Science+Business Media Dordrecht 4 Abstract A fractional wave equation replaces the second time derivative by a Caputo derivative of order between one and two. In this paper, we show that the fractional wave equation governs a stochastic model for wave propagation, with deterministic time replaced by the inverse of a stable subordinator whose index is one-half the order of the fractional time derivative. Keywords Fractional calculus Wave equation Inverse subordinator Stochastic solution Introduction The traditional wave equation p(x, t = x p(x, t (. M. M. Meerschaert (B A. Sikorskii Department of Statistics & Probability, Michigan State University, East Lansing, MI 4884, USA mcubed@stt.msu.edu URL: A. Sikorskii sikorska@stt.msu.edu URL: R. L. Schilling Institut für Mathematische Stochastik, Technische Universität, Dresden, Germany rene.schilling@tu-dresden.de URL: models wave propagation in an ideal conducting medium. Assume a plane wave solution p(x, t = e iωt+ikx for some frequency ω > and substitute into (. to see that we must have ( iω = (ik or in other words k =±ω. This solution is a traveling wave at speed one (justified by a suitable choice of units, and the general solution to (. can be written as a linear combination of plane waves. In a complex inhomogeneous conducting medium, experimental evidence reveals that sound waves exhibit power-law attenuation, with an amplitude that falls off at an exponential rate α = α(ω ω p for some powerlaw index p (e.g., see Duck [] for applications to medical ultrasound. A variety of modified wave equations have been proposed to model wave conduction in complex media [,4,,3]. We note here, apparently for the first time, that a simple time-fractional wave equation γ p(x, t = x p(x, t, <γ <, (. which replaces the second time derivative by a fractional derivative, also exhibits power-law attenuation. Assuming the same plane wave solution, and using the well-known formula dt dγ γ [e at ]=t γ e at (e.g., see [8, Example.6], we now have ( iω γ = k, and a little algebra yields k = β(ω + iα(ω with attenuation coefficient α(ω = α ω γ/. Hence, solutions to the time-fractional wave equation (. also exhibit powerlaw attenuation with power-law index p = γ/. The goal of this paper is to develop a new stochastic solution to the time-fractional wave equation (.. Our
2 M. M. Meerschaert et al. stochastic solution is based on limit theory for random walks and, therefore, provides a simple and illuminating statistical physics model for wave conduction in complex media. Background In one spatial dimension, the general solution to (. with initial conditions p(x, = φ(x and p(x, = ψ(x is given by the d Alembert formula p(x, t = x+t [φ(x + t + φ(x t] + ψ(ydy. x t (. In fact, for all continuous and exponentially bounded functions φ,ψ the unique solution to the equivalent integral equation p(x, t = φ(x + tψ(x + (t s x p(x, sds (. is given by the d Alembert formula. In this paper, we consider the following integral form of the fractional wave equation: t γ/ p(x, t = φ(x + Ɣ( + γ/ ψ(x + Ɣ(γ (t s γ x p(x, sds. (.3 The differential form of equation (.3 employs the Riemann Liouville fractional derivative. The Riemann Liouville fractional integral of non-integer order γ> is defined by I γ t f (t = Ɣ(γ (t s γ f (sds. (.4 The Riemann Liouville fractional derivative of noninteger order γ> is defined by d n D γ t f (t = Ɣ(n γ dt n (t s n γ f (sds where n is the smallest integer greater than γ. (.5 Equation (.3 corresponds to the following fractional differential equation: t γ D γ t p(x, t φ(x Ɣ( γ t γ/ ψ(x Ɣ( γ/ = x p(x, t (.6 with initial conditions p(x, = φ(x and / p(x, = ψ(x. γ/ In this paper, we will show that the solution to (.3 is p(x, t = E [φ(x + E t + φ(x E t ] + x+e E ψ(ydy (.7 x E t where E t is the inverse (hitting time or first passage time of a standard stable subordinator with index γ/. Then, using the general theory of second-order Cauchy problems, we will extend this result to a wide variety of fractional partial differential equations that model wave-like motions. Finally, we will develop random walk models that provide a physical explanation for these fractional wave equations. 3 Fractional wave equations Let D u be a standard stable subordinator with D = a.s. and Laplace transform E[e sd u ]=e usβ for some <β<. The random variable D has a smooth density function g β (u. Define the inverse subordinator (generalized inverse, first passage time, or hitting time E t = inf{u : D u > t} (3. for t. Then, a simple computation [6, Corollary 3.] shows that E t has a smooth density u h(u, t = t β u /β g β (tu /β, u >, t >. (3. Write R + =[,, letb(r R + denote the set of real-valued continuous functions p(x, t on R R + such that p(x, t Ae B( x +t for some constants A, B >, and denote by B(R the set of real-valued continuous functions φ(x on R such that φ(x Ae B x for some constants A, B >. Denote by B m (R
3 Stochastic solutions for fractional wave equations the set of real-valued continuously differentiable functions φ(x such that j φ B(R for all integers d dx j j m and by B m, (R R + the set of realvalued functions p(x, t on R R + continuously differentiable in x with d j p(x, t B(R R dx j + for all integers j m. Theorem 3. Let (x = x ψ(ydy. For any φ, ψ such that φ, B (R, the unique solution to the fractional wave equation (.3 in B, (R R + with p(x, = φ(x and β p(x, = ψ(x is given by β the formula (.7, where E t is the inverse stable subordinator (3. with index β = γ/. Proof The proof uses a result of Fujita [3] together with a duality result from [3]. Fujita considers a stable Lévy process X γ (t with characteristic function [ e ikx γ (t ] = exp [ t k /γ e i(π/( /γ sgn(k]. E (3.3 with index < /γ < along with its supremum process Y γ (t = sup u t X γ (u. (3.4 Fujita shows that for φ, ψ such that φ, B (R, the unique solution to (.3 inb, (R R + is p(x, t = E [ φ(x + Y γ (t + φ(x Y γ (t ] x+y + γ (t E ψ(y dy. (3.5 x Y γ (t Using the parameterization of Samorodnitsky and Taqqu [], the characteristic function of a generic stable process ξ(t = ξ μ,α,σ,θ (t is [ E e ikξ(t] = exp [ ikμ σ α k α { iθ sgn(k tan(πα/} ]. An elementary calculation (e.g., see [3, p. ] shows that the process X γ (t has stability index α = /γ, skewness θ =, scale σ α = tcos(πα/ >, and centering constant μ =. Hence, X γ (t is a spectrally negative stable Lévy process, with no positive jumps. Use the elementary formula (e.g., see [8, Eq. 5.5] (ik α = k α cos(πα/( + i sgn(k tan(πα/ to write E [ e ikx γ (t ] = e t(ikα and then set k = is to see that [ ] E e sx γ (t = e tsα for all s and t. Now, it follows from [7, Theorem, p. 89] that the first-passage time process D u = inf{t : X γ (t >u} is a stable subordinator with Laplace transform E [ e sd ] u = e usβ for all u and s, where the stability index β = /α = γ/. Then, the inverse β- stable subordinator E t in (3. is the generalized inverse of X γ (t, which equals the supremum of X γ (t. Hence, we have E t = Y γ (t pathwise, see also Proposition in [9]. Then, the form of the solution follows from (3.5. The integral form of the fractional wave equation (.3 corresponds to the differential form (. withthe initial conditions p(x, = φ( and ψ = β p(x, t, x β β = γ/. The first initial condition follows directly from (.3. As for the second one, note that Ɣ(γ (t s γ x p(x, sds = I γ t x p(x, t and apply the Caputo derivative to both sides of (.3 to get β β p(x, t = ψ(x + β β Iγ t x p(x, t. Since Caputo and Riemann Liouville derivatives of order <β< are related by (e.g., see [8, p. 39] β t β β f (t = Dβ t f (x f ( Ɣ( β, and I γ t x p(x, t evaluated at t = is zero, we have β β Iγ t x p(x, t = D β t I γ t x p(x, t = I β t x p(x, t. Since p(x, t B, (R R +, x p(x, s Ae B( x +t for some constants A, B > and s t. Therefore, I β t x p(x, t AeB( x +t Ɣ(β (t s β ds = AeB( x +t Ɣ(β + tβ as t. Thus, the initial conditions corresponding to (.3 are p(x, = φ(x and β p(x, = ψ(x. β
4 M. M. Meerschaert et al. In the remainder of this section, we discuss related results in the literature and give some alternative stochastic representations of the solution. Remark 3. Define the reflected stable process Z t = X γ (t inf{x γ (s : s t}, (3.6 where X γ (t is the spectrally negative stable process with index < γ and characteristic function (3.3. Apply [4, Lemma 4.5] to see that Z t has the same one-dimensional distributions as the inverse (3. of a standard β-stable subordinator with β = γ/. Then, it follows from Theorem 3. that for any φ, ψ such that φ, B (R, the unique solution to the fractional wave equation (. inb, (R R + with p(x, = φ(x and β p(x, = ψ(x is given by β the formula p(x, t = E [φ(x + Z t + φ(x Z t ] + x+z E ψ(y dy. (3.7 x Z t The advantage to this representation is that Z t is a Markov process. Remark 3.3 Mainardi [5, Sect. 6.3] considers a version of (. that employs the Caputo fractional derivative γ p(x, t = t Ɣ( γ u p(x, u(t u γ du (3.8 of order <γ <. Mainardi [5, Sect. 6.4] derives the fractional wave equation (. from a viscoelastic model with a power-law stress strain relationship. He notes that the Green s function solution to the fractional wave equation (. can be also expressed in terms of stable densities. He considers the fractional wave equation (. subject to the initial conditions p(x, = δ(x and p(x, =. Since the Caputo and Riemann-Liouville fractional derivatives of order <γ < are related by (e.g., see [5, p. ] t γ γ f (t = Dγ t f (t f ( Ɣ( γ t γ f ( Ɣ( γ, (3.9 when p(x, = Eq.(.6 with initial conditions p(x, = φ(x and β p(x, = has the same β integral form as Eq. (. with p(x, = φ(x and p(x, =. Letting L η α(x be the stable probability density function with characteristic function e ikx L η α (xdx = exp [ k α e i sgn(kηπ/] Mainardi shows that L α α (x = α /α(x for any x R and <α, where β (z = n= ( z n, n!ɣ( nβ β <β< (3. is the Wright function. It follows [5, Eq. 6.37] that the solution to the fractional wave equation (. with initial conditions p(x, = δ(x and p(x, = can be written in the form p(x, t = ( x βt β L(/β /β t β with β = γ/. Hence, the solution to the fractional wave equation (. with initial conditions p(x, = φ(x and p(x, = is given by p(x, t = [ ] φ(x u + φ(x + u βt β L(/β /β ( u t β du. (3. The solution (3. to the fractional wave equation (. involves a stable density with index α = /β = /γ, whereas the solution in Theorem 3. uses the inverse of a stable law with index β = γ/. This can be explained using the Zolotarev duality formula for stable densities [3, Theorem.]. Baeumer et al. [3, Theorem 4.]use Zolotarev duality to prove that βh(u, t = q(u, t for all t > and u, where h(u, t is the density (3. of the standard inverse β-stable subordinator on the set u, q(u, t is the density of the spectrally negative stable process X γ (t with index <γ onthe set < u <, and α = /β. Foru >, the self-similarity argument shows that the function q(u, t = t β L(/β /β ( u t β.
5 Stochastic solutions for fractional wave equations Then, the solution (3. reduces to a special case of (.7 with ψ(x = since p(x, t = = = [φ(x u + φ(x + u] ( u βt β L(/β /β t β du q(u, t [φ(x u + φ(x + u] du β [φ(x u + φ(x + u] h(u, tdu. Theorem 4. in [3] also shows that the conditional distribution of X γ (t given X γ (t > is identical to the distribution of E t. Hence, for any φ, ψ such that φ, B (R, the unique solution to the fractional wave equation (.6 inb, (R R + with p(x, = φ(x and β p(x, = ψ(x can be written as β p(x, t = E φ(x + X γ (t + φ(x X γ (t + x+x γ (t x X γ (t ψ(ydy X γ (t >. (3. An extension of the well-known D. André reflection principle (see Appendix shows that P[Y γ (t x] = P[X γ (t x X γ (t ], and this together with (3.5 gives another proof of (3.. 4 General wave equations Given a closed operator L on a Banach space X of functions, consider the second-order Cauchy problem p(x, t = Lp(x, t; p(x, = φ(x, p(x, = ψ(x. (4. The traditional wave equation (. is a special case where L = x. Bajlekova [5,6] developed the theory of fractional order Cauchy problems p(x, t = Lp(x, t; γ p(x, = φ(x, p(x, = (4. using a Caputo fractional derivative of order <γ <. The general theory of second-order Cauchy problems is laid out in [, Sects ]. A strongly continuous (i.e., continuous in the Banach space norm family of linear operators (Cos(t t is called a cosine family if Cos( = I and Cos(t Cos(s = Cos(t + s + Cos(t s for all s, t. The generator L of the cosine family is defined by Lf(x = lim [Cos(t f (x f (x], t and the domain Dom(L is the set of functions f X for which this limit exists strongly. If the operator L in (4. is a generator of a cosine family (Cos(t t, then for any φ,ψ X, the unique mild solution to the second-order Cauchy problem (4. is given by t p(x, t = Cos(tφ(x + That is, we have p(x, t = φ(x + tψ(x + L Cos(sψ(xds. (4.3 (t sp(x, sds forall t, the integrated version of the second-order Cauchy problem. Furthermore, (4.3 is the unique (classical solution to the second-order Cauchy problem (4. for any φ,ψ Dom(L [, Theorem 3.4.]. Theorem 4. Suppose that the operator L in (4. is a generator of a cosine family (Cos(t t. Then for any φ X, the unique mild solution to the fractional Cauchy problem (4. is given by the formula p(x, t = E [Cos(E t φ(x], (4.4 where Cos(tφ(x is the unique mild solution to the second-order Cauchy problem (4. with ψ =, and E t is the inverse (3. of the standard stable subordinator with index β = γ/. Furthermore, equation (4.4 gives the unique classical solution to (4. for any φ Dom(L. Proof Bajlekova [5, Theorem 3.] proves that if Cos(t φ(x = S (tφ(x solves the second-order Cauchy problem (4. with ψ =, then the unique solution to the fractional Cauchy problem (4. isp(x, t = S γ (tφ(x where the family of solution operators S γ (t is given by the subordination formula S γ (t = S (st γ/ γ/ (st γ/ ds, (4.5
6 M. M. Meerschaert et al. using the Wright function defined in (3.. Recall the identity (e.g., see [5, Eq..3] e zt β (zdz = E β (t (4.6 where the Mittag-Leffler function E β (z = n= z n Ɣ( + βn for β> and z C. Bingham [8] and Bondesson et al. [] show that the inverse E t of a β-stable subordinator has a Mittag-Leffler distribution with E(e se t = e su h(u, tdu = E β ( st β. But, it follows from (4.6 along with a substitution z = u/t β that we also have e su t β β ( u t β du = e sztβ β (z dz = E β ( st β where β = γ/. Then, it follows from the uniqueness of the Laplace transform that the standard inverse β- stable density (3. is related to the Wright function by h(u, t = ( u t β β t β. (4.7 Hence, Bajlekova s solution (4.5 to the fractional wave equation is equivalent to the formula (4.4. Remark 4. Mainardi [5, Sect. 6.3] shows that the solution to the fractional wave equation (. with initial conditions p(x, = φ(x and p(x, = is given by the convolution formula p(x, t = ( u [φ(x u + φ(x +u] t β β t β du. Using (4.7, this reduces to (.7 with ψ(x. Example 4.3 Given an open subset D of R d, consider the Laplace operator L = x on L (D with Dirichlet boundary conditions [, Example 7..]. For any φ Dom(L there exists a unique solution p(x, t to the wave equation p(x, t = x p(x, t; p(x, = φ(x; p(x, = ; p(x, t = x / D by [, Theorem 7..]. Then, it follows from Theorem 4. that the function p γ (x, t = E[p(x, E t ] solves the corresponding fractional wave equation γ p(x, t = x p(x, t; p(x, = φ(x; p(x, = ; p(x, t = x / D on this bounded domain for <γ <, where E t is the inverse stable subordinator (3. with index β = γ/. Example 4.4 If L = B, where B is a generator of a C -semigroup (A(t t on a Banach space of functions, then L is a generator of a cosine family given by Cos(t = (A(t + A( t, t R, see [, Example 3.4.5]. When B = x, (A(t t is a shift semigroup, and equation (4. becomes the traditional wave equation (.. Equation (4.3 giving the solution becomes the d Alembert formula (.. Theorem 4. gives the solution to the fractional wave equation (. with the initial conditions p(x, =, p(x, =. Example 4.5 If L is a self-adjoint linear operator on some Hilbert space such that (Lx, x H ω x H for some ω>and all x Dom(L, then L generates a cosine family [, Example 3.4.6], and hence, (4.3 is the unique classical solution to the wave equation p(x, t = Lp(x, t; p(x, = φ(x; p(x, t = for any φ Dom(L. Then, Theorem 4. implies that the function p γ (x, t = E[p(x, E t ] solves the corresponding fractional wave equation p(x, t = Lp(x, t; γ p(x, = φ(x; p(x, t =, where E t is the inverse stable subordinator (3. with index β = γ/.
7 Stochastic solutions for fractional wave equations 5 Random walk models In this section, we will develop a random walk model for the fractional wave equation (.. First, we decompose the fractional wave equation into simpler parts. Using the notation (.4 for the Riemann Liouville fractional integral, the integral form (.3 of the fractional wave equation with ψ can be written as p(x, t = φ(x + I γ t x p(x, t. (5. Using the property D γ t = D n t In γ t for the Riemann Liouville factional derivative and integral, and the semigroup property I α t Iβ t = I α+β t, it follows easily that D γ t I γ t f (t = f (t [5, Theorem.5]. Apply the operator D γ t to both sides of (5. to get the equivalent form D γ t p(x, t = x p(x, t + D γ t φ(x. (5. An easy computation [8, Example.8] shows that D γ t = t γ /Ɣ( γ, and then, (5. becomes t γ D γ t p(x, t φ(x Ɣ( γ = x p(x, t. (5.3 Since the Caputo and Riemann Liouville fractional derivatives of order <γ <are related by (3.9, Eq. (5.3 is equivalent to the fractional wave equation (. with initial conditions p(x, = φ(x and p(x, =. Now, consider the one way fractional wave equations D γ/ t p + (x, t = x p + (x, t + φ(x t γ/ Ɣ( γ/ D γ/ t p (x, t = x p (x, t + φ(x t γ/ Ɣ( γ/ (5.4 where again <γ <. Apply the operator I γ/ t to both sides to obtain the integral forms (I + I γ/ t x p + (x, t = φ(x (I I γ/ t x p (x, t = φ(x (5.5 where I is the identity operator. Theorem 5. For any φ B (R, the unique solutions to the one way fractional wave equations (5.5 in B, (R R + are given by the formulae p + (x, t = E [φ(x E t] p (x, t = E [φ(x + E t] (5.6 where E t is the generalized inverse (3. of the standard stable subordinator with index β = γ/. Furthermore, the unique solution to the fractional wave equation (.3 in B, (R R + with p(x, = φ(x and p(x, = ψ(x is then given by p(x, t = p + (x, t + p (x, t. Proof Fujita [3] proves the same result with E t replaced by the supremum process (3.4. As noted in the proof of Theorem 3., these two processes have the same one-dimensional distributions. Then, the result follows. Remark 5. A direct proof of Theorem 5. uses an idea from Fujita [3]. Apply [7, Theorem 4.] to see that the density (3. of the inverse stable subordinator with index β = γ/ solves equation D γ/ t t γ/ h(x, t = x h(x, t + δ(x Ɣ( γ/. (5.7 It follows using the principle of superposition that p + (x, t = E [φ(x E t] = φ(x uh(u, tdy solves the positive one way fractional wave equation. Then, a simple change of coordinates shows that p (x, t = E [φ(x + E t ] / solves the negative one way fractional wave equation. Now, write (I I γ t x (p + + p = (I I γ/ t x (I + I γ/ t x p + + (I + I γ/ t x (I I γ/ t x p = (I I γ/ t x φ + (I + Iγ/ t x φ = φ which is equivalent to the fractional wave equation (.3 with p(x, = φ(x and p(x, =. One can also prove Theorem 3. in the same manner. Just apply the same argument again with φ(x replaced by the function (x = x ψ(ydy, and then add the two solutions. Remark 5.3 Here, we indicate an alternative proof of Theorem 4. using Riemann Liouville fractional derivatives and an idea from [9]. Suppose that p(x, t solves the second-order Cauchy problem (4. withinitial conditions p(x, = φ(x and p(x, =. Let h(u, t be the density (3. of the inverse stable subordinator with index β = γ/. From (5.7, it follows that D β t h(u, t = u h(u, t on t > and u >. It follows
8 M. M. Meerschaert et al. from (3. and the asymptotic behavior of stable densities that h(+, t = t β /Ɣ( β and h(, t = for all t > (e.g., see [, Eq. ]. Write p γ (x, t = p(x, uh(u, tdu and integrate by parts to get D β t p γ (x, t = = p(x, ud β t h(u, tdu p(x, u [ u h(u, t] du = p(x, h(+, t + u p(x, uh(u, tdu. Integrate by parts again and use (4. to get D β t p γ (x, t = p(x, D β t h(+, t + u p(x, u [ u h(u, t] du = p(x, D β t h(+, t + u p(x, uh(u, tdu t β = Lp γ (x, t + φ(x Ɣ( β using the general formula D β t [t α ]=Ɣ( + αt α β /Ɣ ( + α β [8, Example.7]. This equation for p γ (x, t is equivalent to the second-order Cauchy problem (4. with initial conditions p γ (x, = φ(x and p γ (x, =, since the Caputo and Riemann Liouville fractional derivatives of order γ = β (, are related by (3.9. Finally, we develop a simple particle tracking method for solving the fractional wave equation (.3, using a continuous-time random walk [6,7] that converges to the stochastic solution of the fractional wave equation. The main idea is to construct a random walk model that converges to the inverse stable subordinator E t and use the fact that the density (3. ofe t solves the positive one way fractional wave equation. Theorem 5.4 Given a continuous probability density function φ(x on R, letx be a random variable with density φ(x. Let X be a Bernoulli random variable independent of X such that P[X = ] =P[X = ] =/, setx n = X for n >, and let S(n = X + +X n = nx. Let W n be iid random variables independent of X, X with P[W n > t] =Ct β for t > C /β, where <β<and C = /Ɣ( β. Let T =, T n = W + + W n for n, and N t = max{n : T n t} for t. Then, X + c β S(N ct U t as c (5.8 in D[, with the Skorokhod J topology, where the random variable U t has density p(x, t = E [φ(x + E t + φ(x E t ], (5.9 and E t is the inverse stable subordinator (3. with index β = γ/. Hence, p(x, t is the unique solution to the fractional wave equation (.3 in B, (R R + with p(x, = φ(x and p(x, = ψ(x. Proof For any c >, it follows from [8, Theorem 3.4 and Eq. 4.9] that c /β T [ct] D t as c in J topology in D[,, where D t is a β-stable subordinator with characteristic function E[e ikd t ]= exp[ tcɣ( β( ik β ]. Taking C = /Ɣ( β, the limit is a standard stable subordinator with Laplace transform E[e sd t ]=e tsβ, and then, [6, Theorem 3. and Corollary 3.4] implies that c β N ct E t as c, where E t is the inverse (3. of the standard stable subordinator D t. Since S(n = nx it follows easily that X + c β S(N ct = X + c β N ct X X + E t X as c. Then, (5.8 holds with U t = X + E t X, and a simple conditioning argument yields (5.9. Then, Theorem 3. shows that p(x, t is the unique solution to the fractional wave equation (.3 inb, (R R + with p(x, = φ(x and p(x, = ψ(x. Theorem 5.4 provides a physical model for the fractional wave equation. Each sample path represents a packet of wave energy moving out from its initial position X at unit speed, represented by the process S(n. For the traditional wave equation, this is the correct particle model. In the fractional case, time
9 Stochastic solutions for fractional wave equations delays with a power-law probability distribution occur between movements, and this retards the progress of the wave outward from the starting point. These delays are related to the heterogeneous structure of the conducting medium, see Mainardi [5, Sect. 6.4]. Remark 5.5 Theorem 5.4 implies that the histogram of a large number M of identical continuous-time randomwalk processes X + c β S(N ct gives an approximate solution to the fractional wave equation, which gains accuracy at M and c.itisasimple matter to simulate the waiting times using the formula W n = (U n /C /β where U n are iid uniform random variables on (,. Theorem 5.4 remains true for any iid waiting times W n > in the domain of attraction of the β-stable subordinator, except that the norming constants c β need to be adjusted as in [6, Theorem 3.]. Acknowledgments This research was partially supported by NIH grant R-EB79. Appendix Reflection principle The goal of this appendix is to establish the following variation of the D. André reflection principle for Brownian motion, which may also be useful in other contexts. Since this extension is not completely standard, we include its simple proof. Theorem 6. (Reflection principle Suppose that Y t is a Lévy process started at the origin, with no positive jumps, and let S t = sup{y u : u t}. Assume that P(Y t > = P(Y, for all t >. Then, P(S t x = P(Y t x Y t (6. = P(Y t x P(Y t for all t, x >. P ( τ x t, Y t < Y τx = P ( τx t, Y t Y τx < (6. = P (τ x t P (Y < = P (Y < P (Y P (τ x t P (Y = P (Y < P (Y P ( τ x t, Y t Y τx. Observe that we have {τ x < t} {S t > x} {τ x t} {S t x} for all t and x >. Therefore, P(S t > x P(τ x t = P(τ x t, Y t Y τx + P(τ x t, Y t < Y τx ( = + P (Y < P(τ x t, Y t Y τx P (Y = P (Y P(τ x t, Y t Y τx = P (Y t P(τ x t, Y t Y τx, Y t = P(τ x t, Y t Y τx Y t. On the other hand, and with similar arguments, P(S t > x P(τ x < t = P(τ x < t, Y t Y τx + P(τ x < t, Y t < Y τx = P(τ x < t, Y t Y τx Y t P(τ x < t, Y t > Y τx Y t. Therefore, P(τ x < t, Y t > Y τx Y t P(S t > x P(τ x t, Y t Y τx Y t. Since Y t has no upward jumps, Y τx = x and {τ x < t} {Y t > x} ={Y t > x}. Therefore, P(Y t > x Y t P(S t > x P(Y t x Y t. We can now use standard approximation techniques: { } X > x /n Proof Let τ x := inf{u > : Y u > x} denote the first-passage time process. Since (Y t t has stationary independent increments, it follows that (Y t+τx Y τx t is a Lévy process, which is independent of the σ - algebra generated by (Y t t τx, and it has the same finitedimensional distributions as (Y t t. Consequently, {X x} = n = n and {X > x} = n { X x /n } {X > x + /n}
10 M. M. Meerschaert et al. to get P(Y t x Y t = lim P(Y t > x /n Y t n lim P(S t > x /n n = P(S t x and P(S t x = lim P(S t > x /n n lim P(Y t x /n Y t n = P(Y t x Y t which proves P(S t x = P(Y t x Y t. Remark 6. If (Y t t is a Brownian motion, then (6. becomes the classical reflection principle: P(Y t = /, so that (6. is equivalent to P(S t x = P(Y t x. The proof of Theorem 6. relies essentially on local symmetry and the strong Markov property. Let Y t be a strong Markov process with càdlàg paths and transition function p t (z, dy = P z (Y t dy. Write τx z = inf{u > : Y u z > x} for the first passage time above the level x + z for the process Y t started at z; observe that, in general, Y τ x + z. We can use the strong Markov property in (6. to get for any starting point z P z ( τx z t, Y t < Y τ = P z ( τx z t, Y t Y τ < = P Y τx z (ω( Y t τ (ω Y < P z (dω. {τ t} If we assume, in addition, some local symmetry, i.e., that for some constant c (, we have P z (Y t z < P z = c for all t >, z R, (6.3 (Y t z then we get P Y τ (ω( Y t τ (ω Y < = c P Y τ (ω ( Yt τ (ω Y and, with a similar argument, P z ( τ t, Y t < Y τ = c P z ( τ t, Y t Y τ. This means that we can follow the lines of the proof of Theorem 6. to derive the following general result. Theorem 6.3 (Markov reflection principle Suppose (Y t, P z is a strong Markov process satisfying the local symmetry condition (6.3. Set S t = sup{y u Y : u t}. Then, we have for all t, x > and z R P z (S t > x P z (S t x, Y t Y τ Y t z (6.4 P z (Y t z x Y t z. If Y t has only non-positive jumps, then Y τ = x + z a.s., and we get for all t, x > and z R P z (S t x = P z (Y t z x Y t z. (6.5 Remark 6.4 Itisalsopossibletoprove (6. using relation (3.6 in Alili and Chaumont []. References. Alili, L., Chaumont, L.: A new fluctuation identity for Lévy processes and some applications. Bernoulli 7, (. Arendt, W., Batty, C., Hieber, M., Neubrander, F.: Vector- Valued Laplace Transforms and Cauchy Problems. Monographs in Mathematics. Birkhaeuser, Berlin ( 3. Baeumer, B., Meerschaert, M.M., Nane, E.: Space-time duality for fractional diffusion. J. Appl. Probab. 46, 5 (9 4. Baeumer, B., Kovács, M., Meerschaert, M.M., Schilling, R.L., Straka, P.: Reflected spectrally negative stable proc esses and their governing equations. Trans. Am. Math. Soc. (3 5. Bajlekova, E.G.: Fractional evolution equations in Banach spaces. Ph.D. thesis, Eindhoven University of Technology ( 6. Bazhlekova, E.: Subordination principle for fractional evolution equations. Fract. Calc. Appl. Anal. 3(3, 3 3 ( 7. Bertoin, J.: Lévy Processes. Cambridge University Press, Cambridge ( Bingham, N.H.: Limit theorems for occupation times of Markov processes. Z. Warsch. Verw. Geb. 7, (97 9. Bingham, N.H.: Maxima of sums of random variables and suprema of stable processes. Z. Wahrsch. Verw. Geb. 6, (973. Bondesson, L., Kristiansen, G., Steutel, F.: Infinite divisibility of random variables and their integer parts. Stat. Probab. Lett. 8, 7 78 (996. Chen, W., Holm, S.: Fractional Laplacian time-space models for linear and nonlinear lossy media exhibiting arbitrary frequency power-law dependency. J. Acoust. Soc. Am. 5, (4. Duck, F.A.: Physical Properties of Tissue: A Comprehensive Reference Book. Academic Press, Boston (99
11 Stochastic solutions for fractional wave equations 3. Fujita, Y.: Integrodifferential equation which interpolates the heat equation and the wave equation. Osaka J. Math. 7(, 39 3 (99 4. Kelly, J.F., McGough, R.J., Meerschaert, M.M.: Analytical time-domain Green s functions for power-law media. J. Acoust. Soc. Am. 4(5, (8 5. Mainardi, F.: Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models. World Scientific, Singapore ( 6. Meerschaert, M.M., Scheffler, H.-P.: Limit theorems for continuous time random walks with infinite mean waiting times. J. Appl. Probab. 4, (4 7. Meerschaert, M.M., Scheffler, H.-P.: Triangular array limits for continuous time random walks. Stoch. Proc. Appl. 8, (8 8. Meerschaert, M.M., Sikorskii, A.: Stochastic Models for Fractional Calculus. De Gruyter, Berlin ( 9. Meerschaert, M.M., Straka, P., Zhou, Y., McGough, R.J.: Stochastic solution to a time-fractional attenuated wave equation. Nonlinear Dyn. 7(, 73 8 (. Meerschaert, M.M., Straka, P.: Inverse stable subordinators. Math. Model. Nat. Phenom. 8(, 6 (3. Samorodnitsky, G., Taqqu, M.S.: Stable Non-Gaussian Processes: Stochastic Models with Infinite Variance. Chapman and Hall, New York (994. Szabo, T.L.: Time domain wave equations for lossy media obeying a frequency power law. J. Acoust. Soc. Am. 96, 49 5 ( Treeby, B.E., Cox, B.T.: Modeling power law absorption and dispersion for acoustic propagation using the fractional Laplacian. J. Acoust. Soc. Am. 7, (
Covariance structure of continuous time random walk limit processes
Covariance structure of continuous time random walk limit processes Alla Sikorskii Department of Statistics and Probability Michigan State University Fractional Calculus, Probability and Non-local Operators:
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationCTRW Limits: Governing Equations and Fractal Dimensions
CTRW Limits: Governing Equations and Fractal Dimensions Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 19-23, 2013 Joint work with Z-Q. Chen, M. D Ovidio, M.M. Meerschaert,
More informationReflected spectrally negative stable processes and fractional Cauchy problems
Reflected spectrally negative stable processes and fractional Cauchy problems 7th International Conference on Lévy Processes: Theory and Applications 15 19 July 2013 Mark M. Meerschaert Department of Statistics
More informationCorrelation Structure of Fractional Pearson Diffusions
Correlation Structure of Fractional Pearson Diffusions Alla Sikorskii Department of Statistics and Probability Michigan State University International Symposium on Fractional PDEs: Theory, Numerics and
More informationON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP
Dedicated to Professor Gheorghe Bucur on the occasion of his 7th birthday ON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP EMIL POPESCU Starting from the usual Cauchy problem, we give
More informationContinuous Time Random Walk Limits in Bounded Domains
Continuous Time Random Walk Limits in Bounded Domains Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY May 16, 2012 Acknowledgement Boris Baeumer, University of Otago, New Zealand.
More informationLogarithmic scaling of planar random walk s local times
Logarithmic scaling of planar random walk s local times Péter Nándori * and Zeyu Shen ** * Department of Mathematics, University of Maryland ** Courant Institute, New York University October 9, 2015 Abstract
More informationApplied Mathematics Letters
Applied Mathematics Letters 24 (211) 219 223 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Laplace transform and fractional differential
More informationCauchy Problems in Bounded Domains and Iterated Processes
Cauchy Problems in Bounded Domains and Iterated Processes Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 18, 21 Continuous time random walks R d X 2 X 3 X 1 J 1 J 2 J 3 J
More informationCauchy Problems Solved by Running Subordinate Processes
Cauchy Problems Solved by Running Subordinate Processes Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 18, 21 Continuous time random walks R d X 2 X 3 X 1 J 1 J 2 J 3 J 4
More informationFunctional central limit theorem for super α-stable processes
874 Science in China Ser. A Mathematics 24 Vol. 47 No. 6 874 881 Functional central it theorem for super α-stable processes HONG Wenming Department of Mathematics, Beijing Normal University, Beijing 1875,
More informationFractional Cauchy Problems for time-changed Processes
Fractional Cauchy Problems for time-changed Processes Erkan Nane Department of Mathematics and Statistics Auburn University February 17, 211 Erkan Nane (Auburn University) Fractional Cauchy Problems February
More informationSpace-Time Duality and Medical Ultrasound
Space-Time Duality and Medical Ultrasound James F. Kelly and Mark M. Meerschaert A Workshop on Future Directions in Fractional Calculus Research and Applications Michigan State University October 21, 2016
More informationPotentials of stable processes
Potentials of stable processes A. E. Kyprianou A. R. Watson 5th December 23 arxiv:32.222v [math.pr] 4 Dec 23 Abstract. For a stable process, we give an explicit formula for the potential measure of the
More information1. Stochastic Processes and filtrations
1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S
More informationSemi-Markov approach to continuous time random walk limit processes
Semi-Markov approach to continuous time random walk limit processes Joint Mathematics Meetings San Antonio, TX January 13, 2015 Mark M. Meerschaert Michigan State University East Lansing, Michigan, USA
More informationConvergence rates in weighted L 1 spaces of kernel density estimators for linear processes
Alea 4, 117 129 (2008) Convergence rates in weighted L 1 spaces of kernel density estimators for linear processes Anton Schick and Wolfgang Wefelmeyer Anton Schick, Department of Mathematical Sciences,
More informationHOPF S DECOMPOSITION AND RECURRENT SEMIGROUPS. Josef Teichmann
HOPF S DECOMPOSITION AND RECURRENT SEMIGROUPS Josef Teichmann Abstract. Some results of ergodic theory are generalized in the setting of Banach lattices, namely Hopf s maximal ergodic inequality and the
More informationThe Codimension of the Zeros of a Stable Process in Random Scenery
The Codimension of the Zeros of a Stable Process in Random Scenery Davar Khoshnevisan The University of Utah, Department of Mathematics Salt Lake City, UT 84105 0090, U.S.A. davar@math.utah.edu http://www.math.utah.edu/~davar
More informationCORRELATED CONTINUOUS TIME RANDOM WALKS
CORRELATED CONTINUOUS TIME RANDOM WALKS MARK M. MEERSCHAERT, ERKAN NANE, AND YIMIN XIAO Abstract. Continuous time random walks impose a random waiting time before each particle jump. Scaling limits of
More informationExistence Results for Multivalued Semilinear Functional Differential Equations
E extracta mathematicae Vol. 18, Núm. 1, 1 12 (23) Existence Results for Multivalued Semilinear Functional Differential Equations M. Benchohra, S.K. Ntouyas Department of Mathematics, University of Sidi
More informationMultivariate Generalized Ornstein-Uhlenbeck Processes
Multivariate Generalized Ornstein-Uhlenbeck Processes Anita Behme TU München Alexander Lindner TU Braunschweig 7th International Conference on Lévy Processes: Theory and Applications Wroclaw, July 15 19,
More informationA FEW REMARKS ON THE SUPREMUM OF STABLE PROCESSES P. PATIE
A FEW REMARKS ON THE SUPREMUM OF STABLE PROCESSES P. PATIE Abstract. In [1], Bernyk et al. offer a power series and an integral representation for the density of S 1, the maximum up to time 1, of a regular
More informationLower Tail Probabilities and Related Problems
Lower Tail Probabilities and Related Problems Qi-Man Shao National University of Singapore and University of Oregon qmshao@darkwing.uoregon.edu . Lower Tail Probabilities Let {X t, t T } be a real valued
More informationLecture 17 Brownian motion as a Markov process
Lecture 17: Brownian motion as a Markov process 1 of 14 Course: Theory of Probability II Term: Spring 2015 Instructor: Gordan Zitkovic Lecture 17 Brownian motion as a Markov process Brownian motion is
More informationGARCH processes continuous counterparts (Part 2)
GARCH processes continuous counterparts (Part 2) Alexander Lindner Centre of Mathematical Sciences Technical University of Munich D 85747 Garching Germany lindner@ma.tum.de http://www-m1.ma.tum.de/m4/pers/lindner/
More informationTHE FUNDAMENTAL SOLUTION OF THE DISTRIBUTED ORDER FRACTIONAL WAVE EQUATION IN ONE SPACE DIMENSION IS A PROBABILITY DENSITY
THE FUNDAMENTAL SOLUTION OF THE DISTRIBUTED ORDER FRACTIONAL WAVE EQUATION IN ONE SPACE DIMENSION IS A PROBABILITY DENSITY RUDOLF GORENFLO Free University of Berlin Germany Email:gorenflo@mi.fu.berlin.de
More informationThe Central Limit Theorem: More of the Story
The Central Limit Theorem: More of the Story Steven Janke November 2015 Steven Janke (Seminar) The Central Limit Theorem:More of the Story November 2015 1 / 33 Central Limit Theorem Theorem (Central Limit
More information{σ x >t}p x. (σ x >t)=e at.
3.11. EXERCISES 121 3.11 Exercises Exercise 3.1 Consider the Ornstein Uhlenbeck process in example 3.1.7(B). Show that the defined process is a Markov process which converges in distribution to an N(0,σ
More informationCOMPOSITION SEMIGROUPS AND RANDOM STABILITY. By John Bunge Cornell University
The Annals of Probability 1996, Vol. 24, No. 3, 1476 1489 COMPOSITION SEMIGROUPS AND RANDOM STABILITY By John Bunge Cornell University A random variable X is N-divisible if it can be decomposed into a
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationAn Introduction to Malliavin calculus and its applications
An Introduction to Malliavin calculus and its applications Lecture 3: Clark-Ocone formula David Nualart Department of Mathematics Kansas University University of Wyoming Summer School 214 David Nualart
More informationSmall ball probabilities and metric entropy
Small ball probabilities and metric entropy Frank Aurzada, TU Berlin Sydney, February 2012 MCQMC Outline 1 Small ball probabilities vs. metric entropy 2 Connection to other questions 3 Recent results for
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationLévy Processes and Infinitely Divisible Measures in the Dual of afebruary Nuclear2017 Space 1 / 32
Lévy Processes and Infinitely Divisible Measures in the Dual of a Nuclear Space David Applebaum School of Mathematics and Statistics, University of Sheffield, UK Talk at "Workshop on Infinite Dimensional
More informationOn Optimal Stopping Problems with Power Function of Lévy Processes
On Optimal Stopping Problems with Power Function of Lévy Processes Budhi Arta Surya Department of Mathematics University of Utrecht 31 August 2006 This talk is based on the joint paper with A.E. Kyprianou:
More informationStatistics and Probability Letters
Statistics and Probability Letters 83 (23) 83 93 Contents lists available at SciVerse ScienceDirect Statistics and Probability Letters journal homepage: www.elsevier.com/locate/stapro Fractal dimension
More informationShifting processes with cyclically exchangeable increments at random
Shifting processes with cyclically exchangeable increments at random Gerónimo URIBE BRAVO (Collaboration with Loïc CHAUMONT) Instituto de Matemáticas UNAM Universidad Nacional Autónoma de México www.matem.unam.mx/geronimo
More informationRegular Variation and Extreme Events for Stochastic Processes
1 Regular Variation and Extreme Events for Stochastic Processes FILIP LINDSKOG Royal Institute of Technology, Stockholm 2005 based on joint work with Henrik Hult www.math.kth.se/ lindskog 2 Extremes for
More informationA Note on the Central Limit Theorem for a Class of Linear Systems 1
A Note on the Central Limit Theorem for a Class of Linear Systems 1 Contents Yukio Nagahata Department of Mathematics, Graduate School of Engineering Science Osaka University, Toyonaka 560-8531, Japan.
More informationON THE SHAPE OF THE GROUND STATE EIGENFUNCTION FOR STABLE PROCESSES
ON THE SHAPE OF THE GROUND STATE EIGENFUNCTION FOR STABLE PROCESSES RODRIGO BAÑUELOS, TADEUSZ KULCZYCKI, AND PEDRO J. MÉNDEZ-HERNÁNDEZ Abstract. We prove that the ground state eigenfunction for symmetric
More information16. M. Maejima, Some sojourn time problems for strongly dependent Gaussian processes, Z. Wahrscheinlichkeitstheorie verw. Gebiete 57 (1981), 1 14.
Publications 1. M. Maejima, Some limit theorems for renewal processes with non-identically distributed random variables, Keio Engrg. Rep. 24 (1971), 67 83. 2. M. Maejima, A generalization of Blackwell
More informationFractional Laplacian
Fractional Laplacian Grzegorz Karch 5ème Ecole de printemps EDP Non-linéaire Mathématiques et Interactions: modèles non locaux et applications Ecole Supérieure de Technologie d Essaouira, du 27 au 30 Avril
More informationA Smooth Operator, Operated Correctly
Clark Department of Mathematics Bard College at Simon s Rock March 6, 2014 Abstract By looking at familiar smooth functions in new ways, we can make sense of matrix-valued infinite series and operator-valued
More informationOn the quantiles of the Brownian motion and their hitting times.
On the quantiles of the Brownian motion and their hitting times. Angelos Dassios London School of Economics May 23 Abstract The distribution of the α-quantile of a Brownian motion on an interval [, t]
More informationSpectral functions of subordinated Brownian motion
Spectral functions of subordinated Brownian motion M.A. Fahrenwaldt 12 1 Institut für Mathematische Stochastik Leibniz Universität Hannover, Germany 2 EBZ Business School, Bochum, Germany Berlin, 23 October
More informationA NEW PROOF OF THE WIENER HOPF FACTORIZATION VIA BASU S THEOREM
J. Appl. Prob. 49, 876 882 (2012 Printed in England Applied Probability Trust 2012 A NEW PROOF OF THE WIENER HOPF FACTORIZATION VIA BASU S THEOREM BRIAN FRALIX and COLIN GALLAGHER, Clemson University Abstract
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationThe Subexponential Product Convolution of Two Weibull-type Distributions
The Subexponential Product Convolution of Two Weibull-type Distributions Yan Liu School of Mathematics and Statistics Wuhan University Wuhan, Hubei 4372, P.R. China E-mail: yanliu@whu.edu.cn Qihe Tang
More informationEULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS
Qiao, H. Osaka J. Math. 51 (14), 47 66 EULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS HUIJIE QIAO (Received May 6, 11, revised May 1, 1) Abstract In this paper we show
More informationONE DIMENSIONAL MARGINALS OF OPERATOR STABLE LAWS AND THEIR DOMAINS OF ATTRACTION
ONE DIMENSIONAL MARGINALS OF OPERATOR STABLE LAWS AND THEIR DOMAINS OF ATTRACTION Mark M. Meerschaert Department of Mathematics University of Nevada Reno NV 89557 USA mcubed@unr.edu and Hans Peter Scheffler
More informationSOLUTIONS OF SEMILINEAR WAVE EQUATION VIA STOCHASTIC CASCADES
Communications on Stochastic Analysis Vol. 4, No. 3 010) 45-431 Serials Publications www.serialspublications.com SOLUTIONS OF SEMILINEAR WAVE EQUATION VIA STOCHASTIC CASCADES YURI BAKHTIN* AND CARL MUELLER
More informationA slow transient diusion in a drifted stable potential
A slow transient diusion in a drifted stable potential Arvind Singh Université Paris VI Abstract We consider a diusion process X in a random potential V of the form V x = S x δx, where δ is a positive
More informationOPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS
APPLICATIONES MATHEMATICAE 29,4 (22), pp. 387 398 Mariusz Michta (Zielona Góra) OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS Abstract. A martingale problem approach is used first to analyze
More informationSolution for Problem 7.1. We argue by contradiction. If the limit were not infinite, then since τ M (ω) is nondecreasing we would have
362 Problem Hints and Solutions sup g n (ω, t) g(ω, t) sup g(ω, s) g(ω, t) µ n (ω). t T s,t: s t 1/n By the uniform continuity of t g(ω, t) on [, T], one has for each ω that µ n (ω) as n. Two applications
More informationOn Reflecting Brownian Motion with Drift
Proc. Symp. Stoch. Syst. Osaka, 25), ISCIE Kyoto, 26, 1-5) On Reflecting Brownian Motion with Drift Goran Peskir This version: 12 June 26 First version: 1 September 25 Research Report No. 3, 25, Probability
More informationLIST OF MATHEMATICAL PAPERS
LIST OF MATHEMATICAL PAPERS 1961 1999 [1] K. Sato (1961) Integration of the generalized Kolmogorov-Feller backward equations. J. Fac. Sci. Univ. Tokyo, Sect. I, Vol. 9, 13 27. [2] K. Sato, H. Tanaka (1962)
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 22 12/09/2013. Skorokhod Mapping Theorem. Reflected Brownian Motion
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.7J Fall 213 Lecture 22 12/9/213 Skorokhod Mapping Theorem. Reflected Brownian Motion Content. 1. G/G/1 queueing system 2. One dimensional reflection mapping
More informationOn semilinear elliptic equations with measure data
On semilinear elliptic equations with measure data Andrzej Rozkosz (joint work with T. Klimsiak) Nicolaus Copernicus University (Toruń, Poland) Controlled Deterministic and Stochastic Systems Iasi, July
More informationSome Properties of NSFDEs
Chenggui Yuan (Swansea University) Some Properties of NSFDEs 1 / 41 Some Properties of NSFDEs Chenggui Yuan Swansea University Chenggui Yuan (Swansea University) Some Properties of NSFDEs 2 / 41 Outline
More information1D Quintic NLS with White Noise Dispersion
2011 年 8 月 8 日 1D Quintic NLS with White Noise Dispersion Yoshio TSUTSUMI, Kyoto University, Arnaud DEBUSSCHE, ENS de Cachan, Bretagne 1D quintic NLS with white noise dispersion idu + xu 2 dβ(t) =λ u 4
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationAn Operator Theoretical Approach to Nonlocal Differential Equations
An Operator Theoretical Approach to Nonlocal Differential Equations Joshua Lee Padgett Department of Mathematics and Statistics Texas Tech University Analysis Seminar November 27, 2017 Joshua Lee Padgett
More informationNonlocal Fractional Semilinear Delay Differential Equations in Separable Banach spaces
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn:2319-765x. Volume 10, Issue 1 Ver. IV. (Feb. 2014), PP 49-55 Nonlocal Fractional Semilinear Delay Differential Equations in Separable Banach
More informationBrownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory
More informationStarting from Heat Equation
Department of Applied Mathematics National Chiao Tung University Hsin-Chu 30010, TAIWAN 20th August 2009 Analytical Theory of Heat The differential equations of the propagation of heat express the most
More informationA New Fractional Process: A Fractional Non-homogeneous Poisson Process Enrico Scalas
A New Fractional Process: A Fractional Non-homogeneous Poisson Process Enrico Scalas University of Sussex Joint work with Nikolai Leonenko (Cardiff University) and Mailan Trinh Fractional PDEs: Theory,
More informationGLOBAL WELL-POSEDNESS FOR NONLINEAR NONLOCAL CAUCHY PROBLEMS ARISING IN ELASTICITY
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 55, pp. 1 7. ISSN: 1072-6691. UL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu GLOBAL WELL-POSEDNESS FO NONLINEA NONLOCAL
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationA New Wavelet-based Expansion of a Random Process
Columbia International Publishing Journal of Applied Mathematics and Statistics doi:10.7726/jams.2016.1011 esearch Article A New Wavelet-based Expansion of a andom Process Ievgen Turchyn 1* eceived: 28
More informationON ADDITIVE TIME-CHANGES OF FELLER PROCESSES. 1. Introduction
ON ADDITIVE TIME-CHANGES OF FELLER PROCESSES ALEKSANDAR MIJATOVIĆ AND MARTIJN PISTORIUS Abstract. In this note we generalise the Phillips theorem [1] on the subordination of Feller processes by Lévy subordinators
More informationHCM PROPERTY AND THE HALF-CAUCHY DISTRIBUTION PIERRE B O S C H (LILLE 1)
PROBABILITY AND MATHEMATICAL STATISTICS Vol. 35, Fasc. 2 25), pp. 9 2 HCM PROPERTY AND THE HALF-CAUCHY DISTRIBUTION BY PIERRE B O S C H LILLE ) Abstract. Let Z α and Z α be two independent positive α-stable
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationNonlocal problems for the generalized Bagley-Torvik fractional differential equation
Nonlocal problems for the generalized Bagley-Torvik fractional differential equation Svatoslav Staněk Workshop on differential equations Malá Morávka, 28. 5. 212 () s 1 / 32 Overview 1) Introduction 2)
More informationBeyond the color of the noise: what is memory in random phenomena?
Beyond the color of the noise: what is memory in random phenomena? Gennady Samorodnitsky Cornell University September 19, 2014 Randomness means lack of pattern or predictability in events according to
More informationfor all f satisfying E[ f(x) ] <.
. Let (Ω, F, P ) be a probability space and D be a sub-σ-algebra of F. An (H, H)-valued random variable X is independent of D if and only if P ({X Γ} D) = P {X Γ}P (D) for all Γ H and D D. Prove that if
More informationPotential Theory on Wiener space revisited
Potential Theory on Wiener space revisited Michael Röckner (University of Bielefeld) Joint work with Aurel Cornea 1 and Lucian Beznea (Rumanian Academy, Bukarest) CRC 701 and BiBoS-Preprint 1 Aurel tragically
More informationOn the Finite Caputo and Finite Riesz Derivatives
EJTP 3, No. 1 (006) 81 95 Electronic Journal of Theoretical Physics On the Finite Caputo and Finite Riesz Derivatives A. M. A. El-Sayed 1 and M. Gaber 1 Faculty of Science University of Alexandria, Egypt
More informationA NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES
A NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES STEFAN TAPPE Abstract. In a work of van Gaans (25a) stochastic integrals are regarded as L 2 -curves. In Filipović and Tappe (28) we have shown the connection
More information1D Quintic NLS with White Noise Dispersion
2011 年 1 月 7 日 1D Quintic NLS with White Noise Dispersion Yoshio TSUTSUMI, Kyoto University, Arnaud DEBUSSCHE, ENS de Cachan, Bretagne 1D quintic NLS with white noise dispersion idu + xu 2 dβ(t) =λ u 4
More informationALMOST PERIODIC SOLUTIONS OF HIGHER ORDER DIFFERENTIAL EQUATIONS ON HILBERT SPACES
Electronic Journal of Differential Equations, Vol. 21(21, No. 72, pp. 1 12. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu ALMOST PERIODIC SOLUTIONS
More informationConservative Control Systems Described by the Schrödinger Equation
Conservative Control Systems Described by the Schrödinger Equation Salah E. Rebiai Abstract An important subclass of well-posed linear systems is formed by the conservative systems. A conservative system
More informationITÔ S ONE POINT EXTENSIONS OF MARKOV PROCESSES. Masatoshi Fukushima
ON ITÔ S ONE POINT EXTENSIONS OF MARKOV PROCESSES Masatoshi Fukushima Symposium in Honor of Kiyosi Itô: Stocastic Analysis and Its Impact in Mathematics and Science, IMS, NUS July 10, 2008 1 1. Itô s point
More informationWeak Limits for Multivariate Random Sums
Journal of Multivariate Analysis 67, 398413 (1998) Article No. MV981768 Weak Limits for Multivariate Random Sums Tomasz J. Kozubowski The University of Tennessee at Chattanooga and Anna K. Panorska - The
More informationPacking-Dimension Profiles and Fractional Brownian Motion
Under consideration for publication in Math. Proc. Camb. Phil. Soc. 1 Packing-Dimension Profiles and Fractional Brownian Motion By DAVAR KHOSHNEVISAN Department of Mathematics, 155 S. 1400 E., JWB 233,
More informationCitation Osaka Journal of Mathematics. 41(4)
TitleA non quasi-invariance of the Brown Authors Sadasue, Gaku Citation Osaka Journal of Mathematics. 414 Issue 4-1 Date Text Version publisher URL http://hdl.handle.net/1194/1174 DOI Rights Osaka University
More informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More informationHarnack Inequalities and Applications for Stochastic Equations
p. 1/32 Harnack Inequalities and Applications for Stochastic Equations PhD Thesis Defense Shun-Xiang Ouyang Under the Supervision of Prof. Michael Röckner & Prof. Feng-Yu Wang March 6, 29 p. 2/32 Outline
More informationPACKING-DIMENSION PROFILES AND FRACTIONAL BROWNIAN MOTION
PACKING-DIMENSION PROFILES AND FRACTIONAL BROWNIAN MOTION DAVAR KHOSHNEVISAN AND YIMIN XIAO Abstract. In order to compute the packing dimension of orthogonal projections Falconer and Howroyd 997) introduced
More informationLAN property for ergodic jump-diffusion processes with discrete observations
LAN property for ergodic jump-diffusion processes with discrete observations Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Arturo Kohatsu-Higa (Ritsumeikan University, Japan) &
More informationOn the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem
On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem Koichiro TAKAOKA Dept of Applied Physics, Tokyo Institute of Technology Abstract M Yor constructed a family
More informationThe Fractional Laplacian
The Fabian Seoanes Correa University of Puerto Rico, Río Piedras Campus February 28, 2017 F. Seoanes Wave Equation 1/ 15 Motivation During the last ten years it has been an increasing interest in the study
More informationLarge deviations and averaging for systems of slow fast stochastic reaction diffusion equations.
Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint work with Michael Salins 2, Konstantinos Spiliopoulos 3.) 1. Department of Mathematics
More informationarxiv: v1 [math.ap] 26 Mar 2013
Analytic solutions of fractional differential equations by operational methods arxiv:134.156v1 [math.ap] 26 Mar 213 Roberto Garra 1 & Federico Polito 2 (1) Dipartimento di Scienze di Base e Applicate per
More informationSOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD
SOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD R. C. Mittal 1 and Ruchi Nigam 2 1 Department of Mathematics, I.I.T. Roorkee, Roorkee, India-247667. Email: rcmmmfma@iitr.ernet.in
More informationOther properties of M M 1
Other properties of M M 1 Přemysl Bejda premyslbejda@gmail.com 2012 Contents 1 Reflected Lévy Process 2 Time dependent properties of M M 1 3 Waiting times and queue disciplines in M M 1 Contents 1 Reflected
More informationRichard F. Bass Krzysztof Burdzy University of Washington
ON DOMAIN MONOTONICITY OF THE NEUMANN HEAT KERNEL Richard F. Bass Krzysztof Burdzy University of Washington Abstract. Some examples are given of convex domains for which domain monotonicity of the Neumann
More information