Energy-conserving integrator for conservative Hamiltonian systems with ten-dimensional phase space

Abstract

In this paper, an implicit nonsymplectic exact energy-preserving integrator is specifically designed for a ten-dimensional phase-space conservative Hamiltonian system with five degrees of freedom. It is based on a suitable discretization-averaging of the Hamiltonian gradient, with a second-order accuracy to numerical solutions. A one-dimensional disordered discrete nonlinear Schrödinger equation and a post-Newtonian Hamiltonian system of spinning compact binaries are taken as our two examples. We demonstrate numerically that the proposed algorithm exhibits good long-term performance in the preservation of energy, if roundoff errors are neglected. This result is independent of time steps, initial orbital eccentricities, and regular and chaotic orbital dynamical behavior. In particular, the application of appropriately large time steps to the new algorithm is helpful in reducing time-consuming and roundoff errors. This new method, combined with fast Lyapunov indicators, is well suited related to chaos in the two example problems. It is found that chaos in the former system is mainly responsible for one of the parameters. In the latter problem, a combination of small initial separations and high initial eccentricities can easily induce chaos.

arXiv comment27 pages, 14 figures

References in corpus (44)

• Observation of Gravitational Waves from a Binary Black Hole Merger
• GW190425: Observation of a Compact Binary Coalescence with Total Mass $\sim 3.4M_{\odot }$
• GW190521: A Binary Black Hole Merger with a Total Mass of $150M_{\odot }$
• Effective one-body approach to general relativistic two-body dynamics
• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• On the determination of the last stable orbit for circular general relativistic binaries at the third post-Newtonian approximation
• WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
• Dimensional regularization of the gravitational interaction of point masses
• Transition from inspiral to plunge in precessing binaries of spinning black holes
• Gravitational waves from inspiralling compact binaries: Energy loss and waveform to second--post-Newtonian order
• Third post-Newtonian dynamics of compact binaries: Noetherian conserved quantities and equivalence between the harmonic-coordinate and ADM-Hamiltonian formalisms
• Dynamical invariants for general relativistic two-body systems at the third post-Newtonian approximation
• Equivalence between the ADM-Hamiltonian and the harmonic-coordinates approaches to the third post-Newtonian dynamics of compact binaries
• Gravity Waves, Chaos, and Spinning Compact Binaries
• Lyapunov indices with two nearby trajectories in a curved spacetime
• Post-Newtonian accurate parametric solution to the dynamics of spinning compact binaries in eccentric orbits: The leading order spin-orbit interaction
• Computation of Lyapunov Exponents in General Relativity
• Effective one body Hamiltonian of two spinning black-holes with next-to-next-to-leading order spin-orbit coupling
• Explicit methods in extended phase space for inseparable Hamiltonian problems
• Resurvey of order and chaos in spinning compact binaries
• The dynamics of precessing binary black holes using the post-Newtonian approximation
• Symplectic structure of post-Newtonian Hamiltonian for spinning compact binaries
• Construction of Explicit Symplectic Integrators in General Relativity. I. Schwarzschild Black Holes
• On post-Newtonian orbits and the Galactic-center stars
• Construction of explicit symplectic integrators in general relativity. II. Reissner-Nordstrom black holes
• Revisit on "Ruling out chaos in compact binary systems"
• Chaos and Order in Models of Black Hole Pairs
• The fate of chaotic binaries
• Ruling Out Chaos in Compact Binary Systems
• Symplectic Integration of Post-Newtonian Equations of Motion with Spin
• Characteristics of chaos evolution in one-dimensional disordered nonlinear lattices
• Comment on "Ruling out chaos in compact binary systems"
• Generalized, energy-conserving numerical simulations of particles in general relativity. I. Time-like and null geodesics
• Generalized, energy-conserving numerical simulations of particles in general relativity. II. Test particles in electromagnetic fields and GRMHD
• The deterministic nature of conservative post-Newtonian accurate dynamics of compact binaries with leading order spin-orbit interaction
• Chaos in two black holes with next-to-leading order spin-spin interactions
• Second post-Newtonian Lagrangian dynamics of spinning compact binaries
• Spin supplementary conditions for spinning compact binaries
• High order three part split symplectic integrators: Efficient techniques for the long time simulation of the disordered discrete nonlinear Schroedinger equation
• A novel energy-conserving scheme for eight-dimensional systems
• Are long-term $N$-body simulations reliable?
• Chaos may make black holes bright
• The use of Kepler solver in numerical integrations of quasi-Keplerian orbits
• On the accuracy of symplectic integrators for secularly evolving planetary systems

Cited by in corpus (5)

• Construction of explicit symplectic integrators in general relativity. IV. Kerr black holes
• Construction of explicit symplectic integrators in general relativity. III. Reissner-Nordstrom-(anti)-de Sitter black holes
• Impact of electric charges on chaos in magnetized $Reissner-Nordstr\ddot{o}m$ spacetimes
• Explicit K-symplectic methods for nonseparable non-canonical Hamiltonian systems
• Construction of second-order six-dimensional Hamiltonian-conserving scheme