Free and constrained symplectic integrators for numerical general relativity

Abstract

We consider symplectic time integrators in numerical General Relativity and discuss both free and constrained evolution schemes. For free evolution of ADM-like equations we propose the use of the Stoermer-Verlet method, a standard symplectic integrator which here is explicit in the computationally expensive curvature terms. For the constrained evolution we give a formulation of the evolution equations that enforces the momentum constraints in a holonomically constrained Hamiltonian system and turns the Hamilton constraint function from a weak to a strong invariant of the system. This formulation permits the use of the constraint-preserving symplectic RATTLE integrator, a constrained version of the Stoermer-Verlet method. The behavior of the methods is illustrated on two effectively 1+1-dimensional versions of Einstein's equations, that allow to investigate a perturbed Minkowski problem and the Schwarzschild space-time. We compare symplectic and non-symplectic integrators for free evolution, showing very different numerical behavior for nearly-conserved quantities in the perturbed Minkowski problem. Further we compare free and constrained evolution, demonstrating in our examples that enforcing the momentum constraints can turn an unstable free evolution into a stable constrained evolution. This is demonstrated in the stabilization of a perturbed Minkowski problem with Dirac gauge, and in the suppression of the propagation of boundary instabilities into the interior of the domain in Schwarzschild space-time.

25 pages, 7 figures; This version contains minor clarifications and corrections

References in corpus (18)

• Evolution of Binary Black Hole Spacetimes
• Accurate Evolutions of Orbiting Black-Hole Binaries Without Excision
• Gravitational wave extraction from an inspiraling configuration of merging black holes
• Calibration of Moving Puncture Simulations
• A New Generalized Harmonic Evolution System
• Numerical Investigation of Cosmological Singularities
• Simulation of Binary Black Hole Spacetimes with a Harmonic Evolution Scheme
• A constrained scheme for Einstein equations based on Dirac gauge and spherical coordinates
• Unequal Mass Binary Black Hole Plunges and Gravitational Recoil
• Toward standard testbeds for numerical relativity
• Einstein's Equations with Asymptotically Stable Constraint Propagation
• New Algorithm for Mixmaster Dynamics
• Hyperbolicity of second-order in space systems of evolution equations
• The Midpoint Rule as a Variational--Symplectic Integrator. I. Hamiltonian Systems
• Hamiltonian Time Evolution for General Relativity
• Numerical stability for finite difference approximations of Einstein's equations
• Different canonical formulations of Einstein's theory of gravity
• On the applicability of constrained symplectic integrators in general relativity

Cited by in corpus (7)

• Continuum and Discrete Initial-Boundary-Value Problems and Einstein's Field Equations
• Cosmology using numerical relativity
• Hyperbolic formulations of General Relativity with Hamiltonian structure
• Numerical relativity simulations in the era of the Einstein Telescope
• Strongly hyperbolic Hamiltonian systems in numerical relativity: Formulation and symplectic integration
• Free Versus Constrained Evolution of the 2+1 Equivariant Wave Map
• Constructing of constraint preserving scheme for Einstein equations