A Symmetric Integrator for non-integrable Hamiltonian Relativistic Systems

TL;DR

This paper introduces a symmetric, reversible integrator with a novel step size controller for non-integrable Hamiltonian systems, improving long-term accuracy and efficiency in simulating relativistic geodesic orbits.

Contribution

A new symmetric integrator with a custom step size controller that conserves constants of motion and outperforms existing methods in accuracy and speed for non-integrable Hamiltonian systems.

Findings

Faster than standard symplectic integrators

More accurate than adaptive Runge-Kutta schemes

Effective for long-term simulations of geodesic orbits

Abstract

By combining a standard symmetric, symplectic integrator with a new step size controller, we provide an integration scheme that is symmetric, reversible and conserves the values of the constants of motion. This new scheme is appropriate for long term numerical integrations of geodesic orbits in spacetime backgrounds, whose corresponding Hamiltonian system is non-integrable, and, in general, for any non-integrable Hamiltonian system whose kinetic part depends on the position variables. We show by numerical examples that the new integrator is faster and more accurate i) than the standard symplectic integration schemes with or without standard adaptive step size controllers and ii) than an adaptive step Runge-Kutta scheme.