GWTC-5.0: Constraints on the Cosmic Expansion Rate and Modified Gravitational-wave Propagation

Using 236 gravitational-wave sources from the GWTC-5.0 catalog, this study refines the Hubble constant estimate to $71.0^{+9.0}_{-7.1}$ km s$^{-1}$ Mpc$^{-1}$ with a 25.7% reduction in uncertainty compared to previous results and confirms no deviations from general relativity in gravitational-wave propagation.

The Gist

Imagine the universe as a giant, expanding balloon. For decades, scientists have been trying to measure exactly how fast this balloon is inflating. This speed is called the Hubble Constant. But here's the problem: when they measure it using light from the very beginning of the universe (the Cosmic Microwave Background), they get one answer. When they measure it using light from nearby exploding stars (Supernovae), they get a different, slightly faster answer. This disagreement is known as the "Hubble Tension," and it's one of the biggest mysteries in physics today.

This paper, written by the LIGO, Virgo, and KAGRA collaborations, introduces a new, independent way to measure that expansion speed using gravitational waves—ripples in the fabric of spacetime caused by massive objects crashing into each other.

Here is a simple breakdown of what they did and what they found, using some everyday analogies.

1. The "Standard Siren" Analogy

Usually, to measure distance in space, astronomers use a "cosmic distance ladder." They start with nearby objects they know the size of, then use those to measure further objects, and so on. It's like trying to measure the length of a football field by using a ruler, then a tape measure, then a car odometer, hoping each step is accurate.

Gravitational waves offer a shortcut. When two black holes or neutron stars merge, they create a sound (a "chirp") that travels through space. Because we know the physics of how these objects merge, the "volume" of the sound tells us exactly how far away they are. The scientists call these Standard Sirens.

• The Problem: The sound tells us the distance, but it doesn't tell us the speed at which the universe is expanding. To get that, we need to know the redshift (how much the universe stretched the signal while it traveled).
• The Catch: The gravitational wave signal itself is "degenerate." It's like hearing a siren in a fog; you can tell how loud it is, but you can't tell if it's a loud siren far away or a quiet siren close by. The signal mixes up the mass of the objects with their distance.

2. Two Ways to Solve the Puzzle

To break this "fog," the team used two clever tricks with 236 gravitational wave events from their new catalog (GWTC-5.0):

Method A: The "Spectral Siren" (The Crowd's Voice)
Imagine you walk into a room full of people shouting. You don't know who is where, but you notice a pattern: most people are shouting at a specific pitch, with a few shouting higher or lower.

• How it works: The scientists looked at the "mass spectrum" of all the merging black holes. They know there are specific "favorite" masses where black holes tend to form (like a crowd preferring a certain pitch). By analyzing the pattern of masses across all 236 events, they could statistically figure out how much the universe stretched the signal. It's like deducing the size of the room by listening to the echo patterns of the whole crowd, rather than asking one person.

Method B: The "Dark Siren" (The Map Search)
Imagine you hear a siren but can't see the source. You pull out a map and look for the most likely houses in the direction the sound came from.

• How it works: For each gravitational wave event, the team looked at the "sky map" to see which galaxies were in the area. They used two massive galaxy catalogs (like a phone book for the universe): one called GLADE+ (a wide but shallow list) and one called DES Year 6 (a deep, detailed list of a smaller area). They matched the gravitational wave event to the galaxies in that spot to guess the redshift.
• The Improvement: In this new study, the "sky maps" for the new events are much sharper (better localization) than before, thanks to the Virgo detector joining in. This is like going from a blurry photo of a neighborhood to a high-definition street view, making it much easier to find the right house.

3. The Results: A New Measurement

By combining these methods, the team calculated the Hubble Constant ($H_0$).

• The Result: They found the universe is expanding at 71.0 km/s per Megaparsec.
• The Precision: The uncertainty (the "fuzziness" of the measurement) has dropped by 25.7% compared to their previous study.
• The Comparison: This result sits right in the middle of the two conflicting previous measurements (the "early universe" vs. "local universe" values). It doesn't fully solve the tension yet, but it provides a strong, independent check that leans slightly toward the faster, local measurement.

Key Takeaway: For the first time, the team found that using just the "Dark Sirens" (statistical methods without a visible light counterpart) gave a tighter, more precise constraint on the expansion rate than the single "Bright Siren" event (GW170817) they had previously relied on. It's like finally having enough data points to draw a clear line, rather than guessing based on a single dot.

4. Checking the Rules of Gravity

The paper also asked a second question: Does gravity behave exactly as Einstein predicted?

• The Test: In Einstein's General Relativity, gravity waves and light waves travel at the same speed and lose energy in the same way as they cross the universe. Some alternative theories suggest gravity might get "friction" or change strength over vast distances.
• The Analogy: Imagine running a race. If Einstein is right, you and a light beam should finish at the exact same time and with the same energy. If modified gravity is right, you might arrive slightly tired or slower.
• The Result: The scientists found no evidence that gravity behaves differently than Einstein predicted. The "friction" is zero. The universe is playing by the standard rules of General Relativity, at least on the scales they tested.

Summary

This paper is a major step forward in "Gravitational Wave Cosmology." By listening to the "chirps" of 236 cosmic collisions and cross-referencing them with galaxy maps and statistical patterns, the team has:

1. Measured the expansion rate of the universe with greater precision than ever before using only gravitational waves.
2. Confirmed that Einstein's theory of gravity holds up, with no signs of "friction" slowing down gravitational waves.

They are essentially tuning the universe's "speedometer" with a new, independent tool, helping to resolve one of the biggest debates in modern physics.

Technical Summary

Technical Summary: GWTC-5.0: Constraints on the Cosmic Expansion Rate and Modified Gravitational-wave Propagation

Problem Statement
The measurement of the Hubble constant ($H_0$) remains a critical challenge in cosmology due to the persistent discrepancy between early-Universe measurements from the Cosmic Microwave Background (CMB) and local measurements from standardizable sources like Type Ia supernovae. Gravitational waves (GWs) from compact binary coalescences (CBCs) offer a "standard siren" approach to measure $H_0$ independently of the cosmic distance ladder. However, GW signals alone provide luminosity distance but lack intrinsic redshift information due to the mass-redshift degeneracy. This paper addresses the challenge of breaking this degeneracy using the fifth LIGO–Virgo–KAGRA Collaboration (LVK) Gravitational-Wave Transient Catalog (GWTC-5.0) to constrain $H_0$ and test for deviations from General Relativity (GR) in GW propagation.

Methodology
The analysis employs a hierarchical Bayesian framework to jointly infer cosmological parameters and population-level properties of GW sources from 236 CBC candidates (235 dark sirens and 1 bright siren, GW170817) with a false alarm rate (FAR) < 0.25 yr$^{-1}$.

1. Redshift Inference Strategies:

   • Spectral Sirens: The method breaks the mass-redshift degeneracy by utilizing features in the intrinsic mass spectrum of the CBC population. The analysis assumes that the observed mass distribution is a redshifted version of a source-frame distribution characterized by specific population parameters.
   • Dark Sirens: For events without electromagnetic (EM) counterparts, redshifts are inferred by associating the GW sky localization with galaxy catalogs. Two catalogs were utilized: GLADE+ (Ks-band, nearly full-sky) and the Dark Energy Survey (DES) Year 6 Gold catalog (r-band, deep southern sky). The analysis constructs redshift priors based on the spatial distribution of galaxies, applying luminosity weighting ($\epsilon=1$) or uniform weighting ($\epsilon=0$).

2. Population Models:
   The study tests three phenomenological mass distribution models:

   • Multi Peak (MLTP): A power-law distribution with two Gaussian peaks, applied to binary black hole (BBH) candidates.
   • FullPop-4.0: A unified model covering the full mass spectrum of BNS, NSBH, and BBH mergers, featuring a broken power law, a dip function for the mass gap, and two Gaussian peaks.
   • FullPop 3 Peaks: An extension of FullPop-4.0 adding a third Gaussian peak to capture additional mass spectrum structure.
     Spin distributions were also modeled using Gaussian and "Transition" models (where spin magnitude depends on mass), though the fiducial analysis primarily relies on the FullPop-4.0 mass model with uniform spin priors.

3. Modified Gravity Tests:
   The paper constrains deviations from GR affecting GW propagation by parameterizing the ratio of GW luminosity distance ($D_L^{GW}$) to EM luminosity distance ($D_L^{EM}$). Two parametrizations were tested:

   • $\Xi_0-n$: A phenomenological model where the ratio approaches $\Xi_0$ at high redshift.
   • $\alpha_M$: A model inspired by Horndeski gravity, where the effective Planck mass evolves with redshift, parameterized by $c_M$.

Key Contributions

• Dataset Expansion: The analysis utilizes the full population of GWTC-5.0, comprising 236 candidates. This includes 94 events from the O4b observing run, which feature significantly improved sky localization compared to previous runs due to Virgo's participation.
• Unified Population Modeling: The study advances the use of unified population models (FullPop-4.0 and FullPop 3 Peaks) that simultaneously analyze BNS, NSBH, and BBH populations, allowing for a more robust reconstruction of the mass spectrum and its features.
• Galaxy Catalog Comparison: This work provides the first direct comparison of cosmological constraints derived from the GLADE+ Ks-band catalog versus the deeper DES r-band catalog, evaluating the impact of sky coverage versus depth and completeness.
• Spin Modeling: The paper introduces the modeling of spin magnitude distributions (Gaussian vs. Transition) within the cosmological inference framework, finding evidence favoring the Transition model.

Results

• Hubble Constant ($H_0$):

  • Combining the 235 dark sirens with the bright siren GW170817, the fiducial analysis (FullPop-4.0 model + DES r-band catalog with luminosity weighting) yields:
    $$H_0 = 71.0^{+9.0}_{-7.1} \text{ km s}^{-1} \text{ Mpc}^{-1}$$
  • The dark siren analysis alone provides $H_0 = 67.6^{+13.6}_{-12.7} \text{ km s}^{-1} \text{ Mpc}^{-1}$, which is tighter than the constraint from GW170817 alone ($79.1^{+27.6}_{-12.4}$).
  • The inclusion of the DES r-band catalog improves the $H_0$ constraint by 7.4% compared to the spectral siren result alone. The total uncertainty reduction compared to the GWTC-4.0 measurement is 25.7%.
  • The results are statistically consistent with both Planck (CMB) and SH0ES (local) measurements at the 68% credible interval.

• Modified Gravity:

  • No departures from GR were found. The constraints on the $\Xi_0-n$ parametrization are $\Xi_0 = 1.1^{+0.6}_{-0.3}$ (wide $H_0$ prior) and $\Xi_0 = 1.0^{+0.3}_{-0.2}$ (narrow $H_0$ prior).
  • The constraint on the Horndeski parameter is $c_M = -0.4^{+1.6}_{-1.3}$ (wide prior) and $c_M = -0.1^{+1.0}_{-0.8}$ (narrow prior).
  • These results represent a 35–50% improvement in constraints on modified gravity parameters compared to GWTC-4.0, driven primarily by the increased number of high-redshift events.

• Population Parameters:

  • The FullPop 3 Peaks model shows a mild preference over FullPop-4.0 in Bayesian evidence, but FullPop-4.0 is retained as the fiducial model for direct comparability with previous work.
  • The Transition spin model is strongly favored over the Gaussian model ($\log_{10} \mathcal{B} \approx 3.16$), though the impact of spin modeling on the $H_0$ constraint is currently modest.

Significance
The paper claims a pivotal milestone in dark siren cosmology: for the first time, the constraint on $H_0$ derived from a sample of $\sim$200 dark sirens is tighter than the constraint derived from the single bright siren event GW170817. This demonstrates that statistical methods utilizing population features and galaxy catalogs can now rival the precision of individual multi-messenger events. Furthermore, the work establishes that GW observations provide an independent probe of cosmological modified gravity, yielding constraints on GW propagation that are competitive with, and complementary to, large-scale structure and CMB analyses. The authors note that while galaxy catalog information currently provides a subdominant improvement (7.4%) over spectral siren constraints, future Stage IV surveys with deeper and wider coverage are expected to significantly enhance the constraining power of the dark siren method.