Abstract:In the absence of numerous gravitational-wave detections with confirmed electromagnetic counterparts, the "dark siren" method has emerged as a leading technique of gravitational-wave cosmology. The method allows redshift information of such events to be inferred statistically from a catalogue of potential host galaxies. Due to selection effects, dark siren analyses necessarily depend on the mass distribution of compact objects and the evolution of their merger rate with redshift. Informative priors on these quantities will impact the inferred posterior constraints on the Hubble constant ($H_0$). It is thus crucial to vary these unknown distributions during an $H_0$ inference. This was not possible in earlier analyses due to the high computational cost, restricting them to either excluding galaxy catalogue information, or fixing the gravitational-wave population mass distribution and risking introducing bias to the $H_0$ measurement. This paper introduces a significantly enhanced version of the Python package GWCOSMO, which allows joint estimation of cosmological and compact binary population parameters. This thereby ensures the analysis is now robust to a major source of potential bias. The gravitational-wave events from the Third Gravitational-Wave Transient Catalogue are reanalysed with the GLADE+ galaxy catalogue, and an updated, more reliable measurement of $H_0=69^{+12}_{-7}$ km s$^{-1}$ Mpc$^{-1}$ is found (maximum a posteriori probability and 68% highest density interval). This improved method will enable cosmological analyses with future gravitational-wave detections to make full use of the information available (both from galaxy catalogues and the compact binary population itself), leading to promising new independent bounds on the Hubble constant.

GWPopulation: Hardware agnostic population inference for compact binaries and beyond

Non-parametric inference of the population of compact binaries from gravitational wave observations using binned Gaussian processes

Massively parallel Bayesian inference for transient gravitational-wave astronomy

ICAROGW: A python package for inference of astrophysical population properties of noisy, heterogeneous and incomplete observations

Projections of the uncertainty on the compact binary population background using popstock

pygwb: Python-based library for gravitational-wave background searches

A physically modelled selection function for compact binary mergers in the LIGO-Virgo O3 run and beyond

Backward Population Synthesis: Mapping the Evolutionary History of Gravitational-Wave Progenitors

Parallelized Inference for Gravitational-Wave Astronomy

A novel approach to infer population and cosmological properties with gravitational waves standard sirens and galaxy surveys

Inferring the Astrophysical Population of Gravitational Wave Sources in the Presence of Noise Transients

PixelPop: High Resolution Nonparameteric Inference of Gravitational-Wave Populations in Multiple Dimensions

What to do when things get crowded? Scalable joint analysis of overlapping gravitational wave signals

Observed Gravitational-Wave Populations

Joint cosmological and gravitational-wave population inference using dark sirens and galaxy catalogues

Probing the large scale structure using gravitational-wave observations of binary black holes

Gravitational wave populations and cosmology with neural posterior estimation

Non-Stationary Astrophysical Stochastic Gravitational-Wave Background: A New Probe to the High Redshift Population of Binary Black Holes

gwforge: A user-friendly package to generate gravitational-wave mock data

Probing multiple populations of compact binaries with third-generation gravitational-wave detectors