A Technique for Studying Strong and Weak Local Errors of Splitting Stochastic Integrators

We present a technique, based on so-called word series, to write down in a systematic way expansions of the strong and weak local errors of splitting algorithms for the integration of Stratonovich stochastic differential equations. Those expansions immediately lead to the corresponding order conditions. Word series are similar to, but simpler than, the B-series used to analyze Runge-Kutta and other one-step integrators. The suggested approach makes it unnecessary to use the Baker-Campbell-Hausdorff formula. As an application, we compare two splitting algorithms recently considered by Leimkuhler and Matthews to integrate the Langevin equations. The word series method bears out clearly reasons for the advantages of one algorithm over the other.

[1]  B. Leimkuhler,et al.  Rational Construction of Stochastic Numerical Methods for Molecular Sampling , 2012, 1203.5428.

[2]  Fernando Casas,et al.  Splitting and composition methods in the numerical integration of differential equations , 2008, 0812.0377.

[3]  J. M. Sanz-Serna,et al.  Averaging and Computing Normal Forms with Word Series Algorithms , 2018 .

[4]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[5]  J. Gaines,et al.  The algebra of iterated stochastic integrals , 1994 .

[6]  Simon J. A. Malham,et al.  Algebraic structure of stochastic expansions and efficient simulation , 2011, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[7]  Jacques Laskar,et al.  New families of symplectic splitting methods for numerical integration in dynamical astronomy , 2012, 1208.0689.

[8]  Kuo-Tsai Ciren,et al.  INTEGRATION OF PATHS, GEOMETRIC INVARIANTS AND A GENERALIZED BAKER-HAUSDORFF FORMULA , 2016 .

[9]  Jesús María Sanz-Serna,et al.  Higher-Order Averaging, Formal Series and Numerical Integration I: B-series , 2010, Found. Comput. Math..

[10]  Jesús María Sanz-Serna,et al.  Word Series for Dynamical Systems and Their Numerical Integrators , 2015, Foundations of Computational Mathematics.

[11]  E. Hairer,et al.  Geometric Numerical Integration , 2022, Oberwolfach Reports.

[12]  Ander Murua,et al.  The Hopf Algebra of Rooted Trees, Free Lie Algebras, and Lie Series , 2006, Found. Comput. Math..

[13]  Ernst Hairer,et al.  On the Butcher group and general multi-value methods , 1974, Computing.

[14]  Free Lie algebras,et al.  Free Lie algebras , 2015 .

[15]  Fabrice Baudoin.  Diffusion Processes and Stochastic Calculus , 2014 .

[16]  Kevin Burrage,et al.  Order Conditions of Stochastic Runge-Kutta Methods by B-Series , 2000, SIAM J. Numer. Anal..

[17]  Houman Owhadi,et al.  Long-Run Accuracy of Variational Integrators in the Stochastic Context , 2007, SIAM J. Numer. Anal..

[18]  J. M. Sanz-Serna,et al.  Formal series and numerical integrators: some history and some new techniques , 2015, 1503.06976.

[19]  Jesús María Sanz-Serna,et al.  Higher-Order Averaging, Formal Series and Numerical Integration II: The Quasi-Periodic Case , 2012, Foundations of Computational Mathematics.

[20]  B. Leimkuhler,et al.  Robust and efficient configurational molecular sampling via Langevin dynamics. , 2013, The Journal of chemical physics.

[21]  J. M. Sanz-Serna,et al.  Order conditions for numerical integrators obtained by composing simpler integrators , 1999, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[22]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[23]  Shouchuan Zhang,et al.  On the Hopf Algebra of Rooted Trees , 2007, 0710.5646.

[24]  B. Leimkuhler,et al.  The computation of averages from equilibrium and nonequilibrium Langevin molecular dynamics , 2013, 1308.5814.

[26]  M. Fliess,et al.  Fonctionnelles causales non linaires et indtermines non commutatives , 1981 .

[27]  R. Ree,et al.  Lie Elements and an Algebra Associated With Shuffles , 1958 .

[28]  G. Milstein.  Numerical Integration of Stochastic Differential Equations , 1994 .

[29]  J. M. Sanz-Serna,et al.  Canonical B-series , 1994 .

[30]  M. V. Tretyakov,et al.  Stochastic Numerics for Mathematical Physics , 2004, Scientific Computation.

[31]  J. M. Sanz-Serna,et al.  A formal series approach to averaging: Exponentially small error estimates , 2012 .

[32]  J. Butcher.  Coefficients for the study of Runge-Kutta integration processes , 1963, Journal of the Australian Mathematical Society.

[33]  Jesús María Sanz-Serna,et al.  Higher-Order Averaging, Formal Series and Numerical Integration III: Error Bounds , 2015, Found. Comput. Math..

[34]  J. M. Sanz-Serna,et al.  Symplectic Methods Based on Decompositions , 1997 .

[35]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .