We present single and binary black-hole simulations that follow the ''moving-puncture'' paradigm of simulating black-hole spacetimes without excision, and use ''moving boxes'' mesh refinement. Focusing on binary black-hole configurations where the simulations cover roughly two orbits, we address five major issues determining the quality of our results: numerical discretization error, finite extraction radius of the radiation signal, physical appropriateness of initial data, gauge choice, and computational performance. We also compare results we have obtained with the BAM code described here with the independent LEAN code.

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.

We study the inspiral, merger and ringdown of unequal mass black hole binaries by analyzing a catalogue of numerical simulations for seven different values of the mass ratio (from q = M2/M1 = 1 to q = 4). We compare numerical and Post-Newtonian results by projecting the waveforms onto spin-weighted spherical harmonics, characterized by angular indices (l , m). We find that the PostNewtonian equations predict remarkably well the relation between the wave amplitude and the orbital frequency for each (l , m), and that the convergence of the Post-Newtonian series to the numerical results is non-monotonic. To leading order the total energy emitted in the merger phase scales like η 2 and the spin of the final black hole scales like η, where η = q/(1 + q) 2 is the symmetric mass ratio. We study the multipolar distribution of the radiation, finding that odd-l multipoles are suppressed in the equal mass limit. Higher multipoles carry a larger fraction of the total energy as q increases. We introduce and compare three different definitions for the ringdown starting time. Applying linear estimation methods (the so-called Prony methods) to the ringdown phase, we find resolution-dependent time variations in the fitted parameters of the final black hole. By crosscorrelating information from different multipoles we show that ringdown fits can be used to obtain precise estimates of the mass and spin of the final black hole, which are in remarkable agreement with energy and angular momentum balance calculations.

Recent progress in numerical relativity has enabled us to model the nonperturbative merger phase of the binary black-hole coalescence problem. Based on these results, we propose a phenomenological family of waveforms which can model the inspiral, merger and ring-down stages of black-hole coalescence. We also construct a template bank using this family of waveforms and discuss its implementation in the search for signatures of gravitational waves produced by black-hole coalescences in the data of ground-based interferometers. This template bank might enable us to extend the present inspiral searches to highermass binary black-hole systems, i.e., systems with total mass greater than about 80 solar masses, thereby increasing the reach of the current generation of ground-based detectors.

Gravitational waveforms from the inspiral and ring-down stages of the binary black-hole coalescences can be modeled accurately by approximation/perturbation techniques in general relativity. Recent progress in numerical relativity has enabled us to model also the nonperturbative merger phase of the binary black-hole coalescence problem. This enables us to coherently search for all three stages of the coalescence of nonspinning binary black holes using a single template bank. Taking our motivation from these results, we propose a family of template waveforms which can model the inspiral, merger, and ringdown stages of the coalescence of nonspinning binary black holes that follow quasicircular inspiral. This two-dimensional template family is explicitly parametrized by the physical parameters of the binary. We show that the template family is not only effectual in detecting the signals from black-hole coalescences, but also faithful in estimating the parameters of the binary. We compare the sensitivity of a search (in the context of different ground-based interferometers) using all three stages of the black-hole coalescence with other template-based searches which look for individual stages separately. We find that the proposed search is significantly more sensitive than other template-based searches for a substantial mass range, potentially bringing about remarkable improvement in the event rate of ground-based interferometers. As part of this work, we also prescribe a general procedure to construct interpolated template banks using nonspinning black-hole waveforms produced by numerical relativity.