High order splitting methods for SDEs satisfying a commutativity condition

Abstract

In this paper, we introduce a new simple approach to developing and establishing the convergence of splitting methods for a large class of stochastic differential equations (SDEs), including additive, diagonal and scalar noise types. The central idea is to view the splitting method as a replacement of the driving signal of an SDE, namely Brownian motion and time, with a piecewise linear path that yields a sequence of ODEs $-$ which can be discretized to produce a numerical scheme. This new way of understanding splitting methods is inspired by, but does not use, rough path theory. We show that when the driving piecewise linear path matches certain iterated stochastic integrals of Brownian motion, then a high order splitting method can be obtained. We propose a general proof methodology for establishing the strong convergence of these approximations that is akin to the general framework of Milstein and Tretyakov. That is, once local error estimates are obtained for the splitting method, then a global rate of convergence follows. This approach can then be readily applied in future research on SDE splitting methods. By incorporating recently developed approximations for iterated integrals of Brownian motion into these piecewise linear paths, we propose several high order splitting methods for SDEs satisfying a certain commutativity condition. In our experiments, which include the Cox-Ingersoll-Ross model and additive noise SDEs (noisy anharmonic oscillator, stochastic FitzHugh-Nagumo model, underdamped Langevin dynamics), the new splitting methods exhibit convergence rates of $O(h^{3/2})$ and outperform schemes previously proposed in the literature.

arXiv comment49 pages, 9 figures

References in corpus (10)

• Accurate sampling using Langevin dynamics
• Strong convergence of an explicit numerical method for SDEs with nonglobally Lipschitz continuous coefficients
• Weak approximation of stochastic differential equations and application to derivative pricing
• A new weak approximation scheme of stochastic differential equations and the Runge-Kutta method
• On arbitrarily slow convergence rates for strong numerical approximations of Cox-Ingersoll-Ross processes and squared Bessel processes
• An optimal polynomial approximation of Brownian motion
• A randomized and fully discrete Galerkin finite element method for semilinear stochastic evolution equations
• Brownian bridge expansions for Lévy area approximations and particular values of the Riemann zeta function
• Optimal friction matrix for underdamped Langevin sampling
• An asymptotic radius of convergence for the Loewner equation and simulation of $SL{E}_k$ traces via splitting

Cited by in corpus (1)

• Generative Modelling of Lévy Area for High Order SDE Simulation