Christopher Philip Luke Berry

Gravitational waves can be gravitationally lensed by massive objects along their path. Depending on the lens mass and the lens–source geometry, this can lead to the observation of a single distorted signal, multiple repeated events, or a mixture of both. We present the results for gravitational-wave lensing searches on the data from the first part of the fourth LIGO–Virgo–KAGRA observing run (O4a). We search for strongly lensed events in the newly acquired data by (1) searching for an overall phase shift present in an image formed at a saddle point of the lens potential, (2) looking for pairs of detected candidates with consistent frequency evolution, and (3) identifying subthreshold counterpart candidates to the detected signals. Beyond strong lensing, we also look for lensing-induced distortions in all detected signals using an isolated point-mass model. We do not find evidence for strongly lensed gravitational-wave signals and use this result to constrain the rate of detectable strongly lensed events and the merger rate density of binary black holes at high redshift. In the search for single distorted lensed signals, we find one outlier: GW231123_135430, for which we report more detailed investigations. While this event is interesting, the associated waveform uncertainties make its interpretation complicated, and future observations of the populations of binary black holes and of gravitational lenses will help determine the probability that this event could be lensed.

We analyze data from 142 of the 218 gravitational-wave (GW) sources in the fourth LIGO–Virgo–KAGRA Collaboration (LVK) Gravitational-Wave Transient Catalog (GWTC-4.0) to estimate the Hubble constant H ₀ jointly with the population properties of merging compact binaries. We measure the luminosity distance and redshifted masses of GW sources directly; in contrast, we infer GW source redshifts statistically through (i) location of features in the compact object mass spectrum and merger rate evolution, and (ii) identifying potential host galaxies in the GW localization volume. Probing the relationship between source luminosity distances and redshifts obtained in this way yields constraints on cosmological parameters. We also constrain parameterized deviations from general relativity that affect GW propagation, specifically those modifying the dependence of a GW signal on the source luminosity distance. Assuming our fiducial model for the source-frame mass distribution and using GW candidates detected up to the end of the fourth observing run (O4a), together with the GLADE+ all-sky galaxy catalog, we estimate H 0 = 75 . 4 − 9.1 + 12.8 ( 75 . 4 − 13.3 + 24.5 ) km s − 1 Mpc − 1 . This value is reported as a median with 68.3% (90%) symmetric credible interval, and includes combination with the H ₀ measurement from GW170817 and its electromagnetic counterpart. Using a parameterization of modified GW propagation in terms of the magnitude parameter Ξ ₀ , we estimate Ξ 0 = 1 . 4 − 0.4 + 1.0 ( 1 . 4 − 0.6 + 2.8 ) , where Ξ ₀ = 1 recovers the behavior of general relativity.

We detail the population properties of binary neutron star, neutron star–black hole binary, and binary black hole mergers using 158 events from the cumulative Gravitational-Wave Transient Catalog 4.0. The black hole primary mass distribution consists of a power-law-like continuum that steepens above 35 M _⊙ with overdensities at 10 M _⊙ and 35 M _⊙ . Binary black holes with primary masses near 10 M _⊙ are more likely to have less massive secondaries, with a mass ratio distribution peaking at q = 0.7 1 − 0.11 + 0.12 , potentially a signature of stable mass transfer during binary evolution. Black hole spins are inferred to be nonextremal, with 90% of black holes having χ < 0.6, and preferentially aligned with binary orbits, implying many merging binaries form in isolation. However, we find that a significant fraction, 0.23–0.41, of binaries have negative effective inspiral spins, suggesting many could be formed dynamically in gas-free environments. We find evidence for correlation between effective inspiral spin and mass ratio, driven either by variation in the mode of the distribution or the width. The binary black hole merger rate increases with redshift, consistent with the cosmic star formation density. While there is no evidence of the mass spectrum evolving with redshift, the distribution of effective inspiral spin is found to broaden as redshift increases out to z ≈ 1. We infer the local merger rates (at z = 0) to be 8.8–250 Gpc ⁻³ yr ⁻¹ for binary neutron stars, 9.3–86 Gpc ⁻³ yr ⁻¹ for neutron star–black hole binaries, and 13–26 Gpc ⁻³ yr ⁻¹ for binary black holes; all values reflect central 90% credible intervals.

Gravitational waves from extreme mass-ratio inspirals (EMRIs) will enable subpercent measurements of massive black hole parameters and provide access to the demographics of compact objects in galactic nuclei. During the LISA mission, multiple EMRIs are expected to be detected, allowing statistical studies of massive black hole populations and their formation pathways. We perform hierarchical Bayesian inference on simulated EMRI catalogs to assess how well LISA could constrain the astrophysical population using parameterized population models. We test our inference framework on a variety of populations, including heterogeneous and homogeneous mixtures of parameterized subpopulations, and scenarios in which the assumed model is deliberately misspecified. Our results show that population parameters governing distributions with sharp features can be tightly constrained. Mixed populations can be disentangled with as few as ∼20 detections, and even with model misspecification, the inference retains sensitivity to key population features. These results demonstrate the capabilities and limitations of EMRI population inference, providing guidance for constructing realistic astrophysical population models for LISA analysis.

LIGO, Virgo, and KAGRA form a network of gravitational-wave observatories. Data and analysis results from this network are made publicly available through the Gravitational Wave Open Science Center. This paper describes open data from this network, including the addition of data from the first part of the fourth observing run and selected periods from the preceding engineering run, collected from 2023 May to 2024 January. The public dataset includes calibrated strain time series for each instrument, data from additional channels used for noise subtraction and detector characterization, and analysis data products from version 4.0 of the Gravitational-Wave Transient Catalog.

Version 4.0 of the Gravitational-Wave Transient Catalog (GWTC-4.0) adds new candidates detected by the LIGO, Virgo, and KAGRA observatories through the first part of the fourth observing run (O4a: 2023 May 24 15:00:00 to 2024 January 16 16:00:00 UTC) and a preceding engineering run. In these new data, we find 128 compact binary coalescence candidates that are identified by at least one of our search algorithms with a probability of astrophysical origin p _astro ≥ 0.5 and that are not vetoed during event validation. We also provide detailed source property measurements for 86 of these that have a false-alarm rate −1 . Based on the inferred component masses, these candidates are consistent with signals from binary black holes (BBHs) and neutron star–black hole binaries (GW230518_125908 and GW230529_181500). Median-inferred component masses of BBHs in the catalog now range from 5.79M⊙ (GW230627_015337) to 137 M ⊙ (GW231123_135430), while GW231123_135430 was probably produced by the most massive binary observed in the catalog. For the first time, we have discovered BBH signals with network signal-to-noise ratio exceeding 30, GW230814_230901 and GW231226_101520, enabling high-fidelity studies of the waveforms and astrophysical properties of these systems. Combined with the 90 candidates included in GWTC-3.0, the catalog now contains 218 candidates with p _astro ≥ 0.5 and not otherwise vetoed, more than doubling the size of the catalog and further opening our view of the gravitational-wave Universe.

GW230814, detected by the LIGO Livingston observatory with a signal-to-noise ratio of 42.4, represents the loudest gravitational-wave signal in the GWTC-4.0 catalog. Its source is consistent with a binary black hole coalescence with component masses m 1 = 33 . 7 − 2.2 + 2.9 M ⊙ and m 2 = 28 . 2 − 3.1 + 2.2 M ⊙ and a small effective inspiral spin χ eff = − 0.0 1 − 0.07 + 0.06 . The high signal-to-noise ratio enabled us to detect an ℓ = ∣ m ∣ = 4 mode in the inspiral–merger–ringdown signal for the first time (with Bayes factor ≈10), as well as enabling a range of tests of consistency between theoretical predictions and the observed waveform. While most of these tests show agreement with theoretical predictions, there are suggestions of minor deviations in the ringdown phase. Simulations that incorporate general relativity and realistic detector noise reproduce similar deviations, suggesting that they do not constitute evidence for a breakdown of general relativity. The observation of GW230814 demonstrates that the unprecedented sensitivity of the detectors enables highly significant detections with a single observatory. However, without corroborating data from a multidetector network, the ability to draw rigorous conclusions about fundamental physics remains severely limited.

The Gravitational-Wave Transient Catalog (GWTC) is a collection of candidate gravitational-wave transient signals identified and characterized by the LIGO–Virgo–KAGRA Collaboration. Producing the contents of the GWTC from detector data requires complex analysis methods. These comprise techniques to model the signal; identify the transients in the data; evaluate the quality of the data and mitigate possible instrumental issues; infer the parameters of each transient; compare the data with the waveform models for compact binary coalescences; and handle the large amount of results associated with all of these different analyses. In this paper, we describe the methods employed to produce the catalog’s fourth release, GWTC-4.0, focusing on the analysis of the first part of the fourth observing run of Advanced LIGO, Advanced Virgo and KAGRA.

The discovery of joint sources of high-energy neutrinos and gravitational waves has been a primary target for the LIGO, Virgo, KAGRA, and IceCube observatories. The joint detection of high-energy neutrinos and gravitational waves would provide insight into cosmic processes, from the dynamics of compact object mergers and stellar collapses to the mechanisms driving relativistic outflows. The joint detection of multiple cosmic messengers can also elevate the significance of the common observation even when some or all of the constituent messengers are subthreshold, i.e., not significant enough to declare their detection individually. Using data from the LIGO, Virgo, and IceCube observatories, including subthreshold events, we searched for common sources of gravitational waves and high-energy neutrinos during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Our search did not identify significant joint sources. We derive constraints on the rate densities of joint sources. Our results constrain the isotropic neutrino emission from gravitational-wave sources for very high values of the total energy emitted in neutrinos (>10 ⁵² –10 ⁵⁴ erg).