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 spectroscopic confirmation of three new two-image gravitationally lensed quasars, compiled from existing strong lens and X-ray catalogs. Images of HSC J091843.27-022007.5 show a red galaxy with two blue point sources at either side, separated by 2.26 arcsec. This system has a source and a lens redshifts $z_s = 0.804$ and $z_l = 0.459$, respectively, as obtained by our follow-up spectroscopic data. CXCO J100201.50+020330.0 shows two point sources separated by 0.85 arcsec on either side of an early-type galaxy. The follow-up spectroscopic data confirm the fainter quasar has the same redshift with the brighter quasar from the SDSS fiber spectrum at $z_s = 2.016$. The deflecting foreground galaxy is a typical early-type galaxy at a redshift of $z_l = 0.439$. SDSS J135944.21+012809.8 has two point sources with quasar spectra at the same redshift $z_s = 1.096$, separated by 1.05 arcsec, and fits to the HSC images confirm the presence of a galaxy between these. These discoveries demonstrate the power of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP)'s deep imaging and wide sky coverage. Combined with existing X-ray source catalogues and follow-up spectroscopy, the HSC-SSP provides us unique opportunities to find multiple-image quasars lensed by a foreground galaxy.

We use KiDS weak lensing data to measure variations in mean halo mass as a function of several key galaxy properties (namely: stellar colour, specific star formation rate, Sersic index, and effective radius) for a volume-limited sample of GAMA galaxies in a narrow stellar mass range ($M_* \sim 2$-$5 \times 10^{10}$ Msol). This mass range is particularly interesting, inasmuch as it is where bimodalities in galaxy properties are most pronounced, and near to the break in both the galaxy stellar mass function and the stellar-to-halo mass relation (SHMR). In this narrow mass range, we find that both size and Sersic index are better predictors of halo mass than colour or SSFR, with the data showing a slight preference for Sersic index. In other words, we find that mean halo mass is more strongly correlated with galaxy structure than either stellar populations or star formation rate. These results lead to an approximate lower bound on the dispersion in halo masses among $ \log M_* \approx {10.5}$ galaxies: we find that the dispersion is $\gtrsim 0.3$ dex. This would imply either that offsets from the mean SHMR are very tightly coupled to size/structure, or that the dispersion in the SHMR is larger than past results have suggested. Our results thus provide new empirical constraints on the relationship between stellar and halo mass assembly at this particularly interesting mass range.

The evolution of the structural and kinematic properties of early-type galaxies (ETGs), their scaling relations, as well as their stellar metallicity and age contain precious information on the assembly history of these systems. We present results on the evolution of the stellar mass-velocity dispersion relation of ETGs, focusing in particular on the effects of some selection criteria used to define ETGs. We also try to shed light on the role that in-situ and ex-situ stellar populations have in massive ETGs, providing a possible explanation of the observed metallicity distributions.

We report the largest sample of candidate strong gravitational lenses belonging to the Survey of Gravitationally-lensed Objects in HSC Imaging for group-to-cluster scale (SuGOHI-c) systems. These candidates are compiled from the S18A data release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Survey. We visually inspect $\sim39,500$ galaxy clusters, selected from several catalogs, overlapping with the Wide, Deep, and UltraDeep fields, spanning the cluster redshift range $0.05<z_{cl}<1.38$. We discover 641 candidate lens systems, of which 536 are new. From the full sample, 47 are almost certainly bonafide lenses, 181 of them are highly probable lenses and 413 are possible lens systems. Additionally, we present 131 lens candidates at galaxy-scale serendipitously discovered during the inspection. We obtained spectroscopic follow-up of 10 candidates using the X-shooter. With this follow-up, we confirm 8 systems as strong gravitational lenses. Of the remaining two, one of the sources is too faint to detect any emission, and the other has a tentative redshift close to the lens redshift, but additional arcs in this system are yet to be observed spectroscopically. Since the HSC-SSP is an ongoing survey, we expect to find $\sim600$ definite or probable lenses using this procedure and even more if combined with other lens finding methods.

Strong lenses are extremely useful probes of the distribution of matter on galaxy and cluster scales at cosmological distances, but are rare and difficult to find. The number of currently known lenses is on the order of 1,000. We wish to use crowdsourcing to carry out a lens search targeting massive galaxies selected from over 442 square degrees of photometric data from the Hyper Suprime-Cam (HSC) survey. We selected a sample of $\sim300,000$ galaxies with photometric redshifts in the range $0.2 < z_{phot} < 1.2$ and photometrically inferred stellar masses $\log{M_*} > 11.2$. We crowdsourced lens finding on this sample of galaxies on the Zooniverse platform, as part of the Space Warps project. The sample was complemented by a large set of simulated lenses and visually selected non-lenses, for training purposes. Nearly 6,000 citizen volunteers participated in the experiment. In parallel, we used YattaLens, an automated lens finding algorithm, to look for lenses in the same sample of galaxies. Based on a statistical analysis of classification data from the volunteers, we selected a sample of the most promising $\sim1,500$ candidates which we then visually inspected: half of them turned out to be possible (grade C) lenses or better. Including lenses found by YattaLens or serendipitously noticed in the discussion section of the Space Warps website, we were able to find 14 definite lenses, 129 probable lenses and 581 possible lenses. YattaLens found half the number of lenses discovered via crowdsourcing. Crowdsourcing is able to produce samples of lens candidates with high completeness and purity, compared to currently available automated algorithms. A hybrid approach, in which the visual inspection of samples of lens candidates pre-selected by discovery algorithms and/or coupled to machine learning is crowdsourced, will be a viable option for lens finding in the 2020s.

We report the serendipitous discovery of HSC J0904$-$0102, a quadruply-lensed Lyman break galaxy (LBG) in the Survey of Gravitationally-lensed Objects in Hyper Suprime-Cam Imaging (SuGOHI). Owing to its point-like appearance, the source was thought to be a lensed active galactic nucleus. We obtained follow-up spectroscopic data with the Gemini Multi-Object Spectrographs on the Gemini South Telescope, which confirmed this to be a lens system. The deflecting foreground galaxy is a typical early-type galaxy at a high redshift of $z_{\ell} = 0.957$ with stellar velocity dispersion $\sigma_v=259\pm56$ km~s$^{-1}$. The lensed source is identified as an LBG at $z_{\rm s} = 3.403$, based on the sharp drop bluewards of Ly$\alpha$ and other absorption features. A simple lens mass model for the system, assuming a singular isothermal ellipsoid, yields an Einstein radius of $\theta_{\rm Ein} = 1. 23^{\prime\prime}$ and a total mass within the Einstein radius of $M_{\rm Ein} = (5.55\pm 0.24) \times 10^{11}M_{\odot}$ corresponding to a velocity dispersion of $\sigma_{\rm SIE}= 283\pm 3$ km~s$^{-1}$, which is in good agreement with the value derived spectroscopically. The most isolated lensed LBG image has a magnification of $\sim 6.5$. In comparison with other lensed LBGs and typical $z\sim4$ LBG populations, HSC J0904$-$0102 is unusually compact, an outlier at $>2\sigma$ confidence. Together with a previously discovered SuGOHI lens, HSC J1152$+$0047, that is similarly compact, we believe that the HSC Survey is extending LBG studies down to smaller galaxy sizes.

The stellar density profile a galaxy is typically summarised with two numbers: total stellar mass and half-light radius. The total mass of a galaxy, however, is not a well-defined quantity, due to the finite depth of photometric observations and the arbitrariness of the distinction between galaxy and diffuse intra-group light. This limits our ability to make accurate comparisons between models and observations. I wish to provide a more robust two-parameter description of the stellar density distribution of elliptical galaxies, in terms of quantities that can be measured unambiguously. I propose to use the stellar mass enclosed within 10 kpc in projection, $M_{*,10}$, and the mass-weighted stellar density slope within the same aperture, $\Gamma_{*,10}$, for this purpose. I measured the distribution in $M_{*,10}$ and $\Gamma_{*,10}$ of a sample of elliptical galaxies from the Sloan Digital Sky Survey and the Galaxy And Mass Assembly survey, using photometry from the Hyper Suprime-Cam survey. The pair of values of $(M_{*,10},\Gamma_{*,10})$ can be used to predict the stellar density profile in the inner 10 kpc of a galaxy with better than 20% accuracy. Similarly, $M_{*,10}$ and $\Gamma_{*,10}$ can be combined to obtain a proxy for stellar velocity dispersion at least as good as the stellar mass fundamental plane. As a first application, I then compared the observed $M_{*,10}-\Gamma_{*,10}$ relation of elliptical galaxies with that of similarly selected galaxies in the EAGLE Reference simulation. Observed and simulated galaxies match at $M_{*,10}=10^{11}M_\odot$, but the EAGLE $M_{*,10}-\Gamma_{*,10}$ relation is shallower and has a larger intrinsic scatter compared to observations. This new parameterisation of the stellar density profile of massive elliptical galaxies provides a more robust way of comparing results from different photometric surveys and from hydrodynamical simulations.

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).

A double source plane (DSP) system is a precious probe for the density profile of distant galaxies and cosmological parameters. However, these measurements could be affected by the surrounding environment of the lens galaxy. Thus, it is important to evaluate the cluster-scale mass for detailed mass modeling. We observed the {\it Eye of Horus}, a DSP system discovered by the Subaru HSC--SSP, with XMM--Newton. We detected two X-ray extended emissions, originating from two clusters, one centered at the {\it Eye of Horus}, and the other located $\sim100$ arcsec northeast to the {\it Eye of Horus}. We determined the dynamical mass assuming hydrostatic equilibrium, and evaluated their contributions to the lens mass interior of the Einstein radius. The contribution of the former cluster is $1.1^{+1.2}_{-0.5}\times10^{12}~M_{\odot}$, which is $21-76\%$ of the total mass within the Einstein radius. The discrepancy is likely due to the complex gravitational structure along the line of sight. On the other hand, the contribution of the latter cluster is only $\sim2\%$ on the {\it Eye of Horus}. Therefore, the influence associated with this cluster can be ignored.

Strong gravitationally lensed quasars provide powerful means to study galaxy evolution and cosmology. We use Chitah to hunt for new lens systems in the Hyper Suprime$-$Cam Subaru Strategic Program (HSC SSP) S16A. We present 46 lens candidates, of which 3 are previously known. Including 2 additional lenses found by YattaLens, we obtain X-shooter spectra of 6 promising candidates for lens confirmation and redshift measurements. We report new spectroscopic redshift measurements for both the lens and source galaxies in 4 lens systems. We apply the lens modeling software Glee to model our 6 X-shooter lenses uniformly. Through our analysis of the HSC images, we find that HSCJ022622$-$042522, HSCJ115252$+$004733, and HSCJ141136$-$010216 have point-like lensed images, and that the lens light distribution is well aligned with mass distribution within 6 deg. Thanks to the X-shooter spectra, we estimate fluxes on the Baldwin- Phillips-Terlevich (BPT) diagram, and find that HSCJ022622$-$042522 has a probable quasar source, based on the upper limit of the Nii flux intensity. We also measure the FWHM of Ly$\alpha$ emission of HSCJ141136$-$010216 to be $\sim$254 km/s, showing that it is a probable Lyman-$\alpha$ emitter.

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.

Taking advantage of a Bayesian hierarchical inference formalism, we study the evolution of the observed correlation between central stellar velocity dispersion $\sigma_0$ and stellar mass $M_*$ of massive ($M_*\gtrsim 10^{10.5}\,\mathrm{M_\odot}$) early-type galaxies (ETGs) out to redshift $z\approx 2.5$. Collecting ETGs from state-of-the-art literature samples, we consider both a fiducial ($0\lesssim z\lesssim 1.2$) and an extended ($0\lesssim z\lesssim 2.5$) samples. Based on the fiducial sample, we find that at $z\lesssim 1.2$ the $M_*-\sigma_0$ relation is well represented by $\sigma_0\propto M_*^{\beta}(1+z)^{\zeta}$, with $\beta\simeq 0.23$ (at given $z$, $\sigma_0$ increases with $M_*$ as a power law with slope similar to the classical Faber-Jackson relation) and $\zeta\simeq 0.26$ (at given $M_*$, $\sigma_0$ decreases for decreasing $z$, for instance by $\simeq 20\%$ from $z=1$ to $z=0$). The analysis of the extended sample suggests that over the wider redshift range $0\lesssim z\lesssim 2.5$ the slope $\beta$ could be slightly steeper at higher $z$ ($\beta\simeq 0.2+0.2\log(1+z)$ describes the data as well as constant $\beta\simeq 0.23$) and that the redshift-dependence of the normalisation could be stronger ($\zeta\simeq 0.42$). The intrinsic scatter in $\sigma_0$ at given $M_*$ is $\simeq 0.07$ dex for the fiducial sample at $z\lesssim 1.2$, but is found to increase with redshift in the extended sample, being $\simeq 0.11$ dex at $z=2$. Our results suggest that, on average, the velocity dispersion of individual massive ($M_*\gtrsim 5\times 10^{11}M_\odot$) ETGs decreases with time while they evolve from $z\approx 1$ to $z\approx 0$.

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).

Context: The determination of the stellar initial mass function (IMF) of massive galaxies is one of the open problems in cosmology. Strong gravitational lensing is one of the few methods that allow us to constrain the IMF outside of the Local Group. Aims: The goal of this study is to statistically constrain the distribution in the IMF mismatch parameter, defined as the ratio between the true stellar mass of a galaxy and that inferred assuming a reference IMF, of massive galaxies from the BOSS CMASS sample. Methods: We take 23 strong lenses drawn from the CMASS sample, measure their Einstein radii and stellar masses using multi-band photometry from the Hyper Suprime-Cam survey, then fit a model distribution for the IMF mismatch parameter and dark matter halo mass to the whole sample. We use a prior on halo mass from weak lensing measurements and account for strong lensing selection effects in our model. Results: Assuming an NFW density profile for the dark matter distribution, we infer a value $\mu_{\mathrm{IMF}} = -0.04\pm0.11$ for the average base-10 logarithm of the IMF mismatch parameter, defined with respect to a Chabrier IMF. A Salpeter IMF is in tension with our measurements. Conclusions: Our results are consistent with a scenario in which the region of massive galaxies where the IMF normalization is significantly heavier than that of the Milky Way is much smaller than the scales $5\sim10$~kpc probed by the Einstein radius of the lenses in our sample, as recent spatially resolved studies of the IMF in massive galaxies suggest. The MCMC chains describing the posterior probability distribution of the model are available online, together with the code used to obtain them.

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 study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$\sigma$ significance. PSF also slightly flattens the measured colour profiles.

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

We present a detailed analysis of the inner mass structure of the Cosmic Horseshoe (J1148+1930) strong gravitational lens system observed with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3). In addition to the spectacular Einstein ring, this systems shows a radial arc. We obtained the redshift of the radial arc counter image $z_\text{s,r} = 1.961 \pm 0.001$ from Gemini observations. To disentangle the dark and luminous matter, we consider three different profiles for the dark matter distribution: a power-law profile, the NFW, and a generalized version of the NFW profile. For the luminous matter distribution, we base it on the observed light distribution that is fitted with three components: a point mass for the central light component resembling an active galactic nucleus, and the remaining two extended light components scaled by a constant M/L. To constrain the model further, we include published velocity dispersion measurements of the lens galaxy and perform a self-consistent lensing and axisymmetric Jeans dynamical modeling. Our model fits well to the observations including the radial arc, independent of the dark matter profile. Depending on the dark matter profile, we get a dark matter fraction between 60 % and 70 %. With our composite mass model we find that the radial arc helps to constrain the inner dark matter distribution of the Cosmic Hoseshoe independently of the dark matter profile.

Using the two large cosmological hydrodynamical simulations, Horizon-AGN (H-AGN) and Horizon-noAGN (H-noAGN, no AGN feedback), we investigate how a typical sub-grid model for AGN feedback affects the evolution of the total density profiles (dark matter + stars) at the effective radius of massive early-type galaxies (M*>10^11 Msun). We have studied the dependencies of the mass-weighted density slope gamma'_tot with the effective radius, the galaxy mass and the host halo mass at z~0.3 and found that the inclusion of AGN feedbackalways leads to a much better agreement with observational values and trends. Our analysis suggests also that the inclusion of AGN feedback favours a strong correlation between gamma'_tot and the density slope of the dark matter component while, in the absence of AGN activity, gamma'_tot is rather strongly correlated with the density slope of the stellar component. Finally, we find that gamma'_tot derived from our samples of galaxies increases from z=2 to z=0,in good agreement with the expected observational trend. The derived slopes are slightly lower than in the data when AGN is included because the simulated galaxies tend to be too extended, especially the least massive ones. However, the simulated compact galaxies without AGN feedback have gamma'_tot values that are significantly too high compared to observations.

Aims. We wish to determine the distribution of dark matter halo masses as a function of the stellar mass and the stellar mass profile, for massive galaxies in the BOSS CMASS sample. Methods. We use grizy photometry from HSC to obtain S\'ersic fits and stellar masses of CMASS galaxies for which HSC weak lensing data is available, visually selected to have spheroidal morphology. We apply a cut in stellar mass, $\log{M_*/M_\odot} > 11.0$,selecting $\sim$10, 000 objects. Using a Bayesian hierarchical inference method, we first investigate the distribution of S\'ersic index and size as a function of stellar mass. Then, making use of shear measurements from HSC, we measure the distribution of halo mass as a function of stellar mass, size and S\'ersic index. Results. Our data reveals a steep stellar mass-size relation $R_e \propto M_*^{\beta_R}$, with $\beta_R$ larger than unity, and a positive correlation between S\'ersic index and stellar mass: $n \propto M_*^{0.46}$. Halo mass scales approximately with the 1.7 power of the stellar mass. We do not find evidence for an additional dependence of halo mass on size or S\'ersic index at fixed stellar mass. Conclusions. Our results disfavour galaxy evolution models that predict significant differences in the size growth efficiency of galaxies living in low and high mass halos.

We investigate the local and line-of-sight overdensities of strong gravitational lens galaxies using wide-area multiband imaging from the Hyper Suprime-Cam Subaru Strategic Program. We present 41 new definite or probable lens candidates discovered in Data Release 2 of the survey. Using a combined sample of 87 galaxy-scale lenses out to a lens redshift of $z_{\mathrm{L}} \sim 0.8$, we compare galaxy number counts in lines of sight toward known and newly-discovered lenses in the survey to those of a control sample consisting of random lines of sight. We also compare the local overdensity of lens galaxies to a sample of "twin" galaxies that have a similar redshift and velocity dispersion to test whether lenses lie in different environments from similar non-lens galaxies. We find that lens fields contain higher number counts of galaxies compared to the control fields, but this effect arises from the local environment of the lens. Once galaxies in the lens plane are removed, the lens lines of sight are consistent with the control sample. The local environments of the lenses are overdense compared to the control sample, and are slightly overdense compared to those of the twin sample, although the significance is marginal. There is no significant evidence of the evolution of the local overdensity of lens environments with redshift.

The Initial Mass Function (IMF) for massive galaxies can be constrained by combining stellar dynamics with strong gravitational lensing. However, this method is limited by degeneracies between the density profile of dark matter and the stellar mass-to-light ratio. In this work we reduce this degeneracy by combining weak lensing together with strong lensing and stellar kinematics. Our analysis is based on two galaxy samples: 45 strong lenses from the SLACS survey and 1,700 massive quiescent galaxies from the SDSS main spectroscopic sample with weak lensing measurements from the Hyper Suprime-Cam survey. We use a Bayesian hierarchical approach to jointly model all three observables. We fit the data with models of varying complexity and show that a model with a radial gradient in the stellar mass-to-light ratio is required to simultaneously describe both galaxy samples. This result is driven by a subset of strong lenses with very steep total density profile, that cannot be fitted by models with no gradient. Our measurements are unable to determine whether $M_*/L$ gradients are due to variations in stellar population parameters at fixed IMF, or to gradients in the IMF itself. The inclusion of $M_*/L$ gradients decreases dramatically the inferred IMF normalisation, compared to previous lensing-based studies, with the exact value depending on the assumed dark matter profile. The main effect of strong lensing selection is to shift the stellar mass distribution towards the high mass end, while the halo mass and stellar IMF distribution at fixed stellar mass are not significantly affected.

We present a Bayesian hierarchical inference formalism to study the relation between the properties of dark matter halos and those of their central galaxies using weak gravitational lensing. Unlike traditional methods, this technique does not resort to stacking the weak lensing signal in bins, and thus allows for a more efficient use of the information content in the data. Our method is particularly useful for constraining scaling relations between two or more galaxy properties and dark matter halo mass, and can also be used to constrain the intrinsic scatter in these scaling relations. We show that, if observational scatter is not properly accounted for, the traditional stacking method can produce biased results when exploring correlations between multiple galaxy properties and halo mass. For example, this bias can affect studies of the joint correlation between galaxy mass, halo mass, and galaxy size, or galaxy color. In contrast, our method easily and efficiently handles the intrinsic and observational scatter in multiple galaxy properties and halo mass. We test our method on mocks with varying degrees of complexity. We find that we can recover the mean halo mass and concentration, each with a $0.1$ dex accuracy, and the intrinsic scatter in halo mass with a $0.05$ dex accuracy. In its current version, our method will be most useful for studying the weak lensing signal around central galaxies in groups and clusters, as well as massive galaxies samples with $\log{M_*} > 11$, which have low satellite fractions.

We present a catalogue of 38 spectroscopically detected strong galaxy-galaxy gravitational lens candidates identified in the Sloan Digital Sky Survey IV (SDSS-IV). We were able to simulate narrow-band images for 8 of them demonstrating evidence of multiple images. Two of our systems are compound lens candidates, each with 2 background source-planes. One of these compound systems shows clear lensing features in the narrow-band image. Our sample is based on 2812 galaxies observed by the Mapping Nearby Galaxies at APO (MaNGA) integral field unit (IFU). This Spectroscopic Identification of Lensing Objects (SILO) survey extends the methodology of the Sloan Lens ACS Survey (SLACS) and BOSS Emission-Line Survey (BELLS) to lower redshift and multiple IFU spectra. We searched ~ 1.5 million spectra, of which 3065 contained multiple high signal-to-noise background emission-lines or a resolved [OII] doublet, that are included in this catalogue. Upon manual inspection, we discovered regions with multiple spectra containing background emission-lines at the same redshift, providing evidence of a common source-plane geometry which was not possible in previous SLACS and BELLS discovery programs. We estimate more than half of our candidates have an Einstein radius > 1.7", which is significantly greater than seen in SLACS and BELLS. These larger Einstein radii produce more extended images of the background galaxy increasing the probability that a background emission-line will enter one of the IFU spectroscopic fibres, making detection more likely.

Time delay lensing is a mature and competitive cosmological probe. However, it is limited in accuracy by the well-known problem of the mass-sheet degeneracy: too rigid assumptions on the density profile of the lens can potentially bias the inference on cosmological parameters. I investigate the degeneracy between the choice of the lens density profile and the inference on the Hubble constant, focusing on double image systems. By expanding lensing observables in terms of the local derivatives of the lens potential around the Einstein radius, and assuming circular symmetry, I show that three degrees of freedom in the radial direction are necessary to achieve a few percent accuracy in the time-delay distance. Additionally, while the time delay is strongly dependent on the second derivative of the potential, observables typically used to constrain lens models in time-delay studies, such as image position and radial magnification information, are mostly sensitive to the first and third derivatives, making it very challenging to accurately determine time-delay distances with lensing data alone. Tests on mock observations show that the assumption of a power-law density profile results in a 5% average bias on $H_0$, with a 6% scatter. Using a more flexible model and adding unbiased velocity dispersion constraints allows to obtain an inference with 1% accuracy. A power-law model can still provide 3% accuracy if velocity dispersion measurements are used to constrain its slope. Although this work is based on the assumption of axisymmetry, its main findings can be generalized to cases with moderate ellipticity.

The Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) is an excellent survey for the search for strong lenses, thanks to its area, image quality and depth. We use three different methods to look for lenses among 43,000 luminous red galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS) sample with photometry from the S16A internal data release of the HSC SSP. The first method is a newly developed algorithm, named YATTALENS, which looks for arc-like features around massive galaxies and then estimates the likelihood of an object being a lens by performing a lens model fit. The second method, CHITAH, is a modeling-based algorithm originally developed to look for lensed quasars. The third method makes use of spectroscopic data to look for emission lines from objects at a different redshift from that of the main galaxy. We find 15 definite lenses, 36 highly probable lenses and 282 possible lenses. Among the three methods, YATTALENS, which was developed specifically for this problem, performs best in terms of both completeness and purity. Nevertheless five highly probable lenses were missed by YATTALENS but found by the other two methods, indicating that the three methods are highly complementary. Based on these numbers we expect to find $\sim$300 definite or probable lenses by the end of the HSC SSP.

We carry out a systematic investigation of the total mass density profile of massive (Mstar~3e11 Msun) early-type galaxies and its dependence on redshift, specifically in the range 0<z<1. We start from a large sample of SDSS early-type galaxies with stellar masses and effective radii measured assuming two different profiles, de Vaucouleurs and S\'{e}rsic. We assign dark matter haloes to galaxies via abundance matching relations with standard LCDM profiles and concentrations. We then compute the total, mass-weighted density slope at the effective radius gamma', and study its redshift dependence at fixed stellar mass. We find that a necessary condition to induce an increasingly flatter gamma' at higher redshifts, as suggested by current strong lensing data, is to allow the intrinsic stellar profile of massive galaxies to be S\'{e}rsic and the input S\'{e}rsic index n to vary with redshift approximately as n(z)~(1+z)^(-1). This conclusion holds irrespective of the input Mstar-Mhalo relation, the assumed stellar initial mass function, or even the chosen level of adiabatic contraction in the model. Secondary contributors to the observed redshift evolution of gamma' may come from an increased contribution at higher redshifts of adiabatic contraction and/or bottom-light stellar initial mass functions. The strong lensing selection effects we have simulated seem not to contribute to this effect. A steadily increasing S\'{e}rsic index with cosmic time is supported by independent observations, though it is not yet clear whether cosmological hierarchical models (e.g., mergers) are capable of reproducing such a fast and sharp evolution.

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.

We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and S\'{e}rsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a S\'{e}rsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations.

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.

The stellar initial mass function (IMF) of early-type galaxies is the combination of the IMF of the stellar population formed in-situ and that of accreted stellar populations. Using as an observable the effective IMF $\alpha_{IMF}$, defined as the ratio between the true stellar mass of a galaxy and the stellar mass inferred assuming a Salpeter IMF, we present a theoretical model for its evolution as a result of dry mergers. We use a simple dry merger evolution model, based on cosmological $N$-body simulations, together with empirically motivated prescriptions for the IMF to make predictions for how the effective IMF of massive early-type galaxies changes from $z=2$ to $z=0$. We find that the IMF normalization of individual galaxies becomes lighter with time. At fixed velocity dispersion, $\alpha_{IMF}$ is predicted to be constant with redshift. Current constraints on the evolution of the IMF are in slight tension with this prediction, even though systematic uncertainties prevent a conclusive statement. The correlation of $\alpha_{IMF}$ with stellar mass becomes shallower with time, while the correlation between $\alpha_{IMF}$ and velocity dispersion is mostly preserved by dry mergers. We also find that dry mergers can mix the dependence of the IMF on stellar mass and velocity dispersion, making it challenging to infer, from $z=0$ observations of global galactic properties, what is the quantity that is originally coupled with the IMF.

We report the serendipitous discovery of HSC J142449-005322, a double source plane lens system in the Hyper Suprime-Cam Subaru Strategic Program. We dub the system Eye of Horus. The lens galaxy is a very massive early-type galaxy with stellar mass of ~7x10^11 Msun located at z_L=0.795. The system exhibits two arcs/rings with clearly different colors, including several knots. We have performed spectroscopic follow-up observations of the system with FIRE on Magellan. The outer ring is confirmed at z_S2=1.988 with multiple emission lines, while the inner arc and counterimage is confirmed at z_S1=1.302. This makes it the first double source plane system with spectroscopic redshifts of both sources. Interestingly, redshifts of two of the knots embedded in the outer ring are found to be offset by delta_z=0.002 from the other knots, suggesting that the outer ring consists of at least two distinct components in the source plane. We perform lens modeling with two independent codes and successfully reproduce the main features of the system. However, two of the lensed sources separated by ~0.7 arcsec cannot be reproduced by a smooth potential, and the addition of substructure to the lens potential is required to reproduce them. Higher-resolution imaging of the system will help decipher the origin of this lensing feature and potentially detect the substructure.

We present a detailed analysis of an individual case of gravitational lensing of a $z\sim8$ Lyman-Break galaxy (LBG) in a blank field, identified in Hubble Space Telescope imaging obtained as part of the Brightest of Reionizing Galaxies survey. To investigate the close proximity of the bright ($m_{AB}=25.8$) $Y_{098}$-dropout to a small group of foreground galaxies, we obtained deep spectroscopy of the dropout and two foreground galaxies using VLT/X-Shooter. We detect H-$\alpha$, H-$\beta$, [OIII] and [OII] emission in the brightest two foreground galaxies (unresolved at the natural seeing of $0.8$ arcsec), placing the pair at $z=1.327$. We can rule out emission lines contributing all of the observed broadband flux in $H_{160}$ band at $70\sigma$, allowing us to exclude the $z\sim8$ candidate as a low redshift interloper with broadband photometry dominated by strong emission lines. The foreground galaxy pair lies at the peak of the luminosity, redshift and separation distributions for deflectors of strongly lensed $z\sim8$ objects, and we make a marginal detection of a demagnified secondary image in the deepest ($J_{125}$) filter. We show that the configuration can be accurately modelled by a singular isothermal ellipsoidal deflector and a S\'{e}rsic source magnified by a factor of $\mu=4.3\pm0.2$. The reconstructed source in the best-fitting model is consistent with luminosities and morphologies of $z\sim8$ LBGs in the literature. The lens model yields a group mass of $9.62\pm0.31\times10^{11} M_{\odot}$ and a stellar mass-to-light ratio for the brightest deflector galaxy of $M_{\star}/L_{B}=2.3^{+0.8}_{-0.6} M_{\odot}/L_{\odot}$ within its effective radius. The foreground galaxies' redshifts would make this one of the few strong lensing deflectors discovered at $z>1$.

We present spectroscopy and laser guide star adaptive optics (LGSAO) images of the doubly imaged lensed quasar SDSS J1206+4332. We revise the deflector redshift proposed previously to $z_{d}=0.745,$ and measure for the first time its velocity dispersion $\sigma=(290\pm30)$ km/s. The LGSAO data show the lensed quasar host galaxy stretching over the astroid caustic thus forming an extra pair of merging images, which was previously thought to be an unrelated galaxy in seeing limited data. Owing to the peculiar geometry, the lens acts as a natural coronagraph on the broad-line region of the quasar so that only narrow [O III] emission is found in the fold arc. We use the data to reconstruct the source structure and deflector potential, including nearby perturbers. We reconstruct the point-spread function (PSF) from the quasar images themselves, since no additional point source is present in the field of view. From gravitational lensing and stellar dynamics, we find the slope of the total mass density profile to be $\gamma^{\prime}=-\log\rho/\log r =1.93\pm0.09.$ We discuss the potential of SDSS J1206+4332 for measuring time delay distance (and thus H$_0$ and other cosmological parameters), or as a standard ruler, in combination with the time delay published by the COSMOGRAIL collaboration. We conclude that this system is very promising for cosmography. However, in order to achieve competitive precision and accuracy, an independent characterization of the PSF is needed. Spatially resolved kinematics of the deflector would reduce the uncertainties further. Both are within the reach of current observational facilities.

In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging, we have found four images of a single supernova forming an Einstein cross configuration around a redshift z=0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The cluster's gravitational potential also creates multiple images of the z=1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses.

We investigate the cosmic evolution of the internal structure of massive early-type galaxies over half of the age of the Universe. We perform a joint lensing and stellar dynamics analysis of a sample of 81 strong lenses from the SL2S and SLACS surveys and combine the results with a hierarchical Bayesian inference method to measure the distribution of dark matter mass and stellar IMF across the population of massive early-type galaxies. Lensing selection effects are taken into account. We find that the dark matter mass projected within the inner 5 kpc increases for increasing redshift, decreases for increasing stellar mass density, but is roughly constant along the evolutionary tracks of early-type galaxies. The average dark matter slope is consistent with that of an NFW profile, but is not well constrained. The stellar IMF normalization is close to a Salpeter IMF at $\log{M_*} = 11.5$ and scales strongly with increasing stellar mass. No dependence of the IMF on redshift or stellar mass density is detected. The anti-correlation between dark matter mass and stellar mass density supports the idea of mergers being more frequent in more massive dark matter halos.

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

Several studies have suggested that the observed size evolution of massive early-type galaxies (ETGs) can be explained as a combination of dry mergers and progenitor bias, at least since z\sim1. In this paper we carry out a new test of the dry-merger scenario based on recent lensing measurements of the evolution of the mass density profile of ETGs. We construct a theoretical model for the joint evolution of the size and mass density profile slope \gamma' driven by dry mergers occurring at rates given by cosmological simulations. Such dry-merger model predicts a strong decrease of \gamma' with cosmic time, inconsistent with the almost constant \gamma' inferred from observations in the redshift range 0<z<1. We then show with a simple toy model that a modest amount of cold gas in the mergers -- consistent with the upper limits on recent star formation in ETGs -- is sufficient to reconcile the model with measurements of \gamma'. By fitting for the amount of gas accreted during mergers, we find that models with dissipation are consistent with observations of the evolution in both size and density slope, if \sim4% of the total final stellar mass arises from the gas accreted since z\sim1. Purely dry merger models are ruled out at > 99% CL. We thus suggest a scenario where the outer regions of massive ETGs grow by accretion of stars and dark matter, while small amounts of dissipation and nuclear star formation conspire to keep the mass density profile constant and approximately isothermal.

We present RingFinder, a tool for finding galaxy-scale strong gravitational lenses in multiband imaging data. By construction, the method is sensitive to configurations involving a massive foreground early-type galaxy and a faint, background, blue source. RingFinder detects the presence of blue residuals embedded in an otherwise smooth red light distribution by difference imaging in two bands. The method is automated for efficient application to current and future surveys, having originally been designed for the 150-deg2 Canada France Hawaii Telescope Legacy Survey (CFHTLS). We describe each of the steps of RingFinder. We then carry out extensive simulations to assess completeness and purity. For sources with magnification mu>4, RingFinder reaches 42% (resp. 25%) completeness and 29% (resp. 86%) purity before (resp. after) visual inspection. The completeness of RingFinder is substantially improved in the particular range of Einstein radii 0.8 < REin < 2. and lensed images brighter than g = 22.5, where it can be as high as 70%. RingFinder does not introduce any significant bias in the source or deflector population. We conclude by presenting the final catalog of RingFinder CFHTLS galaxy-scale strong lens candidates. Additional information obtained with Hubble Space Telescope and Keck Adaptive Optics high resolution imaging, and with Keck and Very Large Telescope spectroscopy, is used to assess the validity of our classification, and measure the redshift of the foreground and the background objects. From an initial sample of 640,000 early type galaxies, RingFinder returns 2500 candidates, which we further reduce by visual inspection to 330 candidates. We confirm 33 new gravitational lenses from the main sample of candidates, plus an additional 16 systems taken from earlier versions of RingFinder. First applications are presented in the SL2S galaxy-scale Lens Sample paper series.

We present optical and near infrared spectroscopy obtained at Keck, VLT, and Gemini for a sample of 36 secure strong gravitational lens systems and 17 candidates identified as part of the SL2S survey. The deflectors are massive early-type galaxies in the redshift range z_d=0.2-0.8, while the lensed sources are at z_s=1-3.5. We combine this data with photometric and lensing measurements presented in the companion paper III and with lenses from the SLACS and LSD surveys to investigate the cosmic evolution of the internal structure of massive early-type galaxies over half the age of the universe. We study the dependence of the slope of the total mass density profile \gamma' (\rho(r)\propto r^{-\gamma'}) on stellar mass, size, and redshift. We find that two parameters are sufficent to determine \gamma' with less than 6% residual scatter. At fixed redshift, \gamma' depends solely on the surface stellar mass density \partial \gamma'/ \partial \Sigma_*=0.38\pm 0.07, i.e. galaxies with denser stars also have steeper slopes. At fixed M_* and R_{eff}, \gamma' depends on redshift, in the sense that galaxies at a lower redshift have steeper slopes (\partial \gamma' / \partial z = -0.31\pm 0.10). However, the mean redshift evolution of \gamma' for an individual galaxy is consistent with zero d\gamma'/dz=-0.10\pm0.12. This result is obtained by combining our measured dependencies of \gamma' on z,M_*,R_{eff} with the evolution of the R_{eff}-M_* taken from the literature, and is broadly consistent with current models of the formation and evolution of massive early-type galaxies. Detailed quantitative comparisons of our results with theory will provide qualitatively new information on the detailed physical processes at work.

We present Hubble Space Telescope (HST) imaging data and CFHT Near IR ground-based images for the final sample of 56 candidate galaxy-scale lenses uncovered in the CFHT Legacy Survey as part of the Strong Lensing in the Legacy Survey (SL2S) project. The new images are used to perform lens modeling, measure surface photometry, and estimate stellar masses of the deflector early-type galaxies. Lens modeling is performed on the HST images (or CFHT when HST is not available) by fitting the spatially extended light distribution of the lensed features assuming a singular isothermal ellipsoid mass profile and by reconstructing the intrinsic source light distribution on a pixelized grid. Based on the analysis of systematic uncertainties and comparison with inference based on different methods we estimate that our Einstein Radii are accurate to \sim3%. HST imaging provides a much higher success rate in confirming gravitational lenses and measuring their Einstein radii than CFHT imaging does. Lens modeling with ground-based images however, when successful, yields Einstein radius measurements that are competitive with spaced-based images. Information from the lens models is used together with spectroscopic information from the companion paper IV to classify the systems, resulting in a final sample of 39 confirmed (grade-A) lenses and 17 promising candidates. The redshifts of the main deflector span a range 0.3<zd< 0.8, providing an excellent sample for the study of the cosmic evolution of the mass distribution of early-type galaxies over the second half of the history of the Universe.

We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6(+1.2,-0.9) light days, a width of the BLR of 6.9(+1.2,-1.1) light days, and a disk opening angle of 25\pm10 degrees above the plane. We also constrain the inclination angle to be 9(+7,-5) degrees, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log10(M(BH)/Msun) = 7.57(+0.44,-0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log10(f) = 0.78(+0.44,-0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M(BH)-{\sigma}* relation for AGN and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios.

Stars and dark matter account for most of the mass of early-type galaxies, but uncertainties in the stellar population and the dark matter profile make it challenging to distinguish between the two components. Nevertheless, precise observations of stellar and dark matter are extremely valuable for testing the many models of structure formation and evolution. We present a measurement of the stellar mass and inner slope of the dark matter halo of a massive early-type galaxy at $z=0.222$. The galaxy is the foreground deflector of the double Einstein ring gravitational lens system SDSSJ0946+1006, also known as the "Jackpot". By combining the tools of lensing and dynamics, we first constrain the mean slope of the total mass density profile ($\rho_{\rm{tot}}\propto r^{-\gamma'}$) within the radius of the outer ring to be $\gamma' = 1.98\pm0.02\pm0.01$. Then we obtain a bulge-halo decomposition, assuming a power-law form for the dark matter halo. Our analysis yields $\gamma_{\rm{DM}} = 1.7\pm0.2$ for the inner slope of the dark matter profile, in agreement with theoretical findings on the distribution of dark matter in ellipticals, and a stellar mass from lensing and dynamics $M_*^{\rm{LD}} = 5.5_{-1.3}^{+0.4}\times10^{11}M_\Sun$. By comparing this measurement with stellar masses inferred from stellar population synthesis fitting we find that a Salpeter IMF provides a good description of the stellar population of the lens while the probability of the IMF being heavier than Chabrier is 95%. Our data suggest that growth by accretion of small systems from a compact red nugget is a plausible formation scenario for this object.

Some active galactic nuclei (AGN) are surrounded by extended emission-line regions (EELRs), which trace both the illumination pattern of escaping radiation and its history over the light-travel time from the AGN to the gas. From a new set of such EELRs, we present evidence that the AGN in many Seyfert galaxies undergo luminous episodes 20,000-200,000 years in duration. Motivated by the discovery of the spectacular nebula known as Hanny's Voorwerp, ionized by a powerful AGN which has apparently faded dramatically within ~ 100,000 years, Galaxy Zoo volunteers have carried out both targeted and serendipitous searches for similar emission-line clouds around low-redshift galaxies.We present the resulting list of candidates and describe spectroscopy identifying 19 galaxies with AGN-ionized regions at projected radii > 10 kpc. This search recovered known EELRs and identified additional previously unknown cases, one with detected emission to r = 37 kpc. At least 14/19 are in interacting or merging systems; tidal tails are a prime source of extraplanar ionized gas. We see a mix of one- and two-sided structures, with observed cone angles from 23-112 degrees. We consider the energy balance in the ionized clouds, with lower and upper bounds on ionizing luminosity from recombination and ionization-parameter arguments, and estimate the luminosity of the core from the far-infrared data. The implied ratio of ionizing radiation seen by the clouds to that emitted by the nucleus, for a constant nuclear source, ranges from 0.02 to > 12; 7/19 exceed unity. Small values imply heavily obscured AGN. However, large values may require that the AGN has faded over tens of thousands of years, giving us several examples of systems in which such dramatic long-period variation has occurred; this is the only current technique for addressing these timescales in AGN history. (Abridged)

The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in Spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our dataset covers a total duration of 4.5 months. During this time, Mrk 50 was highly variable, exhibiting a maximum variability amplitude of a factor of 4 in the U-band continuum and a factor of 2 in the H-beta line. Using standard cross-correlation techniques, we find that H-beta and H-gamma lag the V-band continuum by tau_cen = 10.64(-0.93,+0.82) and 8.43(-1.28,+1.30) days, respectively, while the lag of He II 4686 is unresolved. The H-beta line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the high-velocity wings than in the line core, consistent with an origin in a broad-line region dominated by orbital motion rather than infall or outflow. Assuming a virial normalization factor of f=5.25, the virial estimate of the black hole mass is (3.2+-0.5)*10^7 solar masses. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics.

Weak lensing is one of the best available diagnostic tools to measure the total density profiles of distant clusters of galaxies. Unfortunately, it suffers from the well-known mass-sheet degeneracy, so that weak lensing analyses cannot lead to fully reliable determinations of the total mass of the clusters. One possible way to set the relevant scale of the density profile would be to make a direct measurement of the magnification produced by the clusters as gravitational lenses. In this paper we revisit a suggestion made a few years ago for this general purpose, based on the use of the Fundamental Plane as a standard rod for early-type galaxies. Here we move one step further, beyond the simple outline of the idea given earlier, and quantify some statistical properties of this innovative diagnostic tool, with the final goal of identifying clear guidelines for a future observational test on concrete cases, which turns out to be well within the current instrument capabilities. The study is carried out by discussing the statistical properties of Fundamental Plane measurements for a sample of early-type source galaxies behind a massive cluster, for which a weak lensing analysis is assumed to be available. We proceed to study the best strategy to use Fundamental Plane measurements to determine the mass scale of a given cluster and find that the optimal choice is that of a sample of early-type galaxies behind the cluster distributed approximately uniformly in the sky. We also show that substructures do not contribute much to the magnification signal that we are looking for, but only add a modest amount of scatter. We find that for a massive cluster (M > 10^15 M_Sun) located at redshift 0.3 +/- 0.1, a set of about 20 Fundamental Plane measurements, combined with a good weak lensing analysis, should be able to lead to a mass determination with a precision of 20% or better.

The Seyfert 1 galaxy Zw 229-015 is among the brightest active galaxies being monitored by the Kepler mission. In order to determine the black hole mass in Zw 229-015 from H-beta reverberation mapping, we have carried out nightly observations with the Kast Spectrograph at the Lick 3m telescope during the dark runs from June through December 2010, obtaining 54 spectroscopic observations in total. We have also obtained nightly V-band imaging with the Katzman Automatic Imaging Telescope at Lick Observatory and with the 0.9m telescope at the Brigham Young University West Mountain Observatory over the same period. We detect strong variability in the source, which exhibited more than a factor of 2 change in broad H-beta flux. From cross-correlation measurements, we find that the H-beta light curve has a rest-frame lag of 3.86(+0.69,-0.90) days with respect to the V-band continuum variations. We also measure reverberation lags for H-alpha and H-gamma and find an upper limit to the H-delta lag. Combining the H-beta lag measurement with a broad H-beta width of sigma = 1590+/-47 km/s measured from the root-mean-square variability spectrum, we obtain a virial estimate of M_BH = 1.00(-0.24,+0.19)*10^7 solar masses for the black hole in Zw 229-015. As a Kepler target, Zw 229-015 will eventually have one of the highest-quality optical light curves ever measured for any active galaxy, and the black hole mass determined from reverberation mapping will serve as a benchmark for testing relationships between black hole mass and continuum variability characteristics in active galactic nuclei.