Symplectic Integration of Post-Newtonian Equations of Motion with Spin
arXiv:1003.5122 · doi:10.1103/PhysRevD.81.104025
Abstract
We present a non-canonically symplectic integration scheme tailored to numerically computing the post-Newtonian motion of a spinning black-hole binary. Using a splitting approach we combine the flows of orbital and spin contributions. In the context of the splitting, it is possible to integrate the individual terms of the spin-orbit and spin-spin Hamiltonians analytically, exploiting the special structure of the underlying equations of motion. The outcome is a symplectic, time-reversible integrator, which can be raised to arbitrary order by composition. A fourth-order version is shown to give excellent behavior concerning error growth and conservation of energy and angular momentum in long-term simulations. Favorable properties of the integrator are retained in the presence of weak dissipative forces due to radiation damping in the full post-Newtonian equations.
9 pages, 6 figures
References in corpus (3)
Cited by in corpus (14)
- 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
- A Symmetric Integrator for non-integrable Hamiltonian Relativistic Systems
- Applying explicit symplectic integrator to study chaos of charged particles around magnetized Kerr black hole
- A Note on the Construction of Explicit Symplectic Integrators for Schwarzschild Spacetimes
- Explicit symplectic methods in black hole spacetimes
- Dynamics of charged particles moving around Kerr black hole with inductive charge and external magnetic field
- Charged particle motions near non-Schwarzschild black holes with external magnetic fields in modified theories of gravity
- Energy-conserving integrator for conservative Hamiltonian systems with ten-dimensional phase space
- Adjustment of force-gradient operator in symplectic methods
- Free Versus Constrained Evolution of the 2+1 Equivariant Wave Map
- Contrasting the Implicit Method in Incoherent Lagrangian and the Correction Map Method in Hamiltonian
- Dissipated Correction Map Method with Trapezoidal Rule for the Simulations of Gravitational Waves from Spinning Compact Binary