High-order integrator for sampling the invariant distribution of a class of parabolic SPDEs with additive space-time noise

Abstract

We introduce a time-integrator to sample with high order of accuracy the invariant distribution for a class of semilinear SPDEs driven by an additive space-time noise. Combined with a postprocessor, the new method is a modification with negligible overhead of the standard linearized implicit Euler-Maruyama method. We first provide an analysis of the integrator when applied for SDEs (finite dimension), where we prove that the method has order $2$ for the approximation of the invariant distribution, instead of $1$. We then perform a stability analysis of the integrator in the semilinear SPDE context, and we prove in a linear case that a higher order of convergence is achieved. Numerical experiments, including the semilinear heat equation driven by space-time white noise, confirm the theoretical findings and illustrate the efficiency of the approach.

25 pages

Cited by in corpus (11)

• A full-discrete exponential Euler approximation of invariant measure for parabolic stochastic partial differential equations
• Optimal explicit stabilized integrator of weak order one for stiff and ergodic stochastic differential equations
• Exotic aromatic B-series for the study of long time integrators for a class of ergodic SDEs
• Numerical Unique Ergodicity of Monotone SDEs driven by Nondegenerate Multiplicative Noise
• Convergence analysis of explicit stabilized integrators for parabolic semilinear stochastic PDEs
• Dimension-free convergence rates for gradient Langevin dynamics in RKHS
• Numerical Ergodicity and Uniform Estimate of Monotone SPDEs Driven by Multiplicative Noise
• Finite-Volume approximation of the invariant measure of a viscous stochastic scalar conservation law
• Parareal exponential $\theta$-scheme for longtime simulation of stochastic Schrödinger equations with weak damping
• Numerical Analysis on Ergodic Limit of Approximations for Stochastic NLS Equation via Multi-symplectic Scheme
• Convergence in probability of an ergodic and conformal multi-symplectic numerical scheme for a damped stochastic NLS equation