We present a measurement of the Hubble constant (H0) 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 ΛCDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}~\mathrm{km~s^{-1}~Mpc^{-1}}$, a $2.4{{\ \rm per\ cent}}$ precision measurement, in agreement with local measurements of H0 from type Ia supernovae calibrated by the distance ladder, but in 3.1σ 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.3σ tension with Planck CMB determinations of H0 in flat Λ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 Λ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 H0 inference to cosmological model assumptions. For six different cosmological models, our combined inference on H0 ranges from ∼73–78 km s−1 Mpc−1, which is consistent with the local distance ladder constraints.

Strong gravitational lenses with measured time delays between the multiple images and models of the lens mass distribution allow a one-step determination of the time-delay distance, and thus a measure of cosmological parameters. We present a blind analysis of the gravitational lens RXJ1131−1231 incorporating (1) the newly measured time delays from COSMOGRAIL, the COSmological MOnitoring of GRAvItational Lenses, (2) archival Hubble Space Telescope imaging of the lens system, (3) a new velocity-dispersion measurement of the lens galaxy of 323 ± 20 km s −1 based on Keck spectroscopy, and (4) a characterization of the line-of-sight structures via observations of the lens' environment and ray tracing through the Millennium Simulation. Our blind analysis is designed to prevent experimenter bias. The joint analysis of the data sets allows a time-delay distance measurement to 6% precision that takes into account all known systematic uncertainties. Time-delay lenses constrain especially tightly the Hubble constant H 0 (5.7% and 4.0% respectively in wCDM and open ΛCDM) and curvature of the universe. The overall information content is similar to that of Baryon Acoustic Oscillation experiments. Thus, they complement well other cosmological probes, and provide an independent check of unknown systematics. Our measurement of the Hubble constant is completely independent of those based on the local distance ladder method, providing an important consistency check of the standard cosmological model and of general relativity.

We present a new measurement of the Hubble Constant H 0 and other cosmological parameters based on the joint analysis of three multiply-imaged quasar systems with measured gravitational time delays. First, we measure the time delay of HE 0435−1223 from 13-year light curves obtained as part of the COSMOGRAIL project. Companion papers detail the modeling of the main deflectors and line of sight effects, and how these data are combined to determine the time-delay distance of HE 0435−1223. Crucially, the measurements are carried out blindly with respect to cosmological parameters in order to avoid confirmation bias. We then combine the timedelay distance of HE 0435−1223 with previous measurements from systems B1608+656 and RXJ1131−1231 to create a Time Delay Strong Lensing probe (TDSL). In flat ΛCDM with free matter and energy density, we find H 0 = 71.9 +2.4 −3.0 km s −1 Mpc −1 and Ω Λ = 0.62 +0.24 −0.35 . This measurement is completely independent of, and in agreement with, the local distance ladder measurements of H 0 . We explore more general cosmological models combining TDSL with other probes, illustrating its power to break degeneracies inherent to other methods. The joint constraints from TDSL and Planck are H 0 = 69.2 +1.4 −2.2 km s −1 Mpc −1 , Ω Λ = 0.70 +0.01 −0.01 and Ω k = 0.003 +0.004 −0.006 in open ΛCDM and H 0 = 79.0 +4.4 −4.2 km s −1 Mpc −1 , Ω de = 0.77 +0.02 −0.03 and w = −1.38 +0.14 −0.16 in flat w CDM. In combination with Planck and Baryon Acoustic Oscillation data, when relaxing the constraints on the numbers of relativistic species we find N eff = 3.34 +0.21 −0.21 in N eff ΛCDM and when relaxing the total mass of neutrinos we find Σm ν ≤ 0.182 eV in m ν ΛCDM. Finally, in an open w CDM in combination with Planck and CMB lensing we find H 0 = 77.9 +5.0 −4.2 km s −1 Mpc −1 , Ω de = 0.77 +0.03 −0.03 , Ω k = −0.003 +0.004 −0.004 and w = −1.37 +0.18 −0.23 .

Under the assumption of a flat ΛCDM cosmology, recent data from the Planck satellite point toward a Hubble constant that is in tension with that measured by gravitational lens time delays and by the local distance ladder. Prosaically, this difference could arise from unknown systematic uncertainties in some of the measurements. More interestingly -if systematics were ruled out -resolving the tension would require a departure from the flat ΛCDM cosmology, introducing for example a modest amount of spatial curvature, or a non-trivial dark energy equation of state. To begin to address these issues, we present here an analysis of the gravitational lens RXJ1131−1231 that is improved in one particular regard: we examine the issue of systematic error introduced by an assumed lens model density profile. We use more flexible gravitational lens models with baryonic and dark matter components, and find that the exquisite Hubble Space Telescope image with thousands of intensity pixels in the Einstein ring and the stellar velocity dispersion of the lens contain sufficient information to constrain these more flexible models. The total uncertainty on the time-delay distance is 6.6% for a single system. We proceed to combine our improved time-delay distance measurements with the WMAP9 and Planck posteriors. In an open ΛCDM model, the data for RXJ1131−1231 in combination with Planck favor a flat universe with Ω k = 0.00 +0.01 −0.02 (68% CI). In a flat wCDM model, the combination of RXJ1131−1231 and Planck yields w = −1.52 +0.19 −0.20 (68% CI).

We present a blind time-delay strong lensing (TDSL) cosmographic analysis of the doubly imaged quasar SDSS 1206+4332 . We combine the relative time delay between the quasar images, Hubble Space Telescope imaging, the Keck stellar velocity dispersion of the lensing galaxy, and wide-field photometric and spectroscopic data of the field to constrain two angular diameter distance relations. The combined analysis is performed by forward modelling the individual data sets through a Bayesian hierarchical framework, and it is kept blind until the very end to prevent experimenter bias. After unblinding, the inferred distances imply a Hubble constant H 0 = 68.8 +5.4 −5.1 km s −1 Mpc −1 , assuming a flat Λ cold dark matter cosmology with uniform prior on Ω m in [0.05, 0.5]. The precision of our cosmographic measurement with the doubly imaged quasar SDSS 1206+4332 is comparable with those of quadruply imaged quasars and opens the path to perform on selected doubles the same analysis as anticipated for quads. Our analysis is based on a completely independent lensing code than our previous three H0LiCOW systems and the new measurement is fully consistent with those. We provide the analysis scripts paired with the publicly available software to facilitate independent analysis. The consistency between blind measurements with independent codes provides an important sanity check on lens modelling systematics. By combining the likelihoods of the four systems under the same prior, we obtain H 0 = 72.5 +2.1 −2.3 km s −1 Mpc −1 . This measurement is independent of the distance ladder and other cosmological probes.

Strong gravitational lens systems with time delays between the multiple images allow measurements of time-delay distances, which are primarily sensitive to the Hubble constant that is key to probing dark energy, neutrino physics, and the spatial curvature of the Universe, as well as discovering new physics. We present H0LiCOW (H 0 Lenses in COSMOGRAIL's Wellspring), a program that aims to measure H 0 with < 3.5% uncertainty from five lens systems (B1608+656, RXJ1131−1231, HE 0435−1223, WFI2033−4723 and HE 1104−1805). We have been acquiring (1) time delays through COSMOGRAIL and Very Large Array monitoring, (2) high-resolution Hubble Space Telescope imaging for the lens mass modeling, (3) wide-field imaging and spectroscopy to characterize the lens environment, and (4) moderate-resolution spectroscopy to obtain the stellar velocity dispersion of the lenses for mass modeling. In cosmological models with one-parameter extension to flat ΛCDM, we expect to measure H 0 to < 3.5% in most models, spatial curvature Ω k to 0.004, w to 0.14, and the effective number of neutrino species to 0.2 (1σ uncertainties) when combined with current CMB experiments. These are, respectively, a factor of ∼ 15, ∼ 2, and ∼ 1.5 tighter than CMB alone. Our data set will further enable us to study the stellar initial mass function of the lens galaxies, and the co-evolution of supermassive black holes and their host galaxies. This program will provide a foundation for extracting cosmological distances from the hundreds of time-delay lenses that are expected to be discovered in current and future surveys.

A new method for improving the resolution of astronomical images is presented. It is based on the principle that sampled data cannot be fully deconvolved without violating the sampling theorem. Thus, the sampled image should not be deconvolved by the total Point Spread Function, but by a narrower function chosen so that the resolution of the deconvolved image is compatible with the adopted sampling.Our deconvolution method gives results which are, in at least some cases, superior to those of other commonly used techniques : in particular, it does not produce ringing around point sources superimposed on a smooth background. Moreover, it allows to perform accurate astrometry and photometry of crowded fields. These improvements are a consequence of both the correct treatment of sampling and the recognition that the most probable astronomical image is not a flat one.The method is also well adapted to the optimal combination of different images of the same object, as can be obtained, e.g., from infrared observations or via adaptive optics techniques.

Measuring time delays between the multiple images of gravitationally lensed quasars is now recognized as a competitive way to constrain the cosmological parameters, and it is complementary with other cosmological probes. This requires long and well sampled optical light curves of numerous lensed quasars, such as those obtained by the COSMOGRAIL collaboration. High-quality data from our monitoring campaign call for novel numerical techniques to robustly measure the delays, as well as the associated random and systematic uncertainties, even in the presence of microlensing variations. We propose three different point estimators to measure time delays, which are explicitly designed to handle light curves with extrinsic variability. These methods share a common formalism, which enables them to process data from n-image lenses. Since the estimators rely on significantly contrasting ideas, we expect them to be sensitive to different bias sources. For each method and data set, we empirically estimate both the precision and accuracy (bias) of the time delay measurement using simulated light curves with known time delays that closely mimic the observations. Finally, we test the self-consistency of our approach, and we demonstrate that our bias estimation is serviceable. These new methods, including the empirical uncertainty estimator, will represent the standard benchmark for analyzing the COSMOGRAIL light curves.