Explicit methods in extended phase space for inseparable Hamiltonian problems

Abstract

We present a method for explicit leapfrog integration of inseparable Hamiltonian systems by means of an extended phase space. A suitably defined new Hamiltonian on the extended phase space leads to equations of motion that can be numerically integrated by standard symplectic leapfrog (splitting) methods. When the leapfrog is combined with coordinate mixing transformations, the resulting algorithm shows good long term stability and error behaviour. We extend the method to non-Hamiltonian problems as well, and investigate optimal methods of projecting the extended phase space back to original dimension. Finally, we apply the methods to a Hamiltonian problem of geodesics in a curved space, and a non-Hamiltonian problem of a forced non-linear oscillator. We compare the performance of the methods to a general purpose differential equation solver LSODE, and the implicit midpoint method, a symplectic one-step method. We find the extended phase space methods to compare favorably to both for the Hamiltonian problem, and to the implicit midpoint method in the case of the non-linear oscillator.

Accepted to Celestial Mechanics and Dynamical Astronomy. 22 pages, 4 figures, 3 tables

References in corpus (2)

• Implementing Few-Body Algorithmic Regularization with Post-Newtonian Terms
• Algorithmic regularization with velocity-dependent forces

Cited by in corpus (38)

• Explicit symplectic approximation of nonseparable Hamiltonians: algorithm and long time performance
• Post-Newtonian dynamical modeling of supermassive black holes in galactic-scale simulations
• The formation of extremely diffuse galaxy cores by merging supermassive black holes
• 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
• MSTAR -- a fast parallelised algorithmically regularised integrator with minimum spanning tree coordinates
• General purpose ray-tracing and polarized radiative transfer in General Relativity
• A Note on the Construction of Explicit Symplectic Integrators for Schwarzschild Spacetimes
• Gravitational Waves from the Inspiral of Supermassive Black Holes in Galactic-scale Simulations
• Explicit symplectic methods in black hole spacetimes
• KETJU -- resolving small-scale supermassive black hole dynamics in GADGET-4
• Extended phase-space symplectic-like integrators for coherent post-Newtonian Euler-Lagrange equations
• A novel energy-conserving scheme for eight-dimensional systems
• Introducing an Explicit Symplectic Integration Scheme for Riemannian Manifold Hamiltonian Monte Carlo
• The use of Kepler solver in numerical integrations of quasi-Keplerian orbits
• JANUS: A bit-wise reversible integrator for N-body dynamics
• Energy-conserving integrator for conservative Hamiltonian systems with ten-dimensional phase space
• A general model and toolkit for the ionization of three or more electrons in strongly driven atoms using an effective Coulomb potential for the interaction between bound electrons
• Semiexplicit Symplectic Integrators for Non-separable Hamiltonian Systems
• Explicit symplectic integrators with adaptive time steps in curved spacetimes
• FANTASY: User-Friendly Symplectic Geodesic Integrator for Arbitrary Metrics with Automatic Differentiation
• Preservation of Quadratic Invariants by Semiexplicit Symplectic Integrators for Non-separable Hamiltonian Systems
• Signatures of magnetic field effects in non-sequential double ionization manifesting as back-scattering for molecules versus forward-scattering for atoms
• Performance of different correction maps in the extended phase-space method for spinning compact binaries
• Adjustment of force-gradient operator in symplectic methods
• Explicit K-symplectic methods for nonseparable non-canonical Hamiltonian systems
• Modeling the Time Evolution of Compact Binary Systems with Machine Learning
• Numerical Methods and Closed Orbits in the Kepler-Heisenberg Problem
• Quantum versus semi-classical signatures of correlated triple ionization in Dalitz plots
• Adaptive Hamiltonian Variational Integrators and Symplectic Accelerated Optimization
• Generalized flow-composed symplectic methods for post-Newtonian Hamiltonian systems
• A new correction method for quasi-Keplerian orbits
• Contrasting the Implicit Method in Incoherent Lagrangian and the Correction Map Method in Hamiltonian
• Geometric numerical integration of Liénard systems via a contact Hamiltonian approach
• Explicit Symplectic Integrators for Massive Point Vortex Dynamics in Binary Mixture of Bose--Einstein Condensates
• Extended phase-space symplectic integration for electron dynamics
• Dissipated Correction Map Method with Trapezoidal Rule for the Simulations of Gravitational Waves from Spinning Compact Binary