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.

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 a study of the stellar host galaxy, CO (1$-$0) molecular gas distribution and AGN emission on 50 to 500 pc-scales of the gravitationally lensed dust-obscured AGN MG J0751+2716 and JVAS B1938+666 at redshifts 3.200 and 2.059, respectively. By correcting for the lensing distortion using a grid-based lens modelling technique, we spatially locate the different emitting regions in the source plane for the first time. Both AGN host galaxies have 300 to 500 pc-scale size and surface brightness consistent with a bulge/pseudo-bulge, and 2 kpc-scale AGN radio jets that are embedded in extended molecular gas reservoirs that are 5 to 20 kpc in size. The CO (1$-$0) velocity fields show structures possibly associated with discs (elongated velocity gradients) and interacting objects (off-axis velocity components). There is evidence for a decrement in the CO (1$-$0) surface brightness at the location of the host galaxy, which may indicate radiative feedback from the AGN, or offset star formation.We find CO-H$_2$ conversion factors of around $\alpha_{\rm CO} = 1.5\pm0.5$ (K km s$^{-1}$ pc$^2$)$^{-1}$, molecular gas masses of $> 3\times10^{10}$ M$_{\odot}$, dynamical masses of $\sim 10^{11}$ M$_{\odot}$ and gas fractions of around 60 per cent. The intrinsic CO line luminosities are comparable to those of unobscured AGN and dusty star-forming galaxies at similar redshifts, but the infrared luminosities are lower, suggesting that the targets are less efficient at forming stars. Therefore, they may belong to the AGN feedback phase predicted by galaxy formation models, because they are not efficiently forming stars considering their large amount of molecular gas.

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.

Minor mergers have been proposed as the driving mechanism for the size growth of quiescent galaxies with decreasing redshift. The process whereby large star-forming galaxies quench and join the quiescent population at the large size end has also been suggested as an explanation for this size growth. Given the clear association of quenching with clusters, we explore this mechanism by studying the structural properties of 23 spectroscopically identified recently quenched (or "poststarburst" (PSB)) cluster galaxies at $z\sim1$. Despite clear PSB spectral signatures implying rapid and violent quenching, 87\% of these galaxies have symmetric, undisturbed morphologies in the stellar continuum. Remarkably, they follow a mass-size relation lying midway between the star-forming and quiescent field relations, with sizes $0.1$ dex smaller than $z\sim1$ star-forming galaxies at log$(M_{*}/M_{\odot})=10.5$. This implies a rapid change in the light profile without directly effecting the stellar distribution, suggesting changes in the mass-to-light ratio gradients across the galaxy are responsible. We develop fading toy models to explore how star-forming galaxies move across the mass-size plane as their stellar populations fade to match those of the PSBs. "Outside-in" fading has the potential to reproduce the contraction in size and increase in bulge-dominance observed between star-forming and PSB cluster galaxies. Since cluster PSBs lie on the large size end of the quiescent mass-size relation, and our previous work shows cluster galaxies are smaller than field galaxies, the sizes of quiescent galaxies must grow both from the quenching of star-forming galaxies and dry minor mergers.

We present near-infrared observations of 42 gravitationally lensed galaxies obtained in the framework of the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey, a program aimed at investigating the spatially resolved properties of the ionised gas in 1.2<z<2.5 galaxies by means of a full coverage of the YJ, H and K near-infrared bands. Detailed metallicity maps and gradients are derived for a sub-sample of 28 galaxies from reconstructed source plane emission line maps, exploiting the variety of different emission line diagnostics provided by the broad wavelength coverage of the survey. About 85% of these galaxies are characterised by metallicity gradients shallower than 0.05 dex/kpc and 89% are consistent with flat slope within 3$\sigma$ (67% within 1$\sigma$), suggesting a mild evolution with cosmic time. In the context of cosmological simulations and chemical evolution models, the presence of efficient feedback mechanisms and/or extended star formation profiles on top of the classical "inside-out" scenario of mass assembly is generally required to reproduce the observed flatness of the metallicity gradients beyond z$\sim$1 . Three galaxies with significantly (> 3$\sigma$) "inverted" gradients are also found, showing an anti-correlation between metallicity and star formation rate density on local scales, possibly suggesting recent episodes of pristine gas accretion or strong radial flows in place. Nevertheless, the individual metallicity maps are characterised by a variety of different morphologies, with flat radial gradients sometimes hiding non-axisymmetric variations on kpc scales which are washed out by azimuthal averages, especially in interacting systems or in those undergoing local episodes of recent star formation.

We report the results of the STRong lensing Insights from the Dark Energy Survey (STRIDES) follow-up campaign of the late 2017/early 2018 season. We obtained spectra of 65 lensed quasar candidates either with EFOSC2 on the NTT or ESI on Keck, which confirm 10 new gravitationally lensed quasars and 10 quasar pairs with similar spectra, but which do not show a lensing galaxy in DES images. Eight lensed quasars are doubly imaged with source redshifts between 0.99 and 2.90, one is triply imaged by a group (DESJ0345-2545, $z=1.68$), and one is quadruply imaged (quad: DESJ0053-2012, $z=3.8$). Singular isothermal ellipsoid models for the doubles, based on high-resolution imaging from SAMI on SOAR or NIRC2 on Keck, give total magnifications between 3.2 and 5.6, and Einstein radii between 0.49 and 1.97 arcseconds. After spectroscopic follow-up, we extract multi-epoch $grizY$ photometry of confirmed lensed quasars and contaminant quasar+star pairs from the first 4 years of DES data using parametric multi-band modelling, and compare variability in each system's components. By measuring the reduced ${\chi}^2$ associated with fitting all epochs to the same magnitude, we find a simple cut on the less variable component that retains all confirmed lensed quasars, while removing 94 per cent of contaminant systems with stellar components. Based on our spectroscopic follow-up, this variability information can improve selection of lensed quasars and quasar pairs from 34-45 per cent to 51-70 per cent, with the majority of remaining contaminants being compact star-forming galaxies. Using mock lensed quasar lightcurves we demonstrate that selection based only on variability will over-represent the quad fraction by 10 per cent over a complete DES magnitude-limited sample (excluding microlensing differences), explained by the magnification bias and hence lower luminosity (more variable) sources in quads.

We present an analysis of seven strongly gravitationally lensed quasars and the corresponding constraints on the properties of dark matter. Our results are derived by modelling the lensed image positions and flux-ratios using a combination of smooth macro models and a population of low-mass haloes within the mass range 10^6 to 10^9 Msun. Our lens models explicitly include higher-order complexity in the form of stellar discs and luminous satellites, as well as low-mass haloes located along the observed lines of sight for the first time. Assuming a Cold Dark Matter (CDM) cosmology, we infer an average total mass fraction in substructure of f_sub = 0.012^{+0.007}_{-0.004} (68 per cent confidence limits), which is in agreement with the predictions from CDM hydrodynamical simulations to within 1 sigma. This result is closer to the predictions than those from previous studies that did not include line-of-sight haloes. Under the assumption of a thermal relic dark matter model, we derive a lower limit on the particle relic mass of m th > 5.58 keV (95 per cent confidence limits), which is consistent with a value of m_th > 5.3 keV from the recent analysis of the Ly-alpha forest. We also identify two main sources of possible systematic errors and conclude that deeper investigations in the complex structure of lens galaxies as well as the size of the background sources should be a priority for this field.

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

We present two catalogues of active galactic nucleus (AGN) candidates selected from the latest data of two all-sky surveys -- Data Release 2 (DR2) of the \emph{Gaia} mission and the unWISE catalogue of the \emph{Wide-field Infrared Survey Explorer} (\emph{WISE}). We train a random forest classifier to predict the probability of each source in the \emph{Gaia}-unWISE joint sample being an AGN, $P_{\rm RF}$, based on \emph{Gaia} astrometric and photometric measurements and unWISE photometry. The two catalogues, which we designate C75 and R85, are constructed by applying different $P_{\rm RF}$ threshold cuts to achieve an overall completeness of 75\% ($\approx$90\% at \emph{Gaia} $G\leq20$ mag) and reliability of 85\% respectively. The C75 (R85) catalogue contains 2,734,464 (2,182,193) AGN candidates across the effective 36,000 deg$^2$ sky, of which $\approx$0.91 (0.52) million are new discoveries. Photometric redshifts of the AGN candidates are derived by a random forest regressor using \emph{Gaia} and \emph{WISE} magnitudes and colours. The estimated overall photometric redshift accuracy is 0.11. Cross-matching the AGN candidates with a sample of known bright cluster galaxies, we identify a high-probability strongly-lensed AGN candidate system, SDSS\,J1326$+$4806, with a large image separation of 21\farcs06. All the AGN candidates in our catalogues will have $\sim$5-year long light curves from \emph{Gaia} by the end of the mission, and thus will be a great resource for AGN variability studies. Our AGN catalogues will also be helpful in AGN target selections for future spectroscopic surveys, especially ones in the southern hemisphere. The C75 catalogue can be downloaded at https://www.ast.cam.ac.uk/~ypshu/AGN_Catalogues.html.

Low mass, low metallicity galaxies at low to moderate ($z\lesssim3$) redshifts offer the best opportunity for detailed examination of the interplay between massive stars, ionizing radiation and gas in sources similar to those that likely reionized the universe. We present new narrowband Hubble Space Telescope observations of Ly$\alpha$ emission and the adjacent ultraviolet (UV) continuum in the low mass ($M_{\star} = 2 \times 10^8$ M$_{\odot}$), low metallicity ($Z\sim1/20$ Z$_{\odot}$) and highly ionized gravitationally lensed galaxy SL2S J02176$-$0513 at $z=1.844$. The galaxy has strong Ly$\alpha$ emission with photometric equivalent width $W^{\rm phot}_{\rm{Ly}\alpha} = 218 \pm 12$ \AA, at odds with the Ly$\alpha$ escape fraction of 10%. However, the spectroscopic Ly$\alpha$ profile suggests the presence of broad absorption underlying the emission, and the total equivalent width is consistent with the escape fraction once this underlying absorption is included. The Ly$\alpha$ emission is more spatially extended than the UV continuum, and the 0.14" spatial resolution of HST coupled with the magnification of gravitational lensing enables us to examine the distribution of Ly$\alpha$ and the UV continuum on sub-kiloparsec scales. We find that the peaks of the Ly$\alpha$ emission and the UV continuum are offset by 650 pc, and there is no Ly$\alpha$ emission arising from the region with the strongest UV light. Our combined spectroscopic and imaging data imply a significant range in neutral hydrogen column density across the object. These observations offer indirect support for a model in which ionizing radiation escapes from galaxies through channels with low column density of neutral gas.

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 introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

Strong-gravitational lens systems with quadruply-imaged quasars (quads) are unique probes to address several fundamental problems in cosmology and astrophysics. Although they are intrinsically very rare, ongoing and planned wide-field deep-sky surveys are set to discover thousands of such systems in the next decade. It is thus paramount to devise a general framework to model strong-lens systems to cope with this large influx without being limited by expert investigator time. We propose such a general modelling framework (implemented with the publicly available software Lenstronomy) and apply it to uniformly model three-band Hubble Space Telescope Wide Field Camera 3 images of 13 quads. This is the largest uniformly modelled sample of quads to date and paves the way for a variety of studies. To illustrate the scientific content of the sample, we investigate the alignment between the mass and light distribution in the deflectors. The position angles of these distributions are well-aligned, except when there is strong external shear. However, we find no correlation between the ellipticity of the light and mass distributions. We also show that the observed flux-ratios between the images depart significantly from the predictions of simple smooth models. The departures are strongest in the bluest band, consistent with microlensing being the dominant cause in addition to millilensing. Future papers will exploit this rich dataset in combination with ground based spectroscopy and time delays to determine quantities such as the Hubble constant, the free streaming length of dark matter, and the normalization of the initial stellar mass function.

We report the discovery and spectroscopic confirmation of 22 new gravitationally lensed quasars found using $\textit{Gaia}$ data release 2. The selection was made using several techniques: multiple $\textit{Gaia}$ detections around objects in quasar candidate catalogues, modelling of unWISE coadd pixels using $\textit{Gaia}$ astrometry, and $\textit{Gaia}$ detections offset from photometric and spectroscopic galaxies. Spectra of 33 candidates were obtained with the William Herschel Telescope, 22 of which are lensed quasars, 2 highly probably lensed quasars, 5 nearly identical quasar pairs, 1 inconclusive system, and 3 contaminants. Of the 3 confirmed quadruply imaged systems, J2145+6345 is a 2.1 arcsecond separation quad with 4 bright images ($G$=16.86, 17.26, 18.34, 18.56), making it ideal for time delay monitoring. Analysing this new sample alongside known lenses in the Pan-STARRS footprint, and comparing to expected numbers of lenses, we show that, as expected, we are biased towards systems with bright lensing galaxies and low source redshifts. We discuss possible techniques to remove this bias from future searches. A |b|>20 complete sample of lensed quasars detected by $\textit{Gaia}$ and with image separations above 1 arcsecond will provide a valuable statistical sample of around 350 systems. Currently only 96 known lenses satisfy these criteria, yet promisingly, our unWISE modelling technique is able to recover all of these with simple WISE-$\textit{Gaia}$ colour cuts that remove $\sim$80 per cent of previously followed-up contaminants. Finally, we provide an online database of known lenses, quasar pairs, and contaminant systems.

Minor mergers are thought to be responsible for the size growth of quiescent field galaxies with decreasing redshift. We test this hypothesis using the cluster environment as a laboratory. Satellite galaxies in clusters move at high velocities, making mergers between them rare. The stellar mass-size relation in ten clusters and in the field is measured and compared at $z~\mathtt{\sim}~1$. Our cluster sample contains 344 spectroscopically-confirmed cluster members with Gemini/GMOS and 182 confirmed with HST WFC3 G141 grism spectroscopy. On average, quiescent and star-forming cluster galaxies are smaller than their field counterparts by ($0.08\pm0.04$) dex and ($0.07\pm0.01$) dex respectively. These size offsets are consistent with the average sizes of quiescent and star-forming field galaxies between $1.2\leqslant z\leqslant1.5$, implying the cluster environment has inhibited size growth between this period and $z~\mathtt{\sim}~1$. The negligible differences measured between the $z~\mathtt{\sim}~0$ field and cluster quiescent mass-size relations in other works imply that the average size of quiescent cluster galaxies must rise with decreasing redshift. Using a toy model, we show that the disappearance of the compact cluster galaxies might be explained if, on average, $\mathtt{\sim}40\%$ of them merge with their brightest cluster galaxies (BCGs) and $\mathtt{\sim}60\%$ are tidally destroyed into the intra-cluster light (ICL) between $0\leqslant z\leqslant1$. This is in agreement with the observed stellar mass growth of BCGs between $0\leqslant z\leqslant1$ and the observed ICL stellar mass fraction at $z~\mathtt{\sim}~0$. Our results support minor mergers as the cause for the size growth in quiescent field galaxies, with cluster-specific processes responsible for the similarity between the field and cluster quiescent mass-size relations at low redshift.

We use a sample of 17 strong gravitational lens systems from the BELLS GALLERY survey to quantify the amount of low-mass dark matter haloes within the lensing galaxies and along their lines of sight, and to constrain the properties of dark matter. Based on a detection criterion of 10$\sigma$, we report no significant detection in any of the lenses. Using the sensitivity function at the 10-$\sigma$ level, we have calculated the predicted number of detectable cold dark matter (CDM) line-of-sight haloes to be $\mu_{l} = 1.17\pm1.08$, in agreement with our null detection. Assuming a detection sensitivity that improved to the level implied by a 5-$\sigma$ threshold, the expected number of detectable line-of-sight haloes rises to $\mu_l = 9.0\pm3.0$. Whilst the current data find zero detections at this sensitivity level (which has a probability of P$^{{\rm5}\sigma}_{{\rm CDM}}(n_{\rm det}=0)$=0.0001 and would be in strong tension with the CDM framework), we find that such a low detection threshold leads to many spurious detections and non-detections and therefore the current lack of detections is unreliable and requires data with improved sensitivity. Combining this sample with a subsample of 11 SLACS lenses, we constrain the half-mode mass to be $\log$(M$_{\rm hm}) < 12.26$ at the 2-$\sigma$ level. The latter is consistent with resonantly produced sterile neutrino masses m$_{\rm s} < 0.8$ keV at any value of the lepton asymmetry at the 2-$\sigma$ level.

We present a study of seventeen LAEs at redshift 2$<z<$3 gravitationally lensed by massive early-type galaxies (ETGs) at a mean redshift of approximately 0.5. Using a fully Bayesian grid-based technique, we model the gravitational lens mass distributions with elliptical power-law profiles and reconstruct the UV-continuum surface brightness distributions of the background sources using pixellated source models. We find that the deflectors are close to, but not consistent with isothermal models in almost all cases, at the $2\sigma$-level. We take advantage of the lensing magnification (typically $\mu\simeq$ 20) to characterise the physical and morphological properties of these LAE galaxies. From reconstructing the ultra-violet continuum emission, we find that the star-formation rates range from 0.3 to 8.5 M$_{\odot}$ yr$^{-1}$ and that the galaxies are typically composed of several compact and diffuse components, separated by 0.4 to 4 kpc. Moreover, they have peak star-formation rate intensities that range from 2.1 to 54.1 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$. These galaxies tend to be extended with major axis ranging from 0.2 to 1.8 kpc (median 561 pc), and with a median ellipticity of 0.49. This morphology is consistent with disk-like structures of star-formation for more than half of the sample. However, for at least two sources, we also find off-axis components that may be associated with mergers. Resolved kinematical information will be needed to confirm the disk-like nature and possible merger scenario for the LAEs in the sample.

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.

The primary goals of the STRong lensing Insights into the Dark Energy Survey (STRIDES) collaboration are to measure the dark energy equation of state parameter and the free streaming length of dark matter. To this aim, STRIDES is discovering strongly lensed quasars in the imaging data of the Dark Energy Survey and following them up to measure time delays, high resolution imaging, and spectroscopy sufficient to construct accurate lens models. In this paper, we first present forecasts for STRIDES. Then, we describe the STRIDES classification scheme, and give an overview of the Fall 2016 follow-up campaign. We continue by detailing the results of two selection methods, the Outlier Selection Technique and a morphological algorithm, and presenting lens models of a system, which could possibly be a lensed quasar in an unusual configuration. We conclude with the summary statistics of the Fall 2016 campaign. Including searches presented in companion papers (Anguita et al.; Ostrovski et al.), STRIDES followed up 117 targets identifying 7 new strongly lensed systems, and 7 nearly identical quasars (NIQs), which could be confirmed as lenses by the detection of the lens galaxy. 76 candidates were rejected and 27 remain otherwise inconclusive, for a success rate in the range 6-35\%. This rate is comparable to that of previous searches like SQLS even though the parent dataset of STRIDES is purely photometric and our selection of candidates cannot rely on spectroscopic information.

For a sample of $\sim$80 local ($0.02 \leq z \leq 0.1$) Seyfert-1 galaxies with high-quality long-slit Keck spectra and spatially-resolved stellar-velocity dispersion ($\sigma_{\star}$) measurements, we study the profile of the [OIII]$\lambda$5007A emission line to test the validity of using its width as a surrogate for $\sigma_{\star}$. Such an approach has often been used in the literature, since it is difficult to measure $\sigma_{\star}$ for type-1 active galactic nuclei (AGNs) due to the AGN continuum outshining the stellar-absorption lines. Fitting the [OIII] line with a single Gaussian or Gauss-Hermite polynomials overestimates $\sigma_{\star}$ by 50-100%. When line asymmetries from non-gravitational gas motion are excluded in a double Gaussian fit, the average ratio between the core [OIII] width ($\sigma_{\rm {[OIII],D}}$) and $\sigma_{\star}$ is $\sim$1, but with individual data points off by up to a factor of two. The resulting black-hole-mass-$\sigma_{\rm {[OIII],D}}$ relation scatters around that of quiescent galaxies and reverberation-mapped AGNs. However, a direct comparison between $\sigma_{\star}$ and $\sigma_{\rm {[OIII],D}}$ shows no close correlation, only that both quantities have the same range, average and standard deviation, probably because they feel the same gravitational potential. The large scatter is likely due to the fact that line profiles are a luminosity-weighted average, dependent on the light distribution and underlying kinematic field. Within the range probed by our sample (80-260 km s$^{-1}$), our results strongly caution against the use of [OIII] width as a surrogate for $\sigma_{\star}$ on an individual basis. Even though our sample consists of radio-quiet AGNs, FIRST radio-detected objects have, on average, a $\sim$10% larger [OIII] core width.

We present milliarcsecond (mas) angular resolution observations of the gravitationally lensed radio source MG J0751+2716 (at z=3.2) obtained with global Very Long Baseline Interferometry (VLBI) at 1.65 GHz. The background object is highly resolved in the tangential and radial directions, showing evidence of both compact and extended structure across several gravitational arcs that are 200 to 600~mas in size. By identifying compact sub-components in the multiple images, we constrain the mass distribution of the foreground z=0.35 gravitational lens using analytic models for the main deflector [power-law elliptical mass model; $\rho(r) \propto r^{-\gamma}$, where $\gamma=2$ corresponds to isothermal] and for the members of the galaxy group. Moreover, our mass models with and without the group find an inner mass-density slope steeper than isothermal for the main lensing galaxy, with $\gamma_1 = 2.08 \pm 0.02$ and $\gamma_2 = 2.16 \pm 0.02$ at the 4.2$\sigma$ level and 6.8$\sigma$ level, respectively, at the Einstein radius ($b_1 = 0.4025 \pm 0.0008$ and $b_2 = 0.307 \pm 0.002$ arcsec, respectively). We find randomly distributed image position residuals of about 3 mas, which are much larger that the measurement errors ($40$ $\mu$as on average). This suggests that at the mas level, the assumption of a smooth mass distribution fails, requiring additional structure in the model. However, given the environment of the lensing galaxy, it is not clear whether this extra mass is in the form of sub-haloes within the lens or along the line of sight, or from a more complex halo for the galaxy group.

We report the discovery, spectroscopic confirmation and preliminary characterisation of 24 gravitationally lensed quasars identified using Gaia observations. Candidates were selected in the Pan-STARRS footprint with quasar-like WISE colours or as photometric quasars from SDSS, requiring either multiple detections in Gaia or a single Gaia detection near a morphological galaxy. The Pan-STARRS grizY images were modelled for the most promising candidates and 60 candidate systems were followed up with the William Herschel Telescope. 13 of the lenses were discovered as Gaia multiples and 10 as single Gaia detections near galaxies. We also discover 1 lens identified through a quasar emission line in an SDSS galaxy spectrum. The lenses have median image separation 2.13 arcsec and the source redshifts range from 1.06 to 3.36. 4 systems are quadruply-imaged and 20 are doubly-imaged. Deep CFHT data reveal an Einstein ring in one double system. We also report 12 quasar pairs, 10 of which have components at the same redshift and require further follow-up to rule out the lensing hypothesis. We compare the properties of these lenses and other known lenses recovered by our search method to a complete sample of simulated lenses to show the lenses we are missing are mainly those with small separations and higher source redshifts. The initial Gaia data release only catalogues all images of ~ 30% of known bright lensed quasars, however the improved completeness of Gaia data release 2 will help find all bright lensed quasars on the sky.

Einstein's theory of gravity, General Relativity, has been precisely tested on Solar System scales, but the long-range nature of gravity is still poorly constrained. The nearby strong gravitational lens, ESO 325-G004, provides a laboratory to probe the weak-field regime of gravity and measure the spatial curvature generated per unit mass, $\gamma$. By reconstructing the observed light profile of the lensed arcs and the observed spatially resolved stellar kinematics with a single self-consistent model, we conclude that $\gamma = 0.97 \pm 0.09$ at 68% confidence. Our result is consistent with the prediction of 1 from General Relativity and provides a strong extragalactic constraint on the weak-field metric of gravity.

We report upon the follow up of 34 candidate lensed quasars found in the Dark Energy Survey using NTT-EFOSC, Magellan-IMACS, KECK-ESI and SOAR-SAMI. These candidates were selected by a combination of double component fitting, morphological assessment and color analysis. Most systems followed up are indeed composed of at least one quasar image and 13 with two or more quasar images: two lenses, four projected binaries and seven Nearly Identical Quasar Pairs (NIQs). The two systems confirmed as genuine gravitationally lensed quasars are one quadruple at $z_s=1.713$ and one double at $z_s=1.515$. Lens modeling of these two systems reveals that both systems require very little contribution from the environment to reproduce the image configuration. Nevertheless, small flux anomalies can be observed in one of the images of the quad. Further observations of 9 inconclusive systems (including 7 NIQs) will allow to confirm (or not) their gravitational lens nature.

We report new observations of SL2SJ021737-051329, a lens system consisting of a bright arc at z=1.84435, magnified ~17x by a massive galaxy at z=0.65. SL2SJ0217 is a low-mass (M <10^9 M*), low-metallicity (Z~1/20 Z*) galaxy, with extreme star-forming conditions that produce strong nebular UV emission lines in the absence of any apparent outflows. Here we present several notable features from rest-frame UV Keck/LRIS spectroscopy: (1) Very strong narrow emission lines are measured for CIV 1548,1550, HeII 1640, OIII] 1661,1666, SiIII] 1883,1892, and CIII] 1907,1909. (2) Double-peaked LyA emission is observed with a dominant blue peak and centered near the systemic velocity. (3) The low- and high-ionization absorption features indicate very little or no outflowing gas along the sightline to the lensed galaxy. The relative emission line strengths can be reproduced with a very high-ionization, low-metallicity starburst with binaries, with the exception of He \ii, which indicates an additional ionization source is needed. We rule out large contributions from AGN and shocks to the photoionization budget, suggesting that the emission features requiring the hardest radiation field likely result from extreme stellar populations that are beyond the capabilities of current models. Therefore, SL2S0217 serves as a template for the extreme conditions that are important for reionization and thought to be more common in the early Universe.

We present models for the dark and luminous mass structure of 12 strong lensing early-type galaxies (ETGs). We combine pixel-based modelling of multiband HST/ACS imaging with Jeans modelling of kinematics obtained from Keck/ESI spectra to disentangle the dark and luminous contributions to the mass. Assuming a gNFW profile for the dark matter halo and a spatially constant stellar-mass-to-light ratio $\Upsilon_{\star}$ for the baryonic mass, we infer distributions for $\Upsilon_{\star}$ consistent with IMFs that are heavier than the Milky Way's (with a global mean mismatch parameter relative to a Chabrier IMF $\mu_{\alpha c} = 1.80 \pm 0.14$) and halo inner density slopes which span a large range but are generally cuspier than the dark-matter-only prediction ($\mu_{\gamma'} = 2.01_{-0.22}^{+0.19}$). We investigate possible reasons for overestimating the halo slope, including the neglect of spatially varying stellar-mas-to-light ratios and/or stellar orbital anisotropy, and find that a quarter of the systems prefer radially declining stellar-mass-to-light ratio gradients, but that the overall effect on our inference on the halo slope is small. We suggest a coherent explanation of these results in the context of inside-out galaxy growth, and that the relative importance of different baryonic processes in shaping the dark halo may depend on halo environment.

We present the first spatially-resolved rest-frame UV study of the gravitationally lensed galaxy, the 'Cosmic Horseshoe' (J1148+1930) at z=2.38. Our gravitational lens model shows that the system is made up of four star-forming regions, each ~4-8 kpc^2 in size, from which we extract four spatially exclusive regional spectra. We study the interstellar and wind absorption lines, along with CIII] doublet emission lines, in each region to investigate any variation in emission/absorption line properties. The mapped CIII] emission shows distinct kinematical structure, with velocity offsets of ~+/-50 km/s between regions suggestive of a merging system, and a variation in equivalent width that indicates a change in ionisation parameter and/or metallicity between the regions. Absorption line velocities reveal a range of outflow strengths, with gas outflowing between -200<v(km/s)<-50 relative to the systemic velocity of that region. Interestingly, the strongest gas outflow appears to emanate from the most diffuse star-forming region. The star-formation rates remain relatively constant (~8-16 M_sol/yr), mostly due to large uncertainties in reddening estimates. As such, the outflows appear to be 'global' rather than 'locally' sourced. We measure electron densities with a range of log(Ne)=3.92-4.36 cm^-3, and point out that such high densities may be common when measured using the CIII] doublet due to its large critical density. Overall, our observations demonstrate that while it is possible to trace variations in large scale gas kinematics, detecting inhomogeneities in physical gas properties and their effects on the outflowing gas may be more difficult. This study provides important lessons for the spatially-resolved rest-frame UV studies expected with future observatories, such as JWST.

We report the discovery of a new gravitationally lensed quasar (double) from the Optical Gravitational Lensing Experiment (OGLE) identified inside the $\sim$670 sq. deg area encompassing the Magellanic Clouds. The source was selected as one of $\sim$60 "red W1-W2" mid-IR objects from WISE and having a significant amount of variability in OGLE for both two (or more) nearby sources. This is the first detection of a gravitational lens, where the discovery is made "the other way around", meaning we first measured the time delay between the two lensed quasar images of $-132<t_{\rm AB}<-76$ days (90% CL), with the median $t_{\rm AB}\approx-102$ days (in the observer frame), and where the fainter image B lags image A. The system consists of the two quasar images separated by 1.5" on the sky, with $I\approx20.0$ mag and $I\approx19.6$ mag, respectively, and a lensing galaxy that becomes detectable as $I \approx 21.5$ mag source, 1.0" from image A, after subtracting the two lensed images. Both quasar images show clear AGN broad emission lines at $z=2.16$ in the NTT spectra. The SED fitting with the fixed source redshift provided the estimate of the lensing galaxy redshift of $z \approx 0.9 \pm 0.2$ (90% CL), while its type is more likely to be elliptical (the SED-inferred and lens-model stellar mass is more likely present in ellipticals) than spiral (preferred redshift by the lens model).

We present the full sample of 118 galaxy-scale strong-lens candidates in the Sloan Lens ACS (SLACS) Survey for the Masses (S4TM) Survey, which are spectroscopically selected from the final data release of the Sloan Digital Sky Survey. Follow-up Hubble Space Telescope (HST) imaging observations confirm that 40 candidates are definite strong lenses with multiple lensed images. The foreground lens galaxies are found to be early-type galaxies (ETGs) at redshifts 0.06 to 0.44, and background sources are emission-line galaxies at redshifts 0.22 to 1.29. As an extension of the SLACS Survey, the S4TM Survey is the first attempt to preferentially search for strong-lens systems with relatively lower lens masses than those in the pre-existing strong-lens samples. By fitting HST data with a singular isothermal ellipsoid model, we find total projected mass within the Einstein radius of the S4TM strong-lens sample ranges from $3 \times10^{10} M_{\odot}$ to $2 \times10^{11} M_{\odot}$. In [Shu15], we have derived the total stellar mass of the S4TM lenses to be $5 \times10^{10} M_{\odot}$ to $1 \times10^{12} M_{\odot}$. Both total enclosed mass and stellar mass of the S4TM lenses are on average almost a factor of 2 smaller than those of the SLACS lenses, which also represent typical mass scales of the current strong-lens samples. The extended mass coverage provided by the S4TM sample can enable a direct test, with the aid of strong lensing, for transitions in scaling relations, kinematic properties, mass structure, and dark-matter content trends of ETGs at intermediate-mass scales as noted in previous studies.

We present the first constraints on stellar mass-to-light ratio gradients in an early-type galaxy (ETG) using multiple dynamical tracer populations to model the dark and luminous mass structure simultaneously. We combine the kinematics of the central starlight, two globular cluster populations and satellite galaxies in a Jeans analysis to obtain new constraints on M87's mass structure, employing a flexible mass model which allows for radial gradients in the stellar mass-to-light ratio. We find that, in the context of our model, a radially declining stellar-mass-to-light ratio is strongly favoured. Modelling the stellar mass-to-light ratio as following a power law, $\Upsilon_{\star} \sim R^{-\mu}$, we infer a power-law slope $\mu = -0.54 \pm 0.05$; equally, parameterising the stellar-mass-to-light ratio via a central mismatch parameter relative to a Salpeter IMF, $\alpha$, and scale radius $R_M$, we find $\alpha > 1.48$ at $95\%$ confidence and $R_M = 0.35 \pm 0.04$ kpc. We use stellar population modelling of high-resolution 11-band HST photometry to show that such a steep gradient cannot be achieved by variations in only the metallicity, age, dust extinction and star formation history if the stellar initial mass function (IMF) remains spatially constant. On the other hand, the stellar mass-to-light ratio gradient that we find is consistent with an IMF whose inner slope changes such that it is Salpeter-like in the central $\sim 0.5$ kpc and becomes Chabrier-like within the stellar effective radius. This adds to recent evidence that the non-universality of the IMF in ETGs may be confined to their core regions, and points towards a picture in which the stars in these central regions may have formed in fundamentally different physical conditions.

We report the discovery, spectroscopic confirmation and first lens models of the first two, strongly lensed quasars from a combined search in WISE and Gaia over the DES footprint. The four-image lensWGD2038-4008 (r.a.=20:38:02.65, dec.=-40:08:14.64) has source- and lens-redshifts $z_{s}=0.777 \pm 0.001$ and $z_l = 0.230 \pm 0.002$ respectively. Its deflector has effective radius $R_{\rm eff} \approx 3.4^{\prime\prime}$, stellar mass $\log(M_{\star}/M_{\odot}) = 11.64^{+0.20}_{-0.43}$, and shows extended isophotal shape variation. Simple lens models yield Einstein radii $R_{\rm E}=(1.30\pm0.04)^{\prime\prime},$ axis ratio $q=0.75\pm0.1$ (compatible with that of the starlight) and considerable shear-ellipticity degeneracies. The two-image lensWGD2021-4115 (r.a.=20:21:39.45, dec.=--41:15:57.11) has $z_{s}=1.390\pm0.001$ and $z_l = 0.335 \pm 0.002$, and Einstein radius $R_{\rm E} = (1.1\pm0.1)^{\prime\prime},$ but higher-resolution imaging is needed to accurately separate the deflector and faint quasar image. We also show high-rank candidate doubles selected this way, some of which have been independently identified with different techniques, and discuss a DES+WISE quasar multiplet selection.

We report the discovery, spectroscopic confirmation, and mass modelling of the gravitationally lensed quasar system PS J0630-1201. The lens was discovered by matching a photometric quasar catalogue compiled from Pan-STARRS and WISE photometry to the Gaia DR1 catalogue, exploiting the high spatial resolution of the latter (FWHM $\sim $0.1") to identify the three brightest components of the lens. Follow-up spectroscopic observations with the WHT confirm the multiple objects are quasars at redshift $z_{q}=3.34$. Further follow-up with Keck AO high-resolution imaging reveals that the system is composed of two lensing galaxies and the quasar is lensed into a $\sim$2.8" separation four-image cusp configuration with a fifth image clearly visible, and a 1.0" arc due to the lensed quasar host galaxy. The system is well-modelled with two singular isothermal ellipsoids, reproducing the position of the fifth image. We discuss future prospects for measuring time delays between the images and constraining any offset between mass and light using the faintly detected Einstein arcs associated with the quasar host galaxy.

Gaia's exceptional resolution (FWHM $\sim$ 0.1$^{\prime\prime}$) allows identification and cataloguing of the multiple images of gravitationally lensed quasars. We investigate a sample of 49 known lensed quasars in the SDSS footprint, with image separations less than 2$^{\prime\prime}$, and find that 8 are detected with multiple components in the first Gaia data release. In the case of the 41 single Gaia detections, we generally are able to distinguish these lensed quasars from single quasars when comparing Gaia flux and position measurements to those of Pan-STARRS and SDSS. This is because the multiple images of these lensed quasars are typically blended in ground-based imaging and therefore the total flux and a flux-weighted centroid are measured, which can differ significantly from the fluxes and centroids of the individual components detected by Gaia. We compare the fluxes through an empirical fit of Pan-STARRS griz photometry to the wide optical Gaia bandpass values using a sample of isolated quasars. The positional offsets are calculated from a recalibrated astrometric SDSS catalogue. Applying flux and centroid difference criteria to spectroscopically confirmed quasars, we discover 4 new sub-arcsecond-separation lensed quasar candidates which have two distinct components of similar colour in archival CFHT or HSC data. Our method based on single Gaia detections can be used to identify the $\sim$ 1400 lensed quasars with image separation above 0.5$^{\prime\prime}$, expected to have only one image bright enough to be detected by Gaia.

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.

Flux ratio anomalies in quasar lenses can be attributed to dark matter substructure surrounding the lensing galaxy and, thus, used to constrain the substructure mass fraction. Previous applications of this approach infer a substructure abundance that potentially in tension with the predictions of a $\Lambda$CDM cosmology. However, the assumption that all flux ratio anomalies are due to substructure is a strong one, and alternative explanations have not been fully investigated. Here, we use new high-resolution near-IR Keck~II adaptive optics imaging for the lens system CLASS B0712+472 to perform pixel-based lens modelling for this system and, in combination with new VLBA radio observations, show that the inclusion of the disc in the lens model can explain the flux ratio anomalies without the need for dark matter substructures. The projected disc mass comprises 16% of the total lensing mass within the Einstein radius and the total disc mass is $1.79 \times 10^{10} M_{sun}$. The case of B0712+472 adds to the evidence that not all flux ratio anomalies are due to dark subhaloes, and highlights the importance of taking the effects of baryonic structures more fully into account in order to obtain an accurate measure of the substructure mass fraction.

We present gravitational lens models of the multiply imaged quasar DES J0408-5354, recently discovered in the Dark Energy Survey (DES) footprint, with the aim of interpreting its remarkable quad-like configuration. We first model the DES single-epoch $grizY$ images as a superposition of a lens galaxy and four point-like objects, obtaining spectral energy distributions (SEDs) and relative positions for the objects. Three of the point sources (A,B,D) have SEDs compatible with the discovery quasar spectra, while the faintest point-like image (G2/C) shows significant reddening and a `grey' dimming of $\approx0.8$mag. In order to understand the lens configuration, we fit different models to the relative positions of A,B,D. Models with just a single deflector predict a fourth image at the location of G2/C but considerably brighter and bluer. The addition of a small satellite galaxy ($R_{\rm E}\approx0.2$") in the lens plane near the position of G2/C suppresses the flux of the fourth image and can explain both the reddening and grey dimming. All models predict a main deflector with Einstein radius between $1.7"$ and $2.0",$ velocity dispersion $267-280$km/s and enclosed mass $\approx 6\times10^{11}M_{\odot},$ even though higher resolution imaging data are needed to break residual degeneracies in model parameters. The longest time-delay (B-A) is estimated as $\approx 85$ (resp. $\approx125$) days by models with (resp. without) a perturber near G2/C. The configuration and predicted time-delays of J0408-5354 make it an excellent target for follow-up aimed at understanding the source quasar host galaxy and substructure in the lens, and measuring cosmological parameters. We also discuss some lessons learnt from J0408-5354 on lensed quasar finding strategies, due to its chromaticity and morphology.

We report the discovery and spectroscopic confirmation of the quad-like lensed quasar system DES J0408-5354 found in the Dark Energy Survey (DES) Year 1 (Y1) data. This system was discovered during a search for DES Y1 strong lensing systems using a method that identified candidates as red galaxies with multiple blue neighbors. DES J0408-5354 consists of a central red galaxy surrounded by three bright (i < 20) blue objects and a fourth red object. Subsequent spectroscopic observations using the Gemini South telescope confirmed that the three blue objects are indeed the lensed images of a quasar with redshift z = 2.375, and that the central red object is an early-type lensing galaxy with redshift z = 0.597. DES J0408-5354 is the first quad lensed quasar system to be found in DES and begins to demonstrate the potential of DES to discover and dramatically increase the sample size of these very rare objects.

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.

We present the discovery and spectroscopic confirmation with the ESO NTT and Gemini South telescopes of eight new 6.0 < z < 6.5 quasars with z$_{AB}$ < 21.0. These quasars were photometrically selected without any star-galaxy morphological criteria from 1533 deg$^{2}$ using SED model fitting to photometric data from the Dark Energy Survey (g, r, i, z, Y), the VISTA Hemisphere Survey (J, H, K) and the Wide-Field Infrared Survey Explorer (W1, W2). The photometric data was fitted with a grid of quasar model SEDs with redshift dependent Lyman-{\alpha} forest absorption and a range of intrinsic reddening as well as a series of low mass cool star models. Candidates were ranked using on a SED-model based $\chi^{2}$-statistic, which is extendable to other future imaging surveys (e.g. LSST, Euclid). Our spectral confirmation success rate is 100% without the need for follow-up photometric observations as used in other studies of this type. Combined with automatic removal of the main types of non-astrophysical contaminants the method allows large data sets to be processed without human intervention and without being over run by spurious false candidates. We also present a robust parametric redshift estimating technique that gives comparable accuracy to MgII and CO based redshift estimators. We find two z $\sim$ 6.2 quasars with HII near zone sizes < 3 proper Mpc which could indicate that these quasars may be young with ages < 10$^6$ - 10$^7$ years or lie in over dense regions of the IGM. The z = 6.5 quasar VDESJ0224-4711 has J$_{AB}$ = 19.75 is the second most luminous quasar known with z > 6.5.

We present an exploration of the mass structure of a sample of 12 strongly lensed massive, compact early-type galaxies at redshifts $z\sim0.6$ to provide further possible evidence for their inside-out growth. We obtain new ESI/Keck spectroscopy and infer the kinematics of both lens and source galaxies, and combine these with existing photometry to construct (a) the fundamental plane (FP) of the source galaxies and (b) physical models for their dark and luminous mass structure. We find their FP to be tilted towards the virial plane relative to the local FP, and attribute this to their unusual compactness, which causes their kinematics to be totally dominated by the stellar mass as opposed to their dark matter; that their FP is nevertheless still inconsistent with the virial plane implies that both the stellar and dark structure of early-type galaxies is non-homologous. We also find the intrinsic scatter of their FP to be comparable to the local value, indicating that variations in the stellar mass structure outweight variations in the dark halo in the central regions of early-type galaxies. Finally, we show that inference on the dark halo structure -- and, in turn, the underlying physics -- is sensitive to assumptions about the stellar initial mass function (IMF), but that physically-motivated assumptions about the IMF imply haloes with sub-NFW inner density slopes, and may present further evidence for the inside-out growth of compact early-type galaxies via minor mergers and accretion.

We present the first results of searches for axions and axion-like-particles with the XENON100 experiment. The axion-electron coupling constant, $g_{Ae}$, has been probed by exploiting the axio-electric effect in liquid xenon. A profile likelihood analysis of 224.6 live days $\times$ 34 kg exposure has shown no evidence for a signal. By rejecting $g_{Ae}$, larger than $7.7 \times 10^{-12}$ (90\% CL) in the solar axion search, we set the best limit to date on this coupling. In the frame of the DFSZ and KSVZ models, we exclude QCD axions heavier than 0.3 eV/c$^2$ and 80 eV/c$^2$, respectively. For axion-like-particles, under the assumption that they constitute the whole abundance of dark matter in our galaxy, we constrain $g_{Ae}$, to be lower than $1 \times 10^{-12}$ (90\% CL) for mass range from 1 to 40 keV/c$^2$, and set the best limit to date as well.

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

We present the discovery and preliminary characterization of a gravitationally lensed quasar with a source redshift $z_{s}=2.74$ and image separation of $2.9"$ lensed by a foreground $z_{l}=0.40$ elliptical galaxy. Since the images of gravitationally lensed quasars are the superposition of multiple point sources and a foreground lensing galaxy, we have developed a morphology independent multi-wavelength approach to the photometric selection of lensed quasar candidates based on Gaussian Mixture Models (GMM) supervised machine learning. Using this technique and $gi$ multicolour photometric observations from the Dark Energy Survey (DES), near IR $JK$ photometry from the VISTA Hemisphere Survey (VHS) and WISE mid IR photometry, we have identified a candidate system with two catalogue components with $i_{AB}=18.61$ and $i_{AB}=20.44$ comprised of an elliptical galaxy and two blue point sources. Spectroscopic follow-up with NTT and the use of an archival AAT spectrum show that the point sources can be identified as a lensed quasar with an emission line redshift of $z=2.739\pm0.003$ and a foreground early type galaxy with $z=0.400\pm0.002$. We model the system as a single isothermal ellipsoid and find the Einstein radius $\theta_E \sim 1.47"$, enclosed mass $M_{enc} \sim 4 \times 10^{11}$M$_{\odot}$ and a time delay of $\sim$52 days. The relatively wide separation, month scale time delay duration and high redshift make this an ideal system for constraining the expansion rate beyond a redshift of 1.

We present a new sample of strong gravitational lens systems where both the foreground lenses and background sources are early-type galaxies. Using imaging from HST/ACS and Keck/NIRC2, we model the surface brightness distributions and show that the sources form a distinct population of massive, compact galaxies at redshifts $0.4 \lesssim z \lesssim 0.7$, lying systematically below the size-mass relation of the global elliptical galaxy population at those redshifts. These may therefore represent relics of high-redshift red nuggets or their partly-evolved descendants. We exploit the magnifying effect of lensing to investigate the structural properties, stellar masses and stellar populations of these objects with a view to understanding their evolution. We model these objects parametrically and find that they generally require two S\'ersic components to properly describe their light profiles, with one more spheroidal component alongside a more envelope-like component, which is slightly more extended though still compact. This is consistent with the hypothesis of the inside-out growth of these objects via minor mergers. We also find that the sources can be characterised by red-to-blue colour gradients as a function of radius which are stronger at low redshift -- indicative of ongoing accretion -- but that their environments generally appear consistent with that of the general elliptical galaxy population, contrary to recent suggestions that these objects are predominantly associated with clusters.

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.

Accurate and precise measurements of the Hubble constant are critical for testing our current standard cosmological model and revealing possibly new physics. With Hubble Space Telescope (HST) imaging, each strong gravitational lens system with measured time delays can allow one to determine the Hubble constant with an uncertainty of $\sim 7\%$. Since HST will not last forever, we explore adaptive-optics (AO) imaging as an alternative that can provide higher angular resolution than HST imaging but has a less stable point spread function (PSF) due to atmospheric distortion. To make AO imaging useful for time-delay-lens cosmography, we develop a method to extract the unknown PSF directly from the imaging of strongly lensed quasars. In a blind test with two mock data sets created with different PSFs, we are able to recover the important cosmological parameters (time-delay distance, external shear, lens mass profile slope, and total Einstein radius). Our analysis of the Keck AO image of the strong lens system RXJ1131-1231 shows that the important parameters for cosmography agree with those based on HST imaging and modeling within 1-$\sigma$ uncertainties. Most importantly, the constraint on the model time-delay distance by using AO imaging with $0.045"$resolution is tighter by $\sim 50\%$ than the constraint of time-delay distance by using HST imaging with $0.09"$when a power-law mass distribution for the lens system is adopted. Our PSF reconstruction technique is generic and applicable to data sets that have multiple nearby point sources, enabling scientific studies that require high-precision models of the PSF.

Gravitational lens flux-ratio anomalies provide a powerful technique for measuring dark matter substructure in distant galaxies. However, before using these flux-ratio anomalies to test galaxy formation models, it is imperative to ascertain that the given anomalies are indeed due to the presence of dark matter substructure and not due to some other component of the lensing galaxy halo or to propagation effects. Here we present the case of CLASS~B1555+375, which has a strong radio-wavelength flux-ratio anomaly. Our high-resolution near-infrared Keck~II adaptive optics imaging and archival Hubble Space Telescope data reveal the lensing galaxy in this system to have a clear edge-on disc component that crosses directly over the pair of images that exhibit the flux-ratio anomaly. We find that simple models that include the disc can reproduce the cm-wavelength flux-ratio anomaly without requiring additional dark matter substructure. Although further studies are required, our results suggest the assumption that all flux-ratio anomalies are due to a population of dark matter sub-haloes may be incorrect, and analyses that do not account for the full complexity of the lens macro-model may overestimate the substructure mass fraction in massive lensing galaxies.

We discuss spatially resolved emission line spectroscopy secured for a total sample of 15 gravitationally lensed star-forming galaxies at a mean redshift of $z\simeq2$ based on Keck laser-assisted adaptive optics observations undertaken with the recently-improved OSIRIS integral field unit (IFU) spectrograph. By exploiting gravitationally lensed sources drawn primarily from the CASSOWARY survey, we sample these sub-L$^{\ast}$ galaxies with source-plane resolutions of a few hundred parsecs ensuring well-sampled 2-D velocity data and resolved variations in the gas-phase metallicity. Such high spatial resolution data offers a critical check on the structural properties of larger samples derived with coarser sampling using multiple-IFU instruments. We demonstrate how kinematic complexities essential to understanding the maturity of an early star-forming galaxy can often only be revealed with better sampled data. Although we include four sources from our earlier work, the present study provides a more representative sample unbiased with respect to emission line strength. Contrary to earlier suggestions, our data indicates a more diverse range of kinematic and metal gradient behavior inconsistent with a simple picture of well-ordered rotation developing concurrently with established steep metal gradients in all but merging systems. Comparing our observations with the predictions of hydrodynamical simulations suggests that gas and metals have been mixed by outflows or other strong feedback processes, flattening the metal gradients in early star-forming galaxies.

A sample of 102 local (0.02 < z < 0.1) Seyfert galaxies with black hole masses MBH > 10^7 M_sun was selected from the Sloan Digital Sky Survey (SDSS) and observed using the Keck 10-m telescope to study the scaling relations between MBH and host galaxy properties. We study profile changes of the broad Hbeta emission line within the ~3-9 year time-frame between the two sets of spectra. The variability of the broad Hbeta emission line is of particular interest, not only since it is used to estimate MBH, but also since its strength and width is used to classify Seyfert galaxies into different types. At least some form of broad-line variability (in either width or flux) is observed in the majority (~66%) of the objects, resulting in a Seyfert-type change for ~38% of the objects, likely driven by variable accretion and/or obscuration. The broad Hbeta line virtually disappears in 3/102 (~3%) extreme cases. We discuss potential causes for these changing-look AGNs. While similar dramatic transitions have previously been reported in the literature, either on a case-by-case basis or in larger samples focusing on quasars at higher redshifts, our study provides statistical information on the frequency of H$\beta$ line variability in a sample of low-redshift Seyfert galaxies.

Following a number of conflicting studies of M87's mass profile, we undertake a dynamical analysis of multiple tracer populations to constrain its mass over a large radius range. We combine stellar kinematics in the central regions with the dynamics of 612 globular clusters out to 200 kpc and satellite galaxies extending to scales comparable with the virial radius. Using a spherical Jeans analysis, we are able to disentangle the mass contributions from the dark and baryonic components and set constraints on the structure of each. Assuming isotropy, we explore four different models for the dark matter halo and find that a centrally-cored dark matter distribution is preferred. We infer a stellar mass-to-light ratio $\Upsilon_{\star,v} = 6.9 \pm 0.1$ -- consistent with a Salpeter-like IMF -- and a core radius $r_c = 67 \pm 20$ kpc. We then introduce anisotropy and find that, while the halo remains clearly cored, the radial stellar anisotropy has a strong impact on both $\Upsilon_{\star,v}$ and the core's radius; here we find $\Upsilon_{\star,v} = 3.50_{-0.36}^{+0.32}$ -- consistent with a Chabrier-like IMF -- and $r_c = 19.00_{-8.34}^{+8.38}$ kpc. Thus the presence of a core at the centre of the dark halo is robust against anisotropy assumptions, while the stellar mass and core size are not. We are able to reconcile previously discrepant studies by showing that modelling the globular cluster data alone leads to the very different inference of a super-NFW cusp, thus highlighting the value of multiple-population modelling, and we point to the possible role of M87's AGN and the cluster environment in forming the central dark matter core.

We report the observation and confirmation of the first group- and cluster-scale strong gravitational lensing systems found in Dark Energy Survey (DES) data. Through visual inspection of data from the Science Verification (SV) season, we identified 53 candidate systems. We then obtained spectroscopic follow-up of 21 candidates using the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South telescope and the Inamori-Magellan Areal Camera and Spectrograph (IMACS) at the Magellan/Baade telescope. With this follow-up, we confirmed six candidates as gravitational lenses: Three of the systems are newly discovered, and the remaining three were previously known. Of the 21 observed candidates, the remaining 15 were either not detected in spectroscopic observations, were observed and did not exhibit continuum emission (or spectral features), or were ruled out as lensing systems. The confirmed sample consists of one group-scale and five galaxy cluster-scale lenses. The lensed sources range in redshift z ~ 0.80-3.2, and in i-band surface brightness i_{SB} ~ 23-25 mag/sq.-arcsec. (2" aperture). For each of the six systems, we estimate the Einstein radius and the enclosed mass, which have ranges ~ 5.0 - 8.6" and ~ 7.5 x 10^{12} - 6.4 x 10^{13} solar masses, respectively.

We present a new photometric catalogue of the rich globular cluster (GC) system around M87, the brightest cluster galaxy in Virgo. Using archival Next Generation Virgo cluster Survey (NGVS) images in the ugriz bands, observed with CFHT/MegaPrime, we perform a careful subtraction of the galaxy's halo light in order to detect objects at small galactocentric radii as well as in the wider field, and find 17620 GC candidates over a radius range from 1.3 kpc to 445 kpc with g < 24 magnitudes. By inferring their colour, radial and magnitude distributions in a Bayesian way, we find that they are well described as a mixture of two GC populations and two distinct contaminant populations, but confirm earlier findings of radius-dependent colour gradients in both GC populations. This is consistent with a picture in which the more enriched GCs reside deeper in the galaxy's potential well, indicating a role for dissipative collapse in the formation of both the red and the blue GCs.

Blue Compact Dwarf (BCD) Galaxies in the nearby Universe provide a means for studying feedback mechanisms and star-formation processes in low-metallicity environments in great detail. Due to their vicinity, these local analogues to young galaxies are well suited for high-resolution studies that would be unfeasible for primordial galaxies in the high-redshift universe. Here we present HST-WFC3 observations of one such BCD, Mrk 71, one of the most powerful local starbursts known, in the light of [O II], He II, Hb, [O III], Ha, and [S II]. At D=3.44 Mpc, this extensive suite of emission line images enables us to explore the chemical and physical conditions of Mrk 71 on ~2 pc scales. Using these high spatial-resolution observations, we use emission line diagnostics to distinguish ionisation mechanisms on a pixel-by-pixel basis and show that despite the previously reported hypersonic gas and super-bubble blow out, the gas in Mrk 71 is photoionised, with no sign of shock-excited emission. Using strong-line metallicity diagnostics, we present the first 'metallicity image' of a galaxy, revealing chemically inhomogeneity on scales of <50 pc. We additionally demonstrate that while chemical structure can be lost at large spatial scales, metallicity-diagnostics can break down on spatial scales smaller than a HII region. HeII emission line images are used to identify up to six Wolf-Rayet stars in Mrk 71, three of which lie on the edge of blow-out region. This study not only demonstrates the benefits of high-resolution spatially-resolved observations in assessing the effects of feedback mechanisms, but also the limitations of fine spatial scales when employing emission-line diagnostics. Both aspects are especially relevant as we enter the era of extremely large telescopes, when observing structure on ~10 pc scales will no longer be limited to the local universe.

Lens time delays are a powerful probe of cosmology, provided that the gravitational potential of the main deflector can be modeled with sufficient precision. Recent work has shown that this can be achieved by detailed modeling of the host galaxies of lensed quasars, which appear as "Einstein Rings" in high resolution images. We carry out a systematic exploration of the high resolution imaging required to exploit the thousands of lensed quasars that will be discovered by current and upcoming surveys with the next decade. Specifically, we simulate realistic lens systems as imaged by the Hubble Space Telescope (HST), James Webb Space Telescope (JWST), and ground based adaptive optics images taken with Keck or the Thirty Meter Telescope (TMT). We compare the performance of these pointed observations with that of images taken by the Euclid (VIS), Wide-Field Infrared Survey Telescope (WFIRST) and Large Synoptic Survey Telescope (LSST) surveys. We use as our metric the precision with which the slope $\gamma'$ of the total mass density profile $\rho_{tot}\propto r^{-\gamma'}$ for the main deflector can be measured. Ideally, we require that the statistical error on $\gamma'$ be less than 0.02, such that it is subdominant to other sources of random and systematic uncertainties. We find that survey data will likely have sufficient depth and resolution to meet the target only for the brighter gravitational lens systems, comparable to those discovered by the SDSS survey. For fainter systems, that will be discovered by current and future surveys, targeted follow-up will be required. However, the exposure time required with upcoming facilitites such as JWST, the Keck Next Generation Adaptive Optics System, and TMT, will only be of order a few minutes per system, thus making the follow-up of hundreds of systems a practical and efficient cosmological probe.

We present spectroscopic confirmation of two new lensed quasars via data obtained at the 6.5m Magellan/Baade Telescope. The lens candidates have been selected from the Dark Energy Survey (DES) and WISE based on their multi-band photometry and extended morphology in DES images. Images of DES J0115-5244 show two blue point sources at either side of a red galaxy. Our long-slit data confirm that both point sources are images of the same quasar at $z_{s}=1.64.$ The Einstein Radius estimated from the DES images is $0.51$". DES J2200+0110 is in the area of overlap between DES and the Sloan Digital Sky Survey (SDSS). Two blue components are visible in the DES and SDSS images. The SDSS fiber spectrum shows a quasar component at $z_{s}=2.38$ and absorption compatible with Mg II and Fe II at $z_{l}=0.799$, which we tentatively associate with the foreground lens galaxy. The long-slit Magellan spectra show that the blue components are resolved images of the same quasar. The Einstein Radius is $0.68$" corresponding to an enclosed mass of $1.6\times10^{11}\,M_{\odot}.$ Three other candidates were observed and rejected, two being low-redshift pairs of starburst galaxies, and one being a quasar behind a blue star. These first confirmation results provide an important empirical validation of the data-mining and model-based selection that is being applied to the entire DES dataset.

Laboratory experiments searching for galactic dark matter particles scattering off nuclei have so far not been able to establish a discovery. We use data from the XENON100 experiment to search for dark matter interacting with electrons. With no evidence for a signal above the low background of our experiment, we exclude a variety of representative dark matter models that would induce electronic recoils. For axial-vector couplings to electrons, we exclude cross-sections above 6x10^(-35) cm^2 for particle masses of m_chi = 2 GeV/c^2. Independent of the dark matter halo, we exclude leptophilic models as explanation for the long-standing DAMA/LIBRA signal, such as couplings to electrons through axial-vector interactions at a 4.4 sigma confidence level, mirror dark matter at 3.6 sigma, and luminous dark matter at 4.6 sigma.

We create a baseline of the black hole (BH) mass (MBH) - stellar-velocity dispersion (sigma) relation for active galaxies, using a sample of 66 local (0.02<z<0.09) Seyfert-1 galaxies, selected from the Sloan Digital Sky Survey (SDSS). Analysis of SDSS images yields AGN luminosities free of host-galaxy contamination and morphological classification. 51/66 galaxies have spiral morphology. 28 bulges have Sersic index n<2 and are considered candidate pseudo bulges, with eight being definite pseudo bulges based on multiple classification criteria met. Only 4/66 galaxies show sign of interaction/merging. High signal-to-noise ratio Keck spectra provide the width of the broad Hbeta emission line free of FeII emission and stellar absorption. AGN luminosity and Hbeta line widths are used to estimate MBH. The Keck-based spatially-resolved kinematics is used to determine stellar-velocity dispersion within the spheroid effective radius. We find that sigma can vary on average by up to 40% across definitions commonly used in the literature, emphasizing the importance of using self-consistent definitions in comparisons and evolutionary studies. The MBH-sigma relation for our Seyfert-1 galaxies has the same intercept and scatter as that of reverberation-mapped AGNs as well as quiescent galaxies, consistent with the hypothesis that our single epoch MBH estimator and sample selection do not introduce significant biases. Barred galaxies, merging galaxies, and those hosting pseudo bulges do not represent outliers in the MBH-sigma relation. This is in contrast with previous work, although no firm conclusion can be drawn due to the small sample size and limited resolution of the SDSS images.

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.

The low-background, VUV-sensitive 3-inch diameter photomultiplier tube R11410 has been developed by Hamamatsu for dark matter direct detection experiments using liquid xenon as the target material. We present the results from the joint effort between the XENON collaboration and the Hamamatsu company to produce a highly radio-pure photosensor (version R11410-21) for the XENON1T dark matter experiment. After introducing the photosensor and its components, we show the methods and results of the radioactive contamination measurements of the individual materials employed in the photomultiplier production. We then discuss the adopted strategies to reduce the radioactivity of the various PMT versions. Finally, we detail the results from screening 216 tubes with ultra-low background germanium detectors, as well as their implications for the expected electronic and nuclear recoil background of the XENON1T experiment.

XENON is a dark matter direct detection project, consisting of a time projection chamber (TPC) filled with liquid xenon as detection medium. The construction of the next generation detector, XENON1T, is presently taking place at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy. It aims at a sensitivity to spin-independent cross sections of $2 \cdot 10^{-47} ~ \mathrm{cm}^{\mathrm{2}}$ for WIMP masses around 50 GeV/c$^{2}$, which requires a background reduction by two orders of magnitude compared to XENON100, the current generation detector. An active system that is able to tag muons and muon-induced backgrounds is critical for this goal. A water Cherenkov detector of $\sim$10 m height and diameter has been therefore developed, equipped with 8 inch photomultipliers and cladded by a reflective foil. We present the design and optimization study for this detector, which has been carried out with a series of Monte Carlo simulations. The muon veto will reach very high detection efficiencies for muons ($>99.5%$) and showers of secondary particles from muon interactions in the rock ($>70%$). Similar efficiencies will be obtained for XENONnT, the upgrade of XENON1T, which will later improve the WIMP sensitivity by another order of magnitude. With the Cherenkov water shield studied here, the background from muon-induced neutrons in XENON1T is negligible.

We present observational results from a new Hubble Space Telescope (HST) Snapshot program to extend the methods of the Sloan Lens ACS (SLACS) Survey to lower lens-galaxy masses. We discover 40 new galaxy-scale strong lenses, which we supplement with 58 previously discovered SLACS lenses. In addition, we determine the posterior PDFs of the Einstein radius for 33 galaxies (18 new and 15 from legacy SLACS data) based on single lensed images. We find a less-than-unity slope of $0.64\pm0.06$ for the $\log_{10} {\sigma}_*$-$\log_{10} {\sigma}_{\rm SIE}$ relation, which corresponds to a 6-$\sigma$ evidence that the total mass-density profile of early-type galaxies varies systematically in the sense of being shallower at higher lens-galaxy velocity dispersions. The trend is only significant when single-image systems are considered, highlighting the importance of including both "lenses" and "non-lenses" for an unbiased treatment of the lens population when extending to lower mass ranges. By scaling simple stellar population models to the HST I-band data, we identify a strong trend of increasing dark-matter fraction at higher velocity dispersions, which can be alternatively interpreted as a trend in the stellar initial mass function (IMF) normalization. Consistent with previous findings and the suggestion of a non-universal IMF, we find that a Salpeter IMF is ruled out for galaxies with velocity dispersion less than $180$ km/s. Considered together, our mass-profile and dark-matter-fraction trends with increasing galaxy mass could both be explained by an increasing relative contribution on kiloparsec scales from a dark-matter halo with a spatial profile more extended than that of the stellar component.

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.

We investigate the cosmic evolution of the black hole (BH) mass -- bulge luminosity relation using a sample of 52 active galaxies at $z \sim 0.36$ and $z \sim 0.57$ in the BH mass range of $10^{7.4-9.1} M_{\odot}$. By consistently applying multi-component spectral and structural decomposition to high-quality Keck spectra and high-resolution HST images, BH masses ($M_{\rm BH}$) are estimated using the H$\beta$ broad emission line combined with the 5100 \AA\ nuclear luminosity, and bulge luminosities ($L_{\rm bul}$) are derived from surface photometry. Comparing the resulting $M_{\rm BH}-L_{\rm bul}$ relation to local active galaxies and taking into account selection effects, we find evolution of the form $M_{\rm BH} / L_{\rm bul} \propto (1+z)^{\gamma}$ with $\gamma=1.8\pm0.7$, consistent with BH growth preceding that of the host galaxies. Including an additional sample of 27 active galaxies with $0.5<z<1.9$ taken from the literature and measured in a consistent way, we obtain $\gamma=0.9\pm0.7$ for the $M_{\rm BH}-L_{\rm bul}$ relation and $\gamma=0.4\pm0.5$ for the $M_{\rm BH}$--total host galaxy luminosity ($L_{\rm host}$) relation. The results strengthen the findings from our previous studies and provide additional evidence for host-galaxy bulge growth being dominated by disk-to-bulge transformation via minor mergers and/or disk instabilities.

Central Galaxies (CGs) in massive halos live in unique environments with formation histories closely linked to that of the host halo. In local clusters they have larger sizes ($R_e$) and lower velocity dispersions (sigma) at fixed stellar mass M_star, and much larger R_e at a fixed $\sigma$ than field and satellite galaxies (non-CGs). Using spectroscopic observations of group galaxies selected from the COSMOS survey, we compare the dynamical scaling relations of early-type CGs and non-CGs at z~0.6, to distinguish possible mechanisms that produce the required evolution. CGs are systematically offset towards larger R_e at fixed $\sigma$ compared to non-CGs with similar M_star. The CG R_e-M_star relation also shows differences, primarily driven by a sub-population (~15%) of galaxies with large $R_e$, while the M_star-sigma relations are indistinguishable. These results are accentuated when double Sersic profiles, which better fit light in the outer regions of galaxies, are adopted. They suggest that even group-scale CGs can develop extended components by these redshifts that can increase total $R_e$ and M_star estimates by factors of ~2. To probe the evolutionary link between our sample and cluster CGs, we also analyze two cluster samples at z~0.6 and z~0. We find similar results for the more massive halos at comparable z, but much more distinct CG scaling relations at low-z. Thus, the rapid, late-time accretion of outer components, perhaps via the stripping and accretion of satellites, would appear to be a key feature that distinguishes the evolutionary history of CGs.

We observed six strongly lensed, radio-loud quasars (MG 0414+0534, CLASS B0712+472, JVAS B1030+074, CLASS B1127+385, CLASS B1152+199, and JVAS B1938+666) in order to identify systems suitable for measuring cosmological parameters using time delays between their multiple images. These systems are in standard two and four image configurations, with B1938 having a faint secondary pair of images. Two separate monitoring campaigns were carried out using the VLA and upgraded JVLA. Lightcurves were extracted for each individual lensed image and analyzed for signs of intrinsic variability. While it was not possible to measure time delays from these data, $\chi^2$-based and structure function tests found evidence for variability in a majority of the lightcurves. B0712 and B1030 had particularly strong variations, exhibiting linear flux trends. These results show that most of these systems should be targeted with followup monitoring campaigns, especially B0712 and B1030.

We present constraints on the equation of state of dark energy, $w$, and the total matter density, $\Omega_{\mathrm{M}}$, derived from the double-source-plane strong lens SDSSJ0946+1006, the first cosmological measurement with a galaxy-scale double-source-plane lens. By modelling the primary lens with an elliptical power-law mass distribution, and including perturbative lensing by the first source, we are able to constrain the cosmological scaling factor in this system to be $\beta^{-1}=1.404 \pm 0.016$, which implies $\Omega_{\mathrm{M}}= 0.33_{-0.26}^{+0.33}$ for a flat $\Lambda$ cold dark matter ($\Lambda$CDM) cosmology. Combining with a cosmic microwave background prior from Planck, we find $w$ = $-1.17^{+0.20}_{-0.21}$ assuming a flat $w$CDM cosmology. This inference shifts the posterior by 1${\sigma}$ and improves the precision by 30 per cent with respect to Planck alone, and demonstrates the utility of combining simple, galaxy-scale multiple-source-plane lenses with other cosmological probes to improve precision and test for residual systematic biases.

Radio mode feedback, in which most of the energy of an active galactic nucleus (AGN) is released in a kinetic form via radio-emitting jets, is thought to play an important role in the maintenance of massive galaxies in the present-day Universe. We study the link between radio emission and the properties of the central black hole in a large sample of local radio galaxies drawn from the Sloan Digital Sky Survey (SDSS), based on the catalogue of Best and Heckman (2012). Our sample is mainly dominated by massive black holes (mostly in the range $10^8-10^9 M_{\odot}$) accreting at very low Eddington ratios (typically $\lambda < 0.01$). In broad agreement with previously reported trends, we find that radio galaxies are preferentially associated with the more massive black holes, and that the radio loudness parameter seems to increase with decreasing Eddington ratio. We compare our results with previous studies in the literature, noting potential biases. The majority of the local radio galaxies in our sample are currently in a radiatively inefficient accretion regime, where kinetic feedback dominates over radiative feedback. We discuss possible physical interpretations of the observed trends in the context of a two-stage feedback process involving a transition in the underlying accretion modes.

We explore the properties of lensing galaxies and lensed faint sources at redshifts between 1.5 and 3.0. Our sample consists of 9 carefully selected strongly-lensed galaxies discovered by the CASSOWARY survey in the Sloan Digital Sky Survey (SDSS) data. We find that, despite some limitations of the original SDSS data, the homogeneous sample of lensing systems can provide a useful insight into lens and source properties. We also explore the limitations of using low-resolution data to model and analyse galaxy-galaxy lensing. We derive the relative alignment of mass and light in fitted lens profiles. The range of magnification extends above 5, hence we are able to analyse potentially small and low-mass galaxies at high redshifts. We confirm the likely evolution of the size-luminosity relation for blue star-forming galaxies as a function of redshift

We present the results of a search for galaxy substructures in a sample of 11 gravitational lens galaxies from the Sloan Lens ACS Survey. We find no significant detection of mass clumps, except for a luminous satellite in the system SDSS J0956+5110. We use these non-detections, in combination with a previous detection in the system SDSS J0946+1006, to derive constraints on the substructure mass function in massive early-type host galaxies with an average redshift z ~ 0.2 and an average velocity dispersion of 270 km/s. We perform a Bayesian inference on the substructure mass function, within a median region of about 32 kpc squared around the Einstein radius (~4.2 kpc). We infer a mean projected substructure mass fraction $f = 0.0076^{+0.0208}_{-0.0052}$ at the 68 percent confidence level and a substructure mass function slope $\alpha$ < 2.93 at the 95 percent confidence level for a uniform prior probability density on alpha. For a Gaussian prior based on Cold Dark Matter (CDM) simulations, we infer $f = 0 .0064^{+0.0080}_{-0.0042}$ and a slope of $\alpha$ = 1.90$^{+0.098}_{-0.098}$ at the 68 percent confidence level. Since only one substructure was detected in the full sample, we have little information on the mass function slope, which is therefore poorly constrained (i.e. the Bayes factor shows no positive preference for any of the two models).The inferred fraction is consistent with the expectations from CDM simulations and with inference from flux ratio anomalies at the 68 percent confidence level.

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

Strong gravitational lenses can be used to detect low mass subhalos, based on deviations in image fluxes and positions from what can be achieved with a smooth mass distribution. So far, this method has been limited by the small number of (radio-loud, microlensing free) systems which can be analysed for the presence of substructure. Using the gravitational lens B1422+231, we demonstrate that adaptive optics integral field spectroscopy can also be used to detect dark substructures. We analyse data obtained with OSIRIS on the Keck I Telescope, using a Bayesian method that accounts for uncertainties relating to the point spread function and image positions in the separate exposures. The narrow-line [OIII] fluxes measured for the lensed images are consistent with those measured in the radio, and show a significant deviation from what would be expected in a smooth mass distribution, consistent with the presence of a perturbing low mass halo. Detailed lens modelling shows that image fluxes and positions are fit significantly better when the lens is modelled as a system containing a single perturbing subhalo in addition to the main halo, rather than by the main halo on its own, indicating the significant detection of substructure.The inferred mass of the subhalo depends on the subhalo mass density profile: the 68% confidence interval for the perturber mass within 600 pc are: 8.2$^{+0.6}_{-0.8}$, 8.2$^{+0.6}_{-1}$ and 7.6$\pm0.3$ $\log_{10}$[M$_{\rm{sub}}$/M$_\odot$] respectively for a singular isothermal sphere, a pseudo-Jaffe, and an NFW mass profile. This method can extend the study of flux ratio anomalies to virtually all quadruply imaged quasars, and therefore offers great potential to improve the determination of the subhalo mass function in the near future.

We present X-shooter observations of CASSOWARY 20 (CSWA 20), a star-forming (SFR ~6 Msol/yr) galaxy at z=1.433, magnified by a factor of 11.5 by the gravitational lensing produced by a massive foreground galaxy at z=0.741. We analysed the integrated physical properties of the HII regions of CSWA 20 using temperature- and density-sensitive emission lines. We find the abundance of oxygen to be ~1/7 of solar, while carbon is ~50 times less abundant than in the Sun. The unusually low C/O ratio may be an indication of a particularly rapid timescale of chemical enrichment. The wide wavelength coverage of X-shooter gives us access to five different methods for determining the metallicity of CSWA 20, three based on emission lines from HII regions and two on absorption features formed in the atmospheres of massive stars. All five estimates are in agreement, within the factor of ~2 uncertainty of each method. The interstellar medium of CSWA 20 only partially covers the star-forming region as viewed from our direction; in particular, absorption lines from neutrals and first ions are exceptionally weak. We find evidence for large-scale outflows of the interstellar medium (ISM) with speeds of up 750 km/s, similar to the values measured in other high-z galaxies sustaining much higher rates of star formation.

The long-standing assumption that the stellar initial mass function (IMF) is universal has recently been challenged by a number of observations. Several studies have shown that a "heavy" IMF (e.g., with a Salpeter-like abundance of low mass stars and thus normalisation) is preferred for massive early-type galaxies, while this IMF is inconsistent with the properties of less massive, later-type galaxies. These discoveries motivate the hypothesis that the IMF may vary (possibly very slightly) across galaxies and across components of individual galaxies (e.g. bulges vs discs). In this paper we use a sample of 19 late-type strong gravitational lenses from the SWELLS survey to investigate the IMFs of the bulges and discs in late-type galaxies. We perform a joint analysis of the galaxies' total masses (constrained by strong gravitational lensing) and stellar masses (constrained by optical and near-infrared colours in the context of a stellar population synthesis [SPS] model, up to an IMF normalisation parameter). Using minimal assumptions apart from the physical constraint that the total stellar mass within any aperture must be less than the total mass within the aperture, we find that the bulges of the galaxies cannot have IMFs heavier (i.e. implying high mass per unit luminosity) than Salpeter, while the disc IMFs are not well constrained by this data set. We also discuss the necessity for hierarchical modelling when combining incomplete information about multiple astronomical objects. This modelling approach allows us to place upper limits on the size of any departures from universality. More data, including spatially resolved kinematics (as in paper V) and stellar population diagnostics over a range of bulge and disc masses, are needed to robustly quantify how the IMF varies within galaxies.

We investigate the relationship between the halo mass, M_200, and concentration, c, for a sample of 26 group- and cluster-scale strong gravitational lenses. In contrast with previous results, we find that these systems are only ~ 0.1 dex more over-concentrated than similar-mass halos from dark matter simulations; the concentration of a halo with M_200 = 10^14 M_sun is log c = 0.78\pm0.05, while simulations of halos with this mass at similar redshifts (z ~ 0.4) predict log c ~ 0.56 - 0.71. We also find that we are unable to make informative inference on the slope of the M_200-c relation in spite of our large sample size; we note that the steep slopes found in previous studies tend to follow the slope in the covariance between M_200 and c, indicating that these results may be measuring the scatter in the data rather than the intrinsic signal. Furthermore, we conclude that our inability to constrain the M_200-c slope is due to a limited range of halo masses, as determined by explicitly modelling our halo mass distribution, and we suggest that other studies may be producing biased results by using an incorrect distribution for their halo masses.

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 exploit the group environment of the CAmbridge Sloan Survey Of Wide ARcs in the skY (CASSOWARY) z=0.3 lens J2158+0257 to measure the group dynamical mass as a complement to the central dynamical and lensing mass constraints. Follow-up spectroscopy of candidate group members is performed using VLT/FORS2. From the resulting N=21 confirmed members we measure the group dynamical mass by calibrating an analytic tracer mass estimator with cosmological simulations. The luminosity weighted line-of-sight velocity dispersion and the Einstein radius of the lens are used as mass probes in the inner regions of the galaxy. Combining these three observational probes allows us to independently constrain the mass and concentration of the dark matter halo, in addition to the total stellar mass of the central galaxy. We find a dark matter halo in remarkably good agreement with simulations (log_10 M_200/M_sun = 14.2 +/- 0.2, c_200 = 4.4 (+1.6, -1.4)) and a stellar mass-to-light ratio which favors a Salpeter initial mass function ((M/L)* = 5.7 +/- 1.2). Our measurement of a normal halo concentration suggests that there is no discrepancy between simulations and observations at the group mass scale. This is in contrast to the cluster mass scale for which a number of studies have claimed over-concentrated halos. While the halo mass is robustly determined, and the halo concentration is not significantly affected by systematics, the resulting stellar mass-to-light ratio is sensitive to the choice of stellar parameters, such as density profile and velocity anisotropy.

High precision cosmological distance measurements towards individual objects such as time delay gravitational lenses or type Ia supernovae are affected by weak lensing perturbations by galaxies and groups along the line of sight. In time delay gravitational lenses, "external convergence," kappa, can dominate the uncertainty in the inferred distances and hence cosmological parameters. In this paper we attempt to reconstruct kappa, due to line of sight structure, using a simple halo model. We use mock catalogues from the Millennium Simulation, and calibrate and compare our reconstructed P(kappa) to ray-traced kappa "truth" values; taking into account realistic observational uncertainties. We find that the reconstruction of kappa provides an improvement in precision of ~50% over galaxy number counts. We find that the lowest-kappa lines of sight have the best constrained P(kappa). In anticipation of large future samples of lenses, we find that selecting the third of the systems with the highest precision kappa estimates gives a sample of unbiased time delay distance measurements with just ~1% uncertainty due to line of sight external convergence effects. Photometric data are sufficient to pre-select the best-constrained lines of sight, and can be done before investment in light-curve monitoring. Conversely, we show that selecting lines of sight with high external shear could, with the reconstruction model presented, induce biases of up to 1% in time delay distance. We find that a major potential source of systematic error is uncertainty in the high mass end of the stellar mass-halo mass relation; this could introduce ~2% biases on the time-delay distance if completely ignored. We suggest areas for the improvement of this general analysis framework (including more sophisticated treatment of high mass structures) that should allow yet more accurate cosmological inferences to be made.

In order to use strong gravitational lens time delays to measure precise and accurate cosmological parameters the effects of mass along the line of sight must be taken into account. We present a method to achieve this by constraining the probability distribution function of the effective line of sight convergence k_ext. The method is based on matching the observed overdensity in the weighted number of galaxies to that found in mock catalogs with k_ext obtained by ray-tracing through structure formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, and redshift. This additional information reduces the uncertainty of k_ext from sigma_k $0.06 to ~0.04 for very overdense lines of sight like that of the system B1608+656. For more common lines of sight, sigma_k is reduced to ~<0.03, corresponding to an uncertainty of ~<3% on distance. This uncertainty has comparable effects on cosmological parameters to that arising from the mass model of the deflector and its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to constrain k_ext almost as well as with spectroscopic redshifts. As an illustration, we apply our method to the system B1608+656. Our most reliable k_ext estimator gives sigma_k=0.047 down from 0.065 using only galaxy counts. Although deeper multi-band observations of the field of B1608+656 are necessary to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to characterize the immediate environment, is an effective way to increase the precision of time-delay cosmography.

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.

Bright gravitationally lensed galaxies provide our most detailed view of galaxies at high redshift. Yet as a result of the small number of ultra-bright z~2 lensed systems with confirmed redshifts, most detailed spectroscopic studies have been limited in their scope. With the goal of increasing the number of bright lensed galaxies available for detailed follow-up, we have undertaken a spectroscopic campaign targeting wide separation (>3 arcsec) galaxy-galaxy lens candidates within the Sloan Digital Sky Survey (SDSS). Building on the earlier efforts of our CASSOWARY survey, we target a large sample of candidate galaxy-galaxy lens systems in SDSS using a well-established search algorithm which identifies blue arc-like structures situated around luminous red galaxies. In this paper, we present a new redshift catalog containing 25 lensed sources in SDSS confirmed through spectroscopic follow-up of candidate galaxy-galaxy lens systems. Included in this new sample are two of the brightest galaxies (r=19.6 and 19.7) galaxies known at z~2, a low metallicity (12 + log (O/H)~8.0) extreme nebular line emitting galaxy at z=1.43, and numerous systems for which detailed follow-up will be possible. The source redshifts span 0.9<z<2.5 (median redshift of 1.9), and their optical magnitudes are in the range 19.6<r<22.3. We present a brief source-by-source discussion of the spectroscopic properties extracted from our confirmatory spectra and discuss some initial science results. With more than 50 gravitationally lensed z>1 galaxies now confirmed within SDSS, it will soon be possible for the first time to develop generalized conclusions from detailed spectroscopic studies of the brightest lensed systems at high redshift.

Recent work has suggested that the stellar initial mass function (IMF) is not universal, but rather is correlated with galaxy stellar mass, stellar velocity dispersion, or morphological type. In this paper, we investigate variations of the IMF within individual galaxies. For this purpose, we use strong lensing and gas kinematics to measure independently the normalisation of the IMF of the bulge and disk components of a sample of 5 massive spiral galaxies with substantial bulge components taken from the SWELLS survey. We find that the stellar mass of the bulges are tightly constrained by the lensing and kinematic data. A comparison with masses based on stellar population synthesis models fitted to optical and near infrared photometry favors a Salpeter-like normalisation of the IMF. Conversely, the disk masses are less well constrained due to degeneracies with the dark matter halo, but are consistent with Milky Way type IMFs in agreement with previous studies. The disks are submaximal at 2.2 disk scale lengths, but due to the contribution of the bulges, the galaxies are baryon dominated at 2.2 disk scale lengths. Globally, our inferred IMF normalisation is consistent with that found for early-type galaxies of comparable stellar mass (> 10^11 M_sun). Our results suggest a non-universal IMF within the different components of spiral galaxies, adding to the well-known differences in stellar populations between disks and bulges.

We present a ready-to-use method to constrain the density distribution in early-type galaxy lenses. Assuming a power-law density profile, then joint use of the virial theorem and the lens equation yields simple formulae for the power-law index (or logarithmic density gradient). Any dependence on orbital anisotropy can be tightly constrained or even erased completely. Our results rely just on surface brightnesses and line-of-sight kinematics, making deprojection unnecessary. We revisit three systems that have already been examined in the literature (the Cosmic Horseshoe, the Jackpot and B1608+656) and provide our estimates. Finally, we show that the method yields a good approximation for the density profile even when the true profile is a broken power-law, albeit with a mild bias towards isothermality.

We present new mass models for the gravitational lens system B1938+666, using multi-wavelength data acquired from Keck adaptive optics (AO) and Hubble Space Telescope (HST) observations. These models are the first results from the Strong-lensing at High Angular Resolution Program (SHARP), a project designed to study known quadruple-image and Einstein ring lenses using high-resolution imaging, in order to probe their mass distributions in unprecedented detail. Here, we specifically highlight differences between AO- and HST-derived lens models, finding that -- at least when the lens and source galaxies are both bright and red, and the system has a high degree of circular symmetry -- AO-derived models place significantly tighter constraints on model parameters. Using this improved precision, we infer important physical properties about the B1938+666 system, including the mass density slope of the lensing galaxy (gamma = 2.045), the projected dark matter mass fraction within the Einstein radius (M_dark/M_lens = 0.55), and the total magnification factor of the source galaxy (~ 13). Additionally, we measure an upper-limit constraint on luminous substructure (M_V > 16.2), based on the non-detection of bright satellite galaxies in all data sets. Finally, we utilize the improved image resolution of the AO data to reveal the presence of faint arcs outside of the primary Einstein ring. The positions and orientations of these arcs raise the intriguing possibility that B1938+666 has a second source galaxy, located at a more distant redshift. However, future work is needed to verify this hypothesis.

We investigate the possibility of constraining the dark energy equation of state by measuring the ratio of Einstein radii in a strong gravitational lens system with two source planes. This quantity is independent of the Hubble parameter and directly measures the growth of angular diameter distances as a function of redshift. We investigate the prospects for a single double source plane system and for a forecast population of systems discovered by re-observing a population of single source lenses already known from a photometrically selected catalogue such as CASSOWARY or from a spectroscopically selected catalogue such as SLACS. We find that constraints comparable to current data-sets (15% uncertainty on the dark equation of state at 68%CL) are possible with a handful of double source plane systems. We also find that the method's degeneracy between Omega_M and w is almost orthogonal to that of CMB and BAO measurements, making this method highly complimentary to current probes.

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.

We infer the normalization and the radial and angular distributions of the number density of satellites of massive galaxies ($\log_{10}[M_{h}^*/M\odot]>10.5$) between redshifts 0.1 and 0.8 as a function of host stellar mass, redshift, morphology and satellite luminosity. Exploiting the depth and resolution of the COSMOS HST images, we detect satellites up to eight magnitudes fainter than the host galaxies and as close as 0.3 (1.4) arcseconds (kpc). Describing the number density profile of satellite galaxies to be a projected power law such that $P(R)\propto R^{\rpower}$, we find $\rpower=-1.1\pm 0.3$. We find no dependency of $\rpower$ on host stellar mass, redshift, morphology or satellite luminosity. Satellites of early-type hosts have angular distributions that are more flattened than the host light profile and are aligned with its major axis. No significant average alignment is detected for satellites of late-type hosts. The number of satellites within a fixed magnitude contrast from a host galaxy is dependent on its stellar mass, with more massive galaxies hosting significantly more satellites. Furthermore, high-mass late-type hosts have significantly fewer satellites than early-type galaxies of the same stellar mass, likely a result of environmental differences. No significant evolution in the number of satellites per host is detected. The cumulative luminosity function of satellites is qualitatively in good agreement with that predicted using subhalo abundance matching techniques. However, there are significant residual discrepancies in the absolute normalization, suggesting that properties other than the host galaxy luminosity or stellar mass determine the number of satellites.

Early-type galaxies (ETGs) are observed to be more compact at z>2 than in the local Universe. Remarkably, much of this size evolution appears to take place in a short (1.8 Gyr) time span between z=2.2 and z=1.3, which poses a serious challenge to hierarchical galaxy formation models where mergers occurring on a similar timescale are the main mechanism for galaxy growth. We compute the merger-driven redshift evolution of stellar mass Mstar\propto(1+z)^aM, half-mass radius Re\propto(1+z)^aR and velocity-dispersion sigma0\propto(1+z)^asigma predicted by concordance Lambda cold dark matter for a typical massive ETG in the redshift range z=1.3-2.2. Neglecting dissipative processes, and thus maximizing evolution in surface density, we find -1.5<aM<-0.6, -1.9<aR<-0.7 and 0.06<asigma<0.22, under the assumption that the accreted satellites are spheroids. It follows that the predicted z=2.2 progenitors of z=1.3 ETGs are significantly less compact (on average a factor of 2 larger Re at given Mstar) than the quiescent galaxies observed at z>2. Furthermore, we find that the scatter introduced in the size-mass correlation by the predicted merger-driven growth is difficult to reconcile with the tightness of the observed scaling law. We conclude that - barring unknown systematics or selection biases in the current measurements - minor and major mergers with spheroids are not sufficient to explain the observed size growth of ETGs within the standard model.

Measuring the relative mass contributions of luminous and dark matter in spiral galaxies is important for understanding their formation and evolution. The combination of a galaxy rotation curve and strong lensing is a powerful way to break the disk-halo degeneracy that is inherent in each of the methods individually. We present an analysis of the 10-image radio spiral lens B1933+503 at z_l=0.755, incorporating (1) new global VLBI observations, (2) new adaptive-optics assisted K-band imaging, (3) new spectroscopic observations for the lens galaxy rotation curve and the source redshift. We construct a three-dimensionally axisymmetric mass distribution with 3 components: an exponential profile for the disk, a point mass for the bulge, and an NFW profile for the halo. The mass model is simultaneously fitted to the kinematics and the lensing data. The NFW halo needs to be oblate with a flattening of a/c=0.33^{+0.07}_{-0.05} to be consistent with the radio data. This suggests that baryons are effective at making the halos oblate near the center. The lensing and kinematics analysis probe the inner ~10 kpc of the galaxy, and we obtain a lower limit on the halo scale radius of 16 kpc (95% CI). The dark matter mass fraction inside a sphere with a radius of 2.2 disk scale lengths is f_{DM,2.2}=0.43^{+0.10}_{-0.09}. The contribution of the disk to the total circular velocity at 2.2 disk scale lengths is 0.76^{+0.05}_{-0.06}, suggesting that the disk is marginally submaximal. The stellar mass of the disk from our modeling is log_{10}(M_{*}/M_{sun}) = 11.06^{+0.09}_{-0.11} assuming that the cold gas contributes ~20% to the total disk mass. In comparison to the stellar masses estimated from stellar population synthesis models, the stellar initial mass function of Chabrier is preferred to that of Salpeter by a probability factor of 7.2.

We construct a fully self-consistent mass model for the lens galaxy J2141 at z=0.14, and use it to improve on previous studies by modelling its gravitational lensing effect, gas rotation curve and stellar kinematics simultaneously. We adopt a very flexible axisymmetric mass model constituted by a generalized NFW dark matter halo and a stellar mass distribution obtained by deprojecting the MGE fit to the high-resolution K'-band LGSAO imaging data of the galaxy, with the (spatially constant) M/L ratio as a free parameter. We model the stellar kinematics by solving the anisotropic Jeans equations. We find that the inner logarithmic slope of the dark halo is weakly constrained (gamma = 0.82^{+0.65}_{-0.54}), and consistent with an unmodified NFW profile. We infer the galaxy to have (i) a dark matter fraction within 2.2 disk radii of 0.28^{+0.15}_{-0.10}, independent of the galaxy stellar population, implying a maximal disk for J2141; (ii) an apparently uncontracted dark matter halo, with concentration c_{-2} = 7.7_{-2.5}^{+4.2} and virial velocity v_{vir} = 242_{-39}^{+44} km/s, consistent with LCDM predictions; (iii) a slightly oblate halo (q_h = 0.75^{+0.27}_{-0.16}), consistent with predictions from baryon-affected models. Comparing the stellar mass inferred from the combined analysis (log_{10} Mstar/Msun = 11.12_{-0.09}^{+0.05}) with that inferred from SPS modelling of the galaxies colours, and accounting for a cold gas fraction of 20+/-10%, we determine a preference for a Chabrier IMF over Salpeter IMF by a Bayes factor of 5.7 (substantial evidence). We infer a value beta_{z} = 1 - sigma^2_{z}/sigma^2_{R} = 0.43_{-0.11}^{+0.08} for the orbital anisotropy parameter in the meridional plane, in agreement with most studies of local disk galaxies, and ruling out at 99% CL that the dynamics of this system can be described by a two-integral distribution function. [Abridged]

We present gravitational lens models for 20 strong gravitational lens systems observed as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS) project. Fifteen of the lenses are taken from paper I while five are newly discovered systems. The systems are galaxy-galaxy lenses where the foreground deflector has an inclined disc, with a wide range of morphological types, from late-type spiral to lenticular. For each system, we compare the total mass inside the critical curve inferred from gravitational lens modelling to the stellar mass inferred from stellar population synthesis (SPS) models, computing the stellar mass fraction f* = M(SPS)/M(lens). We find that, for the lower mass SWELLS systems, adoption of a Salpeter stellar initial mass function (IMF) leads to estimates of f* that exceed 1. This is unphysical, and provides strong evidence against the Salpeter IMF being valid for these systems. Taking the lower mass end of the SWELLS sample sigma(SIE) < 230 km/s, we find that the IMF is lighter (in terms of stellar mass-to-light ratio) than Salpeter with 98% probability, and consistent with the Chabrier IMF and IMFs between the two. This result is consistent with previous studies of spiral galaxies based on independent techniques. In combination with recent studies of massive early-type galaxies that have favoured a heavier Salpeter-like IMF, this result strengthens the evidence against a universal stellar IMF.

The mass-function of dwarf satellite galaxies that are observed around Local Group galaxies substantially differs from simulations based on cold dark matter: the simulations predict many more dwarf galaxies than are seen. The Local Group, however, may be anomalous in this regard. A massive dark satellite in an early-type lens galaxy at z = 0.222 was recently found using a new method based on gravitational lensing, suggesting that the mass fraction contained in substructure could be higher than is predicted from simulations. The lack of very low mass detections, however, prohibited any constraint on their mass function. Here we report the presence of a 1.9 +/- 0.1 x 10^8 M_sun dark satellite in the Einstein-ring system JVAS B1938+666 at z = 0.881, where M_sun denotes solar mass. This satellite galaxy has a mass similar to the Sagittarius galaxy, which is a satellite of the Milky Way. We determine the logarithmic slope of the mass function for substructure beyond the local Universe to be alpha = 1.1^+0.6_-0.4, with an average mass-fraction of f = 3.3^+3.6_-1.8 %, by combining data on both of these recently discovered galaxies. Our results are consistent with the predictions from cold dark matter simulations at the 95 per cent confidence level, and therefore agree with the view that galaxies formed hierarchically in a Universe composed of cold dark matter.

Spectroscopic confirmation of galaxies at z~7 and above has been extremely difficult, owing to a drop in intensity of Ly-alpha emission in comparison with samples at z~6. This crucial finding could potentially signal the ending of cosmic reionization. However it is based on small datasets, often incomplete and heterogeneous in nature. We introduce a flexible Bayesian framework, useful to interpret such evidence. Within this framework, we implement two simple phenomenological models: a smooth one, where the distribution of Ly-alpha is attenuated by a factor \es with respect to z~6; a patchy one where a fraction \ep is absorbed/non-emitted while the rest is unabsorbed. From a compilation of 39 observed z~7 galaxies we find \es=0.69+-0.12 and \ep=0.66+-0.16. The models can be used to compute fractions of emitters above any equivalent width W. For W>25\AA, we find X^{25}_{z=7}=0.37+-0.11 (0.14+-0.06) for galaxies fainter (brighter) than M_{UV}=-20.25 for the patchy model, consistent with previous work, but with smaller uncertainties by virtue of our full use of the data. At z~8 we combine new deep (5-\sigma flux limit 10^{-17}ergs^{-1}cm^{-2}) Keck-NIRSPEC observations of a bright Y-dropout identified by our BoRG Survey, with those of three objects from the literature and find that the inference is inconclusive. We compute predictions for future near-infrared spectroscopic surveys and show that it is challenging but feasible to constrain the distribution of Ly-alpha emitters at z~8 and distinguish between models.

We combine Hubble Space Telescope images from the Great Observatories Origins Deep Survey with archival Very Large Telescope and Keck spectra of a sample of 11 X-ray selected broad-line active galactic nuclei in the redshift range 1<z<2 to study the black hole mass - stellar mass relation out to a lookback time of 10 Gyrs. Stellar masses of the spheroidal component are derived from multi-filter surface photometry. Black hole masses are estimated from the width of the broad MgII emission line and the 3000A nuclear luminosity. Comparing with a uniformly measured local sample and taking into account selection effects, we find evolution in the form M_BH/M_spheroid ~ (1+z)^(1.96+/-0.55), in agreement with our earlier studies based on spheroid luminosity. However, this result is more accurate because it does not require a correction for luminosity evolution and therefore avoids the related and dominant systematic uncertainty. We also measure total stellar masses. Combining our sample with data from the literature, we find M_BH/M_host ~ (1+z)^(1.15+/-0.15), consistent with the hypothesis that black holes (in the range M_BH ~ 10^8-9 M_sun) predate the formation of their host galaxies. Roughly one third of our objects reside in spiral galaxies; none of the host galaxies reveal signs of interaction or major merger activity. Combined with the slower evolution in host stellar masses compared to spheroid stellar masses, our results indicate that secular evolution or minor mergers play a non-negligible role in growing both BHs and spheroids.

The relative contribution of baryons and dark matter to the inner regions of spiral galaxies provides critical clues to their formation and evolution, but it is generally difficult to determine. For spiral galaxies that are strong gravitational lenses, however, the combination of lensing and kinematic observations can be used to break the disk-halo degeneracy. In turn, such data constrain fundamental parameters such as i) the mass density profile slope and axis ratio of the dark matter halo, and by comparison with dark matter-only numerical simulations the modifications imposed by baryons; ii) the mass in stars and therefore the overall star formation efficiency, and the amount of feedback; iii) by comparison with stellar population synthesis models, the normalization of the stellar initial mass function. In this first paper of a series, we present a sample of 16 secure, 1 probable, and 6 possible strong lensing spiral galaxies, for which multi-band high-resolution images and rotation curves were obtained using the Hubble Space Telescope and Keck-II Telescope as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS). The sample includes 8 newly discovered secure systems. [abridged] We find that the SWELLS sample of secure lenses spans a broad range of morphologies (from lenticular to late-type spiral), spectral types (quantified by Halpha emission), and bulge to total stellar mass ratio (0.22-0.85), while being limited to M_*>10^{10.5} M_sun. The SWELLS sample is thus well-suited for exploring the relationship between dark and luminous matter in a broad range of galaxies. We find that the deflector galaxies obey the same size-mass relation as that of a comparison sample of elongated non-lens galaxies selected from the SDSS survey. We conclude that the SWELLS sample is consistent with being representative of the overall population of high-mass high-inclination disky galaxies.

We exploit the strong lensing effect to explore the properties of intrinsically faint and compact galaxies at intermediate redshift, at the highest possible resolution at optical wavelengths. Our sample consists of 46 strongly-lensed emission line galaxies discovered by the Sloan Lens ACS (SLACS) Survey. The galaxies have been imaged at high resolution with HST in three bands (V_HST, I_814 and H_160), allowing us to infer their size, luminosity, and stellar mass using stellar population synthesis models. Lens modeling is performed using a new fast and robust code, klens, which we test extensively on real and synthetic non-lensed galaxies, and also on simulated galaxies multiply-imaged by SLACS- like galaxy-scale lenses. Our tests show that our measurements of galaxy size, flux, and Sersic index are robust and accurate, even for objects intrinsically smaller than the HST point spread function. The median magnification is 8.8, with a long tail that extends to magnifications above 40. Modeling the SLACS sources reveals a population of galaxies with colors and Sersic indices (median n ~ 1) consistent with the objects detected in the field with HST in the GEMS survey, but that are (typically) ~ 2 magnitudes fainter and ~ 5 times smaller in apparent size. The closest analog are ultracompact emission line galaxies identified by HST grism surveys. The lowest mass galaxies in our sample are comparable to the brightest Milky Way satellites in stellar mass (10^7 solar masses) and have well-determined half light radii of 0."05 (~0.3 kpc).

The degeneracy among the disk, bulge and halo contributions to galaxy rotation curves prevents an understanding of the distribution of baryons and dark matter in disk galaxies. In an attempt to break this degeneracy, we present an analysis of the spiral galaxy strong gravitational lens SDSS J2141-0001, discovered as part of the SLACS survey. We present new Hubble Space Telescope multicolor imaging, gas and stellar kinematics data derived from long-slit spectroscopy, and K-band LGS adaptive optics imaging, both from the Keck telescopes. We model the galaxy as a sum of concentric axisymmetric bulge, disk and halo components and infer the contribution of each component, using information from gravitational lensing and gas kinematics. This analysis yields a best-fitting total (disk plus bulge) stellar mass of log_{10}(Mstar/Msun) = 10.99(+0.11,-0.25). The photometric data combined with stellar population synthesis models yield log_{10}(Mstar/Msun) = 10.97\pm0.07, and 11.21\pm0.07 for the Chabrier and Salpeter IMFs, respectively. Accounting for the expected gas fraction of \simeq 20% reduces the lensing plus kinematics stellar mass by 0.10\pm0.05 dex, resulting in a Bayes factor of 11.9 in favor of a Chabrier IMF. The dark matter halo is roughly spherical, with minor to major axis ratio q_{halo}=0.91(+0.15,-0.13). The dark matter halo has a maximum circular velocity of V_{max}=276(+17,-18) km/s, and a central density parameter of log_{10}\Delta_{V/2}=5.9(+0.9,-0.5). This is higher than predicted for uncontracted dark matter haloes in LCDM cosmologies, log_{10}\Delta_{V/2}=5.2, suggesting that either the halo has contracted in response to galaxy formation, or that the halo has a higher than average concentration. At 2.2 disk scale lengths the dark matter fraction is f_{DM}=0.55(+0.20,-0.15), suggesting that SDSS J2141-0001 is sub-maximal.

We study the spatial distribution of faint satellites of intermediate redshift (0.1<z<0.8), early-type galaxies, selected from the GOODS fields. We combine high resolution HST images and state-of-the-art host subtraction techniques to detect satellites of unprecedented faintness and proximity to intermediate redshift host galaxies (up to 5.5 magnitudes fainter and as close as 0."5/2.5 kpc to the host centers). We model the spatial distribution of objects near the hosts as a combination of an isotropic, homogenous background/foreground population and a satellite population with a power law radial profile and an elliptical angular distribution. We detect a significant population of satellites, Ns =1.7 (+0.9,-0.8) that is comparable to the number of Milky Way satellites with similar host-satellite contrast.The average projected radial profile of the satellite distribution is isothermal, gamma_p= -1.0(+0.3,-0.4), which is consistent with the observed central mass density profile of massive early-type galaxies. Furthermore, the satellite distribution is highly anisotropic (isotropy is ruled out at a >99.99% confidence level). Defining phi to be the offset between the major axis of the satellite spatial distribution and the major axis of the host light profile, we find a maximum posterior probability of phi = 0 and |phi| less than 42 degrees at the 68% confidence level. The alignment of the satellite distribution with the light of the host is consistent with simulations, assuming that light traces mass for the host galaxy as observed for lens galaxies. The anisotropy of the satellite population enhances its ability to produce the flux ratio anomalies observed in gravitationally lensed quasars.

We present high-quality Keck/LRIS longslit spectroscopy of a pilot sample of 25 local active galaxies selected from the SDSS (0.02<z<0.1; MBH>10^7 M_sun) to study the relations between black hole mass (MBH) and host-galaxy properties. We determine stellar kinematics of the host galaxy, deriving stellar-velocity dispersion profiles and rotation curves from three spectral regions (including CaH&K, MgIb triplet, and CaII triplet). In addition, we perform surface photometry on SDSS images, using a newly developed code for joint multi-band analysis. BH masses are estimated from the width of the Hbeta emission line and the host-galaxy free 5100A AGN luminosity. Combining results from spectroscopy and imaging allows us to study four MBH scaling relations: MBH-sigma, MBH-L(sph), MBH-M(sph,*), MBH-M(sph,dyn). We find the following results. First, stellar-velocity dispersions determined from aperture spectra (e.g. SDSS fiber spectra or unresolved data from distant galaxies) can be biased, depending on aperture size, AGN contamination, and host-galaxy morphology. However, such a bias cannot explain the offset seen in the MBH-sigma relation at higher redshifts. Second, while the CaT region is the cleanest to determine stellar-velocity dispersions, both the MgIb region, corrected for FeII emission, and the CaHK region, although often swamped by the AGN powerlaw continuum and emission lines, can give results accurate to within a few percent. Third, the MBH scaling relations of our pilot sample agree in slope and scatter with those of other local active and inactive galaxies. In the next papers of the series we will quantify the scaling relations, exploiting the full sample of ~100 objects.

We report new results obtained in calibrations of superheated liquid droplet detectors used in dark matter searches with different radiation sources (n,$\alpha$,$\gamma$). In particular, detectors were spiked with alpha-emitters located inside and outside the droplets. It is shown that the responses are different, depending on whether alpha particles or recoil nuclei create the signals. The energy thresholds for $\alpha$-emitters are compared with test beam measurements using mono-energetic neutrons, as well as with theoretical predictions. Finally a model is presented which describes how the observed intensities of particle induced acoustic signals can be related to the dynamics of bubble growth in superheated liquids. An improved understanding of the bubble dynamics is an important first step in obtaining better discrimination between particle types interacting in detectors of this kind.

We use Keck laser guide star adaptive optics imaging and exploit the magnifying effects of strong gravitational lensing (the effective resolution is FWHM ~ 200 pc) to investigate the sub-kpc scale of an intermediate-redshift (z = 0.63) massive early-type galaxy being lensed by a foreground early-type galaxy; we dub this class of strong gravitational lens systems EELs, e.g., early-type/early-type lenses. We find that the background source is massive (M* = 10^{10.9} M_sun) and compact (r_e = 1.1 kpc), and a two-component fit is required to model accurately the surface brightness distribution, including an extended low-surface-brightness component. This extended component may arise from the evolution of higher-redshift `red nuggets' or may already be in place at z ~ 2 but is unobservable due to cosmological surface brightness dimming.

We present a joint gravitational lensing and stellar-dynamical analysis of 11 early-type galaxies (median deflector redshift $\zd=0.5$) from Strong Lenses in the Legacy Survey (SL2S). Using newly measured redshifts and stellar velocity dispersions from Keck spectroscopy with lens models from Paper I, we derive the total mass density slope inside the Einstein radius for each of the 11 lenses. The average total density slope is found to be $<\gamma'>= 2.16^{+0.09}_{-0.09}$ ($\rho_{\rm tot}\propto r^{-\gamma'}$), with an intrinsic scatter of $0.25^{+0.10}_{-0.07}$. We also determine the dark matter fraction for each lens within half the effective radius, and find the average projected dark matter mass fraction to be $0.42^{+0.08}_{-0.08}$ with a scatter of $0.20^{+0.09}_{-0.07}$ for a Salpeter IMF. By combining the SL2S results with those from the Sloan Lens ACS Survey (median $\zd=0.2$) and the Lenses Structure and Dynamics survey (median $\zd=0.8$), we investigate cosmic evolution of $\gamma'$ and find a mild trend $\partial<\gamma'>/\partial\zd = -0.25^{+0.10}_{-0.12}$. This suggests that the total density profile of massive galaxies has become slightly steeper over cosmic time. If this result is confirmed by larger samples, it would indicate that dissipative processes played some role in the growth of massive galaxies since $z\sim1$.

We use stellar masses, photometry, lensing, and velocity dispersions to investigate empirical correlations for the final sample of 73 early-type lens galaxies (ETGs) from the SLACS survey. The traditional correlations (Fundamental Plane [FP] and its projections) are consistent with those found for non-lens galaxies, supporting the thesis that SLACS lens galaxies are representative of massive ETGs. The addition of strong lensing estimates of the total mass allows us to gain further insights into their internal structure: i) the mean slope of the total mass density profile is <gamma'> = 2.078+/-0.027 with an intrinsic scatter of 0.16+/-0.02; ii) gamma' correlates with effective radius and central mass density, in the sense that denser galaxies have steeper profiles; iii) the dark matter fraction within reff/2 is a monotonically increasing function of galaxy mass and size; iv) the dimensional mass M_dim is proportional to the total mass, and both increase more rapidly than stellar mass M*; v) the Mass Plane (MP), obtained by replacing surface brightness with surface mass density in the FP, is found to be tighter and closer to the virial relation than the FP and the M*P, indicating that the scatter of those relations is dominated by stellar population effects; vi) we construct the Fundamental Hyper-Plane by adding stellar masses to the MP and find the M* coefficient to be consistent with zero and no residual intrinsic scatter. Our results demonstrate that the dynamical structure of ETGs is not scale invariant and that it is fully specified by the total mass, r_eff, and sigma. Although the basic trends can be explained qualitatively in terms of varying star formation efficiency as a function of halo mass and as the result of dry and wet mergers, reproducing quantitatively the observed correlations and their tightness may be a significant challenge for galaxy formation models.

We use stellar dynamics, strong lensing, stellar population synthesis models, and weak lensing shear measurements to constrain the dark matter (DM) profile and stellar mass in a sample of 53 massive early-type galaxies. We explore three DM halo models (unperturbed Navarro Frenk & White [NFW] halos and the adiabatic contraction models of Blumenthal and Gnedin) and impose a model for the relationship between the stellar and virial mass (i.e., a relationship for the star-formation efficiency as a function of halo mass). We show that, given our model assumptions, the data clearly prefer a Salpeter-like initial mass function (IMF) over a lighter IMF (e.g., Chabrier or Kroupa), irrespective of the choice of DM halo. In addition, we find that the data prefer at most a moderate amount of adiabatic contraction (Blumenthal adiabatic contraction is strongly disfavored) and are only consistent with no adiabatic contraction (i.e., a NFW halo) if a mass-dependent IMF is assumed, in the sense of a more massive normalization of the IMF for more massive halos.

We measure the average mass properties of a sample of 41 strong gravitational lenses at moderate redshift (z ~ 0.4 - 0.9), and present the lens redshift for 6 of these galaxies for the first time. Using the techniques of strong and weak gravitational lensing on archival data obtained from the Hubble Space Telescope, we determine that the average mass overdensity profile of the lenses can be fit with a power-law profile (Delta_Sigma prop. to R^{-0.86 +/- 0.16}) that is within 1-sigma of an isothermal profile (Delta_Sigma prop. to R^{-1}) with velocity dispersion sigma_v = 260 +/- 20 km/s. Additionally, we use a two-component de Vaucouleurs+NFW model to disentangle the total mass profile into separate luminous and dark matter components, and determine the relative fraction of each component. We measure the average rest frame V-band stellar mass-to-light ratio (Upsilon_V = 4.0 +/- 0.6 h M_sol/L_sol) and virial mass-to-light ratio (tau_V = 300 +/- 90 h M_sol/L_sol) for our sample, resulting in a virial-to-stellar mass ratio of M_vir/M_* = 75 +/- 25. Finally, we compare our results to a previous study using low redshift lenses, to understand how galaxy mass profiles evolve over time. We investigate the evolution of M_vir/M_*(z) = alpha(1+z)^{beta}, and find best fit parameters of alpha = 51 +/- 36 and beta = 0.9 +/- 1.8, constraining the growth of virial to stellar mass ratio over the last ~7 Gigayears. We note that, by using a sample of strong lenses, we are able to constrain the growth of M_vir/M_*(z) without making any assumptions about the IMF of the stellar population.

We use high-resolution adaptive optics (AO) imaging on the Keck II telescope to study the gravitational lens B0128+437 in unprecedented detail, allowing us to resolve individual lensed quasar components and, for the first time, detect and measure properties of the lensing galaxy. B0128+437 is a small separation lens with known flux-ratio and astrometric anomalies. We discuss possible causes for these anomalies, including the presence of substructure in the lensing galaxy, propagation effects due to dust and a turbulent interstellar medium, and gravitational microlensing. This work on B0128 demonstrates that AO will be an essential tool for studying the many new small-separation lenses expected from future surveys.

We study the internal dynamical structure of the early-type lens galaxy SDSS J0728+3835 at z = 0.206. The analysis is based on two-dimensional kinematic maps extending out to 1.7 effective radii obtained from Keck spectroscopy, on lensing geometry and on stellar mass estimates obtained from multiband Hubble Space Telescope imaging. The data are modelled under the assumptions of axial symmetry supported by a two-integral distribution function (DF), by applying the combined gravitational lensing and stellar dynamics code CAULDRON, and yielding high-quality constraints for an early-type galaxy at cosmological redshifts. Modelling the total density profile as a power-law of the form rho_tot ~ 1/r^{gamma}, we find that it is nearly isothermal (logarithmic slope gamma = 2.08^{+0.04}_{-0.02}), and quite flattened (axial ratio q = 0.60^{+0.08}_{-0.03}). The galaxy is mildly anisotropic (delta = 0.08 +/- 0.02) and shows a fair amount of rotational support, in particular towards the outer regions. We determine a dark matter fraction lower limit of 28 per cent within the effective radius. The stellar contribution to the total mass distribution is close to maximal for a Chabrier initial mass function (IMF), whereas for a Salpeter IMF the stellar mass exceeds the total mass within the galaxy inner regions. We find that the combination of a NFW dark matter halo with the maximally rescaled luminous profile provides a remarkably good fit to the total mass distribution over a broad radial range. Our results confirm and expand the findings of the SLACS survey for early-type galaxies of comparable velocity dispersion (sigma_SDSS = 214 +/- 11 km/s). The internal structure of J0728 is consistent with that of local early-type galaxies of comparable velocity dispersion as measured by the SAURON project, suggesting lack of evolution in the past two billion years.

The gravitational lens system CLASS B2108+213 has two radio-loud lensed images separated by 4.56 arcsec. The relatively large image separation implies that the lensing is caused by a group of galaxies. In this paper, new optical imaging and spectroscopic data for the lensing galaxies of B2108+213 and the surrounding field galaxies are presented. These data are used to investigate the mass and composition of the lensing structure. The redshift and stellar velocity dispersion of the main lensing galaxy (G1) are found to be z = 0.3648 +/- 0.0002 and sigma_v = 325 +/- 25 km/s, respectively. The optical spectrum of the lensed quasar shows no obvious emission or absorption features and is consistent with a BL Lac type radio source. However, the tentative detection of the G-band and Mg-b absorption lines, and a break in the spectrum of the host galaxy of the lensed quasar gives a likely source redshift of z = 0.67. Spectroscopy of the field around B2108+213 finds 51 galaxies at a similar redshift to G1, thus confirming that there is a much larger structure at z ~ 0.365 associated with this system. The width of the group velocity distribution is 694 +/- 93 km/s, but is non-Gaussian, implying that the structure is not yet viralized. The main lensing galaxy is also the brightest group member and has a surface brightness profile consistent with a typical cD galaxy. A lensing and dynamics analysis of the mass distribution, which also includes the newly found group members, finds that the logarithmic slope of the mass density profile is on average isothermal inside the Einstein radius, but steeper at the location of the Einstein radius. This apparent change in slope can be accounted for if an external convergence gradient, representing the underlying parent halo of the galaxy group, is included in the mass model.

Strong gravitational lens systems with measured time delays between the multiple images provide a method for measuring the "time-delay distance" to the lens, and thus the Hubble constant. We present a Bayesian analysis of the strong gravitational lens system B1608+656, incorporating (i) new, deep Hubble Space Telescope (HST) observations, (ii) a new velocity dispersion measurement of 260+/-15 km/s for the primary lens galaxy, and (iii) an updated study of the lens' environment. When modeling the stellar dynamics of the primary lens galaxy, the lensing effect, and the environment of the lens, we explicitly include the total mass distribution profile logarithmic slope gamma' and the external convergence kappa_ext; we marginalize over these parameters, assigning well-motivated priors for them, and so turn the major systematic errors into statistical ones. The HST images provide one such prior, constraining the lens mass density profile logarithmic slope to be gamma'=2.08+/-0.03; a combination of numerical simulations and photometric observations of the B1608+656 field provides an estimate of the prior for kappa_ext: 0.10 +0.08/-0.05. This latter distribution dominates the final uncertainty on H_0. Compared with previous work on this system, the new data provide an increase in precision of more than a factor of two. In combination with the WMAP 5-year data set, we find that the B1608+656 data set constrains the curvature parameter to be -0.031 < Omega_k < 0.009 (95% CL), a level of precision comparable to that afforded by the current Type Ia SNe sample. Asserting a flat spatial geometry, we find that, in combination with WMAP, H_0 = 69.7 +4.9/-5.0 km/s/Mpc and w=-0.94 +0.17/-0.19 (68% CL), suggesting that the observations of B1608+656 constrain w as tightly as do the current Baryon Acoustic Oscillation data. (abridged)

From high-resolution images of 23 Seyfert-1 galaxies at z=0.36 and z=0.57 obtained with the Near Infrared Camera and Multi-Object Spectrometer on board the Hubble Space Telescope (HST), we determine host-galaxy morphology, nuclear luminosity, total host-galaxy luminosity and spheroid luminosity. Keck spectroscopy is used to estimate black hole mass (M_BH). We study the cosmic evolution of the M_BH-spheroid luminosity (L_sph) relation. In combination with our previous work, totaling 40 Seyfert-1 galaxies, the covered range in BH mass is substantially increased, allowing us to determine for the first time intrinsic scatter and correct evolutionary trends for selection effects. We re-analyze archival HST images of 19 local reverberation-mapped active galaxies to match the procedure adopted at intermediate redshift. Correcting spheroid luminosity for passive luminosity evolution and taking into account selection effects, we determine that at fixed present-day V-band spheroid luminosity, M_BH/L_sph \propto (1+z)^(2.8+/-1.2). When including a sample of 44 quasars out to z=4.5 taken from the literature, with luminosity and BH mass corrected to a self-consistent calibration, we extend the BH mass range to over two orders of magnitude, resulting in M_BH/L_sph \propto (1+z)^(1.4+/-0.2). The intrinsic scatter of the relation, assumed constant with redshift, is 0.3+/-0.1 dex (<0.6 dex at 95% CL). The evolutionary trend suggests that BH growth precedes spheroid assembly. Interestingly, the M_BH-total host-galaxy luminosity relation is apparently non-evolving. It hints at either a more fundamental relation or that the spheroid grows by a redistribution of stars. However, the high-z sample does not follow this relation, indicating that major mergers may play the dominant role in growing spheroids above z~1.

We determine an absolute calibration of the initial mass function (IMF) of early-type galaxies, by studying a sample of 56 gravitational lenses identified by the SLACS Survey. Under the assumption of standard Navarro, Frenk & White dark matter halos, a combination of lensing, dynamical, and stellar population synthesis models is used to disentangle the stellar and dark matter contribution for each lens. We define an "IMF mismatch" parameter \alpha=M*(L+D)/M*(SPS) as the ratio of stellar mass inferred by a joint lensing and dynamical models (M*(L+D)) to the current stellar mass inferred from stellar populations synthesis models (M*(SPS)). We find that a Salpeter IMF provides stellar masses in agreement with those inferred by lensing and dynamical models (<\log \alpha>=0.00+-0.03+-0.02), while a Chabrier IMF underestimates them (<\log \alpha>=0.25+-0.03+-0.02). A tentative trend is found, in the sense that \alpha appears to increase with galaxy velocity dispersion. Taken at face value, this result would imply a non universal IMF, perhaps dependent on metallicity, age, or abundance ratios of the stellar populations. Alternatively, the observed trend may imply non-universal dark matter halos with inner density slope increasing with velocity dispersion. While the degeneracy between the two interpretations cannot be broken without additional information, the data imply that massive early-type galaxies cannot have both a universal IMF and universal dark matter halos.

We present the current photometric dataset for the Sloan Lens ACS (SLACS) Survey, including HST photometry from ACS, WFPC2, and NICMOS. These data have enabled the confirmation of an additional 15 grade `A' (certain) lens systems, bringing the number of SLACS grade `A' lenses to 85; including 13 grade `B' (likely) systems, SLACS has identified nearly 100 lenses and lens candidates. Approximately 80% of the grade `A' systems have elliptical morphologies while ~10% show spiral structure; the remaining lenses have lenticular morphologies. Spectroscopic redshifts for the lens and source are available for every system, making SLACS the largest homogeneous dataset of galaxy-scale lenses to date. We have developed a novel Bayesian stellar population analysis code to determine robust stellar masses with accurate error estimates. We apply this code to deep, high-resolution HST imaging and determine stellar masses with typical statistical errors of 0.1 dex; we find that these stellar masses are unbiased compared to estimates obtained using SDSS photometry, provided that informative priors are used. The stellar masses range from 10^10.5 to 10^11.8 M$_\odot$ and the typical stellar mass fraction within the Einstein radius is 0.4, assuming a Chabrier IMF. The ensemble properties of the SLACS lens galaxies, e.g. stellar masses and projected ellipticities, appear to be indistinguishable from other SDSS galaxies with similar stellar velocity dispersions. This further supports that SLACS lenses are representative of the overall population of massive early-type galaxies with M* >~ 10^11 M$_\odot$, and are therefore an ideal dataset to investigate the kpc-scale distribution of luminous and dark matter in galaxies out to z ~ 0.5.

The characteristic size of early-type galaxies (ETGs) of given stellar mass is observed to increase significantly with cosmic time, from redshift z>2 to the present. A popular explanation for this size evolution is that ETGs grow through dissipationless ("dry") mergers, thus becoming less compact. Combining N-body simulations with up-to-date scaling relations of local ETGs, we show that such an explanation is problematic, because dry mergers do not decrease the galaxy stellar-mass surface-density enough to explain the observed size evolution, and also introduce substantial scatter in the scaling relations. Based on our set of simulations, we estimate that major and minor dry mergers increase half-light radius and projected velocity dispersion with stellar mass (M) as M^(1.09+/-0.29) and M^(0.07+/-0.11), respectively. This implies that: 1) if the high-z ETGs are indeed as dense as estimated, they cannot evolve into present-day ETGs via dry mergers; 2) present-day ETGs cannot have assembled more than ~45% of their stellar mass via dry mergers. Alternatively, dry mergers could be reconciled with the observations if there was extreme fine tuning between merger history and galaxy properties, at variance with our assumptions. Full cosmological simulations will be needed to evaluate whether this fine-tuned solution is acceptable.

Based on 58 SLACS strong-lens early-type galaxies with direct total-mass and stellar-velocity dispersion measurements, we find that inside one effective radius massive elliptical galaxies with M_eff >= 3x10^10 M_sun are well-approximated by a power-law ellipsoid with an average logaritmic density slope of <gamma'_LD> = -dlog(rho_tot)/dlog(r)=2.085^{+0.025}_{-0.018} (random error on mean) for isotropic orbits with beta_r=0, +-0.1 (syst.) and sigma_gamma' <= 0.20^{+0.04}_{-0.02} intrinsic scatter (all errors indicate the 68 percent CL). We find no correlation of gamma'_LD with galaxy mass (M_eff), rescaled radius (i.e. R_einst/R_eff) or redshift, despite intrinsic differences in density-slope between galaxies. Based on scaling relations, the average logarithmic density slope can be derived in an alternative manner, fully independent from dynamics, yielding <gamma'_SR>=1.959 +- 0.077. Agreement between the two values is reached for <beta_r> =0.45 +- 0.25, consistent with mild radial anisotropy. This agreement supports the robustness of our results, despite the increase in mass-to-light ratio with total galaxy mass: M_eff ~ L_{V,eff}^(1.363+-0.056). We conclude that massive early-type galaxies are structurally close-to homologous with close-to isothermal total density profiles (<=10 percent intrinsic scatter) and have at most some mild radial anisotropy. Our results provide new observational limits on galaxy formation and evolution scenarios, covering four Gyr look-back time.

Upcoming ground and space based observatories such as the DES, the LSST, the JDEM concepts and the SKA, promise to dramatically increase the size of strong gravitational lens samples. A significant fraction of the systems are expected to be time delay lenses. Many of the existing lensing degeneracies become less of an issue with large samples since the distributions of a number of parameters are predictable, and can be incorporated into an analysis, thus helping to lessen the degeneracy. Assuming a mean galaxy density profile that does not evolve with redshift, a Lambda-CDM cosmology, and Gaussian distributions for bulk parameters describing the lens and source populations, we generate synthetic lens catalogues and examine the relationship between constraints on the Omega_m - Omega_Lambda plane and H_0 with increasing lens sample size. We find that, with sample sizes of ~400 time delay lenses, useful constraints can be obtained for Omega_m and Omega_Lambda with approximately similar levels of precision as from the best of other methods. In addition, sample sizes of ~100 time delay systems yield estimates of H_0 with errors of only a couple of percent, exceeding the level of precision from current best estimates such as the HST Key Project. We note that insufficient prior knowledge of the lens samples employed in the analysis, via under or overestimates in the mean values of the sample distributions, results in broadening of constraints. This highlights the need for sound prior knowledge of the sample before useful cosmological constraints can be obtained from large time delay samples (abridged).

In this paper we address the consumption of energy by astronomers while performing their professional duties. Although we find that astronomy uses a negligible fraction of the US energy budget, the rate at which energy is consumed by an average astronomer is similar to that of a typical high-flying businessperson. We review some of the ways in which astronomers are already acting to reduce their energy consumption. In the coming decades, all citizens will have to reduce their energy consumption to conserve fossil fuel reserves and to help avert a potentially catastrophic change in the Earth's climate. The challenges are the same for astronomers as they are for everyone: decreasing the distances we travel and investing in energy-efficient infrastructure. The high profile of astronomy in the media, and the great public interest in our field, can play a role in promoting energy-awareness to the wider population. Our specific recommendations are therefore to 1) reduce travel when possible, through efficient meeting organization, and by investing in high-bandwidth video conference facilities and virtual-world software, 2) create energy-efficient observatories, computing centers and workplaces, powered by sustainable energy resources, and 3) actively publicize these pursuits.

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

The Cold Dark Matter theory of gravitationally-driven hierarchical structure formation has earned its status as a paradigm by explaining the distribution of matter over large spans of cosmic distance and time. However, its central tenet, that most of the matter in the universe is dark and exotic, is still unproven; the dark matter hypothesis is sufficiently audacious as to continue to warrant a diverse battery of tests. While local searches for dark matter particles or their annihilation signals could prove the existence of the substance itself, studies of cosmological dark matter in situ are vital to fully understand its role in structure formation and evolution. We argue that gravitational lensing provides the cleanest and farthest-reaching probe of dark matter in the universe, which can be combined with other observational techniques to answer the most challenging and exciting questions that will drive the subject in the next decade: What is the distribution of mass on sub-galactic scales? How do galaxy disks form and bulges grow in dark matter halos? How accurate are CDM predictions of halo structure? Can we distinguish between a need for a new substance (dark matter) and a need for new physics (departures from General Relativity)? What is the dark matter made of anyway? We propose that the central tool in this program should be a wide-field optical imaging survey, whose true value is realized with support in the form of high-resolution, cadenced optical/infra-red imaging, and massive-throughput optical spectroscopy.

The optical afterglow of long-duration GRB 071003 is among the brightest yet to be detected from any GRB, with R ~ 12 mag in KAIT observations starting 42 s after the GRB trigger, including filtered detections during prompt emission. However, our high S/N ratio afterglow spectrum displays only extremely weak absorption lines at what we argue is the host redshift of z = 1.60435 - in contrast to the three other, much stronger Mg II absorption systems observed at lower redshifts. Together with Keck adaptive optics observations which fail to reveal a host galaxy coincident with the burst position, our observations suggest a halo progenitor and offer a cautionary tale about the use of Mg II for GRB redshift determination. We present early through late-time observations spanning the electromagnetic spectrum, constrain the connection between the prompt emission and early variations in the light curve (we observe no correlation), and discuss possible origins for an unusual, marked rebrightening that occurs a few hours after the burst: likely either a late-time refreshed shock or a wide-angle secondary jet. Analysis of the late-time afterglow is most consistent with a wind environment, suggesting a massive star progenitor. Together with GRB 070125, this may indicate that a small but significant portion of star formation in the early universe occurred far outside what we consider a normal galactic disk.

We present results from a systematic investigation of the X-ray properties of a sample of moderate redshift (0.3<z<0.6) galaxy groups. These groups were selected not by traditional X-ray or optical search methods, but rather by an association, either physical or along the line of sight, with a strong gravitational lens. We calculate the properties of seven galaxy groups in the fields of six lens systems. Diffuse X-ray emission from the intragroup medium is detected in four of the groups. All of the detected groups have X-ray luminosities greater than 10^{42} h^{-2} erg/s, and lie on the L_X-sigma_v relations defined by local groups and clusters. The upper limits for the non-detections are also consistent with the local L_X-sigma_v relationships. Although the sample size is small and deeper optical and X-ray data are needed, these results suggest that lens-selected groups are similar to X-ray selected samples and, thus, are more massive than the typical poor-group environments of local galaxies.

In the local Universe, high-power radio galaxies live in lower density environments than low-luminosity radio galaxies. If this trend continues to higher redshifts, powerful radio galaxies would serve as efficient probes of moderate redshift groups and poor clusters. Photometric studies of radio galaxies at 0.3 < z < 0.5 suggest that the radio luminosity-environment correlation disappears at moderate redshifts, though this could be the result of foreground/background contamination affecting the photometric measures of environment. We have obtained multi-object spectroscopy of in the fields of 14 lower luminosity (L_1.4GHz < 4x10^24 W/Hz) and higher luminosity (L_1.4GHz > 1.2x10^25 W/Hz) radio galaxies at z ~ 0.3 to spectroscopically investigate the link between the environment and the radio luminosity of radio galaxies at moderate redshifts. Our results support the photometric analyses; there does not appear to be a correlation between the luminosity of a radio galaxy and its environment at moderate redshifts. Hence, radio galaxies are not efficient signposts for group environments at moderate redshifts.

We report on an investigation of the environments of the SLACS sample of gravitational lenses. The local and global environments of the lenses are characterized using SDSS photometry and, when available, spectroscopy. We find that the lens systems that are best modelled with steeper than isothermal density profiles are more likely to have close companions than lenses with shallower than isothermal profiles. This suggests that the profile steepening may be caused by interactions with a companion galaxy as indicated by N-body simulations of group galaxies. The global environments of the SLACS lenses are typical of non-lensing SDSS galaxies with comparable properties to the lenses, and the richnesses of the lens groups are not as strongly correlated with the lens density profiles as the local environments. Furthermore, we investigate the possibility of line-of-sight contamination affecting the lens models but do not find a significant over-density of sources compared to lines of sight without lenses.

We report on an investigation of the SBS 1520+530 gravitational lens system and its environment using archival HST imaging, Keck spectroscopic data, and Keck adaptive-optics imaging. The AO imaging has allowed us to fix the lens galaxy properties with a high degree of precision when performing the lens modeling, and the data indicate that the lens has an elliptical morphology and perhaps a disk. The new spectroscopic data suggest that previous determinations of the lens redshift may be incorrect, and we report an updated, though inconclusive, value z_lens = 0.761. We have also spectroscopically confirmed the existence of several galaxy groups at approximately the redshift of the lens system. We create new models of the lens system that explicitly account for the environment of the lens, and we also include improved constraints on the lensing galaxy from our adaptive-optics imaging. Lens models created with these new data can be well-fit with a steeper than isothermal mass slope (alpha = 2.29, with the density proportional to r^-alpha) if H_0 is fixed at 72 km/s/Mpc; isothermal models require H_0 ~ 50 km/s/Mpc. The steepened profile may indicate that the lens is in a transient perturbed state caused by interactions with a nearby galaxy.

We report the serendipitous discovery of two strong gravitational lens candidates (ACS J160919+6532 and ACS J160910+6532) in deep images obtained with the Advanced Camera for Surveys on the Hubble Space Telescope, each less than 40 arcsec from the previously known gravitational lens system CLASS B1608+656. The redshifts of both lens galaxies have been measured with Keck and Gemini: one is a member of a small galaxy group at z~0.63, which also includes the lensing galaxy in the B1608+656 system, and the second is a member of a foreground group at z~0.43. By measuring the effective radii and surface brightnesses of the two lens galaxies, we infer their velocity dispersions based on the passively evolving Fundamental Plane (FP) relation. Elliptical isothermal lens mass models are able to explain their image configurations within the lens hypothesis, with a velocity dispersion compatible with that estimated from the FP for a reasonable source-redshift range. Based on the large number of massive early-type galaxies in the field and the number-density of faint blue galaxies, the presence of two additional lens systems around CLASS B1608+656 is not unlikely in hindsight. Gravitational lens galaxies are predominantly early-type galaxies, which are clustered, and the lensed quasar host galaxies are also clustered. Therefore, obtaining deep high-resolution images of the fields around known strong lens systems is an excellent method of enhancing the probability of finding additional strong gravitational lens systems.

We report on the results of a spectroscopic survey of the environments of the gravitational lens systems CLASS B1600+434 (z_l = 0.41, z_s = 1.59) and CLASS B2319+051 (z_l = 0.62). The B1600+434 system has a time delay measured for it, and we find the system to lie in a group with a velocity dispersion of 100 km/s and at least six members. B2319+051 has a large group in its immediate foreground with at least 10 members and a velocity dispersion of 460 km/s and another in the background of the lens with a velocity dispersion of 190 km/s. There are several other small groups in the fields of these lens systems, and we describe the properties of these moderate redshift groups. Furthermore, we quantify the effects of these group structures on the gravitational lenses and find a ~5% correction to the derived value of H_0 for B1600+434.