Unbiased likelihood-free inference of the Hubble constant from light standard sirens

Abstract

Multi-messenger observations of binary neutron star mergers offer a promising path towards resolution of the Hubble constant ($H_0$) tension, provided their constraints are shown to be free from systematics such as the Malmquist bias. In the traditional Bayesian framework, accounting for selection effects in the likelihood requires calculation of the expected number (or fraction) of detections as a function of the parameters describing the population and cosmology; a potentially costly and/or inaccurate process. This calculation can, however, be bypassed completely by performing the inference in a framework in which the likelihood is never explicitly calculated, but instead fit using forward simulations of the data, which naturally include the selection. This is Likelihood-Free Inference (LFI). Here, we use density-estimation LFI, coupled to neural-network-based data compression, to infer $H_0$ from mock catalogues of binary neutron star mergers, given noisy redshift, distance and peculiar velocity estimates for each object. We demonstrate that LFI yields statistically unbiased estimates of $H_0$ in the presence of selection effects, with precision matching that of sampling the full Bayesian hierarchical model. Marginalizing over the bias increases the $H_0$ uncertainty by only $6\%$ for training sets consisting of $O(10^4)$ populations. The resulting LFI framework is applicable to population-level inference problems with selection effects across astrophysics.

19 pages, 8 figures, comments welcome

References in corpus (15)

• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• Multi-messenger Observations of a Binary Neutron Star Merger
• In the Realm of the Hubble tension $-$ a Review of Solutions
• Exploring the Sensitivity of Next Generation Gravitational Wave Detectors
• A gravitational-wave standard siren measurement of the Hubble constant
• Cosmic Distances Calibrated to 1% Precision with Gaia EDR3 Parallaxes and Hubble Space Telescope Photometry of 75 Milky Way Cepheids Confirm Tension with LambdaCDM
• The trouble with $H_0$
• Scientific Objectives of Einstein Telescope
• Fast likelihood-free cosmology with neural density estimators and active learning
• Characterization of systematic error in Advanced LIGO calibration
• Nuisance hardened data compression for fast likelihood-free inference
• Accuracy Requirements for Empirically-Measured Selection Functions
• How Many Kilonovae Can Be Found in Past, Present, and Future Survey Datasets?
• Heavy double neutron stars: birth, mid-life and death
• Prospects of the local Hubble parameter measurement using gravitational waves from double neutron stars

Cited by in corpus (6)

• Challenges for $\Lambda$CDM: An update
• Forward modeling of galaxy populations for cosmological redshift distribution inference
• Characterizing the observation bias in gravitational-wave detections and finding structured population properties
• Fast and robust Bayesian Inference using Gaussian Processes with GPry
• Calibrating cosmological simulations with implicit likelihood inference using galaxy growth observables
• Likelihood-free Forward Modeling for Cluster Weak Lensing and Cosmology