skip to main content

This content will become publicly available on June 22, 2024

Title: The Hitchhiker’s Guide to the Galaxy Catalog Approach for Dark Siren Gravitational-wave Cosmology
Abstract We outline the “dark siren” galaxy catalog method for cosmological inference using gravitational wave (GW) standard sirens, clarifying some common misconceptions in the implementation of this method. When a confident transient electromagnetic counterpart to a GW event is unavailable, the identification of a unique host galaxy is in general challenging. Instead, as originally proposed by Schutz, one can consult a galaxy catalog and implement a dark siren statistical approach incorporating all potential host galaxies within the localization volume. Trott & Huterer recently claimed that this approach results in a biased estimate of the Hubble constant, H 0 , when implemented on mock data, even if optimistic assumptions are made. We demonstrate explicitly that, as previously shown by multiple independent groups, the dark siren statistical method leads to an unbiased posterior when the method is applied to the data correctly. We highlight common sources of error possible to make in the generation of mock data and implementation of the statistical framework, including the mismodeling of selection effects and inconsistent implementations of the Bayesian framework, which can lead to a spurious bias.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
The Astronomical Journal
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Gravitational-wave (GW) radiation from a coalescing compact binary is a standard siren, as the luminosity distance of each event can be directly measured from the amplitude of the signal. One possibility to constrain cosmology using the GW siren is to perform statistical inference on a population of binary black hole (BBH) events. In essence, this statistical method can be viewed as follows. We can modify the shape of the distribution of observed BBH events by changing the cosmological parameters until it eventually matches the distribution constructed from an astrophysical population model, thereby allowing us to determine the cosmological parameters. In this work, we derive the Cramér–Rao bound for both cosmological parameters and those governing the astrophysical population model from this statistical dark siren method by examining the Fisher information contained in the event distribution. Our study provides analytical insights and enables fast yet accurate estimations of the statistical accuracy of dark siren cosmology. Furthermore, we consider the bias in cosmology due to unmodeled substructures in the merger rate and mass distribution. We find that a 1% deviation in the astrophysical model can lead to a more than 1% error in the Hubble constant. This could limit the accuracy of dark siren cosmology when there are more than 104BBH events detected.

    more » « less

    The current and next observation seasons will detect hundreds of gravitational waves (GWs) from compact binary systems coalescence at cosmological distances. When combined with independent electromagnetic measurements, the source redshift will be known, and we will be able to obtain precise measurements of the Hubble constant H0 via the distance–redshift relation. However, most observed mergers are not expected to have electromagnetic counterparts, which prevents a direct redshift measurement. In this scenario, one possibility is to use the dark sirens method that statistically marginalizes over all the potential host galaxies within the GW location volume to provide a probabilistic source redshift. Here we presented H0 measurements using two new dark sirens compared to previous analyses using DECam data: GW190924$\_$021846 and GW200202$\_$154313. The photometric redshifts of the possible host galaxies of these two events are acquired from the DECam Local Volume Exploration Survey (DELVE) carried out on the Blanco telescope at Cerro Tololo. The combination of the H0 posterior from GW190924$\_$021846 and GW200202$\_$154313 together with the bright siren GW170817 leads to $H_{0} = 68.84^{+15.51}_{-7.74}\, \rm {km\, s^{-1}\, Mpc^{-1}}$. Including these two dark sirens improves the 68  per cent confidence interval (CI) by 7  per cent over GW170817 alone. This demonstrates that the addition of well-localized dark sirens in such analysis improves the precision of cosmological measurements. Using a sample containing 10 well-localized dark sirens observed during the third LIGO/Virgo observation run, without the inclusion of GW170817, we determine a measurement of $H_{0} = 76.00^{+17.64}_{-13.45}\, \rm {km\, s^{-1}\, Mpc^{-1}}$.

    more » « less
  3. Abstract

    We present a measurement of the Hubble ConstantH0using the gravitational wave event GW190412, an asymmetric binary black hole merger detected by LIGO/Virgo, as a dark standard siren. This event does not have an electromagnetic counterpart, so we use the statistical standard siren method and marginalize over potential host galaxies from the Dark Energy Spectroscopic Instrument (DESI) survey. GW190412 is well-localized to 12 deg2(90% credible interval), so it is promising for a dark siren analysis. The dark siren value forH0=85.433.9+29.1km s−1 Mpc−1, with a posterior shape that is consistent with redshift overdensities. When combined with the bright standard siren measurement from GW170817 we recoverH0=77.965.03+23.0km s−1 Mpc−1, consistent with both early and late-time Universe measurements ofH0. This work represents the first standard siren analysis performed with DESI data, and includes the most complete spectroscopic sample used in a dark siren analysis to date.

    more » « less
  4. Abstract We use 47 gravitational wave sources from the Third LIGO–Virgo–Kamioka Gravitational Wave Detector Gravitational Wave Transient Catalog (GWTC–3) to estimate the Hubble parameter H ( z ), including its current value, the Hubble constant H 0 . Each gravitational wave (GW) signal provides the luminosity distance to the source, and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog. Using the binary black hole (BBH) redshifted masses, we simultaneously infer the source mass distribution and H ( z ). The source mass distribution displays a peak around 34 M ⊙ , followed by a drop-off. Assuming this mass scale does not evolve with the redshift results in a H ( z ) measurement, yielding H 0 = 68 − 8 + 12 km s − 1 Mpc − 1 (68% credible interval) when combined with the H 0 measurement from GW170817 and its electromagnetic counterpart. This represents an improvement of 17% with respect to the H 0 estimate from GWTC–1. The second method associates each GW event with its probable host galaxy in the catalog GLADE+ , statistically marginalizing over the redshifts of each event’s potential hosts. Assuming a fixed BBH population, we estimate a value of H 0 = 68 − 6 + 8 km s − 1 Mpc − 1 with the galaxy catalog method, an improvement of 42% with respect to our GWTC–1 result and 20% with respect to recent H 0 studies using GWTC–2 events. However, we show that this result is strongly impacted by assumptions about the BBH source mass distribution; the only event which is not strongly impacted by such assumptions (and is thus informative about H 0 ) is the well-localized event GW190814. 
    more » « less

    Accurate measurement of the Hubble constant from standard sirens such as the gravitational wave (GW) sources with electromagnetic counterparts relies on the robust peculiar velocity correction of the redshift of the host galaxy. We show in this work that the peculiar velocity of the host galaxies exhibits a correlation with the properties of the host galaxy primarily such as its stellar mass and this correlation also evolves with redshift. As the galaxies of higher stellar mass tend to form in galaxies with higher halo masses which are located in spatial regions having a non-linear fluctuation in the density field of the matter distribution, the root mean square peculiar velocity of more massive galaxies is higher. As a result, depending on the formation channel of the binary compact objects, the peculiar velocity contamination to the galaxies will be different. The variation in the peculiar velocity of the host galaxies can lead to a significant variation in the estimation of the Hubble constant inferred using sources such as binary neutron stars. For the network of GW detectors such as LIGO-Virgo-KAGRA (LVK), LVK+LIGO-India, and Cosmic Explorer + Einstein Telescope, the variation in the precision of Hubble constant inferred from 10 bright siren events can vary from $\sim 5.4 - 6~{{\ \rm per \, cent}}$, $\sim 4.5 - 5.3~{{\ \rm per \, cent}}$, and $\sim 1.1 - 2.7~{{\ \rm per \, cent}}$, respectively. The impact of such a correlation between peculiar velocity and stellar mass on the inference of the Hubble constant is not only limited to GW sources but also applicable to type-Ia supernovae.

    more » « less