Multi-messenger constraints on the Hubble constant

Title: Multi-messenger constraints on the Hubble constant through combination of gravitational waves, gamma-ray bursts and kilonovae from neutron star mergers

Authors: Mattia Bulla, Michael W. Coughlin, Suhail Dhawan, Tim Dietrich

First Author’s Institution: The Oskar Klein Centre, Department of Astronomy, Stockholm University, AlbaNova, SE-10691 Stockholm, Sweden

Astrobite written by Alexandra Wernersson. Alexandra is a PhD student at the University of Amsterdam doing research in multi-messenger astrophysics.

Introduction

The Hubble constant, \(H_0\), measures the expansion rate of the universe. In our local universe \(H_0\) can be approximated by the famous linear equation

\(v_H = cz = H_0 D_L\) (1)

where \(D_L\) is the luminosity distance and \(v_H\) the Hubble flow velocity. The direct measurement of the Hubble constant, \(H_0\), is one the biggest questions in contemporary cosmology. As measurements have become more precise, a significant uncertainty in the \(H_0\) measurement have emerged due to the discrepancy of the value from late and early universe measurements. The cosmic microwave background gives us an estimate of \(H_0\) from the early universe which is in tension with late time measurements such as the standard candle Type Ia supernovae. Events producing gravitational waves (GWs) detectable by LIGO/Virgo have been proposed as promising independent measurements of \(H_0\). This is often called using GWs as ”standard sirens” and utilizes the fact that the GW amplitude strain, \(h\), is dependent on the luminosity distance \(D_L\). Furthermore, for GW binaries one measures an effective mass known as the chirp mass, \(M\), compromised of the two compact objects masses. Both the GW amplitude strain, \(h\), and the rate of change of the GW frequency, \(dfGW/dt\), are dependent on the chirp mass. Thus, the measurement of \(h\) and \(dfGW/dt\) provide an absolute measurement of \(D_L\).

The next puzzle piece is the redshift \(z\) which can be obtained either from the electromagnetic (EM) counterpart and its host galaxy, so called ”bright sirens”, or from a statistical analysis of galaxy clustering, ”dark sirens”. \(H_0\) is then obtained from Eq. (1) by the independent measurements of \(D_L\) and \(z\). The benefit of the above procedure is that the value of \(H_0\) is obtained without relying on the cosmic distance ladder or a specific cosmological model. The main uncertainty in measuring \(H_0\) from GWs lies in the difficulty of constraining the inclination angle \(i\) of the binary. The constraints on \(i\) can be improved by taking advantage of the viewing-angle, \(\theta_{obs}\), dependence on the EM signals. The authors of today’s paper review the current constraints on \(H_0\) using gamma-ray-bursts (GRBs) and kilonova (KN) observables from the GW170817 event and comment on the ability and future of this technique to solve the Hubble tension.

Inclination constraints of GRBs

Gamma ray bursts launch relativistic jets from the compact object formed after a black hole-neutron star (BH-NS) or binary neutron star (BNS) merger. The interaction of the jet with the surrounding environment produces an ”afterglow” consisting of synchrotron radiation of the accelerated electrons from the shocked medium in the X-ray and radio wavelengths. For a certain value of the jet angle \(\theta_j\) , the constraint on \(H_0\) is improved by a factor of two from the GW-only analysis. Furthermore, a jet moving close to the speed of light may appear to move with superluminal speed i.e. ”chasing” the emitted radiation along the line-of-sight of the observer. Therefore, a measurement of the apparent velocity can constrain the \(\theta_{obs}\) − \(\theta_j\) hence improving the constraint on \(H_0\). The superluminal motion of the jet associated with GW170817 in combination with the radio light-curves lead to an improvement of the degeneracy between \(\theta_{obs}\) and \(\theta_j\).

Inclination constraints from the kilonova

The kilonova event associated with BNS and BH-NS mergers is a thermal emission powered by the r-process in the neutron rich outflow of the merger. Because the outflows consists of different components with a variety of compositions and geometries the KN has a viewing angle dependence. This fact has been exploited in the case of GW170817 to put constraints on the inclination angle and thereby an improvement on the \(H_0\). One of the first analyses of this type is summarized by the pink curve in Figure 1 which lead to a 34% improvement on the uncertainties on \(H_0\) compared to only GW data. In contrast, a more recent study reduced the errors in \(H_0\) even further by the use of spectroscopy, which lead to an 54% improvement on the uncertainties on \(H_0\) showcased by the orange curve in Fig. 1.

Figure 1: Constraints on \(H_0\) from the BNS merger GW170817 using GW and EM data. The figure shows the marginalized 1D posterior density distributions for \(H_0\) with the inferred values labeled in the legend. Color scheme: Constraints from the GRB170817A with (cyan) or without (light green, red) including jet superluminal motion. Constraints from the KN AT 2017gfo from model fitting with broad-band photometry (pink) and spectroscopy (orange). The black curve correspond to the inferred value from the original standard siren approach. The inferred \(H_0\) value from Planck and SHOES are shown with the 1\(\sigma\) interval in brown and green respectively. Figure 2 in the paper.

Summary and outlook

Gravitational wave events are promising standard sirens for \(H_0\) measurements. In order to reach a ∼1% precision on \(H_0\) one would need about 50 − 200 GW events with an EM counterpart. Today’s paper has shown the constraints may be improved sooner and with fewer events if one improves the constraint on the viewing angle using the GRB afterglow, superluminal motion and KN. Doing so would significantly reduce the degeneracy between the distance and inclination angle of the GW data and improve the accuracy of the inferred \(H_0\) value.

Even so, in order to have sufficient accuracy on the inferred values from GW and EM events one still requires the systematics from the standard siren approach and the GRB/KN modeling to be under control. The systematic errors that may contribute include calibration errors that may affect the inferred value of \(D_L\). Further systematic uncertainties may arise from the models used to fit the GRB and KN data which may introduce a bias on the \(H_0\) value.

In conclusion, the coming decade will provide an increase in GW detectors and sensitivity. This will enable us to provide an understanding of the systematic uncertainties and showcase the potential of the standard siren approach. The combined effort of the multi-messenger community will have the potential to provide precise and accurate values of the Hubble constant.

Astrobite edited by Jason Hinkle & Ali Crisp

Featured image credit: National Science Foundation/LIGO/Sonoma State University/A. Simonnet

Author

Submit a Comment

Your email address will not be published. Required fields are marked *

Astrobites wants to learn more about you and how you use our site! Take the quick survey:

X