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.

Stability estimates for initial data in general Ornstein-Uhlenbeck equations

Stability estimate for a partial data inverse problem for the convection-diffusion equation

Stability and moment estimates for the stochastic singular $Φ$-Laplace equation

Measure data systems with Orlicz growth

On log-concave approximations of high-dimensional posterior measures and stability properties in non-linear inverse problems

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

Large deviations and concentration inequalities for the Ornstein-Uhlenbeck process without tears

Invariant Measures and Maximal L^2 Regularity for Nonautonomous Ornstein-Uhlenbeck Equations

Stability Estimates for the Inverse Fractional Conductivity Problem

Hölder stability for a semilinear elliptic inverse problem

Large-time behavior of solutions of parabolic equations on the real line with convergent initial data III:unstable limit at infinity

Invariant Measures of Stochastic Lattice Plate Equations: Stability, Ergodicity and Mixing

Stability Estimates for Some Parabolic Inverse Problems With the Final Overdetermination via a New Carleman Estimate

Weighted $L^1$-semigroup approach for nonlinear Fokker--Planck equations and generalized Ornstein--Uhlenbeck processes

Regularity of semi-stable solutions to fourth order nonlinear eigenvalue problems on general domains

Schauder estimates for degenerate stable Kolmogorov equations

On ergodic invariant measures for the stochastic Landau-Lifschitz-Gilbert equation in 1D

Exponential stability of the nonlinear Schrödinger equation with locally distributed damping on compact Riemannian manifold

Heavy tailed large deviations for time averages of diffusions: the Ornstein-Uhlenbeck case

Stable nearly self-similar blowup of the 2D Boussinesq and 3D Euler equations with smooth data II: Rigorous Numerics