Convergence in Density of Splitting AVF Scheme for Stochastic Langevin Equation
Jianbo Cui, Jialin Hong, Derui Sheng
TL;DR
The paper investigates convergence in density for numerical approximations of a stochastic Langevin equation with non-globally monotone drift using a splitting AVF scheme. It develops a Malliavin-calculus framework to prove exponential integrability and the existence and smoothness of densities for both the exact and numerical solutions, and establishes uniform non-degeneracy of the Malliavin covariance. The authors obtain an optimal strong convergence rate of order $h^{1/2}$ in the Malliavin–Sobolev sense and, leveraging Donsker’s delta, prove a density-convergence rate of order $O(h)$ for the numerical scheme. These results bridge strong-path convergence with probabilistic density convergence under Hörmander-type nondegeneracy, enabling reliable density approximations for degenerate SDEs with super-linear drift in practical simulations.
Abstract
In this article, we study the density function of the numerical solution of the splitting averaged vector field (AVF) scheme for the stochastic Langevin equation. To deal with the non-globally monotone coefficient in the considered equation, we first present the exponential integrability properties of the exact and numerical solutions. Then we show the existence and smoothness of the density function of the numerical solution by proving its uniform non-degeneracy in Malliavin sense. In order to analyze the approximate error between the density function of the exact solution and that of the numerical solution, we derive the optimal strong convergence rate in every Malliavin--Sobolev norm of the numerical scheme via Malliavin calculus. Combining the approximation result of Donsker's delta function and the smoothness of the density functions, we prove that the convergence rate in density coincides with the optimal strong convergence rate of the numerical scheme.
