Symmetric Partitioned Runge-Kutta Methods for Differential Equations on Lie-Groups
Abstract
Computer simulations in QCD are based on the discretization of the theory on a Euclidean lattice. To compute the mean value of an observable, usually the Hybrid Monte Carlo method is applied. Here equations of motion, derived from an Hamiltonian, have to be solved numerically. Commonly the Leapfrog (Stoermer-Verlet) method or splitting methods with multiple timescales à la Sexton-Weingarten are used to integrate the dynamical system, defined on a Lie group. Here we formulate time-reversible higher order integrators based on implicit partitioned Runge-Kutta schemes and show that they allow for larger step-sizes than the Leapfrog method. Since these methods are based on an infinite series of exponential functions, we concentrate on the truncation of this series with respect to the global error and accuracy. Finally, we see that the global error of a SPRK scheme is always even such that a convergence order of one is gained for methods with odd convergence order.
- Publication:
-
Proceedings of the XXIX International Symposium on Lattice Field Theory (Lattice 2011). July 10-16
- Pub Date:
- 2011
- DOI:
- 10.22323/1.139.0049
- arXiv:
- arXiv:1112.4336
- Bibcode:
- 2011slft.confE..49S
- Keywords:
-
- High Energy Physics - Lattice;
- High Energy Physics - Theory
- E-Print:
- 7 pages, 2 figures, talk given at the XXIX International Symposium on Lattice Field Theory, July 10-16, 2011, Squaw Valley, Lake Tahoe, CA