Abstract:Let $(X_t)_{t \ge 0}$ be the solution of the stochastic differential equation $$dX_t = b(X_t) dt+A dZ_t, \quad X_{0}=x,$$ where $b: \mathbb{R}^d \rightarrow \mathbb R^d$ is a Lipschitz function, $A \in \mathbb R^{d \times d}$ is a positive definite matrix, $(Z_t)_{t\geq 0}$ is a $d$-dimensional rotationally invariant $\alpha$-stable Lévy process with $\alpha \in (1,2)$ and $x\in\mathbb{R}^{d}$. We use two Euler-Maruyama schemes with decreasing step sizes $\Gamma = (\gamma_n)_{n\in \mathbb{N}}$ to approximate the invariant measure of $(X_t)_{t \ge 0}$: one with i.i.d. $\alpha$-stable distributed random variables as its innovations and the other with i.i.d. Pareto distributed random variables as its innovations. We study the convergence rate of these two approximation schemes in the Wasserstein-1 distance. For the first scheme, when the function $b$ is Lipschitz and satisfies a certain dissipation condition, we show that the convergence rate is $\gamma^{1/\alpha}_n$. Under an additional assumption on the second order directional derivatives of $b$, this convergence rate can be improved to $\gamma^{1+\frac 1 {\alpha}-\frac{1}{\kappa}}_n$ for any $\kappa \in [1,\alpha)$. For the second scheme, when the function $b$ is twice continuously differentiable, we obtain a convergence rate of $\gamma^{\frac{2-\alpha}{\alpha}}_n$. We show that the rate $\gamma^{\frac{2-\alpha}{\alpha}}_n$ is optimal for the one dimensional stable Ornstein-Uhlenbeck process. Our theorems indicate that the recent remarkable result about the unadjusted Langevin algorithm with additive innovations can be extended to the SDEs driven by an $\alpha$-stable Lévy process and the corresponding convergence rate has a similar behaviour. Compared with the previous result, we have relaxed the second order differentiability condition to the Lipschitz condition for the first scheme.

High order weak approximation of Stochastic Differential Equations for bounded and measurable test functions

Weak Backward Error Analysis for SDEs

Weak approximation schemes for SDEs with super-linearly growing coefficients

Boundary-preserving weak approximations of some semilinear stochastic partial differential equations

Weak Approximations of Stochastic Partial Differential Equations with Fractional Noise

Optimal Weak Order and Approximation of the Invariant Measure with a Fully-Discrete Euler Scheme for Semilinear Stochastic Parabolic Equations with Additive Noise

Strongly asymptotically optimal schemes for the strong approximation of stochastic differential equations with respect to the supremum error

Weak Convergence of Euler Scheme for SDEs with Low Regular Drift

Approximation of the invariant measure for stable SDE by the Euler-Maruyama scheme with decreasing step-sizes

Approximation of SDEs: a stochastic sewing approach

On the convergence of adaptive approximations for stochastic differential equations

A new weak approximation scheme of stochastic differential equations and the Runge-Kutta method

Weak well-posedness and weak discretization error for stable-driven SDEs with Lebesgue drift

Weak convergence of Galerkin approximations of stochastic partial differential equations driven by additive Lévy noise

Approximation of the invariant measure of stable SDEs by an Euler--Maruyama scheme

Analysis of a Modified Regularity-Preserving Euler Scheme for Parabolic Semilinear SPDEs: Total Variation Error Bounds for the Numerical Approximation of the Invariant Distribution

A generalized coupling approach for the weak approximation of stochastic functional differential equations

Weak convergence and invariant measure of a full discretization for non-globally Lipschitz parabolic SPDE

Higher order numerical methods for SDEs without globally monotone coefficients

Weak Convergence of Euler-Maruyama's Approximation for SDEs under Integrability Condition

Strong Approximation of Monotone Stochastic Partial Differential Equations Driven by Multiplicative Noise