Backward error analysis of stochastic Poisson integrators

TL;DR

This paper performs backward error analysis on stochastic Poisson integrators to evaluate their long-term structure-preserving properties, providing theoretical error estimates and confirming results with numerical experiments.

Contribution

It introduces stochastic modified equations for Poisson integrators and rigorously analyzes their long-term accuracy and structure preservation.

Findings

Error estimates show long-term accuracy of integrators

Numerical experiments confirm theoretical predictions

Integrators effectively preserve stochastic Hamiltonian structure

Abstract

We address our attention to the numerical time discretization of stochastic Poisson systems via Poisson integrators. The aim of the investigation regards the backward error analysis of such integrators to reveal their ability of being structure-preserving, for long times of integration. In particular, we first provide stochastic modified equations suitable for such integrators and then we rigorously study them to prove accurate estimates on the long-term numerical error along the dynamics generated by stochastic Poisson integrators, with reference to the preservation of the random Hamiltonian conserved along the exact flow of the approximating Wong-Zakai Poisson system. Finally, selected numerical experiments confirm the effectiveness of the theoretical analysis.