Systematic bias in dark siren statistical methods and its impact on Hubble constant measurement
Phys. Rev. D 112, 063561 – Published 29 September, 2025
DOI: https://doi.org/10.1103/vd36-3mys
Abstract
The advent of multimessenger cosmology, marked by the detection of GW170817, demonstrated that standard sirens are a valuable cosmological probe. In the absence of an electromagnetic counterpart identification, gravitational waves (GWs) carry valuable information through the dark siren (DS) approach, where the source redshift is estimated using galaxy catalogs of potential hosts within the localization volume. However, the DS analysis can be affected by galaxy catalog incompleteness at the limits of gravitational-wave detectability, potentially introducing biases in the constraints on cosmological parameters. Focusing on GWs from binary black holes (BBH) detected by the LIGO, Virgo, and KAGRA (LVK) collaboration, we explore the possible systematic biases in the measurement of the Hubble constant (). These biases may arise from (1) the incompleteness of catalogs due to the apparent magnitude thresholds of optical telescope sensitivity, and (2) the use of incorrect weighting schemes (e.g., using star formation or stellar mass as tracers of the host galaxy) for each potential host. We found that an unbiased estimate of can be obtained when the corrected weighting scheme is applied to a complete or volume-limited catalog. We use a complete galaxy catalog covering 90% of the localization probability () for each GW detection, employing stellar mass as a tracer. Our results show that a sample of 100 binary black hole events with and measured luminosity distances below 1600(2500) Mpc, detected by the LVK at O4(O5) sensitivity, can provide a percent-level measurement of , with a precision of approximately 3%(1%). This number of detections is expected to be accumulated after approximately 8 and 3 years of observations with the LVK at O4 and O5 sensitivity, respectively. The O5 run provides a reduction in the uncertainty by compared to the O4-like configuration. The precision increases to approximately 6% when it is assumed that every galaxy has an equal probability of being the host.