No Arabic abstract
We present the measurement of the Hubble Constant, $H_0$, with three strong gravitational lens systems. We describe a blind analysis of both PG1115+080 and HE0435-1223 as well as an extension of our previous analysis of RXJ1131-1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP 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 $H_{0}$. 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 PG1115+080, $H_{0}=70.1^{+5.3}_{-4.5}~rm km,s^{-1},Mpc^{-1}$ for HE0435-1223, and $H_{0}=77.0^{+4.0}_{-4.6}~rm km,s^{-1},Mpc^{-1}$ for RXJ1131-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 the above results assume a flat $Lambda$ cold dark matter cosmology with a uniform prior on $Omega_{textrm{m}}$ in [0.05, 0.5] and $H_{0}$ 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 $H_0$. The full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper.
Strongly lensed quasars can provide measurements of the Hubble constant ($H_{0}$) independent of any other methods. One of the key ingredients is exquisite high-resolution imaging data, such as Hubble Space Telescope (HST) imaging and adaptive-optics (AO) imaging from ground-based telescopes, which provide strong constraints on the mass distribution of the lensing galaxy. In this work, we expand on the previous analysis of three time-delay lenses with AO imaging (RXJ1131-1231, HE0435-1223, and PG1115+080), and perform a joint analysis of J0924+0219 by using AO imaging from the Keck Telescope, obtained as part of the SHARP (Strong lensing at High Angular Resolution Program) AO effort, with HST imaging to constrain the mass distribution of the lensing galaxy. Under the assumption of a flat $Lambda$CDM model with fixed $Omega_{rm m}=0.3$, we show that by marginalizing over two different kinds of mass models (power-law and composite models) and their transformed mass profiles via a mass-sheet transformation, we obtain $Delta t_{rm BA}hhat{sigma}_{v}^{-2}=6.89substack{+0.8-0.7}$ days, $Delta t_{rm CA}hhat{sigma}_{v}^{-2}=10.7substack{+1.6-1.2}$ days, and $Delta t_{rm DA}hhat{sigma}_{v}^{-2}=7.70substack{+1.0-0.9}$ days, where $h=H_{0}/100~rm km,s^{-1},Mpc^{-1}$ is the dimensionless Hubble constant and $hat{sigma}_{v}=sigma^{rm ob}_{v}/(280~rm km,s^{-1})$ is the scaled dimensionless velocity dispersion. Future measurements of time delays with 10% uncertainty and velocity dispersion with 5% uncertainty would yield a $H_0$ constraint of $sim15$% precision.
Accurate and precise measurements of the Hubble constant are critical for testing our current standard cosmological model and revealing possibly new physics. With Hubble Space Telescope (HST) imaging, each strong gravitational lens system with measured time delays can allow one to determine the Hubble constant with an uncertainty of $sim 7%$. Since HST will not last forever, we explore adaptive-optics (AO) imaging as an alternative that can provide higher angular resolution than HST imaging but has a less stable point spread function (PSF) due to atmospheric distortion. To make AO imaging useful for time-delay-lens cosmography, we develop a method to extract the unknown PSF directly from the imaging of strongly lensed quasars. In a blind test with two mock data sets created with different PSFs, we are able to recover the important cosmological parameters (time-delay distance, external shear, lens mass profile slope, and total Einstein radius). Our analysis of the Keck AO image of the strong lens system RXJ1131-1231 shows that the important parameters for cosmography agree with those based on HST imaging and modeling within 1-$sigma$ uncertainties. Most importantly, the constraint on the model time-delay distance by using AO imaging with $0.045$resolution is tighter by $sim 50%$ than the constraint of time-delay distance by using HST imaging with $0.09$when a power-law mass distribution for the lens system is adopted. Our PSF reconstruction technique is generic and applicable to data sets that have multiple nearby point sources, enabling scientific studies that require high-precision models of the PSF.
Time-delay strong lensing provides a unique way to directly measure the Hubble constant ($H_{0}$). The precision of the $H_{0}$ measurement depends on the uncertainties in the time-delay measurements, the mass distribution of the main deflector(s), and the mass distribution along the line of sight. Tie and Kochanek (2018) have proposed a new microlensing effect on time delays based on differential magnification of the coherent accretion disc variability of the lensed quasar. If real, this effect could significantly broaden the uncertainty on the time delay measurements by up to $30%$ for lens systems such as PG1115+080, which have relatively short time delays and monitoring over several different epochs. In this paper we develop a new technique that uses the time-delay ratios and simulated microlensing maps within a Bayesian framework in order to limit the allowed combinations of microlensing delays and thus to lessen the uncertainties due to the proposed effect. We show that, under the assumption of Tie and Kochanek (2018), the uncertainty on the time-delay distance ($D_{Delta t}$, which is proportional to 1/$H_{0}$) of short time-delay ($sim18$ days) lens, PG1115+080, increases from $sim7%$ to $sim10%$ by simultaneously fitting the three time-delay measurements from the three different datasets across twenty years, while in the case of long time-delay ($sim90$ days) lens, the microlensing effect on time delays is negligible as the uncertainty on $D_{Delta t}$ of RXJ1131-1231 only increases from $sim2.5%$ to $sim2.6%$.
We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $Lambda$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1sigma$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3sigma$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $Lambda$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $Lambda$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $sim73$-$78~mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints.
In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated datasets. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ($ |A| < 2%$) and precision ($< 6%$ per system), even in the presence of poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix, and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges.