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.

We use 47 gravitational wave sources from the Third LIGO–Virgo–Kamioka Gravitational Wave Detector Gravitational Wave Transient Catalog (GWTC–3) to estimate the Hubble parameter H(z), including its current value, the Hubble constant H ₀. Each gravitational wave (GW) signal provides the luminosity distance to the source, and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog. Using the binary black hole (BBH) redshifted masses, we simultaneously infer the source mass distribution and H(z). The source mass distribution displays a peak around 34 M _⊙, followed by a drop-off. Assuming this mass scale does not evolve with the redshift results in a H(z) measurement, yielding H 0 = 68 − 8 + 12 km s − 1 Mpc − 1 (68% credible interval) when combined with the H ₀ measurement from GW170817 and its electromagnetic counterpart. This represents an improvement of 17% with respect to the H ₀ estimate from GWTC–1. The second method associates each GW event with its probable host galaxy in the catalog GLADE+, statistically marginalizing over the redshifts of each event’s potential hosts. Assuming a fixed BBH population, we estimate a value of H 0 = 68 − 6 + 8 km s − 1 Mpc − 1 with the galaxy catalog method, an improvement of 42% with respect to our GWTC–1 result and 20% with respect to recent H ₀ studies using GWTC–2 events. However, we show that this result is strongly impacted by assumptions about the BBH source mass distribution; the only event which is not strongly impacted by such assumptions (and is thus informative about H ₀) is the well-localized event GW190814.

We report the observation of a coalescing compact binary with component masses 2.5–4.5 M _⊙ and 1.2–2.0 M _⊙ (all measurements quoted at the 90% credible level). The gravitational-wave signal GW230529_181500 was observed during the fourth observing run of the LIGO–Virgo–KAGRA detector network on 2023 May 29 by the LIGO Livingston observatory. The primary component of the source has a mass less than 5 M _⊙ at 99% credibility. We cannot definitively determine from gravitational-wave data alone whether either component of the source is a neutron star or a black hole. However, given existing estimates of the maximum neutron star mass, we find the most probable interpretation of the source to be the coalescence of a neutron star with a black hole that has a mass between the most massive neutron stars and the least massive black holes observed in the Galaxy. We provisionally estimate a merger rate density of 55 − 47 + 127 Gpc − 3 yr − 1 for compact binary coalescences with properties similar to the source of GW230529_181500; assuming that the source is a neutron star–black hole merger, GW230529_181500-like sources may make up the majority of neutron star–black hole coalescences. The discovery of this system implies an increase in the expected rate of neutron star–black hole mergers with electromagnetic counterparts and provides further evidence for compact objects existing within the purported lower mass gap.

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the l = m = 2 mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the l = 2, m = 1, 2 modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found, so we present 95% credible upper limits on the strain amplitudes h ₀ for the single-harmonic search along with limits on the pulsars’ mass quadrupole moments Q ₂₂ and ellipticities ε. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437−4715 and J0711−6830, which have spin-down ratios of 0.87 and 0.57, respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars, our limits are factors of ∼100 and ∼20 more constraining than their spin-down limits, respectively. For the dual-harmonic searches, new limits are placed on the strain amplitudes C ₂₁ and C ₂₂. For 23 pulsars, we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

Gravitational lensing by massive objects along the line of sight to the source causes distortions to gravitational wave (GW) signals; such distortions may reveal information about fundamental physics, cosmology, and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO-Virgo network. We search for repeated signals from strong lensing by (1) performing targeted searches for subthreshold signals, (2) calculating the degree of overlap among the intrinsic parameters and sky location of pairs of signals, (3) comparing the similarities of the spectrograms among pairs of signals, and (4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by (1) frequency-independent phase shifts in strongly lensed images, and (2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the nondetection of GW lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects.