Splitting AVF method for generalized Langevin equations: probability density function and geometric ergodicity

Abstract

The generalized Langevin equation (GLE) constitutes a fundamental model for describing nonequilibrium dynamics with memory effects. To overcome the numerical challenges arising from superquadratically growing potentials and degenerate noise, we propose and analyze a structure-preserving splitting averaged vector field (AVF) method for a quasi-Markovian GLE. The core advantage of this method lies in its ability to simultaneously preserve the exponential integrability, Malliavin differentiability, and ergodicity of the underlying continuous system. Notably, by leveraging exponential integrability, Malliavin differentiability, and uniform non-degeneracy of the numerical solution, we obtain the existence and smoothness of its probability density function, which converges to that of the exact solution with first-order accuracy. Furthermore, by validating the Lyapunov condition and the minorization condition using a localized technique, we establish the geometric ergodicity of the numerical solution. Finally, numerical experiments are conducted to confirm the theoretical results.