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 .

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.

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.

Based on a new sample of 355 quasars with significant optical polarization and using complementary statistical methods, we confirm that quasar polarization vectors are not randomly oriented over the sky with a probability often in excess of 99.9%. The polarization vectors appear coherently oriented or aligned over huge (∼1 Gpc) regions of the sky located at both low (z ∼ 0.5) and high (z ∼ 1.5) redshifts and characterized by different preferred directions of the quasar polarization. In fact, there seems to exist a regular alternance along the line of sight of regions of randomly and aligned polarization vectors with a typical comoving length scale of 1.5 Gpc. Furthermore, the mean polarization angleθ appears to rotate with redshift at the rate of ∼30• per Gpc. The symmetry of the theθ − z relation is mirror-like, the mean polarization angle rotating clockwise with increasing redshift in North Galactic hemisphere and counter-clockwise in the South Galactic one. These characteristics make the alignment effect difficult to explain in terms of local mechanisms, namely a contamination by interstellar polarization in our Galaxy. While interpretations like a global rotation of the Universe can potentially explain the effect, the properties we observe qualitatively correspond to the dichroism and birefringence predicted by photon-pseudoscalar oscillation within a magnetic field. Interestingly, the alignment effect seems to be prominent along an axis not far from preferred directions tentatively identified in the Cosmic Microwave Background maps. Although many questions and more particularly the interpretation of the effect remain open, alignments of quasar polarization vectors appear as a promising new way to probe the Universe and its dark components at extremely large scales.

The light travel time differences in strong gravitational lensing systems allows an independent determination of the Hubble constant. This method has been successfully applied to several lens systems. The formally most precise measurements are, however, in tension with the recent determination of H 0 from the Planck satellite for a spatially flat six-parameters ΛCDM cosmology. We reconsider the uncertainties of the method, concerning the mass profile of the lens galaxies, and show that the formal precision relies on the assumption that the mass profile is a perfect power law. Simple analytical arguments and numerical experiments reveal that masssheet like transformations yield significant freedom in choosing the mass profile, even when exquisite Einstein rings are observed. Furthermore, the characterization of the environment of the lens does not break that degeneracy which is not physically linked to extrinsic convergence. We present an illustrative example where the multiple imaging properties of a composite (baryons + dark matter) lens can be extremely well reproduced by a power-law model having the same velocity dispersion, but with predictions for the Hubble constant that deviate by ∼20%. Hence we conclude that the impact of degeneracies between parametrized models have been underestimated in current H 0 measurements from lensing, and need to be carefully reconsidered.

When an image of a strongly lensed quasar is microlensed, the different components of its spectrum are expected to be differentially magnified owing to the different sizes of the corresponding emitting region. Chromatic changes are expected to be observed in the continuum while the emission lines should be deformed as a function of the size, geometry and kinematics of the regions from which they originate. Microlensing of the emission lines has been reported only in a handful of systems so far. In this paper we search for microlensing deformations of the optical spectra of pairs of images in 17 lensed quasars with bolometric luminosities between 10 44.7−47.4 erg/s and black hole masses 10 7.6−9.8 M . This sample is composed of 13 pairs of previously unpublished spectra and four pairs of spectra from literature. Our analysis is based on a simple spectral decomposition technique which allows us to isolate the microlensed fraction of the flux independently of a detailed modeling of the quasar emission lines. Using this technique, we detect microlensing of the continuum in 85% of the systems. Among them, 80% show microlensing of the broad emission lines. Focusing on the most common emission lines in our spectra (C iii] and Mg ii) we detect microlensing of either the blue or the red wing, or of both wings with the same amplitude. This observation implies that the broad line region is not in general spherically symmetric. In addition, the frequent detection of microlensing of the blue and red wings independently but not simultaneously with a different amplitude, does not support existing microlensing simulations of a biconical outflow. Our analysis also provides the intrinsic flux ratio between the lensed images and the magnitude of the microlensing affecting the continuum. These two quantities are particularly relevant for the determination of the fraction of matter in clumpy form in galaxies and for the detection of dark matter substructures via the identification of flux ratio anomalies.