Search

We study the approximation of the ergodic measure of the following stochastic differential equation (SDE) on $\mathbb{R}^d$: \begin{eqnarray}\label{e:SDEE} d X_t &=& (b_1(X_t)+b_2(X_t)) d t+\sigma(X_t) d W_t, \end{eqnarray} where $W_t$ is a $d$-dimensional standard Brownian motion, and $b_1: \mathbb{R}^d \mapsto \mathbb{R}^d$, $b_2: \mathbb{R}^d \mapsto \mathbb{R}^d$ and $\sigma: \mathbb{R}^d \mapsto \mathbb{R}^{d\times d}$ are the functions to be specified in Assumption 2.1 below. In particular, $b_1$ satisfies $b_1\in \mathbb{L}^\infty(\mathbb{R}^d)\cap \mathbb{L}^1(\mathbb{R}^d)$ or $b_1 \in \mathcal{C}_b^{\alpha}(\mathbb{R}^d)$ with $\alpha\in (0,1)$, which makes the standard numerical schemes not work or fail to give a good convergence rate. In order to overcome these two difficulties, we first apply a Zvonkin's transform to SDE and obtain a new SDE which has coefficients with nice properties and admits a unique ergodic measure $\widehat \mu$, then discretize the new equation by Euler-Maruyama scheme to approximate $\widehat \mu$, and finally use the inverse Zvonkin's transform to get an approximation of the ergodic measure of SDE, denoted by $\mu$. Our approximation method is inspired by Xie and Zhang [22]. The proof of our main result is based on the method of introducing a stationary Markov chain, a key ingredient in this method is establishing the regularity of a Poisson equation, which is done by combining the classical PDE local regularity and a nice extension trick introduced by Gurvich [10].

We are interested in the following two $\mathbb{R}^d$-valued stochastic differential equations (SDEs): \begin{gather*} d X_t=b(X_t)\, d t + \sigma\,d L_t, \quad X_0=x, %\label{BM-SDE} d Y_t=b(Y_t)\,dt + \sigma\,d B_t, \quad Y_0=y, \end{gather*} where $\sigma$ is an invertible $d\times d$ matrix, $L_t$ is a rotationally symmetric $\alpha$-stable L\'evy process, and $B_t$ is a $d$-dimensional standard Brownian motion. We show that for any $\alpha_0 \in (1,2)$ the Wasserstein-$1$ distance $W_1$ satisfies for $\alpha \in [\alpha_0,2)$ \begin{gather*} W_{1}\left(X_{t}^x, Y_{t}^y\right) \leq C e^{-Ct}|x-y| +C_{\alpha_0}d\cdot\log(1+d)(2-\alpha)\log\frac{1}{2-\alpha}, \end{gather*} which implies, in particular, \normal \begin{equation} \label{e:W1Rate} W_1(\mu_\alpha, \mu) \leq C_{\alpha_0} d \cdot \log(1+d)(2-\alpha) \log \frac{1}{2-\alpha}, \end{equation} where $\mu_\alpha$ and $\mu$ are the ergodic measures of $(X^x_t)_{t \ge 0}$ and $(Y^y_t)_{t \ge 0}$ respectively. The term $d\cdot\log(1+d)$ appearing in this estimate seems to be optimal. For the special case of a $d$-dimensional Ornstein--Uhlenbeck system, we show that $W_1(\mu_\alpha, \mu) \geq C_{\alpha_0,d} (2-\alpha)$; this indicates that the convergence rate with respect to $\alpha$ in \eqref{e:W1Rate} is optimal up to a logarithmic correction. We conjecture that the sharp rate with respect to $\alpha$ and $d$ is $d\cdot\log (1+d) (2-\alpha)$.

In this paper, we develop a stochastic algorithm based on the Euler--Maruyama scheme to approximate the invariant measure of the limiting multidimensional diffusion of $G/Ph/n+GI$ queues in the Halfin-Whitt regime. Specifically, we prove a non-asymptotic error bound between the invariant measures of the approximate model from the algorithm and the limiting diffusion. To establish the error bound, we employ the recently developed Stein's method for multi-dimensional diffusions, in which the regularity of Stein's equation developed by Gurvich (2014, 2022) plays a crucial role. We further prove the central limit theorem (CLT) and the moderate deviation principle (MDP) for the occupation measures of the limiting diffusion of $G/Ph/n+GI$ queues and its Euler-Maruyama scheme. In particular, the variances in the CLT and MDP associated with the limiting diffusion are determined by Stein's equation and Malliavin calculus, in which properties of a mollified diffusion and an associated weighted occupation time play a crucial role., Comment: arXiv admin note: substantial text overlap with arXiv:2109.03623

We propose two Euler-Maruyama (EM) type numerical schemes in order to approximate the invariant measure of a stochastic differential equation (SDE) driven by an $\alpha$-stable L\'evy process ($1, Comment: Accepted by SPA. This version is a little different version from the published one, in which we add in the appendix a shorter proof for a Bismut-Elworthy-Li formula suggested by a very careful and knowledgeable referee

We consider a stochastic 2D Navier-Stokes equation in a bounded domain. The random force is assumed to be non-degenerate and periodic in time, its law has a support localised with respect to both time and space. Slightly strengthening the conditions in the pioneering work about exponential ergodicity by Shirikyan [Shi15], we prove that the stochastic system satisfies Donsker-Varadhan typle large deviations principle. Our proof is based on a criterion of [JNPS15] in which we need to verify uniform irreducibility and uniform Feller property for the related Feynman-Kac semigroup.

Inspired by \citet{Berkes14} and \citet{Wu07}, we prove an almost sure invariance principle for stationary $β-$mixing stochastic processes defined on Hilbert space. Our result can be applied to Markov chain satisfying Meyn-Tweedie type Lyapunov condition and thus generalises the contraction condition in \citet[Example 2.2]{Berkes14}. We prove our main theorem by the big and small blocks technique and an embedding result in \citet{gotze2011estimates}. Our result is further applied to the ergodic Markov chain and functional autoregressive processes.

Letting~$N=\left\{N(t), t\geq0\right\}$ be a standard Poisson process, Stroock~ \cite{Stroock-1981} constructed a family of continuous processes by $$\Theta_{\epsilon}(t)=\int_0^t\theta_{\epsilon}(r)dr, \ \ \ \ \ 0 \le t \le 1,$$ where $\theta_{\epsilon}(r)=\frac{1}{\epsilon}(-1)^{N(\epsilon^{-2}r)}$, and proved that it weakly converges to a standard Brownian motion under the continuous function topology. We establish the functional large deviations principle (LDP) for the approximations of a class of Gaussian processes constructed by integrals over $\Theta_{\epsilon}(t)$, and find the explicit form for rate function. As an application, we consider the following (non-Markovian) stochastic differential equation \begin{equation*} \begin{aligned} X^{\epsilon}(t) &=x_{0}+\int^{t}_{0}b(X^{\epsilon}(s))ds+\lambda(\epsilon)\int^{t}_{0}\sigma(X^{\epsilon}(s))d\Theta_{\epsilon}(s), \end{aligned} \end{equation*} where $b$ and $\sigma$ are both Lipschitz functions, and establish its Freidlin-Wentzell type LDP as $\epsilon \rightarrow 0$. The rate function indicates a phase transition phenomenon as $\lambda(\epsilon)$ moves from one region to the other.

In this paper, we establish normalized and self-normalized Cramér-type moderate deviations for Euler-Maruyama scheme for SDE. As a consequence of our results, Berry-Esseen's bounds and moderate deviation principles are also obtained. Our normalized Cramér-type moderate deviations refines the recent work of [Lu, J., Tan, Y., Xu, L., 2022. Central limit theorem and self-normalized Cramér-type moderate deviation for Euler-Maruyama scheme. Bernoulli 28(2): 937--964]., 21 pages

There has been a surge of interest in developing robust estimators for models with heavy-tailed and bounded variance data in statistics and machine learning, while few works impose unbounded variance. This paper proposes two type of robust estimators, the ridge log-truncated M-estimator and the elastic net log-truncated M-estimator. The first estimator is applied to convex regressions such as quantile regression and generalized linear models, while the other one is applied to high dimensional non-convex learning problems such as regressions via deep neural networks. Simulations and real data analysis demonstrate the {robustness} of log-truncated estimations over standard estimations., Comment: 44 pages

We generalize the { ${\rm M}$-estimator} put forward by Catoni in his seminal paper [C12] to the case in which samples can have finite $��$-th moment with $��\in (1,2)$ rather than finite variance, our approach is by slightly modifying the influence function $��$ therein. The choice of the new influence function is inspired by the Taylor-like expansion developed in [C-N-X]. We obtain a deviation bound of the estimator, as $��\rightarrow 2$, this bound is the same as that in [C12]. Experiment shows that our generalized ${\rm M}$-estimator performs better than the empirical mean estimator, the smaller the $��$ is, the better the performance will be. As an application, we study an $\ell_{1}$ regression considered by Zhang et al. [Z-Z] who assumed that samples have finite variance, and relax their assumption to be finite {$��$-th} moment with $��\in (1,2)$.

We study stochastic integrals driven by a general subordinator and establish a zero-one law for the finiteness of the resulting integral as well as moment estimates. As an application, we use these results to obtain structural properties of SPDEs driven by multiplicative pure jump noise, which include (1) a maximal inequality for a multiplicative stochastic convolution $Z_t$, (2) a small ball probability of $Z_t$, (3) the existence of invariant measures and accessibility to zero of SPDEs, and (4) a Galerkin approximation of solutions to SPDEs., Comment: Accepted by TAMS

We consider a stochastic differential equation and its Euler-Maruyama (EM) scheme, under some appropriate conditions, they both admit a unique invariant measure, denoted by $��$ and $��_��$ respectively ($��$ is the step size of the EM scheme). We construct an empirical measure $��_��$ of the EM scheme as a statistic of $��_��$, and use Stein's method developed in \citet{FSX19} to prove a central limit theorem of $��_��$. The proof of the self-normalized Cram��r-type moderate deviation (SNCMD) is based on a standard decomposition on Markov chain, splitting $��^{-1/2}(��_��(.)-��(.))$ into a martingale difference series sum $\mcl H_��$ and a negligible remainder $\mcl R_��$. We handle $\mcl H_��$ by the time-change technique for martingale, while prove that $\mcl R_��$ is exponentially negligible by concentration inequalities, which have their independent interest. Moreover, we show that SNCMD holds for $x = o(��^{-1/6})$, which has the same order as that of the classical result in \citet{shao1999cramer,JSW03}., To appear in Bernoulli, correct a small error in the proof of Lemma 4.1

Stein's method has been widely used for probability approximations. However, in the multi-dimensional setting, most of the results are for multivariate normal approximation or for test functions with bounded second- or higher-order derivatives. For a class of multivariate limiting distributions, we use Bismut's formula in Malliavin calculus to control the derivatives of the Stein equation solutions by the first derivative of the test function. Combined with Stein's exchangeable pair approach, we obtain a general theorem for multivariate approximations with near optimal error bounds on the Wasserstein distance.We apply the theorem to the unadjusted Langevin algorithm., Comment: We changed the title and corrected some small errors and typos

This paper is concerned with the Stein's method associated with a (possibly) asymmetric $��$-stable distribution $Z$, in dimension one. More precisely, its goal is twofold. In the first part, we exhibit a genuine bound for the Wasserstein distance between $Z$ and any integrable random variable $X$, in terms of an operator that reduces to the classical fractional Laplacian in the symmetric case. Then, in the second part we apply the aforementioned bound to compute error rates in the stable central limit theorem, when the entries are in the domain $\mathcal{D}_��$ of normal attraction of a stable law of exponent $��$. To conclude, we study the specific case where the entries are Pareto like multiplied by a slowly varying function, which provides an example of random variables that do not belong to $\mathcal{D}_��$, but for which our approach continues to apply., 14 pages. A new example where the entries are not in the domain of normal attraction of a stable law of exponent $\alpha$ has been added

By the Lindeberg principle, we develop in this paper an approximation to one dimensional (possibly) asymmetric $\alpha$-stable distributions with $\alpha \in (0,2)$ in the smooth Wasserstein distance. It is the first time that the general stable central limit theorem is proved by the Lindeberg principle, and that this theorem with $\alpha \in (0,1]$ is proved by a new method other than Fourier analysis. Our main tools are a Taylor-like expansion and a Kolmogorov forward equation.

Using heat kernel estimates, we prove the pathwise uniqueness for strong solutions of irregular stochastic differential equation driven by a family of Markov process, whose generator is a non-local and non-symmetric L\'evy type operator. Due to the extra term $1_{[0,\sigma(X_{s-},z)]}(r)$ in multiplicative noise, we need to derive some new regularity results for the generator and use a trick of mixing $L_1$ and $L_2$-estimates by Kurtz and Protter \cite{Ku-Po}.

By a coupling method, we prove that a family of stochastic partial differential equations (SPDEs) driven by highly degenerate pure jump L��vy noises are exponential mixing. These pure jump L��vy noises include $��$-stable process with $��\in (0,2)$

We establish a large deviation principle for the occupation measure of the stochastic real Ginzburg-Landau equation driven by $��$-stable noises. The proof is based on a hyper-exponential recurrence criterion. Our result indicates a phenomenon that strong dissipation beats heavy tailed noises to produce a large deviation, it seems to us that this phenomenon has not been reported in the known literatures., This paper fills the gap in the paper arXiv:1410.7247 which has been withdrawn. We complete the proof in arXiv:1410.7247. arXiv admin note: text overlap with arXiv:1510.01904

Let $(P_t)$ be the transition semigroup of a L��vy process $L$ taking values in a Hilbert space $H$. Let $��$ be the L��vy measure of $L$. It is shown that for any bounded and measurable function $f$, $$ \int_H\left\vert P_tf(x+y)-P_tf(x)\right\vert ^2 ��(\dif y)\le \frac 1 t P_tf^2(x) \qquad \text{for all $t>0$, $x\in H$.} $$ As $��$ can be infinite this formula establishes some smoothening effect of the semigroup $(P_t)$. In the paper some applications of the formula will be presented as well.

Using the Bismut's approach to Malliavin calculus, we introduce a simplified Malliavin matrix ([11]) for stochastic differential equations (SDEs) force by degenerate stable like noises. For the degenerate SDEs driven by Wiener noises, one can derive a Norris type lemma and use it \emph{iteratively} to prove the smoothness of density functions. Unfortunately, Norris type lemma is very hard to be iteratively applied to SDEs with stable like noises. In this paper, we derive a simple inequality as a replacement and use it to show that two families of degenerate SDEs with stable like noises admit smooth density functions. One family is the linear SDEs studied by Priola and Zabczyk ([13]), under some additional assumption we can iteratively use the inequality to get the smoothness of the density. The other family is the general SDEs with stable like noises, we can apply this inequality only \emph{one} time and thus derive that the SDEs admit smooth density if the \emph{first} order Lie brackets span $\R^d$. The crucial step in this paper is estimating the smallest eigenvalue of the simplified Malliavin matrix, which only uses some elementary facts of Poisson processes and undergraduate level ordinary differential equations., Comment: Added some references by Bismut which include integration by parts of SDEs forced by stable like noises, Rewrited the abstract, Corrected some very small errors

Gradient estimates are derived, for the first time, for the semigroup associated to a class of stochastic differential equations driven by multiplicative L��vy noise. In particular, the estimates are sharp for $��$-stable type noises. To derive these estimates, a new derivative formula of Bismut-Elworthy-Li's type is established for the semigroup by using the Malliavin calculus and a finite-jump approximation argument., 21pp

Our simple but useful technique is using an integration by parts to split the stochastic convolution into two terms. We develop five applications for this technique. The first one is getting a uniform estimate of stochastic convolution of $��$-stable processes. Since $��$-stable noises only have $p, 12 pages

We study the ergodicity of stochastic real Ginzburg-Landau equation driven by additive $��$-stable noises, showing that as $��\in (3/2,2)$, this stochastic system admits a unique invariant measure. After establishing the existence of invariant measures by the same method as in [9], we prove that the system is strong Feller and accessible to zero. These two properties imply the ergodicity by a simple but useful criterion in [16]. To establish the strong Feller property, we need to truncate the nonlinearity and apply a gradient estimate established in [26] (or see [24]} for a general version for the finite dimension systems). Because the solution has discontinuous trajectories and the nonlinearity is not Lipschitz, we can not solve a control problem to get irreducibility. Alternatively, we use a replacement, i.e., the fact that the system is accessible to zero. In section 3, we establish a maximal inequality for stochastic $��$-stable convolution, which is crucial for studying the well-posedness, strong Feller property and the accessibility of the mild solution. We hope this inequality will also be useful for studying other SPDEs forced by $��$-stable noises., 25 pages. Gave a new estimate in Lemma 3.1 which is more delicate than the original, corrected some typos in the proof of Lemma 3.1

In this work, we prove the strong Feller property and the exponential ergodicity of stochastic Burgers equations driven by $��/2$-subordinated cylindrical Brownian motions with $��\in(1,2)$. To prove the results, we truncate the nonlinearity and use the derivative formula for SDEs driven by $��$-stable noises established in Zhang (arXiv:1204.2630v2)., 17pp

In this note we prove the well-posedness for stochastic 2D Navier-Stokes equation driven by general L��vy processes (in particular, $��$-stable processes), and obtain the existence of invariant measures., 9 pages

We modify the coupling method established in [22, 20] and develop a technique to prove the exponential mixing of a 2D stochastic system forced by degenerate Levy noises. In particular, these Levy noises include $\alpha$-stable noises (0 < $\alpha$ < 2). Thanks to the stimulating discussion [14], this technique is promising to study the exponential mixing problem of SPDEs driven by degenerate symmetric $\alpha$-stable noises., Comment: Some small errors in the proof of Theorem 2.4 were corrected

We show how gradient estimates for transition semigroups can be used to establish exponential mixing for a class of Markov processes in infinite dimensions. We concentrate on semilinear systems driven by cylindrical $��$-stable noises, $��\in (0,2)$, introduced in Priola-Zabczyk "Structural properties of semilinear SPDEs driven by cylindrical stable processes" (PTRF to appear). We first prove that if the nonlinearity is bounded, then the system is ergodic and strong mixing. Then we show that the system is exponentially mixing provided that the nonlinearity, or its Lipschitz constant, are sufficiently small.

We study an infinite white $��$-stable systems with unbounded interactions, proving the existence by Galerkin approximation and exponential mixing property by an $��$-stable version of gradient bounds., 28 pages (Corrected several errors)

We proved the existence of an infinite dimensional stochastic system driven by white $��$-stable noises ($1, 20 pages