Gravitational waves emitted by black-hole binary systems have the highest signal-to-noise ratio in LIGO and Virgo detectors when black-hole spins are aligned with the orbital angular momentum and extremal. For such systems, we extend the effective-one-body inspiral-merger-ringdown waveforms to generic mass ratios and spins calibrating them to 38 numerical-relativity nonprecessing waveforms produced by the SXS Collaboration. The numerical-relativity simulations span mass ratios from 1 to 8, spin magnitudes up to 98% of extremality, and last for 40 to 60 gravitational-wave cycles. When the total mass of the binary is between 20 and 200M ⊙ , the effective-one-body nonprecessing (dominant mode) waveforms have overlap above 99% (using the advanced-LIGO design noise spectral density) with all of the 38 nonprecessing numerical waveforms, when maximizing only on initial phase and time. This implies a negligible loss in event rate due to modeling. We also show that-without further calibration-the precessing effective-onebody (dominant mode) waveforms have overlap above 97% with two very long, strongly precessing numerical-relativity waveforms, when maximizing only on the initial phase and time.

Numerical simulations of 15 orbits of an equal-mass binary black-hole system are presented. Gravitational waveforms from these simulations, covering more than 30 cycles and ending about 1.5 cycles before merger, are compared with those from quasicircular zero-spin post-Newtonian (PN) formulae. The cumulative phase uncertainty of these comparisons is about 0.05 radians, dominated by effects arising from the small residual spins of the black holes and the small residual orbital eccentricity in the simulations. Matching numerical results to PN waveforms early in the run yields excellent agreement (within 0.05 radians) over the first 15 cycles, thus validating the numerical simulation and establishing a regime where PN theory is accurate. In the last 15 cycles to merger, however, generic time-domain Taylor approximants build up phase differences of several radians. But, apparently by coincidence, one specific post-Newtonian approximant, TaylorT4 at 3.5PN order, agrees much better with the numerical simulations, with accumulated phase differences of less than 0.05 radians over the 30-cycle waveform. Gravitational-wave amplitude comparisons are also done between numerical simulations and postNewtonian, and the agreement depends on the post-Newtonian order of the amplitude expansion: the amplitude difference is about 6%-7% for zeroth order and becomes smaller for increasing order. A newly derived 3.0PN amplitude correction improves agreement significantly ( < 1% amplitude difference throughout most of the run, increasing to 4% near merger) over the previously known 2.5PN amplitude terms.

We describe a general procedure to generate spinning, precessing waveforms that include inspiral, merger, and ringdown stages in the effective-one-body (EOB) approach. The procedure uses a precessing frame in which precession-induced amplitude and phase modulations are minimized, and an inertial frame, aligned with the spin of the final black hole, in which we carry out the matching of the inspiral-plunge to merger-ringdown waveforms. As a first application, we build spinning, precessing EOB waveforms for the gravitational modes l ¼ 2 such that in the nonprecessing limit those waveforms agree with the EOB waveforms recently calibrated to numerical-relativity waveforms. Without recalibrating the EOB model, we then compare EOB and post-Newtonian precessing waveforms to two numerical-relativity waveforms produced by the Caltech-Cornell-CITA collaboration. The numerical waveforms are strongly precessing and have 35 and 65 gravitational-wave cycles. We find a remarkable agreement between EOB and numerical-relativity precessing waveforms and spins' evolutions. The phase difference is ∼0.2 rad rad at merger, while the mismatches, computed using the advanced-LIGO noise spectral density, are below 2% when maximizing only on the time and phase at coalescence and on the polarization angle.