Terrestrial Exoplanet Simulator (TES): an error optimal planetary systems integrator that permits close encounters

We present TES, a new n-body integration code for the accurate and rapid propagation of planetary systems in the presence of close encounters. TES builds upon the classic Encke method and integrates only the perturbations to Keplerian trajectories to reduce both the error and runtime of simulations. Variable step size is used throughout to enable close encounters to be precisely handled. A suite of numerical improvements are presented that together makes TES optimal in terms of energy error. Lower runtimes are found in all test problems considered when compared to direct integration using ias15. TES is freely available.