skip to main content


Title: TDCOSMO: I. An exploration of systematic uncertainties in the inference of H 0 from time-delay cosmography
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
Award ID(s):
1906976 1907396
PAR ID:
10175672
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; « less
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
639
ISSN:
0004-6361
Page Range / eLocation ID:
A101
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. null (Ed.)
    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
  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. Time delay cosmography uses the arrival time delays between images in strong gravitational lenses to measure cosmological parameters, in particular the Hubble constant H 0 . The lens models used in time delay cosmography omit dark matter subhalos and line-of-sight halos because their effects are assumed to be negligible. We explicitly quantify this assumption by analyzing mock lens systems that include full populations of dark matter subhalos and line-of-sight halos, applying the same modeling assumptions used in the literature to infer H 0 . We base the mock lenses on six quadruply imaged quasars that have delivered measurements of the Hubble constant, and quantify the additional uncertainties and/or bias on a lens-by-lens basis. We show that omitting dark substructure does not bias inferences of H 0 . However, perturbations from substructure contribute an additional source of random uncertainty in the inferred value of H 0 that scales as the square root of the lensing volume divided by the longest time delay. This additional source of uncertainty, for which we provide a fitting function, ranges from 0.7 − 2.4%. It may need to be incorporated in the error budget as the precision of cosmographic inferences from single lenses improves, and it sets a precision limit on inferences from single lenses. 
    more » « less
  5. ABSTRACT

    Approximate methods to populate dark-matter haloes with galaxies are of great utility to galaxy surveys. However, the limitations of simple halo occupation models (HODs) preclude a full use of small-scale galaxy clustering data and call for more sophisticated models. We study two galaxy populations, luminous red galaxies (LRGs) and star-forming emission-line galaxies (ELGs), at two epochs, z = 1 and z = 0, in the large-volume, high-resolution hydrodynamical simulation of the MillenniumTNG project. In a partner study we concentrated on the small-scale, one-halo regime down to r ∼ 0.1 h−1 Mpc, while here we focus on modelling galaxy assembly bias in the two-halo regime, r ≳ 1 h−1 Mpc. Interestingly, the ELG signal exhibits scale dependence out to relatively large scales (r ∼ 20 h−1 Mpc), implying that the linear bias approximation for this tracer is invalid on these scales, contrary to common assumptions. The 10–15 per cent discrepancy is only reconciled when we augment our halo occupation model with a dependence on extrinsic halo properties (‘shear’ being the best-performing one) rather than intrinsic ones (e.g. concentration, peak mass). We argue that this fact constitutes evidence for two-halo galaxy conformity. Including tertiary assembly bias (i.e. a property beyond mass and ‘shear’) is not an essential requirement for reconciling the galaxy assembly bias signal of LRGs, but the combination of external and internal properties is beneficial for recovering ELG the clustering. We find that centrals in low-mass haloes dominate the assembly bias signal of both populations. Finally, we explore the predictions of our model for higher order statistics such as nearest neighbour counts. The latter supplies additional information about galaxy assembly bias and can be used to break degeneracies between halo model parameters.

     
    more » « less