Strong backward error analysis of symplectic integrators for stochastic Hamiltonian systems

Abstract

Backward error analysis is a powerful tool in order to detect the long-term conservative behavior of numerical methods. In this work, we present a long-term analysis of symplectic stochastic numerical integrators, applied to Hamiltonian systems with multiplicative noise. We first compute and analyze the associated stochastic modified differential equations. Then, suitable bounds for the coefficients of such equations are provided towards the computation of long-term estimates for the Hamiltonian deviations occurring along the aforementioned numerical dynamics. This result generalizes Benettin-Giorgilli Theorem to the scenario of stochastic symplectic methods. Finally, specific numerical methods are considered, in order to provide a numerical evidence confirming the effectiveness of the theoretical investigation.