The main goal of numerical relativity is the long time simulation of highly nonlinear spacetimes that cannot be treated by perturbation theory. This involves analytic, computational and physical issues. At present, the major impasses to achieving global simulations of physical usefulness are of an analytic/computational nature. We present here some examples of how analytic insight can lend useful guidance for the improvement of numerical approaches.
We present a number of open problems within general relativity. After a brief introduction to some technical mathematical issues and the famous singularity theorems, we discuss the cosmic censorship hypothesis and the Penrose inequality, the uniqueness of black hole solutions and the stability of Kerr spacetime and the final state conjecture, critical phenomena and the Einstein-Yang--Mills equations, and a number of other problems in classical general relativity. We then broaden the scope and discuss some mathematical problems motivated by quantum gravity, including AdS/CFT correspondence and problems in higher dimensions and, in particular, the instability of anti-de Sitter spacetime, and in cosmology, including the cosmological constant problem and dark energy, the stability of de Sitter spacetime and cosmological singularities and spikes. Finally, we briefly discuss some problems in numerical relativity and relativistic astrophysics.
We produce the first astrophysically-relevant numerical binary black hole gravitational waveform in a higher-curvature theory of gravity beyond general relativity. We simulate a system with parameters consistent with GW150914, the first LIGO detection, in order-reduced dynamical Chern-Simons gravity, a theory with motivations in string theory and loop quantum gravity. We present results for the leading-order corrections to the merger and ringdown waveforms, as well as the ringdown quasi-normal mode spectrum. We estimate that such corrections may be discriminated in detections with signal to noise ratio $gtrsim 180-240$, with the precise value depending on the dimension of the GR waveform family used in data analysis.
A numerical-relativity calculation yields in general a solution of the Einstein equations including also a radiative part, which is in practice computed in a region of finite extent. Since gravitational radiation is properly defined only at null infinity and in an appropriate coordinate system, the accurate estimation of the emitted gravitational waves represents an old and non-trivial problem in numerical relativity. A number of methods have been developed over the years to extract the radiative part of the solution from a numerical simulation and these include: quadrupole formulas, gauge-invariant metric perturbations, Weyl scalars, and characteristic extraction. We review and discuss each method, in terms of both its theoretical background as well as its implementation. Finally, we provide a brief comparison of the various methods in terms of their inherent advantages and disadvantages.
This document proposes data formats to exchange numerical relativity results, in particular gravitational waveforms. The primary goal is to further the interaction between gravitational-wave source modeling groups and the gravitational-wave data-analysis community. We present a simple and extendable format which is applicable to various kinds of gravitational wave sources including binaries of compact objects and systems undergoing gravitational collapse, but is nevertheless sufficiently general to be useful for other purposes.
We discuss a general formalism for numerically evolving initial data in general relativity in which the (complex) Ashtekar connection and the Newman-Penrose scalars are taken as the dynamical variables. In the generic case three gauge constraints and twelve reality conditions must be solved. The analysis is applied to a Petrov type {1111} planar spacetime where we find a spatially constant volume element to be an appropriate coordinate gauge choice.