Forward and non-forward symplectic integrators in solving classical dynamics problems

Abstract

Forward time step integrators are splitting algorithms with only positive splitting coefficients. When used in solving physical evolution equations, these positive coefficients correspond to positive time steps. Forward algorithms are essential for solving time-irreversible equations that cannot be evolved using backward time steps. However, forward integrators are also better in solving time-reversible equations of classical dynamics by tracking as closely as possible the physical trajectory. This work compares in detail various forward and non-forward fourth-order integrators using three, fourth, five and six force evaluations. In the case of solving the 2D Kepler orbit, all non-forward integrators are optimized by simply minimizing the size of their backward time steps.

26 pages, 17 figures, to be published in a "Special Issue on Splitting Methods for Differential Equation", Intl. J. of Comp. Math

References in corpus (10)

• Worm Algorithm for Continuous-space Path Integral Monte Carlo Simulations
• On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems
• Fourth-Order Algorithms for Solving the Imaginary Time Gross-Pitaevskii Equation in a Rotating Anisotropic Trap
• Quantum Statistical Calculations and Symplectic Corrector Algorithms
• Fourth Order Algorithms for Solving the Multivariable Langevin Equation and the Kramers Equation
• Higher Order Force Gradient Symplectic Algorithms
• Forward Symplectic Integrators and the Long Time Phase Error in Periodic Motions
• The physics of symplectic integrators: perihelion advances and symplectic corrector algorithms
• Short time evolved wave functions for solving quantum many-body problems
• A Fundamental Theorem on the Structure of Symplectic Integrators