Journal
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY
Volume 446, Issue 2, Pages 1424-1437Publisher
OXFORD UNIV PRESS
DOI: 10.1093/mnras/stu2164
Keywords
gravitation; methods: numerical; planets and satellites: dynamical evolution and stability
Categories
Funding
- NSF [AST-0807444]
- Keck Fellowship
- Association of Members of the Institute for Advanced
Ask authors/readers for more resources
We present IAS15, a 15th-order integrator to simulate gravitational dynamics. The integrator is based on a Gauss-Radau quadrature and can handle conservative as well as non-conservative forces. We develop a step-size control that can automatically choose an optimal timestep. The algorithm can handle close encounters and high-eccentricity orbits. The systematic errors are kept well below machine precision, and long-term orbit integrations over 109 orbits show that IAS15 is optimal in the sense that it follows Brouwer's law, i.e. the energy error behaves like a random walk. Our tests show that IAS15 is superior to a mixed-variable symplectic integrator and other popular integrators, including high-order ones, in both speed and accuracy. In fact, IAS15 preserves the symplecticity of Hamiltonian systems better than the commonly used nominally symplectic integrators to which we compared it. We provide an open-source implementation of IAS15. The package comes with several easy-to-extend examples involving resonant planetary systems, Kozai-Lidov cycles, close encounters, radiation pressure, quadrupole moment and generic damping functions that can, among other things, be used to simulate planet-disc interactions. Other non-conservative forces can be added easily.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available