High precision symplectic integrators for the Solar System

Abstract

Using a Newtonian model of the Solar System with all 8 planets, we perform extensive tests on various symplectic integrators of high orders, searching for the best splitting scheme for long term studies in the Solar System. These comparisons are made in Jacobi and heliocentric coordinates and the implementation of the algorithms is fully detailed for practical use. We conclude that high order integrators should be privileged, with a preference for the new (10, 6, 4) method of Blanes et al. (2013). © 2013 Springer Science+Business Media Dordrecht.This work was supported by GTSNext project. The work of SB, FC, JM and AM has been partially supported by Ministerio de Ciencia e Innovacion (Spain) under project MTM2010-18246-C03 (co-financed by FEDER Funds of the European Union).Farrés, A.; Laskar, J.; Blanes Zamora, S.; Casas Perez, F.; Makazaga, J.; Murua, A. (2013). High precision symplectic integrators for the Solar System. Celestial Mechanics and Dynamical Astronomy. 116:141-174. doi:10.1007/s10569-013-9479-6S14117411