The H0LiCOW collaboration inferred via gravitational lensing time delays a Hubble constant $H_0=73.3^{+1.7}_{-1.8}$ km s$^{-1}{\rm Mpc}^{-1}$, describing deflector mass density profiles by either a power-law or stars plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in $H_0$. We quantify any potential effect of the MST with flexible mass models that are maximally degenerate with H0. Our calculation is based on a new hierarchical approach in which the MST is only constrained by stellar kinematics. The approach is validated on hydrodynamically simulated lenses. We apply the method to the TDCOSMO sample of 7 lenses (6 from H0LiCOW) and measure $H_0=74.5^{+5.6}_{-6.1}$ km s$^{-1}{\rm Mpc}^{-1}$. In order to further constrain the deflector mass profiles, we then add imaging and spectroscopy for 33 strong gravitational lenses from the SLACS sample. For 9 of the SLAC lenses we use resolved kinematics to constrain the stellar anisotropy. From the joint analysis of the TDCOSMO+SLACS sample, we measure $H_0=67.4^{+4.1}_{-3.2}$ km s$^{-1}{\rm Mpc}^{-1}$, assuming that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. While our new analysis does not statistically invalidate the mass profile assumptions by H0LiCOW, and thus their $H_0$ measurement relying on those, it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on $H_0$ derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data, chiefly spatially resolved kinematics of lens galaxies.

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

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

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

In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated datasets presented in paper I. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20\%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ($ |A| < 2\%$) and precision ($< 6\%$ per system), even in the presence of poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix, and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges.

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

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

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

We present a 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.

In time-delay cosmography, three of the key ingredients are 1) determining the velocity dispersion of the lensing galaxy, 2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and 3) estimating the external convergence $\kappa_\mathrm{ext}$ from less massive structures that are not included in the mass model. We present results on all three of these ingredients for two time-delay lensed quasar systems, DES J0408-5354 and WGD 2038-4008. We use the Gemini, Magellan and VLT telescopes to obtain spectra to both measure the stellar velocity dispersions of the main lensing galaxies and to identify the line-of-sight galaxies in these systems. Next, we identify 10 groups in DES J0408-5354 and 2 groups in WGD 2038-4008using a group-finding algorithm. We then identify the most significant galaxy and galaxy-group perturbers using the "flexion shift" criterion. We determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for both of these systems based on our spectroscopy and on the DES-only multiband wide-field observations. Using weighted galaxy counts, calibrated based on the Millennium Simulation, we find that DES J0408-5354 is located in a significantly underdense environment, leading to a tight (width $\sim3\%$), negative-value $\kappa_\mathrm{ext}$ distribution. On the other hand, WGD 2038-4008 is located in an environment of close to unit density, and its low source redshift results in a much tighter $\kappa_\mathrm{ext}$ of $\sim1\%$, as long as no external shear constraints are imposed.

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.

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

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.

In this study we investigate 179 radio-IR galaxies drawn from a sample of spectroscopically-confirmed galaxies that are detected in radio and mid-infrared (MIR) in the redshift range of $0.55 \leq z \leq 1.30$ in the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey. We constrain the Active Galactic Nuclei (AGN) contribution in the total IR luminosity (f$_{\text{AGN}}$), and estimate the AGN luminosity (L$_{\text{AGN}}$) and the star formation rate (SFR) using the CIGALE Spectral Energy Distribution (SED) fitting routine. Based on the f$_{\text{AGN}}$ and radio luminosity, radio-IR galaxies are split into: galaxies that host either high or low f$_{\text{AGN}}$ AGN (high-/low-f$_{\text{AGN}}$), and star forming galaxies with little to no AGN activity (SFGs). We study the colour, stellar mass, radio luminosity, L$_{\text{AGN}}$ and SFR properties of the three radio-IR sub-samples, comparing to a spec-IR sample drawn from spectroscopically-confirmed galaxies that are also detected in MIR. No significant difference between radio luminosity of these sub-samples was found, which could be due to the combined contribution of radio emission from AGN and star formation. We find a positive relationship between L$_{\text{AGN}}$and specific SFR (sSFR) for both AGN sub-samples, strongly suggesting a co-evolution scenario of AGN and SF in these galaxies. A toy model is designed to demonstrate this co-evolution scenario, where we find that, in almost all cases, a rapid quenching timescale is required, which we argue is a signature of AGN quenching. The environmental preference for intermediate/infall regions of clusters/groups remains across the co-evolution scenario, which suggests that galaxies might be in an orbital motion around the cluster/group during the scenario.

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

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 the measurement of the Hubble Constant, $H_0$, with three strong gravitational lens systems. We describe a blind analysis of both PG1115+080 and HE0435-1223 as well as an extension of our previous analysis of RXJ1131-1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP AO effort, with Hubble Space Telescope (HST) imaging, velocity dispersion measurements, and a description of the line-of-sight mass distribution to build an accurate and precise lens mass model. This mass model is then combined with the COSMOGRAIL measured time delays in these systems to determine $H_{0}$. We do both an AO-only and an AO+HST analysis of the systems and find that AO and HST results are consistent. After unblinding, the AO-only analysis gives $H_{0}=82.8^{+9.4}_{-8.3}~\rm km\,s^{-1}\,Mpc^{-1}$ for PG1115+080, $H_{0}=70.1^{+5.3}_{-4.5}~\rm km\,s^{-1}\,Mpc^{-1}$ for HE0435-1223, and $H_{0}=77.0^{+4.0}_{-4.6}~\rm km\,s^{-1}\,Mpc^{-1}$ for RXJ1131-1231. The joint AO-only result for the three lenses is $H_{0}=75.6^{+3.2}_{-3.3}~\rm km\,s^{-1}\,Mpc^{-1}$. The joint result of the AO+HST analysis for the three lenses is $H_{0}=76.8^{+2.6}_{-2.6}~\rm km\,s^{-1}\,Mpc^{-1}$. All of the above results assume a flat $\Lambda$ cold dark matter cosmology with a uniform prior on $\Omega_{\textrm{m}}$ in [0.05, 0.5] and $H_{0}$ in [0, 150] $\rm km\,s^{-1}\,Mpc^{-1}$. This work is a collaboration of the SHARP and H0LiCOW teams, and shows that AO data can be used as the high-resolution imaging component in lens-based measurements of $H_0$. The full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper.

We present a blind time-delay strong lensing (TDSL) cosmographic analysis of the doubly imaged quasar SDSS 1206+4332. We combine the relative time delay between the quasar images, Hubble Space Telescope imaging, the Keck stellar velocity dispersion of the lensing galaxy, and wide-field photometric and spectroscopic data of the field to constrain two angular diameter distance relations. The combined analysis is performed by forward modelling the individual data sets through a Bayesian hierarchical framework, and it is kept blind until the very end to prevent experimenter bias. After unblinding, the inferred distances imply a Hubble constant $H_0 = 68.8^{+5.4}_{-5.1}$ kms$^{-1}$Mpc$^{-1}$, assuming a flat Lambda cold dark matter cosmology with uniform prior on $\Omega_{\rm m}$ in [0.05, 0.5]. The precision of our cosmographic measurement with the doubly imaged quasar SDSS 1206+4332 is comparable with those of quadruply imaged quasars and opens the path to perform on selected doubles the same analysis as anticipated for quads. Our analysis is based on a completely independent lensing code than our previous three H0LiCOW systems and the new measurement is fully consistent with those. We provide the analysis scripts paired with the publicly available software to facilitate independent analysis. The consistency between blind measurements with independent codes provides an important sanity check on lens modelling systematics. By combining the likelihoods of the four systems under the same prior, we obtain $H_0 = 72.5^{+2.1}_{-2.3}$kms$^{-1}$Mpc$^{-1}$. This measurement is independent of the distance ladder and other cosmological probes.

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

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

We present results of a systematic search for gravitationally lensed quasars in Pan-STARRS1. Our final sample of candidates comprises of 91 systems, not including 25 rediscovered lensed quasars and quasar pairs. In the absence of spectroscopy to verify the lensing nature of the candidates, the main sources of contaminants are likely to be quasar pairs, which we consider to be a byproduct of our work, and a smaller number of quasar$+$star associations. Amongst the independently discovered quads is 2M1134$-$2103, for which we obtained spectroscopy for the first time, finding a redshift of 2.77 for the quasar. There is evidence for microlensing in at least one image. We perform detailed mass modeling of this system using archival imaging data, and find that the unusually large shear responsible for the diamond-like configuration can be attributed mainly to a faint companion $\sim 4''$ away, and to a galaxy group/cluster $\sim 30''$ away. We also set limits of $z\sim0.5-1.5$ on the redshift of the lensing galaxy, based on its brightness, the image separation of the lensed images, and an analysis of the observed photometric flux ratios.

For nearly 40 years, dark matter has been widely assumed to be cold and collisionless. Cold dark matter models make fundamental predictions for the behavior of dark matter on small (<10 kpc) scales. These predictions include cuspy density profiles at the centers of dark matter halos and a halo mass function that increases as dN/dM ~ M^-1.9 down to very small masses. We suggest two observational programs relying on extremely large telescopes to critically test these predictions, and thus shed new light on the nature of dark matter. (1) Combining adaptive optics-enabled imaging with deep spectroscopy to measure the three-dimensional motions of stars within a sample of Local Group dwarf galaxies that are the cleanest dark matter laboratories known in the nearby universe. From these observations the inner slope of the dark matter density profile can be determined with an accuracy of 0.20 dex, enabling a central cusp to be distinguished from a core at 5 sigma significance. (2) Diffraction-limited AO imaging and integral field spectroscopy of gravitationally lensed galaxies and quasars to quantify the abundance of dark substructures in the halos of the lens galaxies and along the line of sight. Observations of 50 lensed arcs and 50 multiply-imaged quasars will be sufficient to measure the halo mass function over the range 10^7 < M < 10^10 Msun at cosmological scales, independent of the baryonic and stellar composition of those structures. These two observational probes provide complementary information about the small scale structure, with a joint self-consistent analysis mitigating limitations of either probe. This program will produce the strongest existing constraints on the properties of dark matter on small scales, allowing conclusive tests of alternative warm, fuzzy, and self-interacting dark matter models.

Many of the fundamental physical constants in Physics, as a discipline, are measured to exquisite levels of precision. The fundamental constants that define Cosmology, however, are largely determined via a handful of independent techniques that are applied to even fewer datasets. The history of the measurement of the Hubble Constant (H0), which serves to anchor the expansion history of the Universe to its current value, is an exemplar of the difficulties of cosmological measurement; indeed, as we approach the centennial of its first measurement, the quest for H0 still consumes a great number of resources. In this white paper, we demonstrate how the approaching era of Extremely Large Telescopes (ELTs) will transform the astrophysical measure of H0 from the limited and few into a fundamentally new regime where (i) multiple, independent techniques are employed with modest use of large aperture facilities and (ii) 1% or better precision is readily attainable. This quantum leap in how we approach H0 is due to the unparalleled sensitivity and spatial resolution of ELT's and the ability to use integral field observations for simultaneous spectroscopy and photometry, which together permit both familiar and new techniques to effectively by-pass the conventional 'ladder' framework to minimize total uncertainty. Three independent techniques are discussed -- (i) standard candles via a two-step distance ladder applied to metal, poor stellar populations, (ii) standard clocks via gravitational lens cosmography, and (iii) standard sirens via gravitational wave sources -- each of which can reach 1% with relatively modest investment from 30-m class facilities.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

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.

Using 57 Radio Active Galactic nuclei (RAGN) at 0.55 $\leq$ z $\leq$ 1.3 drawn from five fields of the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey, we study the effect of injection of energy from outbursts of RAGN on their spectroscopically-confirmed neighboring galaxies (SNGs). We observe an elevated fraction of quenched neighbors (fq) within 500 kpc projected radius of RAGN in the most dense local environments compared to those of non-RAGN control samples matched to the RAGN population in colour, stellar mass, and local environment at 2$\sigma$ significance. Further analyses show that there are offsets at similar significance between fqs of RAGN-SNGs and the appropriate control samples for galaxies specifically in cluster environments and those hosted by most massive cluster galaxies, which tentatively suggests that some negative feedback from the RAGN is occurring in these dense environments. In addition, we find that the median radio power of RAGN increases with increasing local overdensity, an effect which may lend itself to the quenching of neighboring galaxies. Furthermore, we find that, in the highest local overdensities, the fq of the sub-sample of lower stellar mass RAGN-SNGs is larger than that of the higher stellar mass RAGN-SNGs sub-sample, which indicates a more pronounced effect from RAGN on lower stellar mass galaxies. We propose a scenario in which RAGN residing within clusters might heat the intracluster medium (ICM) affecting both in situ star formation and any inflowing gas that remains in their neighboring galaxies.

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.

Time-delay strong lensing provides a unique way to directly measure the Hubble constant ($H_{0}$). The precision of the $H_{0}$ measurement depends on the uncertainties in the time-delay measurements, the mass distribution of the main deflector(s), and the mass distribution along the line of sight. Tie and Kochanek (2018) have proposed a new microlensing effect on time delays based on differential magnification of the coherent accretion disc variability of the lensed quasar. If real, this effect could significantly broaden the uncertainty on the time delay measurements by up to $30\%$ for lens systems such as PG1115+080, which have relatively short time delays and monitoring over several different epochs. In this paper we develop a new technique that uses the time-delay ratios and simulated microlensing maps within a Bayesian framework in order to limit the allowed combinations of microlensing delays and thus to lessen the uncertainties due to the proposed effect. We show that, under the assumption of Tie and Kochanek (2018), the uncertainty on the time-delay distance ($D_{\Delta t}$, which is proportional to 1/$H_{0}$) of short time-delay ($\sim18$ days) lens, PG1115+080, increases from $\sim7\%$ to $\sim10\%$ by simultaneously fitting the three time-delay measurements from the three different datasets across twenty years, while in the case of long time-delay ($\sim90$ days) lens, the microlensing effect on time delays is negligible as the uncertainty on $D_{\Delta t}$ of RXJ1131-1231 only increases from $\sim2.5\%$ to $\sim2.6\%$.

We present 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 evaluated the effectiveness of different indicators of cluster virialization using 12 large-scale structures in the ORELSE survey spanning from $0.7<z<1.3$. We located diffuse X-ray emission from 16 galaxy clusters using Chandra observations. We studied the properties of these clusters and their members, using Chandra data in conjunction with optical and near-IR imaging and spectroscopy. We measured X-ray luminosities and gas temperatures of each cluster, as well as velocity dispersions of their member galaxies. We compared these results to scaling relations derived from virialized clusters, finding significant offsets of up to 3-4$\sigma$ for some clusters, which could indicate they are disturbed or still forming. We explored if other properties of the clusters correlated with these offsets by performing a set of tests of virialization and substructure on our sample, including Dressler-Schectman tests, power ratios, analyses of the velocity distributions of galaxy populations, and centroiding differences. For comparison to a wide range of studies, we used two sets of tests: ones that did and did not use spectral energy distribution fitting to obtain rest-frame colours, stellar masses, and photometric redshifts of galaxies. Our results indicated that the difference between the stellar mass or light mean-weighted center and the X-ray center, as well as the projected offset of the most-massive/brightest cluster galaxy from other cluster centroids had the strongest correlations with scaling relation offsets, implying they are the most robust indicators of cluster virialization and can be used for this purpose when X-ray data is insufficiently deep for reliable $L_X$ and $T_X$ measurements.

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 present time-delay estimates for the quadruply imaged quasar PG 1115+080. Our resuls are based on almost daily observations for seven months at the ESO MPIA 2.2m telescope at La Silla Observatory, reaching a signal-to-noise ratio of about 1000 per quasar image. In addition, we re-analyse existing light curves from the literature that we complete with an additional three seasons of monitoring with the Mercator telescope at La Palma Observatory. When exploring the possible source of bias we consider the so-called microlensing time delay, a potential source of systematic error so far never directly accounted for in previous time-delay publications. In fifteen years of data on PG 1115+080, we find no strong evidence of microlensing time delay. Therefore not accounting for this effect, our time-delay estimates on the individual data sets are in good agreement with each other and with the literature. Combining the data sets, we obtain the most precise time-delay estimates to date on PG 1115+080, with Dt(AB) = 8.3+1.5-1.6 days (18.7% precision), Dt(AC) = 9.9+1.1-1.1 days (11.1%) and Dt(BC) = 18.8+1.6-1.6 days (8.5%). Turning these time delays into cosmological constraints is done in a companion paper that makes use of ground-based Adaptive Optics (AO) with the Keck telescope.

Constraining the sub-galactic matter-power spectrum on 1-10 kpc scales would make it possible to distinguish between the concordance $\Lambda$CDM model and various alternative dark-matter models due to the significantly different levels of predicted mass structure. Here, we demonstrate a novel approach to observationally constrain the population of overall law-mass density fluctuations in the inner regions of massive elliptical lens galaxies, based on the power spectrum of the associated surface-brightness perturbations observable in highly magnified galaxy-scale Einstein rings and gravitational arcs. The application of our method to the SLACS lens system SDSS J0252+0039 results in the following limits (at the 99 per cent confidence level) on the dimensionless convergence-power spectrum (and the associated standard deviation in aperture mass): $\Delta^{2}_{\delta\kappa}<1$ ($\sigma_{AM}< 0.8 \times 10^8 M_\odot$) on 0.5-kpc scale, $\Delta^{2}_{\delta\kappa}<0.1$ ($\sigma_{AM}< 1 \times 10^8 M_\odot$) on 1-kpc scale and $\Delta^{2}_{\delta\kappa}<0.01$ ($\sigma_{AM}< 3 \times 10^8 M_\odot$) on 3-kpc scale. The estimated effect of CDM sub-haloes lies considerably below these first observational upper-limit constraints on the level of inhomogeneities in the projected total mass distribution of galactic haloes. Future analysis for a larger sample of galaxy-galaxy strong lens systems will narrow down these constraints and rule out all cosmological models predicting a significantly larger level of clumpiness on these critical sub-galactic scales.

The flux ratios in the multiple images of gravitationally lensed quasars can provide evidence for dark matter substructure in the halo of the lensing galaxy if the flux ratios differ from those predicted by a smooth model of the lensing galaxy mass distribution. However, it is also possible that baryonic structures in the lensing galaxy, such as edge-on discs, can produce flux-ratio anomalies. In this work, we present the first statistical analysis of flux-ratio anomalies due to baryons from a numerical simulation perspective. We select galaxies with various morphological types in the Illustris simulation and ray-trace through the simulated halos, which include baryons in the main lensing galaxies but exclude any substructures, in order to explore the pure baryonic effects. Our ray-tracing results show that the baryonic components can be a major contribution to the flux-ratio anomalies in lensed quasars and that edge-on disc lenses induce the strongest anomalies. We find that the baryonic components increase the probability of finding high flux-ratio anomalies in the early-type lenses by about 8% and by about 10 - 20% in the disc lenses. The baryonic effects also induce astrometric anomalies in 13% of the mock lenses. Our results indicate that the morphology of the lens galaxy becomes important in the analysis of flux-ratio anomalies when considering the effect of baryons, and that the presence of baryons may also partially explain the discrepancy between the observed (high) anomaly frequency and what is expected due to the presence of subhalos as predicted by the CDM simulations.

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

According to the current understanding of cosmic structure formation, the precursors of the most massive structures in the Universe began to form shortly after the Big Bang, in regions corresponding to the largest fluctuations in the cosmic density field. Observing these structures during their period of active growth and assembly - the first few hundred million years of the Universe - is challenging because it requires surveys that are sensitive enough to detect the distant galaxies that act as signposts for these structures and wide enough to capture the rarest objects. As a result, very few such objects have been detected so far. Here we report observations of a far-infrared-luminous object at redshift 6.900 (less than 800 Myr after the Big Bang) that was discovered in a wide-field survey. High-resolution imaging reveals this source to be a pair of extremely massive star-forming galaxies. The larger of these galaxies is forming stars at a rate of 2900 solar masses per year, contains 270 billion solar masses of gas and 2.5 billion solar masses of dust, and is more massive than any other known object at a redshift of more than 6. Its rapid star formation is probably triggered by its companion galaxy at a projected separation of just 8 kiloparsecs. This merging companion hosts 35 billion solar masses of stars and has a star-formation rate of 540 solar masses per year, but has an order of magnitude less gas and dust than its neighbor and physical conditions akin to those observed in lower-metallicity galaxies in the nearby Universe. These objects suggest the presence of a dark-matter halo with a mass of more than 400 billion solar masses, making it among the rarest dark-matter haloes that should exist in the Universe at this epoch.

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

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.

We present time-delay measurements for the new quadruply imaged quasar DES J0408-5354, the first quadruply imaged quasar found in the Dark Energy Survey (DES). Our result is made possible by implementing a new observational strategy using almost daily observations with the MPIA 2.2m telescope at La Silla observatory and deep exposures reaching a signal-to-noise ratio of about 1000 per quasar image. This data quality allows us to catch small photometric variations (a few mmag rms) of the quasar, acting on temporal scales much shorter than microlensing, hence making the time delay measurement very robust against microlensing. In only 7 months we measure very accurately one of the time delays in DES J0408-5354: Dt(AB) = -112.1 +- 2.1 days (1.8%) using only the MPIA 2.2m data. In combination with data taken with the 1.2m Euler Swiss telescope, we also measure two delays involving the D component of the system Dt(AD) = -155.5 +- 12.8 days (8.2%) and Dt(BD) = -42.4 +- 17.6 days (41%), where all the error bars include systematics. Turning these time delays into cosmological constraints will require deep HST imaging or ground-based Adaptive Optics (AO), and information on the velocity field of the lensing galaxy.

In this study we investigate 89 radio galaxies that are spectroscopically-confirmed to be members of five large scale structures in the redshift range of $0.65 \le z \le 0.96$. Based on a two-stage classification scheme, the radio galaxies are classified into three sub-classes: active galactic nucleus (AGN), hybrid, and star-forming galaxy (SFG). We study the properties of the three radio sub-classes and their global and local environmental preferences. We find AGN hosts are the most massive population and exhibit quiescence in their star-formation activity. The SFG population has a comparable stellar mass to those hosting a radio AGN but are unequivocally powered by star formation. Hybrids, though selected as an intermediate population in our classification scheme, were found in almost all analyses to be a unique type of radio galaxies rather than a mixture of AGN and SFGs. They are dominated by a high-excitation radio galaxy (HERG) population. We discuss environmental effects and scenarios for each sub-class. AGN tend to be preferentially located in locally dense environments and in the cores of clusters/groups, with these preferences persisting when comparing to galaxies of similar colour and stellar mass, suggesting that their activity may be ignited in the cluster/group virialized core regions. Conversely, SFGs exhibit a strong preference for intermediate-density global environments, suggesting that dusty starbursting activity in LSSs is largely driven by galaxy-galaxy interactions and merging.

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

We present the discovery of 3 quasar lenses in the Sloan Digital Sky Survey (SDSS), selected using two novel photometry-based selection techniques. The J0941+0518 system, with two point sources separated by 5.46" on either side of a galaxy, has source and lens redshifts $z_s = 1.54$ and $z_l = 0.343$. The AO-assisted images of J2211+1929 show two point sources separated by 1.04", corresponding to the same quasar at $z_s = 1.07,$ besides the lens galaxy and Einstein ring. Images of J2257+2349 show two point sources separated by 1.67" on either side of an E/S0 galaxy. The extracted spectra show two images of the same quasar at redshift $z_s = 2.10$. In total, the two selection techniques identified 309 lens candidates, including 47 known lenses, and 6 previously ruled out candidates. 55 of the remaining candidates were observed using NIRC2 and ESI at Keck Observatory, EFOSC2 at the ESO-NTT (La Silla), and SAM and the Goodman spectrograph at SOAR. Of the candidates observed, 3 were confirmed as lenses, 36 were ruled out, and 16 remain inconclusive. Taking into account that we recovered known lenses, this gives us a success rate of at least 50/309 (16%). This initial campaign demonstrates the power of purely photometric selection techniques in finding lensed quasars. Developing and refining these techniques is essential for efficient identification of these rare lenses in ongoing and future photometric surveys.

We report the discovery and constrain the physical conditions of the interstellar medium of the highest-redshift millimeter-selected dusty star-forming galaxy (DSFG) to date, SPT-S J031132-5823.4 (hereafter SPT0311-58), at $z=6.900 +/- 0.002$. SPT0311-58 was discovered via its 1.4mm thermal dust continuum emission in the South Pole Telescope (SPT)-SZ survey. The spectroscopic redshift was determined through an ALMA 3mm frequency scan that detected CO(6-5), CO(7-6) and [CI](2-1), and subsequently confirmed by detections of CO(3-2) with ATCA and [CII] with APEX. We constrain the properties of the ISM in SPT0311-58 with a radiative transfer analysis of the dust continuum photometry and the CO and [CI] line emission. This allows us to determine the gas content without ad hoc assumptions about gas mass scaling factors. SPT0311-58 is extremely massive, with an intrinsic gas mass of $M_{\rm gas} = 3.3 \pm 1.9 \times10^{11}\,M_{\odot}$. Its large mass and intense star formation is very rare for a source well into the Epoch of Reionization.

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

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

We present a detection in pre-explosion Hubble Space Telescope (HST) imaging of a point source consistent with being the progenitor star of the Type IIb supernova (SN IIb) 2016gkg. Post-explosion imaging from the Keck Adaptive Optics system was used to perform relative astrometry between the Keck and HST imaging. We identify a single point source in the HST images coincident with the SN position to 0.89-sigma. The HST photometry is consistent with the progenitor star being an A0Ia star with T=9500 K and log (L/Lsun)=5.15. We find that the SN 2016gkg progenitor star appears more consistent with binary than single-star evolutionary models. In addition, early-time light curve data from SN 2016gkg revealed a rapid rise in luminosity within ~0.4 days of non-detection limits, consistent with models of the cooling phase after shock break-out. We use these data to determine an explosion date of 20.15 September 2016 and progenitor star radius of log (R/Rsun)=2.41, which agrees with photometry from the progenitor star. Our findings are also consistent with detections of other SNe IIb progenitor stars, although more luminous and bluer than most other examples.

We studied AGN activity in twelve LSSs in the ORELSE survey, at 0.65<z<1.28, using a combination of Chandra observations, optical and NIR imaging and spectroscopy. We located a total of 61 AGNs across our sample that were successfully matched to optical counterparts in the LSSs. Seeking to study AGN triggering mechanisms, we examined the spatial distribution of the AGNs and their average spectral properties. We found that AGN populations across our sample had less time since the last starburst than the overall galaxy populations. We did not find any relation between AGN activity and location within the LSSs, suggesting triggering mechanisms which depend on global environment are at most sub-dominant. To focus on differences between our AGNs, we grouped them into four sub-samples based on the spectral properties of their parents LSSs. We found one of the sub-samples, SG0023 & SC1604, stood out from the others: AGNs in this sample were disproportionately luminous, their average time since the last starburst event was the smallest, despite the fact that this was not true of the overall galaxy population in those LSSs, and both the AGNs and the overall galaxy population had the largest fraction of close kinematic pairs, which indicates a higher rate of galaxy mergers and interactions. These results suggest that major mergers are driving AGN activity in SG0023 & SC1604, while other processes are likely triggering less luminous AGNs in the rest of our sample. Additionally, minor mergers are unlikely to play a significant role, since the same conditions that lead to more major mergers should should also lead to more minor mergers, which is not observed in SG0023 & SC1604.

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.

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

The South Pole Telescope has discovered one hundred gravitationally lensed, high-redshift, dusty, star-forming galaxies (DSFGs). We present 0.5" resolution 870um Atacama Large Millimeter/submillimeter Array imaging of a sample of 47 DSFGs spanning z=1.9-5.7, and construct gravitational lens models of these sources. Our visibility-based lens modeling incorporates several sources of residual interferometric calibration uncertainty, allowing us to properly account for noise in the observations. At least 70% of the sources are strongly lensed by foreground galaxies (mu_870um > 2), with a median magnification mu_870um = 6.3, extending to mu_870um > 30. We compare the intrinsic size distribution of the strongly lensed sources to a similar number of unlensed DSFGs and find no significant differences in spite of a bias between the magnification and intrinsic source size. This may indicate that the true size distribution of DSFGs is relatively narrow. We use the source sizes to constrain the wavelength at which the dust optical depth is unity and find this wavelength to be correlated with the dust temperature. This correlation leads to discrepancies in dust mass estimates of a factor of 2 compared to estimates using a single value for this wavelength. We investigate the relationship between the [CII] line and the far-infrared luminosity and find that the same correlation between the [CII]L_FIR ratio and Sigma_FIR found for low-redshift star-forming galaxies applies to high-redshift galaxies and extends at least two orders of magnitude higher in Sigma_FIR. This lends further credence to the claim that the compactness of the IR-emitting region is the controlling parameter in establishing the "[CII] deficit."

We use the Atacama Large Millimeter/submillimeter Array (ALMA) in Cycle 1 to determine spectroscopic redshifts of high-redshift dusty star-forming galaxies (DSFGs) selected by their 1.4mm continuum emission in the South Pole Telescope (SPT) survey. We present ALMA 3mm spectral scans between 84-114GHz for 15 galaxies and targeted ALMA 1mm observations for an additional eight sources. Our observations yield 30 new line detections from CO, [CI] , [NII] , H_2O and NH_3. We further present APEX [CII] and CO mid-J observations for seven sources for which only a single line was detected in spectral-scan data from ALMA Cycle 0 or Cycle 1. We combine the new observations with previously published and new mm/submm line and photometric data of the SPT-selected DSFGs to study their redshift distribution. The combined data yield 39 spectroscopic redshifts from molecular lines, a success rate of >85%. Our sample represents the largest data set of its kind today and has the highest spectroscopic completeness among all redshift surveys of high-z DSFGs. The median of the redshift distribution is z=3.9+/-0.4, and the highest-redshift source in our sample is at z=5.8. We discuss how the selection of our sources affects the redshift distribution, focusing on source brightness, selection wavelength, and strong gravitational lensing. We correct for the effect of gravitational lensing and find the redshift distribution for 1.4mm-selected sources with a median redshift of z=3.1+/-0.3. Comparing to redshift distributions selected at shorter wavelengths from the literature, we show that selection wavelength affects the shape of the redshift distribution.

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 study the abundance of substructure in the matter density near galaxies using ALMA Science Verification observations of the strong lensing system SDP.81. We present a method to measure the abundance of subhalos around galaxies using interferometric observations of gravitational lenses. Using simulated ALMA observations, we explore the effects of various systematics, including antenna phase errors and source priors, and show how such errors may be measured or marginalized. We apply our formalism to ALMA observations of SDP.81. We find evidence for the presence of a $M=10^{8.96\pm 0.12} M_{\odot}$ subhalo near one of the images, with a significance of $6.9\sigma$ in a joint fit to data from bands 6 and 7; the effect of the subhalo is also detected in both bands individually. We also derive constraints on the abundance of dark matter subhalos down to $M\sim 2\times 10^7 M_{\odot}$, pushing down to the mass regime of the smallest detected satellites in the Local Group, where there are significant discrepancies between the observed population of luminous galaxies and predicted dark matter subhalos. We find hints of additional substructure, warranting further study using the full SDP.81 dataset (including, for example, the spectroscopic imaging of the lensed carbon monoxide emission). We compare the results of this search to the predictions of $\Lambda$CDM halos, and find that given current uncertainties in the host halo properties of SDP.81, our measurements of substructure are consistent with theoretical expectations. Observations of larger samples of gravitational lenses with ALMA should be able to improve the constraints on the abundance of galactic substructure.

To understand cosmic mass assembly in the Universe at early epochs, we primarily rely on measurements of stellar mass and star formation rate of distant galaxies. In this paper, we present stellar masses and star formation rates of six high-redshift ($2.8\leq z \leq 5.7$) dusty, star-forming galaxies (DSFGs) that are strongly gravitationally lensed by foreground galaxies. These sources were first discovered by the South Pole Telescope (SPT) at millimeter wavelengths and all have spectroscopic redshifts and robust lens models derived from ALMA observations. We have conducted follow-up observations, obtaining multi-wavelength imaging data, using {\it HST}, {\it Spitzer}, {\it Herschel} and the Atacama Pathfinder EXperiment (APEX). We use the high-resolution {\it HST}/WFC3 images to disentangle the background source from the foreground lens in {\it Spitzer}/IRAC data. The detections and upper limits provide important constraints on the spectral energy distributions (SEDs) for these DSFGs, yielding stellar masses, IR luminosities, and star formation rates (SFRs). The SED fits of six SPT sources show that the intrinsic stellar masses span a range more than one order of magnitude with a median value $\sim$ 5 $\times 10^{10}M_{\Sun}$. The intrinsic IR luminosities range from 4$\times 10^{12}L_{\Sun}$ to 4$\times 10^{13}L_{\Sun}$. They all have prodigious intrinsic star formation rates of 510 to 4800 $M_{\Sun} {\rm yr}^{-1}$. Compared to the star-forming main sequence (MS), these six DSFGs have specific SFRs that all lie above the MS, including two galaxies that are a factor of 10 higher than the MS. Our results suggest that we are witnessing the ongoing strong starburst events which may be driven by major mergers.

We present spatially-resolved imaging obtained with the Australia Telescope Compact Array (ATCA) of three CO lines in two high-redshift gravitationally lensed dusty star-forming galaxies, discovered by the South Pole Telescope. Strong lensing allows us to probe the structure and dynamics of the molecular gas in these two objects, at z=2.78 and z=5.66, with effective source-plane resolution of less than 1kpc. We model the lensed emission from multiple CO transitions and the dust continuum in a consistent manner, finding that the cold molecular gas as traced by low-J CO always has a larger half-light radius than the 870um dust continuum emission. This size difference leads to up to 50% differences in the magnification factor for the cold gas compared to dust. In the z=2.78 galaxy, these CO observations confirm that the background source is undergoing a major merger, while the velocity field of the other source is more complex. We use the ATCA CO observations and comparable resolution Atacama Large Millimeter/submillimeter Array dust continuum imaging of the same objects to constrain the CO-H_2 conversion factor with three different procedures, finding good agreement between the methods and values consistent with those found for rapidly star-forming systems. We discuss these galaxies in the context of the star formation - gas mass surface density relation, noting that the change in emitting area with observed CO transition must be accounted for when comparing high-redshift galaxies to their lower redshift counterparts.

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.

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 TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

Strong gravitational lenses provide an important tool to measure masses in the distant Universe, thus testing models for galaxy formation and dark matter; to investigate structure at the Epoch of Reionization; and to measure the Hubble constant and possibly w as a function of redshift. However, the limiting factor in all of these studies has been the currently small samples of known gravitational lenses (~10^2). The era of the SKA will transform our understanding of the Universe with gravitational lensing, particularly at radio wavelengths where the number of known gravitational lenses will increase to ~10^5. Here we discuss the technical requirements, expected outcomes and main scientific goals of a survey for strong gravitational lensing with the SKA. We find that an all-sky (3pi sr) survey carried out with the SKA1-MID array at an angular resolution of 0.25-0.5 arcsec and to a depth of 3 microJy / beam is required for studies of galaxy formation and cosmology with gravitational lensing. In addition, the capability to carryout VLBI with the SKA1 is required for tests of dark matter and studies of supermassive black holes at high redshift to be made using gravitational lensing.

We present [CII] observations of 20 strongly lensed dusty star forming galaxies at 2.1 < z < 5.7 using APEX and Herschel. The sources were selected on their 1.4 mm flux (S_1.4mm > 20 mJy) from the South Pole Telescope survey, with far-infrared (FIR) luminosities determined from extensive photometric data. The [CII] line is robustly detected in 17 sources, all but one being spectrally resolved. Eleven out of 20 sources observed in [CII] also have low-J CO detections from ATCA. A comparison with mid- and high-J CO lines from ALMA reveals consistent [CII] and CO velocity profiles, suggesting that there is little differential lensing between these species. The [CII], low-J CO and FIR data allow us to constrain the properties of the interstellar medium. We find [CII] to CO(1-0) luminosity ratios in the SPT sample of 5200 +- 1800, with significantly less scatter than in other samples. This line ratio can be best described by a medium of [CII] and CO emitting gas with a higher [CII] than CO excitation temperature, high CO optical depth tau_CO >> 1, and low to moderate [CII] optical depth tau_CII ~< 1. The geometric structure of photodissociation regions allows for such conditions.

The time delays between point-like images in gravitational lens systems can be used to measure cosmological parameters. The number of lenses with measured time delays is growing rapidly; the upcoming \emph{Large Synoptic Survey Telescope} (LSST) will monitor $\sim10^3$ strongly lensed quasars. In an effort to assess the present capabilities of the community to accurately measure the time delays, and to provide input to dedicated monitoring campaigns and future LSST cosmology feasibility studies, we have invited the community to take part in a "Time Delay Challenge" (TDC). The challenge is organized as a set of "ladders," each containing a group of simulated datasets to be analyzed blindly by participating teams. Each rung on a ladder consists of a set of realistic mock observed lensed quasar light curves, with the rungs' datasets increasing in complexity and realism. The initial challenge described here has two ladders, TDC0 and TDC1. TDC0 has a small number of datasets, and is designed to be used as a practice set by the participating teams. The (non-mandatory) deadline for completion of TDC0 was the TDC1 launch date, December 1, 2013. The TDC1 deadline was July 1 2014. Here we give an overview of the challenge, we introduce a set of metrics that will be used to quantify the goodness-of-fit, efficiency, precision, and accuracy of the algorithms, and we present the results of TDC0. Thirteen teams participated in TDC0 using 47 different methods. Seven of those teams qualified for TDC1, which is described in the companion paper II.

We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $\chi^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $\chi^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f = $20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.

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.

Herschel Space Observatory photometry and extensive multiwavelength followup have revealed that the powerful radio galaxy 3C 220.3 at z=0.685 acts as a gravitational lens for a background submillimeter galaxy (SMG) at z=2.221. At an observed wavelength of 1mm, the SMG is lensed into three distinct images. In the observed near infrared, these images are connected by an arc of 1.8" radius forming an Einstein half-ring centered near the radio galaxy. In visible light, only the arc is apparent. 3C 220.3 is the only known instance of strong galaxy-scale lensing by a powerful radio galaxy not located in a galaxy cluster and therefore it offers the potential to probe the dark matter content of the radio galaxy host. Lens modeling rejects a single lens, but two lenses centered on the radio galaxy host A and a companion B, separated by 1.5", provide a fit consistent with all data and reveal faint candidates for the predicted fourth and fifth images. The model does not require an extended common dark matter halo, consistent with the absence of extended bright X-ray emission on our Chandra image. The projected dark matter fractions within the Einstein radii of A (1.02") and B (0.61") are about 0.4 +/- 0.3 and 0.55 +/- 0.3. The mass to i-band light ratios of A and B, M/L ~ 8 +/- 4 Msun/Lsun, appear comparable to those of radio-quiet lensing galaxies at the same redshift in the CASTLES, LSD, and SL2S samples. The lensed SMG is extremely bright with observed f(250um) = 440mJy owing to a magnification factor mu~10. The SMG spectrum shows luminous, narrow CIV 154.9nm emission, revealing that the SMG houses a hidden quasar in addition to a violent starburst. Multicolor image reconstruction of the SMG indicates a bipolar morphology of the emitted ultraviolet (UV) light suggestive of cones through which UV light escapes a dust-enshrouded nucleus.

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.

Obtaining lensing time delay measurements requires long-term monitoring campaigns with a high enough resolution (< 1 arcsec) to separate the multiple images. In the radio, a limited number of high-resolution interferometer arrays make these observations difficult to schedule. To overcome this problem, we propose a technique for measuring gravitational time delays which relies on monitoring the total flux density with low-resolution but high-sensitivity radio telescopes to follow the variation of the brighter image. This is then used to trigger high-resolution observations in optimal numbers which then reveal the variation in the fainter image. We present simulations to assess the efficiency of this method together with a pilot project observing radio lens systems with the Westerbork Synthesis Radio Telescope (WSRT) to trigger Very Large Array (VLA) observations. This new method is promising for measuring time delays because it uses relatively small amounts of time on high-resolution telescopes. This will be important because instruments that have high sensitivity but limited resolution, together with an optimum usage of followup high-resolution observations from appropriate radio telescopes may in the future be useful for gravitational lensing time delay measurements by means of this new method.

We present the average rest-frame spectrum of high-redshift dusty, star-forming galaxies from 250-770GHz. This spectrum was constructed by stacking ALMA 3mm spectra of 22 such sources discovered by the South Pole Telescope and spanning z=2.0-5.7. In addition to multiple bright spectral features of 12CO, [CI], and H2O, we also detect several faint transitions of 13CO, HCN, HNC, HCO+, and CN, and use the observed line strengths to characterize the typical properties of the interstellar medium of these high-redshift starburst galaxies. We find that the 13CO brightness in these objects is comparable to that of the only other z>2 star-forming galaxy in which 13CO has been observed. We show that the emission from the high-critical density molecules HCN, HNC, HCO+, and CN is consistent with a warm, dense medium with T_kin ~ 55K and n_H2 >~ 10^5.5 cm^-3. High molecular hydrogen densities are required to reproduce the observed line ratios, and we demonstrate that alternatives to purely collisional excitation are unlikely to be significant for the bulk of these systems. We quantify the average emission from several species with no individually detected transitions, and find emission from the hydride CH and the linear molecule CCH for the first time at high redshift, indicating that these molecules may be powerful probes of interstellar chemistry in high-redshift systems. These observations represent the first constraints on many molecular species with rest-frame transitions from 0.4-1.2mm in star-forming systems at high redshift, and will be invaluable in making effective use of ALMA in full science operations.

From a set of adaptive optics (AO) observations collected with the W.M. Keck telescope between August and September 2009, we derived the orbital parameters of the most recently discovered satellites of the large C-type asteroid (93) Minerva. The satellites of Minerva, which are approximately 3 and 4 km in diameter, orbit very close to the primary $\sim$5 & $\sim$8 $\times$ Rp and $\sim$1% & $\sim$2% $\times$ RHill) in a circular manner, sharing common characteristics with most of the triple asteroid systems in the main-belt. Combining these AO observations with lightcurve data collected since 1980 and two stellar occultations in 2010 & 2011, we removed the ambiguity of the pole solution of Minerva's primary and showed that it has an almost regular shape with an equivalent diameter Deq = 154 $\pm$ 6 km in agreement with IRAS observations. The surprisingly high bulk density of 1.75 $\pm$ 0.30 g/cm$\^3$ for this C-type asteroid, suggests that this taxonomic class is composed of asteroids with different compositions, For instance, Minerva could be made of the same material as dry CR, CO, and CV meteorites. We discuss possible scenarios on the origin of the system and conclude that future observations may shine light on the nature and composition of this fifth known triple main-belt asteroid.

We present observations of SPT-S J053816-5030.8, a gravitationally-lensed dusty star forming galaxy (DSFG) at z = 2.7817, first discovered at millimeter wavelengths by the South Pole Telescope. SPT 0538-50 is typical of the brightest sources found by wide-field millimeter-wavelength surveys, being lensed by an intervening galaxy at moderate redshift (in this instance, at z = 0.441). We present a wide array of multi-wavelength spectroscopic and photometric data on SPT 0538-50, including data from ALMA, Herschel PACS and SPIRE, Hubble, Spitzer, VLT, ATCA, APEX, and the SMA. We use high resolution imaging from HST to de-blend SPT 0538-50, separating DSFG emission from that of the foreground lens. Combined with a source model derived from ALMA imaging (which suggests a magnification factor of 21 +/- 4), we derive the intrinsic properties of SPT 0538-50, including the stellar mass, far-IR luminosity, star formation rate, molecular gas mass, and - using molecular line fluxes - the excitation conditions within the ISM. The derived physical properties argue that we are witnessing compact, merger-driven star formation in SPT 0538-50, similar to local starburst galaxies, and unlike that seen in some other DSFGs at this epoch.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

We report the detection of CO(1-0) line emission in the bright, lensed star-forming galaxies SPT-S 233227-5358.5 (z=2.73) and SPT-S 053816-5030.8 (z=2.78), using the Australia Telescope Compact Array (ATCA). Both galaxies were discovered in a large-area millimeter survey with the South Pole Telescope (SPT) and found to be gravitationally lensed by intervening structures. The measured CO intensities imply galaxies with molecular gas masses of (3.2 \pm 0.5)x10^10 (mu/15)^{-1}(X_CO/0.8) M_sun and (1.7 \pm 0.3)x10^10 (mu/20)^{-1}(X_CO/0.8) M_sun, and gas depletion timescales of 4.9x10^7 (X_CO/0.8) yr and 2.6x10^7 (X_CO/0.8) yr, respectively, where mu corresponds to the lens magnification and X_CO is the CO luminosity to gas mass conversion factor. In the case of SPT-S 053816-5030.8, we also obtained significant detections of the rest-frame 115.7 and 132.4 GHz radio continuum. Based on the radio to infrared spectral energy distribution and an assumed synchrotron spectral index, we find that 42 \pm 10 % and 55 \pm 13 % of the flux at rest-frame 115.7 and 132.4 GHz arise from free-free emission. We find a radio-derived intrinsic star formation rate (SFR) of 470 \pm 170 M_sun/yr, consistent within the uncertainties with the infrared estimate. Based on the morphology of this object in the source plane, the derived gas mass and the possible flattening of the radio spectral index towards low frequencies, we argue that SPT-S 053816-5030.8 exhibits properties compatible with a scaled-up local ultra-luminous infrared galaxy.

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.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 860 micrometer imaging of four high-redshift (z=2.8-5.7) dusty sources that were detected using the South Pole Telescope (SPT) at 1.4 mm and are not seen in existing radio to far-infrared catalogs. At 1.5 arcsec resolution, the ALMA data reveal multiple images of each submillimeter source, separated by 1-3 arcsec, consistent with strong lensing by intervening galaxies visible in near-IR imaging of these sources. We describe a gravitational lens modeling procedure that operates on the measured visibilities and incorporates self-calibration-like antenna phase corrections as part of the model optimization, which we use to interpret the source structure. Lens models indicate that SPT0346-52, located at z=5.7, is one of the most luminous and intensely star-forming sources in the universe with a lensing corrected FIR luminosity of 3.7 X 10^13 L_sun and star formation surface density of 4200 M_sun yr^-1 kpc^-2. We find magnification factors of 5 to 22, with lens Einstein radii of 1.1-2.0 arcsec and Einstein enclosed masses of 1.6-7.2x10^11 M_sun. These observations confirm the lensing origin of these objects, allow us to measure the their intrinsic sizes and luminosities, and demonstrate the important role that ALMA will play in the interpretation of lensed submillimeter sources.

Using the Atacama Large Millimeter/submillimeter Array (ALMA), we have conducted a blind redshift survey in the 3 mm atmospheric transmission window for 26 strongly lensd dusty star-forming galaxies (DSFGs) selected with the South Pole Telescope (SPT). The sources were selected to have S_1.4mm>20 mJy and a dust-like spectrum and, to remove low-z sources, not have bright radio (S_843MHz<6mJy) or far-infrared counterparts (S_100um<1 Jy, S_60um<200mJy). We robustly detect 44 line features in our survey, which we identify as redshifted emission lines of 12CO, 13CO, [CI], H2O, and H2O+. We find one or more spectral features in 23 sources yielding a ~90% detection rate for this survey; in 12 of these sources we detect multiple lines, while in 11 sources we detect only a single line. For the sources with only one detected line, we break the redshift degeneracy with additional spectroscopic observations if available, or infer the most likely line identification based on photometric data. This yields secure redshifts for ~70% of the sample. The three sources with no lines detected are tentatively placed in the redshift desert between 1.7<z<2.0. The resulting mean redshift of our sample is <z>=3.5. This finding is in contrast to the redshift distribution of radio-identified DSFGs, which have a significantly lower mean redshift of <z>=2.3 and for which only 10-15% of the population is expected to be at z>3. We discuss the effect of gravitational lensing on the redshift distribution and compare our measured redshift distribution to that of models in the literature.

In the past decade, our understanding of galaxy evolution has been revolutionized by the discovery that luminous, dusty, starburst galaxies were 1,000 times more abundant in the early Universe than at present. It has, however, been difficult to measure the complete redshift 2 distribution of these objects, especially at the highest redshifts (z > 4). Here we report a redshift survey at a wavelength of three millimeters, targeting carbon monoxide line emission from the star-forming molecular gas in the direction of extraordinarily bright millimetrewave-selected sources. High-resolution imaging demonstrates that these sources are strongly gravitationally lensed by foreground galaxies. We detect spectral lines in 23 out of 26 sources and multiple lines in 12 of those 23 sources, from which we obtain robust, unambiguous redshifts. At least 10 of the sources are found to lie at z > 4, indicating that the fraction of dusty starburst galaxies at high redshifts is greater than previously thought. Models of lens geometries in the sample indicate that the background objects are ultra-luminous infrared galaxies, powered by extreme bursts of star formation.

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.

We use Chandra observations of nine optically and X-ray selected clusters in five different structures at z ~ 0.7-1.1 from the Observations of Redshift Evolution in Large-Scale Environments (ORELSE) survey to study diffuse X-ray emission from galaxy clusters. X-ray gas temperatures and bolometric rest-frame luminosities are measured for each cluster in the sample. We present new redshift measurements, derived from dataobtained using the Deep Imaging Multi-Object Spectrograph on the Keck 10-m telescope, for two clusters in the RX J0910 supercluster at z ~ 1.1, from which velocity dispersions are measured. Dispersions for all clusters are combined with X-ray luminosities and gas temperatures to evaluate how the cluster properties compare to low-redshift scaling relations. We also measure the degree of substructure in each cluster by examining the velocity histograms, performing Dressler-Shectman tests, and computing the offsets between the X-ray emission center and optically-derived centroids. We find that only two clusters show clear indications of being unrelaxed, based on their scaling relations and other dynamical state diagnostics. Using our sample, we evaluate the redshift evolution of the L_x-T relation and investigate the implications of our results for precision cosmology surveys.

We present APEX SABOCA 350micron and LABOCA 870micron observations of 11 representative examples of the rare, extremely bright (S_1.4mm > 15mJy), dust-dominated millimeter-selected galaxies recently discovered by the South Pole Telescope (SPT). All 11 sources are robustly detected with LABOCA with 40 < S_870micron < 130mJy, approximately an order of magnitude higher than the canonical submillimeter galaxy (SMG) population. Six of the sources are also detected by SABOCA at >3sigma, with the detections or upper limits providing a key constraint on the shape of the spectral energy distribution (SED) near its peak. We model the SEDs of these galaxies using a simple modified blackbody and perform the same analysis on samples of SMGs of known redshift from the literature. These calibration samples inform the distribution of dust temperature for similar SMG populations, and this dust temperature prior allows us to derive photometric redshift estimates and far infrared luminosities for the sources. We find a median redshift of <z> = 3.0, higher than the <z> = 2.2 inferred for the normal SMG population. We also derive the apparent size of the sources from the temperature and apparent luminosity, finding them to appear larger than our unlensed calibration sample, which supports the idea that these sources are gravitationally magnified by massive structures along the line of sight.

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

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 investigate the properties of the 525 spectroscopically confirmed members of the Cl1604 supercluster at z~0.9 as part of the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey. Using extensive Keck LRIS/DEIMOS spectroscopy in conjunction with ten-band ground-based, Spitzer, and HST imaging, we investigate the buildup of the red sequence in groups and clusters at high redshift. Nearly all of the brightest and most massive red-sequence galaxies are found within the bounds of the clusters and groups. Despite the prevalence of these red-sequence galaxies, we find that the average cluster galaxy has a spectrum indicative of a star-forming galaxy, with a star formation rate between those of z~1 field galaxies and moderate redshift cluster galaxies. The average group galaxy is even more active, exhibiting properties indicative of a starburst. The presence of massive, red galaxies and the high fraction of starbursting galaxies suggest that significant processing is occurring in the group environment at z~1 and earlier. There is a deficit of low-luminosity red-sequence galaxies in all Cl1604 clusters and groups, suggesting that such galaxies transition to the red sequence at later times. Extremely massive (10^12) red sequence galaxies are also absent from the Cl1604 clusters and groups. We suggest that such galaxies form at later times through merging processes. There are also large populations of transition galaxies at intermediate stellar masses present in the groups and clusters, suggesting that such masses are important in the buildup of the red-sequence mass function at z~1. Through a comparison of the transitional populations present in the Cl1604 clusters and groups, we find evidence that massive blue cloud galaxies are quenched earliest in the most dynamically relaxed systems and at progressively later times in dynamically unrelaxed systems.

We use deep Chandra imaging and an extensive optical spectroscopy campaign on the Keck 10-m telescopes to study the properties of X-ray point sources in five large-scale structures at redshifts of z ~ 0.7-0.9. We first study X-ray point sources using the statistical measure of cumulative source counts, finding that the measured overdensities are consistent with previous results, but we recommend caution in overestimating the precision of the technique. Optical spectroscopy of objects matched to X-ray point sources confirms a total of 27 AGN within the five structures, and we find that their host galaxies tend to be located away from dense cluster cores. More than 36% of host galaxies are located in the `green valley', which suggests they are a transitional population. Based on analysis of OII and Hd line strengths, the average spectral properties of the AGN host galaxies in all structures indicate either on-going star formation or a starburst within ~ 1 Gyr, and the host galaxies are younger than the average galaxy in the parent population. These results indicate a clear connection between starburst and nuclear activity. We use composite spectra of the spectroscopically confirmed members in each structure to separate them based on a measure of the overall evolutionary state of their constituent galaxies. We define structures as having more evolved populations if their average galaxy has lower EW(OII) and EW(Hd). The AGN in the more evolved structures have lower rest-frame 0.5-8 keV X-ray luminosities (all below 10^43.3 erg s^-1) and longer times since a starburst than those in the less evolved structures, suggesting that the peak of both star formation and AGN activity has occurred at earlier times. With the wide range of evolutionary states and timeframes in the structures, we use our results to analyze the evolution of X-ray AGN and evaluate potential triggering mechanisms.

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.

We present the discovery and characterisation of PTF10iya, a short-lived (dt ~ 10 d, with an optical decay rate of ~ 0.3 mag per d), luminous (M_g ~ -21 mag) transient source found by the Palomar Transient Factory. The ultraviolet/optical spectral energy distribution is reasonably well fit by a blackbody with T ~ 1-2 x 10^4 K and peak bolometric luminosity L_BB ~ 1-5 x 10^44 erg per s (depending on the details of the extinction correction). A comparable amount of energy is radiated in the X-ray band that appears to result from a distinct physical process. The location of PTF10iya is consistent with the nucleus of a star-forming galaxy (z = 0.22405 +/- 0.00006) to within 350 mas (99.7 per cent confidence radius), or a projected distance of less than 1.2 kpc. At first glance, these properties appear reminiscent of the characteristic "big blue bump" seen in the near-ultraviolet spectra of many active galactic nuclei (AGNs). However, emission-line diagnostics of the host galaxy, along with a historical light curve extending back to 2007, show no evidence for AGN-like activity. We therefore consider whether the tidal disruption of a star by an otherwise quiescent supermassive black hole may account for our observations. Though with limited temporal information, PTF10iya appears broadly consistent with the predictions for the early "super-Eddington" phase of a solar-type star disrupted by a ~ 10^7 M_sun black hole. Regardless of the precise physical origin of the accreting material, the large luminosity and short duration suggest that otherwise quiescent galaxies can transition extremely rapidly to radiate near the Eddington limit; many such outbursts may have been missed by previous surveys lacking sufficient cadence.

The X-ray source CXOXBJ142607.6+353351 (CXOJ1426+35), which was identified in a 172 ks Chandra image in the Bootes field, shows double-peaked rest-frame optical/UV emission lines, separated by 0.69" (5.5 kpc) in the spatial dimension and by 690 km s^-1 in the velocity dimension. The high excitation lines and emission line ratios indicate both systems are ionized by an AGN continuum, and the double-peaked profile resembles that of candidate dual AGN. At a redshift of z=1.175, this source is the highest redshift candidate dual AGN yet identified. However, many sources have similar emission line profiles for which other interpretations are favored. We have analyzed the substantial archival data available in this field, as well as acquired near-infrared (NIR) adaptive optics (AO) imaging and NIR slit spectroscopy. The X-ray spectrum is hard, implying a column density of several 10^23 cm^-2. Though heavily obscured, the source is also one of the brightest in the field, with an absorption-corrected 2-10 keV luminosity of ~10^45 erg s^-1. Outflows driven by an accretion disk may produce the double-peaked lines if the central engine accretes near the Eddington limit. However, we may be seeing the narrow line regions of two AGN following a galactic merger. While the AO image reveals only a single source, a second AGN would easily be obscured by the significant extinction inferred from the X-ray data. Understanding the physical processes producing the complex emission line profiles seen in CXOJ1426+35 and related sources is important for interpreting the growing population of dual AGN candidates.

We compare galaxy number counts in HST/ACS fields containing moderate-redshift (0.2<z<1.0) strong gravitational lenses with those in two control samples: (1) the first square degree of the COSMOS survey, comprising 259 ACS fields and (2) 20 "pure parallel" fields randomly located on the sky. Through a Bayesian analysis we determine the expectation values (mu_0) and confidence levels of the underlying number counts for a range of apertures and magnitude bins. Our analysis has produced the following results: (i) We infer that our control samples are not consistent, with the number counts in the COSMOS sample being significantly higher than in the pure parallel sample for 21 <= F814W <= 23. This result matches those found in previous analyses of COSMOS data using different techniques. (ii) We find that small-size apertures, centered on strong lenses, are overdense compared with randomly placed apertures in the control samples, even compared to the COSMOS sample. Correcting for the local clustering of elliptical galaxies, based on the average two-point correlation function reduces this overdensity to the 1-2 sigma level. Thus, the overdensity of galaxies seen along a typical line of sight to a lens can be explained mostly by the natural clustering of galaxies, rather than being due to lenses lying along otherwise biased lines of sight. However, a larger sample of lenses is needed to see if the remaining bias persists when the lens-field uncertainties are smaller. (iii) There is considerable scatter in the lines of sight to _individual_ lens systems, but quantities that are linearly dependent on the external convergence (e.g., H_0) should become unbiased if the extra galaxies that cause the bias can be accounted for in the lens models either through direct modeling or via an informed prior on the external convergence. The number counts can used to set such an informed prior.

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.

The Hubble constant value is currently known to 10% accuracy unless assumptions are made for the cosmology (Sandage et al. 2006). Gravitational lens systems provide another probe of the Hubble constant using time delay measurements. However, current investigations of ~20 time delay lenses, albeit of varying levels of sophistication, have resulted in different values of the Hubble constant ranging from 50-80 km/s/Mpc. In order to reduce uncertainties, more time delay measurements are essential together with better determined mass models (Oguri 2007, Saha et al. 2006). We propose a more efficient technique for measuring time delays which does not require regular monitoring with a high-resolution interferometer array. The method uses double image and long-axis quadruple lens systems in which the brighter component varies first and dominates the total flux density. Monitoring the total flux density with low-resolution but high sensitivity radio telescopes provides the variation of the brighter image and is used to trigger high-resolution observations which can then be used to see the variation in the fainter image. We present simulations of this method together with a pilot project using the WSRT (Westerbork Radio Synthesis Telescope) to trigger VLA (Very Large Array) observations. This new method is promising for measuring time delays because it uses relatively small amounts of time on high-resolution telescopes. This will be important because many SKA pathfinder telescopes, such as MeerKAT (Karoo Array Telescope) and ASKAP (Australian Square Kilometre Array Pathfinder), have high sensitivity but limited resolution.

Using Spitzer-MIPS 24um imaging and Keck spectroscopy we examine the nature of the obscured star forming population in three clusters and three groups at z~0.9. These six systems are components of the Cl1604 supercluster, the largest structure imaged by Spitzer at redshifts near unity. We find that the average density of 24um-detected galaxies within the Cl1604 clusters is nearly twice that of the surrounding field and that this overdensity scales with the cluster's dynamical state. The 24um-bright members often appear optically unremarkable and exhibit only moderate [OII] line emission due to severe obscuration. Their spatial distribution suggests they are an infalling population, but an examination of their spectral properties, morphologies and optical colors indicate they are not simply analogs of the field population that have yet to be quenched. Using stacked composite spectra, we find the 24um-detected cluster and group galaxies exhibit elevated levels of Balmer absorption compared to galaxies undergoing normal, continuous star formation. A similar excess is not observed in field galaxies with equivalent infrared luminosities, indicating a greater fraction of the detected cluster and group members have experienced a burst of star formation in the recent past compared to their counterparts in the field. Our results suggest that gas-rich galaxies at high redshift experience a temporary increase in their star formation activity as they assemble into denser environments. Using HST-ACS imaging we find that disturbed morphologies are common among the 24um-detected cluster and group members and become more prevalent in regions of higher galaxy density. We conclude that mergers are the dominant triggering mechanism responsible for the enhanced star formation found in the Cl1604 groups, while a mix of harassment and mergers are likely driving the activity of the cluster galaxies.

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 present a tomographic cosmological weak lensing analysis of the HST COSMOS Survey. Applying our lensing-optimized data reduction, principal component interpolation for the ACS PSF, and improved modelling of charge-transfer inefficiency, we measure a lensing signal which is consistent with pure gravitational modes and no significant shape systematics. We carefully estimate the statistical uncertainty from simulated COSMOS-like fields obtained from ray-tracing through the Millennium Simulation. We test our pipeline on simulated space-based data, recalibrate non-linear power spectrum corrections using the ray-tracing, employ photometric redshifts to reduce potential contamination by intrinsic galaxy alignments, and marginalize over systematic uncertainties. We find that the lensing signal scales with redshift as expected from General Relativity for a concordance LCDM cosmology, including the full cross-correlations between different redshift bins. For a flat LCDM cosmology, we measure sigma_8(Omega_m/0.3)^0.51=0.75+-0.08 from lensing, in perfect agreement with WMAP-5, yielding joint constraints Omega_m=0.266+0.025-0.023, sigma_8=0.802+0.028-0.029 (all 68% conf.). Dropping the assumption of flatness and using HST Key Project and BBN priors only, we find a negative deceleration parameter q_0 at 94.3% conf. from the tomographic lensing analysis, providing independent evidence for the accelerated expansion of the Universe. For a flat wCDM cosmology and prior w in [-2,0], we obtain w<-0.41 (90% conf.). Our dark energy constraints are still relatively weak solely due to the limited area of COSMOS. However, they provide an important demonstration for the usefulness of tomographic weak lensing measurements from space. (abridged)

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)

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We report on the results of a Chandra search for evidence of triggered nuclear activity within the Cl0023+0423 four-way group merger at z ~ 0.84. The system consists of four interacting galaxy groups in the early stages of hierarchical cluster formation and, as such, provides a unique look at the level of processing and evolution already under way in the group environment prior to cluster assembly. We present the number counts of X-ray point sources detected in a field covering the entire Cl0023 structure, as well as a cross-correlation of these sources with our extensive spectroscopic database. Both the redshift distribution and cumulative number counts of X-ray sources reveal little evidence to suggest that the system contains X-ray luminous active galactic nuclei (AGNs) in excess to what is observed in the field population. If preprocessing is under way in the Cl0023 system, our observations suggest that powerful nuclear activity is not the predominant mechanism quenching star formation and driving the evolution of Cl0023 galaxies. We speculate that this is due to a lack of sufficiently massive nuclear black holes required to power such activity, as previous observations have found a high late-type fraction among the Cl0023 population. It may be that disruptive AGN-driven outflows become an important factor in the preprocessing of galaxy populations only during a later stage in the evolution of such groups and structures when sufficiently massive galaxies (and central black holes) have built up, but prior to hydrodynamical processes stripping them of their gas reservoirs.

Compact Symmetric Objects (CSOs) are small (less than 1 kpc) radio sources which have symmetric double lobes or jets. The dominant theory for the small size of these objects is that they are young radio sources which could grow into larger radio galaxies, but the currently small population of known CSOs makes it difficult to definitively determine whether or not this is the case. While a greater number of Gigahertz peaked sources can be identified by sifting through spectral surveys, this yields none of the dynamics of the sources, and also brings Quasars into the sample, which although interesting are peaked around 1 Gigahertz for very different reasons. We have used the 5 GHz VLBA Imaging and Polarization Survey (VIPS) to identify 103 CSO candidates morphologically, and are following up on these sources with multifrequency VLBA observations to confirm CSO identifications and to study their dynamics. The identification of candidates from within the survey will be discussed, as well as preliminary results from the follow-up observations.

We present results of a spectroscopic search for Lyman-alpha emitters (LAEs) in the Cl1604 supercluster field using the extensive spectroscopic Keck/DEIMOS database taken as part of the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey. A total of 12 slitmasks were observed and inspected in the Cl1604 field, spanning a survey volume of 1.365x10^4 co-moving Mpc^3. We find a total of 17 high redshift (4.39 < z < 5.67) LAE candidates down to a limiting flux of 1.9x10^(-18) ergs/s/cm (~0.1L* at z~5), 13 of which we classify as high quality. The resulting LAE number density is nearly double that of LAEs found in the Subaru deep field at z~4.9 and nearly an order of magnitude higher than in other surveys of LAEs at similar redshifts, an excess that is essentially independent of LAE luminosity. We also report on the discovery of two possible LAE group structures at z~4.4 and z~4.8 and investigate the effects of cosmic variance of LAEs on our results. Fitting a simple truncated single Gaussian model to a composite spectrum of the 13 high quality LAE candidates, we find a best-fit stellar velocity dispersion of 136 km/s. Additionally, we see modest evidence of a second peak in the composite spectrum, possibly caused by galactic outflows, as well as evidence for a non-trivial Lyman-alpha escape fraction. We find an average LAE star formation rate density (SFRD) of ~5x10^(-3) M_solar/yr/Mpc^3 with moderate evidence for negative evolution in the LAE SFRD from z~4.6 to z~5.7. We measure a best-fit luminosity function generally consistent with measurements from other surveys at similar epochs. Finally, we investigate any possible effects from weak or strong gravitational lensing induced by the foreground supercluster, finding that our LAE candidates are minimally affected by lensing processes.

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

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.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

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.

Strong gravitational lensing is a powerful technique for probing galaxy mass distributions and for measuring cosmological parameters. We present a pixelated approach to modeling simultaneously the lens potential and source intensity of strong gravitational lens systems with extended source-intensity distributions. For systems with sources of sufficient extent such that the separate lensed images are connected by intensity measurements, the accuracy in the reconstructed potential is solely limited by the quality of the data. We apply this potential reconstruction technique to deep HST observations of B1608+656, a four-image gravitational lens system formed by a pair of interacting lens galaxies. We present a comprehensive Bayesian analysis of the system that takes into account the extended source-intensity distribution, dust extinction, and the interacting lens galaxies. Our approach allows us to compare various models of the components of the lens system, which include the point-spread function (PSF), dust, lens galaxy light, source-intensity distribution, and lens potential. Using optimal combinations of the PSF, dust, and lens galaxy light models, we successfully reconstruct both the lens potential and the extended source-intensity distribution of B1608+656. The resulting reconstruction can be used as the basis of a measurement of the Hubble constant. We use our reconstruction of the gravitational potential to study the relative distribution of mass and light in the lensing galaxies. We find that the mass-to-light ratio for the primary lens galaxy is (2.0+/-0.2)h M_{\sun} L_{B,\sun}^{-1} within the Einstein radius 3.9 h^{-1} kpc, in agreement with what is found for noninteracting lens galaxies at the same scales. (Abridged)

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.

To investigate the role of feedback from Active Galactic Nuclei (AGN) in driving the evolution of their host galaxies, we have carried out a study of the environments and optical properties of galaxies harboring X-ray luminous AGN in the Cl1604 supercluster at z~0.9. Making use of Chandra, HST/ACS and Keck/DEIMOS observations, we examine the integrated colors, morphologies and spectral properties of nine moderate-luminosity (L_x ~ 10^43 erg s^-1) type 2 Seyferts detected in the Cl1604 complex. We find that the AGN are predominantly hosted by luminous spheroids and/or bulge dominated galaxies which have colors that place them in the valley between the blue cloud and red sequence in color-magnitude space, consistent with predictions that AGN hosts should constitute a transition population. Half of the hosts have bluer overall colors as a result of blue resolved cores in otherwise red spheroids and a majority show signs of recent or pending interactions. We also find a substantial number exhibit strong Balmer absorption features indicative of post-starburst galaxies, despite the fact that we detect narrow [OII] emission lines in all of the host spectra. If the [OII] lines are due in part to AGN emission, as we suspect, then this result implies that a significant fraction of these galaxies (44%) have experienced an enhanced level of star formation within the last ~1 Gyr which was rapidly suppressed. Overall we find that the properties of the nine host galaxies are generally consistent with a scenario in which recent interactions have triggered both increased levels of nuclear activity and an enhancement of centrally concentrated star formation, followed by a rapid truncation of the latter, possibly as a result of feedback from the AGN itself. [Abridged]

We present multi-frequency multi-epoch Very Long Baseline Array (VLBA) observations of J11584+2450. These observations clearly show this source, previously classified as a core-jet, to be a compact symmetric object (CSO). Comparisons between these new data and data taken over the last 9 years shows the edge brightened hot spots retreating towards the core (and slightly to the west) at approximately 0.3c. Whether this motion is strictly apparent or actually physical in nature is discussed, as well as possible explanations, and what implications a physical contraction of J11584+2450 would have for current CSO models.

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.

(abridged) The Observatory for Multi-Epoch Gravitational Lens Astrophysics (OMEGA) is a mission concept for a 1.5-m near-UV through near-IR space observatory that will be dedicated to frequent imaging and spectroscopic monitoring of ~100 multiply-imaged active galactic nuclei over the whole sky. Using wavelength-tailored dichroics with extremely high transmittance, efficient imaging in six channels will be done simultaneously during each visit to each target. The separate spectroscopic mode, engaged through a flip-in mirror, uses an image slicer spectrograph. After a period of many visits to all targets, the resulting multidimensional movies can then be analyzed to a) measure the mass function of dark matter substructure; b) measure precise masses of the accreting black holes as well as the structure of their accretion disks and their environments over several decades of physical scale; and c) measure a combination of the local Hubble expansion and cosmological distances to unprecedented precision.

Lens modeling is the key to successful and meaningful automated strong galaxy-scale gravitational lens detection. We have implemented a lens-modeling "robot" that treats every bright red galaxy (BRG) in a large imaging survey as a potential gravitational lens system. Using a simple model optimized for "typical" galaxy-scale lenses, we generate four assessments of model quality that are used in an automated classification. The robot infers the lens classification parameter H that a human would have assigned; the inference is performed using a probability distribution generated from a human-classified training set, including realistic simulated lenses and known false positives drawn from the HST/EGS survey. We compute the expected purity, completeness and rejection rate, and find that these can be optimized for a particular application by changing the prior probability distribution for H, equivalent to defining the robot's "character." Adopting a realistic prior based on the known abundance of lenses, we find that a lens sample may be generated that is ~100% pure, but only ~20% complete. This shortfall is due primarily to the over-simplicity of the lens model. With a more optimistic robot, ~90% completeness can be achieved while rejecting ~90% of the candidate objects. The remaining candidates must be classified by human inspectors. We are able to classify lens candidates by eye at a rate of a few seconds per system, suggesting that a future 1000 square degree imaging survey containing 10^7 BRGs, and some 10^4 lenses, could be successfully, and reproducibly, searched in a modest amount of time. [Abridged]

We present the results of Chandra observations of the Cl1604 supercluster at z~0.9. The system is the largest structure mapped at redshifts approaching unity, containing at least eight spectroscopically confirmed galaxy clusters and groups. Using two 50-ksec ACIS-I pointings we examine both the X-ray point source population and the diffuse emission from individual clusters in the system. We find a 2.5\sigma excess of point sources detected in the hard band (2-10 keV) relative to the number of sources found in blank fields observed by Chandra. No such excess is observed in the soft band (0.5-2 keV). The hard-band source density is 1.47 times greater than that of a blank field, in agreement with the previously reported correlation between overdensity amplitude and cluster redshift. Using a maximum likelihood technique we have matched 112 of the 161 detected X-ray point sources to optical counterparts and found 15 sources that are associated with the supercluster. All 15 sources have rest-frame luminosities consistent with emission from active galactic nuclei (AGN). We find that the supercluster AGN largely avoid the densest regions of the system and are instead distributed on the outskirts of massive clusters or within poorer clusters and groups. We have also detected diffuse emission from two of the eight clusters and groups in the system, clusters Cl1604+4304 and Cl1604+4314. The systems have bolometric luminosities of 1.43x10^44 and 8.20x10^43 h70^-2 erg s^-1 and gas temperatures of 3.50 (+1.82-1.08) and 1.64 (+0.65-0.45) keV, respectively. Using updated velocity dispersions, we compare the properties of these systems to the cluster scaling relations followed by other X-ray and optically selected galaxy clusters at high redshift.

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.

The gravitational lens system CLASS B2108+213 has two lensed images separated by 4.56 arcsec. Such a wide image separation suggests that the lens is either a massive galaxy, or is composed of a group of galaxies. To investigate the structure of the lensing potential we have carried out new high resolution imaging of the two lensed images at 1.7 GHz with the VLBA and at 5 GHz with global VLBI. Compact and extended emission is detected from the two lensed images, which provides additional constraints to the lensing mass model. We find that the data are consistent with either a single lensing galaxy, or a two galaxy lens model that takes account of a nearby companion to the main lensing galaxy within the Einstein radius of the system. However, for an ensemble of global power-law mass models, those with density profiles steeper than isothermal are a better fit. The best-fitting profile for a single spherical mass model has a slope of $\gamma=$~2.45$_{-0.18}^{+0.19}$. The system also has a third radio component which is coincident with the main lensing galaxy. This component is detected at milli-arcsecond scales for the first time by the 1.7 GHz VLBA and 5 GHz global VLBI imaging. However, the third radio component is found not to be consistent with a core lensed image because the radio spectrum differs from the two lensed images, and its flux-density is too high when compared to what is expected from simple mass models with a variable power-law density profile and/or a reasonable core radius. Furthermore, 1.4 GHz imaging of the system with the MERLIN finds extended lobe emission on either side of the main lensing galaxy. Therefore, the radio emission from the third radio component is almost certainly from an AGN within the main lensing galaxy, which is classified as an FR I type radio source.

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 examine the radio properties of EGRET-detected blazars observed as part of the VLBA Imaging and Polarimetry Survey (VIPS). VIPS has a flux limit roughly an order of magnitude below the MOJAVE survey and most other samples that have been used to study the properties of EGRET blazars. At lower flux levels, radio flux density does not directly correlate with gamma-ray flux density. We do find that the EGRET-detected blazars tend to have higher brightness temperatures, greater core fractions, and possibly larger than average jet opening angles. A weak correlation is also found with jet length and with polarization. All of the well-established trends can be explained by systematically larger Doppler factors in the gamma-ray loud blazars, consistent with the measurements of higher apparent velocities found in monitoring programs carried out at radio frequencies above 10 GHz.

We discuss gravitational lensing by elliptical galaxies with some particular mass distributions. Using simple techniques from the theory of quadrature domains and the Schwarz function (cf. \cite{Sh}) we show that when the mass density is constant on confocal ellipses, the total number of lensed images of a point source cannot exceed 5 (4 bright images and 1 dim image). Also, using the Dive--Nikliborc converse of the celebrated Newton's theorem concerning the potentials of ellipsoids, we show that ``Einstein rings'' must always be either circles (in the absence of a tidal shear), or ellipses.

We investigate the galaxy populations in seven X-ray selected, intermediate-redshift groups (0.2<z<0.6). Overall, the galaxy populations in these systems are similar to those in clusters at the same redshift; they have large fractions of early-type galaxies (f_e~70%) and small fractions of galaxies with significant star formation (f_[OII]~30%). We do not observe a strong evolution in the galaxy populations from those seen in X-ray luminous groups at low-redshift. Both f_e and f_[OII] are correlated with radius but do not reach the field value out to ~r_500. However, we find significant variation in the galaxy populations between groups with some groups having field-like populations. Comparisons between the morphological and spectral properties of group galaxies reveal both gas-poor mergers and a population of passive spirals. Unlike low-redshift, X-ray emitting groups, in some of these groups the brightest galaxy does not lie at the center of the X-ray emission, and in several of the groups that do have a central BGG, the BGG has multiple components. These groups appear to represent a range of evolutionary stages in the formation of the BGG. Some groups have relatively large central galaxy densities, and one group contains a string of seven bright galaxies within a radius of 200 kpc that have a lower velocity dispersion than the rest of the system. None of the central galaxies, including those with multiple components, have significant [OII] emission. These observations support a scenario in which BGGs are formed relatively late through gas-poor mergers.

The existence of flux-ratio anomalies between fold and cusp images in galaxy-scale strong-lens systems has led to an interpretation based on the presence of a high mass-fraction of cold-dark-matter (CDM) substructures around galaxies, as predicted by numerical N-body simulations. The flux-ratio anomaly is particularly evident in the radio-loud quadruple gravitational lens system CLASS B2045+265. In this paper, new high-resolution radio, optical, and infrared imaging of B2045+265 is presented which sheds more light on this anomaly and its possible causes. First, deep Very Long Baseline Array observations show very compact images, possibly with a hint of a jet, but with no evidence for differential scattering or scatter broadening. Second, optical and infrared observations with the Hubble Space Telescope and through Adaptive-Optics imaging with the W. M. Keck Telescope, show a previously undiscovered object -- interpreted as a (tidally disrupted) dwarf satellite based on its colours and slight extension -- between the main lens galaxy and the three anomalous flux-ratio images. Third, colour variations in the early-type lens galaxy indicate recent star-formation, possibly the result of secondary infall of gas-rich satellites. A population of such galaxies around the lens system could explain the previously discovered strong [O II] emission. However, spiral structure and/or normal star formation in the lens galaxy cannot be excluded. In light of these new data, we propose a lens model for the system, including the observed dwarf satellite, which reproduces all positional and flux-ratio constraints, without the need for additional CDM substructure. [abridged]

We present the selection of the Jodrell Bank Flat-spectrum (JBF) radio source sample, which is designed to reduce the uncertainties in the Cosmic Lens All-Sky Survey (CLASS) gravitational lensing statistics arising from the lack of knowledge about the parent population luminosity function. From observations at 4.86 GHz with the Very Large Array, we have selected a sample of 117 flat-spectrum radio sources with flux densities greater than 5 mJy. These sources were selected in a similar manner to the CLASS complete sample and are therefore representative of the parent population at low flux densities. The vast majority (~90 per cent) of the JBF sample are found to be compact on the arcsecond scales probed here and show little evidence of any extended radio jet emission. Using the JBF and CLASS complete samples we find the differential number counts slope of the parent population above and below the CLASS 30 mJy flux density limit to be -2.07+/-0.02 and -1.96+/-0.12, respectively.

We present the first results of the VLBA Imaging and Polarimetry Survey (VIPS), a 5 GHz VLBI survey of 1,127 sources with flat radio spectra. Through automated data reduction and imaging routines, we have produced publicly available I, Q, and U images and have detected polarized flux density from 37% of the sources. We have also developed an algorithm to use each source's I image to automatically classify it as a point-like source, a core-jet, a compact symmetric object (CSO) candidate, or a complex source. The mean ratio of the polarized to total 5 GHz flux density for VIPS sources with detected polarized flux density ranges from 1% to 20% with a median value of about 5%. We have also found significant evidence that the directions of the jets in core-jet systems tend to be perpendicular to the electric vector position angles (EVPAs). The data is consistent with a scenario in which ~24% of the polarized core-jets have EVPAs that are anti-aligned with the directions of their jet components and which have a substantial amount of Faraday rotation. In addition to these initial results, plans for future follow-up observations are discussed.

We report on the results of a visual search for galaxy-scale strong gravitational lenses over 650 arcmin^2 of HST/ACS imaging in the DEEP2-EGS field. In addition to a previously-known Einstein Cross (the "Cross," HST J141735+52264, with z_lens=0.8106 and a published z_source=3.40), we identify two new strong galaxy-galaxy lenses with multiple extended arcs. The first, HST J141820+52361 (the ``Dewdrop''; z_lens=0.5798, lenses two distinct extended sources into two pairs of arcs z_source=0.while), 9818 the second, HST J141833+52435 (the ``Anchor''; z_lens=0.4625), produces a single pair of arcs (source redshift not yet known). All three definite lenses are fit well by simple singular isothermal ellipsoid models including external shear. Using the three-dimensional line-of-sight (LOS) information on galaxies from the DEEP2 data, we calculate the convergence and shear contributions, assuming singular isothermal sphere halos truncated at 200 h^-1 kpc. These are also compared against three-dimensional local-density estimates. We find that even strong lenses in demonstrably underdense local environments may be considerably affected by LOS contributions, which in turn, may be underestimates of the effect of large scale structure.

We present spectroscopic confirmation of nine moderate redshift galaxy groups and poor clusters selected from the ROSAT Deep Cluster Survey. The groups span the redshift range z=0.23-0.59 and have between 4 and 20 confirmed members. The velocity dispersions of these groups range from 125 to 650 km/s. Similar to X-ray groups at low redshift, these systems contain a significant number of early-type galaxies. Therefore, the trend for X-ray luminous groups to have high early-type fractions is already in place by at least z=0.5. In 4 of the 9 groups, the X-ray emission is clearly peaked on the most luminous early-type galaxy in the group. However, in several cases the central galaxy is composed of multiple luminous nuclei, suggesting that the brightest group galaxy may still be undergoing major mergers. In at least 3 of the groups in our sample, a dominant early-type galaxy is not found at the center of the group potential. This suggests that many of our groups are not dynamically evolved despite their high X-ray luminosities. While similar systems have been identified at low redshift, the X-ray luminosities of the intermediate redshift examples are one to three orders of magnitude higher than those of their low redshift counterparts. We suggest that this may be evidence for group downsizing: while massive groups are still in the process of collapsing and virializing at intermediate redshifts, only low-mass groups are in the process of forming at the present day.

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.

We report the discovery of four groups of galaxies along the line of sight to the B1608+656 gravitational lens system. One group is at the redshift of the primary lensing galaxy (z = 0.631) and appears to have a low mass, with eight spectroscopically-confirmed members and an estimated velocity dispersion of 150+/-60 km/s. The three other groups are in the foreground of the lens. These groups contain ~10 confirmed members each, and are located at redshifts of 0.265, 0.426, and 0.520. Two of the three additional groups are centered roughly on the lens system, while the third is centered ~1 arcminute south of the lens. We find that each group provides an external convergence (kappa_ext) of ~0.005--0.060, depending on the assumptions made in calculating the convergence. For lens systems where no additional observables that can break the mass-sheet degeneracy exist, the determination of H_0 will be biased high by a factor of (1-kappa_ext)^{-1} if the external convergence is not properly included in the model. For the B1608+656 system, the stellar velocity dispersion of the lensing galaxy has been measured, thus breaking the mass-sheet degeneracy due to the group that is physically associated with the lens. The effect of the other groups along the line of sight can be folded into the overall uncertainties due to large-scale structure (LSS) along the line of sight. Because the B1608+656 system appears to lie along an overdense line of sight, the LSS will cause the measurement of H_0 to be biased high for this system. The systematic bias introduced by LSS could be 5% or greater. (Abridged)

As a step toward investigating the parsec-scale properties of faint extragalactic radio sources, the Very Long Baseline Array (VLBA) was used at 5.0 GHz to obtain phase-referenced images of 76 sources in the NOAO Bootes field. These 76 sources were selected from the FIRST catalog to have peak flux densities above 10 mJy at 5 arcsec resolution and deconvolved major diameters of less than 3 arcsec at 1.4 GHz. Fifty-seven of these faint radio sources were identified with accretion-powered radio galaxies and quasars brighter than 25.5 mag in the optical I band. On VLA scales at 1.4 GHz, a measure of the compactness of the faint sources (the ratio of the peak flux density from FIRST to the integrated flux density from the NVSS catalog) spans the full range of possibilites arising from source-resolution effects. Thirty of the faint radio sources, or 39 +9/-7%, were detected with the VLBA at 5.0 GHz with peak flux densities above 6 sigma ~2 mJy at 2 mas resolution. The VLBA detections occur through the full range of compactness ratios. The stronger VLBA detections can themselves serve as phase-reference calibrators, boding well for opening up much of the radio sky to VLBA imaging. For the adopted cosmology, the VLBA resolution correponds to 17 pc or finer. Most VLBA detections are unresolved or slightly resolved but one is diffuse and five show either double or core-jet structures; the properties of these latter six are discussed in detail. Three VLBA detections are unidentified and fainter than 25.5 mag in the optical I band; their properties are highlighted because they likely mark optically-obscured active nuclei at high redshift.

We report the discovery of a new gravitational lens system from the CLASS survey, CLASS B0631+519. VLA, MERLIN and VLBA observations show a doubly-imaged radio core, a doubly-imaged lobe and a second lobe that is probably quadruply-imaged. The maximum image separation is 1.16 arcseconds. The VLBA resolves the most magnified image of the flat-spectrum radio core into a number of sub-components spread across approximately 20 milli-arcseconds. Optical and near-infrared imaging with the ACS and NICMOS cameras on the HST show that there are two galaxies along the line of sight to the lensed source, as was previously discovered by optical spectroscopy. The nearer galaxy at z=0.0896 is a small blue irregular, while the more distant galaxy at z=0.6196 is an elliptical type and appears to contribute most of the lensing effect. The host galaxy of the lensed source is visible in NICMOS imaging as a set of arcs that form an almost complete Einstein ring. Mass modelling using non-parametric techniques can reproduce the ring and indicates that the irregular galaxy has a (localised) effect on the flux density distribution in the Einstein ring at the 5-10% level.

We present full polarization Very Long Baseline Array (VLBA) observations at 5 GHz and 15 GHz of 24 compact active galactic nuclei (AGN). These sources were observed as part of a pilot project to demonstrate the feasibility of conducting a large VLBI survey to further our understanding of the physical properties and temporal evolution of AGN jets. The sample is drawn from the Cosmic Lens All-Sky Survey (CLASS) where it overlaps with the Sloan Digital Sky Survey at declinations north of 15 degrees. There are 2100 CLASS sources brighter than 50 mJy at 8.4 GHz, of which we have chosen 24 for this pilot study. All 24 sources were detected and imaged at 5 GHz with a typical dynamic range of 500:1, and 21 of 24 sources were detected and imaged at 15 GHz. Linear polarization was detected in 8 sources at both 5 and 15 GHz, allowing for the creation of Faraday rotation measure (RM) images. The core RMs for the sample were found to have an average absolute value of 390 +/- 100 rad/m^2. We also present the discovery of a new Compact Symmetric Object, J08553+5751. All data were processed automatically using pipelines created or adapted for the survey.

We present observations of CLASS B2108+213, the widest separation gravitational lens system discovered by the Cosmic Lens All-Sky Survey. Radio imaging using the VLA at 8.46 GHz and MERLIN at 5 GHz shows two compact components separated by 4.56 arcsec with a faint third component in between which we believe is emission from a lensing galaxy. 5-GHz VLBA observations reveal milliarcsecond-scale structure in the two lensed images that is consistent with gravitational lensing. Optical emission from the two lensed images and two lensing galaxies within the Einstein radius is detected in Hubble Space Telescope imaging. Furthermore, an optical gravitational arc, associated with the strongest lensed component, has been detected. Surrounding the system are a number of faint galaxies which may help explain the wide image separation. A plausible mass distribution model for CLASS B2108+213 is also presented.

Compact groups of galaxies recently have been discovered in association with several strong gravitational lens systems. These groups provide additional convergence to the lensing potential and thus affect the value of H_0 derived from the systems. Lens system time delays are now being measured with uncertainties of only a few percent or better. Additionally, vast improvements are being made in incorporating observational constraints such as Einstein ring structures and stellar velocity dispersions into the lens models. These advances are reducing the uncertainties on H_0 to levels at which the the effects of associated galaxy groups may contribute significantly to the overall error budget. We describe a dedicated multiwavelength program, using Keck, HST, and Chandra, to find such groups and measure their properties. We present, as a case study, results obtained from observations of the CLASS lens system B1608+656 and discuss the implications for the value of H_0 derived from this system.

We are using the NRAO VLBA at 5.0 GHz to image about 200 FIRST sources stronger than 10 mJy at 1.4 GHz in the NDWFS.

We present the optical spectra of four newly discovered gravitational lenses from the Cosmic Lens All-Sky Survey (CLASS). These observations were carried out using the Low Resolution Imaging Spectrograph on the W. M. Keck-I Telescope as part of a program to study galaxy-scale gravitational lenses. From our spectra we found the redshift of the background source in CLASS B0128+437 (z_s=3.1240+-0.0042) and the lensing galaxy redshifts in CLASS B0445+123 (z_l=0.5583+-0.0003) and CLASS B0850+054 (z_l=0.5883+-0.0006). Intriguingly, we also discovered that CLASS B0631+519 may have two lensing galaxies (z_l,1=0.0896+-0.0001, z_l,2=0.6196+-0.0004). We also found a single unidentified emission line from the lensing galaxy in CLASS B0128+437 and the lensed source in CLASS B0850+054. We find the lensing galaxies in CLASS B0445+123 and CLASS B0631+519 (l,2) to be early-type galaxies with Einstein Radii of 2.8-3.0 h^{-1} kpc. The deflector in CLASS B0850+054 is a late-type galaxy with an Einstein Radius of 1.6 h^{-1} kpc.

We present results from a systematic search for strong gravitational lenses in the GOODS ACS data. The search technique involves creating a sample of likely lensing galaxies, which we define as massive early-type galaxies in a redshift range 0.3 < z <1.3. The target galaxies are selected by color and magnitude, giving a sample of 1092 galaxies. For each galaxy in the sample, we subtract a smooth description of the galaxy light from the z_{850}-band data. The residuals are examined, along with true-color images created from the B_{435}V_{606}i_{775} data, for morphologies indicative of strong lensing. We present our six most promising lens candidates, as well as our full list of candidates.

We present a refined gravitational lens model of the four-image lens system B1608+656 based on new and improved observational constraints: (i) the three independent time-delays and flux-ratios from VLA observations, (ii) the radio-image positions from VLBA observations, (iii) the shape of the deconvolved Einstein Ring from optical and infrared HST images, (iv) the extinction-corrected lens-galaxy centroids and structural parameters, and (v) a stellar velocity dispersion, sigma_ap=247+-35 km/s, of the primary lens galaxy (G1), obtained from an echelle spectrum taken with the Keck--II telescope. The lens mass model consists of two elliptical mass distributions with power-law density profiles and an external shear, totaling 22 free parameters, including the density slopes which are the key parameters to determine the value of H_0 from lens time delays. This has required the development of a new lens code that is highly optimized for speed. The minimum-chi^2 model reproduces all observations very well, including the stellar velocity dispersion and the shape of the Einstein Ring. A combined gravitational-lens and stellar dynamical analysis leads to a value of the Hubble Constant of H_0=75(+7/-6) km/s/Mpc (68 percent CL; Omega_m=0.3, Omega_Lambda=0.7. The non-linear error analysis includes correlations between all free parameters, in particular the density slopes of G1 and G2, yielding an accurate determination of the random error on H_0. The lens galaxy G1 is ~5 times more massive than the secondary lens galaxy (G2), and has a mass density slope of gamma_G1=2.03(+0.14/-0.14) +- 0.03 (68 percent CL) for rho~r^-gamma', very close to isothermal (gamma'=2). (Abridged)

We present photometric and spectroscopic observations of supernova (SN) 2002cx, which reveal it to be unique among all observed type Ia supernovae (SNe Ia). SN 2002cx exhibits a SN 1991T-like premaximum spectrum, a SN 1991bg-like luminosity, and expansion velocities roughly half those of normal SNe Ia. Photometrically, SN 2002cx has a broad peak in the $R$ band and a plateau phase in the $I$ band, and slow late-time decline. The $(B - V)$ color evolution is nearly normal, but the $(V - R)$ and $(V - I)$ colors are very red. Early-time spectra of SN 2002cx evolve very quickly and are dominated by lines from Fe-group elements; features from intermediate-mass elements (Ca, S, Si) are weak or absent. Mysterious emission lines are observed around 7000 \AA\ at about 3 weeks after maximum brightness. The nebular spectrum of SN 2002cx is also unique, consisting of narrow iron and cobalt lines. The observations of SN 2002cx are inconsistent with the observed spectral/photometric sequence, and provide a major challenge to our understanding of SNe Ia. No existing theoretical model can successfully explain all observed aspects of SN 2002cx.

We report the final results of the search for gravitationally lensed flat-spectrum radio sources found in the combination of CLASS (Cosmic Lens All Sky Survey) and JVAS (Jodrell-Bank VLA Astrometric Survey). VLA observations of 16,503 sources have been made, resulting in the largest sample of arcsec-scale lens systems available. Contained within the 16,503 sources is a complete sample of 11,685 sources having two-point spectral indices between 1.4 and 5 GHz flatter than -0.5 and 5 GHz flux densities $\geq$30 mJy. A subset of 8,958 sources form a well-defined statistical sample suitable for analysis of the lens statistics. We describe the systematic process by which 149 candidate lensed sources were picked from the statistical sample on the basis of possessing multiple compact components in the 0.2 arcsec-resolution VLA maps. Candidates were followed up with 0.05 arcsec resolution MERLIN and 0.003 arcsec VLBA observations at 5 GHz and rejected as lens systems if they failed well-defined surface brightness and/or morphological tests. Maps for all the candidates can be found on the World Wide Web at http://www.jb.man.ac.uk/research/gravlens/index.html We summarize the properties of each of the 22 gravitational lens systems in JVAS/CLASS. Twelve are double-image systems, nine are four-image systems and one is a six-image system. Thirteen constitute a statistically well-defined sample giving a point-source lensing rate of 1:690$\pm$190. The interpretation of the results in terms of the properties of the lensing galaxy population and cosmological parameters will be published elsewhere. (Abridged)

The Cosmic Lens All-Sky Survey (CLASS) is an international collaborative program which has obtained high-resolution radio images of over 10000 flat-spectrum radio sources in order to create the largest and best studied statistical sample of radio-loud gravitationally lensed systems. With this survey, combined with detailed studies of the lenses found therein, constraints can be placed on the expansion rate, matter density, and dark energy (e.g. cosmological constant, quintessence) content of the Universe that are complementary to and independent of those obtained through other methods. CLASS is aimed at identifying lenses where multiple images are formed from compact flat-spectrum radio sources, which should be easily identifiable in the radio maps. Because CLASS is radio-based, dust obscuration in lensing galaxies is not a factor, and the relative insensitivity of the instrument to environmental conditions leads to nearly uniform sensitivity and resolution over the entire survey. In four observing seasons from 1994-1999, CLASS has observed 13783 radio sources with the VLA at 8.4 GHz at 0.2 arcsecond resolution. When combined with the JVAS survey, the CLASS sample contains over 16,000 images. A complete sample of 11685 flat-spectrum sources was observed, selected from GB6 catalogue at 4.85 GHz and the NVSS at 1.4 GHz. So far, CLASS has found 16 new gravitational lens systems, and the JVAS/CLASS survey contains a total of 22 lenses. (Abridged)

We report the discovery of a new gravitational lens system from the CLASS survey. Radio observations with the VLA, the WSRT and MERLIN show that the radio source B0850+054 is comprised of two compact components with identical spectra, a separation of 0.7 arcsec and a flux density ratio of 6:1. VLBA observations at 5 GHz reveal structures that are consistent with the gravitational lens hypothesis. The brighter of the two images is resolved into a linear string of at least six sub-components whilst the weaker image is radially stretched towards the lens galaxy. UKIRT K-band imaging detects an 18.7 mag extended object, but the resolution of the observations is not sufficient to resolve the lensed images and the lens galaxy. Mass modelling has not been possible with the present data and the acquisition of high-resolution optical data is a priority for this system.

A new two-image gravitational lens system has been discovered as a result of the Cosmic Lens All-Sky Survey (CLASS). Radio observations with the VLA, MERLIN and the VLBA at increasingly higher resolutions all show two components with a flux density ratio of ~7:1 and a separation of 1.34". Both components are compact and have the same spectral index. Followup observations made with the VLA at 8.4 GHz show evidence of a feature to the south-east of the brighter component and a corresponding extension of the weaker component to the north-west. Optical observations with the WHT show ~1.7" extended emission aligned in approximately the same direction as the separation between the radio components with an R-band magnitude of 21.8 +/- 0.4.

We derive constraints on cosmological parameters and the properties of the lensing galaxies from gravitational lens statistics based on the final Cosmic Lens All Sky Survey (CLASS) data. For a flat universe with a classical cosmological constant, we find that the present matter fraction of the critical density is $\Omega_{\rm m}=0.31^{+0.27}_{-0.14}$ (68%) $^{+0.12}_{-0.10}$ (systematic). For a flat universe with a constant equation of state for dark energy $w = p_x({pressure})/\rho_x({energy density})$, we find $w < -0.55^{+0.18}_{-0.11}$ (68%).

The gravitational lens CLASS B1608+656 is the only four-image lens system for which all three independent time delays have been measured. This makes the system an excellent candidate for a high-quality determination of H_0 at cosmological distances. However, the original measurements of the time delays had large (12-20%) uncertainties, due to the low level of variability of the background source during the monitoring campaign. In this paper, we present results from two additional VLA monitoring campaigns. In contrast to the ~5% variations seen during the first season of monitoring, the source flux density changed by 25-30% in each of the subsequent two seasons. We analyzed the combined data set from all three seasons of monitoring to improve significantly the precision of the time delay measurements; the delays are consistent with those found in the original measurements, but the uncertainties have decreased by factors of two to three. We combined the delays with revised isothermal mass models to derive a measurement of H_0. Depending on the positions of the galaxy centroids, which vary by up to 0.1 arcsec in HST images obtained with different filters, we obtain H_0 = 61-65 km/s/Mpc, for (Omega_M,Omega_L) = (0.3,0.7). The value of H_0 decreases by 6% if (Omega_M,Omega_L) = (1.0,0.0). The formal uncertainties on H_0 due to the time delay measurements are +/- 1 (+/- 2) km/s/Mpc for the 1-sigma (2-sigma) confidence limits. Thus, the systematic uncertainties due to the lens model, which are on the order of +/- 15 km/s/Mpc, now dominate the error budget for this system. In order to improve the measurement of H_0 with this lens, new models that incorporate the constraints provided by stellar dynamics and the optical/infrared Einstein ring seen in HST images must be developed.

Previous observations of the environments of the lensing galaxies in gravitational lens systems suggest that many of the lensing galaxies are associated with small groups of galaxies. As a result, we have begun a coordinated program to study the local environments of all known gravitational lens systems. In this paper, we present results on the gravitational lens system CLASS B0712+472, which has previously measured source and lens redshifts of (z_l,z_s) = (0.4060,1.339). Although we have not found a galaxy group associated with the primary lensing galaxy, we have found a foreground group which is spatially coincident with the lens system. Based on multi-object spectroscopy taken at the Keck 10-m telescope, we have confirmed ten group members with a mean redshift of <z> = 0.2909. The resulting velocity dispersion and estimated virial mass are 306(+110/-58) km/s and 3.0(+2.2/-1.2) x 10^13 h^-1 M_sun, respectively, for (Omega_m,Omega_Lambda) = (0.2,0.0). The dynamical properties of this moderate-redshift group are completely consistent with the range of values found in nearby groups of galaxies. Five of the group members are red, elliptical-like galaxies, while the remaining five are active, star-forming galaxies. Based on the spectroscopic results and the publically-available HST imaging of nine group members, we find that the early-type fraction is 40%. We estimate that the effect of this foreground group on the gravitational lensing potential of B0712+472 is small, producing an external shear which is only a few percent, although the shear could be larger if the group centroid is significantly closer to the lens system than it appears to be.

We present results from the first intensive monitoring campaign of a sample of Compact Symmetric Objects (CSOs). We observed seven CSOs at 8.5 GHz over a period of eight months, with an average spacing between observations of 2.7 days. Our results show that, as predicted, the flux densities of the CSOs are extremely stable; the mean RMS variability of the sample was 0.7% in flux density. The low variability of the CSOs makes them excellent flux density calibrators at this frequency. We recommend that at least four CSOs be included in any VLA monitoring campaign which requires precise epoch-to-epoch calibration, such as those to measure gravitational lens time delays. The CSO data enable the correction of small systematic errors in the primary flux calibration.

We report the discovery of a new two-image gravitational lens system from the Cosmic Lens All-Sky Survey, CLASS B2319+051. Radio imaging with the Very Large Array (VLA) and Multi-Element Radio-Linked Interferometer Network (MERLIN) shows two compact components with a flux density ratio of 5:1, separated by 1.36 arcsec. Observations with the Very Long Baseline Array (VLBA) resolve each of the radio components into a pair of parity-reversed subcomponents. Hubble Space Telescope (HST) observations with the Near-Infrared Camera and Multi-Object Spectrometer (NICMOS) show a bright elliptical galaxy (G1) coincident with the radio position, and a second irregular galaxy (G2) 3.4 arcsec to the northwest. Previous spectroscopic studies have indicated that these galaxies are at different redshifts: z(G1) = 0.624, z(G2) = 0.588. Infrared counterparts to the lensed radio components are not detected in the NICMOS image, and the source redshift has not yet been determined. Preliminary mass modeling based on the VLBA subcomponent data indicates that the lensing potential includes a strong external shear contribution. A VLA monitoring program is currently being undertaken to measure the differential time delay.

The NRAO VLBA was used at 1.4 GHz to image the inner 550 mas (23 lt-years) of the nearest known Seyfert 1 nucleus, in the Sm galaxy NGC 4395. One continuum source was detected, with flux density $530\pm130 \mu$Jy, diameter d<11 mas (0.46 lt-year), and brightness temperature $T_{\rm b} > 2.0\times10^6$ K. The spectral power $P$ of the VLBA source is intermediate between those of Sagittarius A and Sagittarius A$^*$ in the Galactic Center. For the VLBA source in NGC 4395, the constraints on $T_b$, $d$, and $P$ are consistent with an origin from a black hole but exclude an origin from a compact starburst or a supernova remnant like Cassiopeia A. Moreover, the spectral powers of NGC 4395 at 1.4 and 4.9 GHz appear to be too low and too constant to allow analogy with SN 1988Z, a suggested prototype for models of compact supernova remnants. The variable and warm X-ray absorber in NGC 4395 has a free-free optical depth much larger than unity at 1.4 GHz and, therefore, cannot fully cover the VLBA source.

Strongly correlated non-intrinsic variability between 5 and 8.5 GHz has been observed in one of the lensed images of the gravitational lens B1600+434. These non-intrinsic (i.e. `external') variations are interpreted as `radio-microlensing' of relativistic micro-arsec-scale jet components in the source at a redshift of z=1.59 by massive compact objects in the halo of the edge-on disk lens galaxy at z=0.41. We shortly summarize these observations and discuss several new observational and theoretical programs to investigate this new phenomenon in more detail.

High resolution MERLIN observations of a newly-discovered four-image gravitational lens system, B0128+437, are presented. The system was found after a careful re-analysis of the entire CLASS dataset. The MERLIN observations resolve four components in a characteristic quadruple-image configuration; the maximum image separation is 542 mas and the total flux density is 48 mJy at 5 GHz. A best-fit lens model with a singular isothermal ellipsoid results in large errors in the image positions. A significantly improved fit is obtained after the addition of a shear component, suggesting that the lensing system is more complex and may consist of multiple deflectors. The integrated radio spectrum of the background source indicates that it is a GigaHertz-Peaked Spectrum (GPS) source. It may therefore be possible to resolve structure within the radio images with deep VLBI observations and thus better constrain the lensing mass distribution.

We present the discovery of CLASS B0739+366, a new gravitational lens system from the Cosmic Lens All-Sky Survey. Radio imaging of the source with the Very Large Array (VLA) shows two compact components separated by $0\farcs54$, with a flux density ratio of $\sim$ 6:1. High-resolution follow-up observations using the Very Long Baseline Array (VLBA) at 1.7 GHz detect weak, parity-reversed jet emission from each of the radio components. Hubble Space Telescope NICMOS F160W observations detect infrared counterparts to the lensed images, as well as an extended object between them which we identify as the lensing galaxy. Redshifts for the galaxy and lensed source have not yet been obtained. For typical lens and source redshifts of $z=0.5$ and $z=1.5$, respectively, preliminary mass modeling predicts a time delay of $\sim7h^{-1}$ days in a flat $\Omega_{M}=1.0$ universe. The small predicted time delay and weak radio components will make CLASS B0739+366 a challenging target for Hubble constant determination.

We present radio, optical, near-infrared and spectroscopic observations of the source B0827+525. We consider this source as the best candidate from the Cosmic Lens All-Sky Survey (CLASS) for a `dark lens' system or binary radio-loud quasar. The system consists of two radio components with somewhat different spectral indices, separated by 2.815 arcsec. VLBA observations show that each component has substructure on a scale of a few mas. A deep K-band exposure with the W.M.Keck-II Telescope reveals emission near both radio components. The K-band emission of the weaker radio component appears extended, whereas the emission from the brighter radio component is consistent with a point source. Hubble Space Telescope F160W-band observations with the NICMOS instrument confirms this. A redshift of 2.064 is found for the brighter component, using the LRIS instrument on the W.M.Keck-II Telescope. The probability that B0827+525 consists of two unrelated compact flat-spectrum radio sources is ~3%, although the presence of similar substructure in both component might reduce this. We discuss two scenarios to explain this system: (i) CLASS B0827+525 is a `dark lens' system or (ii) B0827+525 is a binary radio-loud quasar. B0827+525 has met all criteria that thus far have in 100% of the cases confirmed a source as an indisputable gravitational lens system. Despite this, no lens galaxy has been detected with m_F160W<=23 mag. Hence, we might have found the first binary radio-loud quasar. At this moment, however, we feel that the `dark lens' hypothesis cannot yet be fully excluded.

The complex ten-component gravitational lens system B1933+503 has been monitored with the VLA during the period February to June 1998 with a view to measuring the time delay between the four compact components and hence to determine the Hubble parameter. Here we present the results of an `A' configuration 8.4-GHz monitoring campaign which consists of 37 epochs with an average spacing of 2.8 days. The data have yielded light curves for the four flat-spectrum radio components (components 1, 3, 4 and 6). We observe only small flux density changes in the four flat-spectrum components which we do not believe are predominantly intrinsic to the source. Therefore the variations do not allow us to determine the independent time delays in this system. However, the data do allow us to accurately determine the flux density ratios between the four flat-spectrum components. These will prove important as modelling constraints and could prove crucial in future monitoring observations should these data show only a monotonic increase or decrease in the flux densities of the flat-spectrum components.

First, we review the current status of the detection of strong `external' variability in the CLASS gravitational B1600+434, focusing on the 1998 VLA 8.5-GHz and 1998/9 WSRT multi-frequency observations. We show that this data can best be explained in terms of radio-microlensing. We then proceed to show some preliminary results from our new multi-frequency VLA monitoring program, in particular the detection of a strong feature (~30%) in the light curve of the lensed image which passes predominantly through the dark-matter halo of the lens galaxy. We tentatively interpret this event, which lasted for several weeks, as a radio-microlensing caustic crossing, i.e. the superluminal motion of a micro-arcsec-scale jet-component in the lensed source over a single caustic in the magnification pattern, that has been created by massive compact objects along the line-of-sight to the lensed image.

We present Very Large Array (VLA) 8.5-GHz light curves of the two lens images of the Cosmic Lens All Sky Survey (CLASS) gravitational lens B1600+434. We find a nearly linear decrease of 18-19% in the flux densities of both lens images over a period of eight months (February-October) in 1998. Additionally, the brightest image A shows modulations up to 11% peak-to-peak on scales of days to weeks over a large part of the observing period. Image B varies significantly less on this time scale. We conclude that most of the short-term variability in image A is not intrinsic source variability, but is most likely caused by microlensing in the lens galaxy. The alternative, scintillation by the ionized Galactic ISM, is shown to be implausible based on its strong opposite frequency dependent behavior compared with results from multi-frequency WSRT monitoring observations (Koopmans & de Bruyn 1999). From these VLA light curves we determine a median time delay between the lens images of 47^{+5}_{-6} d (68%) or 47^{+12}_{-9} d (95%). We use two different methods to derive the time delay; both give the same result within the errors. We estimate an additional systematic error between -8 and +7 d. If the mass distribution of lens galaxy can be described by an isothermal model (Koopmans, de Bruyn & Jackson 1998), this time delay would give a value for the Hubble parameter, H_0=57^{+14}_{-11} (95% statistical) ^{+26}_{-15} (systematic) km/s/Mpc (Omega_m=1 and Omega_Lambda=0). Similarly, the Modified-Hubble-Profile mass model would give H_0=74^{+18}_{-15} (95% statistical) ^{+22}_{-22} (systematic) km/s/Mpc. For Omega_m=0.3 and Omega_Lambda=0.7, these values increase by 5.4%. ... (ABRIDGED)

Compact Symmetric Objects (CSOs) are young radio galaxies whose jet axes lie close to the plane of the sky, and whose appearance is therefore not dominated by relativistic beaming effects. The small linear sizes of CSOs make them valuable for studies of both the evolution of radio galaxies and testing unified schemes for active galactic nuclei (AGN). A parsec-scale region of gas surrounding the central engine is predicted by both accretion and obscuration scenarios. Working surfaces, or ``hot spots,'' and the radio jets of CSOs are close enough to the central engines that this circumnuclear gas can be seen in absorption. The CSOs Observed in the Northern Sky (COINS) sample is comprised of 52 CSO candidates identified in three VLBI surveys. Of these, 27 have now been confirmed as CSOs. Optical redshifts are available in the literature for 28 of the CSO candidates, and HI absorption has been detected toward four. We present new optical spectroscopic redshifts for three of the candidates and summarize the current status of optical identifications. We further report on the discovery of HI in absorption towards the CSO J1816+3457 and summarize the results of neutral hydrogen absorption studies of the sources in this sample.

We present new results from a continuing Keck program to study gravitational lens systems. We have obtained redshifts for three lens systems, SBS 0909+532, HST 1411+5211, and CLASS B2319+051. For all of these systems, either the source or lens redshift (or both) has been previously unidentified. We find (z_l, z_s) = (0.830, 1.377) for SBS 0909+532; (z_l, z_s) = (0.465, 2.811) for HST 1411+5211, although the source redshift is still tentative; and (z_l1, z_l2) = (0.624, 0.588) for the two lensing galaxies in CLASS B2319+051. The background radio source in B2319+051 has not been detected optically; its redshift is, therefore, still unknown. We find that the spectral features of the central lensing galaxy in all three systems are typical of an early-type galaxy. The observed image splittings in SBS 0909+532 and HST 1411+5211 imply that the masses within the Einstein ring radii of the lensing galaxies are 1.4 x 10^{11} and 2.0 x 10^{11} h^{-1} M_sun, respectively. The resulting B band mass-to-light ratio for HST 1411+5211 is 41.3 +/- 1.2 h (M/L)_sun, a factor of 5 times higher than the average early-type lensing galaxy. This large mass-to-light is almost certainly the result of the additional mass contribution from the cluster CL 3C295 at z = 0.46. For the lensing galaxy in SBS 0909+532, we measure (M/L)_B = 4^{+11}_{-3} h (M/L)_sun where the large errors are the result of significant uncertainty in the galaxy luminosity. While we cannot measure directly the mass-to-light ratio of the lensing galaxy in B2319+051, we estimate that (M/L)_B is between 3-7 h (M/L)_sun.

The four-component gravitational lens CLASS B1608+656 has been monitored with the VLA for two seasons in order to search for time delays between the components. These time delays can be combined with mass models of the lens system to yield a measurement of H_0. The component light curves show significantly different behavior in the two observing seasons. In the first season the light curves have maximum variations of ~5%, while in the second season the components experienced a nearly monotonic ~40% decrease in flux. We present the time delays derived from a joint analysis of the light curves from the two seasons.

We present the results of a program to monitor the four-image gravitational lens B1608+656 with the VLA. The system was observed over a seven month period from 1996 October to 1997 May. The 64 epochs of observation have an average spacing of 3.6~d. The light curves of the four images of the background source show that the flux density of the background source has varied at the ~5% level. We measure time delays in the system based on common features that are seen in all four light curves. The three independent time delays in the system are found to be Delta t_{BA} = 31 +/- 7~d, Delta t_{BC} = 36 +/- 7~d, and Delta t_{BD} = 76^{+9}_{-10}~d at 95% confidence. This is the first gravitational lens system for which three independent time delays have been measured. A companion paper presents a mass model for the lensing galaxy which correctly reproduces the observed image positions, flux density ratios, and time delay ratios. The last condition is crucial for determining H_0 with a four-image lens. We combine the time delays with the model to obtain a value for the Hubble constant of H_0 = 59^{+8}_{-7} km/s/Mpc at 95% confidence (statistical) for (Omega_M, Omega_{Lambda}) = (1,0). In addition, there is an estimated systematic uncertainty of +/- 15 km/s/Mpc from uncertainties in modeling the radial mass profiles of the lensing galaxies. The value of H_0 presented in this paper is comparable to recent measurements of H_0 from the gravitational lenses 0957+561, PG1115+080, B0218+357, and PKS1830-211.

EDITED FROM PAPER: We present mass models of the four-image gravitational lens system B1608+656. A mass model for the lens galaxies has been determined that reproduces the image positions, two out of three flux-density ratios and the model time delays. Using the time delays determined by Fassnacht et al. (1999a), we find that the best isothermal mass model gives H_0=59^{+7}_{-6} km/s/Mpc for Omega_m=1 and Omega_l=0.0, or H_0=(65-63)^{+7}_{-6} km/s/Mpc for Omega_m=0.3 and Omega_l = 0.0-0.7 (95.4% statistical confidence). A systematic error of +/-15 km/s/Mpc is estimated. This cosmological determination of H_0 agrees well with determinations from three other gravitational lens systems (i.e. B0218+357, Q0957+561 and PKS1830-211), SNe Ia, the S-Z effect and local determinations. The current agreement on H_0 from four out of five gravitational lens systems (i) emphasizes the reliability of its determination from isolated gravitational lens systems and (ii) suggests that a close-to-isothermal mass profile can describe disk galaxies, ellipticals and central cluster ellipticals. The average of H_0 from B0218+357, Q0957+561, B1608+656 and PKS1830-211, gives H_0(GL)=69 +/-7 km/s/Mpc for a flat universe with Omega_m=1 or H_0(GL)=74 +/-8 km/s/Mpc for Omega_m=0.3 and Omega_l=0.0-0.7. When including PG1115+080, these values decrease to 64 +/-11 km/s/Mpc and 68 +/-13 km/s/Mpc (2-sigma errors), respectively.

The third phase of the Cosmic Lens All-Sky Survey (CLASS) has recently been completed, bringing the total number of sources imaged to over 15000 in the CLASS and JVAS combined survey. In the VLA observations carried out in March and April of 1998, two new candidate lensed systems were discovered: CLASS B1152+199 and B1359+154. B1152+199 is a 1.6 arcsecond double, with a background quasar at z=1.019 lensed by a foreground galaxy at z=0.439. The relatively flat radio spectra of the lensed images, combined with a previous ROSAT detection of the source, make B1152+199 a strong candidate for time delay studies at both radio and X-ray wavelengths. B1359+154 is a quadruply lensed quasar at z=3.235, with a maximum image separation of 1.7 arcseconds. As yet, the redshift of the lensing object in this system is undetermined. The steep spectral index of the source suggests that B1359+154 will not exhibit strong variability, and is therefore unlikely to be useful for determining the Hubble constant from measured time delays.

The JVAS (Jodrell Bank-VLA Astrometric Survey) and CLASS (Cosmic Lens All-Sky Survey) are well-defined surveys containing about ten thousand flat-spectrum radio sources. For many reasons, flat-spectrum radio sources are particularly well-suited as a population from which one can obtain unbiased samples of gravitational lenses. These are by far the largest gravitational (macro)lens surveys, and particular attention was paid to constructing a cleanly-defined sample for the survey itself and for the underlying luminosity function. Here we present the constraints on cosmological parameters, particularly the cosmological constant, derived from JVAS and combine them with constraints from optical gravitational lens surveys, `direct' measurements of $\Omega_{0}$, $H_{0}$ and the age of the universe, and constraints derived from CMB anisotropies, before putting this final result into the context of the latest results from other, independent cosmological tests.

We present the discovery of a new gravitational lens system with two compact radio images separated by 0.701+-0.001 arcsec. The lens system was discovered in the Cosmic Lens All Sky Survey (CLASS) as a flat spectrum radio source. Both radio components show structure in a VLBA 8.4 GHz radio image. No further extended structure is seen in either the VLA, MERLIN or VLBA images. Hubble Space Telescope (HST) WFPC2 images in F555W and F814W show two extended objects close to the radio components, which we identify as two lens galaxies. Their colours and mass-to-light ratios seem to favour two late-type spiral galaxies at relatively high redshifts (z_d>~0.5). Faint emission is also detected at positions corresponding to the radio images. A two-lens mass model can explain the observed VLBA structure. The best fit model has a reduced chi^2 of 1.1. The relative positions of the VLBA subcomponents are reproduced within 0.08 mas, the flux density ratios within 0.19. We also reproduce the position angle and separation of the two VLBA subcomponents in A and B within the observational errors, which we consider strong evidence for the validity of the lens model. Moreover, we find a surface density axis ratio of 0.74^{+0.10}_{-0.12} for the primary lens (G1), consistent with the surface brightness axis ratio of 0.69+-0.15. Also, the surface density position angle of 4.9^{+28.2}_{-22.4} degree of G1 compares well with the -6+-13 degree position angle of the surface brightness distribution. The errors indicate the 99 per cent confidence interval.

We have discovered a new gravitational lens in the Cosmic Lens All-Sky Survey (CLASS). B2045+265 is a four-image system with a maximum separation of 1.9 arcsec. A fifth radio component is detected, but its radio spectrum and its positional coincidence with infrared emission from the lensing galaxy strongly suggests that it is the radio core of the lensing galaxy. This implies that the B2045+265 system consists of a flat-spectrum radio source being lensed by another flat-spectrum radio source. Infrared images detect the lensed images of the background source and the lensing galaxy. The lensed images have relative positions and flux densities that are consistent with those seen at radio wavelengths. Spectra of the system reveal a lens redshift of z_l=0.8673 and a source redshift of z_s=1.28. The image splitting and system redshifts imply that the projected mass inside the Einstein radius of the lensing galaxy is M_E=4.7x10^{11} h^{-1} M_sun. An estimate of the light emitted inside the Einstein radius gives a mass-to-light ratio in the rest frame B band of (M/L_B)_E = 20 h (M/L_B)_sun. Both the mass and mass-to-light ratio are higher than what is seen in nearby Sa galaxies. In fact, the implied rotation velocity for the lensing galaxy is two to three times higher than what is seen in nearby spirals. The large projected mass inside the Einstein ring radius may be the result of a significant amount of dark matter in the system, perhaps from a compact group of galaxies; however, it may also arise from a misidentification of the source redshift. A simple model of the gravitational potential of the lens reproduces the image positions well, but further modeling is required to satisfy the constraints from the image flux density ratios. With further observations and modeling, this lens may yield an estimate of H_0.

In a recent paper, Hawkins (1997) argues on the basis of statistical studies of double-image gravitational lenses and lens candidates that a large population of dark lenses exists and that these outnumber galaxies with more normal mass-to-light ratios by a factor of 3:1. If correct, this is a very important result for many areas of astronomy including galaxy formation and cosmology. In this paper we discuss our new radio-selected gravitational lens sample, JVAS/CLASS, in order to test and constrain this proposition. We have obtained ground-based and HST images of all multiple-image lens systems in our sample and in 12 cases out of 12 we find the lensing galaxies in the optical and/or near infrared. Our success in finding lensing galaxies creates problems for the dark lens hypothesis. If it is to survive, ad hoc modifications seem to be necessary: only very massive galaxies (more than about one trillion solar masses) can be dark, and the cutoff in mass must be sharp. Our finding of lens galaxies in all the JVAS/CLASS systems is complementary evidence which supports the conclusion of Kochanek et al. (1997) that many of the wide-separation optically-selected pairs are physically distinct quasars rather than gravitational lens systems.

We present spectra of three gravitational lens systems taken with the Low Resolution Imaging Spectrograph on the W. M. Keck Telescopes. All of the systems were discovered in the JVAS and CLASS radio surveys, which were designed to find lenses suitable for measuring $H_0$. Previous spectra of these systems had low signal-to-noise ratios, and only one of the source redshifts was secure. Our observations give unambiguous lens and source redshifts for all of the systems, with ($z_l$, $z_s$) = (0.4060,1.339), (0.5990,1.535) and (0.4144,1.589) for B0712+472, B1030+074 and B1600+434, respectively. The observed image splittings in the systems imply that the masses of the lensing galaxies within their Einstein rings are 5.4$\times 10^{10}$, 1.2$\times 10^{11}$, and 6.3$\times 10^{10} h^{-1} M_{\sun}$. The resulting V-band mass-to-light ratios for B0712+472 and B1030+074, measured inside their Einstein ring radii, are $\sim 10h (M/L)_{\sun, V}$, slightly higher than values observed in nearby ellipticals. For B1600+434, the mass-to-light ratio is $48h (M/L)_{\sun, V}$. This high value can be explained, at least in part, by the prominent dust lane running through the galaxy. Two of the three lens systems show evidence of variability, so monitoring may yield a time delay and thus a measurement of $H_0$.

We report the discovery of the most complex arcsec-scale radio gravitational lens system yet known. B1933+503 was found during the course of the CLASS survey and MERLIN and VLA radio maps reveal up to 10 components. Four of these are compact and have flat spectra; the rest are more extended and have steep spectra. The background lensed object appears to consist of a flat spectrum core (quadruply imaged) and two compact "lobes" symmetrically disposed relative to the core. One of the lobes is quadruply imaged while the other is doubly imaged. An HST observation of the system with the WFPC2 shows a galaxy with an axial ratio of 0.5, but none of the images of the background object are detected. A redshift of 0.755 has been measured for the lens galaxy.

We report the discovery, using NICMOS on the Hubble Space Telescope, of an arcsecond-diameter Einstein ring in the gravitational lens system B1938+666. The lensing galaxy is also detected, and is most likely an early-type. Modelling of the ring is presented and compared with the radio structure from MERLIN maps. We show that the Einstein ring is consistent with the gravitational lensing of an extended infrared component, centred between the two radio components.

The discovery of HDF~J123652+621227, a candidate gravitational lens in the HDF, is reported. This lens may be multiply imaging several optical sources at different redshifts. If follow-up spectroscopy of the lens and the brightest image confirms this hypothesis, observations of this system alone can be used to obtain an estimate of the redshift distribution at extremely faint flux levels.