Slimplectic Integrators: Variational Integrators for General Nonconservative Systems

Abstract

Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative interactions. In this Letter, we develop the "slimplectic" integrator, a new type of numerical integrator that shares many of the benefits of traditional symplectic integrators yet is applicable to general nonconservative systems. We utilize a fixed time-step variational integrator formalism applied to the principle of stationary nonconservative action developed in Galley, 2013; Galley, Tsang & Stein, 2014. As a result, the generalized momenta and energy (Noether current) evolutions are well-tracked. We discuss several example systems, including damped harmonic oscillators, Poynting-Robertson drag, and gravitational radiation reaction, by utilizing our new publicly available code to demonstrate the slimplectic integrator algorithm. Slimplectic integrators are well-suited for integrations of systems where nonconservative effects play an important role in the long-term dynamical evolution. As such they are particularly appropriate for cosmological or celestial N-body dynamics problems where nonconservative interactions, e.g. gas interactions or dissipative tides, can play an important role.

6 pages, 5 Figures; Accepted to ApJL; code repository at http://github.com/davtsang/slimplectic (typo fixed in repository URL)

References in corpus (1)

• Variational Integrators for Almost-Integrable Systems

Cited by in corpus (16)

• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• 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
• Integrability of eccentric, spinning black hole binaries up to second post-Newtonian order
• Explicit symplectic methods in black hole spacetimes
• Variational order for forced Lagrangian systems
• Numerical integration in celestial mechanics: a case for contact geometry
• Explicit symplectic integrators with adaptive time steps in curved spacetimes
• EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour
• 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
• A new variational discretization technique for initial value problems bypassing governing equations
• Dissipated Correction Map Method with Trapezoidal Rule for the Simulations of Gravitational Waves from Spinning Compact Binary
• Contrasting the Implicit Method in Incoherent Lagrangian and the Correction Map Method in Hamiltonian