We present the results from three gravitational-wave searches for coalescing compact binaries with component masses above $1M_{\odot }$ during the first and second observing runs of the advanced gravitational-wave detector network. During the first observing run ($O1$), from September 12, 2015 to January 19, 2016, gravitational waves from three binary black hole mergers were detected. The second observing run ($O2$), which ran from November 30, 2016 to August 25, 2017, saw the first detection of gravitational waves from a binary neutron star inspiral, in addition to the observation of gravitational waves from a total of seven binary black hole mergers, four of which we report here for the first time: GW170729, GW170809, GW170818, and GW170823. For all significant gravitational-wave events, we provide estimates of the source properties. The detected binary black holes have total masses between ${18.6}_{-0.7}^{+3.2}{M}_{\odot }$ and ${84.4}_{-11.1}^{+15.8}{M}_{\odot }$ and range in distance between ${320}_{-110}^{+120}$ and ${2840}_{-1360}^{+1400}\mathrm{Mpc}$. No neutron star–black hole mergers were detected. In addition to highly significant gravitational-wave events, we also provide a list of marginal event candidates with an estimated false-alarm rate less than 1 per 30 days. From these results over the first two observing runs, which include approximately one gravitational-wave detection per 15 days of data searched, we infer merger rates at the 90% confidence intervals of $110-3840{\mathrm{Gpc}}^{-3}{\mathrm{y}}^{-1}$ for binary neutron stars and $9.7-101{\mathrm{Gpc}}^{-3}{\mathrm{y}}^{-1}$ for binary black holes assuming fixed population distributions and determine a neutron star–black hole merger rate 90% upper limit of $610{\mathrm{Gpc}}^{-3}{\mathrm{y}}^{-1}$.

On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from to ( – if we restrict the dimensionless component spin magnitudes to be smaller than 0.05). These mass parameters are consistent with the individual binary components being neutron stars. However, both the source-frame chirp mass and the total mass of this system are significantly larger than those of any other known binary neutron star (BNS) system. The possibility that one or both binary components of the system are black holes cannot be ruled out from gravitational-wave data. We discuss possible origins of the system based on its inconsistency with the known Galactic BNS population. Under the assumption that the signal was produced by a BNS coalescence, the local rate of neutron star mergers is updated to 250–2810 .

On August 17, 2017, the Advanced LIGO and Advanced Virgo gravitational-wave detectors observed a low-mass compact binary inspiral. The initial sky localization of the source of the gravitational-wave signal, GW170817, allowed electromagnetic observatories to identify NGC 4993 as the host galaxy. In this work, we improve initial estimates of the binary’s properties, including component masses, spins, and tidal parameters, using the known source location, improved modeling, and recalibrated Virgo data. We extend the range of gravitational-wave frequencies considered down to 23 Hz, compared to 30 Hz in the initial analysis. We also compare results inferred using several signal models, which are more accurate and incorporate additional physical effects as compared to the initial analysis. We improve the localization of the gravitational-wave source to a 90% credible region of $16{\mathrm{deg}}^2$. We find tighter constraints on the masses, spins, and tidal parameters, and continue to find no evidence for nonzero component spins. The component masses are inferred to lie between 1.00 and $1.89{\mathrm{M}}_{\odot }$ when allowing for large component spins, and to lie between 1.16 and $1.60{\mathrm{M}}_{\odot }$ (with a total mass ${2.73}_{-0.01}^{+0.04}{\mathrm{M}}_{\odot }$) when the spins are restricted to be within the range observed in Galactic binary neutron stars. Using a precessing model and allowing for large component spins, we constrain the dimensionless spins of the components to be less than 0.50 for the primary and 0.61 for the secondary. Under minimal assumptions about the nature of the compact objects, our constraints for the tidal deformability parameter $\tilde{\mathrm{\Lambda }}$ are (0,630) when we allow for large component spins, and ${300}_{-230}^{+420}$ (using a 90% highest posterior density interval) when restricting the magnitude of the component spins, ruling out several equation-of-state models at the 90% credible level. Finally, with LIGO and GEO600 data, we use a Bayesian analysis to place upper limits on the amplitude and spectral energy density of a possible postmerger signal.

The detection of gravitational waves by Advanced LIGO and Advanced Virgo provides an opportunity to test general relativity in a regime that is inaccessible to traditional astronomical observations and laboratory tests. We present four tests of the consistency of the data with binary black hole gravitational waveforms predicted by general relativity. One test subtracts the best-fit waveform from the data and checks the consistency of the residual with detector noise. The second test checks the consistency of the low- and high-frequency parts of the observed signals. The third test checks that phenomenological deviations introduced in the waveform model (including in the post-Newtonian coefficients) are consistent with 0. The fourth test constrains modifications to the propagation of gravitational waves due to a modified dispersion relation, including that from a massive graviton. We present results both for individual events and also results obtained by combining together particularly strong events from the first and second observing runs of Advanced LIGO and Advanced Virgo, as collected in the catalog GWTC-1. We do not find any inconsistency of the data with the predictions of general relativity and improve our previously presented combined constraints by factors of 1.1 to 2.5. In particular, we bound the mass of the graviton to be $m_g\le 4.7\times {10}^{-23}\mathrm{eV}/c^2$ (90% credible level), an improvement of a factor of 1.6 over our previously presented results. Additionally, we check that the four gravitational-wave events published for the first time in GWTC-1 do not lead to stronger constraints on alternative polarizations than those published previously.

We report on the population of 47 compact binary mergers detected with a false-alarm rate of < 1 yr − 1 in the second LIGO–Virgo Gravitational-Wave Transient Catalog. We observe several characteristics of the merging binary black hole (BBH) population not discernible until now. First, the primary mass spectrum contains structure beyond a power law with a sharp high-mass cutoff; it is more consistent with a broken power law with a break at 39.7 − 9.1 + 20.3 M ⊙ or a power law with a Gaussian feature peaking at 33.1 − 5.6 + 4.0 M ⊙ (90% credible interval). While the primary mass distribution must extend to ∼ 65 M ⊙ or beyond, only 2.9 − 1.7 + 3.5 % of systems have primary masses greater than 45 M ⊙ . Second, we find that a fraction of BBH systems have component spins misaligned with the orbital angular momentum, giving rise to precession of the orbital plane. Moreover, 12 %– 44 % of BBH systems have spins tilted by more than 90°, giving rise to a negative effective inspiral spin parameter, χ eff . Under the assumption that such systems can only be formed by dynamical interactions, we infer that between 25% and 93% of BBHs with nonvanishing ∣ χ eff ∣ > 0.01 are dynamically assembled. Third, we estimate merger rates, finding R BBH = 23.9 − 8.6 + 14.3 Gpc − 3 yr − 1 for BBHs and R BNS = 320 − 240 + 490 Gpc − 3 yr − 1 for binary neutron stars. We find that the BBH rate likely increases with redshift ( 85 % credibility) but not faster than the star formation rate ( 86 % credibility). Additionally, we examine recent exceptional events in the context of our population models, finding that the asymmetric masses of GW190412 and the high component masses of GW190521 are consistent with our models, but the low secondary mass of GW190814 makes it an outlier.

We present results on the mass, spin, and redshift distributions with phenomenological population models using the 10 binary black hole (BBH) mergers detected in the first and second observing runs completed by Advanced LIGO and Advanced Virgo. We constrain properties of the BBH mass spectrum using models with a range of parameterizations of the BBH mass and spin distributions. We find that the mass distribution of the more massive BH in such binaries is well approximated by models with no more than 1% of BHs more massive than 45 M ⊙ and a power-law index of α = 1.3 − 1.7 + 1.4 (90% credibility). We also show that BBHs are unlikely to be composed of BHs with large spins aligned to the orbital angular momentum. Modeling the evolution of the BBH merger rate with redshift, we show that it is flat or increasing with redshift with 93% probability. Marginalizing over uncertainties in the BBH population, we find robust estimates of the BBH merger rate density of R = 53.2 − 28.2 + 55.8 Gpc⁻³ yr⁻¹ (90% credibility). As the BBH catalog grows in future observing runs, we expect that uncertainties in the population model parameters will shrink, potentially providing insights into the formation of BHs via supernovae, binary interactions of massive stars, stellar cluster dynamics, and the formation history of BHs across cosmic time.