Extending the lifetime of 3D black hole computations with a new hyperbolic system of evolution equations

The evolution of black-hole binaries in vacuum spacetimes constitutes the two-body problem in general relativity. The solution of this problem in the framework of the Einstein field equations is a substantially more complex exercise than that of the dynamics of two point masses in Newtonian gravity, but it also presents us with a wealth of new exciting physics. Numerical methods are likely the only method to compute the dynamics of black-hole systems in the fully non-linear r...

The 3+1 Hamiltonian formulation in the gauge $D_tN=-K$ on the lapse function fixes the direction of time associated with the trace $K$ of the extrinsic curvature tensor. The Hamiltonian equations hereby become hyperbolic. We study this new system for black hole spacetimes that are asymptotically quiescent, which introduces analyticity properties that can be exploited for numerical calculations by compactification in spherical coordinates with complex radius following a M\"obi...

This is the first paper in a series aimed to implement boundary conditions consistent with the constraints' propagation in 3D numerical relativity. Here we consider spherically symmetric black hole spacetimes in vacuum or with a minimally coupled scalar field, within the Einstein-Christoffel symmetric hyperbolic formulation of Einstein's equations. By exploiting the characteristic propagation of the main variables and constraints, we are able to single out the only free modes...

This article begins with a brief introduction to numerical relativity aimed at readers who have a background in applied mathematics but not necessarily in general relativity. I then introduce and summarise my work on the problem of treating asymptotically flat spacetimes of infinite extent with finite computational resources. Two different approaches are considered. The first approach is the standard one and is based on evolution on Cauchy hypersurfaces with artificial timeli...

We present a new fully first order strongly hyperbolic representation of the BSSN formulation of Einstein's equations with optional constraint damping terms. We describe the characteristic fields of the system, discuss its hyperbolicity properties, and present two numerical implementations and simulations: one using finite differences, adaptive mesh refinement and in particular binary black holes, and another one using the discontinuous Galerkin method in spherical symmetry. ...

I describe approaches to the study of black hole spacetimes via numerical relativity. After a brief review of the basic formalisms and techniques used in numerical black hole simulations, I discuss a series of calculations from axisymmetry to full 3D that can be seen as stepping stones to simulations of the full 3D coalescence of two black holes. In particular, I emphasize the interplay between perturbation theory and numerical simulation that build both confidence in present...

Since the breakthroughs in 2005 which have led to long term stable solutions of the binary black hole problem in numerical relativity, much progress has been made. I present here a short summary of the state of the field, including the capabilities of numerical relativity codes, recent physical results obtained from simulations, and improvements to the methods used to evolve and analyse binary black hole spacetimes.

It is often the case in numerical relativity that schemes that are known to be convergent for well posed systems are used in evolutions of weakly hyperbolic (WH) formulations of Einstein's equations. Here we explicitly show that with several of the discretizations that have been used through out the years, this procedure leads to non-convergent schemes. That is, arbitrarily small initial errors are amplified without bound when resolution is increased, independently of the amo...

We consider unconstrained evolution schemes for the hyperboloidal initial value problem in numerical relativity as a promising candidate for the optimally efficient numerical treatment of radiating compact objects. Here, spherical symmetry already poses nontrivial problems and constitutes an important first step to regularize the resulting singular PDEs. We evolve the Einstein equations in their generalized BSSN and Z4 formulations coupled to a massless self-gravitating scala...

We present a new evolution system for Einstein's field equations which is based on tetrad fields and conformally compactified hyperboloidal spatial hypersurfaces which reach future null infinity. The boost freedom in the choice of the tetrad is fixed by requiring that its timelike leg be orthogonal to the foliation, which consists of constant mean curvature slices. The rotational freedom in the tetrad is fixed by the 3D Nester gauge. With these conditions, the field equations...