🔭 astrophysics

Improved constraint on the Hubble constant from dark sirens with LIGO/Virgo/KAGRA O4a

This paper presents an improved measurement of the Hubble constant ( $H_0 = 69.9^{+4.1}_{-4.0}$ km/s/Mpc) by combining statistical dark siren analyses of 17 well-localized gravitational-wave events from LIGO/Virgo/KAGRA with the bright siren GW170817, utilizing methodological enhancements such as $r$ -band luminosity weighting and refined selection effect treatments to reduce uncertainty by approximately 11% compared to the GW170817-only result.

Original authors: V. Alfradique, C. R. Bom, G. Teixeira, A. Santos

Published 2026-03-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: V. Alfradique, C. R. Bom, G. Teixeira, A. Santos

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe is 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 ( $H_0$ ).

Here's the problem: When scientists look at the "baby pictures" of the universe (the Cosmic Microwave Background), they get one speed. When they look at the "adult pictures" (nearby exploding stars called supernovae), they get a different speed. They disagree by about 5–6 standard deviations, which in science is like two people arguing over the price of a coffee where one says $3 and the other says $5, and both are absolutely sure they are right. This is the famous "Hubble Tension."

This paper is about a new, clever way to settle the argument using Gravitational Waves (ripples in spacetime caused by colliding black holes).

The Two Types of "Sirens"

The authors use a method called "Standard Sirens." Think of a gravitational wave event like a lighthouse or a siren on a ship.

The Bright Siren: Sometimes, when two black holes or neutron stars collide, they flash light (like a camera flash). We can see the flash, identify the host galaxy, and know exactly how far away it is. This is like seeing a lighthouse and knowing its exact brightness, so you can calculate the distance. The most famous one is GW170817.
The Dark Siren: Most of the time, black holes collide in total darkness. There is no flash. We hear the "siren" (the gravitational wave) but can't see the ship. We know the sound tells us the distance, but we don't know which galaxy it came from. It's like hearing a siren in a foggy city and trying to guess which street it's on.

How They Solved the "Dark" Mystery

The authors took 17 "Dark Sirens" (collisions with no light) from the LIGO/Virgo/KAGRA detectors. Since they couldn't see the host galaxy, they had to play a game of probability:

The Map: They looked at the "foggy" area of the sky where the sound came from.
The Guest List: They pulled up a massive digital catalog of every galaxy in that area (from the DESI Legacy Imaging Survey).
The Guessing Game: They asked, "Which of these galaxies is the most likely host?"
- Old Method: Assume every galaxy has an equal chance.
- New Method (The Upgrade): They realized that heavy black holes are more likely to be born in bright, massive galaxies (like how big parties happen in big houses, not tiny apartments). So, they gave "weight" to the brighter galaxies. It's like betting more on the big house being the party location.

The New Tricks They Used

The paper introduces three major upgrades to their detective work:

The "Heavy" Filter: They used a new AI model to estimate the distance to these galaxies more accurately.
The "Mass" Clue: They didn't just listen to the sound; they analyzed the weight of the colliding black holes. Heavier black holes have a specific "signature" that helps narrow down the distance. It's like knowing the siren belongs to a specific type of truck, which helps you guess how far away it is.
Fixing the "Missing" Data: Their galaxy map wasn't perfect; it missed some faint, distant galaxies. Instead of ignoring this, they mathematically "filled in the blanks" with fake galaxies that followed the same rules as the real ones, ensuring their sample wasn't biased toward only the bright, easy-to-see ones.

The Results: A Better Estimate

After crunching the numbers on these 17 dark events, they got a result:

Dark Sirens Alone: The universe is expanding at about 78 km/s/Mpc (with a margin of error).
The Grand Finale: They combined their 17 dark sirens with the one famous "Bright Siren" (GW170817).

The Final Answer: By mixing the "loud" event (Bright Siren) with the "quiet" events (Dark Sirens), they refined the measurement to 69.9 km/s/Mpc.

Why This Matters

The Improvement: This new result is about 11% more precise than using the single bright event alone.
The Verdict: This number sits right in the middle of the two conflicting measurements (Planck vs. SH0ES). It doesn't fully solve the Hubble Tension yet, but it proves that Gravitational Waves are a powerful, independent referee that can help settle the debate.
The Future: As we detect more of these "dark sirens" in the coming years, and as our galaxy maps get deeper and clearer, we expect to pin down the expansion rate of the universe with incredible precision, potentially revealing if our current understanding of the universe (Dark Energy, General Relativity) needs a rewrite.

In short: The authors used a clever statistical trick to turn 17 invisible black hole collisions into a precise ruler for the universe, helping us measure the cosmic expansion rate with greater confidence than ever before.

1. Problem Statement

The paper addresses the "Hubble tension," a persistent 5–6 $\sigma$ discrepancy between the value of the Hubble constant ( $H_0$ ) inferred from early-universe Cosmic Microwave Background (CMB) measurements (Planck: $\sim 67.4$ km/s/Mpc) and late-universe local distance ladder measurements (SH0ES: $\sim 73.0$ km/s/Mpc).

Gravitational Waves (GWs) from compact binary coalescences offer an independent "standard siren" method to measure $H_0$ without relying on the cosmic distance ladder. While "bright sirens" (events with electromagnetic counterparts, like GW170817) provide direct redshift measurements, they are rare. "Dark sirens" (events without counterparts) rely on statistical methods to associate the GW source with potential host galaxies in a catalog. This paper aims to improve the precision of dark siren $H_0$ constraints by incorporating the latest data from the fourth LIGO–Virgo–KAGRA (LVK) observing run (O4a) and refining the statistical methodology.

2. Methodology

The authors employ the statistical dark siren method, combining GW luminosity distance measurements with photometric redshifts from galaxy catalogs. The analysis involves several key methodological advancements over previous works (specifically Bom et al. [51]):

Data Sample:
- GW Events: A total of 17 well-localized binary black hole (BBH) events. This includes 7 new events from the O4a run (GWTC-4) and 10 events from previous runs (O1–O3).
- Selection Criteria: Events were selected based on luminosity distance ( $d_L < 1500$ Mpc), sky localization area ( $A_{50\%} < 1000$ deg $^2$ ), and sky coverage by the galaxy catalog ( $>70\%$ ).
- Galaxy Catalog: The DESI Legacy Imaging Survey (LS) (Data Releases 9 and 10). The authors constructed a volume-limited sample by applying apparent and absolute magnitude cuts based on a Schechter function fit to the galaxy luminosity distribution.
Photometric Redshift (Photo-z) Estimation:
- A Mixture Density Network (MDN), a deep learning model, was used to estimate probability density functions (PDFs) for galaxy redshifts rather than single point estimates.
- For the northern sky (declination $>30^\circ$ ), a new MDN model was trained exclusively on LS-DR9 data to improve performance in that region, differing from previous combined training sets.
Likelihood Construction & Weighting:
- Luminosity Weighting: The authors introduced $r$ -band luminosity weighting for potential host galaxies. This assumes BBH mergers occur preferentially in more luminous (and thus more massive) galaxies, assigning higher statistical weight to brighter galaxies.
- Extended GW Likelihood: Instead of a simple Gaussian approximation dependent only on luminosity distance and sky position, the authors used a full GW likelihood derived from a Kernel Density Estimate (KDE) of posterior samples. This incorporates information from the binary black hole component masses and reweights them using a population prior (Broken Power Law + 2 Peaks model from GWTC-4).
- Selection Effects: The analysis explicitly accounts for the incompleteness of the magnitude-limited LS catalog and the GW detector selection function ( $\beta(H_0)$ ) to correct for biases.
Synthetic Galaxies: For events where the LS catalog coverage was incomplete (approaching the 70% threshold), synthetic galaxies were injected into uncovered regions to maintain statistical completeness, following established procedures.

3. Key Contributions

Inclusion of O4a Data: The first application of the statistical dark siren method to the seven well-localized BBH events from the O4a observing run (GWTC-4).
Methodological Refinements:
- Implementation of $r$ -band luminosity weighting to better model the host galaxy population.
- Adoption of an extended GW likelihood that utilizes BBH mass population priors, moving beyond the standard Gaussian distance approximation.
- Development of a region-specific MDN for photo-z estimation in the northern sky.
Rigorous Selection Bias Treatment: A consistent treatment of selection effects for both the GW detection network and the magnitude-limited galaxy catalog, ensuring the inference remains unbiased by catalog incompleteness.

4. Results

The paper presents $H_0$ constraints under various configurations:

Dark Sirens Only (17 events):
- Without Luminosity Weighting: $H_0 = 78.8^{+14.6}_{-12.2}$ km/s/Mpc (using full GW likelihood).
- With $r$ -band Luminosity Weighting: $H_0 = 78.2^{+12.0}_{-11.0}$ km/s/Mpc.
- Observation: Luminosity weighting reduced the uncertainty by approximately 14% when using the full GW likelihood. The results are consistent with both Planck and SH0ES within $1\sigma$ .
Combined Bright + Dark Sirens:
- The authors combined the luminosity-weighted dark siren sample with the bright siren GW170817.
- They applied corrections for the host galaxy's peculiar velocity and included constraints on the jet viewing angle from VLBI data.
- Final Constraint: $H_0 = 69.9^{+4.1}_{-4.0}$ km/s/Mpc.
Impact Analysis:
- The addition of the 17 dark sirens to the GW170817 measurement improved the uncertainty by approximately 11% compared to GW170817 alone.
- The use of the full GW likelihood (incorporating mass information) shifted the posterior peak slightly higher compared to the Gaussian approximation but remained statistically consistent.
- The inclusion of synthetic galaxies for incomplete sky coverage had a negligible impact on the final uncertainty ( $<1\%$ ).

5. Significance

Precision Improvement: The combined result ( $H_0 = 69.9 \pm 4.0$ km/s/Mpc) represents a significant step toward percent-level precision in GW cosmology. The uncertainty is reduced by ~11% relative to the bright siren-only result, demonstrating the complementarity of bright and dark sirens.
Validation of Methodology: The study validates that luminosity weighting and the inclusion of BBH mass population priors can tighten constraints without introducing significant systematic bias, provided the galaxy catalog is sufficiently complete and volume-limited.
Future Outlook: The paper highlights that as the number of well-localized events increases in future observing runs (O4b, O5) and with deeper galaxy surveys (e.g., Vera C. Rubin Observatory), GW cosmology will become a critical independent probe for resolving the Hubble tension, provided systematic effects like host galaxy weighting and selection biases are carefully controlled.

In summary, this work refines the statistical dark siren framework using the latest LVK data and advanced machine learning techniques, yielding a robust and precise measurement of the Hubble constant that bridges the gap between early and late-universe cosmological probes.