Efficient high-order explicit symplectic splitting methods for post-Newtonian Hamiltonian systems

Abstract

The nonseparability of post-Newtonian (PN) Hamiltonian systems typically necessitates the use of computationally expensive implicit integrators. Recent research overcomes this limitation by embedding the dynamics into a doubled phase space, which enables the development of explicit symplectic methods. However, existing specially designed explicit integrators suffer from order reduction for high-order methods when the time stepsize is small, i.e., h <3. In this paper, we propose a novel extension and splitting approach for the doubled Hamiltonian, under which specially designed explicit symplectic integrators can be constructed. It is shown that the proposed integrators achieve genuine high-order convergence without order reduction and take advantage of the small PN parameter . Numerical results from simulations with 2PN spinning binaries demonstrate superior long-term conservation of invariants and significantly higher computational efficiency compared to both implicit methods and existing explicit splitting techniques.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…