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, more » 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. « less
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; « less
Award ID(s):
1906976 1907396
Publication Date:
Journal Name:
Astronomy & Astrophysics
Page Range or eLocation-ID:
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 velocitymore »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.« less
  2. 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 themore »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).« less
  3. 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 ofmore »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).« less
  4. 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 inmore »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.« less
  5. ABSTRACT Core formation and runaway core collapse in models with self-interacting dark matter (SIDM) significantly alter the central density profiles of collapsed haloes. Using a forward modelling inference framework with simulated data-sets, we demonstrate that flux ratios in quadruple image strong gravitational lenses can detect the unique structural properties of SIDM haloes, and statistically constrain the amplitude and velocity dependence of the interaction cross-section in haloes with masses between 106 and 1010 M⊙. Measurements on these scales probe self-interactions at velocities below $30 \ \rm {km} \ \rm {s^{-1}}$, a relatively unexplored regime of parameter space, complimenting constraints at highermore »velocities from galaxies and clusters. We cast constraints on the amplitude and velocity dependence of the interaction cross-section in terms of σ20, the cross-section amplitude at $20 \ \rm {km} \ \rm {s^{-1}}$. With 50 lenses, a sample size available in the near future, and flux ratios measured from spatially compact mid-IR emission around the background quasar, we forecast $\sigma _{20} \lt 11\rm {\small {--}}23 \ \rm {cm^2} \rm {g^{-1}}$ at $95 {{\ \rm per\ cent}}$ CI, depending on the amplitude of the subhalo mass function, and assuming cold dark matter (CDM). Alternatively, if $\sigma _{20} = 19.2 \ \rm {cm^2}\rm {g^{-1}}$ we can rule out CDM with a likelihood ratio of 20:1, assuming an amplitude of the subhalo mass function that results from doubly efficient tidal disruption in the Milky Way relative to massive elliptical galaxies. These results demonstrate that strong lensing of compact, unresolved sources can constrain SIDM structure on sub-galactic scales across cosmological distances, and the evolution of SIDM density profiles over several Gyr of cosmic time.« less