We report the observation of a compact binary coalescence involving a 22.2–24.3 M _⊙ black hole and a compact object with a mass of 2.50–2.67 M _⊙ (all measurements quoted at the 90% credible level). The gravitational-wave signal, GW190814, was observed during LIGO’s and Virgo’s third observing run on 2019 August 14 at 21:10:39 UTC and has a signal-to-noise ratio of 25 in the three-detector network. The source was localized to 18.5 deg² at a distance of Mpc; no electromagnetic counterpart has been confirmed to date. The source has the most unequal mass ratio yet measured with gravitational waves, , and its secondary component is either the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system. The dimensionless spin of the primary black hole is tightly constrained to ≤0.07. Tests of general relativity reveal no measurable deviations from the theory, and its prediction of higher-multipole emission is confirmed at high confidence. We estimate a merger rate density of 1–23 Gpc⁻³ yr⁻¹ for the new class of binary coalescence sources that GW190814 represents. Astrophysical models predict that binaries with mass ratios similar to this event can form through several channels, but are unlikely to have formed in globular clusters. However, the combination of mass ratio, component masses, and the inferred merger rate for this event challenges all current models of the formation and mass distribution of compact-object binaries.

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.

Gravitational waves detected from the binary neutron star (NS) merger GW170817 constrained the NS equation of state by placing an upper bound on certain parameters, describing the binary’s tidal interactions. We show that the interpretation of the UV/optical/infrared counterpart of GW170817 with kilonova models, combined with new numerical-relativity results, imply a complementary lower bound on the tidal deformability parameter. The joint constraints tentatively rule out both extremely stiff and soft NS equations of state.

We present a systematic numerical relativity study of the mass ejection and the associated electromagnetic transients and nucleosynthesis from binary neutron star (NS) mergers. We find that a few 10⁻³ M _⊙ of material is ejected dynamically during the mergers. The amount and the properties of these outflows depend on binary parameters and on the NS equation of state (EOS). A small fraction of these ejecta, typically ∼10⁻⁶ M _⊙, is accelerated by shocks formed shortly after merger to velocities larger than 0.6c and produces bright radio flares on timescales of weeks, months, or years after merger. Their observation could constrain the strength with which the NSs bounce after merger and, consequently, the EOS of matter at extreme densities. The dynamical ejecta robustly produce second and third r-process peak nuclei with relative isotopic abundances close to solar. The production of light r-process elements is instead sensitive to the binary mass ratio and the neutrino radiation treatment. Accretion disks of up to ∼0.2 M _⊙ are formed after merger, depending on the lifetime of the remnant. In most cases, neutrino- and viscously driven winds from these disks dominate the overall outflow. Finally, we generate synthetic kilonova light curves and find that kilonovae depend on the merger outcome and could be used to constrain the NS EOS.

The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in 2019 April and lasting six months, O3b starting in 2019 November and lasting five months, and O3GK starting in 2020 April and lasting two weeks. In this paper we describe these data and various other science products that can be freely accessed through the Gravitational Wave Open Science Center at https://gwosc.org. The main data set, consisting of the gravitational-wave strain time series that contains the astrophysical signals, is released together with supporting data useful for their analysis and documentation, tutorials, as well as analysis software packages.