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.1\sigma$ 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\sigma$ 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.

The local expansion rate of the Universe is parametrized by the Hubble constant, $H_0$, the ratio between recession velocity and distance. Different techniques lead to inconsistent estimates of $H_0$. Observations of Type Ia supernovae (SNe) can be used to measure $H_0$, but this requires an external calibrator to convert relative distances to absolute ones. We use the angular diameter distance to strong gravitational lenses as a suitable calibrator, which is only weakly sensitive to cosmological assumptions. We determine the angular diameter distances to two gravitational lenses, $810^{+160}_{-130}$ and $1230^{+180}_{-150}$~Mpc, at redshifts of $z=0.295$ and $0.6304$. Using these absolute distances to calibrate 740 previously-measured relative distances to SNe, we measure the Hubble constant to be $H_0=82.4^{+8.4}_{-8.3} ~{\rm km\,s^{-1}\,Mpc^{-1}}$.

Context. The precise determination of the present-day expansion rate of the Universe, expressed through the Hubble constant $H_0$, is one of the most pressing challenges in modern cosmology. Assuming flat $\Lambda$CDM, $H_0$ inference at high redshift using cosmic-microwave-background data from Planck disagrees at the 4.4$\sigma$ level with measurements based on the local distance ladder made up of parallaxes, Cepheids and Type Ia supernovae (SNe Ia), often referred to as "Hubble tension". Independent, cosmological-model-insensitive ways to infer $H_0$ are of critical importance. Aims. We apply an inverse-distance-ladder approach, combining strong-lensing time-delay-distance measurements with SN Ia data. By themselves, SNe Ia are merely good relative distance indicators, but by anchoring them to strong gravitational lenses one can obtain an $H_0$ measurement that is relatively insensitive to other cosmological parameters. Methods. A cosmological parameter estimate is performed for different cosmological background models, both for strong-lensing data alone and for the combined lensing + SNe Ia data sets. Results. The cosmological-model dependence of strong-lensing $H_0$ measurements is significantly mitigated through the inverse distance ladder. In combination with SN Ia data, the inferred $H_0$ consistently lies around 73-74 km s$^{-1}$ Mpc$^{-1}$, regardless of the assumed cosmological background model. Our results agree nicely with those from the local distance ladder, but there is a >2$\sigma$ tension with Planck results, and a ~1.5$\sigma$ discrepancy with results from an inverse distance ladder including Planck, Baryon Acoustic Oscillations and SNe Ia. Future strong-lensing distance measurements will reduce the uncertainties in $H_0$ from our inverse distance ladder.

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.

Strong lensing time-delay systems constrain cosmological parameters via the so-called time-delay distance and the angular diameter distance to the lens. In previous studies, only the former information was used. In this paper, we show that the cosmological constraints improve significantly when the latter information is also included. Specifically, the angular diameter distance plays a crucial role in breaking the degeneracy between the curvature of the Universe and the time-varying equation of state of dark energy. Using a mock sample of 55 bright quadruple lens systems based on expectations for ongoing/future imaging surveys, we find that adding the angular diameter distance information to the time-delay distance information and the cosmic microwave background data of Planck improves the constraint on the constant equation of state by 30%, on the time variation in the equation of state by a factor of two, and on the Hubble constant in the flat $\Lambda$CDM model by a factor of two. Therefore, previous forecasts for the statistical power of time-delay systems were significantly underestimated, i.e., time-delay systems are more powerful than previously appreciated.

The distance-redshift relation plays a fundamental role in constraining cosmological models. In this paper, we show that measurements of positions and time delays of strongly lensed images of a background galaxy, as well as those of the velocity dispersion and mass profile of a lens galaxy, can be combined to extract the angular diameter distance of the lens galaxy. Physically, as the velocity dispersion and the time delay give a gravitational potential ($GM/r$) and a mass ($GM$) of the lens, respectively, dividing them gives a physical size ($r$) of the lens. Comparing the physical size with the image positions of a lensed galaxy gives the angular diameter distance to the lens. A mismatch between the exact locations at which these measurements are made can be corrected by measuring a local slope of the mass profile. We expand on the original idea put forward by Paraficz and Hjorth, who analyzed singular isothermal lenses, by allowing for an arbitrary slope of a power-law spherical mass density profile, an external convergence, and an anisotropic velocity dispersion. We find that the effect of external convergence cancels out when dividing the time delays and velocity dispersion measurements. We derive a formula for the uncertainty in the angular diameter distance in terms of the uncertainties in the observables. As an application, we use two existing strong lens systems, B1608+656 ($z_{\rm L}=0.6304$) and RXJ1131$-$1231 ($z_{\rm L}=0.295$), to show that the uncertainty in the inferred angular diameter distances is dominated by that in the velocity dispersion, $\sigma^2$, and its anisotropy. We find that the current data on these systems should yield about 16% uncertainty in $D_A$ per object. This improves to 13% when we measure $\sigma^2$ at the so-called sweet-spot radius. Achieving 7% is possible if we can determine $\sigma^2$ with 5% precision.

Survey observations of the three-dimensional locations of galaxies are a powerful approach to measure the distribution of matter in the universe, which can be used to learn about the nature of dark energy, physics of inflation, neutrino masses, etc. A competitive survey, however, requires a large volume (e.g., Vsurvey is roughly 10 Gpc3) to be covered, and thus tends to be expensive. A "sparse sampling" method offers a more affordable solution to this problem: within a survey footprint covering a given survey volume, Vsurvey, we observe only a fraction of the volume. The distribution of observed regions should be chosen such that their separation is smaller than the length scale corresponding to the wavenumber of interest. Then one can recover the power spectrum of galaxies with precision expected for a survey covering a volume of Vsurvey (rather than the volume of the sum of observed regions) with the number density of galaxies given by the total number of observed galaxies divided by Vsurvey (rather than the number density of galaxies within an observed region). We find that regularly-spaced sampling yields an unbiased power spectrum with no window function effect, and deviations from regularly-spaced sampling, which are unavoidable in realistic surveys, introduce calculable window function effects and increase the uncertainties of the recovered power spectrum. While we discuss the sparse sampling method within the context of the forthcoming Hobby-Eberly Telescope Dark Energy Experiment, the method is general and can be applied to other galaxy surveys.

We present an analytic model for the local bias of dark matter halos in a LCDM universe. The model uses the halo mass density instead of the halo number density and is searched for various halo mass cuts, smoothing lengths, and redshift epoches. We find that, when the logarithmic density is used, the second-order polynomial can fit the numerical relation between the halo mass distribution and the underlying matter distribution extremely well. In this model the logarithm of the dark matter density is expanded in terms of log halo mass density to the second order. The model remains excellent for all halo mass cuts (from M_{cut}=3\times10^{11}$ to $3\times10^{12}h^{-1}M_{\odot}$), smoothing scales (from $R=5h^{-1}$Mpc to $50h^{-1}$Mpc), and redshift ranges (from z=0 to 1.0) considered in this study. The stochastic term in the relation is found not entirely random, but a part of the term can be determined by the magnitude of the shear tensor.