A combined multimessenger analysis: Neutron stars are stellar remnants with densities greater than that of an atomic nucleus. The properties of matter under such extreme conditions are poorly understood and inaccessible to terrestrial laboratories. Dietrich et al. developed a framework to combine multiple constraints on the masses and radii of neutron stars, including data from gravitational waves, electromagnetic observations, and theoretical nuclear physics calculations. They used this analysis to constrain the neutron-star equation of state and also improved the precision on the gravitational wave (standard siren) measurement of the Hubble constant—the expansion rate of the Universe. Science , this issue p. 1450

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.

Accurate and fast gravitational waveform (GW) models are essential to extract information about the properties of compact binary systems that generate GWs. Building on previous work, we present an extension of the NRTidal model for binary neutron star (BNS) waveforms. The upgrades are (i) a new closed-form expression for the tidal contribution to the GW phase which includes further analytical knowledge and is calibrated to more accurate numerical relativity data than previously available; (ii) a tidal correction to the GW amplitude; and (iii) an extension of the spin-sector incorporating equation-of-state-dependent finite size effects at quadrupolar and octupolar order; these appear in the spin-spin tail terms and cubic-in-spin terms, both at 3.5 PN. We add the new description to the precessing binary black hole waveform model IMRPhenomPv2 to obtain a frequency-domain precessing binary neutron star model. In addition, we extend the SEOBNRv4_ROM and IMRPhenomD aligned-spin binary black hole waveform models with the improved tidal phase corrections. Focusing on the new IMRPhenomPv2_NRTidalv2 approximant, we test the model by comparing with numerical relativity waveforms as well as hybrid waveforms combining tidal effective-one-body and numerical relativity data. We also check consistency against a tidal effective-one-body model across large regions of the BNS parameter space.

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.