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.

We examine the effects of spin-orbit and spin-spin coupling on the inspiral of a coalescing binary system of spinning compact objects and on the gravitational radiation emitted therefrom. Using a formalism developed by Blanchet, Damour, and Iyer, we calculate the contributions due to the spins of the bodies to the symmetric trace-free radiative multipole moments which are used to calculate the waveform, energy loss, and angular momentum loss from the inspiralling binary. Using equations of motion which include terms due to spin-orbit and spin-spin coupling, we evolve the orbit of a coalescing binary and use the orbit to calculate the emitted gravitational waveform. We find the spins of the bodies affect the waveform in several ways: 1) The spin terms contribute to the orbital decay of the binary, and thus to the accumulated * Present address: Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208. 1 phase of the gravitational waveform. 2) The spins cause the orbital plane to precess, which changes the orientation of the orbital plane with respect to an observer, thus causing the shape of the waveform to be modulated. 3) The spins contribute directly to the amplitude of the waveform. We discuss the size and importance of spin effects for the case of two coalescing neutron stars, and for the case of a neutron star orbiting a rapidly rotating 10M ⊙ black hole.

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.