The LIGO detection of GW150914 provides an unprecedented opportunity to study the two-body motion of a compact-object binary in the large-velocity, highly nonlinear regime, and to witness the final merger of the binary and the excitation of uniquely relativistic modes of the gravitational field. We carry out several investigations to determine whether GW150914 is consistent with a binary black-hole merger in general relativity. We find that the final remnant's mass and spin, as determined from the low-frequency (inspiral) and high-frequency (postinspiral) phases of the signal, are mutually consistent with the binary black-hole solution in general relativity. Furthermore, the data following the peak of GW150914 are consistent with the least-damped quasinormal mode inferred from the mass and spin of the remnant black hole. By using waveform models that allow for parametrized general-relativity violations during the inspiral and merger phases, we perform quantitative tests on the gravitational-wave phase in the dynamical regime and we determine the first empirical bounds on several high-order post-Newtonian coefficients. We constrain the graviton Compton wavelength, assuming that gravitons are dispersed in vacuum in the same way as particles with mass, obtaining a 90%-confidence lower bound of 10 13 km. In conclusion, within our statistical uncertainties, we find no evidence for violations of general relativity in the genuinely strong-field regime of gravity.

On September 14, 2015, the Laser Interferometer Gravitational-Wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in general relativity. GW150914 was produced by a nearly equal mass binary black hole of masses 36 þ5 −4 M ⊙ and 29 þ4 −4 M ⊙ ; for each parameter we report the median value and the range of the 90% credible interval. The dimensionless spin magnitude of the more massive black hole is bound to be < 0.7 (at 90% probability). The luminosity distance to the source is 410 −0.07 . This black hole is significantly more massive than any other inferred from electromagnetic observations in the stellar-mass regime.

We improve the accuracy of the effective-one-body (EOB) waveforms that were employed during the first observing run of Advanced LIGO for binaries of spinning, nonprecessing black holes by calibrating them to a set of 141 numerical-relativity (NR) waveforms. The NR simulations expand the domain of calibration towards larger mass ratios and spins, as compared to the previous EOBNR model. Merger-ringdown waveforms computed in black-hole perturbation theory for Kerr spins close to extremal provide additional inputs to the calibration. For the inspiral-plunge phase, we use a Markov-chain Monte Carlo algorithm to efficiently explore the calibration space. For the merger-ringdown phase, we fit the NR signals with phenomenological formulae. After extrapolation of the calibrated model to arbitrary mass ratios and spins, the (dominant-mode) EOBNR waveforms have faithfulness -at design Advanced-LIGO sensitivity -above 99% against all the NR waveforms, including 16 additional waveforms used for validation, when maximizing only on initial phase and time. This implies a negligible loss in event rate due to modeling for these binary configurations. We find that future NR simulations at mass ratios 4 and double spin 0.8 will be crucial to resolve discrepancies between different ways of extrapolating waveform models. We also find that some of the NR simulations that already exist in such region of parameter space are too short to constrain the low-frequency portion of the models. Finally, we build a reduced-order version of the EOBNR model to speed up waveform generation by orders of magnitude, thus enabling intensive data-analysis applications during the upcoming observation runs of Advanced LIGO.

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.

Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of numerical simulations for merging black holes. The catalog contains 2018 distinct configurations (a factor of 11 increase compared to the 2013 SXS catalog), including 1426 spin-precessing configurations, with mass ratios between 1 and 10, and spin magnitudes up to 0.998. The median length of a waveform in the catalog is 39 cycles of the dominant = m = 2 gravitational-wave mode, with the shortest waveform containing 7.0 cycles and the longest 351.3 cycles. We discuss improvements such as correcting for moving centers of mass and extended coverage of the parameter space. We also present a thorough analysis of numerical errors, finding typical truncation errors corresponding to a waveform mismatch of ∼ 10 −4 . The simulations provide remnant masses and spins with uncertainties of 0.03% and 0.1% (90 th percentile), about an order of magnitude better than analytical models for remnant properties. The full catalog is publicly available at https://www.black-holes.org/waveforms . black holes and of the surrounding spacetime [31,32]. Simulations have also been used for visualizations of curved spacetime [33][34][35][36][37][38][39][40], investigations of spin quantities [41], and the relaxation of spacetime to the Kerr solution following merger [42][43][44]. The motion of the black hole horizons and horizon curvature quantities have been used to explore eccentric dynamics [45][46][47][48], spin precession [49][50][51][52], and the first law of binary black hole mechanics [53][54][55][56][57]. These in turn have been compared to analytic post-Newtonian and self-force approximations (see also [58][59][60]), mapping out the bounds of validity of these approximations.A key application of BBH simulations is the accurate modeling of gravitational waves emitted by these systems during their late inspiral, merger, and final ringdown. Waveforms extracted from BBH simulations are essential for analyzing observed gravitational-wave signals from black hole binaries. Indeed, all BBH observations by LIGO and Virgo were analyzed using waveform families that rely on numerical relativity for their construction, most notably effective-one-body waveform models [61-65] and phenomenological waveform models [66][67][68]. Numerical simulations are also central in validating such waveform models [69][70][71][72][73][74][75][76], and were used to validate GW searches [77][78][79]. Waveforms from numerical relativity are also used directly in parameter estimation [80,81], to construct template banks [82], and to construct waveform families without intermediate analytical models, using methods such as reduced order modeling [83][84][85][86]. Today's simula...