The H0LiCOW collaboration inferred via gravitational lensing time delays a Hubble constant $H_0=73.3^{+1.7}_{-1.8}$ km s$^{-1}{\rm Mpc}^{-1}$, describing deflector mass density profiles by either a power-law or stars 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 flexible mass models that are maximally degenerate with H0. Our calculation is based on a new hierarchical approach in which the MST is only constrained by stellar kinematics. The approach is validated on hydrodynamically simulated lenses. We apply the method to the TDCOSMO sample of 7 lenses (6 from H0LiCOW) and measure $H_0=74.5^{+5.6}_{-6.1}$ km s$^{-1}{\rm Mpc}^{-1}$. In order to further constrain the deflector mass profiles, we then add imaging and spectroscopy for 33 strong gravitational lenses from the SLACS sample. For 9 of the SLAC lenses we use resolved kinematics to constrain the stellar anisotropy. From the joint analysis of the TDCOSMO+SLACS sample, we measure $H_0=67.4^{+4.1}_{-3.2}$ km s$^{-1}{\rm Mpc}^{-1}$, assuming 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. While our new analysis does not statistically invalidate the mass profile assumptions by H0LiCOW, and thus their $H_0$ measurement relying on those, 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, chiefly spatially resolved kinematics of lens galaxies.

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 velocity 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.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

We present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the MPIA 2.2m telescope at La Silla observatory between October 2016 and February 2020. The lensed quasars HE 0047-1756, WG 0214-2105, DES 0407-5006, 2M 1134-2103, PSJ 1606-2333 and DES 2325-5229 were observed almost daily at high signal-to-noise ratio to obtain high-quality light curves where we can record fast and small-amplitude variations of the quasars. We measure time delays between all pairs of multiple images with only one or two seasons of monitoring at the exception of the time delays relative to image D of PSJ 1606-2333. The most precise estimate is obtained for the delay between image A and image B of DES 0407-5006, with $\tau_{AB} = -128.4^{+3.5}_{-3.8}$ days (2.8% precision) including systematics due to extrinsic variability in the light curves. For HE 0047-1756, we combine our high-cadence data with measurements from decade-long light-curves from previous COSMOGRAIL campaigns, and reach a precision of 0.9 day on the final measurement. The present work demonstrates the feasibility of measuring time delays in lensed quasars in only one or two seasons, provided high signal-to-noise ratio data are obtained at a cadence close to daily.

Strong gravitational lenses with measured time delay are a powerful tool to measure cosmological parameters, especially the Hubble constant ($H_0$). Recent studies show that by combining just three multiply-imaged AGN systems, one can determine $H_0$ to 2.4% precision. Furthermore, the number of time-delay lens systems is growing rapidly, enabling the determination of $H_0$ to 1% precision in the near future. However, as the precision increases it is important to ensure that systematic errors and biases remain subdominant. For this purpose, challenges with simulated datasets are a key component in this process. Following the experience of the past challenge on time delay, where it was shown that time delays can indeed be measured precisely and accurately at the sub-percent level, we now present the "Time Delay Lens Modeling Challenge" (TDLMC). The goal of this challenge is to assess the present capabilities of lens modeling codes and assumptions and test the level of accuracy of inferred cosmological parameters given realistic mock datasets. We invite scientists to model a set of simulated HST observations of 50 mock lens systems. The systems are organized in rungs, with the complexity and realism increasing going up the ladder. The goal of the challenge is to infer $H_0$ for each rung, given the HST images, the time delay, and stellar velocity dispersion of the deflector for a fixed background cosmology. The TDLMC challenge starts with the mock data release on 2018 January 8th. The deadline for blind submission is different for each rung. The deadline for Rung 0-1 is 2018 September 8; the deadline for Rung 2 is 2019 April 8 and the one for Rung 3 is 2019 September 8. This first paper gives an overview of the challenge including the data design, and a set of metrics to quantify the modeling performance and challenge details.

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 presented in paper I. 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.

One of the main challenges in using high redshift active galactic nuclei to study the correlations between the mass of the supermassive Black Hole (MBH) and the properties of their active host galaxies is instrumental resolution. Strong lensing magnification effectively increases instrumental resolution and thus helps to address this challenge. In this work, we study eight strongly lensed active galactic nuclei (AGN) with deep Hubble Space Telescope imaging, using the lens modelling code Lenstronomy to reconstruct the image of the source. Using the reconstructed brightness of the host galaxy, we infer the host galaxy stellar mass based on stellar population models. MBH are estimated from broad emission lines using standard methods. Our results are in good agreement with recent work based on non-lensed AGN, providing additional evidence that the correlation evolves over cosmic time. At the moment, the sample size of lensed AGN is small and thus they provide mostly a consistency check on systematic errors related to resolution for the non-lensed AGN. However, the number of known lensed AGN is expected to increase dramatically in the next few years, through dedicated searches in ground and space based wide field surveys, and they may become a key diagnostic of black hole and galaxy co-evolution.

We present Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimetre Array (ALMA) observations of SDSS J0924+0219, a z = 1.524 radio-quiet lensed quasar with an intrinsic radio flux density of about 3 micro-Jy. The four lensed images are clearly detected in the radio continuum and the CO(5-4) line, whose centroid is at z = 1.5254 +/- 0.0001, with a marginal detection in the submillimetre continuum. The molecular gas displays ordered motion, in a structure approximately 1--2.5 kpc in physical extent, with typical velocities of 50-100 km/s. Our results are consistent with the radio emission being emitted from the same region, but not with a point source of radio emission. SDSS J0924+0219 shows an extreme anomaly in the flux ratios of the two merging images in the optical continuum and broad emission lines, suggesting the influence of microlensing by stars in the lensing galaxy. We find the flux ratio in the radio, submillimetre continuum and CO lines to be slightly greater than 1 but much less than that in the optical, which can be reproduced with a smooth galaxy mass model and an extended source. Our results, supported by a microlensing simulation, suggest that the most likely explanation for the optical flux anomaly is indeed microlensing.

We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data is less sensitive to small-scale structures than the space-based observations. Marginalising over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, $\kappa_{\rm ext}$ at the position of B1608+656 is $\kappa_{\rm ext} = 0.11^{+0.06}_{-0.04}$, where the error bars corresponds respectively to the 16th and 84th quartiles. This result is compatible with previous estimates using the number-counts technique, suggesting that B1608+656 resides in an over-dense line of sight, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favor truncated halo models.

Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Calibrating the Baryon Acoustic Oscillations (BAO) and Type Ia supernovae with combinations of the distance indicators, we obtain a joint and self-consistent measurement of $H_0$ and $r_d$ at low redshift, independent of cosmological models and CMB inference. In an attempt to alleviate the tension between late-time and CMB-based measurements, we consider four extensions of the standard $\Lambda$CDM model. The sound horizon from our different measurements is $r_d=(137\pm3^{stat.}\pm2^{syst.})$~Mpc. Depending on the adopted distance indicators, the $combined$ tension in $H_0$ and $r_d$ ranges between 2.3 and 5.1 $\sigma$. We find that modifications of $\Lambda$CDM that change the physics after recombination fail to solve the problem, for the reason that they only resolve the tension in $H_0$, while the tension in $r_d$ remains unchanged. Pre-recombination extensions (with early dark energy or the effective number of neutrinos $\rm{N}_{\rm{eff}}=3.24 \pm 0.16$) are allowed by the data, unless the calibration from Cepheids is included. Results from time-delay lenses are consistent with those from distance-ladder calibrations and point to a discrepancy between absolute distance scales measured from the CMB (assuming the standard cosmological model) and late-time observations. New proposals to resolve this tension should be examined with respect to reconciling not only the Hubble constant but also the sound horizon derived from the CMB and other cosmological probes.

We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for this lens system using three band Hubble Space Telescope imaging. We combine the measured time delays, line-of-sight central velocity dispersion of the deflector, and statistically constrained external convergence with our lens models to estimate two cosmological distances. We measure the "effective" time-delay distance corresponding to the redshifts of the deflector and the lensed quasar $D_{\Delta t}^{\rm eff}=3382^{+146}_{-115}$ Mpc and the angular diameter distance to the deflector $D_{\rm d}=1711^{+376}_{-280}$ Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant $H_0 = 74.2^{+2.7}_{-3.0}$ km s$^{-1}$ Mpc$^{-1}$ assuming a flat $\Lambda$CDM cosmology and a uniform prior for $\Omega_{\rm m}$ as $\Omega_{\rm m} \sim \mathcal{U}(0.05, 0.5)$. This measurement gives the most precise constraint on $H_0$ to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analyzed by the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of $H_0$ based on the local distance ladder, reinforcing the tension with the inference from early Universe probes, for example, with 2.2$\sigma$ discrepancy from the cosmic microwave background measurement.

We present the HOLISMOKES project on strong gravitational lensing of supernovae (SNe) as a probe of SN physics and cosmology. We investigate the effects of microlensing on early-phase SN Ia spectra, and find that within 10 rest-frame days after SN explosion, distortions of SN Ia spectra due to microlensing are typically negligible (<1% distortion within the 68% credible region, and ~10% distortion within the 95% credible region). This shows great prospects of using lensed SNe Ia to obtain intrinsic early-phase SN spectra for deciphering SN Ia progenitors. As a demonstration of the usefulness of lensed SNe Ia for cosmology, we simulate a sample of mock lensed SN Ia systems that are expected to have accurate and precise time-delay measurements in the era of the Legacy Survey of Space and Time (LSST). Adopting realistic yet conservative uncertainties on their time-delay distances and lens angular diameter distances (of 6.6% and 5%, respectively), we find that a sample of 20 lensed SNe Ia would allow a constraint on the Hubble constant ($H_0$) with 1.3% uncertainty in the flat $\Lambda$CDM cosmology. We find a similar constraint on $H_0$ in an open $\Lambda$CDM cosmology, while the constraint degrades to 3% in a flat wCDM cosmology. We anticipate lensed SNe to be an independent and powerful probe of SN physics and cosmology in the upcoming LSST era.

We present the results of 15 years of monitoring of lensed quasars conducted by the COSMOGRAIL program at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed QSOs are presented for the first time. We complement our data set with other monitoring data available in the literature to measure the time delays in 18 systems, among which 9 reach a relative precision better than 15% for at least one time delay. To achieve this, we develop an automated version of the curve-shifting toolbox PyCS to ensure robust estimation of the time delay in the presence of microlensing while accounting for the errors due to imperfect representation of microlensing. We also reanalyze the previously published time delays of RX J1131$-$1231 and HE 0435$-$1223, adding respectively 6 and 2 new seasons of monitoring and confirming the previous time-delay measurements. When the time delay measurement is possible, we correct the light curves of the lensed images from their time delay and present the difference curves to highlight the microlensing signal contained in the data. This is to date the largest sample of decade-long lens monitoring data, useful to measure $H_0$, to measure the size of quasar accretion disks with microlensing, and to study quasar variability.

The discovery of multiply-imaged gravitationally lensed QSOs is fundamental to many astronomical and cosmological studies. However, these objects are rare and challenging to discover due to requirements of high-angular resolution astrometric, multiwavelength photometric and spectroscopic data. This has limited the number of known systems to a few hundred objects. We aim to reduce the constraints on angular resolution and discover multiply-imaged QSO candidates by using new candidate selection principles based on unresolved photometric time-series and ground-based images from public surveys. We selected candidates for multiply-imaged QSOs based on low levels of entropy computed from Catalina unresolved photometric time-series or Euclidean similarity to known lenses in a space defined by the wavelet power spectra of Pan-STARSS DR2 or DECaLS DR7 images, combined with multiple {\it Gaia} DR2 sources or large astrometric errors and supervised and unsupervised learning methods. We then confirmed spectroscopically some candidates with the Palomar Hale, Keck-I, and ESO/NTT telescopes. Here we report the discovery and confirmation of seven doubly-imaged QSOs and one likely double quasar. This demonstrates the potential of combining space-astrometry, even if unresolved, with low spatial-resolution photometric time-series and/or low-spatial resolution multi-band imaging to discover multiply-imaged lensed QSOs.

The magnifications of compact-source lenses are extremely sensitive to the presence of low mass dark matter halos along the entire sight line from the source to the observer. Traditionally, the study of dark matter structure in compact-source strong gravitational lenses has been limited to radio-loud systems, as the radio emission is extended and thus unaffected by microlensing which can mimic the signal of dark matter structure. An alternate approach is to measure quasar nuclear-narrow line emission, which is free from microlensing and present in virtually all quasar lenses. In this paper, we double the number of systems which can be used for gravitational lensing analyses by presenting measurements of narrow-line emission from a sample of 8 quadruply imaged quasar lens systems, WGD J0405-3308, HS 0810+2554, RX J0911+0551, SDSS J1330+1810, PS J1606-2333, WFI 2026-4536, WFI 2033-4723 and WGD J2038-4008. We describe our updated grism spectral modelling pipeline, which we use to measure narrow-line fluxes with uncertainties of 2-10\%, presented here. We fit the lensed image positions with smooth mass models and demonstrate that these models fail to produce the observed distribution of image fluxes over the entire sample of lenses. Furthermore, typical deviations are larger than those expected from macromodel uncertainties. This discrepancy indicates the presence of perturbations caused by small-scale dark matter structure. The interpretation of this result in terms of dark matter models is presented in a companion paper.

Gravitational microlensing is a powerful tool for probing the inner structure of distant quasars. In this context, we have obtained spectropolarimetric observations of the two images of the broad absorption line (BAL) quasar SDSSJ081830.46+060138.0 (J0818+0601) at redshift $z \simeq$ 2.35. We first show that J0818+0601 is actually gravitationally lensed, and not a binary quasar. A strong absorption system detected at $z$ = 1.0065$\pm$0.0002 is possibly due to the lensing galaxy. Microlensing is observed in one image and it magnifies the emission lines, the continuum, and the BALs differently. By disentangling the part of the spectrum that is microlensed from the part that is not microlensed, we unveil two sources of continuum that must be spatially separated: a compact one, which is microlensed, and an extended one, which is not microlensed and contributes to two thirds of the total continuum emission. J0818+0601 is the second BAL quasar in which an extended source of rest-frame ultraviolet continuum is found. We also find that the images are differently polarized, suggesting that the two continua might be differently polarized. Our analysis provides constraints on the BAL flow. In particular, we find that the outflow is seen with a nonzero onset velocity, and stratified according to ionization.

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 $\textit{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 of 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\%$. This translates to a Hubble constant $H_{0} = 71.6_{-4.9}^{+3.8}~\mathrm{km~s^{-1}~Mpc^{-1}}$, assuming a flat $\Lambda$CDM cosmology with a uniform prior on $\Omega_\mathrm{m}$ in the range [0.05, 0.5]. This work is part of the $H_0$ 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).

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.

Galaxies and galaxy groups located along the line of sight towards gravitationally lensed quasars produce high-order perturbations of the gravitational potential at the lens position. When these perturbation are too large, they can induce a systematic error on $H_0$ of a few-percent if the lens system is used for cosmological inference and the perturbers are not explicitly accounted for in the lens model. In this work, we present a detailed characterization of the environment of the lens system WFI2033-4723 ($z_{\rm src} = 1.662$, $z_{\rm lens}$ = 0.6575), one of the core targets of the H0LICOW project for which we present cosmological inferences in a companion paper (Rusu et al. 2019). We use the Gemini and ESO-Very Large telescopes to measure the spectroscopic redshifts of the brightest galaxies towards the lens, and use the ESO-MUSE integral field spectrograph to measure the velocity-dispersion of the lens ($\sigma_{\rm {los}}= 250^{+15}_{-21}$ km/s) and of several nearby galaxies. In addition, we measure photometric redshifts and stellar masses of all galaxies down to $i < 23$ mag, mainly based on Dark Energy Survey imaging (DR1). Our new catalog, complemented with literature data, more than doubles the number of known galaxy spectroscopic redshifts in the direct vicinity of the lens, expanding to 116 (64) the number of spectroscopic redshifts for galaxies separated by less than 3 arcmin (2 arcmin) from the lens. Using the flexion-shift as a measure of the amplitude of the gravitational perturbation, we identify 2 galaxy groups and 3 galaxies that require specific attention in the lens models. The ESO MUSE data enable us to measure the velocity-dispersions of three of these galaxies. These results are essential for the cosmological inference analysis presented in Rusu et al. (2019).

The quadruply lensed quasar HE0435-1223 shows a clear microlensing effect that affects differently the blue and red wings of the H$\alpha$ line profile in its image D. To interpret these observations, and constrain the broad emission line region (BLR) properties, the effect of gravitational microlensing on quasar broad emission line profiles and their underlying continuum has been simulated considering representative BLR models and microlensing magnification maps. The amplification and distortion of the H$\alpha$ line profile, characterized by a set of four indices, can be reproduced by the simulations. Although the constraints on the BLR models set by the observed single-epoch microlensing signal are not very robust, we found that flattened geometries (Keplerian disk and equatorial wind) can more easily reproduce the observed line profile deformations than a biconical polar wind. With an additional independent constraint on the size of the continuum source, the Keplerian disk model of the H$\alpha$ BLR is slightly favored.

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.

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}$ kms$^{-1}$Mpc$^{-1}$, assuming a flat Lambda cold dark matter cosmology with uniform prior on $\Omega_{\rm 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}$kms$^{-1}$Mpc$^{-1}$. This measurement is independent of the distance ladder and other cosmological probes.

We present new measurements of the time delays of WFI2033-4723. The data sets used in this work include 14 years of data taken at the 1.2m Leonhard Euler Swiss telescope, 13 years of data from the SMARTS 1.3m telescope at Las Campanas Observatory and a single year of high-cadence and high-precision monitoring at the MPIA 2.2m telescope. The time delays measured from these different data sets, all taken in the R-band, are in good agreement with each other and with previous measurements from the literature. Combining all the time-delay estimates from our data sets results in Dt_AB = 36.2-0.8+0.7 days (2.1% precision), Dt_AC = -23.3-1.4+1.2 days (5.6%) and Dt_BC = -59.4-1.3+1.3 days (2.2%). In addition, the close image pair A1-A2 of the lensed quasars can be resolved in the MPIA 2.2m data. We measure a time delay consistent with zero in this pair of images. We also explore the prior distributions of microlensing time-delay potentially affecting the cosmological time-delay measurements of WFI2033-4723. There is however no strong indication in our measurements that microlensing time delay is neither present nor absent. This work is part of a H0LiCOW series focusing on measuring the Hubble constant from WFI2033-4723.

Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide information from distinct emission regions in the quasar. Through monitoring of these strong gravitational lenses, a full tomographic view can emerge with Astronomical-Unit scale resolution. Work to date has demonstrated the potential of this technique in about a dozen systems. In the 2020s there will be orders of magnitude more systems to work with. Monitoring of lens systems for caustic-crossing events to enable triggering of multi-platform, multi-wavelength observations in the 2020s will fulfill the potential of quasar microlensing as a unique and comprehensive probe of active black hole structure and dynamics.

Optical and X-ray observations of strongly gravitationally lensed quasars (especially when four separate images of the quasar are produced) determine not only the amount of matter in the lensing galaxy but also how much is in a smooth component and how much is composed of compact masses (e.g., stars, stellar remnants, primordial black holes, CDM sub-halos, and planets). Future optical surveys will discover hundreds to thousands of quadruply lensed quasars, and sensitive X-ray observations will unambiguously determine the ratio of smooth to clumpy matter at specific locations in the lensing galaxies and calibrate the stellar mass fundamental plane, providing a determination of the stellar $M/L$. A modest observing program with a sensitive, sub-arcsecond X-ray imager, combined with the planned optical observations, can make those determinations for a large number (hundreds) of the lensing galaxies, which will span a redshift range of $\sim$$0.25<z<1.5$

If the disappearance of the broad emission lines observed in changing-look quasars originates from the obscuration of the quasar core by dusty clouds moving in the torus, high linear optical polarization would be expected in those objects. We then measured the rest-frame UV-blue linear polarization of a sample of 13 changing-look quasars, 7 of them being in a type 1.9-2 state. For all quasars but one the polarization degree is lower than 1%. This suggests that the disappearance of the broad emission lines cannot be attributed to dust obscuration, and supports the scenario in which changes of look are caused by a change in the rate of accretion onto the supermassive black hole. Such low polarization degrees also indicate that these quasars are seen under inclinations close to the system axis. One type 1.9-2 quasar in our sample shows a high polarization degree of 6.8%. While this polarization could be ascribed to obscuration by a moving dusty cloud, we argue that this is unlikely given the very long time needed for a cloud from the torus to eclipse the broad emission line region of that object. We propose that the high polarization is due to the echo of a past bright phase seen in polar-scattered light. This interpretation raises the possibility that broad emission lines observed in the polarized light of some type 2 active galactic nuclei can be echoes of past type 1 phases and not evidence of hidden broad emission line regions.

This Astro2020 white paper summarizes the unknowns of active galactic nuclei (AGN) physics that could be unveiled thanks to a new, space-born, ultraviolet spectropolarimeter. The unique capabilities of high energy polarimetry would help us to understand the precise mechanisms of matter and energy transfer and supermassive black holes growth, together with the impact of AGN feedback on galaxy evolution.

We present e-MERLIN and EVN observations which reveal unambiguous jet activity within radio quiet quasar HS~0810+2554. With an intrinsic flux density of 880~nJy, this is the faintest radio source ever imaged. The findings present new evidence against the idea that radio loud and radio quiet quasars are powered by different underlying radio emission mechanisms, showing instead that the same AGN mechanism can operate as the dominant source of radio emission even in the very lowest radio luminosity quasars. Thanks to strong gravitational lensing, our source is not only visible, but with VLBI is imaged to a scale of just 0.27~pc: the highest ever resolution image of a radio quiet quasar. Brightness temperatures of at least $8.4\times 10^6$~K are associated with two highly compact components, and subsequent modelling of the lensed system has revealed that the components are linearly aligned on opposing sides of the optical quasar core, with the typical morphology of a compact symmetric object (CSO). Given that this source has been found to fall on the radio--FIR correlation, we suggest that the radio--FIR correlation cannot always be used to rule out AGN activity in favour of star-formation activity. The correlation -- or at least its scatter -- may conceal the coexistence of kinetic and radiative feedback modes in AGN. Modelling of the lensing mass itself points to a non-smooth mass distribution, hinting at the presence of dark matter substructure which has manifested as astrometric perturbations of the VLBI lensed images, posing no threat to the CDM paradigm.

The ultraviolet (UV) polarization spectrum of nearby active galactic nuclei (AGN) is poorly known. The Wisconsin Ultraviolet Photo-Polarimeter Experiment and a handful of instruments on board the Hubble Space Telescope were able to probe the near- and mid-UV polarization of nearby AGN, but the far-UV band (from 1200 angs down to the Lyman limit at 912 angs) remains completely uncharted. In addition, the linewidth resolution of previous observations was at best 1.89 angs. Such a resolution is not sufficient to probe in detail quantum mechanical effects, synchrotron and cyclotron processes, scattering by electrons and dust grains, and dichroic extinction by asymmetric dust grains. Exploring those physical processes would require a new, high-resolution, broadband polarimeter with full ultraviolet-band coverage. In this context, we discuss the AGN science case for POLLUX, a high-resolution UV spectropolarimeter, proposed for the 15-meter primary mirror option of LUVOIR (a multi-wavelength space observatory concept being developed by the Goddard Space Flight Center and proposed for the 2020 Decadal Survey Concept Study).

We report 87 previously unpublished optical linear polarization measurements of 86 quasars obtained in May and October 2008, and from April to July 2015 with the Focal Reducer and low-dispersion Spectrographs FORS1 and FORS2 attached to the Very Large Telescope at the Paranal Observatory. Of the 86 quasars, 37 have $p \geq 0.6\%$, 9 have $p \geq 2\%$, and 1 has $p \geq 10\%$.

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 ALMA 2-mm continuum and CO (2-1) spectral line imaging of the gravitationally lensed z=0.654 star-forming/quasar composite RX J1131-1231 at 240-400 mas angular resolution. The continuum emission is found to be compact and coincident with the optical emission, whereas the molecular gas forms a complete Einstein ring, which shows strong differential magnification. The de-lensed source structure is determined on 400-pc resolution using a visibility-fitting lens modelling technique. The reconstructed molecular gas velocity-field is consistent with a rotating disk with a maximum rotational velocity of 280 km/s. From dynamical model fitting we find an enclosed mass M(r<5 kpc)=(1.46+/-0.31)*10^11 M_sol. The molecular gas distribution is highly structured, with clumps that are co-incident with higher gas velocity dispersion regions 40-50 km/s and with the intensity peaks in the optical emission, which are associated with sites of on-going turbulent star-formation. The peak in the CO (2-1) distribution is not co-incident with the AGN, where there is a paucity of molecular gas emission, possibly due to radiative feedback from the central engine. The intrinsic molecular gas luminosity is L'_CO=(1.2+/-0.3)*10^10 K km/s pc^2 and the inferred gas mass is M(H2)=(8.3+/-3.0)*10^10 M_sol, which given its dynamical mass is consistent with a CO-H2 conversion factor of alpha = 5.5+/-2.0 M_solar(K km/s pc^2)^-1. This suggests that the star-formation efficiency is dependent on the host galaxy morphology as opposed to the nature of the AGN. The far-infrared continuum spectral energy distribution shows evidence for heated dust, equivalent to an obscured star-formation rate of SFR=69^(+41)_(-25)*(7.3/u_IR)M_sol/yr, which demonstrates the composite star-forming/AGN nature of this system. RX J1131-1231

We present a weak gravitational lensing measurement of the external convergence along the line of sight to the quadruply lensed quasar HE$\,$0435$-$1223. Using deep r-band images from Subaru-Suprime-Cam we observe galaxies down to a 3$\sigma$ limiting magnitude of $\sim 26$ mags resulting in a source galaxy density of 14 galaxies / arcmin$^2$ after redshift-based cuts. Using an inpainting technique and Multi-Scale Entropy filtering algorithm, we find that the region in close proximity to the lens has an estimated external convergence of $\kappa=-0.012^{+0.020}_{-0.013}$ and is hence marginally under-dense. We also rule out the presence of any halo with a mass greater than $M_{\rm vir}=1.6\times10^{14}h^{-1}M_\odot$ (68$\%$ confidence limit). Our results, consistent with previous studies of this lens, confirm that the intervening mass along the line of sight to HE$\,$0435$-$1223 does not affect significantly the cosmological results inferred from the time delay measurements of that specific object.

Gravitational microlensing is a powerful tool allowing one to probe the structure of quasars on sub-parsec scale. We report recent results, focusing on the broad absorption and emission line regions. In particular microlensing reveals the intrinsic absorption hidden in the P Cygni-type line profiles observed in the broad absorption line quasar H1413+117, as well as the existence of an extended continuum source. In addition, polarization microlensing provides constraints on the scattering region. In the quasar Q2237+030, microlensing differently distorts the H$\alpha$ and CIV broad emission line profiles, indicating that the low- and high-ionization broad emission lines must originate from regions with distinct kinematical properties. We also present simulations of the effect of microlensing on line profiles considering simple but representative models of the broad emission line region. Comparison of observations to simulations allows us to conclude that the H$\alpha$ emitting region in Q2237+030 is best represented by a Keplerian disk.

We report 192 previously unpublished optical linear polarization measurements of quasars obtained in April 2003, April 2007, and October 2007 with the European Southern Observatory Faint Object Spectrograph and Camera (EFOSC2) instrument attached to the 3.6m telescope at the La Silla Observatory. Each quasar was observed once. Among the 192 quasars, 89 have a polarization degree $p \geq 0.6\%$, 18 have $p \geq 2\%$, and two have $p \geq 10\%$.

Strongly lensed active galactic nuclei (AGN) provide a unique opportunity to make progress in the study of the evolution of the correlation between the mass of supermassive black holes ($\mathcal M_{BH}$) and their host galaxy luminosity ($L_{host}$). We demonstrate the power of lensing by analyzing two systems for which state-of-the-art lens modelling techniques have been applied to Hubble Space Telescope imaging data. We use i) the reconstructed images to infer the total and bulge luminosity of the host and ii) published broad-line spectroscopy to estimate $\mathcal M_{BH}$ using the so-called virial method. We then enlarge our sample with new calibration of previously published measurements to study the evolution of the correlation out to z~4.5. Consistent with previous work, we find that without taking into account passive luminosity evolution, the data points lie on the local relation. Once passive luminosity evolution is taken into account, we find that BHs in the more distant Universe reside in less luminous galaxies than today. Fitting this offset as $\mathcal M_{BH}$/$L_{host}$ $\propto$ (1+z)$^{\gamma}$, and taking into account selection effects, we obtain $\gamma$ = 0.6 $\pm$ 0.1 and 0.8$\pm$ 0.1 for the case of $\mathcal M_{BH}$-$L_{bulge}$ and $\mathcal M_{BH}$-$L_{total}$, respectively. To test for systematic uncertainties and selection effects we also consider a reduced sample that is homogeneous in data quality. We find consistent results but with considerably larger uncertainty due to the more limited sample size and redshift coverage ($\gamma$ = 0.7 $\pm$ 0.4 and 0.2$\pm$ 0.5 for $\mathcal M_{BH}$-$L_{bulge}$ and $\mathcal M_{BH}$-$L_{total}$, respectively), highlighting the need to gather more high-quality data for high-redshift lensed quasar hosts. Our result is consistent with a scenario where the growth of the black hole predates that of the host galaxy.

Recent studies have shown that line profile distortions are commonly observed in gravitationally lensed quasar spectra. We investigate the effect of gravitational microlensing on quasar broad emission line profiles and their underlying continuum, combining the emission from simple representative BLR models with generic microlensing magnification maps. Specifically, we considered Keplerian disk, polar, and equatorial wind BLR models of various sizes. The effect of microlensing has been quantified with four observables: $\mu^{BLR}$, the total magnification of the broad emission line; $\mu^{cont}$, the magnification of the underlying continuum; as well as red/blue, RBI and wings/core, WCI, indices that characterize the line profile distortions. The simulations showed that distortions of line profiles, such as those recently observed in lensed quasars, can indeed be reproduced and attributed to the differential effect of microlensing on spatially separated regions of the BLR. While the magnification of the emission line $\mu^{BLR}$ sets an upper limit on the BLR size and, similarly, the magnification of the continuum $\mu^{cont}$ sets an upper limit on the size of the continuum source, the line profile distortions mainly depend on the BLR geometry and kinematics. We thus built (WCI,RBI) diagrams that can serve as diagnostic diagrams to discriminate between the various BLR models on the basis of quantitative measurements. It appears that a strong microlensing effect puts important constraints on the size of the BLR and on its distance to the high-magnification caustic. In that case, BLR models with different geometries and kinematics are more prone to produce distinctive line profile distortions for a limited number of caustic configurations, which facilitates their discrimination.

If the disappearance of the broad emission lines observed in changing-look quasars were caused by the obscuration of the quasar core through moving dust clouds in the torus, high linear polarization typical of type 2 quasars would be expected. We measured the polarization of the changing-look quasar J1011+5442 in which the broad emission lines have disappeared between 2003 and 2015. We found a polarization degree compatible with null polarization. This measurement suggests that the observed change of look is not due to a change of obscuration hiding the continuum source and the broad line region, and that the quasar is seen close to the system axis. Our results thus support the idea that the vanishing of the broad emission lines in J1011+5442 is due to an intrinsic dimming of the ionizing continuum source that is most likely caused by a rapid decrease in the rate of accretion onto the supermassive black hole.

Strong gravitational lensing provides a powerful test of Cold Dark Matter (CDM) as it enables the detection and mass measurement of low mass haloes even if they do not contain baryons. Compact lensed sources such as Active Galactic Nuclei (AGN) are particularly sensitive to perturbing subhalos, but their use as a test of CDM has been limited by the small number of systems which have significant radio emission which is extended enough avoid significant lensing by stars in the plane of the lens galaxy, and red enough to be minimally affected by differential dust extinction. Narrow-line emission is a promising alternative as it is also extended and, unlike radio, detectable in virtually all optically selected AGN lenses. We present first results from a WFC3 grism narrow-line survey of lensed quasars, for the quadruply lensed AGN HE0435-1223. Using a forward modelling pipeline which enables us to robustly account for spatial blending, we measure the [OIII] 5007 \AA~ flux ratios of the four images. We find that the [OIII] fluxes and positions are well fit by a simple smooth mass model for the main lens. Our data rule out a $M_{600}>10^{8} (10^{7.2}) M_\odot$ NFW perturber projected within $\sim$1\farcs0 (0\farcs1) arcseconds of each of the lensed images, where $M_{600}$ is the perturber mass within its central 600 pc. The non-detection is broadly consistent with the expectations of $\Lambda$CDM for a single system. The sensitivity achieved demonstrates that powerful limits on the nature of dark matter can be obtained with the analysis of $\sim20$ narrow-line lenses.

The radio-quiet quasar Q2059-360 at redshift $z=3.08$ is known to be close to a small Lyman $\alpha$ blob (LAB) and to be absorbed by a proximate damped Ly$\alpha$ (PDLA) system. Here, we present the Multi Unit Spectroscopic Explorer (MUSE) integral field spectroscopy follow-up of this quasi-stellar object (QSO). Our primary goal is to characterize this LAB in detail by mapping it both spatially and spectrally using the Ly$\alpha$ line, and by looking for high-ionization lines to constrain the emission mechanism. Combining the high sensitivity of the MUSE integral field spectrograph mounted on the Yepun telescope at ESO-VLT with the natural coronagraph provided by the PDLA, we map the LAB down to the QSO position, after robust subtraction of QSO light in the spectral domain. In addition to confirming earlier results for the small bright component of the LAB, we unveil a faint filamentary emission protruding to the south over about 80 pkpc (physical kpc); this results in a total size of about 120 pkpc. We derive the velocity field of the LAB (assuming no transfer effects) and map the Ly$\alpha$ line width. Upper limits are set to the flux of the N V $\lambda 1238-1242$, C IV $\lambda 1548-1551$, He II $\lambda 1640$, and C III] $\lambda 1548-1551$ lines. We have discovered two probable Ly$\alpha$ emitters at the same redshift as the LAB and at projected distances of 265 kpc and 207 kpc from the QSO; their Ly$\alpha$ luminosities might well be enhanced by the QSO radiation. We also find an emission line galaxy at $z=0.33$ near the line of sight to the QSO. This LAB shares the same general characteristics as the 17 others surrounding radio-quiet QSOs presented previously. However, there are indications that it may be centered on the PDLA galaxy rather than on the QSO.

Early-type galaxies provide unique tests for the predictions of the cold dark matter cosmology and the baryonic physics assumptions entering models for galaxy formation. In this work, we use the Illustris simulation to study correlations of three main properties of early-type galaxies, namely, the stellar orbital anisotropies, the central dark matter fractions and the central radial density slopes, as well as their redshift evolution since $z=1.0$. We find that lower-mass galaxies or galaxies at higher redshift tend to be bluer in rest-frame colour, have higher central gas fractions, and feature more tangentially anisotropic orbits and steeper central density slopes than their higher-mass or lower-redshift counterparts, respectively. The projected central dark matter fraction within the effective radius shows a very mild mass dependence but positively correlates with galaxy effective radii due to the aperture effect. The central density slopes obtained by combining strong lensing measurements with single aperture kinematics are found to differ from the true density slopes. We identify systematic biases in this measurement to be due to two common modelling assumptions, isotropic stellar orbital distributions and power-law density profiles. We also compare the properties of early-type galaxies in Illustris to those from existing galaxy and strong lensing surveys, we find in general broad agreement but also some tension, which poses a potential challenge to the stellar formation and feedback models adopted by the simulation.

A power-law density model, i.e., $\rho(r) \propto r^{-\gamma'}$ has been commonly employed in strong gravitational lensing studies, including the so-called time-delay technique used to infer the Hubble constant $H_0$. However, since the radial scale at which strong lensing features are formed corresponds to the transition from the dominance of baryonic matter to dark matter, there is no known reason why galaxies should follow a power law in density. The assumption of a power law artificially breaks the mass-sheet degeneracy, a well-known invariance transformation in gravitational lensing which affects the product of Hubble constant and time delay and can therefore cause a bias in the determination of $H_0$ from the time-delay technique. In this paper, we use the Illustris hydrodynamical simulations to estimate the amplitude of this bias, and to understand how it is related to observational properties of galaxies. Investigating a large sample of Illustris galaxies that have velocity dispersion $\sigma_{SIE}$>160 km/s at redshifts below $z=1$, we find that the bias on $H_0$ introduced by the power-law assumption can reach 20%-50%, with a scatter of $10\%-30\%$ (rms). However, we find that by selecting galaxies with an inferred power-law model slope close to isothermal, it is possible to reduce the bias on $H_0$ to <5%, and the scatter to <10%. This could potentially be used to form less biased statistical samples for $H_0$ measurements in the upcoming large survey era.

We present a new measurement of the Hubble Constant H0 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 HE0435-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 time-delay distance of HE0435-1223 with previous measurements from systems B1608+656 and RXJ1131-1231 to create a Time Delay Strong Lensing probe (TDSL). In flat $\Lambda$CDM with free matter and energy density, we find $H_0$ = 71.9 +2.4 -3.0 km/s/Mpc and $\Omega_{\Lambda}$ = 0.62 +0.24 -0.35 . This measurement is completely independent of, and in agreement with, the local distance ladder measurements of H0. We explore more general cosmological models combining TDSL with other probes, illustrating its power to break degeneracies inherent to other methods. The TDSL and Planck joint constraints are $H_0$ = 69.2 +1.4 -2.2 km/s/Mpc, $\Omega_{\Lambda}$ = 0.70 +0.01 -0.01 and $\Omega_k$ = 0.003 +0.004 -0.006 in open $\Lambda$CDM and $H_0$ = 79.0 +4.4 -4.2 km/s/Mpc, $\Omega_{de}$ = 0.77 +0.02 -0.03 and $w$ = -1.38 +0.14 -0.16 in flat $w$CDM. Combined 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 and when relaxing the total mass of neutrinos we find 0.182 eV. In an open $w$CDM in combination with Planck and CMB lensing we find $H_0$ = 77.9 +5.0 -4.2 km/s/Mpc, $\Omega_{de}$ = 0.77 +0.03 -0.03, $\Omega_k$ = -0.003 +0.004 -0.004 and $w$ = -1.37 +0.18 -0.23.

Galaxies located in the environment or on the line of sight towards gravitational lenses can significantly affect lensing observables, and can lead to systematic errors on the measurement of $H_0$ from the time-delay technique. We present the results of a systematic spectroscopic identification of the galaxies in the field of view of the lensed quasar HE0435-1223, using the W. M. Keck, Gemini and ESO-Very Large telescopes. Our new catalog triples the number of known galaxy redshifts in the vicinity of the lens, expanding to 102 the number of measured redshifts for galaxies separated by less than 3 arcmin from the lens. We complement our catalog with literature data to gather redshifts up to 15 arcmin from the lens, and search for galaxy groups or clusters projected towards HE0435-1223. We confirm that the lens is a member of a small group that includes at least 12 galaxies, and find 8 other group candidates near the line of sight of the lens. The flexion shift, namely the shift of lensed images produced by high order perturbation of the lens potential, is calculated for each galaxy/group and used to identify which objects produce the largest perturbation of the lens potential. This analysis demonstrates that i) at most three of the five brightest galaxies projected within 12 arcsec of the lens need to be explicitly used in the lens models, and ii) the groups can be treated in the lens model as an external tidal field (shear) contribution.

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, HE0435-1223, WFI2033-4723 and HE1104-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 $\Lambda$CDM, we expect to measure $H_0$ to $<3.5\%$ in most models, spatial curvature $\Omega_{\rm k}$ to 0.004, $w$ to 0.14, and the effective number of neutrino species to 0.2 (1$\sigma$ uncertainties) when combined with current CMB experiments. These are, respectively, a factor of $\sim15$, $\sim2$, and $\sim1.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.

Strong gravitational lenses with measured time delays between the multiple images allow a direct measurement of the time-delay distance to the lens, and thus a measure of cosmological parameters, particularly the Hubble constant, $H_{0}$. We present a blind lens model analysis of the quadruply-imaged quasar lens HE 0435-1223 using deep Hubble Space Telescope imaging, updated time-delay measurements from the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL), a measurement of the velocity dispersion of the lens galaxy based on Keck data, and a characterization of the mass distribution along the line of sight. HE 0435-1223 is the third lens analyzed as a part of the $H_{0}$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) project. We account for various sources of systematic uncertainty, including the detailed treatment of nearby perturbers, the parameterization of the galaxy light and mass profile, and the regions used for lens modeling. We constrain the effective time-delay distance to be $D_{\Delta t} = 2612_{-191}^{+208}~\mathrm{Mpc}$, a precision of 7.6%. From HE 0435-1223 alone, we infer a Hubble constant of $H_{0} = 73.1_{-6.0}^{+5.7}~\mathrm{km~s^{-1}~Mpc^{-1}}$ assuming a flat $\Lambda$CDM cosmology. The cosmographic inference based on the three lenses analyzed by H0LiCOW to date is presented in a companion paper (H0LiCOW Paper V).

The empirical correlation between the mass of a super-massive black hole (MBH) and its host galaxy properties is widely considered to be evidence of their co-evolution. A powerful way to test the co-evolution scenario and learn about the feedback processes linking galaxies and nuclear activity is to measure these correlations as a function of redshift. Unfortunately, currently MBH can only be estimated in active galaxies at cosmological distances. At these distances, bright active galactic nuclei (AGN) can outshine the host galaxy, making it extremely difficult to measure the host's luminosity. Strongly lensed AGNs provide in principle a great opportunity to improve the sensitivity and accuracy of the host galaxy luminosity measurements as the host galaxy is magnified and more easily separated from the point source, provided the lens model is sufficiently accurate. In order to measure the MBH-L correlation with strong lensing, it is necessary to ensure that the lens modelling is accurate, and that the host galaxy luminosity can be recovered to at least a precision and accuracy better than that of the typical MBH measurement. We carry out extensive and realistic simulations of deep Hubble Space Telescope observations of lensed AGNs obtained by our collaboration. We show that the host galaxy luminosity can be recovered with better accuracy and precision than the typical uncertainty on MBH(~ 0.5 dex) for hosts as faint as 2-4 magnitudes dimmer than the AGN itself. Our simulations will be used to estimate bias and uncertainties on the actual measurements to be presented in a future paper.

Based on spectroscopy and multiband wide-field observations of the gravitationally lensed quasar HE 0435-1223, we determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for this system. We measure the under/overdensity of the line of sight towards the lens system and compare it to the average line of sight throughout the universe, determined by using the CFHTLenS as a control field. Aiming to constrain $\kappa_\mathrm{ext}$ as tightly as possible, we determine under/overdensities using various combinations of relevant informative weighing schemes for the galaxy counts, such as projected distance to the lens, redshift, and stellar mass. We then convert the measured under/overdensities into a $\kappa_\mathrm{ext}$ distribution, using ray-tracing through the Millennium Simulation. We explore several limiting magnitudes and apertures, and account for systematic and statistical uncertainties relevant to the quality of the observational data, which we further test through simulations. Our most robust estimate of $\kappa_\mathrm{ext}$ has a median value $\kappa^\mathrm{med}_\mathrm{ext} = 0.004$ and a standard deviation of $\sigma_\kappa = 0.025$. The measured $\sigma_\kappa$ corresponds to $2.5\%$ uncertainty on the time delay distance, and hence the Hubble constant $H_0$ inference from this system. The median $\kappa^\mathrm{med}_\mathrm{ext}$ value is robust to $\sim0.005$ (i.e. $\sim0.5\%$ on $H_0$) regardless of the adopted aperture radius, limiting magnitude and weighting scheme, as long as the latter incorporates galaxy number counts, the projected distance to the main lens, and a prior on the external shear obtained from mass modeling. The availability of a well-constrained $\kappa_\mathrm{ext}$ makes \hequad\ a valuable system for measuring cosmological parameters using strong gravitational lens time delays.

Strong gravitational lensing is regarded as the most precise technique to measure the mass in the inner region of galaxies or galaxy clusters. In particular, the mass within one Einstein radius can be determined with an accuracy of order of a few percent or better, depending on the image configuration. For other radii, however, degeneracies exist between galaxy density profiles, precluding an accurate determination of the enclosed mass. The source position transformation (SPT), which includes the well-known mass-sheet transformation (MST) as a special case, describes this degeneracy of the lensing observables in a more general way. In this paper we explore properties of an SPT, removing the MST to leading order, i.e., we consider degeneracies which have not been described before. The deflection field $\boldsymbol{\hat{\alpha}}(\boldsymbol{\theta})$ resulting from an SPT is not curl-free in general, and thus not a deflection that can be obtained from a lensing mass distribution. Starting from a variational principle, we construct lensing potentials that give rise to a deflection field $\boldsymbol{\tilde{\alpha}}$, which differs from $\boldsymbol{\hat{\alpha}}$ by less than an observationally motivated upper limit. The corresponding mass distributions from these 'valid' SPTs are studied: their radial profiles are modified relative to the original mass distribution in a significant and non-trivial way, and originally axi-symmetric mass distributions can obtain a finite ellipticity. These results indicate a significant effect of the SPT on quantitative analyses of lens systems. We show that the mass inside the Einstein radius of the original mass distribution is conserved by the SPT; hence, as is the case for the MST, the SPT does not affect the mass determination at the Einstein radius. [...]

We investigate the kinematics and ionization structure of the broad emission line region of the gravitationally lensed quasar QSO2237+0305 (the Einstein cross) using differential microlensing in the high- and low-ionization broad emission lines. We combine visible and near-infrared spectra of the four images of the lensed quasar and detect a large-amplitude microlensing effect distorting the high-ionization CIV and low-ionization H$\alpha$ line profiles in image A. While microlensing only magnifies the red wing of the Balmer line, it symmetrically magnifies the wings of the CIV emission line. Given that the same microlensing pattern magnifies both the high- and low-ionization broad emission line regions, these dissimilar distortions of the line profiles suggest that the high- and low-ionization regions are governed by different kinematics. Since this quasar is likely viewed at intermediate inclination, we argue that the differential magnification of the blue and red wings of H$\alpha$ favors a flattened, virialized, low-ionization region whereas the symmetric microlensing effect measured in CIV can be reproduced by an emission line formed in a polar wind, without the need of fine-tuned caustic configurations.

COSMOGRAIL is a long-term photometric monitoring of gravitationally lensed QSOs aimed at implementing Refsdal's time-delay method to measure cosmological parameters, in particular H0. Given long and well sampled light curves of strongly lensed QSOs, time-delay measurements require numerical techniques whose quality must be assessed. To this end, and also in view of future monitoring programs or surveys such as the LSST, a blind signal processing competition named Time Delay Challenge 1 (TDC1) was held in 2014. The aim of the present paper, which is based on the simulated light curves from the TDC1, is double. First, we test the performance of the time-delay measurement techniques currently used in COSMOGRAIL. Second, we analyse the quantity and quality of the harvest of time delays obtained from the TDC1 simulations. To achieve these goals, we first discover time delays through a careful inspection of the light curves via a dedicated visual interface. Our measurement algorithms can then be applied to the data in an automated way. We show that our techniques have no significant biases, and yield adequate uncertainty estimates resulting in reduced chi2 values between 0.5 and 1.0. We provide estimates for the number and precision of time-delay measurements that can be expected from future time-delay monitoring campaigns as a function of the photometric signal-to-noise ratio and of the true time delay. We make our blind measurements on the TDC1 data publicly available

Luminosity profiles of galaxies acting as strong gravitational lenses can be tricky to study. Indeed, strong gravitational lensing images display several lensed components, both point-like and diffuse, around the lensing galaxy. Those objects limit the study of the galaxy luminosity to its inner parts. Therefore, the usual fitting methods perform rather badly on such images. Previous studies of strong lenses luminosity profiles using software such as GALFIT or IMFITFITS and various PSF-determining methods have resulted in discrepant results. The present work aims at investigating the causes of those discrepancies, as well as at designing more robust techniques for studying the morphology of early-type lensing galaxies with the ability to subtract a lensed signal from their luminosity profiles. We design a new method to independently measure each shape parameter, namely, the position angle, ellipticity, and half-light radius of the galaxy. Our half-light radius measurement method is based on an innovative scheme for computing isophotes that is well suited to measuring the morphological properties of gravitational lensing galaxies. Its robustness regarding various specific aspects of gravitational lensing image processing is analysed and tested against GALFIT. It is applied to a sample of systems from the CASTLES database. Simulations show that when restricted to small, inner parts of the lensing galaxy, the technique presented here is more trustworthy than GALFIT. It gives more robust results than GALFIT, which shows instabilities regarding the fitting region, the value of the Sersic index, and the signal-to-noise ratio. It is therefore better suited than GALFIT for gravitational lensing galaxies. It is also able to study lensing galaxies that are not much larger than the PSF. New values for the half-light radius of the objects in our sample are presented and compared to previous works.

We have obtained spectropolarimetric observations of the four images of the gravitationally lensed broad absorption line quasar H1413+117. The polarization of the microlensed image D is significantly different both in the continuum and in the broad lines from the polarization of image A which is essentially unaffected by microlensing. The observations suggest that the continuum is scattered off two regions spatially separated and producing roughly perpendicular polarizations. These results are compatible with a model in which the microlensed polarized continuum comes from a compact region located in the equatorial plane close to the accretion disk and the non-microlensed continuum from an extended region located along the polar axis.

Testing the standard Shakura-Sunyaev model of accretion is a challenging task because the central region of quasars where accretion takes place is unresolved with telescopes. The analysis of microlensing in gravitationally lensed quasars is one of the few techniques that can test this model, yielding to the measurement of the size and of temperature profile of the accretion disc. We present spectroscopic observations of the gravitationally lensed broad absorption line quasar H1413+117, which reveal partial microlensing of the continuum emission that appears to originate from two separated regions: a microlensed region, corresponding to the compact accretion disc; and a non-microlensed region, more extended and contributing to at least 30\% of the total UV-continuum flux. Because this extended continuum is occulted by the broad absorption line clouds, it is not associated with the host galaxy, but rather with light scattered in the neighbourhood of the central engine. We measure the amplitude of microlensing of the compact continuum over the rest-frame wavelength range 1000-7000 \AA. Following a Bayesian scheme, we confront our measurements to microlensing simulations of an accretion disc with a temperature varying as $T \propto R^{-1/\nu}$. We find a most likely source half-light radius of $R_{1/2} = 0.61 \times 10^{16}\,$cm (i.e., 0.002\,pc) at 0.18\,$\mu$m, and a most-likely index of $\nu=0.4$. The standard disc ($\nu=4/3$) model is not ruled out by our data, and is found within the 95\% confidence interval associated with our measurements. We demonstrate that, for H1413+117, the existence of an extended continuum in addition to the disc emission only has a small impact on the inferred disc parameters, and is unlikely to solve the tension between the microlensing source size and standard disc sizes, as previously reported in the literature.

We present VLA detections of radio emission in four four-image gravitational lens systems with quasar sources: HS0810+2554, RXJ0911+0511, HE0435$-$1223 and SDSSJ0924+0219, and e-MERLIN observations of two of the systems. The first three are detected at a high level of significance, and SDSS J0924+0219 is detected. HS0810+2554 is resolved, allowing us for the first time to achieve 10-mas resolution of the source frame in the structure of a radio quiet quasar. The others are unresolved or marginally resolved. All four objects are among the faintest radio sources yet detected, with intrinsic flux densities in the range 1-5$\mu$Jy; such radio objects, if unlensed, will only be observable routinely with the Square Kilometre Array. The observations of HS0810+2554, which is also detected with e-MERLIN, strongly suggest the presence of a mini-AGN, with a radio core and milliarcsecond scale jet. The flux densities of the lensed images in all but HE0435-1223 are consistent with smooth galaxy lens models without the requirement for smaller-scale substructure in the model, although some interesting anomalies are seen between optical and radio flux densities. These are probably due to microlensing effects in the optical.

Discrepancies between the observed and model-predicted radio flux ratios are seen in a number of quadruply-lensed quasars. The most favored interpretation of these anomalies is that CDM substructures present in lensing galaxies perturb the lens potentials and alter image magnifications and thus flux ratios. So far no consensus has emerged regarding whether or not the predicted CDM substructure abundance fully accounts for the lensing flux anomaly observations. Accurate modeling relies on a realistic lens sample in terms of both the lens environment and internal structures and substructures. In this paper we construct samples of generalised and specific lens potentials, to which we add (rescaled) subhalo populations from the galaxy-scale Aquarius and the cluster-scale Phoenix simulation suites. We further investigate the lensing effects from subhalos of masses several orders of magnitude below the simulation resolution limit. The resulting flux ratio distributions are compared to the currently best available sample of radio lenses. The observed anomalies in B0128+437, B0712+472 and B1555+375 are more likely to be caused by propagation effects or oversimplified lens modeling, signs of which are already seen in the data. Among the quadruple systems that have closely located image triplets/pairs, the anomalous flux ratios of MG0414+0534 can be reproduced by adding CDM subhalos to its macroscopic lens potential, with a probability of 5%-20%; for B0712+472, B1422+231, B1555+375 and B2045+265, these probabilities are only of a few percent. We hence find that CDM substructures are unlikely to be the whole reason for radio flux anomalies. We discuss other possible effects that might also be at work.

Lensing flux-ratio anomalies are most likely caused by gravitational lensing by small-scale dark matter structures. These anomalies offer the prospect of testing a fundamental prediction of the cold dark matter (CDM) cosmological model: the existence of numerous substructures that are too small to host visible galaxies. In two previous studies we found that the number of subhalos in the six high-resolution simulations of CDM galactic halos of the Aquarius project is not sufficient to account for the observed frequency of flux ratio anomalies seen in selected quasars from the CLASS survey. These studies were limited by the small number of halos used, their narrow range of masses (1-2E12 solar masses) and the small range of lens ellipticities considered. We address these shortcomings by investigating the lensing properties of a large sample of halos with a wide range of masses in two sets of high resolution simulations of cosmological volumes and comparing them to a currently best available sample of radio quasars. We find that, as expected, substructures do not change the flux-ratio probability distribution of image pairs and triples with large separations, but they have a significant effect on the distribution at small separations. For such systems, CDM substructures can account for a substantial fraction of the observed flux-ratio anomalies. For large close-pair separation systems, the discrepancies existing between the observed flux ratios and predictions from smooth halo models are attributed to simplifications inherent in these models which do not take account of fine details in the lens mass distributions.

We have measured the optical linear polarization of quasars belonging to Gpc-scale quasar groups at redshift z ~ 1.3. Out of 93 quasars observed, 19 are significantly polarized. We found that quasar polarization vectors are either parallel or perpendicular to the directions of the large-scale structures to which they belong. Statistical tests indicate that the probability that this effect can be attributed to randomly oriented polarization vectors is of the order of 1%. We also found that quasars with polarization perpendicular to the host structure preferentially have large emission line widths while objects with polarization parallel to the host structure preferentially have small emission line widths. Considering that quasar polarization is usually either parallel or perpendicular to the accretion disk axis depending on the inclination with respect to the line of sight, and that broader emission lines originate from quasars seen at higher inclinations, we conclude that quasar spin axes are likely parallel to their host large-scale structures.

Owing to the advent of large area photometric surveys, the possibility to use broad band photometric data, instead of spectra, to measure the size of the broad line region of active galactic nuclei, has raised a large interest. We describe here a new method using time-delay lensed quasars where one or several images are affected by microlensing due to stars in the lensing galaxy. Because microlensing decreases (or increases) the flux of the continuum compared to the broad line region, it changes the contrast between these two emission components. We show that this effect can be used to effectively disentangle the intrinsic variability of those two regions, offering the opportunity to perform reverberation mapping based on single band photometric data. Based on simulated light curves generated using a damped random walk model of quasar variability, we show that measurement of the size of the broad line region can be achieved using this method, provided one spectrum has been obtained independently during the monitoring. This method is complementary to photometric reverberation mapping and could also be extended to multi-band data. Because the effect described above produces a variability pattern in difference light curves between pairs of lensed images which is correlated with the time-lagged continuum variability, it can potentially produce systematic errors in measurement of time delays between pairs of lensed images. Simple simulations indicate that time-delay measurement techniques which use a sufficiently flexible model for the extrinsic variability are not affected by this effect and produce accurate time delays.

Using infrared spectra of the z = 1.693 quadruply lensed quasar HE0435-1223 acquired in 2009 with the spectrograph SINFONI at the ESO Very Large Telescope, we have detected a clear microlensing effect in images A and D. While microlensing affects the blue and red wings of the H{\alpha} line profile in image D very differently, it de-magnifies the line core in image A. The combination of these different effects sets constraints on the line-emitting region; these constraints suggest that a rotating ring is at the origin of the H{\alpha} line. Visible spectra obtained in 2004 and 2012 indicate that the MgII line profile is microlensed in the same way as the H{\alpha} line. Our results therefore favour flattened geometries for the low-ionization line-emitting region, for example, a Keplerian disk. Biconical models cannot be ruled out but require more fine-tuning. Flux ratios between the different images are also derived and confirm flux anomalies with respect to estimates from lens models with smooth mass distributions.

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 $\Lambda 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 mass-sheet 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 $\sim 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.

We present the results from nine years of optically monitoring the gravitationally lensed z=0.658 quasar RX J1131-1231. The R-band light curves of the four individual images of the quasar were obtained using deconvolution photometry for a total of 707 epochs. Several sharp quasar variability features strongly constrain the time delays between the quasar images. Using three different numerical techniques, we measure these delays for all possible pairs of quasar images while always processing the four light curves simultaneously. For all three methods, the delays between the three close images A, B, and C are compatible with being 0, while we measure the delay of image D to be 91 days, with a fractional uncertainty of 1.5% (1 sigma), including systematic errors. Our analysis of random and systematic errors accounts in a realistic way for the observed quasar variability, fluctuating microlensing magnification over a broad range of temporal scales, noise properties, and seasonal gaps. Finally, we find that our time-delay measurement methods yield compatible results when applied to subsets of the data.

The main obstacle for gravitational lensing to determine accurate masses of deflectors, or to determine precise estimates for the Hubble constant, is the degeneracy of lensing observables with respect to the mass-sheet transformation (MST). The MST is a global modification of the mass distribution which leaves all image positions, shapes and flux ratios invariant, but which changes the time delay. Here we show that another global transformation of lensing mass distributions exists which almost leaves image positions and flux ratios invariant, and of which the MST is a special case. Whereas for axi-symmetric lenses this source position transformation exactly reproduces all strong lensing observables, it does so only approximately for more general lens situations. We provide crude estimates for the accuracy with which the transformed mass distribution can reproduce the same image positions as the original lens model, and present an illustrative example of its performance. This new invariance transformation most likely is the reason why the same strong lensing information can be accounted for with rather different mass models.

Multiply-imaged quasars and AGNs observed in the mid-infrared (MIR) range are commonly assumed to be unaffected by the microlensing produced by the stars in their lensing galaxy. In this paper, we investigate the validity domain of this assumption. Indeed, that premise disregards microlensing of the accretion disc in the MIR range, and does not account for recent progress in our knowledge of the dusty torus. To simulate microlensing, we first built a simplified image of the quasar composed of an accretion disc, and of a larger ring-like torus. The mock quasars are then microlensed using an inverse ray-shooting code. We simulated the wavelength and size dependence of microlensing for different lensed image types and fraction of compact objects projected in the lens. This allows us to derive magnification probabilities as a function of wavelength, as well as to calculate the microlensing-induced deformation of the spectral energy distribution of the lensed images. We find that microlensing variations as large as 0.1 mag are very common at 11 microns (observer-frame). The main signal comes from microlensing of the accretion disc, which may be significant even when the fraction of flux from the disc is as small as 5 % of the total flux. We also show that the torus of sources with Lbol <~ 10^45 erg/s is expected to be noticeably microlensed. Microlensing may thus be used to get insight into the rest near-infrared inner structure of AGNs. Finally, we investigate whether microlensing in the mid-infrared can alter the so-called Rcusp relation that links the fluxes of the lensed images triplet produced when the source lies close to a cusp macro-caustic. This relation is commonly used to identify massive (dark-matter) substructures in lensing galaxies. We find that significant deviations from Rcusp may be expected, which means that microlensing can explain part of the flux ratio problem.

We report the discovery of a pair of quasars at $z=1.487$, with a separation of $8\farcs585\pm0\farcs002$. Subaru Telescope infrared imaging reveals the presence of an elliptical and a disk-like galaxy located almost symmetrically between the quasars, creating a cross-like configuration. Based on absorption lines in the quasar spectra and the colors of the galaxies, we estimate that both galaxies are located at redshift $z=0.899$. This, as well as the similarity of the quasar spectra, suggests that the system is a single quasar multiply imaged by a galaxy group or cluster acting as a gravitational lens, although the possibility of a binary quasar cannot be fully excluded. We show that the gravitational lensing hypothesis implies these galaxies are not isolated, but must be embedded in a dark matter halo of virial mass $\sim 4 \times 10^{14}\ h_{70}^{-1}\ {M}_\odot$ assuming an NFW model with a concentration parameter of $c_{vir}=6$, or a singular isothermal sphere profile with a velocity dispersion of $\sim 670$ km s$^{-1}$. We place constraints on the location of the dark matter halo, as well as the velocity dispersions of the galaxies. In addition, we discuss the influence of differential reddening, microlensing and intrinsic variability on the quasar spectra and broadband photometry.

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, (2) archival HST imaging of the lens system, (3) a new velocity-dispersion measurement of the lens galaxy of 323+/-20km/s 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. In combination with the WMAP7 data set in flat wCDM cosmology, our unblinded cosmological constraints for RXJ1131-1231 are: H_0=80.0+5.8/-5.7km/s/Mpc, OmegaDE=0.79+/-0.03 and w=-1.25+0.17/-0.21. We find the results to be statistically consistent with those from the analysis of the gravitational lens B1608+656. The joint constraints from the two lenses and WMAP7 are H_0=75.2+4.4/-4.2km/s/Mpc, OmegaDE=0.76+0.02/-0.03 and w=-1.14+0.17/-0.20 in flat wCDM, and H_0=73.1+2.4/-3.6km/s/Mpc, OmegaL=0.75+0.01/-0.02 and OmegaK=0.003+0.005/-0.006 in open LCDM. Time-delay lenses constrain especially tightly the Hubble constant (5.7% and 4.0% respectively in wCDM and open LCDM) and curvature of the universe. 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.

Microlensing is a powerful technique which can be used to study the continuum and the broad line emitting regions in distant AGNs. After a brief description of the methods and required data, we present recent applications of this technique. We show that microlensing allows one to measure the temperature profile of the accretion disc, estimate the size and study the geometry of the region emitting the broad emission lines.

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. This sample is composed of 13 pairs of previously unpublished spectra and four pairs of spectra from literature. Our analysis is based on a 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 lines in our spectra (CIII] and MgII) 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.

(abridged) Gravitationally lensed quasars can be used as powerful cosmological and astrophysical probes. We can (i) infer the Hubble constant based on the time-delay technique, (ii) unveil substructures along the l.o.s. toward distant galaxies, and (iii) compare the shape and the slope of baryons and dark matter distributions in galaxies. To reach these goals, we need high-accuracy astrometry and morphology measurements of the lens. In this work, we first present new astrometry for 11 lenses with measured time delays. Using MCS deconvolution on NIC2 HST images, we reached an astrometric accuracy of about 1-2.5 mas and an accurate shape measurement of the lens galaxy. Second, we combined these measurements with those of 14 other systems to present new mass models of these lenses. This led to the following results: 1) In 4 double-image quasars, we show that the influence of the lens environment on the time delay can easily be quantified and modeled, hence putting these lenses with high priority for time-delay determination. 2) For quadruple-image quasars, the difficulty often encountered in reproducing the image positions to milli-arcsec accuracy (astrometric anomaly) is overcome by explicitly including the nearest visible galaxy in the model. However, one anomalous system (J1131-1231) does not show any luminous perturber in its vicinity, and three others (WFI2026-4536, WFI2033-4723, and B2045+265) have problematic modeling. These 4 systems are the best candidates for a pertubation by a dark matter substructure. 3) We find a significant correlation between the PA of the light and of the mass distributions in lensing galaxies. In contrast with other studies, we find that the ellipticity of the light and of the mass also correlate well, suggesting that the overall spatial distribution of matter is not very different from the baryon distribution in the inner \sim 5 kpc of lensing galaxies.

We report the discovery of three new cases of QSOs acting as strong gravitational lenses on background emission line galaxies: SDSS J0827+5224 (zQSO = 0.293, zs = 0.412), SDSS J0919+2720 (zQSO = 0.209, zs = 0.558), SDSS J1005+4016 (zQSO = 0.230, zs = 0.441). The selection was carried out using a sample of 22,298 SDSS spectra displaying at least four emission lines at a redshift beyond that of the foreground QSO. The lensing nature is confirmed from Keck imaging and spectroscopy, as well as from HST/WFC3 imaging in the F475W and F814W filters. Two of the QSOs have face-on spiral host galaxies and the third is a QSO+galaxy pair. The velocity dispersion of the host galaxies, inferred from simple lens modeling, is between \sigma_v = 210 and 285 km/s, making these host galaxies comparable in mass with the SLACS sample of early-type strong lenses.

We present accurate time delays for the quadruply imaged quasar HE 0435-1223. The delays were measured from 575 independent photometric points obtained in the R-band between January 2004 and March 2010. With seven years of data, we clearly show that quasar image A is affected by strong microlensing variations and that the time delays are best expressed relative to quasar image B. We measured Delta_t(BC) = 7.8+/-0.8 days, Delta_t(BD) = -6.5+/-0.7 days and Delta_t_CD = -14.3+/-0.8 days. We spacially deconvolved HST NICMOS2 F160W images to derive accurate astrometry of the quasar images and to infer the light profile of the lensing galaxy. We combined these images with a stellar population fitting of a deep VLT spectrum of the lensing galaxy to estimate the baryonic fraction, $f_b$, in the Einstein radius. We measured f_b = 0.65+0.13-0.10 if the lensing galaxy has a Salpeter IMF and f_b = 0.45+0.04-0.07 if it has a Kroupa IMF. The spectrum also allowed us to estimate the velocity dispersion of the lensing galaxy, sigma_ap = 222+/-34 km/s. We used f_b and sigma_ap to constrain an analytical model of the lensing galaxy composed of an Hernquist plus generalized NFW profile. We solve the Jeans equations numerically for the model and explored the parameter space under the additional requirement that the model must predict the correct astrometry for the quasar images. Given the current error bars on f_b and sigma_ap, we did not constrain H0 yet with high accuracy, i.e., we found a broad range of models with chi^2 < 1. However, narrowing this range is possible, provided a better velocity dispersion measurement becomes available. In addition, increasing the depth of the current HST imaging data of HE 0435-1223 will allow us to combine our constraints with lens reconstruction techniques that make use of the full Einstein ring that is visible in this object.

Anomalous flux ratios between quasar images are suspected to be caused by substructures in lens galaxies. We present new deep and high resolution H and Ks imaging of the strongly lensed quasar SDSS J0924+0219 obtained using the ESO VLT with adaptive optics and the Laser Guide Star system. SDSS J0924+0219 is particularly interesting as the observed flux ratio between the quasar images vastly disagree with the predictions from smooth mass models. With our adaptive optics observations we find a luminous object, Object L, located ~0.3" to the North of the lens galaxy, but we show that it can not be responsible for the anomalous flux ratios. Object L as well as a luminous extension of the lens galaxy to the South are seen in the archival HST/ACS image in the F814W filter. This suggests that Object L is part of a bar in the lens galaxy, as also supported by the presence of a significant disk component in the light profile of the lens galaxy. Finally, we do not find evidence for any other luminous substructure that may explain the quasar images flux ratios. However, owe to the persistence of the flux ratio anomaly over time (~ 7 years) a combination of microlensing and milli-lensing is the favorite explanation for the observations.

We aim to use microlensing taking place in the lensed quasar Q2237+0305 to study the structure of the broad line region and measure the size of the region emitting the CIV and CIII] lines. Methods: Based on 39 spectrophotometric monitoring data points obtained between Oct. 2004 and Dec. 2007, we derived lightcurves for the CIV and CIII] emission lines. We used three different techniques to analyse the microlensing signal. Different components of the lines (narrow, broad and very broad) are identified and studied. We built a library of simulated microlensing lightcurves that reproduce the signal observed in the continuum and in the lines provided only the source size is changed. A Bayesian analysis scheme is then developed to derive the size of the various components of the BLR. Results: 1. The half-light radius of the region emitting the CIV line is found to be R_CIV ~ 66^{+110}_{-46} lt-days = 0.06$^{+0.09}_{-0.04}$ pc = 1.7$^{+2.8}_{-1.1}$\,10$^{17}$ cm (at 68.3% CI). Similar values are obtained for CIII]. Relative sizes of the carbon-line and V-band continuum emitting-regions are also derived with median values of R(line)/R(cont) in the range 4 to 29, depending of the FWHM of the line component. 2. The size of the CIV emitting region agrees with the Radius-Luminosity relationship derived from reverberation mapping. Using the virial theorem we derive the mass of the black hole in Q2237+0305 to be M_BH ~ 10^{8.3+/-0.3} M_sun. 3. We find that the CIV and CIII] lines are produced in at least 2 spatially distinct regions, the most compact one giving rise to the broadest component of the line. The broad and narrow line profiles are slightly different for CIV and CIII]. 4. Our analysis suggests a different structure for the CIV and FeII+III emitting regions, with the latter produced in the inner part of the BLR or in a less extended emitting region than CIV.

We present new optical circular polarization measurements with typical uncertainties < 0.1% for a sample of 21 quasars. All but two objects have null circular polarization. We use this result to constrain the polarization due to photon-pseudoscalar mixing along the line of sight. We detect significant (> 3 sigma) circular polarization in two blazars with high linear polarization and discuss the implications of this result for quasar physics. In particular, the recorded polarization degrees may be indicative of magnetic fields as strong as 1 kG or a significant contribution of inverse Compton scattering to the optical continuum.

We apply the iterative MCS deconvolution method (ISMCS) to near-IR HST archives data of seven gravitationally lensed quasars currently monitored by the COSMOGRAIL collaboration: HE 0047-1756, RX J1131-1231, SDSS J1138+0314, SDSS J1155+6346, SDSS J1226-0006, WFI J2026-4536 and HS 2209+1914. In doing so, we obtain relative positions for the lensed images and shape parameters for the light distribution of the lensing galaxy in each system. The lensed image positions are derived with 1-2 mas accuracy. To predict time delays and to test the ability of simple mass models to reproduce the observed configuration, isothermal and de Vaucouleurs mass models are calculated for the whole sample using state-of-the-art modeling techniques. The effect of the lens environment on the lens mass models is taken into account with a shear term. Doubly imaged quasars are equally well fitted by each of these models. A large amount of shear is necessary to reproduce SDSS J1155+6346 and SDSS J1226-006. In the latter case, we identify a nearby galaxy as the dominant source of shear. The quadruply imaged quasar SDSS J1138+0314 is well reproduced by simple lens models, which is not the case for the two other quads, RX J1131-1231 and WFI J2026-4536. This might be the signature of astrometric perturbations due to massive substructures in the lensing galaxy unaccounted for by the models. Other possible explanations are also presented.

We present the first case of strong gravitational lensing by a QSO : SDSS J0013+1523, at z = 0.120. The discovery is the result of a systematic search for emission lines redshifted behind QSOs, among 22298 spectra of the SDSS data release 7. Apart from the z = 0.120 spectral features of the foreground QSO, the spectrum of SDSS J0013+1523 also displays the OII and Hbeta emission lines and the OIII doublet, all at the same redshift, z = 0.640. Using sharp Keck adaptive optics K-band images obtained using laser guide stars, we unveil two objects within a radius of 2 arcsec from the QSO. Deep Keck optical spectroscopy clearly confirms one of these objects at z = 0.640 and shows traces of the OIII, emission line of the second object, also at z = 0.640. Lens modeling suggests that they represent two images of the same z = 0.640 emission-line galaxy. Our Keck spectra also allow us to measure the redshift of an intervening galaxy at z = 0.394, located 3.2 arcsec away from the line of sight to the QSO. If the z = 0.120 QSO host galaxy is modeled as a singular isothermal sphere, its mass within the Einstein radius is M_E(r < 1 kpc) = 2.16e10 M_Sun and its velocity dispersion is sigma_SIS = 169 km/s. This is about 1 sigma away from the velocity dispersion estimated from the width of the QSO Hbeta emission line, sigma_*(M_BH) = 124 +/- 47 km/s. Deep optical HST imaging will be necessary to constrain the total radial mass profile of the QSO host galaxy using the detailed shape of the lensed source. This first case of a QSO acting as a strong lens on a more distant object opens new directions in the study of QSO host galaxies.

On the basis of 16 years of spectroscopic observations of the four components of the gravitationally lensed broad absorption line (BAL) quasar H1413+117, covering the ultraviolet to visible rest-frame spectral range, we analyze the spectral differences observed in the P Cygni-type line profiles and have used the microlensing effect to derive new clues to the BAL profile formation. We confirm that the spectral differences observed in component D can be attributed to a microlensing effect lasting at least a decade. We show that microlensing magnifies the continuum source in image D, leaving the emission line region essentially unaffected. We interpret the differences seen in the absorption profiles of component D as the result of an emission line superimposed onto a nearly black absorption profile. We also find that the continuum source and a part of the broad emission line region are likely de-magnified in component C, while components A and B are not affected by microlensing. We show that microlensing of the continuum source in component D has a chromatic dependence compatible with the thermal continuum emission of a standard Shakura-Sunyaev accretion disk. Using a simple decomposition method to separate the part of the line profiles affected by microlensing and coming from a compact region from the part unaffected by this effect and coming from a larger region, we disentangle the true absorption line profiles from the true emission line profiles. The extracted emission line profiles appear double-peaked, suggesting that the emission is occulted by a strong absorber, narrower in velocity than the full absorption profile, and emitting little by itself. We propose that the outflow around H1413+117 is constituted by a high-velocity polar flow and a denser, lower velocity disk seen nearly edge-on.

Given an observed gravitational lens mirage produced by a foreground deflector (cf. galaxy, quasar, cluster,...), it is possible via numerical lens inversion to retrieve the real source image, taking full advantage of the magnifying power of the cosmic lens. This has been achieved in the past for several remarkable gravitational lens systems. Instead, we propose here to invert an observed multiply imaged source directly at the telescope using an ad-hoc optical instrument which is described in the present paper. Compared to the previous method, this should allow one to detect fainter source features as well as to use such an optimal gravitational lens telescope to explore even fainter objects located behind and near the lens. Laboratory and numerical experiments illustrate this new approach.

We present the first VLT near-IR observations of a gravitationally lensed quasar, using adaptive optics and laser guide star. These observations can be considered as a test bench for future systematic observations of lensed quasars with adaptive optics, even when bright natural guide stars are not available in the nearby field. With only 14 minutes of observing time, we derived very accurate astrometry of the quasar images and of the lensing galaxy, with 0.05 \arcsec spatial resolution, comparable to the Hubble Space Telescope (HST). In combination with deep VLT optical spectra of the quasar images, we use our adaptive optics images to constrain simple models for the mass distribution of the lensing galaxy. The latter is almost circular and does not need any strong external shear to fit the data. The time delay predicted for SDSS0806+2006, assuming a singular isothermal ellipsoid model and the concordance cosmology, is Delta t \simeq 50 days. Our optical spectra indicate a flux ratio between the quasar images of A/B=1.3 in the continuum and A/B=2.2 in both the MgII and in the CIII] broad emission lines. This suggests that microlensing affects the continuum emission. However, the constant ratio between the two emission lines indicates that the broad emission line region is not microlensed. Finally, we see no evidence of reddening by dust in the lensing galaxy.

Soft gamma-ray repeaters (SGRs) are a rare type of gamma-ray transient sources that are ocasionally detected as bursts in the high-energy sky. They are thought to be produced by magnetars, young neutron stars with very strong magnetic fields of the order of 10^(14-15) G. Only three such objects are known in our Galaxy, and a fourth one is associated with the supernova remnant N49 in the Large Magellanic Cloud. In none of these cases has an optical counterpart to either the gamma-ray flares or the quiescent source been identified. Here we present multi-wavelength observations of a puzzling source, SWIFT J195509+261406, for which we detected more than 40 flaring episodes in the optical band over a time span of 3 days, plus a faint infrared flare 11 days later, after which it returned to quiescence. We propose that SWIFT J195509+261406 is a member of a subgroup of SGRs for which the long-term X-ray emission is transient in nature. Furthermore, it is the first SGR for which bursts have been detected in the optical and near-infrared bands and maybe the link between the "persistent" SGRs and the dim isolated neutron stars.

Based on a sample of 355 quasars with significant optical polarization, we found that quasar polarization vectors are not randomly oriented over the sky as naturally expected. The probability that the observed distribution of polarization angles is due to chance is lower than 0.1%. The polarization vectors of the light from quasars are aligned although the sources span huge regions of the sky (~ 1 Gpc). Groups of quasars located along similar lines of sight but at different redshifts (typically z ~ 0.5 and z ~ 1.5) are characterized by different preferred directions of polarization. These characteristics make the observed alignment effect difficult to explain in terms of a local contamination by interstellar polarization in our Galaxy. Interpreted in terms of a cosmological-size effect, we show that the dichroism and birefringence predicted by a mixing between photons and very light pseudoscalar particles within a magnetic field can qualitatively reproduce the observations. We find that circular polarization measurements could help constrain this mechanism.

We present the main results of the first long-term spectrophotometric monitoring of the ``Einstein cross'' Q2237+0305 and of the single-epoch spectra of the lensed quasar J1131-1231. From October 2004 to December 2006, we find that two prominent microlensing events affect images A & B in Q2237+0305 while images C & D remain grossly unaffected by microlensing on a time scale of a few months. Microlensing in A & B goes with chromatic variations of the quasar continuum. We observe stronger micro-amplification in the blue than in the red part of the spectrum, as expected for continuum emission arising from a standard accretion disk. Microlensing induced variations of the CIII] emission are observed both in the integrated line intensity and profile. Finally, we also find that images C & D are about 0.1-0.3 mag redder than images A & B. The spectra of images A-B-C in J1131-1231 reveal that, in April 2003, microlensing was at work in images A and C. We find that microlensing de-amplifies the continuum emission and the Broad Line Region (BLR) in these images. Contrary to the case of Q2237+0305, we do not find evidence for chromatic microlensing of the continuum emission. On the other hand, we observe that the Balmer and MgII broad line profiles are deformed by microlensing. These deformations imply an anti-correlation between the width of the emission line and the size of the corresponding emitting region. Finally, the differential microlensing of the FeII emission suggests that the bulk of FeII is emitted in the outer parts of the BLR while another fraction of FeII is produced in a compact region.

Gravitationally lensed quasars can be used to map the mass distribution in lensing galaxies and to estimate the Hubble constant H0 by measuring the time delays between the quasar images. Here we report the measurement of two independent time delays in the quadruply imaged quasar WFI J2033-4723 (z = 1.66). Our data consist of R-band images obtained with the Swiss 1.2 m EULER telescope located at La Silla and with the 1.3 m SMARTS telescope located at Cerro Tololo. The light curves have 218 independent epochs spanning 3 full years of monitoring between March 2004 and May 2007, with a mean temporal sampling of one observation every 4th day. We measure the time delays using three different techniques, and we obtain Dt(B-A) = 35.5 +- 1.4 days (3.8%) and Dt(B-C) = 62.6 +4.1/-2.3 days (+6.5%/-3.7%), where A is a composite of the close, merging image pair. After correcting for the time delays, we find R-band flux ratios of F_A/F_B = 2.88 +- 0.04, F_A/F_C = 3.38 +- 0.06, and F_A1/F_A2 = 1.37 +- 0.05 with no evidence for microlensing variability over a time scale of three years. However, these flux ratios do not agree with those measured in the quasar emission lines, suggesting that longer term microlensing is present. Our estimate of H0 agrees with the concordance value: non-parametric modeling of the lensing galaxy predicts H0 = 67 +13/-10 km s-1 Mpc-1, while the Single Isothermal Sphere model yields H0 = 63 +7/-3 km s-1 Mpc-1 (68% confidence level). More complex lens models using a composite de Vaucouleurs plus NFW galaxy mass profile show twisting of the mass isocontours in the lensing galaxy, as do the non-parametric models. As all models also require a significant external shear, this suggests that the lens is a member of the group of galaxies seen in field of view of WFI J2033-4723.

We present the results of the first long-term (2.2 years) spectroscopic monitoring of a gravitationally lensed quasar, namely the Einstein Cross Q2237+0305. We spatially deconvolve deep VLT/FORS1 spectra to accurately separate the spectrum of the lensing galaxy from the spectra of the quasar images. Accurate cross-calibration of the observations at 31 epochs from October 2004 to December 2006 is carried out using foreground stars observed simultaneously with the quasar. The quasar spectra are further decomposed into a continuum component and several broad emission lines. We find prominent microlensing events in the quasar images A and B, while images C and D are almost quiescent on a timescale of a few months. The strongest variations are observed in the continuum, and their amplitude is larger in the blue than in the red, consistent with microlensing of an accretion disk. Variations in the intensity and profile of the broad emission lines are also reported, most prominently in the wings of the CIII] and in the center of the CIV emission lines. During a strong microlensing episode observed in quasar image A, the broad component of the CIII] is more magnified than the narrow component. In addition, the emission lines with higher ionization potentials are more magnified than the lines with lower ionization potentials, consistent with the stratification of the broad line region (BLR) infered from reverberation mapping observations.

Thanks to its sharp view, HST has significantly improved our knowledge of tens of gravitationally lensed quasars in four different respects: (1) confirming their lensed nature; (2) detecting the lensing galaxy responsible for the image splitting; (3) improving the astrometric accuracy on the positions of the unresolved QSO images and of the lens; (4) resolving extended lensed structures from the QSO hosts into faint NIR or optical rings or arcs. These observations have helped to break some degeneracies on the lens potential, to probe the galaxy evolution and to reconstruct the true shape of the QSO host with an increased angular resolution.

We present the results of the first long-term (2.2 years) spectroscopic monitoring of a gravitationally lensed quasar, namely the Einstein Cross Q2237+0305. The goal of this paper is to present the observational facts to be compared in follow-up papers with theoretical models to constrain the inner structure of the source quasar. We spatially deconvolve deep VLT/FORS1 spectra to accurately separate the spectrum of the lensing galaxy from the spectra of the quasar images. Accurate cross-calibration of the 58 observations at 31-epoch from October 2004 to December 2006 is carried out with non-variable foreground stars observed simultaneously with the quasar. The quasar spectra are further decomposed into a continuum component and several broad emission lines to infer the variations of these spectral components. We find prominent microlensing events in the quasar images A and B, while images C and D are almost quiescent on a timescale of a few months. The strongest variations are observed in the continuum of image A. Their amplitude is larger in the blue (0.7 mag) than in the red (0.5 mag), consistent with microlensing of an accretion disk. Variations in the intensity and profile of the broad emission lines are also reported, most prominently in the wings of the CIII] and center of the CIV emission lines. During a strong microlensing episode observed in June 2006 in quasar image A, the broad component of the CIII] is more highly magnified than the narrow component. In addition, the emission lines with higher ionization potentials are more magnified than the lines with lower ionization potentials, consistent with the results obtained with reverberation mapping. Finally, we find that the V-band differential extinction by the lens, between the quasar images, is in the range 0.1-0.3 mag.

We present the analysis of single epoch long slit spectra of the three brightest images of the gravitationally lensed system J1131-1231. These spectra provide one of the clearest observational evidence for differential micro-lensing of broad emission lines (BELs) in a gravitationally lensed quasar. The micro-lensing effect enables us: (1) to confirm that the width of the emission lines is anti-correlated to the size of the emitting region; (2) to show that the bulk of FeII is emitted in the outer parts of the Broad Line Region (BLR) while another fraction of FeII is produced in a compact region; (3) to derive interesting informations on the origin of the narrow intrinsic MgII absorption doublet observed in that system.

(ABRIDGED) Aims: We discuss and characterize micro-lensing among the 3 brightest lensed images (A-B-C) of the gravitational lens system RXS J1131-1231 (a quadruply imaged AGN) by means of long slit optical and NIR spectroscopy. Qualitative constraints on the size of different emission regions are derived. Methods: We decompose the spectra into their individual emission components using a multi-component fitting approach. A complementary decomposition of the spectra enables us to isolate the macro-lensed fraction of the spectra independently of any spectral modelling. Results: -1. The data support micro-lensing de-amplification of images A and C. Not only is the continuum emission microlensed in those images but also a fraction of the Broad Line emitting Region (BLR).-2. Micro-lensing of a very broad component of MgII emission line suggests that the corresponding emission occurs in a region more compact than the other components of the emission line. -3. We find evidence that a large fraction of the FeII emission arises in the outer parts of the BLR. We also find very compact emitting region in the ranges 3080-3540 A and 4630-4800 A that is likely associated with FeII. -4. The OIII narrow emission line regions are partly spatially resolved. This enables us to put a lower limit of 110h^{-1} pc on their intrinsic size. -5. Analysis of MgII absorption found in the spectra indicates that the absorbing medium is intrinsic to the quasar, has a covering factor of 20%, and is constituted of small clouds homogeneously distributed in front of the continuum and BLRs. -6. Two neighbour galaxies are detected at redshifts z=0.10 and z=0.289. These galaxies are possible members of galaxy groups reported at those redshifts.

Aims: Our aim is to measure the time delay between the two gravitationally lensed images of the z = 1.547 quasar SDSS J1650+4251, in order to estimate the Hubble constant H_0. Methods: Our measurement is based on R-band light curves with 57 epochs obtained at Maidanak Observatory, in Uzbekistan, from May 2004 to September 2005. The photometry is performed using simultaneous deconvolution of the data, which provides the individual light curves of the otherwise blended quasar images. The time delay is determined from the light curves using two very different numerical techniques, i.e., polynomial fitting and direct cross-correlation. The time delay is converted into H_0 following analytical modeling of the potential well. Results: Our best estimate of the time delay is Dt = 49.5 +/- 1.9 days, i.e., we reach a 3.8% accuracy. The R-band flux ratio between the quasar images, corrected for the time delay and for slow microlensing, is F_A /F_B = 6.2 +/- 5%. Conclusions: The accuracy reached on the time delay allows us to discriminate well between families of lens models. As for most other multiply imaged quasars, only models of the lensing galaxy that have a de Vaucouleurs mass profile plus external shear give a Hubble constant compatible with the current most popular value (H_0 = 72 +/- 8 km s-1 Mpc-1). A more realistic singular isothermal sphere model plus external shear gives H_0 = 51.7 +4.0 -3.0 km s-1 Mpc-1.

We report the discovery of the first known probable case of a physical triple quasar (not a gravitational lens). A previously known double system, QQ 1429-008 at z = 2.076, is shown to contain a third, fainter QSO component at the same redshift within the measurement errors. Deep optical and IR imaging at the Keck and VLT telescopes has failed to reveal a plausible lensing galaxy group or a cluster, and moreover, we are unable to construct any viable lensing model which could lead to the observed distribution of source positions and relative intensities of the three QSO image components. Furthermore, there are hints of differences in broad-band spectral energy distributions of different components, which are more naturally understood if they are physically distinct AGN. Therefore, we conclude that this system is most likely a physical triple quasar, the first such close QSO grouping known at any redshift. The projected component separations in the restframe are ~ 30 - 50 kpc for the standard concordance cosmology, typical of interacting galaxy systems. The existence of this highly unusual system supports the standard picture in which galaxy interactions lead to the onset of QSO activity.

(Abridged) We present our VLT/FORS1 deep spectroscopic observations of the gravitationally lensed quasar SDSS J0924+0219, as well as archival HST/NICMOS and ACS images of the same object. The two-epoch spectra, obtained in the Multi Object Spectroscopy (MOS) mode, allow for very accurate flux calibration, spatial deconvolution of the data, and provide the redshift of the lensing galaxy z=0.394 +/- 0.001. These spectra, taken 15 days apart, show only slight continuum variations, while the broad emission lines display obvious changes in the red wing of the Mg II line, in the Fe II bands, and in the central part of the C III] line. Even though variations in the line profiles are present, we do not see any significant differences between the continuum and emission line flux ratios of images A and B of the quasar. Spatial deconvolution of the HST images reveals a double Einstein ring. One ring corresponds to the lensed quasar host galaxy at z=1.524 and a second bluer one, is the image either of a star-forming region in the host galaxy, or of another unrelated lower redshift object. We find that a broad range of lens models gives a satisfactory fit to the data. However, they predict very different time delays, making SDSS J0924+0219 an object of particular interest for photometric monitoring. In addition, the lens models reconstructed using exclusively the constraints from the Einstein rings, or using exclusively the astrometry of the quasar images, are not compatible. This suggests that substructures play an important role in SDSS J0924+0219.

High angular resolution images of the complex gravitational lens system RXS J1131-1231 (a quadruply imaged AGN with a bright Einstein ring) obtained with the Advanced Camera for Surveys and NICMOS instruments onboard the Hubble Space Telescope are analysed to determine the lens model and to reconstruct the host galaxy. Results: 1- Precise astrometry and photometry of the four QSO lensed images (A-D) and of the lensing galaxy (G) are obtained. They are found in agreement with an independent study presented in a companion paper. The position and colours of the X object seen in projection close to the lens are found to be only compatible with a satellite galaxy associated with the lens. 2- The Singular Isothermal Ellipsoid plus external shear provides a good fit of the astrometry of images A-D. The positions of extended substructures are also well reproduced. However an octupole (m=4) must be added to the lens potential in order to reproduce the observed lens position, as well as the IB/IC point-like image flux ratio. The ellipticity and orientation of the mass quadrupole are found similar to those of the light distribution, fitted by a Sersic profile. The lens (z=0.295) is found to be a massive elliptical in a rich environment and showing possible evolution with respect to z=0. 3- The host galaxy (z=0.658) is found to be a substantially magnified (M ~ 9) luminous Seyfert 1 spiral galaxy. The angular resolution is sufficient to see regions where stars are intensively forming. Interaction with a closeby companion is also observed. 4- Finally, in the case of RXS J1131-1231, extended lensed structures do not help much in constraining the lens model.

(abridged) Aims: RXS J113155.4-123155 (z=0.66) is a quadruply imaged lensed quasar with a resolved Einstein Ring. The goal of this paper is to provide a full characterization of this system, and more particularly accurate astrometry and photometry. Method: Visible and near-infrared imaging observations of RXS J113155.4-123155 were carried out at various epochs using several ground based telescopes and the HST. The frames have been deconvolved using the MCS algorithm. A Singular Isothermal Ellipsoid (SIE) + external shear has been used to model the lensing galaxy potential. Results: MCS deconvolution enables us to separate the flux of the QSO (point-like images) from that of its host galaxy and to accurately track the flux variations of the point-like images in various filters. The deconvolved frames unveil several multiply imaged structures in the Einstein ring and an unidentified object in the vicinity of the lensing galaxy. We discuss the lightcurves and the chromatic flux ratio variations and deduce that both intrinsic variability and microlensing took place during a span longer than one year. We demonstrate that microlensing may easily account for the so called anomalous flux ratios presented in the discovery paper. However, the observed flux ratios are still poorly reproduced when modeling the lens potential with a SIE+shear. We argue that this disagreement can hardly be explained by milli-lensing caused by substructures in the lensing galaxy. A solution proposed in Paper II consists in a more complex lens model including an octupole term to the lens gravitational potential.

Aims: To predict time delays for a sample of gravitationally lensed quasars and to evaluate the accuracy that can be realistically achieved on the value of H_0. Methods: We consider 14 lensed quasars that are candidates for time-delay monitoring and model them in detail using pixelized lens models. For each system, we provide a mass map, arrival-time surface and the distribution of predicted time-delays in a concordance cosmology, assuming H_0^{-1}=14 Gyr (H_0=70 in local units). Based on the predicted time-delays and on the observational circumstances, we rate each lens as `excellent' or `good' or `unpromising' for time-delay monitoring. Finally, we analyze simulated time delays for the 11 lens rated excellent or good, and show that H_0 can be recovered to a precision of 5%. Results: In combination with COSMOGRAIL paper I on the temporal sampling of lensed quasar light curves, the present work will help design monitoring campaigns of lensed quasars.

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 \bar{\theta} appears to rotate with redshift at the rate of ~ 30\degr per Gpc. The symmetry of the the \bar{\theta} -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 CMB 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.

New linear polarization measurements (mainly in the V band) are presented for 203 quasi-stellar objects (QSOs). The sample is made up of 94 QSOs located in the North Galactic Pole (NGP) region and of 109 QSOs in the South Galactic Pole (SGP) region. First time measurements have been obtained for 184 QSOs. Among them, 109 known radio-emitters, 42 known Broad Absorption Line (BAL) QSOs, and 1 gravitationally lensed quasi-stellar object. We found high polarization levels (p > 3%) for 12 QSOs, including the BAL QSO SDSS J1409+0048. For 10 objects, measurements obtained at different epochs do exist. Two of them show evidence for variability: the highly polarized BL Lac candidate PKS 1216-010 and the radio source PKS 1222+037

A survey measuring quasar polarization vectors has been started in two regions towards the North and South Galactic Poles. Here, We review the discovery of significant correlations of orientations of polarization vectors over huge angular distances. We report new results including a larger sample of the quasars confirming the existence of coherent orientations at redshifts z>1.

We report the discovery of a new quadruply imaged quasar surrounded by an optical Einstein ring candidate. Spectra of the different components of 1RXS J113155.4-123155 reveal a source at z=0.658. Up to now, this object is the closest known gravitationally lensed quasar. The lensing galaxy is clearly detected. Its redshift is measured to be z=0.295. Additionally, the total V magnitude of the system has varied by 0.3 mag between two epochs separated by 33 weeks. The measured relative astrometry of the lensed images is best fitted with an SIS model plus shear. This modeling suggests very high magnification of the source (up to 50 for the total magnification) and predicts flux ratios between the lensed images significantly different from what is actually observed. This suggests that the lensed images may be affected by a combination of micro or milli-lensing and dust extinction effects.