No Arabic abstract
We propose a novel approach to accurately pin down the systematics due to the peculiar velocities of galaxies in measuring the Hubble constant from nearby galaxies in current and future gravitational-wave (GW) standard-siren experiments. Given the precision that future GW standard-siren experiments aim to achieve, the peculiar velocities of nearby galaxies will be a major source of uncertainty. Unlike the conventional backward reconstruction that requires additional redshift-independent distance indicators to recover the peculiar velocity field, we forwardly model the peculiar velocity field by using a high-fidelity mock galaxy catalog built from high-resolution dark matter only (DMO) N-body simulations with a physically motivated subhalo abundance matching technique without introducing any free parameters. Our mock galaxy catalog can impressively well reproduce the observed spectroscopic redshift space distortions (RSDs) in highly non-linear regimes down to very small scales, which is a robust test of the velocity field of our mock galaxy catalog. Based on this mock galaxy catalog, we accurately, for the first time, measure the peculiar velocity probability distributions for the SDSS main galaxy samples. We find that the systematics induced by the peculiar velocities of SDSS like galaxies on the measured Hubble constant can be reduced to below $1%$($1sigma$) for GW host galaxies with a Hubble flow redshift just above $0.13$, a distance that can be well probed by future GW experiments and galaxy surveys.
In this work we investigate the systematic uncertainties that arise from the calculation of the peculiar velocity when estimating the Hubble constant ($H_0$) from gravitational wave standard sirens. We study the GW170817 event and the estimation of the peculiar velocity of its host galaxy, NGC 4993, when using Gaussian smoothing over nearby galaxies. NGC 4993 being a relatively nearby galaxy, at $sim 40 {rm Mpc}$ away, is subject to a significant effect of peculiar velocities. We demonstrate a direct dependence of the estimated peculiar velocity value on the choice of smoothing scale. We show that when not accounting for this systematic, a bias of $sim 200 {rm km s ^{-1}}$ in the peculiar velocity incurs a bias of $sim 4 {rm km s ^{-1} Mpc^{-1}}$ on the Hubble constant. We formulate a Bayesian model that accounts for the dependence of the peculiar velocity on the smoothing scale and by marginalising over this parameter we remove the need for a choice of smoothing scale. The proposed model yields $H_0 = 68.6 ^{+14.0}_{-8.5}~{rm km s^{-1} Mpc^{-1}}$. We demonstrate that under this model a more robust unbiased estimate of the Hubble constant from nearby GW sources is obtained.
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.
The detection of GW170817 in both gravitational waves and electromagnetic waves heralds the age of gravitational-wave multi-messenger astronomy. On 17 August 2017 the Advanced LIGO and Virgo detectors observed GW170817, a strong signal from the merger of a binary neutron-star system. Less than 2 seconds after the merger, a gamma-ray burst (GRB 170817A) was detected within a region of the sky consistent with the LIGO-Virgo-derived location of the gravitational-wave source. This sky region was subsequently observed by optical astronomy facilities, resulting in the identification of an optical transient signal within $sim 10$ arcsec of the galaxy NGC 4993. These multi-messenger observations allow us to use GW170817 as a standard siren, the gravitational-wave analog of an astronomical standard candle, to measure the Hubble constant. This quantity, which represents the local expansion rate of the Universe, sets the overall scale of the Universe and is of fundamental importance to cosmology. Our measurement combines the distance to the source inferred purely from the gravitational-wave signal with the recession velocity inferred from measurements of the redshift using electromagnetic data. This approach does not require any form of cosmic distance ladder; the gravitational wave analysis can be used to estimate the luminosity distance out to cosmological scales directly, without the use of intermediate astronomical distance measurements. We determine the Hubble constant to be $70.0^{+12.0}_{-8.0} , mathrm{km} , mathrm{s}^{-1} , mathrm{Mpc}^{-1}$ (maximum a posteriori and 68% credible interval). This is consistent with existing measurements, while being completely independent of them. Additional standard-siren measurements from future gravitational-wave sources will provide precision constraints of this important cosmological parameter.
The Hubble constant ($H_0$) estimated from the local Cepheid-supernova (SN) distance ladder is in 3-$sigma$ tension with the value extrapolated from cosmic microwave background (CMB) data assuming the standard cosmological model. Whether this tension represents new physics or systematic effects is the subject of intense debate. Here, we investigate how new, independent $H_0$ estimates can arbitrate this tension, assessing whether the measurements are consistent with being derived from the same model using the posterior predictive distribution (PPD). We show that, with existing data, the inverse distance ladder formed from BOSS baryon acoustic oscillation measurements and the Pantheon SN sample yields an $H_0$ posterior near-identical to the Planck CMB measurement. The observed local distance ladder value is a very unlikely draw from the resulting PPD. Turning to the future, we find that a sample of $sim50$ binary neutron star standard sirens (detectable within the next decade) will be able to adjudicate between the local and CMB estimates.
We investigate a recently proposed method for measuring the Hubble constant from gravitational wave detections of binary black hole coalescences without electromagnetic counterparts. In the absence of a direct redshift measurement, the missing information on the left-hand side of the Hubble-Lema^itre law is provided by the statistical knowledge on the redshift distribution of sources. We assume that source distribution in redshift depends on just one unknown hyper-parameter, modeling our ignorance of the astrophysical binary black hole distribution. With tens of thousands of these black sirens -- a realistic figure for the third generation detectors Einstein Telescope and Cosmic Explorer -- an observational constraint on the value of the Hubble parameter at percent level can be obtained. This method has the advantage of not relying on electromagnetic counterparts, which accompany a very small fraction of gravitational wave detections, nor on often unavailable or incomplete galaxy catalogs.