Critical Phenomena in the Collapse of Gravitational Waves.

Fine-tuning generic but smooth spherically symmetric initial data for general relativity to the threshold of dynamical black hole formation creates arbitrarily large curvatures, mediated by a universal self-similar solution that acts as an intermediate attractor. For vacuum gravitational waves, however, these critical phenomena have been elusive. We present, for the first time, excellent agreement among three independent numerical simulations of this collapse. Surprisingly, we find no universality, and observe approximate self-similarity for some families of initial data but not for others.

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.

Geometry and regularity of moving punctures.

Significant advances in numerical simulations of black-hole binaries have recently been achieved using the puncture method. We examine how and why this method works by evolving a single black hole. The coordinate singularity and hence the geometry at the puncture are found to change during evolution, from representing an asymptotically flat end to being a cylinder. We construct an analytic solution for the stationary state of a black hole in spherical symmetry that matches the numerical result and demonstrates that the evolution is not dominated by artefacts at the puncture but indeed finds the analytical result.

Maximum kick from nonspinning black-hole binary inspiral.

When unequal-mass black holes merge, the final black hole receives a kick due to the asymmetric loss of linear momentum in the gravitational radiation emitted during the merger. The magnitude of this kick has important astrophysical consequences. Recent breakthroughs in numerical relativity allow us to perform the largest parameter study undertaken to date in numerical simulations of binary black-hole inspirals. We study nonspinning black-hole binaries with mass ratios from q=M1/M2=1 to q=0.25 (eta=q/(1+q)2 from 0.25 to 0.16). We accurately calculate the velocity of the kick to within 6%, and the final spin of the black holes to within 2%. A maximum kick of 175.2+/-11 km s(-1) is achieved for eta=0.195+/-0.005.

Supermassive recoil velocities for binary black-hole mergers with antialigned spins.

Recent calculations of the recoil velocity in binary black-hole mergers have found the kick velocity to be of the order of a few hundred km/s in the case of nonspinning binaries and about 500 km/s in the case of spinning configurations, and have lead to predictions of a maximum kick of up to 1300 km/s. We test these predictions and demonstrate that kick velocities of at least 2500 km/s are possible for equal-mass binaries with antialigned spins in the orbital plane. Kicks of that magnitude are likely to have significant repercussions for models of black-hole formation, the population of intergalactic black holes, and the structure of host galaxies.

Numerical simulation of orbiting black holes.

We present numerical simulations of binary black hole systems which for the first time last for about one orbital period for close but still separate black holes as indicated by the absence of a common apparent horizon. An important part of the method is the construction of comoving coordinates, in which both the angular and the radial motion are minimized through a dynamically adjusted shift condition. We use fixed mesh refinement for computational efficiency.