---

High weak order methods for stochastic differential equations based on modified equations

---

A. Abdulle, D. Cohen, G. Vilmart and K.C. Zygalakis

---

Abstract. Inspired by recent advances in the theory of modified differential equations, we propose a new methodology for constructing numerical integrators with high weak order for the time integration of stochastic differential equations. This approach is illustrated with the constructions of new methods of weak order two, in particular, semi-implicit integrators well suited for stiff (mean-square stable) stochastic problems, and implicit integrators that exactly conserve all quadratic first integrals of a stochastic dynamical system. Numerical examples confirm the theoretical results and show the versatility of our methodology.

---

Key Words. weak convergence, modified equations, backward error analysis, stiff integrator, invariant preserving integrator

---