A novel energy-conserving scheme for eight-dimensional systems

Abstract

We design a novel, exactly energy-conserving implicit non-symplectic integration method for an eight-dimensional Hamiltonian system with four degrees of freedom. In our algorithm, each partial derivative of the Hamiltonian with respect to one of phase-space variables is discretized by the average of eight Hamiltonian difference terms. Such a discretization form is a second-order approximation to the Hamiltonian gradient. It is shown numerically via simulations of an FPU-$\beta$ system and a post-Newtonian conservative system of compact binaries with one body spinning that the newly proposed method has extremely good energy-conserving performance, compared to the Runge-Kutta, implicit midpoint symplectic method and extended phase-space explicit symplectic-like integrators. The new method is advantageous over very long times and for large time steps compared to state-of-the-art Runge-Kutta method in the accuracy of numerical solutions. Although such an energy-conserving integrator exhibits a higher computational cost than any one of the other three algorithms, the superior results justify its use for satisfying some specific purposes on the preservation of energies in numerical simulations with much longer times, e.g. obtaining a high enough accuracy of the semimajor axis in a Keplerian problem in the Solar System or accurately grasping the frequency of a gravitational wave from a circular orbit in a post-Newtonian system of compact binaries. The new integrator will be potentially applied to model time-varying external electromagnetic fields or time-dependent spacetimes.

arXiv comment17 pages, 6 figures

References in corpus (18)

• Coalescing binary systems of compact objects to (post)$^{5/2}-Newtonian order. V. Spin Effects
• On the determination of the last stable orbit for circular general relativistic binaries at the third post-Newtonian approximation
• Dimensional regularization of the gravitational interaction of point masses
• Transition from inspiral to plunge in precessing binaries of spinning black holes
• Dynamical invariants for general relativistic two-body systems at the third post-Newtonian approximation
• Lyapunov indices with two nearby trajectories in a curved spacetime
• What does a binary black hole merger look like?
• 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
• Detecting chaos, determining the dimensions of tori and predicting slow diffusion in Fermi--Pasta--Ulam lattices by the Generalized Alignment Index method
• Symplectic structure of post-Newtonian Hamiltonian for spinning compact binaries
• Inducing chaos by breaking axial symmetry in a black hole magnetosphere
• 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
• Efficient integration of the variational equations of multi-dimensional Hamiltonian systems: Application to the Fermi-Pasta-Ulam lattice
• Spin supplementary conditions for spinning compact binaries
• An fully implicit scheme for numerical integration of the relativistic particle equation of motion

Cited by in corpus (11)

• Construction of Explicit Symplectic Integrators in General Relativity. I. Schwarzschild Black Holes
• Construction of explicit symplectic integrators in general relativity. II. Reissner-Nordstrom black holes
• 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
• The use of Kepler solver in numerical integrations of quasi-Keplerian orbits
• Energy-conserving integrator for conservative Hamiltonian systems with ten-dimensional phase space
• Performance of different correction maps in the extended phase-space method for spinning compact binaries
• Explicit K-symplectic methods for nonseparable non-canonical Hamiltonian systems
• Construction of second-order six-dimensional Hamiltonian-conserving scheme
• Dissipated Correction Map Method with Trapezoidal Rule for the Simulations of Gravitational Waves from Spinning Compact Binary
• A new correction method for quasi-Keplerian orbits