Continuum and Discrete Initial-Boundary Value Problems and Einstein's Field Equations

Self-Consistent Modeling of Gravitational Theories beyond General Relativity

The majority of extensions to general relativity (GR) display mathematical pathologies-higher derivatives, character change in equations that can be classified within partial differential equation theory, and even unclassifiable ones-that cause severe difficulties to study them, especially in dynamical regimes. We present here an approach that enables their consistent treatment and extraction of physical consequences. We illustrate this method in the context of single and merging black holes in a highly challenging beyond GR theory.

Intrinsic Pathology of Self-Interacting Vector Fields

We show that self-interacting vector field theories exhibit unphysical behavior even when they are not coupled to any external field. This means any theory featuring such vectors is in danger of being unphysical, an alarming prospect for many proposals in cosmology, gravity, high energy physics, and beyond. The problem arises when vector fields with healthy configurations naturally reach a point where time evolution is mathematically ill defined. We develop tools to easily identify this issue, and provide a simple and unifying framework to investigate it.

On ab initio-based, free and closed-form expressions for gravitational waves

We introduce a new approach for finding high accuracy, free and closed-form expressions for the gravitational waves emitted by binary black hole collisions from ab initio models. More precisely, our expressions are built from numerical surrogate models based on supercomputer simulations of the Einstein equations, which have been shown to be essentially indistinguishable from each other. Distinct aspects of our approach are that: (i) representations of the gravitational waves can be explicitly written in a few lines, (ii) these representations are free-form yet still fast to search for and validate and (iii) there are no underlying physical approximations in the underlying model. The key strategy is combining techniques from Artificial Intelligence and Reduced Order Modeling for parameterized systems. Namely, symbolic regression through genetic programming combined with sparse representations in parameter space and the time domain using Reduced Basis and the Empirical Interpolation Method enabling fast free-form symbolic searches and large-scale a posteriori validations. As a proof of concept we present our results for the collision of two black holes, initially without spin, and with an initial separation corresponding to 25-31 gravitational wave cycles before merger. The minimum overlap, compared to ground truth solutions, is 99%. That is, 1% difference between our closed-form expressions and supercomputer simulations; this is considered for gravitational (GW) science more than the minimum required due to experimental numerical errors which otherwise dominate. This paper aims to contribute to the field of GWs in particular and Artificial Intelligence in general.

Well-Posed Formulation of Scalar-Tensor Effective Field Theory

Effective field theory provides a way of parametrizing strong-field deviations from general relativity that might be observable in the gravitational waves emitted in a black hole merger. To perform numerical simulations of mergers in such theories it is necessary that the equations be written in a form that admits a well-posed initial value formulation. We study gravity coupled to a scalar field including the leading (four-derivative) effective field theory corrections. We introduce a new class of "modified harmonic" gauges and gauge-fixed equations of motion, such that, at weak coupling, the equations are strongly hyperbolic and therefore admit a well-posed initial value formulation.

Fundamentals of numerical relativity for gravitational wave sources

Einstein's theory of general relativity affords an enormously successful description of gravity. The theory encodes the gravitational interaction in the metric, a tensor field on spacetime that satisfies partial differential equations known as the Einstein equations. This review introduces some of the fundamental concepts of numerical relativity-solving the Einstein equations on the computer-in simple terms. As a primary example, we consider the solution of the general relativistic two-body problem, which features prominently in the new field of gravitational wave astronomy.

Initial boundary-value problem for the spherically symmetric Einstein equations with fluids with tangential pressure

We prove that, for a given spherically symmetric fluid distribution with tangential pressure on an initial space-like hypersurface with a time-like boundary, there exists a unique, local in time solution to the Einstein equations in a neighbourhood of the boundary. As an application, we consider a particular elastic fluid interior matched to a vacuum exterior.

Exploring New Physics Frontiers Through Numerical Relativity

The demand to obtain answers to highly complex problems within strong-field gravity has been met with significant progress in the numerical solution of Einstein's equations - along with some spectacular results - in various setups. We review techniques for solving Einstein's equations in generic spacetimes, focusing on fully nonlinear evolutions but also on how to benchmark those results with perturbative approaches. The results address problems in high-energy physics, holography, mathematical physics, fundamental physics, astrophysics and cosmology.

Universality, maximum radiation, and absorption in high-energy collisions of black holes with spin

We explore the impact of black hole spins on the dynamics of high-energy black hole collisions. We report results from numerical simulations with γ factors up to 2.49 and dimensionless spin parameter χ=+0.85, +0.6, 0, -0.6, -0.85. We find that the scattering threshold becomes independent of spin at large center-of-mass energies, confirming previous conjectures that structure does not matter in ultrarelativistic collisions. It has further been argued that in this limit all of the kinetic energy of the system may be radiated by fine tuning the impact parameter to threshold. On the contrary, we find that only about 60% of the kinetic energy is radiated for γ=2.49. By monitoring apparent horizons before and after scattering events we show that the "missing energy" is absorbed by the individual black holes in the encounter, and moreover the individual black-hole spins change significantly. We support this conclusion with perturbative calculations. An extrapolation of our results to the limit γ→∞ suggests that about half of the center-of-mass energy of the system can be emitted in gravitational radiation, while the rest must be converted into rest-mass and spin energy.

Characteristic Evolution and Matching

I review the development of numerical evolution codes for general relativity based upon the characteristic initial-value problem. Progress in characteristic evolution is traced from the early stage of 1D feasibility studies to 2D-axisymmetric codes that accurately simulate the oscillations and gravitational collapse of relativistic stars and to current 3D codes that provide pieces of a binary black-hole spacetime. Cauchy codes have now been successful at simulating all aspects of the binary black-hole problem inside an artificially constructed outer boundary. A prime application of characteristic evolution is to extend such simulations to null infinity where the waveform from the binary inspiral and merger can be unambiguously computed. This has now been accomplished by Cauchy-characteristic extraction, where data for the characteristic evolution is supplied by Cauchy data on an extraction worldtube inside the artificial outer boundary. The ultimate application of characteristic evolution is to eliminate the role of this outer boundary by constructing a global solution via Cauchy-characteristic matching. Progress in this direction is discussed.