Statistics and Probability Letters
|
|
- Audra Warren
- 5 years ago
- Views:
Transcription
1 Statistics and Probability Letters 83 (23) Contents lists available at SciVerse ScienceDirect Statistics and Probability Letters journal homepage: Fractal dimension results for continuous time random walks Mark M. Meerschaert a,, rkan Nane b, Yimin Xiao a a Department of Statistics and Probability, Michigan State University, ast Lansing, MI 48824, United States b Department of Mathematics and Statistics, Auburn University, Auburn, AL 36849, United States a r t i c l e i n f o a b s t r a c t Article history: Received 7 December 22 Received in revised form 3 December 22 Accepted January 23 Available online 5 January 23 Continuous time random walks impose random waiting times between particle jumps. This paper computes the fractal dimensions of their process limits, which represent particle traces in anomalous diffusion. 23 lsevier B.V. All rights reserved. Keywords: Continuous time random walk Random fractal Hausdorff dimension Packing dimension. Introduction Given a sequence of space time random vectors {(J n, W n ), n } on R d [, ), a particle arrives at location S(n) = J + + J n at time T(n) = W + + W n. The renewal process N(t) = max{n : T(n) t} counts the number of jumps, and the continuous time random walk (CTRW) S(N t ) with T() = gives the particle location at time t. Under suitable conditions, the normalized sample paths of the space time random walk (S(n), T(n)) converge to a limit {(Y(u), D(u)), u }, and continuous mapping arguments yield a CTRW scaling limit {Y((t)), t }, where (t) = inf{x : D(u) > t}, t (.) is the inverse of the time process D(u). CTRW limits provide a microscopic model for particle motions, leading to macroscopic differential equation models that can involve fractional derivatives in space and time. See Meerschaert and Sikorskii (22) for details. Sample paths of the CTRW limit process Y((t)) model particle traces in statistical physics. This paper investigates the fractal properties of those particle traces. Our results provide physical insight into the microscopic behavior of individual particles undergoing anomalous diffusion. The process (t) models particle resting times. The graph of the inverse process u = (t) flips the axes on a graph of the time process t = D(u). Thus, when D(u) jumps, (t) remains constant, and the particle rests. Unless the resting periods are exponentially distributed, the CTRW limit Y((t)) is non-markovian, and novel methods are required to compute the Hausdorff and packing dimensions of the range and graph. For example, assume i.i.d. particle jumps J n with finite variance and zero mean, and i.i.d. waiting times P(W n > t) Ct β independent of the jumps for some C > and < β <. Then Theorem 4.2 in Meerschaert and Scheffler (24) implies that the CTRW limit Y((t)) is a Brownian motion with an inverse β-stable time change. A typical sample path, shown Corresponding author. -mail addresses: mcubed@stt.msu.edu (M.M. Meerschaert), nane@auburn.edu (. Nane), xiao@stt.msu.edu (Y. Xiao). URLs: (M.M. Meerschaert), ezn (. Nane), (Y. Xiao) /$ see front matter 23 lsevier B.V. All rights reserved. doi:.6/j.spl.23..
2 84 M.M. Meerschaert et al. / Statistics and Probability Letters 83 (23) x t Fig.. Typical sample path of a CTRW limit process Y((t)). Here Y(t) is a Brownian motion and (t) is the inverse of a.8-stable subordinator. in Fig., resembles a Brownian motion interrupted by long resting periods. This process appears in the theory of random conductance models; see Barlow and Cerný (2). Stochastic processes with locally constant paths are also discussed in Davydov (22). Proposition 2.3 shows that the fractal (Hausdorff or packing) dimension is + β/2, less than the dimension 3/2 of a Brownian motion graph, showing that long resting times also affect the geometry of particle traces during motion. 2. CTRW dimension results First we establish a general result concerning a time-changed stochastic process X(t) = Y((t)) for t. We assume that Y(u) is a stochastic process on R d and (t) is a real-valued stochastic process with () = and nondecreasing continuous sample paths. We emphasize that these two processes are not necessarily independent. We use the standard definitions of Hausdorff and packing measures; e.g., see Falconer (99), Kahane (985), Taylor (986) and Xiao (24). Since (t) is continuous, naturally the range of Y((t)) is the same as the range of Y(t) up to a random time. Hence the next result is intuitively obvious, but we include it here for completeness. Proposition 2.. If () > a.s. and there exist constants c and c 2 such that for all constants < a < dim H Y([, a]) = c, dim P Y([, a]) = c 2 a.s., (2.) then almost surely dim H X([, ]) = c and dim P X([, ]) = c 2. (2.2) Proof. Since the process t (t) is nondecreasing and continuous, the range ([, ]) is the random interval [, ()]. Hence X([, ]) = Y([, ()]). It follows from the σ -stability of dim H and (2.) that dim H Y([, )) = c a.s. Hence dim H X([, ]) c almost surely. On the other hand, (2.) implies P dim H Y([, q]) = c, q Q + =, (2.3) where Q + denotes the set of positive rational numbers. Since () > almost surely, we see that there is an event Ω Ω with P(Ω ) = such that for every ω Ω we have (, ω) > and dim H Y([, q], ω) = c for all q Q +. Since for every ω Ω there is a q Q + such that < q < (, ω), we derive that dim H X([, ], ω) = dim H Y([, (, ω)], ω) dim H Y([, q], ω) = c. Combining the upper and lower bounds for dim H X([, ]) yields the first equation in (2.2). The proof of the second equation in (2.2) is similar and is omitted. The graph of the CTRW limit X is closely related to the range of the space time limit Z = {(Y(u), D(u)), u }. In fact, the graph of X is obtained by connecting the points in the range of Y by horizontal line segments, representing the particle resting periods. This motivates the following result. Theorem 2.2. If () > a.s. and there exist constants c 5 and c 6 such that for all constants < a < dim H Z([, a]) = c 5 and dim P Z([, a]) = c 6 a.s., (2.4)
3 M.M. Meerschaert et al. / Statistics and Probability Letters 83 (23) then dim H GrX([, ]) = max, dim H Z([, ]), a.s. (2.5) and dim P GrX([, ]) = max, dim P Z([, ]), a.s. (2.6) Proof. We only prove (2.5), and the proof of (2.6) is similar. The sample function x D(x) is a.s. strictly increasing and we can write the unit interval [, ] in the state space of D as [, ] = D([, ())) I i, i= where for each i, I i is a subinterval on which (t) is a constant. Using D we can express I i = [D(x i ), D(x i )), which is the gap corresponding to the jumping site x i of D, except in the case when x i = (). In the latter case, I i = [D(x i ), ]. Notice that I i (i ) are disjoint intervals and (t) = (s) if and only if s, t I i for some i. Thus, over each interval I i, the graph of X is a horizontal line segment. More precisely, we can decompose the graph set of X as GrX([, ]) = (t, Y((t))) : t [, ] = (t, Y((t))) : t D([, ())) Hence, by the σ -stability of dim H, we have (2.7) (t, Y((t))) : t Ii. (2.8) i= dim H GrX([, ]) = max, dim H (t, Y((t))) : t D([, ())). (2.9) On the other hand, every t D([, ())) can be written as t = D(x) for some x < () and (t) = (D(x)) = x, we see that (t, Y((t))) : t D([, ())) = (D(x), Y(x)) : x [, ()), a.s. (2.) It follows from (2.4) that P ω : dim H (D(x, ω), Y(x, ω)) : x [, q] = c5, q Q + =. (2.) Combining this with the assumption that (, ω) > almost surely, we can find an event Ω 2 such that P(Ω 2 ) = and for every ω Ω 2 we derive from (2.) that dim H (D(x, ω), Y(x, ω)) : x [, (, ω)) = c5, (2.2) since q < (, ω 2 ) < q 2 for some q, q 2 Q +, and U V implies dim H (U) dim H (V). Combining (2.) and (2.2) yields dim H (t, Y((t))) : t D([, ())) = c5, a.s. (2.3) Therefore, (2.5) follows from (2.9) and (2.3). Now we apply these general results to CTRW limits. We assume that the space time random walk (S(n), T(n)) has a process limit {(Y(u), D(u)), u } such that D(u) is a subordinator (a nondecreasing Lévy process) with D() = and [e sd(u) ] = e uσ (s), where the Laplace exponent σ (s) = ( e sy )ν(dy), and the Lévy measure ν of D σ satisfies ν(, ) =, so that the sample paths x D σ (x) are a.s. strictly increasing. Then the sample paths t (t) are a.s. continuous and nondecreasing, with P () > = ; e.g., see Proposition A. in Veillette and Taqqu (2). (2.4)
4 86 M.M. Meerschaert et al. / Statistics and Probability Letters 83 (23) Uncoupled CTRW Consider a CTRW whose i.i.d. waiting times {W n, n } belong to the domain of attraction of a positive β-stable random variable D(), and whose i.i.d. jumps {J n, n } belong to the strict domain of attraction of a d-dimensional stable random vector Y(). We assume that {W n } and {J n } are independent; that is, the CTRW is uncoupled. It follows from Theorem 4.2 in Meerschaert and Scheffler (24) that the scaling limit of this CTRW is a time-changed process X(t) = Y((t)), where (t) is the inverse (.) of a β-stable subordinator D. Since D is self-similar with index /β, its inverse is self-similar with index β. Since Y is independent of, the CTRW scaling limit X is self-similar with index β/α. Proposition 2.3. The uncoupled CTRW limit has a.s. dim H X([, ]) = dim P X([, ]) = min{d, α} (2.5) and max{, α} if α d, dim H GrX([, ]) = dim P GrX([, ]) = + β if α > d =. α (2.6) Proof. The result (2.5) follows from Proposition 2. and the results of Blumenthal and Getoor (96a,b) on the Hausdorff dimension and Pruitt and Taylor (996) on the packing dimension of the range of the stable Lévy process Y. To prove (2.6), recall from Pruitt and Taylor (969) that for any constant a >, β if α β, α dim H Z([, a]) = if β < α d, + β a.s. (2.7) if α > d =, α Theorem 3.2 in Meerschaert and Xiao (25) shows that dim P Z([, a]) also equals the right hand side of (2.7). Then (2.6) follows using Theorem 2.2. Remark 2.4. In increasing generality, Blumenthal and Getoor (962), Jain and Pruitt (968), Pruitt and Taylor (969) and Rezakhanlou and Taylor (988) showed that max{, α} if α d dim H GrY([, ]) = dim P GrY([, ]) = 2 a.s. if α > d = α Compare with (2.6) to see that the inverse stable time change modifies the fractal dimension of the CTRW limit graph in dimension d = when α > Coupled CTRW In a coupled CTRW, the space time jumps {(J n, W n ), n } are i.i.d., but J n can depend on the waiting time W n. Now the CTRW S(N(t)) has scaling limit Y((t )) and the so-called oracle CTRW S(N(t)+) has scaling limit Y((t)); see Henry and Straka (2) or Jurlewicz et al. (22). In the uncoupled case, the two limit processes are the same. The proof of Theorem 2.2 extends immediately to the process Y((t )), with the same dimension results, because the graphs of Z(u) = (D(u), Y(u)) and Z (u) = (D(u), Y(u )) have the same Hausdorff and packing dimensions, as they differ by at most a countable number of discrete points. In the following, we discuss examples for Y((t)), with the understanding that the same dimension results hold for Y((t )). The simplest case is W n = J n, so that X(t) = D((t)). This process is self-similar with index ; see for example Becker- Kern et al. (24). It follows from Proposition 2., Theorem 2.2 and the fact that for any constant a >, that and dim H D([, a]) = dim H {(D(x), D(x)) : x [, a]} = β, a.s. dim H X([, ]) = dim P X([, ]) = β, a.s. (2.8) dim H GrX([, ]) = dim P GrX([, ]) =, a.s. (2.9) These results can also be obtained using uniform Hausdorff and packing dimension results for the β-stable subordinator; see Perkins and Taylor (987).
5 M.M. Meerschaert et al. / Statistics and Probability Letters 83 (23) Shlesinger et al. (982) consider a CTRW where the waiting times W n are i.i.d. with a β-stable random variable D such that (e sd ) = e sβ and, conditional on W n = t, the jump J n is normal with mean zero and variance 2t. Then J n is symmetric stable with index α = 2β. This model was applied to stock market prices by Meerschaert and Scalas (26). Becker-Kern et al. (24) show that the CTRW limit is X(t) = Y((t)), where Y is a real-valued stable Lévy process with index α = 2β and (t) is the inverse of a β-stable subordinator, which is not independent of Y. Here X(t) is self-similar with index /2, the same as Brownian motion. However, the Hausdorff dimensions of the range and graph of X are completely different than those for Brownian motion. Proposition 2. gives that dim H X([, ]) = min{, 2β} a.s. To determine the Hausdorff dimension of the graph of X(t), we first verify that the Fourier Laplace transform of (D(), Y()) is e iξ Y() ηd() = e ηd() e iξ Y() D() = e (η+ξ 2 )D() = e (η+ξ 2 ) β. It follows that the Lévy process Z(u) = (D(u), Y(u)) is operator stable (cf. Meerschaert and Scheffler, 2) with the unique exponent β C =. (2β) (2.2) Now Theorem 3.2 in Meerschaert and Xiao (25) implies that for any a >, 2β if 2β, dim H Z([, a]) = dim P Z([, a]) = 2 + β if 2β >, a.s. (2.2) Consequently, we use Theorem 2.2 to derive dim H GrX([, ]) = dim P GrX([, ]) = max, β +, a.s. (2.22) CTRW triangular arrays Meerschaert and Scheffler (28) prove limit theorems for CTRW triangular arrays. Here Z(u) = (Y(u), D(u)) is a Lévy process with values in R p and characteristic exponent Φ (i.e., (e i ξ,z(u) ) = e uφ(ξ) ). Proposition 2.5. Let X = {Y((t)), t }, where Y is a Lévy process with values in R d and characteristic exponent ψ, and (t) is the inverse (.) of a subordinator {D(u), u } with characteristic exponent σ. Suppose Z = (D(u), Y(u)), u is a Lévy process on R +d whose characteristic exponent Φ satisfies K Re Re K Re (2.23) + σ (η) + ψ(ξ) + Φ(η, ξ) + σ (η) + ψ(ξ) for all (η, ξ) R +d with η + ξ large, where K is a constant. Then a.s. dξ dim H X([, ]) = sup γ < d : + ψ(ξ) ξ < γ Re {ξ R d : ξ } and dim H GrX([, ]) = max{, χ} a.s., where χ = sup γ < + d : + σ (η) + ψ(ξ) Re { η + ξ } dηdξ ( η + ξ ) γ <. Proof. Corollary.8 in Khoshnevisan et al. (23) shows that for any a >, dξ dim H Z([, a]) = sup γ < p : Re {ξ R p : ξ } + Φ(ξ) ξ <, γ a.s. (2.24) where is the uclidean norm on R p. Apply Proposition 2. and Theorem 2.2.
6 88 M.M. Meerschaert et al. / Statistics and Probability Letters 83 (23) Proposition 2.6. Let X, Z be as in Proposition 2.5. If the characteristic exponent Φ of Z satisfies (2.23), then W(r) dim P X([, ]) = sup γ : lim inf = r + r γ, a.s., where W(r) is defined by d W(r) = Re R d + ψ(ξ/r) j= + ξ 2 j dξ, and dim P GrX([, ]) = max{, χ } almost surely, where W(r) χ = sup γ : lim inf = r + r γ and where W(r) = Re R +d + σ (η/r) + ψ(ξ/r) + η 2 d j= + ξ 2 j dη dξ. (2.25) Proof. Theorem. in Khoshnevisan and Xiao (28) shows that for any a >, W(r) dim P Z([, a]) = sup γ : lim inf = r + r γ, a.s., with W(r) given by (2.25). Apply Proposition 2. and Theorem 2.2. Our first example is a stochastic model for ultraslow diffusion from Meerschaert and Scheffler (26); see also Chechkin et al. (22). At each scale c > we are given i.i.d. waiting times (W c) n and i.i.d. jumps (J c n ). Let Sc (n) = J c + + J c n and T c (n) = W c+ +W c n, and suppose Sc (ct) Y(t) and T c (ct) D(t) as c, where Y(t) and D(t) are independent Lévy processes. Letting N c t = max{n : T c (n) t}, the CTRW scaling limit S c (N c) t Y( t) by Theorem 2. in Meerschaert and Scheffler (28). To be specific, take {B i } i.i.d. with < B i < and assume P{W c i > u B i = β} = c u β for u c /β, so that the waiting times are Pareto distributed, conditional on the mixing variables. Then [e sd(t) ] = e tσ (s), where (2.4) holds with ν(t, ) = t β µ(dβ), and µ is the distribution of the mixing variable. Suppose < β < β 2 < < β n < and take µ(dβ) = n k= d β k (Γ k ( β k)) δ βk (dβ), where δ a is the unit mass at a. Then D(u) = n k= d kd k (u), a mixture of independent β k -stable subordinators. Lemma 2.7. Let Y = {Y(x), x } be a strictly stable Lévy motion of index α (, 2] with values in R d, independent of D, and let Φ be the characteristic exponent of the Lévy process Z(u) = (D(u), Y(u)). Then for all (η, ξ) R +d that satisfies η + ξ >, we have K η β n + ξ α Re K + Φ(η, ξ) η, β n + ξ α (2.27) where K is a constant that may depend on n, α, β k, d k. (2.26) Proof. For simplicity, assume that Y has characteristic exponent ψ(ξ) = ξ α. Then Φ(η, ξ) = = n idk η β k + ξ α k= n d k η β k [cos(πβ k /2) i sin(πβ k /2)] + ξ α k= =: f (η, ξ) ig(η). (2.28)
7 M.M. Meerschaert et al. / Statistics and Probability Letters 83 (23) Since β k (, ), we have f (η, ξ) for all η, ξ R +d. Moreover, g(η) Kf (η, ξ) for some constant K >. Hence ( + K 2 )( + f (η, ξ)) Re + Φ(η, ξ) + f (η, ξ). From here it is elementary to verify (2.27). Proposition 2.8. For the triangular array CTRW limit described above, we have a.s. dim H X([, ]) = dim P X([, ]) = min{d, α}, (2.29) and max{, α} if α d, dim H GrX([, ]) = dim P GrX([, ]) = + β n if α > d =. α (2.3) Proof. The proof is similar to Proposition 4. in Meerschaert and Xiao (25), using Lemma 2.7, and Propositions 2.5 and 2.6, hence we omit the details CTRW with correlated jumps Now we consider an uncoupled CTRW whose jumps {J n } form a correlated sequence of random variables, and whose waiting times {W n } are i.i.d. and belong to the domain of attraction of a positive β-stable random variable D(). In this case, Meerschaert et al. (29) show that, under certain conditions on the correlation structure of jumps, the CTRW scaling limit is the (Hβ)-self-similar process X = {Y((t)) : t }, where Y is a fractional Brownian motion with index H (, ), and (t) is the inverse of a β-stable subordinator D, independent of Y. Proposition 2.9. The correlated CTRW limit X described above satisfies a.s. dim H X([, ]) = dim P X([, ]) = min, H, (2.3) and dim H GrX([, ]) = dim P GrX([, ]) = β + ( Hβ). (2.32) The proof of Proposition 2.9 requires a few preliminary results. Let c 7 > be a fixed constant. A collection Λ(b) of intervals of length b in R is called c 7 -nested if no interval of length b in R can intersect more than c 7 intervals of Λ(b). Note that for each integer n, the collection of dyadic intervals I n,j = [j/2 n, (j + )/2 n ] is c 7 -nested with c 7 = 3. Lemma 2.. Let {D(t), t } be a β-stable subordinator and let Λ(b) be a c 7 -nested family. Denote by M u (b, s) the number of intervals in Λ(b) which intersect D([u, u + s]). Then there exists a positive constant c 8 such that for all u and all < b β s, M u (b, s) c 8 sb β. (2.33) If one takes b = s, then we have M u (b, s) c 8. (2.34) Proof. This is an immediate consequence of Lemma 6. in Pruitt and Taylor (969). It can also be derived from Lemma 3.2 in Liu and Xiao (998) where general self-similar Markov processes are considered. Lemma 2.. Let the assumptions of Proposition 2.9 hold and let a > be a constant. Then /H dim P Z([, a]) β + ( Hβ)d if Hd, if > Hd, a.s. (2.35) Proof. The proof is based on a moment argument. Note that for every ε > the function Y(u) ( u a) satisfies a uniform Hölder condition of order H ε. We divide the interval [, a] into ( a + )2 n dyadic intervals I n,j of length 2 n.
8 9 M.M. Meerschaert et al. / Statistics and Probability Letters 83 (23) First we construct a covering of the range Z([, a]) by using balls in R d+ of radius 2 Hn as follows. Define t n,j = j/2 n so that for each I n,j = [t n,j, t n,j + 2 n ], the image Y(I n,j ) is contained in a ball in R d of radius sup s In,j Y(s) Y(t n,j ) and can be covered by at most d sup Y(s) Y(t n,j ) s I n,j N n,j = c 9 (2.36) 2 Hn balls of radius 2 Hn. By the self-similarity and stationarity of increments of Y, we have d N n,j = c9 2 Hdn = c 9 sup Y(s) Y(t n,j ) s I n,j sup Y(s) s [,] d := c <, (2.37) where the last inequality follows from the well known tail probability for the supremum of Gaussian processes (e.g., Fernique s inequality). To get a covering for D(I n,i ), let Γ (2 n ) be the collection of dyadic intervals of order n in R +. Let M n,j be the number of dyadic intervals in Γ (2 n ) which intersect D(I n,j ). Applying (2.34) in Lemma 2. with b n = s n = 2 n, we obtain that (M n,j ) c 8, j ( a + )2 n. (2.38) Since 2 n < 2 Hn, we see that Z(I n,j ) = {(D(x), Y(x)) : x I n,j } can be covered by at most M n,j N n,j balls in R d+ of radius 2 Hn. Denote by N Z([, a]), 2 Hn the smallest number of balls in R d+ of radius 2 Hn that cover Z([, a]), then N Z([, a]), 2 Hn ( a +)2 n j= M n,j N n,j. It follows from (2.37), (2.38) and the independence of Y and D that N Z([, a]), 2 Hn ( a + )c 8 c 2 n. Hence, for any ε >, P N Z([, a]), 2 Hn ( a + )c 8 c 2 n(+ε) 2 nε. It follows from the Borel Cantelli lemma that almost surely N Z([, a]), 2 Hn < ( a + )c 8 c 2 n(+ε) for all n large enough. The upper box-counting dimension of F is defined as dim M F = lim sup ε log N(F, ε). log ε (2.39) Hence dim M Z([, a]) ( + ε)/h a.s. It is well known that for every (bounded) set F R d, dim H F dim P F dim M F d. (2.4) Since ε > is arbitrary, it follows that dim P Z([, a]) /H almost surely. Next we construct a covering for the range Z([, a]) by using balls in R d+ of radius 2 n/β. Let Γ (2 n/β ) be the collection of intervals in R + of the form I = [ k n,k, k+ 2 n/β 2 n/β ], where k is an integer. Then the class Γ (2 n/β ) is 3-nested. Let M n,j be the number of intervals in Γ (2 n/β ) that intersect D(I n,j ). By Lemma 2. with b n = 2 n/β and s n = 2 n, we derive (M ) n,j c 8. Thus D(I n,j ) can almost surely be covered by M n,j intervals of length 2 n/β from Γ (2 n/β ). On the other hand, the image Y(I n,j ) can be covered by at most d sup Y(s) Y(t n,j ) N = n,j c s I n,j 9 2 n/β
9 M.M. Meerschaert et al. / Statistics and Probability Letters 83 (23) balls of radius 2 n/β, where t n,j = j/2 n, and then similar to (2.37) we derive N n,j = c 2 n β H d. (2.4) Denote by N Z([, a]), 2 n/β the smallest number of balls in R d+ of radius 2 n/β that cover Z([, a]), then N Z([, a]), 2 n/β ( a +)2 n j= M n,j N n,j. By (2.4) and the independence of Y and D we have N Z([, a]), 2 n/β ( a + )c 8 c 2 n + β H d Hence, for any ε >, the Borel Cantelli Lemma implies that a.s. N Z([, a]), 2 n/β + < ( a + )c 8 c 2 n β H d+ε. for all n large enough. This and (2.39) imply that dim M Z([, a]) β + ( βh)d + βε almost surely which, in turn, implies dim P Z([, a]) β + ( βh)d a.s. Combining the above we have dim P Z([, a]) min, β + ( βh)d a.s. H This proves (2.35). Lemma 2.2. Under the assumptions of Proposition 2.9, we have /H dim H Z([, ]) β + ( Hβ)d if Hd, if > Hd, a.s. (2.42) Proof. The projection of Z([, ]) into R d is Y([, ]), and dim H Y([, ]) = a.s. when Hd. This implies the first H inequality in (2.42). To prove the inequality in (2.42) for the case > Hd, by Frostman s theorem (cf. Kahane, 985, p. 33) along with the inequality Z(x) Z(y) 2 D(x) D(y) + Y(x) Y(y), it is sufficient to prove that for every constant γ (, β + ( Hβ)d), we have a a dx dy D(x) D(y) + Y(x) Y(y) γ <. (2.43) Since > Hd, we have β + ( Hβ)d > d. We only need to verify (2.43) for every γ (d, β + ( Hβ)d). For this purpose, we will make use of the following easily verifiable fact (see, e.g., Kahane, 985, p. 279): if Ξ is a standard normal vector in R d, then there is a finite constant c > such that for any constants γ > d and ρ, ρ + Ξ γ c ρ (γ d). Fix x, y [, a] such that x y. We use to denote the conditional expectation given the subordinator D, apply the above fact with ρ = D(x) D(y) x y H and use the self-similarity of D to derive γ = x y Hγ D(x) D(y) + Y(x) Y(y) (ρ + Ξ ) γ c x y Hγ H(γ d) x y D(x) D(y) γ d = c 2, (2.44) x y Hd+(γ d)/β
10 92 M.M. Meerschaert et al. / Statistics and Probability Letters 83 (23) where the last equality follows from the /β-self-similarity of D and the constant c 2 Hawkes (97, Lemma ) that, as r +, P(D() r) c 3 r β/(2( β)) exp ( β)β β/( β) r β/( β), = c D() (γ d). Recall from where c 3 = 2π( β)β β/(2( β)) /2. We verify easily c2 <. It follows from Fubini s theorem and (2.44) that a a dx dy D(x) D(y) + Y(x) Y(y) γ c 2 a a dx dy <, (2.45) x y Hd+(γ d)/β the last integral is convergent because Hd + (γ d)/β <. This proves (2.43) and thus the lemma. Proof of Proposition 2.9. q. (2.3) follows from Proposition 2. and the well-known results on the Hausdorff and packing dimensions for the range of a fractional Brownian motion (see, e.g., Chapter 8 of Kahane, 985). In order to prove (2.32), apply Theorem 2.2 and (2.4) along with Lemmas 2. and 2.2. Remark 2.3. For a CTRW with dependent heavy tailed jumps, the outer process Y can be a linear fractional stable motion; see Meerschaert et al. (29). Proposition 2. and Theorem 2.2 are applicable here, but the Hausdorff dimension of the range and graph sets of the processes Y and Z(u) = (D(u), Y(u)) are unknown in general. For some partial results in this direction, see Shieh and Xiao (2) and Xiao and Lin (994). 3. Discussion CTRW limits are random fractals. The fractal dimension is altered by a random time change that represents particle resting times between diffusive movements. An important and useful example is the inverse stable subordinator with index < β <, that corresponds to power law waiting times with an infinite mean. The graph of a time-changed Brownian motion goes from dimension 3/2 to dimension +β/2. The fractal dimension of a fractional Brownian motion graph changes from 2 H to β + ( Hβ). The graph of a stable Lévy motion in one dimension with index < α < 2 is a random fractal with index 2 α. After the time change, the graph has dimension + β( α ). In every case, the inverse β-stable subordinator alters the fractal dimension, such that substituting β = recovers the dimension formula for the original outer process. Acknowledgments The research of M.M. Meerschaert was partially supported by NSF grant DMS and NIH grant R-B279. The research of Y. Xiao was partially supported by NSF grant DMS-693. References Barlow, M.T., Cerný, J., 2. Convergence to fractional kinetics for random walks associated with unbounded conductances. Probab. Theory Related Fields 49, Becker-Kern, P., Meerschaert, M.M., Scheffler, H.P., 24. Limit theorems for coupled continuous time random walks. Ann. Probab. 32, Blumenthal, R.M., Getoor, R., 96a. Some theorems on stable processes. Trans. Amer. Math. Soc. 95, Blumenthal, R.M., Getoor, R., 96b. A dimension theorem for sample functions of stable processes. Illinois J. Math. 4, Blumenthal, R.M., Getoor, R., 962. The dimension of the set of zeros and the graph of a symmetric stable process. Illinois J. Math. 6, Chechkin, A.V., Gorenflo, R., Sokolov, I.M., 22. Retarding subdiffusion and accelerating superdiffusion governed by distributed-order fractional diffusion equations. Phys. Rev. 66, Davydov, Yu., 22. On convex hull of d-dimensional fractional Brownian motion. Statist. Probab. Lett. 82, Falconer, K.J., 99. Fractal Geometry Mathematical Foundations and Applications. Wiley, New York. Hawkes, J., 97. A lower Lipschitz condition for stable subordinator. Z. Wahrscheinlichkeitstheor. Verwandte Geb. 7, Henry, B.I., Straka, P., 2. Lagging and leading coupled continuous time random walks, renewal times and their joint limits. Stochastic Process. Appl. 2, Jain, N.C., Pruitt, W.., 968. The correct measure function for the graph of a transient stable process. Z. Wahrscheinlichkeitstheor. Verwandte Geb. 9, Jurlewicz, A., Kern, P., Meerschaert, M.M., Scheffler, H.-P., 22. Fractional governing equations for coupled random walks. Comput. Math. Appl. 64, Kahane, J.-P., 985. Some Random Series of Functions, second ed. Cambridge University Press, Cambridge. Khoshnevisan, D., Xiao, Y., 28. Packing dimension of the range of a Lévy process. Proc. Amer. Math. Soc. 36, Khoshnevisan, D., Xiao, Y., Zhong, Y., 23. Measuring the range of an additive Lévy processes. Ann. Probab. 3, Liu, L., Xiao, Y., 998. Hausdorff dimension theorems for self-similar Markov processes. Probab. Math. Statist. 8, Meerschaert, M.M., Nane,., Xiao, Y., 29. Correlated continuous time random walks. Statist. Probab. Lett. 79, Meerschaert, M.M., Scalas,., 26. Coupled continuous time random walks in finance. Physica A 37, 4 8. Meerschaert, M.M., Scheffler, H.P., 2. Limit Distributions for Sums of Independent Random Vectors: Heavy Tails in Theory and Practice. Wiley Interscience, New York. Meerschaert, M.M., Scheffler, H.P., 24. Limit theorems for continuous time random walks with infinite mean waiting times. J. Appl. Probab. 4, Meerschaert, M.M., Scheffler, H.P., 26. Stochastic model for ultraslow diffusion. Stochastic Process. Appl. 6,
11 M.M. Meerschaert et al. / Statistics and Probability Letters 83 (23) Meerschaert, M.M., Scheffler, H.P., 28. Triangular array limits for continuous time random walks. Stochastic Process. Appl. 8, Meerschaert, M.M., Sikorskii, A., 22. Stochastic Models for Fractional Calculus. De Gruyter, Berlin. Meerschaert, M.M., Xiao, Y., 25. Dimension results for the sample paths of operator stable processes. Stochastic Process. Appl. 5, Perkins,.A., Taylor, S.J., 987. Uniform measure results for the image of subsets under Brownian motion. Probab. Theory Related Fields 76, Pruitt, W.., Taylor, S.J., 969. Sample path properties of processes with stable components. Z. Wahrscheinlichkeitstheor. Verwandte Geb. 2, Pruitt, W.., Taylor, S.J., 996. Packing and covering indices for a general Lévy process. Ann. Probab. 24, Rezakhanlou, F., Taylor, S.J., 988. The packing measure of the graph of a stable process. Astérisque (57 58), Shieh, N.-R., Xiao, Y., 2. Hausdorff and packing dimensions of the images of random fields. Bernoulli 6, Shlesinger, M., Klafter, J., Wong, Y.M., 982. Random walks with infinite spatial and temporal moments. J. Stat. Phys. 27, Taylor, S.J., 986. The measure theory of random fractals. Math. Proc. Cambridge Philos. Soc., Veillette, M., Taqqu, M.S., 2. Using differential equations to obtain joint moments of first-passage times of increasing Lévy processes. Statist. Probab. Lett. 8, Xiao, Y., 24. Random fractals and Markov processes. In: Lapidus, M.L., van Frankenhuijsen, M. (ds.), Fractal Geometry and Applications: A Jubilee of Benoit Mandelbrot. American Mathematical Society, pp Xiao, Y., Lin, H., 994. Dimension properties of the sample paths of self-similar processes. Acta Math. Sinica (NS),
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 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 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 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 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 informationA Uniform Dimension Result for Two-Dimensional Fractional Multiplicative Processes
To appear in Ann. Inst. Henri Poincaré Probab. Stat. A Uniform Dimension esult for Two-Dimensional Fractional Multiplicative Processes Xiong Jin First version: 8 October 213 esearch eport No. 8, 213, Probability
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 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 informationRandom Fractals and Markov Processes
Proceedings of Symposia in Pure Mathematics Random Fractals and Markov Processes Yimin Xiao Abstract. This is a survey on the sample path properties of Markov processes, especially fractal properties of
More informationAsymptotic Behavior of Semistable Lévy Exponents and Applications to Fractal Path Properties
DOI.7/s959-6-72-6 Asymptotic Behavior of Semistable Lévy Exponents and Applications to Fractal Path Properties Peter Kern Mark M. Meerschaert 2 Yimin Xiao 2 Received: 27 June 26 Springer Science+Business
More informationAn almost sure invariance principle for additive functionals of Markov chains
Statistics and Probability Letters 78 2008 854 860 www.elsevier.com/locate/stapro An almost sure invariance principle for additive functionals of Markov chains F. Rassoul-Agha a, T. Seppäläinen b, a Department
More informationAdditive Lévy Processes
Additive Lévy Processes Introduction Let X X N denote N independent Lévy processes on. We can construct an N-parameter stochastic process, indexed by N +, as follows: := X + + X N N for every := ( N )
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 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 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 information9 Brownian Motion: Construction
9 Brownian Motion: Construction 9.1 Definition and Heuristics The central limit theorem states that the standard Gaussian distribution arises as the weak limit of the rescaled partial sums S n / p n of
More informationSTAT 7032 Probability Spring Wlodek Bryc
STAT 7032 Probability Spring 2018 Wlodek Bryc Created: Friday, Jan 2, 2014 Revised for Spring 2018 Printed: January 9, 2018 File: Grad-Prob-2018.TEX Department of Mathematical Sciences, University of Cincinnati,
More informationWiener Measure and Brownian Motion
Chapter 16 Wiener Measure and Brownian Motion Diffusion of particles is a product of their apparently random motion. The density u(t, x) of diffusing particles satisfies the diffusion equation (16.1) u
More informationGaussian Random Fields: Geometric Properties and Extremes
Gaussian Random Fields: Geometric Properties and Extremes Yimin Xiao Michigan State University Outline Lecture 1: Gaussian random fields and their regularity Lecture 2: Hausdorff dimension results and
More informationFRACTAL BEHAVIOR OF MULTIVARIATE OPERATOR-SELF-SIMILAR STABLE RANDOM FIELDS
Communications on Stochastic Analysis Vol. 11, No. 2 2017) 233-244 Serials Publications www.serialspublications.com FRACTAL BEHAVIOR OF MULTIVARIATE OPERATOR-SELF-SIMILAR STABLE RANDOM FIELDS ERCAN SÖNMEZ*
More informationLIMIT THEOREMS FOR SHIFT SELFSIMILAR ADDITIVE RANDOM SEQUENCES
Watanabe, T. Osaka J. Math. 39 22, 561 63 LIMIT THEOEMS FO SHIFT SELFSIMILA ADDITIVE ANDOM SEQUENCES TOSHIO WATANABE eceived November 15, 2 1. Introduction We introduce shift selfsimilar random sequences,
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
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 informationPacking Dimension Profiles and Lévy Processes
Packing Dimension Profiles and Lévy Processes D. Khoshnevisan, R.L. Schilling and Y. Xiao Abstract We extend the concept of packing dimension profiles, due to Falconer and Howroyd (997 and Howroyd (200,
More information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
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 informationMAJORIZING MEASURES WITHOUT MEASURES. By Michel Talagrand URA 754 AU CNRS
The Annals of Probability 2001, Vol. 29, No. 1, 411 417 MAJORIZING MEASURES WITHOUT MEASURES By Michel Talagrand URA 754 AU CNRS We give a reformulation of majorizing measures that does not involve measures,
More informationBrownian Motion. Chapter Stochastic Process
Chapter 1 Brownian Motion 1.1 Stochastic Process A stochastic process can be thought of in one of many equivalent ways. We can begin with an underlying probability space (Ω, Σ,P and a real valued stochastic
More informationGAUSSIAN PROCESSES; KOLMOGOROV-CHENTSOV THEOREM
GAUSSIAN PROCESSES; KOLMOGOROV-CHENTSOV THEOREM STEVEN P. LALLEY 1. GAUSSIAN PROCESSES: DEFINITIONS AND EXAMPLES Definition 1.1. A standard (one-dimensional) Wiener process (also called Brownian motion)
More informationVolatility and Returns in Korean Futures Exchange Markets
Volatility and Returns in Korean Futures Exchange Markets Kyungsik Kim*,, Seong-Min Yoon and Jum Soo Choi Department of Physics, Pukyong National University, Pusan 608-737, Korea Division of Economics,
More informationHölder regularity for operator scaling stable random fields
Stochastic Processes and their Applications 119 2009 2222 2248 www.elsevier.com/locate/spa Hölder regularity for operator scaling stable random fields Hermine Biermé a, Céline Lacaux b, a MAP5 Université
More informationCorrelation dimension for self-similar Cantor sets with overlaps
F U N D A M E N T A MATHEMATICAE 155 (1998) Correlation dimension for self-similar Cantor sets with overlaps by Károly S i m o n (Miskolc) and Boris S o l o m y a k (Seattle, Wash.) Abstract. We consider
More informationSupermodular ordering of Poisson arrays
Supermodular ordering of Poisson arrays Bünyamin Kızıldemir Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological University 637371 Singapore
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 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 informationPrime numbers and Gaussian random walks
Prime numbers and Gaussian random walks K. Bruce Erickson Department of Mathematics University of Washington Seattle, WA 9895-4350 March 24, 205 Introduction Consider a symmetric aperiodic random walk
More informationINDISTINGUISHABILITY OF ABSOLUTELY CONTINUOUS AND SINGULAR DISTRIBUTIONS
INDISTINGUISHABILITY OF ABSOLUTELY CONTINUOUS AND SINGULAR DISTRIBUTIONS STEVEN P. LALLEY AND ANDREW NOBEL Abstract. It is shown that there are no consistent decision rules for the hypothesis testing problem
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 informationAdditive functionals of infinite-variance moving averages. Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535
Additive functionals of infinite-variance moving averages Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535 Departments of Statistics The University of Chicago Chicago, Illinois 60637 June
More informationPhysica A. Extremal behavior of a coupled continuous time random walk
Physica A 39 (211) 55 511 Contents lists available at ScienceDirect Physica A journal homepage: www.elsevier.com/locate/physa Extremal behavior of a coupled continuous time random walk Rina Schumer a,,
More informationLecture 25: Large Steps and Long Waiting Times
Lecture 25: Large Steps and Long Waiting Times Scribe: Geraint Jones (and Martin Z. Bazant) Department of Economics, MIT Proofreader: Sahand Jamal Rahi Department of Physics, MIT scribed: May 10, 2005,
More informationAsymptotic stability of an evolutionary nonlinear Boltzmann-type equation
Acta Polytechnica Hungarica Vol. 14, No. 5, 217 Asymptotic stability of an evolutionary nonlinear Boltzmann-type equation Roksana Brodnicka, Henryk Gacki Institute of Mathematics, University of Silesia
More information(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS
(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS Svetlana Janković and Miljana Jovanović Faculty of Science, Department of Mathematics, University
More informationMi-Hwa Ko. t=1 Z t is true. j=0
Commun. Korean Math. Soc. 21 (2006), No. 4, pp. 779 786 FUNCTIONAL CENTRAL LIMIT THEOREMS FOR MULTIVARIATE LINEAR PROCESSES GENERATED BY DEPENDENT RANDOM VECTORS Mi-Hwa Ko Abstract. Let X t be an m-dimensional
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 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 informationMath212a1413 The Lebesgue integral.
Math212a1413 The Lebesgue integral. October 28, 2014 Simple functions. In what follows, (X, F, m) is a space with a σ-field of sets, and m a measure on F. The purpose of today s lecture is to develop the
More informationFrom Fractional Brownian Motion to Multifractional Brownian Motion
From Fractional Brownian Motion to Multifractional Brownian Motion Antoine Ayache USTL (Lille) Antoine.Ayache@math.univ-lille1.fr Cassino December 2010 A.Ayache (USTL) From FBM to MBM Cassino December
More informationThe Brownian graph is not round
The Brownian graph is not round Tuomas Sahlsten The Open University, Milton Keynes, 16.4.2013 joint work with Jonathan Fraser and Tuomas Orponen Fourier analysis and Hausdorff dimension Fourier analysis
More informationProbability and Measure
Part II Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 84 Paper 4, Section II 26J Let (X, A) be a measurable space. Let T : X X be a measurable map, and µ a probability
More information2 (Bonus). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due 9/5). Prove that every countable set A is measurable and µ(a) = 0. 2 (Bonus). Let A consist of points (x, y) such that either x or y is
More informationMATHS 730 FC Lecture Notes March 5, Introduction
1 INTRODUCTION MATHS 730 FC Lecture Notes March 5, 2014 1 Introduction Definition. If A, B are sets and there exists a bijection A B, they have the same cardinality, which we write as A, #A. If there exists
More informationWeak Variation of Gaussian Processes
Journal of Theoretical Probability, Vol. 10, No. 4, 1997 Weak Variation of Gaussian Processes Yimin Xiao1,2 Received December 6, 1995; revised June 24, 1996 Let X(l) (ter) be a real-valued centered Gaussian
More informationConvergence of a Generalized Midpoint Iteration
J. Able, D. Bradley, A.S. Moon under the supervision of Dr. Xingping Sun REU Final Presentation July 31st, 2014 Preliminary Words O Rourke s conjecture We begin with a motivating question concerning the
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 informationBrownian Motion and Conditional Probability
Math 561: Theory of Probability (Spring 2018) Week 10 Brownian Motion and Conditional Probability 10.1 Standard Brownian Motion (SBM) Brownian motion is a stochastic process with both practical and theoretical
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 informationHAUSDORFF DIMENSION AND ITS APPLICATIONS
HAUSDORFF DIMENSION AND ITS APPLICATIONS JAY SHAH Abstract. The theory of Hausdorff dimension provides a general notion of the size of a set in a metric space. We define Hausdorff measure and dimension,
More information1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer 11(2) (1989),
Real Analysis 2, Math 651, Spring 2005 April 26, 2005 1 Real Analysis 2, Math 651, Spring 2005 Krzysztof Chris Ciesielski 1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer
More informationLeast Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises
Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,
More informationMetric Spaces and Topology
Chapter 2 Metric Spaces and Topology From an engineering perspective, the most important way to construct a topology on a set is to define the topology in terms of a metric on the set. This approach underlies
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 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 informationSelf-similar fractals as boundaries of networks
Self-similar fractals as boundaries of networks Erin P. J. Pearse ep@ou.edu University of Oklahoma 4th Cornell Conference on Analysis, Probability, and Mathematical Physics on Fractals AMS Eastern Sectional
More informationThe strictly 1/2-stable example
The strictly 1/2-stable example 1 Direct approach: building a Lévy pure jump process on R Bert Fristedt provided key mathematical facts for this example. A pure jump Lévy process X is a Lévy process such
More informationCoupled continuous time random walks in finance
Physica A 37 (26) 114 118 www.elsevier.com/locate/physa Coupled continuous time random walks in finance Mark M. Meerschaert a,,1, Enrico Scalas b,2 a Department of Mathematics and Statistics, University
More informationERRATA: Probabilistic Techniques in Analysis
ERRATA: Probabilistic Techniques in Analysis ERRATA 1 Updated April 25, 26 Page 3, line 13. A 1,..., A n are independent if P(A i1 A ij ) = P(A 1 ) P(A ij ) for every subset {i 1,..., i j } of {1,...,
More informationA NOTE ON THE COMPLETE MOMENT CONVERGENCE FOR ARRAYS OF B-VALUED RANDOM VARIABLES
Bull. Korean Math. Soc. 52 (205), No. 3, pp. 825 836 http://dx.doi.org/0.434/bkms.205.52.3.825 A NOTE ON THE COMPLETE MOMENT CONVERGENCE FOR ARRAYS OF B-VALUED RANDOM VARIABLES Yongfeng Wu and Mingzhu
More informationSYMMETRIC STABLE PROCESSES IN PARABOLA SHAPED REGIONS
SYMMETRIC STABLE PROCESSES IN PARABOLA SHAPED REGIONS RODRIGO BAÑUELOS AND KRZYSZTOF BOGDAN Abstract We identify the critical exponent of integrability of the first exit time of rotation invariant stable
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 informationRecent Developments on Fractal Properties of Gaussian Random Fields
Recent Developments on Fractal Properties of Gaussian Random Fields Yimin Xiao Abstract We review some recent developments in studying fractal and analytic properties of Gaussian random fields. It is shown
More informationDo stochastic processes exist?
Project 1 Do stochastic processes exist? If you wish to take this course for credit, you should keep a notebook that contains detailed proofs of the results sketched in my handouts. You may consult any
More informationSelf-similar fractals as boundaries of networks
Self-similar fractals as boundaries of networks Erin P. J. Pearse ep@ou.edu University of Oklahoma 4th Cornell Conference on Analysis, Probability, and Mathematical Physics on Fractals AMS Eastern Sectional
More informationThe Chaotic Character of the Stochastic Heat Equation
The Chaotic Character of the Stochastic Heat Equation March 11, 2011 Intermittency The Stochastic Heat Equation Blowup of the solution Intermittency-Example ξ j, j = 1, 2,, 10 i.i.d. random variables Taking
More informationKrzysztof Burdzy University of Washington. = X(Y (t)), t 0}
VARIATION OF ITERATED BROWNIAN MOTION Krzysztof Burdzy University of Washington 1. Introduction and main results. Suppose that X 1, X 2 and Y are independent standard Brownian motions starting from 0 and
More informationMultiple points of the Brownian sheet in critical dimensions
Multiple points of the Brownian sheet in critical dimensions Robert C. Dalang Ecole Polytechnique Fédérale de Lausanne Based on joint work with: Carl Mueller Multiple points of the Brownian sheet in critical
More informationProblem set 1, Real Analysis I, Spring, 2015.
Problem set 1, Real Analysis I, Spring, 015. (1) Let f n : D R be a sequence of functions with domain D R n. Recall that f n f uniformly if and only if for all ɛ > 0, there is an N = N(ɛ) so that if n
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 informationExistence and Uniqueness
Chapter 3 Existence and Uniqueness An intellect which at a certain moment would know all forces that set nature in motion, and all positions of all items of which nature is composed, if this intellect
More informationUniformly and strongly consistent estimation for the Hurst function of a Linear Multifractional Stable Motion
Uniformly and strongly consistent estimation for the Hurst function of a Linear Multifractional Stable Motion Antoine Ayache and Julien Hamonier Université Lille 1 - Laboratoire Paul Painlevé and Université
More informationNECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES
NECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES JUHA LEHRBÄCK Abstract. We establish necessary conditions for domains Ω R n which admit the pointwise (p, β)-hardy inequality u(x) Cd Ω(x)
More informationNon-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings. Krzysztof Burdzy University of Washington
Non-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings Krzysztof Burdzy University of Washington 1 Review See B and Kendall (2000) for more details. See also the unpublished
More informationPoisson random measure: motivation
: motivation The Lévy measure provides the expected number of jumps by time unit, i.e. in a time interval of the form: [t, t + 1], and of a certain size Example: ν([1, )) is the expected number of jumps
More informationMATH31011/MATH41011/MATH61011: FOURIER ANALYSIS AND LEBESGUE INTEGRATION. Chapter 2: Countability and Cantor Sets
MATH31011/MATH41011/MATH61011: FOURIER ANALYSIS AND LEBESGUE INTEGRATION Chapter 2: Countability and Cantor Sets Countable and Uncountable Sets The concept of countability will be important in this course
More informationON THE COMPLETE CONVERGENCE FOR WEIGHTED SUMS OF DEPENDENT RANDOM VARIABLES UNDER CONDITION OF WEIGHTED INTEGRABILITY
J. Korean Math. Soc. 45 (2008), No. 4, pp. 1101 1111 ON THE COMPLETE CONVERGENCE FOR WEIGHTED SUMS OF DEPENDENT RANDOM VARIABLES UNDER CONDITION OF WEIGHTED INTEGRABILITY Jong-Il Baek, Mi-Hwa Ko, and Tae-Sung
More informationCourse 212: Academic Year Section 1: Metric Spaces
Course 212: Academic Year 1991-2 Section 1: Metric Spaces D. R. Wilkins Contents 1 Metric Spaces 3 1.1 Distance Functions and Metric Spaces............. 3 1.2 Convergence and Continuity in Metric Spaces.........
More informationExponential martingales: uniform integrability results and applications to point processes
Exponential martingales: uniform integrability results and applications to point processes Alexander Sokol Department of Mathematical Sciences, University of Copenhagen 26 September, 2012 1 / 39 Agenda
More informationIGNACIO GARCIA, URSULA MOLTER, AND ROBERTO SCOTTO. (Communicated by Michael T. Lacey)
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 135, Number 10, October 2007, Pages 3151 3161 S 0002-9939(0709019-3 Article electronically published on June 21, 2007 DIMENSION FUNCTIONS OF CANTOR
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 informationFiltrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition
Filtrations, Markov Processes and Martingales Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition David pplebaum Probability and Statistics Department,
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 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 informationThe main results about probability measures are the following two facts:
Chapter 2 Probability measures The main results about probability measures are the following two facts: Theorem 2.1 (extension). If P is a (continuous) probability measure on a field F 0 then it has a
More informationWittmann Type Strong Laws of Large Numbers for Blockwise m-negatively Associated Random Variables
Journal of Mathematical Research with Applications Mar., 206, Vol. 36, No. 2, pp. 239 246 DOI:0.3770/j.issn:2095-265.206.02.03 Http://jmre.dlut.edu.cn Wittmann Type Strong Laws of Large Numbers for Blockwise
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 informationNotes on uniform convergence
Notes on uniform convergence Erik Wahlén erik.wahlen@math.lu.se January 17, 2012 1 Numerical sequences We begin by recalling some properties of numerical sequences. By a numerical sequence we simply mean
More informationErgodic Theorems. Samy Tindel. Purdue University. Probability Theory 2 - MA 539. Taken from Probability: Theory and examples by R.
Ergodic Theorems Samy Tindel Purdue University Probability Theory 2 - MA 539 Taken from Probability: Theory and examples by R. Durrett Samy T. Ergodic theorems Probability Theory 1 / 92 Outline 1 Definitions
More informationLecture 4 Lebesgue spaces and inequalities
Lecture 4: Lebesgue spaces and inequalities 1 of 10 Course: Theory of Probability I Term: Fall 2013 Instructor: Gordan Zitkovic Lecture 4 Lebesgue spaces and inequalities Lebesgue spaces We have seen how
More informationHEAVY-TRAFFIC EXTREME-VALUE LIMITS FOR QUEUES
HEAVY-TRAFFIC EXTREME-VALUE LIMITS FOR QUEUES by Peter W. Glynn Department of Operations Research Stanford University Stanford, CA 94305-4022 and Ward Whitt AT&T Bell Laboratories Murray Hill, NJ 07974-0636
More information