Initial Data for Numerical Relativity

Hamiltonian formulation of general relativity and post-Newtonian dynamics of compact binaries

Hamiltonian formalisms provide powerful tools for the computation of approximate analytic solutions of the Einstein field equations. The post-Newtonian computations of the explicit analytic dynamics and motion of compact binaries are discussed within the most often applied Arnowitt-Deser-Misner formalism. The obtention of autonomous Hamiltonians is achieved by the transition to Routhians. Order reduction of higher derivative Hamiltonians results in standard Hamiltonians. Tetrad representation of general relativity is introduced for the tackling of compact binaries with spinning components. Configurations are treated where the absolute values of the spin vectors can be considered constant. Compact objects are modeled by use of Dirac delta functions and their derivatives. Consistency is achieved through transition to d-dimensional space and application of dimensional regularization. At the fourth post-Newtonian level, tail contributions to the binding energy show up. The conservative spin-dependent dynamics finds explicit presentation in Hamiltonian form through next-to-next-to-leading-order spin-orbit and spin1-spin2 couplings and to leading-order in the cubic and quartic in spin interactions. The radiation reaction dynamics is presented explicitly through the third-and-half post-Newtonian order for spinless objects, and, for spinning bodies, to leading-order in the spin-orbit and spin1-spin2 couplings. The most important historical issues get pointed out.

Dynamical boson stars

The idea of stable, localized bundles of energy has strong appeal as a model for particles. In the 1950s, John Wheeler envisioned such bundles as smooth configurations of electromagnetic energy that he called geons, but none were found. Instead, particle-like solutions were found in the late 1960s with the addition of a scalar field, and these were given the name boson stars. Since then, boson stars find use in a wide variety of models as sources of dark matter, as black hole mimickers, in simple models of binary systems, and as a tool in finding black holes in higher dimensions with only a single Killing vector. We discuss important varieties of boson stars, their dynamic properties, and some of their uses, concentrating on recent efforts.

The initial value problem as it relates to numerical relativity

Spacetime is foliated by spatial hypersurfaces in the 3+1 split of general relativity. The initial value problem then consists of specifying initial data for all fields on one such a spatial hypersurface, such that the subsequent evolution forward in time is fully determined. On each hypersurface the 3-metric and extrinsic curvature describe the geometry. Together with matter fields such as fluid velocity, energy density and rest mass density, the 3-metric and extrinsic curvature then constitute the initial data. There is a lot of freedom in choosing such initial data. This freedom corresponds to the physical state of the system at the initial time. At the same time the initial data have to satisfy the Hamiltonian and momentum constraint equations of general relativity and can thus not be chosen completely freely. We discuss the conformal transverse traceless and conformal thin sandwich decompositions that are commonly used in the construction of constraint satisfying initial data. These decompositions allow us to specify certain free data that describe the physical nature of the system. The remaining metric fields are then determined by solving elliptic equations derived from the constraint equations. We describe initial data for single black holes and single neutron stars, and how we can use conformal decompositions to construct initial data for binaries made up of black holes or neutron stars. Orbiting binaries will emit gravitational radiation and thus lose energy. Since the emitted radiation tends to circularize the orbits over time, one can thus expect that the objects in a typical binary move on almost circular orbits with slowly shrinking radii. This leads us to the concept of quasi-equilibrium, which essentially assumes that time derivatives are negligible in corotating coordinates for binaries on almost circular orbits. We review how quasi-equilibrium assumptions can be used to make physically well motivated approximations that simplify the elliptic equations we have to solve.

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.

Spectral Methods for Numerical Relativity

Equations arising in general relativity are usually too complicated to be solved analytically and one must rely on numerical methods to solve sets of coupled partial differential equations. Among the possible choices, this paper focuses on a class called spectral methods in which, typically, the various functions are expanded in sets of orthogonal polynomials or functions. First, a theoretical introduction of spectral expansion is given with a particular emphasis on the fast convergence of the spectral approximation. We then present different approaches to solving partial differential equations, first limiting ourselves to the one-dimensional case, with one or more domains. Generalization to more dimensions is then discussed. In particular, the case of time evolutions is carefully studied and the stability of such evolutions investigated. We then present results obtained by various groups in the field of general relativity by means of spectral methods. Work, which does not involve explicit time-evolutions, is discussed, going from rapidly-rotating strange stars to the computation of black-hole-binary initial data. Finally, the evolution of various systems of astrophysical interest are presented, from supernovae core collapse to black-hole-binary mergers.

Recoil velocities from equal-mass binary-black-hole mergers

The final evolution of a binary-black-hole system gives rise to a recoil velocity if an asymmetry is present in the emitted gravitational radiation. Measurements of this effect for nonspinning binaries with unequal masses have pointed out that kick velocities approximately 175 km/s can be reached for a mass ratio approximately 0.36. However, a larger recoil can be obtained for equal-mass binaries if the asymmetry is provided by the spins. Using two independent methods we show that the merger of such binaries yields velocities as large as approximately 440 km/s for black holes having unequal spins that are antialigned and parallel to the orbital angular momentum.

Evolution of binary black-hole spacetimes

We describe early success in the evolution of binary black-hole spacetimes with a numerical code based on a generalization of harmonic coordinates. Indications are that with sufficient resolution this scheme is capable of evolving binary systems for enough time to extract information about the orbit, merger, and gravitational waves emitted during the event. As an example we show results from the evolution of a binary composed of two equal mass, nonspinning black holes, through a single plunge orbit, merger, and ringdown. The resultant black hole is estimated to be a Kerr black hole with angular momentum parameter a approximately 0.70. At present, lack of resolution far from the binary prevents an accurate estimate of the energy emitted, though a rough calculation suggests on the order of 5% of the initial rest mass of the system is radiated as gravitational waves during the final orbit and ringdown.

Metric of a tidally distorted nonrotating black hole

The metric of a tidally distorted, nonrotating black hole is presented in a light-cone coordinate system that penetrates the event horizon and possesses a clear geometrical meaning. The metric is expressed as an expansion in powers of r/R<<1, where r is a measure of distance from the black hole and R is the local radius of curvature of the external spacetime; this is assumed to be much larger than M, the mass of the black hole. The metric is calculated up to a remainder of order (r/R)4, and it depends on a family of tidal gravitational fields which characterize the hole's local environment. The coordinate system allows an easy identification of the event horizon, and expressions are derived for its surface gravity and the rates at which the tidal interaction transfers mass and angular momentum to the black hole.