skip to main content


Title: TDCOSMO: IV. Hierarchical time-delay cosmography – joint inference of the Hubble constant and galaxy density profiles
The H0LiCOW collaboration inferred via strong gravitational lensing time delays a Hubble constant value of H 0 = 73.3 −1.8 +1.7 km s −1 Mpc −1 , describing deflector mass density profiles by either a power-law or stars (constant mass-to-light ratio) plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in H 0 . We quantify any potential effect of the MST with a flexible family of mass models, which directly encodes it, and they are hence maximally degenerate with H 0 . Our calculation is based on a new hierarchical Bayesian approach in which the MST is only constrained by stellar kinematics. The approach is validated on mock lenses, which are generated from hydrodynamic simulations. We first applied the inference to the TDCOSMO sample of seven lenses, six of which are from H0LiCOW, and measured H 0 = 74.5 −6.1 +5.6 km s −1 Mpc −1 . Secondly, in order to further constrain the deflector mass density profiles, we added imaging and spectroscopy for a set of 33 strong gravitational lenses from the Sloan Lens ACS (SLACS) sample. For nine of the 33 SLAC lenses, we used resolved kinematics to constrain the stellar anisotropy. From the joint hierarchical analysis of the TDCOSMO+SLACS sample, we measured H 0 = 67.4 −3.2 +4.1 km s −1 Mpc −1 . This measurement assumes that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. Without relying on the form of the mass density profile used by H0LiCOW, we achieve a ∼5% measurement of H 0 . While our new hierarchical analysis does not statistically invalidate the mass profile assumptions by H0LiCOW – and thus the H 0 measurement relying on them – it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on H 0 derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data.  more » « less
Award ID(s):
1906976 1907396
NSF-PAR ID:
10287447
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
643
ISSN:
0004-6361
Page Range / eLocation ID:
A165
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, H 0 . As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on H 0 from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in H 0 values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the “bulge-halo” conspiracy, H 0 is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields H 0 = 74.0 −1.8 +1.7 km s −1 Mpc −1 , while the power-law model yields 74.2 −1.6 +1.6 km s −1 Mpc −1 . In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO. 
    more » « less
  2. Strong-lensing time delays enable the measurement of the Hubble constant ( H 0 ) independently of other traditional methods. The main limitation to the precision of time-delay cosmography is mass-sheet degeneracy (MSD). Some of the previous TDCOSMO analyses broke the MSD by making standard assumptions about the mass density profile of the lens galaxy, reaching 2% precision from seven lenses. However, this approach could potentially bias the H 0 measurement or underestimate the errors. For this work, we broke the MSD for the first time using spatially resolved kinematics of the lens galaxy in RXJ1131−1231 obtained from the Keck Cosmic Web Imager spectroscopy, in combination with previously published time delay and lens models derived from Hubble Space Telescope imaging. This approach allowed us to robustly estimate H 0 , effectively implementing a maximally flexible mass model. Following a blind analysis, we estimated the angular diameter distance to the lens galaxy D d  = 865 −81 +85 Mpc and the time-delay distance D Δt  = 2180 −271 +472 Mpc, giving H 0  = 77.1 −7.1 +7.3 km s −1 Mpc −1 – for a flat Λ cold dark matter cosmology. The error budget accounts for all uncertainties, including the MSD inherent to the lens mass profile and line-of-sight effects, and those related to the mass–anisotropy degeneracy and projection effects. Our new measurement is in excellent agreement with those obtained in the past using standard simply parametrized mass profiles for this single system ( H 0  = 78.3 −3.3 +3.4 km s −1 Mpc −1 ) and for seven lenses ( H 0  = 74.2 −1.6 +1.6 km s −1 Mpc −1 ), or for seven lenses using single-aperture kinematics and the same maximally flexible models used by us ( H 0  = 73.3 −5.8 +5.8 km s −1 Mpc −1 ). This agreement corroborates the methodology of time-delay cosmography. 
    more » « less
  3. null (Ed.)
    Strong lensing time delays can measure the Hubble constant H 0 independently of any other probe. Assuming commonly used forms for the radial mass density profile of the lenses, a 2% precision has been achieved with seven Time-Delay Cosmography (TDCOSMO) lenses, in tension with the H 0 from the cosmic microwave background. However, without assumptions on the radial mass density profile – and relying exclusively on stellar kinematics to break the mass-sheet degeneracy – the precision drops to 8% with the current data obtained using the seven TDCOSMO lenses, which is insufficient to resolve the H 0 tension. With the addition of external information from 33 Sloan Lens ACS (SLACS) lenses, the precision improves to 5% if the deflectors of TDCOSMO and SLACS lenses are drawn from the same population. We investigate the prospect of improving the precision of time-delay cosmography without relying on mass profile assumptions to break the mass-sheet degeneracy. Our forecasts are based on a previously published hierarchical framework. With existing samples and technology, 3.3% precision on H 0 can be reached by adding spatially resolved kinematics of the seven TDCOSMO lenses. The precision improves to 2.5% with the further addition of kinematics for 50 nontime-delay lenses from SLACS and the Strong Lensing Legacy Survey. Expanding the samples to 40 time-delay and 200 nontime-delay lenses will improve the precision to 1.5% and 1.2%, respectively. Time-delay cosmography can reach sufficient precision to resolve the Hubble tension at 3–5 σ , without assumptions on the radial mass profile of lens galaxies. By obtaining this precision with and without external datasets, we will test the consistency of the samples and enable further improvements based on even larger future samples of time-delay and nontime-delay lenses (e.g., from the Rubin , Euclid , and Roman Observatories). 
    more » « less
  4. Abstract We present the lens mass model of the quadruply-imaged gravitationally lensed quasar WFI2033−4723, and perform a blind cosmographical analysis based on this system. Our analysis combines (1) time-delay measurements from 14 years of data obtained by the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL) collaboration, (2) high-resolution Hubble Space Telescope imaging, (3) a measurement of the velocity dispersion of the lens galaxy based on ESO-MUSE data, and (4) multi-band, wide-field imaging and spectroscopy characterizing the lens environment. We account for all known sources of systematics, including the influence of nearby perturbers and complex line-of-sight structure, as well as the parametrization of the light and mass profiles of the lensing galaxy. After unblinding, we determine the effective time-delay distance to be $4784_{-248}^{+399}~\mathrm{Mpc}$, an average precision of $6.6{{\ \rm per\ cent}}$. This translates to a Hubble constant $H_{0} = 71.6_{-4.9}^{+3.8}~\mathrm{km~s^{-1}~Mpc^{-1}}$, assuming a flat ΛCDM cosmology with a uniform prior on Ωm in the range [0.05, 0.5]. This work is part of the H0 Lenses in COSMOGRAIL’s Wellspring (H0LiCOW) collaboration, and the full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper (H0LiCOW XIII). 
    more » « less
  5. ABSTRACT We present the measurement of the Hubble constant, H0, with three strong gravitational lens systems. We describe a blind analysis of both PG 1115+080 and HE 0435−1223 as well as an extension of our previous analysis of RXJ 1131−1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP (Strong-lensing High Angular Resolution Programme) AO effort, with Hubble Space Telescope (HST) imaging, velocity dispersion measurements, and a description of the line-of-sight mass distribution to build an accurate and precise lens mass model. This mass model is then combined with the COSMOGRAIL-measured time delays in these systems to determine H0. We do both an AO-only and an AO + HST analysis of the systems and find that AO and HST results are consistent. After unblinding, the AO-only analysis gives $H_{0}=82.8^{+9.4}_{-8.3}~\rm km\, s^{-1}\, Mpc^{-1}$ for PG 1115+080, $H_{0}=70.1^{+5.3}_{-4.5}~\rm km\, s^{-1}\, Mpc^{-1}$ for HE 0435−1223, and $H_{0}=77.0^{+4.0}_{-4.6}~\rm km\, s^{-1}\, Mpc^{-1}$ for RXJ 1131−1231. The joint AO-only result for the three lenses is $H_{0}=75.6^{+3.2}_{-3.3}~\rm km\, s^{-1}\, Mpc^{-1}$. The joint result of the AO + HST analysis for the three lenses is $H_{0}=76.8^{+2.6}_{-2.6}~\rm km\, s^{-1}\, Mpc^{-1}$. All of these results assume a flat Λ cold dark matter cosmology with a uniform prior on Ωm in [0.05, 0.5] and H0 in [0, 150] $\rm km\, s^{-1}\, Mpc^{-1}$. This work is a collaboration of the SHARP and H0LiCOW teams, and shows that AO data can be used as the high-resolution imaging component in lens-based measurements of H0. The full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper. 
    more » « less