Tempering stable processes

A tempered stable Lévy process combines both the α$\alpha$-stable and Gaussian trends. In a short time frame it is close to an α$\alpha$-stable process while in a long time frame it approximates a Brownian motion. In this paper we consider a general and robust class of multivariate tempered stable distributions and establish their identifiable parametrization. We prove short and long time behavior of tempered stable Lévy processes and investigate their absolute continuity with respect to the underlying α$\alpha$-stable processes. We find probabilistic representations of tempered stable processes which specifically show how such processes are obtained by cutting (tempering) jumps of stable processes. These representations exhibit α$\alpha$-stable and Gaussian tendencies in tempered stable processes and thus give probabilistic intuition for their study. Such representations can also be used for simulation. We also develop the corresponding representations for Ornstein–Uhlenbeck-type processes.     

Limit theorems for occupation time fluctuations of branching systems I: Long-range dependence

We give a functional limit theorem for the fluctuations of the rescaled occupation time process of a critical branching particle system in R^d${\mathbb{R}}^d$ with symmetric α$\alpha$-stable motion and α<d<2α$\alpha <d<2\alpha$, which leads to a long-range dependence process involving sub-fractional Brownian motion. We also give an analogous result for the system without branching and d<α$d<\alpha$, which involves fractional Brownian motion. We use a space–time random field approach.     

Weak convergence of censored and reflected stable processes

It is shown that if a sequence of open n$n$-sets _Dk$D_k$ increases to an open n$n$-set D$D$ then reflected stable processes in _Dk$D_k$ converge weakly to the reflected stable process in D$D$ for every starting point x$x$ in D$D$. The same result holds for censored α$\alpha$-stable processes for every x$x$ in D$D$ if D$D$ and _Dk$D_k$ satisfy the uniform Hardy inequality. Using the method in the proof of the above results, we also prove the weak convergence of reflected Brownian motions in unbounded domains.     

Ergodicity for time-changed symmetric stable processes

In this paper we study ergodicity and related semigroup property for a class of symmetric Markov jump processes associated with time-changed symmetric α$\alpha$-stable processes. For this purpose, explicit and sharp criteria for Poincaré type inequalities (including Poincaré, super Poincaré and weak Poincaré inequalities) of the corresponding non-local Dirichlet forms are derived. Moreover, our main results, when applied to a class of one-dimensional stochastic differential equations driven by symmetric α$\alpha$-stable processes, yield sharp criteria for their various ergodic properties and corresponding functional inequalities.     

Limit theorems for occupation time fluctuations of branching systems II: Critical and large dimensions

We give functional limit theorems for the fluctuations of the rescaled occupation time process of a critical branching particle system in R^d${\mathbb{R}}^d$ with symmetric α$\alpha$-stable motion in the cases of critical and large dimensions, d=2α$d=2\alpha$ and d>2α$d>2\alpha$. In a previous paper [T. Bojdecki, L.G. Gorostiza, A. Talarczyk, Limit theorems for occupation time fluctuations of branching systems I: long-range dependence, Stochastic Process. Appl., this issue.] we treated the case of intermediate dimensions, α<d<2α$\alpha <d<2\alpha$, which leads to a long-range dependence limit process. In contrast, in the present cases the limits are generalized Wiener processes. We use the same space–time random field method of the previous paper, the main difference being that now the tightness requires a new approach and the proofs are more difficult. We also give analogous results for the system without branching in the cases d=α$d=\alpha$ and d>α$d>\alpha$.     

On functional limits of short- and long-memory linear processes with GARCH(1,1) noises

This paper considers the short- and long-memory linear processes with GARCH (1,1) noises. The functional limit distributions of the partial sum and the sample autocovariances are derived when the tail index α$\alpha$ is in (0,2)$(0,2)$, equal to 2, and in (2,∞)$(2,\infty )$, respectively. The partial sum weakly converges to a functional of α$\alpha$-stable process when α<2$\alpha <2$ and converges to a functional of Brownian motion when α≥2$\alpha \ge 2$. When the process is of short-memory and α<4$\alpha <4$, the autocovariances converge to functionals of α/2$\alpha /2$-stable processes; and if α≥4$\alpha \ge 4$, they converge to functionals of Brownian motions. In contrast, when the process is of long-memory, depending on α$\alpha$ and β$\beta$ (the parameter that characterizes the long-memory), the autocovariances converge to either (i) functionals of α/2$\alpha /2$-stable processes; (ii) Rosenblatt processes (indexed by β$\beta$, 1/2<β<3/4$1/2<\beta <3/4$); or (iii) functionals of Brownian motions. The rates of convergence in these limits depend on both the tail index α$\alpha$ and whether or not the linear process is short- or long-memory. Our weak convergence is established on the space of càdlàg functions on [0,1]$[0,1]$ with either (i) the _J1$J_1$ or the _M1$M_1$ topology (Skorokhod, 1956); or (ii) the weaker form S$S$ topology (Jakubowski, 1997). Some statistical applications are also discussed.     

Simulation of the continuous time random walk of the space-fractional diffusion equations

In this article, we discuss the solution of the space-fractional diffusion equation with and without central linear drift in the Fourier domain and show the strong connection between it and the α$\alpha$-stable Lévy distribution, 0<α<2$0<\alpha <2$. We use some relevant transformations of the independent variables x$x$ and t$t$, to find the solution of the space-fractional diffusion equation with central linear drift which is a special form of the space-fractional Fokker–Planck equation which is useful in studying the dynamic behaviour of stochastic differential equations driven by the non-Gaussian (Lévy) noises. We simulate the continuous time random walk of these models by using the Monte Carlo method.     

Tail behavior of random products and stochastic exponentials

In this paper we study the distributional tail behavior of the solution to a linear stochastic differential equation driven by infinite variance α$\alpha$-stable Lévy motion. We show that the solution is regularly varying with index α$\alpha$. An important step in the proof is the study of a Poisson number of products of independent random variables with regularly varying tail. The study of these products merits its own interest because it involves interesting saddle-point approximation techniques.     

A short proof of the Doob–Meyer theorem

Every submartingale S$S$ of class D$D$ has a unique Doob–Meyer decomposition S=M+A$S=M+A$, where M$M$ is a martingale and A$A$ is a predictable increasing process starting at 0.     

Asymptotic results for the two-parameter Poisson–Dirichlet distribution

The two-parameter Poisson–Dirichlet distribution is the law of a sequence of decreasing nonnegative random variables with total sum one. It can be constructed from stable and gamma subordinators with the two parameters, α$\alpha$ and θ$\theta$, corresponding to the stable component and the gamma component respectively. The moderate deviation principle is established for the distribution when θ$\theta$ approaches infinity, and the large deviation principle is established when both α$\alpha$ and θ$\theta$ approach zero.     

Monotonicity of the reflected Bessel transition density on the diagonal

For α∈R$\alpha \in \mathbb{R}$, let p^R(t,x,x)$p^R(t,x,x)$ denote the diagonal of the transition density of the α$\alpha$-Bessel process in (0,1]$(0,1]$, killed at 0 and reflected at 1. As a function of x$x$, if either α≥3$\alpha \ge 3$ or α=1$\alpha =1$, then for t>0$t>0$, the diagonal is nondecreasing. This monotonicity property fails if 1≠α<3$1\ne \alpha <3$.     

Oscillation of harmonic functions for subordinate Brownian motion and its applications

In this paper, we establish an oscillation estimate of nonnegative harmonic functions for a pure-jump subordinate Brownian motion. The infinitesimal generator of such subordinate Brownian motion is an integro-differential operator. As an application, we give a probabilistic proof of the following form of relative Fatou theorem for such subordinate Brownian motion X$X$ in a bounded κ$\kappa$-fat open set; if u$u$ is a positive harmonic function with respect to X$X$ in a bounded κ$\kappa$-fat open set D$D$ and h$h$ is a positive harmonic function in D$D$ vanishing on D^c$D^c$, then the non-tangential limit of u/h$u/h$ exists almost everywhere with respect to the Martin-representing measure of h$h$.     

Smale’s point estimate theory for Newton’s method on Lie groups

We introduce the notion of the (one-parameter subgroup) γ$\gamma$-condition for a map f$f$ from a Lie group to its Lie algebra and establish α$\alpha$-theory and γ$\gamma$-theory for Newton’s method for a map f$f$ satisfying this condition. Applications to analytic maps are provided, and Smale’s α$\alpha$-theory and γ$\gamma$-theory are extended and developed. Examples arising from initial value problems on Lie group are presented to illustrate applications of our results.     

Phase transition in equilibrium fluctuations of symmetric slowed exclusion

We analyze the equilibrium fluctuations of density, current and tagged particle in symmetric exclusion with a slow bond. The system evolves in the one-dimensional lattice and the jump rate is everywhere equal to one except at the slow bond where it is αn^−β$\alpha n^{-\beta }$, with α>0$\alpha >0$, β∈[0,+∞]$\beta \in [0,+\infty ]$ and n$n$ is the scaling parameter. Depending on the regime of β$\beta$, we find three different behaviors for the limiting fluctuations whose covariances are explicitly computed. In particular, for the critical value β=1$\beta =1$, starting a tagged particle near the slow bond, we obtain a family of Gaussian processes indexed in α$\alpha$, interpolating a fractional Brownian motion of Hurst exponent 1/4$1/4$ and the degenerate process equal to zero.     

The scaling limit of Poisson-driven order statistics with applications in geometric probability

Let _ηt${\eta }_t$ be a Poisson point process of intensity t≥1$t\ge 1$ on some state space Y$\mathbb{Y}$ and let f$f$ be a non-negative symmetric function on Y^k${\mathbb{Y}}^k$ for some k≥1$k\ge 1$. Applying f$f$ to all k$k$-tuples of distinct points of _ηt${\eta }_t$ generates a point process _ξt${\xi }_t$ on the positive real half-axis. The scaling limit of _ξt${\xi }_t$ as t$t$ tends to infinity is shown to be a Poisson point process with explicitly known intensity measure. From this, a limit theorem for the m$m$-th smallest point of _ξt${\xi }_t$ is concluded. This is strengthened by providing a rate of convergence. The technical background includes Wiener–Itô chaos decompositions and the Malliavin calculus of variations on the Poisson space as well as the Chen–Stein method for Poisson approximation. The general result is accompanied by a number of examples from geometric probability and stochastic geometry, such as k$k$-flats, random polytopes, random geometric graphs and random simplices. They are obtained by combining the general limit theorem with tools from convex and integral geometry.     

Multi-dimensional G-Brownian motion and related stochastic calculus under G-expectation

We develop a notion of nonlinear expectation–G$G$-expectation–generated by a nonlinear heat equation with infinitesimal generator G$G$. We first study multi-dimensional G$G$-normal distributions. With this nonlinear distribution we can introduce our G$G$-expectation under which the canonical process is a multi-dimensional G$G$-Brownian motion. We then establish the related stochastic calculus, especially stochastic integrals of Itô’s type with respect to our G$G$-Brownian motion, and derive the related Itô’s formula. We have also obtained the existence and uniqueness of stochastic differential equations under our G$G$-expectation.