Weak Error on the densities for the Euler scheme of stable additive SDEs with Besov drift
Mathis Fitoussi, Elena Issoglio, Stéphane Menozzi
TL;DR
The paper analyzes Euler discretization for multi-dimensional SDEs driven by symmetric α-stable noise with a distributional Besov drift. By leveraging a Besov-space framework, heat-kernel estimates, and Duhamel representations, the authors derive a weak-density convergence rate of (γ−ε)/α, where γ>0 encodes the drift-regularity and Krylov-type conditions. A careful error decomposition into six Δ-terms, combined with controlled Hölder-modulus bounds and a Grönwall-type argument, yields the rate and clarifies how parameter choices (α,β,p,q,r) influence convergence. The results extend weak-error-rate insights from Brownian and Lévy-driven settings to the challenging Besov-drift, multi-dimensional, strictly stable case, with potential extensions toBrownian noise under analogous heat-kernel estimates. The framework provides robust density-weak convergence results for discretizations under rough drifts, relevant for numerical approximation of stable-driven SDEs in high-dimensional settings.
Abstract
We are interested in the Euler-Maruyama dicretization of the formal SDE, $dX_t=b(t,X_t)dt+dZ_t$, where $Z$ is a symmetric isotropic d dimensional stable process of index $α\in (1,2)$, and $b$ is distributional. It belongs to a mix Lebesgue-Besov space. The associated parameters satisfy some constraints which guarantee weak-well posedness. Defining an appropriate Euler scheme, we obtain a convergence rate for the weak error on the densities. The rate depends on the parameters.
