The Hubble Frontier Fields data, along with multiple data sets obtained by other telescopes, have provided some of the most extensive constraints on cluster lenses to date. Multiple lens modeling teams analyzed the fields and made public a number of deliverables. By comparing these results, we can then undertake a unique and vital test of the state of cluster lens modeling. Specifically, we see how well the different teams can reproduce similar magnifications and mass profiles. We find that the circularly averaged mass profiles of the fields are remarkably constrained (scatter <5%) at distances of 1 arcmin from the cluster core, yet magnifications can vary significantly. Averaged across the six fields, we find a bias of -6% (-17%) and a scatter of ~40% (~65%) at a modest magnification of 3 (10). Statistical errors reported by individual teams are often significantly smaller than the differences among all the teams, indicating the importance of continued systematics studies in cluster lensing.

Cluster lensing has become an important tool in the search for high redshift galaxies through its ability to magnify sources. In order to determine the intrinsic properties of these galaxies, lensing mass models must be constructed to determine the magnification of the images. These models are traditionally two-dimensional, focusing on the mass within the cluster and either ignoring or approximating any contribution from line-of-sight galaxies. In this paper, we present the first full set of three-dimensional mass models of the six Hubble Frontier Fields and use them to test for systematic biases in magnifications due to using the traditional 2D approach. We find that omitting foreground or background galaxies causes image position offsets between 0.1-0.4", a non-negligible fraction of the typical 0.3-0.7" residuals of current state-of-the-art models. We also find that median image magnifications can shift by up to 6%, though it is dependent on the field. This can be alleviated in some cases by approximating the mass in the lensing plane, but a 5% magnification bias still exists in other cases; image position offsets are also improved, but are still present at 0.10".

In the cold dark matter (CDM) picture of structure formation, galaxy mass distributions are predicted to have a considerable amount of structure on small scales. Strong gravitational lensing has proven to be a useful tool for studying this small-scale structure. Much of the attention has been given to detecting individual dark matter subhalos through lens modeling, but recent work has suggested that the full population of subhalos could be probed using a power spectrum analysis. In this paper we quantify the power spectrum of small-scale structure in simulated galaxies, with the goal of understanding theoretical predictions and setting the stage for using measurements of the power spectrum to test dark matter models. We use a sample of simulated galaxies generated from the \texttt{Galacticus} semi-analytic model to determine the power spectrum distribution first in the CDM paradigm and then in a warm dark matter scenario. We find that a measurement of the slope and amplitude of the power spectrum on galaxy strong lensing scales ($k\sim 1$ kpc$^{-1}$) could be used to distinguish between CDM and alternate dark matter models, especially if the most massive subhalos can be directly detected via gravitational imaging.

Gravitational lensing has emerged as a powerful probe of the matter distribution on subgalactic scales, which itself may contain important clues about the fundamental origins and properties of dark matter. Broadly speaking, two different approaches have been taken in the literature to map the small-scale structure of the Universe using strong lensing, with one focused on measuring the position and mass of a small number of discrete massive subhalos appearing close in projection to lensed images, and the other focused on detecting the collective effect of all the small-scale structure between the lensed source and the observer. In this paper, we follow the latter approach and perform a detailed study of the sensitivity of galaxy-scale gravitational lenses to the ensemble properties of small-scale structure. As in some previous studies, we adopt the language of the substructure power spectrum to characterize the statistical properties of the small-scale density field. We present a comprehensive theory that treats lenses with extended sources as well as those with time-dependent compact sources (such as quasars) in a unified framework for the first time. Our approach uses mode functions to provide both computational advantages and insights about couplings between the lens and source. The goal of this paper is to develop the theory and gain the intuition necessary to understand how the sensitivity to the substructure power spectrum depends on the source and lens properties, with the eventual aim of identifying the most promising targets for such studies.

We show that the asymptotic scaling of the magnification volume cross section corresponding to an elliptic umbilic caustic surface is $\mu^{-2.5}$ in the two-image region and $\mu^{-2}$ in the four-image region, where $\mu$ is the total unsigned magnification.

We present IRAM/NOEMA, JCMT/SCUBA-2 and VLA observations of the most distant known gravitationally lensed quasar J0439+1634 at z = 6.5. We detect strong dust emission, [CII] 158 $\mu$m, [CI] 369 $\mu$m, [OI] 146 $\mu$m, CO(6-5), CO(7-6), CO(9-8), CO(10-9), H$_{\rm 2}$O $3_{\rm 1,2}-2_{\rm 2,1}$, and H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ lines as well as a weak radio continuum. The strong [CII] line yields a systemic redshift of the host galaxy to be z=6.5188$\pm$0.0001. The magnification makes J0439+1634 the far-infrared (FIR) brightest quasar at z > 6 known, with the brightest [CII] line yet detected at this redshift. The FIR luminosity is (3.4$\pm$0.2)$\times$10$^{13}$ $\mu^{-1}$ $L_{\odot}$, where $\mu$ $\sim$ 2.6 - 6.6 is the magnification of the host galaxy, estimated based on the lensing configuration from HST imaging. We estimate the dust mass to be (2.2$\pm$0.1)$\times$10$^{9}$ $\mu^{-1}$ $M_{\odot}$. The CO Spectral Line Energy Distribution using four CO lines is best fit by a two-component model of the molecular gas excitation. The estimates of molecular gas mass derived from CO lines and atomic carbon mass are consistent, in the range of 3.9 - 8.9 $\times$10$^{10} \mu^{-1}$ $M_{\odot}$. The [CII]/[CI], [CII]/CO, and [OI]/[CII] line luminosity ratios suggest a photodissociation region model with more than one component. The ratio of H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ line luminosity to $L_{\rm TIR}$ is consistent with values in local and high redshift ultra-/hyper-luminous infrared galaxies. The VLA observations reveal an unresolved radio continuum source, and indicate that J0439+1634 is a radio quiet quasar with R = 0.05 - 0.17.

We report Northern Extended Millimeter Array (NOEMA) CO($J = 3 - 2$) observations of the dusty star-forming galaxy ACT-S\,J020941+001557 at $z = 2.5528$, which was detected as an unresolved source in the Atacama Cosmology Telescope (ACT) equatorial survey. Our spatially resolved spectral line data support the derivation of a gravitational lens model from 37 independent velocity channel maps using a pixel-based algorithm, from which we infer a velocity-dependent magnification factor $\mu \approx 7-22$ with a luminosity-weighted mean $\left<\mu\right>\approx 13$. The resulting source-plane reconstruction is consistent with a rotating disk, although other scenarios cannot be ruled out by our data. After correction for lensing, we derive a line luminosity $L^{\prime}_{\rm CO(3-2)}= (5.53\pm 0.69) \times 10^{10}\,{\rm \,K\,km\,s^{-1}\,pc^{2}}$, a cold gas mass $M_{{\rm gas}}= (3.86 \pm 0.33) \times 10^{10}\,M_{\odot}$, a dynamical mass $M_{\rm dyn}\,{\rm sin}^2\,i = 3.9^{+1.8}_{-1.5} \times 10^{10}\,M_{\odot}$, and a gas mass fraction $f_{\rm gas}\,{\rm csc}^2\,i = 1.0^{+0.8}_{-0.4}$. The line brightness temperature ratio of $r_{3,1}\approx 1.6$ relative to a Green Bank Telescope CO($J=1-0$) detection may be elevated by a combination of external heating of molecular clouds, differential lensing, and/or pointing errors.

We present ~1" resolution (~2 kpc in the source plane) observations of the CO(1-0), CO(3-2), Halpha, and [N II] lines in the strongly-lensed z=2.26 star-forming galaxy SDSS J0901+1814. We use these observations to constrain the lensing potential of a foreground group of galaxies, and our source-plane reconstructions indicate that SDSS J0901+1814 is a nearly face-on (i~30 degrees) massive disk with r_{1/2}>~4 kpc for its molecular gas. Using our new magnification factors (mu_tot~30), we find that SDSS J0901+1814 has a star formation rate (SFR) of 268^{+63}_{-61} M_sun/yr, M_gas=(1.6^{+0.3}_{-0.2})x10^11x(alpha_CO/4.6) M_sun, and M_star=(9.5^{+3.8}_{-2.8})x10^10 M_sun, which places it on the star-forming galaxy "main sequence." We use our matched high-angular resolution gas and SFR tracers (CO and Halpha, respectively) to perform a spatially resolved (pixel-by-pixel) analysis of SDSS J0901+1814 in terms of the Schmidt-Kennicutt relation. After correcting for the large fraction of obscured star formation (SFR_Halpha/SFR_TIR=0.054^{+0.015}_{-0.014}), we find SDSS J0901+1814 is offset from "normal" star-forming galaxies to higher star formation efficiencies independent of assumptions for the CO-to-H_2 conversion factor. Our mean best-fit index for the Schmidt-Kennicutt relation for SDSS J0901+1814, evaluated with different CO lines and smoothing levels, is n=1.54+/-0.13; however, the index may be affected by gravitational lensing, and we find n=1.24+/-0.02 when analyzing the source-plane reconstructions. While the Schmidt-Kennicutt index largely appears unaffected by which of the two CO transitions we use to trace the molecular gas, the source-plane reconstructions and dynamical modeling suggest that the CO(1-0) emission is more spatially extended than the CO(3-2) emission.

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

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.

The degeneracy between disc and halo contributions in spiral galaxy rotation curves makes it difficult to obtain a full understanding of the distribution of baryons and dark matter in disc galaxies like our own Milky Way. Using mock data, we study how constraints on dark matter profiles obtained from kinematics, strong lensing, or a combination of the two are affected by assumptions about the halo model. We compare four different models: spherical isothermal and Navarro-Frenk-White halos, along with spherical and elliptical Burkert halos. For both kinematics and lensing we find examples where different models fit the data well but give enclosed masses that are inconsistent with the true (i.e., input) values. This is especially notable when the input and fit models differ in having cored or cuspy profiles (such as fitting an NFW model when the underlying dark matter distribution follows a different profile). We find that mass biases are more pronounced with lensing than with kinematics, and using both methods can help reduce the bias and provide stronger constraints on the dark matter distributions.

Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first such object detected at the epoch of reionization, and the brightest quasar yet known at z>5. High-resolution HST imaging reveals a multiple imaged system with a maximum image separation theta ~ 0.2", best explained by a model of three quasar images lensed by a low luminosity galaxy at z~0.7, with a magnification factor of ~50. The existence of this source suggests that a significant population of strongly lensed, high redshift quasars could have been missed by previous surveys, as standard color selection techniques would fail when the quasar color is contaminated by the lensing galaxy.

The free streaming length of dark matter particles determines the abundance of structure on sub-galactic scales. We present a statistical technique, amendable to any parameterization of subhalo density profile and mass function, to probe dark matter on these scales with quadrupole image lenses. We consider a warm dark matter particle with a mass function characterized by a normalization and free streaming scale $m_{\rm{hm}}$. We forecast bounds on dark matter warmth for 120-180 lenses, attainable with future surveys, at typical lens (source) redshifts of 0.5 (1.5) in early-type galaxies with velocity dispersions of 220-270 km/sec. We demonstrate that limits on $m_{\rm{hm}}$ deteriorate rapidly with increasing uncertainty in image fluxes, underscoring the importance of precise measurements and accurate lens models. For our forecasts, we assume the deflectors in the lens sample do not exhibit complex morphologies, so we neglect systematic errors in their modeling. Omitting the additional signal from line of sight halos, our constraints underestimate the true power of the method. Assuming cold dark matter, for a low normalization, corresponding the destruction of all subhalos within the host scale radius, we forecast $2\sigma$ bounds on $m_{\rm{hm}}$ (thermal relic mass) of $10^{7.5} \ (5.0)$, $10^{8} \ (3.6)$, and $10^{8.5} \ (2.7) \ M_{\odot} \left(\rm{keV}\right)$ for flux errors of $2\%$, $4\%$, and $8\%$. With a higher normalization, these constraints improve to $10^{7.2} \ (6.6)$, $10^{7.5} \ (5.3) $, and $10^{7.8} \ (4.3) \ M_{\odot} \left(\rm{keV}\right)$ with 120 systems. We are also able to measure the normalization of the mass function, which has implications for baryonic feedback models and tidal stripping.

With a large, unique spectroscopic survey in the fields of 28 galaxy-scale strong gravitational lenses, we identify groups of galaxies in the 26 adequately-sampled fields. Using a group finding algorithm, we find 210 groups with at least five member galaxies; the median number of members is eight. Our sample spans redshifts of 0.04 $\le z_{grp} \le$ 0.76 with a median of 0.31, including 174 groups with $0.1 < z_{grp} < 0.6$. Groups have radial velocity dispersions of 60 $\le \sigma_{grp} \le$ 1200 km s$^{-1}$ with a median of 350 km s$^{-1}$. We also discover a supergroup in field B0712+472 at $z =$ 0.29 consisting of three main groups. We recover groups similar to $\sim$ 85% of those previously reported in these fields within our redshift range of sensitivity and find 187 new groups with at least five members. The properties of our group catalog, specifically 1) the distribution of $\sigma_{grp}$, 2) the fraction of all sample galaxies that are group members, and 3) the fraction of groups with significant substructure, are consistent with those for other catalogs. The distribution of group virial masses agrees well with theoretical expectations. Of the lens galaxies, 12 of 26 (46%) (B1422+231, B1600+434, B2114+022, FBQS J0951+2635, HE0435-1223, HST J14113+5211, MG0751+2716, MGJ1654+1346, PG 1115+080, Q ER 0047-2808, RXJ1131-1231, and WFI J2033-4723) are members of groups with at least five galaxies, and one more (B0712+472) belongs to an additional, visually identified group candidate. There are groups not associated with the lens that still are likely to affect the lens model; in six of 25 (24%) fields (excluding the supergroup), there is at least one massive ($\sigma_{grp} \ge$ 500 km s$^{-1}$) group or group candidate projected within 2$^{\prime}$ of the lens.

Strong gravitational lensing provides an independent measurement of the Hubble parameter ($H_0$). One remaining systematic is a bias from the additional mass due to a galaxy group at the lens redshift or along the sightline. We quantify this bias for more than 20 strong lenses that have well-sampled sightline mass distributions, focusing on the convergence $\kappa$ and shear $\gamma$. In 23% of these fields, a lens group contributes a $\ge$1% convergence bias; in 57%, there is a similarly significant line-of-sight group. For the nine time delay lens systems, $H_0$ is overestimated by 11$^{+3}_{-2}$% on average when groups are ignored. In 67% of fields with total $\kappa \ge$ 0.01, line-of-sight groups contribute $\gtrsim 2\times$ more convergence than do lens groups, indicating that the lens group is not the only important mass. Lens environment affects the ratio of four (quad) to two (double) image systems; all seven quads have lens groups while only three of 10 doubles do, and the highest convergences due to lens groups are in quads. We calibrate the $\gamma$-$\kappa$ relation: $\log(\kappa_{\rm{tot}}) = (1.94 \pm 0.34) \log(\gamma_{\rm{tot}}) + (1.31 \pm 0.49)$ with a rms scatter of 0.34 dex. Shear, which, unlike convergence, can be measured directly from lensed images, can be a poor predictor of $\kappa$; for 19% of our fields, $\kappa$ is $\gtrsim 2\gamma$. Thus, accurate cosmology using strong gravitational lenses requires precise measurement and correction for all significant structures in each lens field.

Magnetic fields play a pivotal role in the physics of interstellar medium in galaxies, but there are few observational constraints on how they evolve across cosmic time. Spatially resolved synchrotron polarization maps at radio wavelengths reveal well-ordered large-scale magnetic fields in nearby galaxies that are believed to grow from a seed field via a dynamo effect. To directly test and characterize this theory requires magnetic field strength and geometry measurements in cosmologically distant galaxies, which are challenging to obtain due to the limited sensitivity and angular resolution of current radio telescopes. Here, we report the cleanest measurements yet of magnetic fields in a galaxy beyond the local volume, free of the systematics traditional techniques would encounter. By exploiting the scenario where the polarized radio emission from a background source is gravitationally lensed by a foreground galaxy at z = 0.439 using broadband radio polarization data, we detected coherent $\mu$G magnetic fields in the lensing disk galaxy as seen 4.6 Gyrs ago, with similar strength and geometry to local volume galaxies. This is the highest redshift galaxy whose observed coherent magnetic field property is compatible with a mean-field dynamo origin.

We report the serendipitous discovery of the first gravitationally lensed quasar candidate from Pan-STARRS. The grizy images reveal four point-like images with magnitudes between 14.9 mag and 18.1 mag. The colors of the point sources are similar, and they are more consistent with quasars than with stars or galaxies. The lensing galaxy is detected in the izy bands, with an inferred photometric redshift of ~0.6, lower than that of the point sources. We successfully model the system with a singular isothermal ellipsoid with shear, using the relative positions of the five objects as constraints. While the brightness ranking of the point sources is consistent with that of the model, we find discrepancies between the model-predicted and observed fluxes, likely due to microlensing by stars and millilensing due to the dark matter substructure. In order to fully confirm the gravitational lens nature of this system and add it to the small but growing number of the powerful probes of cosmology and astrophysics represented by quadruply lensed quasars, we require further spectroscopy and high-resolution imaging.

We present a combined strong and weak lensing analysis of the J085007.6+360428 (J0850) field, which was selected by its high projected concentration of luminous red galaxies and contains the massive cluster Zwicky 1953. Using Subaru/Suprime-Cam $BVR_{c}I_{c}i^{\prime}z^{\prime}$ imaging and MMT/Hectospec spectroscopy, we first perform a weak lensing shear analysis to constrain the mass distribution in this field, including the cluster at $z = 0.3774$ and a smaller foreground halo at $z = 0.2713$. We then add a strong lensing constraint from a multiply-imaged galaxy in the imaging data with a photometric redshift of $z \approx 5.03$. Unlike previous cluster-scale lens analyses, our technique accounts for the full three-dimensional mass structure in the beam, including galaxies along the line of sight. In contrast with past cluster analyses that use only lensed image positions as constraints, we use the full surface brightness distribution of the images. This method predicts that the source galaxy crosses a lensing caustic such that one image is a highly-magnified "fold arc", which could be used to probe the source galaxy's structure at ultra-high spatial resolution ($< 30$ pc). We calculate the mass of the primary cluster to be $\mathrm{M_{vir}} = 2.93_{-0.65}^{+0.71} \times 10^{15}~\mathrm{M_{\odot}}$ with a concentration of $\mathrm{c_{vir}} = 3.46_{-0.59}^{+0.70}$, consistent with the mass-concentration relation of massive clusters at a similar redshift. The large mass of this cluster makes J0850 an excellent field for leveraging lensing magnification to search for high-redshift galaxies, competitive with and complementary to that of well-studied clusters such as the HST Frontier Fields.

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.

Matter near a gravitational lens galaxy or projected along the line of sight (LOS) can affect strong lensing observables by more than contemporary measurement errors. We simulate lens fields with realistic three-dimensional mass configurations (self-consistently including voids), and then fit mock lensing observables with increasingly complex lens models to quantify biases and uncertainties associated with different ways of treating the lens environment (ENV) and LOS. We identify the combination of mass, projected offset, and redshift that determines the importance of a perturbing galaxy for lensing. Foreground structures have a stronger effect on the lens potential than background structures, due to nonlinear effects in the foreground and downweighting in the background. There is dramatic variation in the net strength of ENV/LOS effects across different lens fields; modeling fields individually yields stronger priors for $H_0$ than ray tracing through $N$-body simulations. Models that ignore mass outside the lens yield poor fits and biased results. Adding external shear can account for tidal stretching from galaxies at redshifts $z \ge z_{\rm lens}$, but it requires corrections for external convergence and cannot reproduce nonlinear effects from foreground galaxies. Using the tidal approximation is reasonable for most perturbers as long as nonlinear redshift effects are included. Even then, the scatter in $H_0$ is limited by the lens profile degeneracy. Asymmetric image configurations produced by highly elliptical lens galaxies are less sensitive to the lens profile degeneracy, so they offer appealing targets for precision lensing analyses in future surveys like LSST and Euclid.

The arrival times, positions, and fluxes of multiple images in strong lens systems can be used to infer the presence of dark subhalos in the deflector, and thus test predictions of cold dark matter models. However, gravitational lensing does not distinguish between perturbations to a smooth gravitational potential arising from baryonic and non-baryonic mass. In this work, we quantify the extent to which the stellar mass distribution of a deflector can reproduce flux ratio and astrometric anomalies typically associated with the presence of a dark matter subhalo. Using Hubble Space Telescope images of nearby galaxies, we simulate strong lens systems with real distributions of stellar mass as they would be observed at redshift $z_d=0.5$. We add a dark matter halo and external shear to account for the smooth dark matter field, omitting dark substructure, and use a Monte Carlo procedure to characterize the distributions of image positions, time delays, and flux ratios for a compact background source of diameter 5 pc. By convolving high-resolution images of real galaxies with a Gaussian PSF, we simulate the most detailed smooth potential one could construct given high quality data, and find scatter in flux ratios of $\approx 10\%$, which we interpret as a typical deviation from a smooth potential caused by large and small scale structure in the lensing galaxy. We demonstrate that the flux ratio anomalies arising from galaxy-scale baryonic structure can be minimized by selecting the most massive and round deflectors, and by simultaneously modeling flux ratio and astrometric data.

In the standard structure formation scenario based on the cold dark matter paradigm, galactic halos are predicted to contain a large population of dark matter subhalos. While the most massive members of the subhalo population can appear as luminous satellites and be detected in optical surveys, establishing the existence of the low mass and mostly dark subhalos has proven to be a daunting task. Galaxy-scale strong gravitational lenses have been successfully used to study mass substructures lying close to lensed images of bright background sources. However, in typical galaxy-scale lenses, the strong lensing region only covers a small projected area of the lens's dark matter halo, implying that the vast majority of subhalos cannot be directly detected in lensing observations. In this paper, we point out that this large population of dark satellites can collectively affect gravitational lensing observables, hence possibly allowing their statistical detection. Focusing on the region of the galactic halo outside the strong lensing area, we compute from first principles the statistical properties of perturbations to the gravitational time delay and position of lensed images in the presence of a mass substructure population. We find that in the standard cosmological scenario, the statistics of these lensing observables are well approximated by Gaussian distributions. The formalism developed as part of this calculation is very general and can be applied to any halo geometry and choice of subhalo mass function. Our results significantly reduce the computational cost of including a large substructure population in lens models and enable the use of Bayesian inference techniques to detect and characterize the distributed satellite population of distant lens galaxies.

We present the spectroscopic redshift catalog from a wide-field survey of the fields of 28 galaxy-mass strong gravitational lenses. We discuss the acquisition and reduction of the survey data, collected over 40 nights of 6.5m MMT and Magellan time, employing four different multi-object spectrographs. We determine that no biases are introduced by combining datasets obtained with different instrument/spectrograph combinations. Special care is taken to determine redshift uncertainties using repeat observations. The redshift catalog consists of 9768 new and unique galaxy redshifts. 82.4% of the catalog redshifts are between z=0.1 and z=0.7, and the catalog median redshift is z=0.36. The data from this survey will be used to study the lens environments and line-of-sight structures to gain a better understanding of the effects of large scale structure on lens statistics and lens-derived parameters.

In strong gravitational lens systems, the light bending is usually dominated by one main galaxy, but may be affected by other mass along the line of sight (LOS). Shear and convergence can be used to approximate the contributions from less significant perturbers (e.g. those that are projected far from the lens or have a small mass), but higher order effects need to be included for objects that are closer or more massive. We develop a framework for multiplane lensing that can handle an arbitrary combination of tidal planes treated with shear and convergence and planes treated exactly (i.e., including higher order terms). This framework addresses all of the traditional lensing observables including image positions, fluxes, and time delays to facilitate lens modelling that includes the non-linear effects due to mass along the LOS. It balances accuracy (accounting for higher-order terms when necessary) with efficiency (compressing all other LOS effects into a set of matrices that can be calculated up front and cached for lens modelling). We identify a generalized multiplane mass sheet degeneracy, in which the effective shear and convergence are sums over the lensing planes with specific, redshift-dependent weighting factors.

Gravitational lens modeling of spatially resolved sources is a challenging inverse problem with many observational constraints and model parameters. We examine established pixel-based source reconstruction algorithms for de-lensing the source and constraining lens model parameters. Using test data for four canonical lens configurations, we explore statistical and systematic uncertainties associated with gridding, source regularisation, interpolation errors, noise, and telescope pointing. Specifically, we compare two gridding schemes in the source plane: a fully adaptive grid that follows the lens mapping but is irregular, and an adaptive Cartesian grid. We also consider regularisation schemes that minimise derivatives of the source (using two finite difference methods) and introduce a scheme that minimises deviations from an analytic source profile. Careful choice of gridding and regularisation can reduce "discreteness noise" in the $\chi^2$ surface that is inherent in the pixel-based methodology. With a gridded source, some degree of interpolation is unavoidable, and errors due to interpolation need to be taken into account (especially for high signal-to-noise data). Different realisations of the noise and telescope pointing lead to slightly different values for lens model parameters, and the scatter between different "observations" can be comparable to or larger than the model uncertainties themselves. The same effects create scatter in the lensing magnification at the level of a few percent for a peak signal-to-noise ratio of 10, which decreases as the data quality improves.

Certain configurations of massive structures projected along the line of sight maximize the number of detections of gravitationally lensed $z\sim10$ galaxies. We characterize such lines of sight with the \'etendue $\sigma_\mu$, the area in the source plane magnified over some threshold $\mu$. We use the Millennium I and Millennium XXL cosmological simulations to determine the frequency of high $\sigma_\mu$ beams on the sky, their properties, and efficient selection criteria. We define the best beams as having $\sigma_{\mu>3} >2000$ arcsec$^2$, for a $z\sim10$ source plane, and predict $477 \pm 21$ such beams on the sky. The total mass in the beam and $\sigma_{\mu>3}$ are strongly correlated. After controlling for total mass, we find a significant residual correlation between $\sigma_{\mu>3}$ and the number of cluster-scale halos ($>10^{14} M_\odot h^{-1}$) in the beam. Beams with $\sigma_{\mu>3} >2000$ arcsec$^2$, which should be best at lensing $z\sim10$ galaxies, are ten times more likely to contain multiple cluster-scale halos than a single cluster-scale halo. Beams containing an Abell 1689-like massive cluster halo often have additional structures along the line of sight, including at least one additional cluster-scale ($M_{200}>10^{14}M_\odot h^{-1}$) halo 28% of the time. Selecting beams with multiple, massive structures will lead to enhanced detection of the most distant and intrinsically faint galaxies.

Lines of sight with multiple, projected, cluster-scale halos have high total masses and complex lens plane interactions that can boost the area of magnification, or etendue, making detection of faint background sources more likely than elsewhere. To identify these new "compound" cosmic telescopes, we have found lines-of-sight with the highest integrated mass densities, as traced by the projected concentrations of Luminous Red Galaxies (LRGs). We use 1151 MMT Hectospec spectra to derive preliminary magnification maps for two such lines of sight with total mass exceeding ~ 3 x 10$^{15}$ Msun -- J0850+3604 (0850) and J1306+4632 (1306). We identify 2-3 group- and cluster-scale halos in each beam over 0.1 < z < 0.7, all of which are well-traced by LRGs. In Subaru Suprime-Cam imaging of beam 0850, we discover serendipitously a candidate multiply-imaged V-dropout source at z = 5.03, whose location is consistent with the critical curves for a source plane of $z_s$ = 5.03 predicted by our mass model. Incorporating the position of the candidate multiply-imaged galaxy as a constraint on the critical curve location in 0850 narrows the 68% confidence band on lens plane area with mu > 10 for a source plane of $z_s$ = 10 to [1.8, 4.2] square arcminutes, comparable to that of MACS 0717+3745 and El Gordo, two of the most powerful known single cluster lenses. The 68% confidence intervals on the lens plane area with mu > 10 for 1306 are [2.3, 6.7] square arcminutes. The significant lensing power of our beams makes them powerful probes of reionization and galaxy formation in the early Universe.

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.

The best gravitational lenses for detecting distant galaxies are those with the largest mass concentrations and the most advantageous configurations of that mass along the line of sight. Our new method for finding such gravitational telescopes uses optical data to identify projected concentrations of luminous red galaxies (LRGs). LRGs are biased tracers of the underlying mass distribution, so lines of sight with the highest total luminosity in LRGs are likely to contain the largest total mass. We apply this selection technique to the Sloan Digital Sky Survey and identify the 200 fields with the highest total LRG luminosities projected within a 3.5' radius over the redshift range 0.1 < z < 0.7. The redshift and angular distributions of LRGs in these fields trace the concentrations of non-LRG galaxies. These fields are diverse; 22.5% contain one known galaxy cluster and 56.0% contain multiple known clusters previously identified in the literature. Thus, our results confirm that these LRGs trace massive structures and that our selection technique identifies fields with large total masses. These fields contain 2-3 times higher total LRG luminosities than most known strong-lensing clusters and will be among the best gravitational lensing fields for the purpose of detecting the highest redshift galaxies.

We present the discovery of a Type Ia supernova (SN) at redshift $z = 1.914$ from the CANDELS multi-cycle treasury program on the \textit{Hubble Space Telescope (HST)}. This SN was discovered in the infrared using the Wide-Field Camera 3, and it is the highest-redshift Type Ia SN yet observed. We classify this object as a SN\,Ia by comparing its light curve and spectrum with those of a large sample of Type Ia and core-collapse supernovae (SNe). Its apparent magnitude is consistent with that expected from the $\Lambda$CDM concordance cosmology. We discuss the use of spectral evidence for classification of $z > 1.5$ SNe\,Ia using {\it HST} grism simulations, finding that spectral data alone can frequently rule out SNe\,II, but distinguishing between SNe\,Ia and SNe\,Ib/c can require prohibitively long exposures. In such cases, a quantitative analysis of the light curve may be necessary for classification. Our photometric and spectroscopic classification methods can aid the determination of SN rates and cosmological parameters from the full high-redshift CANDELS SN sample.

We investigate the gravitational lensing properties of lines of sight containing multiple cluster-scale halos, motivated by their ability to lens very high-redshift (z ~ 10) sources into detectability. We control for the total mass along the line of sight, isolating the effects of distributing the mass among multiple halos and of varying the physical properties of the halos. Our results show that multiple-halo lines of sight can increase the magnified source-plane region compared to the single cluster lenses typically targeted for lensing studies, and thus are generally better fields for detecting very high-redshift sources. The configurations that result in optimal lensing cross sections benefit from interactions between the lens potentials of the halos when they overlap somewhat on the sky, creating regions of high magnification in the source plane not present when the halos are considered individually. The effect of these interactions on the lensing cross section can even be comparable to changing the total mass of the lens from 10^15 M_sun to 3x10^15 M_sun. The gain in lensing cross section increases as the mass is split into more halos, provided that the lens potentials are projected close enough to interact with each other. A nonzero projected halo angular separation, equal halo mass ratio, and high projected halo concentration are the best mass configurations, whereas projected halo ellipticity, halo triaxiality, and the relative orientations of the halos are less important. Such high mass, multiple-halo lines of sight exist in the SDSS.

In this second paper, we develop an analytical theory of quasi-equatorial lensing by Kerr black holes. In this setting we solve perturbatively our general lens equation with displacement given in Paper I, going beyond weak-deflection Kerr lensing to third order in our expansion parameter epsilon, which is the ratio of the angular gravitational radius to the angular Einstein radius. We obtain new formulas and results for the bending angle, image positions, image magnifications, total unsigned magnification, and centroid, all to third order in epsilon and including the displacement. New results on the time delay between images are also given to second order in epsilon, again including displacement. For all lensing observables we show that the displacement begins to appear only at second order in epsilon. When there is no spin, we obtain new results on the lensing observables for Schwarzschild lensing with displacement.

We develop a unified, analytic framework for gravitational lensing by Kerr black holes. In this first paper we present a new, general lens equation and magnification formula governing lensing by a compact object. Our lens equation assumes that the source and observer are in the asymptotically flat region and does not require a small angle approximation. Furthermore, it takes into account the displacement that occurs when the light ray's tangent lines at the source and observer do not meet on the lens plane. We then explore our lens equation in the case when the compact object is a Kerr black hole. Specifically, we give an explicit expression for the displacement when the observer is in the equatorial plane of the Kerr black hole as well as for the case of spherical symmetry.

We investigate the properties of dark matter substructure in the gravitational lens HE 0435-1223 (z_l=0.455) via its effects on the positions and flux ratios of the quadruply-imaged background quasar (z_s=1.689). We start with a smooth mass model, add individual, truncated isothermal clumps near the lensed images, and use the Bayesian evidence to compare the quality of different models. Compared with smooth models, models with at least one clump near image A are strongly favored. The mass of this clump within its Einstein radius is log(Mein/Msun) = 7.65 +0.87/-0.84. The Bayesian evidence provides weaker support for a second clump near image B, with log(Mein/Msun) = 6.55 +1.01/-1.51. We also examine models with a full population of substructure, and find the mass fraction in substructure at the Einstein radius to be f_sub > 0.00077, assuming the total clump masses follow a mass function dN/dM proportional to M^(-1.9) over the range M = 10^7-10^10 Msun. Few-clump and population models produce similar Bayesian evidence values, so neither type of model is objectively favored.

Nested sampling has emerged as a valuable tool for Bayesian analysis, in particular for determining the Bayesian evidence. The method is based on a specific type of random sampling of the likelihood function and prior volume of the parameter space. I study the statistical uncertainty in the evidence computed with nested sampling. I examine the uncertainty estimator from Skilling (2004, 2006) and introduce a new estimator based on a detailed analysis of the statistical properties of nested sampling. Both perform well in test cases and make it possible to obtain the statistical uncertainty in the evidence with no additional computational cost.

We examine the wavelength dependence of flux ratios for six gravitationally lensed quasars using K and L' images obtained at the Gemini North 8m telescope. We select lenses with source redshifts z_s < 2.8 so that K-band images probe rest-frame optical emission from accretion disks, while L'-band images probe rest-frame near-infrared flux emitted (in part) from the more extended surrounding torus. Since the observations correspond to different source sizes, the K and L' flux ratios are sensitive to structure on different scales and may be useful for studying small-structure in the lens galaxies. Four of the six lenses show differences between K and L' flux ratios. In HE 0435$-1223, SDSS 0246-0825, and HE 2149-2745 the differences may be attributable to known microlensing and/or intrinsic variability. In SDSS 0806+2006 the wavelength dependence is not easily attributed to known variations, and may indicate the presence of substructure. By contrast, in Q0142-100 and SBS 0909+523 the K and L' flux ratios are consistent within the uncertainties. We discuss the utility of the current data for studying chromatic effects related to microlensing, dust extinction, and dark matter substructure.

Using new photometric and spectroscopic data in the fields of nine strong gravitational lenses that lie in galaxy groups, we analyze the effects of both the local group environment and line-of-sight galaxies on the lens potential. We use Monte Carlo simulations to derive the shear directly from measurements of the complex lens environment, providing the first detailed independent check of the shear obtained from lens modeling. We account for possible tidal stripping of the group galaxies by varying the fraction of total mass apportioned between the group dark matter halo and individual group galaxies. The environment produces an average shear of gamma = 0.08 (ranging from 0.02 to 0.17), significant enough to affect quantities derived from lens observables. However, the direction and magnitude of the shears do not match those obtained from lens modeling in three of the six 4-image systems in our sample (B1422, RXJ1131, and WFI2033). The source of this disagreement is not clear, implying that the assumptions inherent in both the environment and lens model approaches must be reconsidered. If only the local group environment of the lens is included, the average shear is gamma = 0.05 (ranging from 0.01 to 0.14), indicating that line-of-sight contributions to the lens potential are not negligible. We isolate the effects of various theoretical and observational uncertainties on our results. Of those uncertainties, the scatter in the Faber-Jackson relation and error in the group centroid position dominate. Future surveys of lens environments should prioritize spectroscopic sampling of both the local lens environment and objects along the line of sight, particularly those bright (I < 21.5) galaxies projected within 5' of the lens.

We present a detailed strong lensing analysis of an HST/ACS legacy dataset for the first gravitational lens, Q0957+561. With deep imaging we identify 24 new strongly lensed features, which we use to constrain mass models. We model the stellar component of the lens galaxy using the observed luminosity distribution, and the dark matter halo using several different density profiles. We draw on the weak lensing analysis by Nakajima et al. (2009) to constrain the mass sheet and environmental terms in the lens potential. Adopting the well-measured time delay, we find H_0 = 85 (+14/-13) km/s/Mpc (68% CL) using lensing constraints alone. The principal uncertainties in H_0 are tied to the stellar mass-to-light ratio (a variant of the radial profile degeneracy in lens models). Adding constraints from stellar population synthesis models, we obtain H_0 = 79.3 (+6.7/-8.5) km/s/Mpc (68% CL). We infer that the lens galaxy has a rising rotation curve and a dark matter distribution with an inner core. Intriguingly, we find the quasar flux ratios predicted by our models to be inconsistent with existing radio measurements, suggesting the presence of substructure in the lens.

We examine the ability of gravitational lens time delays to reveal complex structure in lens potentials. In Congdon, Keeton & Nordgren (2008), we predicted how the time delay between the bright pair of images in a "fold" lens scales with the image separation, for smooth lens potentials. Here we show that the proportionality constant increases with the quadrupole moment of the lens potential, and depends only weakly on the position of the source along the caustic. We use Monte Carlo simulations to determine the range of time delays that can be produced by realistic smooth lens models consisting of isothermal ellipsoid galaxies with tidal shear. We can then identify outliers as "time delay anomalies". We find evidence for anomalies in close image pairs in the cusp lenses RX J1131$-$1231 and B1422+231. The anomalies in RX J1131$-$1231 provide strong evidence for substructure in the lens potential, while at this point the apparent anomalies in B1422+231 mainly indicate that the time delay measurements need to be improved. We also find evidence for time delay anomalies in larger-separation image pairs in the fold lenses, B1608+656 and WFI 2033$-$4723, and the cusp lens RX J0911+0551. We suggest that these anomalies are caused by some combination of substructure and a complex lens environment. Finally, to assist future monitoring campaigns we use our smooth models with shear to predict the time delays for all known four-image lenses.

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.

Motivated by debris disk studies, we investigate the gravitational microlensing of background starlight by a planetesimal disk around a foreground star. We use dynamical survival models to construct a plausible example of a planetesimal disk and study its microlensing properties using established ideas of microlensing by small bodies. When a solar-type source star passes behind a planetesimal disk, the microlensing light curve may exhibit short-term, low-amplitude residuals caused by planetesimals several orders of magnitude below Earth mass. The minimum planetesimal mass probed depends on the photometric sensitivity and the size of the source star, and is lower when the planetesimal lens is located closer to us. Planetesimal lenses may be found more nearby than stellar lenses because the steepness of the planetesimal mass distribution changes how the microlensing signal depends on the lens/source distance ratio. Microlensing searches for planetesimals require essentially continuous monitoring programs that are already feasible and can potentially set constraints on models of debris disks, the supposed extrasolar analogues of Kuiper belts.

Mass clumps in gravitational lens galaxies can perturb lensed images in characteristic ways. Strong lens flux ratios have been used to constrain the amount of dark matter substructure in lens galaxies, and various other observables have been considered as additional probes of substructure. We study the general theory of lensing with stochastic substructure in order to understand how lensing observables depend on the mass function and spatial distribution of clumps. We find that magnification perturbations are mainly sensitive to the total mass in substructure projected near the lensed images; when the source is small, flux ratios are not very sensitive to the shape of the clump mass function. Position perturbations are mainly sensitive to a characteristic clump mass scale, namely m_eff = <m^2>/<m>, with some mild dependence on other mass moments when the spatial distribution is not uniform. They have contributions from both "local" and "global" populations of clumps (i.e., those projected near the images, and those farther away). Time delay perturbations are sensitive to the same characteristic mass, m_eff, and mainly driven by the global population of clumps. While there is significant scatter in all lensing quantities, there are some non-trivial correlations that may contain further information about the clump population. Our results indicate that a joint analysis of multiple lens observables will offer qualitatively new constraints on the mass function and spatial distribution of dark matter substructure in distant galaxies.

[Abridged] Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio, and image separation distribution. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sersic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population.

Studies of strong gravitational lensing in current and upcoming wide and deep photometric surveys, and of stellar kinematics from (integral-field) spectroscopy at increasing redshifts, promise to provide valuable constraints on galaxy density profiles and shapes. However, both methods are affected by various selection and modelling biases, whch we aim to investigate in a consistent way. In this first paper in a series we develop a flexible but efficient pipeline to simulate lensing by realistic galaxy models. These galaxy models have separate stellar and dark matter components, each with a range of density profiles and shapes representative of early-type, central galaxies without significant contributions from other nearby galaxies. We use Fourier methods to calculate the lensing properties of galaxies with arbitrary surface density distributions, and Monte Carlo methods to compute lensing statistics such as point-source lensing cross-sections. Incorporating a variety of magnification bias modes lets us examine different survey limitations in image resolution and flux. We rigorously test the numerical methods for systematic errors and sensitivity to basic assumptions. We also determine the minimum number of viewing angles that must be sampled in order to recover accurate orientation-averaged lensing quantities. We find that for a range of non-isothermal stellar and dark matter density profiles typical of elliptical galaxies, the combined density profile and corresponding lensing properties are surprisingly close to isothermal around the Einstein radius. The converse implication is that constraints from strong lensing and/or stellar kinematics, which are indeed consistent with isothermal models near the Einstein radius, cannot trivially be extrapolated to smaller and larger radii.

We show that dark matter substructure in galaxy-scale halos perturbs the time delays between images in strong gravitational lens systems. The variance of the effect depends on the subhalo mass function, scaling as the product of the substructure mass fraction and a characteristic mass of subhalos (namely <m^2>/<m>). Time delay perturbations therefore complement gravitational lens flux ratio anomalies and astrometric perturbations by measuring a different moment of the subhalo mass function. Unlike flux ratio anomalies, "time delay millilensing" is unaffected by dust extinction or stellar microlensing in the lens galaxy. Furthermore, we show that time delay ratios are immune to the radial profile degeneracy that usually plagues lens modeling. We lay out a mathematical theory of time delay perturbations and find it to be tractable and attractive. We predict that in "cusp" lenses with close triplets of images, substructure may change the arrival-time order of the images (compared with smooth models). We discuss the possibility that this effect has already been observed in RX J1131-1231.

In strong gravitational lensing, the multiple images we see correspond to light rays that leave the source in slightly different directions. If the source emission is anisotropic, the images may differ from conventional lensing predictions (which assume isotropy). To identify scales on which source anisotropy may be important, we study the angle delta between the light rays emerging from the source, for different lensing configurations. If the lens has a power law profile M propto R^gamma, the angle delta initially increases with lens redshift and then either diverges (for a steep profile gamma<1), remains constant (for an isothermal profile gamma=1), or vanishes (for a shallow profile gamma>1) as zl approaches zs. The scaling with lens mass is roughly delta propto M^(1/(2-gamma)). The results for an NFW profile are qualitatively similar to those for a shallow power law, with delta peaking at about half the redshift of the source (not half the distance). In practice, beaming could modify the statistics of beamed sources lensed by massive clusters: for an opening angle theta_jet, there is a probability as high as P ~ 0.02-0.07 (theta_jet/{0.5} deg)^-1 that one of the lensed images may be missed (for 2 < zs < 6). Differential absorption within Active Galactic Nuclei could modify the flux ratios of AGNs lensed by clusters; a sample of AGNs lensed by clusters could provide further constraints on the sizes of absorbing regions. Source anisotropy is not likely to be a significant effect in galaxy-scale strong lensing.

Attempts to constrain the Hubble constant using the strong gravitational lens system Q0957+561 are limited by systematic uncertainties in the mass model, since the time delay is known very precisely. One important systematic effect is the mass sheet degeneracy, which arises because strong lens modeling cannot constrain the presence or absence of a uniform mass sheet $\kappa$, which rescales $H_0$ by the factor $(1-\kappa)$. In this paper we present new constraints on the mass sheet derived from a weak-lensing analysis of the Hubble Space Telescope imaging of a 6 arcmin square region surrounding the lensed quasar. The average mass sheet within a circular aperture (the strong lens model region) is constrained by integrating the tangential weak gravitational shear over the surrounding area. We find the average convergence within a $30"$ radius around the lens galaxy to be $\kappa(<30") = 0.166\pm0.056$ ($1\sigma$ confidence level), normalized to the quasar redshift. This includes contributions from both the lens galaxy and the surrounding cluster. We also constrain a few other low-order terms in the lens potential by applying a multipole aperture mass formalism to the gravitational shear in an annulus around the strong lensing region. Implications for strong lens models and the Hubble constant are discussed in an accompanying paper.

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.

Gravitational lensing provides a unique and powerful probe of the mass distributions of distant galaxies. Four-image lens systems with fold and cusp configurations have two or three bright images near a critical point. Within the framework of singularity theory, we derive analytic relations that are satisfied for a light source that lies a small but finite distance from the astroid caustic of a four-image lens. Using a perturbative expansion of the image positions, we show that the time delay between the close pair of images in a fold lens scales with the cube of the image separation, with a constant of proportionality that depends on a particular third derivative of the lens potential. We also apply our formalism to cusp lenses, where we develop perturbative expressions for the image positions, magnifications and time delays of the images in a cusp triplet. Some of these results were derived previously for a source asymptotically close to a cusp point, but using a simplified form of the lens equation whose validity may be in doubt for sources that lie at astrophysically relevant distances from the caustic. Along with the work of Keeton et al. (2005), this paper demonstrates that perturbation theory plays an important role in theoretical lensing studies.

We present evidence that the optically unidentified radio source, FIRST J121839.7+295325, may be strongly lensing a background galaxy. We estimate the redshift of the assumed gravitational arc, discovered in parallel imaging with HST, from MMT-Blue Channel spectroscopy to be z_{arc}=2.48_-0.05^+0.14. We present lens models with an Einstein radius of R_E=1.3" which contains a mass of M_{dyn}=10^{12 +- 0.5} M_{sol}, where the uncertainty reflects the range of possible lens redshifts. The putative lens is not detected to J_{lim}=22.0 mag and H_{lim}=20.7 mag in our MMT-SWIRC imaging. Using the flux limits from WFPC2 and SWIRC, we estimate that the dynamical mass-to-light ratio of J121839.7+295325 is M_{dyn}/L_B >~ 10 M_sol/L_sol for A_V=1 mag, and this lower limit could be as high as 30 M_sol/L_sol for A_V=0 mag. Since the radio source is optically unidentified (V_{lim}=25.5 mag) and has a radio flux of S_{1.4 GHz}=33 mJy, it is likely a massive early-type galaxy which hosts a radio-loud AGN at 0.8<z<1.5. However, the present data cannot uniquely determine the mass-to-light ratio of the lensing galaxy, and hence the possibility that this system may be a reasonably dark lens is not ruled out.

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

Spectroscopic galaxy-galaxy lens searches are presently the most prolific method of identifying strong lens systems in large data sets. We study the probabilities associated with these lens searches, namely the probability of identifying a candidate with rogue [OII] emission lines in a galaxy's spectrum, and the probability that the candidate will show features of strong lensing in follow-up photometric observations. We include selection effects unique to spectroscopic data, and apply them to the Sloan Lens ACS (SLACS) survey (Bolton et al. 2006). The most significant selection effect is the finite size of the spectroscopic fiber which selects against large separation lenses and results in a non-monotonic dependence of the rogue line probability on velocity dispersion. For example, with the 3 arcsec diameter SDSS fiber and 2 arcsec FWHM seeing, we find that the probability that a given LRG has a rogue [OII] line in its spectrum decreases with velocity dispersion from 150 km/s to 300 km/s and then increases up to 400 km/s for a given source size. The total probability for observing a rogue line in a single survey spectrum is ~0.9-3.0%, and the total lensing rate is ~0.5-1.3%. The range is due to uncertainties in the physical size of [OII] emission regions, and in the evolution of the [OII] luminosity function. Our estimates are a factor of ~5 higher than the results of the SLACS survey, a discrepancy which we attribute to the SLACS requirement that multiple rogue lines be observed simultaneously.

We report the first results of our systematic search for strongly lensed quasars using the spectroscopically confirmed quasars in the Sloan Digital Sky Survey (SDSS). Among 46,420 quasars from the SDSS Data Release 3 (~4188 deg^2), we select a subsample of 22,683 quasars that are located at redshifts between 0.6 and 2.2 and are brighter than the Galactic extinction corrected i-band magnitude of 19.1. We identify 220 lens candidates from the quasar subsample, for which we conduct extensive and systematic follow-up observations in optical and near-infrared wavebands, in order to construct a complete lensed quasar sample at image separations between 1'' and 20'' and flux ratios of faint to bright lensed images larger than 10^{-0.5}. We construct a statistical sample of 11 lensed quasars. Ten of these are galaxy-scale lenses with small image separations (~1''-2'') and one is a large separation (15'') system which is produced by a massive cluster of galaxies, representing the first statistical sample of lensed quasars including both galaxy- and cluster-scale lenses. The Data Release 3 spectroscopic quasars contain an additional 11 lensed quasars outside the statistical sample.

We present cosmological results from the statistics of lensed quasars in the Sloan Digital Sky Survey (SDSS) Quasar Lens Search. By taking proper account of the selection function, we compute the expected number of quasars lensed by early-type galaxies and their image separation distribution assuming a flat universe, which is then compared with 7 lenses found in the SDSS Data Release 3 to derive constraints on dark energy under strictly controlled criteria. For a cosmological constant model (w=-1) we obtain \Omega_\Lambda=0.74^{+0.11}_{-0.15}(stat.)^{+0.13}_{-0.06}(syst.). Allowing w to be a free parameter we find \Omega_M=0.26^{+0.07}_{-0.06}(stat.)^{+0.03}_{-0.05}(syst.) and w=-1.1\pm0.6(stat.)^{+0.3}_{-0.5}(syst.) when combined with the constraint from the measurement of baryon acoustic oscillations in the SDSS luminous red galaxy sample. Our results are in good agreement with earlier lensing constraints obtained using radio lenses, and provide additional confirmation of the presence of dark energy consistent with a cosmological constant, derived independently of type Ia supernovae.

We model the massive dark object at the center of the Galaxy as a Schwarzschild black hole as well as Janis-Newman-Winicour naked singularities, characterized by the mass and scalar charge parameters, and study gravitational lensing (particularly time delay, magnification centroid, and total magnification) by them. We find that the lensing features are qualitatively similar (though quantitatively different) for the Schwarzschild black holes, weakly naked, and marginally strongly naked singularities. However, the lensing characteristics of strongly naked singularities are qualitatively very different from those due the Schwarzschild black holes. The images produced by Schwarzschild black hole lenses and weakly naked and marginally strongly naked singularity lenses always have positive time delays. On the other hand, the strongly naked singularity lenses can give rise to images with positive, zero, or negative time delays. In particular, for a large angular source position the direct image (the outermost image on the same side as the source) due to strongly naked singularity lensing always has negative time delay. We also found that the scalar field decreases the time delay and increases the magnitude of magnifications of images; this result could have important implications for cosmology. As the Janis-Newman-Winicour metric also describes the exterior gravitational field of a scalar star, naked singularities as well as scalar star lenses, if these exist in nature, will serve as more efficient cosmic telescopes than regular gravitational lenses.

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.

Microlensing promises to be a powerful tool for studying distant galaxies and quasars. As the data and models improve, there are systematic effects that need to be explored. Quasar continuum and broad-line regions may respond differently to microlensing due to their different sizes; to understand this effect, we study microlensing of finite sources by a mass function of stars. We find that microlensing is insensitive to the slope of the mass function but does depend on the mass range. For negative parity images, diluting the stellar population with dark matter increases the magnification dispersion for small sources and decreases it for large sources. This implies that the quasar continuum and broad-line regions may experience very different microlensing in negative-parity lensed images. We confirm earlier conclusions that the surface brightness profile and geometry of the source have little effect on microlensing. Finally, we consider non-circular sources. We show that elliptical sources that are aligned with the direction of shear have larger magnification dispersions than sources with perpendicular alignment, an effect that becomes more prominent as the ellipticity increases. Elongated sources can lead to more rapid variability than circular sources, which raises the prospect of using microlensing to probe source shape.

We study microlensing of the faint images that form close to the centers of strong gravitational lens galaxies. These central images, which have finally begun to yield to observations, naturally appear in dense stellar fields and may be particularly sensitive to fine granularity in the mass distribution. The microlensing magnification maps for overfocussed (i.e., demagnified) images differ strikingly from those for magnified images. In particular, the familiar "fold" and "cusp" features of maps for magnified images are only present for certain values of the fraction, f, of the surface mass density contained in stars. For central images, the dispersion in microlensing magnifications is generally larger than for normal (minimum and saddle) images, especially when the source is comparable to or larger than the stellar Einstein radius. The dispersion depends in a complicated way on f; this behaviour may hold the key to using microlensing as a probe of the relative densities of stars and dark matter in the cores of distant galaxies. Quantitatively, we predict that the central image C in PMN J1632-0033 has a magnification dispersion of 0.6 magnitudes for Rsrc/Rein <~ 1, or 0.3 mag for Rsrc/Rein = 10. For comparison, the dispersions are 0.5-0.6 mag for image B and 0.05-0.1 mag for image A, if Rsrc/Rein <~ 1; and just 0.1 mag for B and 0.008 mag for A if Rsrc/Rein = 10. (The dispersions can be extrapolated to larger sources sizes as sigma \propto Rsrc^{-1}.) Thus, central images are more susceptible than other lensed images to microlensing and hence good probes for measuring source sizes.

Using ray tracing for simple analytic profiles, we demonstrate that the lensing cross section for producing giant arcs has distinct contributions due to arcs formed through image distortion only, and arcs form from the merging of two or three images. We investigate the dependence of each of these contributions on halo ellipticity and on the slope of the density profile, and demonstrate that at fixed Einstein radius, the lensing cross section increases as the halo profile becomes steeper. We then compare simulations with and without baryonic cooling of the same cluster for a sample of six clusters, and demonstrate that cooling can increase the overall abundance of giant arcs by factors of a few. The net boost to the lensing probability for individual clusters is mass dependent, and can lower the effective low mass limit of lensing clusters. This last effect can potentially increase the number of lensing clusters by an extra 50%. While the magnitude of these effects may be overestimated due to the well known overcooling problem in simulations, it is evident that baryonic cooling has a non-negligible impact on the expected abundance of giant arcs, and hence cosmological constraints from giant arc abundances may be subject to large systematic errors.

Given the number of recently discovered galaxy-galaxy lens systems, we anticipate that a gravitationally lensed supernova will be observed within the next few years. We explore the possibility that stars in the lens galaxy will produce observable microlensing fluctuations in lensed supernova light curves. For typical parameters, we predict that ~70% of lensed SNe will show microlensing fluctuations > 0.5 mag, while ~25% will have fluctuations > 1 mag. Thus microlensing of lensed supernova will be both ubiquitous and observable. Additionally, we show that microlensing fluctuations will complicate measurements of time delays from multiply imaged supernovae: time delays accurate to better than a few days will be difficult to obtain. We also consider prospects for extracting the lens galaxy's stellar mass fraction and mass function from microlensing fluctuations via a new statistical measure, the time-weighted light curve derivative.

We present 426 epochs of optical monitoring data spanning 1000 days from December 2003 to June 2006 for the gravitationally lensed quasar SDSS J1004+4112. The time delay between the A and B images is 38.4+/-2.0 days in the expected sense that B leads A and the overall time ordering is C-B-A-D-E. The measured delay invalidates all published models. The models failed because they neglected the perturbations from cluster member galaxies. Models including the galaxies can fit the data well, but strong conclusions about the cluster mass distribution should await the measurement of the longer, and less substructure sensitive, delays of the C and D images. For these images, a CB delay of 681+/-15 days is plausible but requires confirmation, while CB and AD delays of >560 days and > 800 days are required. We clearly detect microlensing of the A/B images, with the delay-corrected flux ratios changing from B-A=0.44+/-0.01 mag in the first season to 0.29+/-0.01 mag in the second season and 0.32+/-0.01 mag in the third season.

In the present-day universe, supermassive black hole masses (MBH) appear to be strongly correlated with their galaxy's bulge luminosity, among other properties. In this study, we explore the analogous relationship between MBH, derived using the virial method, and the stellar R-band bulge luminosity (Lr) or stellar bulge mass (M*) at epochs of 1 < z < 4.5 using a sample of 31 gravitationally lensed AGNs and 20 non-lensed AGNs. At redshifts z > 1.7 (10--12 Gyrs ago), we find that the observed MBH--Lr relation is nearly the same (to within ~0.3 mag) as it is today. When the observed Lr are corrected for luminosity evolution, this means that the black holes grew in mass faster than their hosts, with the MBH/M* mass ratio being a factor of > 4(+2)(-1) times larger at z > 1.7 than it is today. By the redshift range 1<z<1.7 (8-10 Gyrs ago), the MBH/M* ratio is at most two times higher than today, but it may be consistent with no evolution. Combining the results, we conclude that the ratio MBH/M* rises with look-back time, although it may saturate at ~6 times the local value. Scenarios in which moderately luminous quasar hosts at z>1.7 were fully formed bulges that passively faded to the present epoch are ruled out.

Braneworld gravity is a model that endows physical space with an extra dimension. In the type II Randall-Sundrum braneworld gravity model, the extra dimension modifies the spacetime geometry around black holes, and changes predictions for the formation and survival of primordial black holes. We develop a comprehensive analytical formalism for far-field black hole lensing in this model, using invariant quantities to compute all geometric optics lensing observables. We then make the first analysis of wave optics in braneworld lensing, working in the semi-classical limit. We show that wave optics offers the only realistic way to observe braneworld effects in black hole lensing. We point out that if primordial braneworld black holes exist, have mass M, and contribute a fraction f of the dark matter, then roughly 3e5 x f (M/1e-18 Msun)^(-1) of them lie within our Solar System. These objects, which we call "attolenses," would produce interference fringes in the energy spectra of gamma-ray bursts at energies ~100 (M/1e-18 Msun)^(-1) MeV (which will soon be accessible with the GLAST satellite). Primordial braneworld black holes spread throughout the universe could produce similar interference effects; the probability for "attolensing" may be non-negligible. If interference fringes were observed, the fringe spacing would yield a simple upper limit on M. Detection of a primordial black hole with M <~ 1e-19 Msun would challenge general relativity and favor the braneworld model. Further work on lensing tests of braneworld gravity must proceed into the physical optics regime, which awaits a description of the full spacetime geometry around braneworld black holes.

We present results from Chandra observations of SDSS J1004+4112, a strongly lensed quasar system with a maximum image separation of 15". All four bright images of the quasar, as well as resolved X-ray emission originating from the lensing cluster, are clearly detected. The emission from the lensing cluster extends out to approximately 1.5 arcmin. We measure the bolometric X-ray luminosity and temperature of the lensing cluster to be 4.7e44 erg s^-1 and 6.4 keV, consistent with the luminosity-temperature relation for distant clusters. The mass estimated from the X-ray observation shows excellent agreement with the mass derived from gravitational lensing. The X-ray flux ratios of the quasar images differ markedly from the optical flux ratios, and the combined X-ray spectrum of the images possesses an unusually strong Fe Kalpha emission line, both of which are indicative of microlensing.

Many strong gravitational lenses lie in complex environments, such as poor groups of galaxies, that significantly bias conclusions from lens analyses. We are undertaking a photometric survey of all known galaxy-mass strong lenses to characterize their environments and include them in careful lens modeling, and to build a large, uniform sample of galaxy groups at intermediate redshifts for evolutionary studies. In this paper we present wide-field photometry of the environments of twelve lens systems with 0.24 < z_lens < 0.5. Using a red-sequence identifying technique, we find that eight of the twelve lenses lie in groups, and that ten group-like structures are projected along the line of sight towards seven of these lenses. Follow-up spectroscopy of a subset of these fields confirms these results. For lenses in groups, the group centroid position is consistent with the direction of the external tidal shear required by lens models. Lens galaxies are not all super-L_* ellipticals; the median lens luminosity is < L_*, and the distribution of lens luminosities extends 3 magnitudes below L_* (in agreement with theoretical models). Only two of the lenses in groups are the brightest group galaxy, in qualitative agreement with theoretical predictions. As in the local Universe, the highest velocity-dispersion groups contain a brightest member spatially coincident with the group centroid, whereas lower-dispersion groups tend to have an offset brightest group galaxy. This suggests that higher-dispersion groups are more dynamically relaxed than lower-dispersion groups and that at least some evolved groups exist by z ~ 0.5.

Gravitational lensing assists in the detection of quasar hosts by amplifying and distorting the host light away from the unresolved quasar core images. We present the results of HST observations of 30 quasar hosts at redshifts 1 < z < 4.5. The hosts are small in size (r_e <~ 6 kpc), and span a range of morphologies consistent with early-types (though smaller in mass) to disky/late-type. The ratio of the black hole mass (MBH, from the virial technique) to the bulge mass (M_bulge, from the stellar luminosity) at 1<z<1.7 is broadly consistent with the local value; while MBH/M_bulge at z>1.7 is a factor of 3--6 higher than the local value. But, depending on the stellar content the ratio may decline at z>4 (if E/S0-like), flatten off to 6--10 times the local value (if Sbc-like), or continue to rise (if Im-like). We infer that galaxy bulge masses must have grown by a factor of 3--6 over the redshift range 3>z>1, and then changed little since z~1. This suggests that the peak epoch of galaxy formation for massive galaxies is above z~1. We also estimate the duty cycle of luminous AGNs at z>1 to be ~1%, or 10^7 yrs, with sizable scatter.

We study gravitational lensing by compact objects in gravity theories that can be written in a Post-Post-Newtonian (PPN) framework: i.e., the metric is static and spherically symmetric, and can be written as a Taylor series in m/r, where m is the gravitational radius of the compact object. Working invariantly, we compute corrections to standard weak-deflection lensing observables at first and second order in the ratio of the angular gravitational radius to the angular Einstein ring radius of the lens. We show that the first-order corrections to the total magnification and centroid position vanish universally for gravity theories that can be written in the PPN framework. This arises from some surprising, fundamental relations among the lensing observables in PPN gravity models. We derive these relations for the image positions, magnifications, and time delays. A deep consequence is that any violation of the universal relations would signal the need for a gravity model outside the PPN framework (provided that some basic assumptions hold). In practical terms, the relations will guide observational programs to test general relativity, modified gravity theories, and possibly the Cosmic Censorship conjecture. We use the new relations to identify lensing observables that are accessible to current or near-future technology, and to find combinations of observables that are most useful for probing the spacetime metric. We give explicit applications to the Galactic black hole, microlensing, and the binary pulsar J0737-3039.

(Abridged) We present the first results from our spectroscopic survey of the environments of strong gravitational lenses. The lens galaxy belongs to a poor group of galaxies in six of the eight systems in our sample. We discover three new groups associated with the lens galaxies of BRI 0952-0115 (five members), MG 1654+1346 (seven members), and B2114+022 (five members). We more than double the number of members for another three previously known groups around the lenses MG 0751+2716 (13 total members), PG 1115+080 (13 total members), and B1422+231 (16 total members). We determine the kinematics of the six groups, including their mean velocities, velocity dispersions, and projected spatial centroids. The velocity dispersions of the groups range from 110 +170, -80 to 470 +100, -90 km/s. In at least three of the lenses -- MG0751, PG1115, and B1422 -- the group environment significantly affects the lens potential. These lenses happen to be the quadruply-imaged ones in our sample, which suggests a connection between image configuration and environment. The lens galaxy is the brightest member in fewer than half of the groups. Our survey also allows us to assess for the first time whether mass structures along the line of sight are important for lensing. We first show that, in principle, the lens potential may be affected by line-of-sight structures over a wide range of spatial and redshift offsets from the lens. We then quantify real line-of-sight effects using our survey and find that at least four of the eight lens fields have substantial interloping structures close in projection to the lens, and at least one of those structures (in the field of MG0751) significantly affects the lens potential.

We are developing a general, unified, and rigorous analytical framework for using gravitational lensing by compact objects to test different theories of gravity beyond the weak-deflection limit. In this paper we present the formalism for computing corrections to lensing observables for static, spherically symmetric gravity theories in which the corrections to the weak-deflection limit can be expanded as a Taylor series in one parameter, namely the gravitational radius of the lens object. We take care to derive coordinate-independent expressions and compute quantities that are directly observable. We compute first- and second-order corrections to the image positions, magnifications, and time delays. Interestingly, we find that the first-order corrections to the total magnification and centroid position vanish in all gravity theories that agree with general relativity in the weak-deflection limit, but they can remain nonzero in modified theories that disagree with GR in the weak-deflection limit. For the Reissner-Nordstrom metric and a related metric from heterotic string theory, our formalism reveals an intriguing connection between lensing observables and the condition for having a naked singularity, which could provide an observational method for testing the existence of such objects. We apply our formalism to the Galactic black hole and predict that the corrections to the image positions are at the level of 10 micro-arcseconds, while the correction to the time delay is a few hundredths of a second. These corrections would be measurable today if a pulsar were found to be lensed by the Galactic black hole; and they should be readily detectable with planned missions like MAXIM.

SDSS J0924+0219 is the most glaring example of a gravitational lens with anomalous flux ratios: optical broad-band photometry shows image D to be a factor of 12 fainter than expected for smooth lens potentials. We report spectroscopy showing that the anomaly is present in the broad emission line flux ratios as well. There are differences between the emission line and continuum flux ratios: the A/D ratio is 10 in the broad Lyman-alpha line and 19 in the associated continuum. Known variability argues for the presence of microlensing. We show that microlensing can account for both the continuum and emission line flux ratios, if the broad emission line region is comparable in size to the Einstein radii of the microlenses. Specifically, we need the half-light radius of the broad-line region to be R_BLR <~ 0.4 Rein ~ 9$ lt-days, which is small but reasonable. If the broad-line region is that large, then stars can contribute only 15-20% of the surface mass density at the positions of the images. While we cannot exclude the possibility that millilensing by dark matter substructure is present as well, we conclude that microlensing is present and sufficient to explain existing data. Under this hypothesis, the A/D flux ratio should return to a value close to unity on a time scale of years rather than millennia.

We study the impact of lens galaxy environments on the image separation distribution of lensed quasars. We account for both environmental convergence and shear, using a joint distribution derived from galaxy formation models calibrated by galaxy-galaxy lensing data and number counts of massive elliptical galaxies. We find that the external field enhances lensing probabilities, particularly at large image separations; the increase is ~30% at \theta=3'' and ~200% at \theta=5'', when we adopt a power-law source luminosity function \Phi(L) \propto L^-2.1. The enhancement is mainly driven by convergence, which boosts both the image separation and magnification bias (for a fixed lens galaxy mass). These effects have been neglected in previous studies of lens statistics. Turning the problem around, we derive the posterior convergence and shear distributions and point out that they are strong functions of image separation; lens systems with larger image separations are more likely to lie in dense environments.

Anomalous radio flux ratios in four-image gravitational lens systems can be explained by adding a clumpy cold dark matter (CDM) component to the potential of the lens galaxy. Evans & Witt (2003) recently suggested that smooth multipole perturbations can provide a reasonable alternative to CDM substructure in some but not all cases. We generalize their method in two ways so as to determine whether multipole models can explain highly anomalous systems. We carry the multipole expansion to higher order, and also include external tidal shear as a free parameter. Fitting for the shear proves crucial to finding a physical (positive-definite density) model. For B1422+231, working to order k = 5 (and including shear) yields a model that is physical but implausible. Going to higher order (k >~ 9) reduces global departures from ellipticity at the cost of introducing small scale wiggles near the bright images. These localized undulations are more pronounced in B2045+265, where k ~ 17 multipoles are required to smooth out large scale deviations from elliptical symmetry. Such modes surely cannot be taken at face value; they must indicate that the models are trying to reproduce some other sort of structure. Our formalism naturally finds models that fit the data exactly, but we use B0712+472 to show that measurement uncertainties have little effect on our results. Finally, we consider the system B1933+503, where two sources are lensed by the same foreground galaxy. The additional constraints provided by the images of the second source render the multipole model unphysical. We conclude that external shear must be taken into account to obtain plausible models, and that a purely smooth angular structure for the lens galaxy does not provide a viable alternative to the prevailing CDM clump hypothesis.

When the source in a 4-image gravitational lens lies close to a "fold" caustic, two of the lensed images lie close together. If the lens potential is smooth on the scale of the separation between the two close images, then the difference between their fluxes should approximately vanish, Rfold = (F1-F2)/(F1+F2) \approx 0. Violations of this "fold relation" in observed lenses are thought to indicate the presence of structure on scales smaller than the separation between the close images. We study the fold relation and find it to be more subtle and rich than was previously realized. The degree to which Rfold can differ from zero for realistic smooth lenses depends not only on the distance of the source from the caustic, but also on its location along the caustic, and on the angular structure of the lens potential. Thus, it is impossible to say from Rfold alone whether observed flux ratios are anomalous or not. Instead, we must consider the full distribution of Rfold values that can be obtained from smooth potentials that reproduce the separation between the two close images and the distance to the next nearest image. We analyze the generic and specific features of this distribution, and then show that 5 of the 12 known lenses with fold configurations have flux ratio anomalies: B0712+472, SDSS 0924+0219, PG 1115+080, B1555+375, and B1933+503. Combining this with our previous analysis revealing anomalies in 3 of the 4 known lenses with cusp configurations, we conclude that at least half (8/16) of all 4-image lenses that admit generic, local analyses exhibit flux ratio anomalies indicative of small-scale structure. The fold and cusp relations do not reveal the nature of the implied small-scale structure, but do provide the formal foundation for substructure studies. [Abridged]

The supermassive black hole at the center of a distant galaxy can be weighed, in rare but realistic cases, when the galaxy acts as a strong gravitational lens. The central image that should be produced by the lens is either destroyed or accompanied by a second central image, depending on the mass of the black hole. We demonstrate that when a central image pair is detected, the mass of the black hole can be determined with an accuracy of < 0.1 dex, if the form of the smooth mass distribution near the galaxy core is known. Uncertainty in the central mass distribution introduces a systematic error in the black hole mass measurement. However, even with nearly complete ignorance of the inner mass distribution, the black hole mass can still be determined to within a factor of 10. Central image pairs should be readily observable with future radio interferometers, allowing this technique to be used for a census of supermassive black holes in inactive galaxies at significant redshift (0.2 < z < 1.0).

We report the discovery of a fifth image in the large separation lensed quasar system SDSS J1004+4112. A faint point source located 0.2'' from the center of the brightest galaxy in the lensing cluster is detected in images taken with the Advanced Camera for Surveys (ACS) and the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) on the Hubble Space Telescope. The flux ratio between the point source and the brightest lensed component in the ACS image is similar to that in the NICMOS image. The location and brightness of the point source are consistent with lens model predictions for a lensed image. We therefore conclude that the point source is likely to be a fifth image of the source quasar. In addition, the NICMOS image reveals the lensed host galaxy of the source quasar, which can strongly constrain the structure of the lensing critical curves and thereby the mass distribution of the lensing cluster.

Flux ratio `anomalies' in quadruply-imaged gravitational lenses can be explained with galactic substructure of the sort predicted by CDM, but the strength and uniqueness of that hypothesis needs to be further assessed. A good way to do that is to use the physical scale associated with the size of the source quasar, and its dependence on wavelength. We develop an analytic toy model to study finite source effects in substructure lensing and identify general principles. We find that image positions and magnifications are basically independent of source size until the source is large enough to intersect a substructure caustic. Even sources that are much larger than the substructure Einstein radius can be perturbed at a detectable level. A tremendous amount could be learned by comparing image positions and magnifications at wavelengths that correspond to different source sizes. In a separate analysis, we carefully study four observed radio lenses to determine which of the images are anomalous. In B0712+472, if the flux ratio anomaly is real it is probably in image C. In B1422+231, the anomaly is in image A. Interestingly, B2045+265 and B1555+375 both appear to have two anomalous images. It remains to be seen whether CDM predicts enough substructure to explain multiple anomalies in multiple lenses. When we join our modeling results and substructure theory, we obtain lower bounds on the masses of the substructures responsible for the observed anomalies. The mass bounds are broadly consistent with expectations for CDM. Perhaps more importantly, we outline various systematic effects in the mass bounds; poor knowledge of whether the substructure lies within the main lens galaxy or elsewhere along the line of sight appears to be the dominant systematic. [Abridged]

We study the effects of ellipticity in lens galaxies and external tidal shear from neighboring objects on the statistics of strong gravitational lenses. For isothermal lens galaxies normalized so that the Einstein radius is independent of ellipticity and shear, ellipticity {\it reduces} the lensing cross section slightly, and shear leaves it unchanged. Ellipticity and shear can significantly enhance the magnification bias, but only if the luminosity function of background sources is steep. Realistic distributions of ellipticity and shear {\it lower} the total optical depth by a few percent for most source luminosity functions, and increase the optical depth only for steep luminosity functions. The boost in the optical depth is noticeable (>5%) only for surveys limited to the brightest quasars (L/L_* > 10). Ellipticity and shear broaden the distribution of lens image separations but do not affect the mean. Ellipticity and shear naturally increase the abundance of quadruple lenses relative to double lenses, especially for steep source luminosity functions, but the effect is not enough (by itself) to explain the observed quadruple-to-double ratio. With such small changes to the optical depth and image separation distribution, ellipticity and shear have a small effect on cosmological constraints from lens statistics: neglecting the two leads to biases of just Delta Omega_M = 0.00 \pm 0.01 and Delta Omega_Lambda = -0.02 \pm 0.01 (where the errorbars represent statistical uncertainties in our calculations).

We combine the Cosmic Lens All-Sky Survey (CLASS) with new Sloan Digital Sky Survey (SDSS) data on the local velocity dispersion distribution function of E/S0 galaxies, $\phi(\sigma)$, to derive lens statistics constraints on $\Omega_\Lambda$ and $\Omega_m$. Previous studies of this kind relied on a combination of the E/S0 galaxy luminosity function and the Faber-Jackson relation to characterize the lens galaxy population. However, ignoring dispersion in the Faber-Jackson relation leads to a biased estimate of $\phi(\sigma)$ and therefore biased and overconfident constraints on the cosmological parameters. The measured velocity dispersion function from a large sample of E/S0 galaxies provides a more reliable method for probing cosmology with strong lens statistics. Our new constraints are in good agreement with recent results from the redshift-magnitude relation of Type Ia supernovae. Adopting the traditional assumption that the E/S0 velocity function is constant in comoving units, we find a maximum likelihood estimate of $\Omega_\Lambda = 0.74$--0.78 for a spatially flat unvierse (where the range reflects uncertainty in the number of E/S0 lenses in the CLASS sample), and a 95% confidence upper bound of $\Omega_\Lambda<0.86$. If $\phi(\sigma)$ instead evolves in accord with extended Press-Schechter theory, then the maximum likelihood estimate for $\Omega_\Lambda$ becomes 0.72--0.78, with the 95% confidence upper bound $\Omega_\Lambda<0.89$. Even without assuming flatness, lensing provides independent confirmation of the evidence from Type Ia supernovae for a nonzero dark energy component in the universe.

Gravitational lensing can amplify the apparent brightness of distant sources. Images that are highly magnified are often part of multiply-imaged systems, but we consider the possibility of having large magnifications without additional detectable images. In rare but non-negligible situations, lensing can produce a singly highly magnified image; this phenomenon is mainly associated with massive cluster-scale halos (>~1e13.5 Msun). Alternatively, lensing can produce multiply-imaged systems in which the extra images are either unresolved or too faint to be detectable. This phenomenon is dominated by galaxies and lower-mass halos (<~1e12 Msun), and is very sensitive to the inner density profile of the halos. Although we study the general problem, we customize our calculations to four quasars at redshift z~6 in the Sloan Digital Sky Survey (SDSS), for which Richards et al. (2004) have ruled out the presence of extra images down to an image splitting of 0.3" and a flux ratio of f=0.01. We predict that 9-29% of all z~6 quasars that are magnified by a factor of mu>10 would lack detectable extra images, with 5-10% being true singly-imaged systems. The maximum of 29% is reached only in the unlikely event that all low-mass (<~1e10 Msun) halos have highly concentrated (isothermal) profiles. In more realistic models where dwarf halos have flatter (NFW) inner profiles, the maximum probability is ~10%. We conclude that the probability that all four SDSS quasars are magnified by a factor of 10 is <~1e-4. The only escape from this conclusion is if there are many (>10) multiply-imaged z~6 quasars in the SDSS database that have not yet been identified, which seems unlikely. In other words, lensing cannot explain the brightnesses of the z~6 quasars, and models that invoke lensing to avoid having billion-Msun black holes in the young universe are not viable.

While many strong gravitational lens galaxies are suspected to lie in groups or clusters of galaxies, environmental effects in lens models are often unconstrained and sometimes ignored. We show that this creates significant biases in a variety of lensing applications, by creating mock lenses associated with each of 13 galaxies in a realistic model group, and then analyzing them with standard techniques. We find that standard models of double lenses, which neglect environment, grossly overestimate the ellipticity of the lens galaxy (de/e~0.5) and the Hubble constant (dh/h~0.22). Standard models of quad lenses, which approximate the environment as a tidal shear, recover the ellipticity reasonably well (|de/e|<~0.24) but overestimate the Hubble constant (dh/h~0.15), and have significant (~30%) errors in the millilensing analyses used to constrain the amount of substructure in dark matter halos. For both doubles and quads, standard models slightly overestimate the velocity dispersion of the lens galaxy (d(sigma)/sigma~0.06), and underestimate the magnifications of the images (d(mu)/mu ~ -0.25). Standard analyses of lens statistics overestimate Omega_Lambda (by 0.05-0.14), and underestimate the ratio of quads to doubles (by a factor of 2). These biases help explain some long-standing puzzles (such as the high observed quad/double ratio), but aggravate others (such as the low value of H_0 inferred from lensing). Most of the biases are caused by neglect of the convergence from the mass associated with the environment, but additional uncertainty is introduced by neglect of higher-order terms. Fortunately, we show that directly observing and modeling lens environments should make it possible to remove the biases and reduce the uncertainties associated with environments to the few percent level. (Abridged)

We study the statistics of large-separation gravitational lens systems produced by non-spherical halos in the Cold Dark Matter (CDM) model. Specifically, we examine how the triaxiality of CDM halos affects the overall lensing probabilities and the relative numbers of different image configurations (double, quadruple, and naked cusp lenses). We find that triaxiality significantly enhances lensing probabilities by a factor of ~2-4, so it cannot be ignored. If CDM halos have central density slopes alpha<~1.5, we predict that a significant fraction (>~20%) of large-separation lenses should have naked cusp image configurations; this contrasts with lensing by isothermal (\alpha~2) galaxies where naked cusp configurations are rare. The image multiplicities depend strongly on the inner density slope \alpha: for \alpha=1, the naked cusp fraction is >~60%; while for \alpha=1.5, quadruple lenses are actually the most probable. Thus, the image multiplicities in large-separation lenses offer a simple new probe of the inner density profiles of dark matter halos. We also compute the expected probabilities and image multiplicities for lensed quasars in the Sloan Digital Sky Survey, and argue that the recent discovery of the large-separation quadruple lens SDSS J1004+4112 is consistent with expectations for CDM.

We present seven epochs of spectroscopy on the quadruply imaged quasar SDSS J1004+4112, spanning observed-frame time delays from 1 to 322 days. The spectra reveal differences in the emission lines between the lensed images. Specifically, component A showed a strong enhancement in the blue wings of several high-ionization lines relative to component B, which lasted at least 28 days (observed frame) then faded. Since the predicted time delay between A and B is <30 days, our time coverage suggests that the event was not intrinsic to the quasar. We attribute these variations to microlensing of part of the broad emission line region of the quasar, apparently resolving structure in the source plane on a scale of ~10^{16} cm at z=1.734. In addition, we observed smaller differences in the emission line profiles between components A and B that persisted throughout the time span, which may also be due to microlensing or millilensing. Further spectroscopic monitoring of this system holds considerable promise for resolving the structure of the broad emission line region in quasars.

We study the recently discovered gravitational lens SDSS J1004+4112, the first quasar lensed by a cluster of galaxies. It consists of four images with a maximum separation of 14.62''. The system has been confirmed as a lensed quasar at z=1.734 on the basis of deep imaging and spectroscopic follow-up observations. We present color-magnitude relations for galaxies near the lens plus spectroscopy of three central cluster members, which unambiguously confirm that a cluster at z=0.68 is responsible for the large image separation. We find a wide range of lens models consistent with the data, but they suggest four general conclusions: (1) the brightest cluster galaxy and the center of the cluster potential well appear to be offset by several kpc; (2) the cluster mass distribution must be elongated in the North--South direction, which is consistent with the observed distribution of cluster galaxies; (3) the inference of a large tidal shear (~0.2) suggests significant substructure in the cluster; and (4) enormous uncertainty in the predicted time delays between the images means that measuring the delays would greatly improve constraints on the models. We also compute the probability of such large separation lensing in the SDSS quasar sample, on the basis of the CDM model. The lack of large separation lenses in previous surveys and the discovery of one in SDSS together imply a mass fluctuation normalization \sigma_8=1.0^{+0.4}_{-0.2} (95% CL), if cluster dark matter halos have an inner slope -1.5. Shallower profiles would require higher values of \sigma_8. Although the statistical conclusion might be somewhat dependent on the degree of the complexity of the lens potential, the discovery is consistent with the predictions of the abundance of cluster-scale halos in the CDM scenario. (Abridged)

Gravitational lensing is a powerful tool for the study of the distribution of dark matter in the Universe. The cold-dark-matter model of the formation of large-scale structures predicts the existence of quasars gravitationally lensed by concentrations of dark matter so massive that the quasar images would be split by over 7 arcsec. Numerous searches for large-separation lensed quasars have, however, been unsuccessful. All of the roughly 70 lensed quasars known, including the first lensed quasar discovered, have smaller separations that can be explained in terms of galaxy-scale concentrations of baryonic matter. Although gravitationally lensed galaxies with large separations are known, quasars are more useful cosmological probes because of the simplicity of the resulting lens systems. Here we report the discovery of a lensed quasar, SDSS J1004+4112, which has a maximum separation between the components of 14.62 arcsec. Such a large separation means that the lensing object must be dominated by dark matter. Our results are fully consistent with theoretical expectations based on the cold-dark-matter model.

Using stellar dynamics and strong gravitational lensing as complementary probes, Sand et al. (2002, 2003) have recently claimed strong evidence for shallow dark matter density profiles in several lensing clusters, which may conflict with predictions of the Cold Dark Matter paradigm. However, systematic uncertainties in the analysis weaken the constraints. By re-analyzing their data, we argue that the tight constraints claimed by Sand et al., were driven by prior assumptions. Relaxing the assumptions, we find that no strong constraints may be derived on the dark matter inner profile from the Sand et al. data; we find satisfactory fits (with reasonable parameters) for a wide range of inner slopes 0 < beta < 1.4. Useful constraints on the mass distributions of lensing clusters can still be obtained, but they require moving beyond mere measurements of lensing critical radii into the realm of detailed lens modeling.

We study lens statistics in flat, low-density universes with different equations of state $w=p_Q/\rho_Q$ for the dark energy component. Dark energy modifies the distance-redshift relation and the mass function of dark matter halos leading to changes in the lensing optical depth as a function of image separation. Those effects must, however, be distinguished from effects associated with the structure of dark matter halos. Baryonic cooling causes galaxy-mass halos to have different central density profiles than group- and cluster-mass halos, which causes the distribution of normal arcsecond-scale lenses to differ from the distribution of ``wide-separation'' ($\Delta\theta \gtrsim 4\arcsec$) lenses. Fortunately, the various parameters related to cosmology and halo structure have very different effects on the overall image separation distribution: (1) the abundance of wide-separation lenses is exremely sensitive (by orders of magnitude) to the distribution of ``concentration'' parameters for massive halos modeled with the Navarro-Frenk-White profile; (2) the transition between normal and wide-separation lenses depends mainly on the mass scale where baryonic cooling ceases to be efficient; and (3) dark energy has effects at all image separation scales. While current lens samples cannot usefully constrain all of the parameters, ongoing and future imaging surveys should discover hundreds or thousands of lenses and make it possible to disentangle the various effects and constrain all of the parameters simultaneously. (abridged)

The inability of standard models to explain the flux ratios in many 4-image gravitational lenses has been cited as evidence for significant small-scale structure in lens galaxies. That claim has generally relied on detailed lens modeling, so it is both model dependent and somewhat difficult to interpret. We present a more robust and generic method for identifying lenses with small-scale structure. For a close triplet of images associated with a source near a cusp caustic, the sum of the signed magnifications should approximately vanish. We derive realistic upper bounds on the sum, and argue that lenses with flux ratios that significiantly violate the bounds can be said to have structure in the lens potential on scales smaller than the image separation. Five observed lenses have such flux ratio ``anomalies'': B2045+265, 1RXS J1131-1231, and SDSS J0924+0219 have strong anomalies; B0712+472 has a strong anomaly at optical/near-IR wavelengths and a marginal anomaly at radio wavelengths; and RX J0911+0551 appears to have an anomaly, but this conclusion is subject to uncertainties about octopole modes in early-type galaxies. Analysis of the cusp relation does not identify the known anomaly in B1422+231, so methods that are more sophisticated (and less generic) than the cusp relation may be necessary to uncover flux ratio anomalies in some systems. Although these flux ratio anomalies might represent milli- or micro-lensing, we cannot identify the cause; we can only conclude that the lenses have significant structure in the potential on scales smaller than the separation between the images. Additional arguments must be invoked to specify the nature of this small-scale structure. [Abridged]

We report the discovery of a new gravitationally lensed quasar from the Sloan Digital Sky Survey, SDSS J090334.92+502819.2. This object was targeted for SDSS spectroscopy as a Luminous Red Galaxy (LRG), but manual examination of the spectrum showed the presence of a quasar at z= 3.6 in addition to a red galaxy at z=0.388, and the SDSS image showed a second possible quasar image nearby. Follow-up imaging and spectroscopy confirmed the lensing hypothesis. In images taken at the ARC 3.5-meter telescope, two quasars are separated by 2.8 arc-seconds; the lensing galaxy is clearly seen and is blended with one of the quasar images. Spectroscopy taken at the Keck II telescope shows that the quasars have identical redshifts of z=3.6 and both show the presence of the same broad absorption line-like troughs. We present simple lens models which account for the geometry and magnifications. The lens galaxy lies near two groups of galaxies and may be a part of them. The models suggest that the groups may contribute considerable shear and may have a strong effect on the lens configuration.

We introduce a self-similar mass model for early-type galaxies, and constrain it using the aperture mass-radius relations determined from the geometries of 22 gravitational lenses. The model consists of two components: a concentrated component which traces the light distribution, and a more extended power-law component (rho propto r^-n) which represents the dark matter. We find that lens galaxies have total mass profiles which are nearly isothermal, or slightly steeper, on the several-kiloparsec radial scale spanned by the lensed images. In the limit of a single-component, power-law radial profile, the model implies n=2.07+/-0.13, consistent with isothermal (n=2). Models in which mass traces light are excluded at >99 percent confidence. An n=1 cusp (such as the Navarro-Frenk-White profile) requires a projected dark matter mass fraction of f_cdm = 0.22+/-0.10 inside 2 effective radii. These are the best statistical constraints yet obtained on the mass profiles of lenses, and provide clear evidence for a small but non-zero dark matter mass fraction in the inner regions of early-type galaxies. In addition, we derive the first strong lensing constraint on the relation between stellar mass-to-light ratio (Upsilon) and galaxy luminosity (L): Upsilon propto L^[0.14 (+0.16)(-0.12)], which is consistent with the relation suggested by the fundamental plane. Finally, we apply our self-similar mass models to current problems regarding the interpretation of time delays and flux ratio anomalies in gravitational lens systems.

Multiply-imaged quasar lenses can be used to constrain the substructure mass fraction in galaxy-sized dark matter halos via anomalous flux ratios in lensed images. The flux ratios, however, can be affected by both the substructure in the lens halo and by isolated small-mass halos along the entire line-of-sight to the lensed source. While lensing by dark matter clumps near the lens galaxy is more efficient than elsewhere, the cumulative effect of all objects along the line-of-sight could be significant. Here we estimate the potential contribution of isolated clumps to the substructure lensing signal using a simple model motivated by cosmological simulations. We find that the contribution of isolated clumps to the total lensing optical depth ranges from a few to tens percent, depending on assumptions and the particular configuration of the lens. Therefore, although the contribution of isolated clumps to the lensing signal is not dominant, it should not be neglected in detailed analyses of substructure lensing. For the currently favored Lambda-CDM model, the total calculated optical depth for lensing is high, tau ~ 0.2-20 and could, therefore, naturally explain the high frequency of anomalous flux ratios in observed lenses. The prediction, however, is highly sensitive to the spatial distribution of substructure halos in the innermost regions of the lens halo, which is still very uncertain. Therefore, constraints on the properties of the substructure population or accurate cosmological constraints, such as the mass of the warm dark matter particle, are difficult if not impossible to derive at this point.

We present results from high-resolution radio and optical observations of PMN J0134-0931, a gravitational lens with a unique radio morphology and an extremely red optical counterpart. Our data support the theory of Keeton & Winn (2003): five of the six observed radio components are multiple images of a single quasar, produced by a pair of lens galaxies. Multi-frequency VLBA maps show that the sixth and faintest component has a different radio spectrum than the others, confirming that it represents a second component of the background source rather than a sixth image. The lens models predict that there should be additional faint images of this second source component, and we find evidence for one of the predicted images. The previously-observed large angular sizes of two of the five bright components are not intrinsic (which would have excluded the possibility that they are lensed images), but are instead due to scatter broadening. Both the extended radio emission observed at low frequencies, and the intrinsic image shapes observed at high frequencies, can be explained by the lens models. The pair of lens galaxies is marginally detected in HST images. The differential extinction of the quasar images suggests that the extreme red color of the quasar is at least partly due to dust in the lens galaxies.

The strange morphology of the six-component gravitational lens PMN J0134-0931 has resisted explanation. We present the first successful quantitative models for the system, based on the idea that there are two lens galaxies and two components of the background source. One source is quintuply imaged and corresponds to the five brightest observed radio components. The other source is triply imaged and corresponds to the sixth component, along with two others too faint to have been detected. The models reproduce the observed image positions and fluxes, and make falsifiable predictions about other properties of the system. Some of these predictions have been confirmed by high-resolution radio and optical observations, as described in the companion paper by Winn et al. (2003). Although we cannot determine the lens model uniquely with current data, we predict that the lens galaxies are spiral galaxies with roughly equal velocity dispersions sigma~120km/s and a projected separation of only 0.4" (2/h kpc at z_l=0.76). This system is the first known lens with five images of a single quasar, and the second with more than four images.

We investigate the evolution of mass-selected early-type field galaxies using a sample of 28 gravitational lenses spanning the redshift range 0 < z < 1. Based on the redshift-dependent intercept of the fundamental plane in the rest frame B band, we measure an evolution rate of d log (M/L)_B / dz = -0.56 +/- 0.04 (all errors are 1 sigma unless noted) if we directly compare to the local intercept measured from the Coma cluster. Re-fitting the local intercept helps minimize potential systematic errors, and yields an evolution rate of d log (M/L)_B / dz = -0.54 +/- 0.09. An evolution analysis of properly-corrected aperture mass-to-light ratios (defined by the lensed image separations) is closely related to the Faber-Jackson relation. In rest frame B band we find an evolution rate of d log (M/L)_B / dz = -0.41 +/- 0.21, a present-day characteristic magnitude of M_{*0} = -19.70 + 5 log h +/- 0.29 (assuming a characteristic velocity dispersion of sigma_{DM*} = 225 km/s), and a Faber-Jackson slope of gamma_{FJ} = 3.29 +/- 0.58. The measured evolution rates favor old stellar populations (mean formation redshift z_f > 1.8 at 2 sigma confidence for a Salpeter initial mass function and a flat Omega_m =0.3 cosmology) among early-type field galaxies, and argue against significant episodes of star formation at z < 1.

The magnifications of the images in a strong gravitational lens system are sensitive to small mass clumps in the lens potential; this effect has been used to infer the amount of substructure in galaxy dark matter halos. I study the theory of substructure lensing to identify important general features, and to compute analytic cross sections that will facilitate further theoretical studies. I show that the problem of a clump anywhere along the line of sight to a lens can be mapped onto an equivalent problem of a clump in a simple convergence and shear field; clumps at arbitrary redshifts are therefore not hard to handle in calculations. For clumps modeled as singular isothermal spheres (SIS), I derive simple analytic estimates of the cross section for magnification perturbations of a given strength. The results yield two interesting conceptual points. First, lensed images with positive parity are always made brighter by SIS clumps; images with negative parity can be brightened but are much more likely to be dimmed. Second, the clumps need not lie within the lens galaxy; they can be moved in redshift by several tenths and still have a significant lensing effect. Isolated small halos are expected to be common in hierarchical structure formation models, but it is not yet known whether they are abundant enough compared with clumps inside lens galaxies to affect the interpretation of substructure lensing.

The Sloan Digital Sky Survey (SDSS) automatically targeted as a quasar candidate the recently discovered, gravitationally lensed, extremely reddened z=2.2 quasar PMN J0134-0931. The SDSS spectrum exhibits Ca II absorption at z=0.76451, which we identify as the redshift of a lensing galaxy. Hubble Space Telescope imaging shows that components CDE of the system are significantly redder than components A or B and detects faint galaxy emission between D and A+B. The redshift of the dust responsible for the reddening remains unconstrained with current data. However, we outline a model wherein lensing and differential reddening by a z=0.76451 galaxy pair can entirely explain this system. Detecting mm-wave molecular line absorption from the lensing galaxy or galaxies may be possible in PMN J0134-0931, just as in the lenses PKS1830-211 and B0218+357. Well-constructed optical quasar surveys like the SDSS can contribute to the detection and study of reddened quasars. (Expanded)

Strong gravitational lensing has traditionally been one of the few phenomena said to oppose a large cosmological constant; many analyses of lens statistics have given upper limits on $\Omega_\Lambda$ that are marginally inconsistent with the concordance cosmology. Those conclusions were based on models where the predicted number counts of galaxies at moderate redshifts (z~0.5-1) increased significantly with $\Omega_\Lambda$. I argue that the models should now be calibrated by counts of distant galaxies. When this is done lens statistics lose most of their sensitivity to the cosmological constant.

Gravitational lensing provides a unique probe of the inner 10-1000 pc of distant galaxies (z=0.2-1). Lens theory predicts that every strong lens system should have a faint image near the center of the lens galaxy, which should be visible in radio lenses but have not been observed. We study these ``core'' images using models derived from the stellar distributions in nearby early-type galaxies. We find that realistic galaxies predict a remarkably wide range of core images, with lensing magnifications spanning some six orders of magnitude. More concentrated galaxies produce fainter core images, although not with any simple, quantitative, model independent relation. Some real galaxies have diffuse cores and predict bright core images (magnification mu>~0.1), but more common are galaxies that predict faint core images (mu<~0.001). Thus, stellar mass distributions alone are probably concentrated enough to explain the lack of observed core images, and may require observational sensitivity to improve by an order of magnitude before detections of core images become common. Two-image lenses will tend to have brighter core images than four-image lenses, so they will be the better targets for finding core images and exploiting these tools for studying the central mass distributions of distant galaxies.

Many different families of mass models are used in modern applications of strong gravitational lensing. I review a wide range of popular models, with two points of emphasis: (1) a discussion of strategies for building models suited to a particular lensing problem; and (2) a summary of technical results for a canonical set of models. All of the models reviewed here are included in publicly-available lensing software called gravlens.

We study the prospects for three-dimensional mapping of the dark matter to high redshift through the shearing of faint galaxies images at multiple distances by gravitational lensing. Such maps could provide invaluable information on the nature of the dark energy and dark matter. While in principle well-posed, mapping by direct inversion introduces exceedingly large, but usefully correlated noise into the reconstruction. By carefully propagating the noise covariance, we show that lensing contains substantial information, both direct and statistical, on the large-scale radial evolution of the density field. This information can be efficiently distilled into low-order signal-to-noise eigenmodes which may be used to compress the data by over an order of magnitude. Such compression will be useful for the statistical analysis of future large data sets. The reconstructed map also contains useful information on the localization of individual massive dark matter halos, and hence the dark energy from halo number counts, but its extraction depends strongly on prior assumptions. We outline a procedure for maximum entropy and point-source regularization of the maps that can identify alternate reconstructions.

We report on a program to obtain HST observations of galaxy-mass gravitational lens systems at optical and infrared wavelengths. Here we discuss the properties of 10 two-image gravitational lens systems (Q0142-100=UM673, B0218+357, SBS0909+532, BRI0952-0115, LBQS1009-0252, Q1017-207=J03.13, B1030+074, HE1104-1805, Q1208+1011, and PKS1830-211). We grouped these 10 systems because they have limited lens model constraints and often show poor contrast between the images and the lens galaxy. Of the 10 lens galaxies, 7 are probably early-type galaxies, 2 are probably late-type galaxies (B0218+357 and PKS1830-211), and one was not detected (Q1208+1011). We detect the host galaxies of the z_s=4.50 lensed quasar in BRI0952-0115, the z_s=2.32 lensed quasar in HE1104-1805, and the unlensed z=1.63 quasar near LBQS1009-0252. We fit a set of four standard lens models to each lens that had sufficient constraints to compare isothermal dark matter and constant mass-to-light lens models, and to explore the effects of local tidal shears.

By directly probing mass distributions, gravitational lensing offers several new tests of the CDM paradigm. Lens statistics place upper limits on the dark matter content of elliptical galaxies. Galaxies built from CDM mass distributions are too concentrated to satisfy these limits, so lensing extends the ``concentration problem'' in CDM to elliptical galaxies. The central densities of the model galaxies are too low on ~10 pc scales to agree with the lack of central images in observed lenses. The flux ratios of four-image lenses imply a substantial population of dark matter clumps with a typical mass ~10^6 Msun. Thus, lensing implies the need for a mechanism that reduces dark matter densities on kiloparsec scales without erasing structure on smaller scales.

Substructure in distant gravitational lens galaxies can be detected because it alters the brightnesses and shapes of the lensed images. The optical and radio flux ratios in the four-image lens B1422+231 imply that there is a \sim 10^4-10^7 Msun/h mass clump in the lens galaxy in front of image A; this is the first constraint on the mass of a particular clump lying in a distant galaxy and detected by its mass. The data also indicate that a small clump, perhaps a star, is passing in front of image B and making the optical flux ratios variable. Both of these hypotheses can be tested with new observations. B1422+231 demonstrates how data at different wavelengths can be used in individual lenses to probe individual mass clumps in distant galaxies.

The statistics of gravitationally lensed arcs, which can be used for a variety of cosmological tests, are sensitive to the intrinsic shapes of the source galaxies. I present an analytic formalism that makes it simple to include elliptical sources in analytic calculations of lens statistics. For cuspy lens models, sources with an axis ratio of 2:1 enhance the total number of arcs longer than 10:1 by a factor of order two, while modestly decreasing the number ratio of radial arcs to tangential arcs. Source ellipticity is therefore an important systematic effect in detailed quantitative studies, but it should not hinder cosmological applications such as attempts to constrain cluster dark matter profiles with arc statistics.

(Abridged) Using the number and sizes of observed gravitational lenses, I derive upper limits on the dark matter content of elliptical galaxies. Galaxies built from Cold Dark Matter (CDM) mass distributions are too concentrated to comfortably satisfy these limits; SCDM is ruled out, and LCDM is only marginally consistent with the data. Thus, lensing adds to the evidence that CDM mass distributions are too concentrated on kiloparsec scales to agree with real galaxies, and extends the argument to elliptical galaxies. By contrast, the lack of central images in radio lenses implies that the central densities of CDM galaxies are too low on ~10 parsec scales, even if supermassive black holes are included. Self-interacting dark matter, or any other modification to regular cold dark matter, must simultaneously reduce the densities on kiloparsec scales and increase the densities on parsec scales in order to satisfy the unique constraints from lensing.

Recent observations of galaxy luminosity profiles and dark matter simulations find luminosity and mass distributions characterized by central cusps rather than finite core radii. We introduce and implement a set of cusped ellipsoidal lens models which include limits similar to the Jaffe, Hernquist, \eta and NFW models and apply them to the gravitational lenses APM 08279+5255 and B 1933+503. A successful model of APM 08279+5255 with its central, odd image requires a very shallow cusp, $\gamma \ltorder 0.4$ where $\rho \propto r^{-\gamma}$ as $r\to 0$, which is similar to a core rather than the favored $1 \ltorder \gamma \ltorder 2$ cusps. B~1933+503, by contrast, is well modeled with a steep density cusp, $1.6 \ltorder \gamma \ltorder 2.0$.

We consider a wide range of parametric mass models for B 1933+503, a ten-image radio lens, and identify shared properties of the models with the best fits. The approximate rotation curves varies by less than 8.5% from the average value between the innermost and the outermost image (1.5h^{-1} kpc to 4.1h^{-1} kpc) for models within 1 \sigma of the best fit, and the radial dependence of the shear strength and angle also have common behavior for the best models. The time delay between images 1 and 6, the longest delay between the radio cores, is \Delta t = (10.6^{+2.4}_{-1.1})h^{-1} days (\Omega_0=0.3, \lambda_0=0.7) including all the modeling uncertainties. Deeper infrared observations, to more precisely register the lens galaxy with the radio images and to measure the properties of the Einstein ring image of the radio source's host galaxy, would significantly improve the model constraints and further reduce the uncertainties in the mass distribution and time delay.

Modern applications of strong gravitational lensing require the ability to use precise and varied observational data to constrain complex lens models. I discuss two sets of computational methods for lensing calculations. The first is a new algorithm for solving the lens equation for general mass distributions. This algorithm makes it possible to apply arbitrarily complicated models to observed lenses. The second is an evaluation of techniques for using observational data including positions, fluxes, and time delays of point-like images, as well as maps of extended images, to constrain models of strong lenses. The techniques presented here are implemented in a flexible and user-friendly software package called gravlens, which is made available to the community.

The statistics of wide-separation (6'' < theta < 15'') gravitational lenses constrain the amount of mass in the cores of dark matter halos on group and cluster mass scales. For a family of halo models with a central cusp $\rho \propto r^{-alpha}$ (1.0 <= alpha <= 1.9), the lack of wide-separation lenses in the large Cosmic Lens All-Sky Survey yields an upper limit on the fraction of the halo mass that is contained within ~4% of the virial radius, f_{core} < 0.023 (95% confidence level, LCDM). This limit offers an important test of the cold dark matter paradigm. While the halo profiles derived from numerical simulations appear to be consistent with this upper limit, larger surveys currently underway such as the 2dF and SDSS should detect wide-separation lenses and thus provide a measurement of the core mass fraction in massive dark matter halos.

HST V and I-band observations show that the gravitational lens B1359+154 consists of six images of a single z_s=3.235 radio source and its star-forming host galaxy, produced by a compact group of galaxies at z_l = 1. VLBA observations at 1.7 GHz strongly support this conclusion, showing six compact cores with similar low-frequency radio spectra. B1359+154 is the first example of galaxy-scale gravitational lensing in which more than four images are observed of the same background source. The configuration is due to the unique lensing mass distribution: three primary lens galaxies lying on the vertices of a triangle separated by 0.7 arcsec (4/h kpc), inside the 1.7 arcsec diameter Einstein ring defined by the radio images. The gravitational potential has additional extrema within this triangle, creating a pair of central images that supplement the ``standard'' four-image geometry of the outer components. Simple mass models consisting of three lens galaxies constrained by HST and VLBA astrometry naturally reproduce the observed image positions but must be finely-tuned to fit the flux densities.

We report the discovery of a new two-image gravitational lens candidate. The system RX J0921+4529 contains two z_s=1.66 quasars separated by 6."93 with an H band magnitude difference of $\Delta m=1.39$. The HST NIC2 H band images reveal an H=18.2 spiral galaxy between the quasar images, which is probably a member of a z_l=0.32 X-ray cluster centered on the field. We detect an extended source near the fainter quasar image but not in the brighter image. If this extended source is the host galaxy of the fainter quasar, then the system is a binary quasar rather than a gravitational lens. VLA observations at 3.6 cm reveal emission from the lens galaxy at the flux level of 1 mJy and a marginal detection of the brighter quasar.

We predict how the observed variations in galaxy populations with environment affect the number and properties of gravitational lenses in different environments. Two trends dominate: lensing strongly favors early-type galaxies, which tend to lie in dense environments, but dense environments tend to have a larger ratio of dwarf to giant galaxies than the field. The two effects nearly cancel, and the distribution of environments for lens and non-lens galaxies are not substantially different (lens galaxies are slightly less likely than non-lens galaxies to lie in groups and clusters). We predict that about 20% of lens galaxies are in bound groups (defined as systems with a line-of-sight velocity dispersion sigma in the range 200 < sigma < 500 km/s), and another roughly 3% are in rich clusters (sigma > 500 km/s). Therefore at least roughly 25% of lenses are likely to have environments that significantly perturb the lensing potential. If such perturbations do not significantly increase the image separation, we predict that lenses in groups have a mean image separation that is about 0.2'' smaller than that for lenses in the field and estimate that 20-40 lenses in groups are required to test this prediction with significance. The tail of the distribution of image separations is already illuminating. Although lensing by galactic potential wells should rarely produce lenses with image separations theta >~ 6'', two such lenses are seen among 49 known lenses, suggesting that environmental perturbations of the lensing potential can be significant. Further comparison of theory and data will offer a direct probe of the dark halos of galaxies and groups and reveal the extent to which they affect lensing estimates of cosmological parameters.

We develop a theory of Einstein rings and demonstrate it using the infrared Einstein ring images of the quasar host galaxies observed in PG1115+080, B1608+656 and B1938+666. The shape of an Einstein ring accurately and independently determines the shape of the lens potential and the shape of the lensed host galaxy. We find that the host galaxies of PG1115+080, B1608+656 and B1938+666 have axis ratios of 0.58+/-0.02, 0.69+/-0.02 and 0.62+/-0.15 including the uncertainties in the lens models. The Einstein rings break the degeneracies in the mass distributions or Hubble constants inferred from observations of gravitational lenses. In particular, the Einstein ring in PG1115+080 rules out the centrally concentrated mass distributions that lead to a high Hubble constant (H_0>60 km/s Mpc) given the measured time delays. Deep, detailed observations of Einstein rings will be revolutionary for constraining mass models and determining the Hubble constant from time delay measurements.

Most gravitational lenses are early-type galaxies in relatively low density environments -- a ``field'' rather than a ``cluster'' population. We show that field early-type galaxies with 0 < z < 1, as represented by the lens galaxies, lie on the same fundamental plane as those in rich clusters at similar redshifts. We then use the fundamental plane to measure the combined evolutionary and K-corrections for early-type galaxies in the V, I and H bands. Only for passively evolving stellar populations formed at z > 2 (H_0=65 km/s Mpc, Omega_0=0.3, Lambda_0=0.7) can the lens galaxies be matched to the local fundamental plane. The high formation epoch and the lack of significant differences between the field and cluster populations contradict many current models of the formation history of early-type galaxies. Lens galaxy colors and the fundamental plane provide good photometric redshift estimates with an empirical accuracy of -0.03 +/- 0.11 for the 17 lenses with known redshifts. A mass model dominated by dark matter is more consistent with the data than either an isotropic or radially anisotropic constant M/L mass model, and a radially anisotropic model is better than an isotropic model.

We study gravitational lens time delays for a general family of lensing potentials, which includes the popular singular isothermal elliptical potential and singular isothermal elliptical density distribution but allows general angular structure. Using a novel approach, we show that the time delay can be cast in a very simple form, depending only on the observed image positions. Including an external shear changes the time delay proportional to the shear strength, and varying the radial profile of the potential changes the time delay approximately linearly. These analytic results can be used to obtain simple estimates of the time delay and the Hubble constant in observed gravitational lenses. The naive estimates for four of five time delay lenses show surprising agreement with each other and with local measurements of H_0; the complicated Q 0957+561 system is the only outlier. The agreement suggests that it is reasonable to use simple isothermal lens models to infer H_0, although it is still important to check this conclusion by examining detailed models and by measuring more lensing time delays.

We study the problem of gravitational lensing by an isothermal elliptical density galaxy in the presence of a tidal perturbation. When the perturbation is fairly strong and oriented near the galaxy's minor axis, the lens can produce image configurations with six or even eight highly magnified images lying approximately on a circle. We classify the caustic structures in the model and identify the range of models that can produce such lenses. Sextuple and octuple lenses are likely to be rare because they require special lens configurations, but a full calculation of the likelihood will have to include both the existence of lenses with multiple lens galaxies and the strong magnification bias that affects sextuple and octuple lenses. At optical wavelengths these lenses would probably appear as partial or complete Einstein rings, but at radio wavelengths the individual images could probably be resolved.

Infrared images of the Q 0957+561 gravitational lens obtained with the Hubble Space Telescope show two large (\sim 5 arcsec) lensed images of the z_s=1.41 quasar host galaxy. Parts of the host galaxy are doubly-imaged like the quasar, while other parts are quadruply-imaged. The distortions of the host galaxy offer the best probe yet of the global structure of the lensing potential, which is essential for determining the Hubble constant from the measured time delay. The distortions are inconsistent with the predictions of previously published lens models, which invalidates those models and their implications for H_0. New models show that the distortions finally break the long-standing degeneracy between the shape of the lens galaxy and the tidal shear contributed by the cluster containing the lens galaxy. The shape of the lens galaxy's mass distribution must be remarkably similar to the shape of its luminosity distribution, and most models that produce reasonable values for the Hubble constant roughly match the observed ellipticity gradient and isophote twist of the lens galaxy. Also, the cluster must be non-spherical and produce a relatively small tidal shear. Although there are still degeneracies in the lens models that lead to a 25% uncertainty in the derived value of the Hubble constant, there are also strong prospects for new observations to further improve the constraints and reduce the uncertainties.

Most gravitational lens galaxies are early-type galaxies in relatively low density environments. We show that they lie on the same fundamental plane as early-type galaxies in both local and distant rich clusters. Their surface brightness evolution requires a typical star formation epoch of z=2-3, almost indistinguishable from that of rich cluster galaxies at comparable redshifts. The restricted galaxy type range of the lenses means that photometric redshifts work well even with only 1-3 filter photometry. We make preliminary measurements of the mass and luminosity functions of the lens galaxies, and find they are consistent with the standard model used for deriving cosmological limits using lens statistics. As expected for a mass-weighted sample, they are more massive and more luminous than the overall early-type galaxy population.

We use observations from the CASTLES survey of gravitational lenses to study extinction in 23 lens galaxies with $0 < z_l < 1$. The median differential extinction between lensed images is $\Delta E(\bv) = 0.05$ mag, and the directly measured extinctions agree with the amount needed to explain the differences between the statistics of radio and (optical) quasar lens surveys. We also measure the first extinction laws outside the local universe, including an $R_V=7.2$ curve for a molecular cloud at $z_l=0.68$ and an $R_V=1.5$ curve for the dust in a redshift $z_l=0.96$ elliptical galaxy.

We describe a series of new applications of gravitational lenses as astrophysical and cosmological tools. Such applications are becoming possible thanks to advances in the quality and quantity of observations. CASTLES (CfA-Arizona-Space-Telescope-LEns-Survey) {cfa-www.harvard.edu/castles} is an ongoing project that exploits the sensitivity and resolution of the Hubble Space Telescope (HST) at optical and infrared wavelengths to study the sample of over 50 known gravitational lenses. The observational goal of CASTLES is a uniform sample of multi-band images of all known galaxy-mass lens systems, to derive precise photometry and astrometry for the lens galaxies, all the known images, and any source or lens components that might have escaped detection. With these measurements we are investigating: (1) the properties of dust and of dark matter in lens galaxies out to z\sim1; (2) the dark matter in lens galaxies and in their environments; (3) the evolution of lens galaxies; and (4) the cosmological model, for instance by refining constraints on the Hubble constant H_0.

The lens candidate RXJ 0921+4529 consists of two z_s=1.66 quasar separated by 6."93 with an H band magnitude difference of \Delta m=1.39. The lens appears to be a z_l=0.31 X-ray cluster, including a m_H=18.5 late-type galaxy lying between the quasar images. We detect an extended source overlapping the faint quasar but not the bright quasar. If this extended source is the host galaxy of the fainter quasar, then the system is a quasar binary rather than a gravitational lens.

We examine models of the mass distribution for the first known case of gravitational lensing. Several new sets of constraints are used, based on recent observations. We remodel the VLBI observations of the radio jets in the two images of Q0957+561, showing that the previously derived positions and uncertainties were incorrect. We use as additional constraints the candidate lensed pairs of galaxies identified recently with the Hubble Space Telescope. We explore a wider range of lens models than before, and find that the Hubble constant is not tightly constrained once elliptical lens models are considered. We also discuss the systematic uncertainties caused by the cluster containing the lens galaxy. We conclude that additional observations of the candidate lensed galaxies as well as direct measurements of the cluster mass profile are needed before an accurate value of the Hubble constant can be derived from this lens.

Deep HST/NICMOS H (F160W) band observations of the z=1.96 quasar pair Q 1634+267A,B reveal no signs of a lens galaxy to a 1 sigma threshold of approximately 22.5 mag. The minimum luminosity for a normal lens galaxy would be a 6L_* galaxy at z > 0.5, which is 650 times greater than our detection threshold. Our observation constrains the infrared mass-to-light ratio of any putative, early-type, lens galaxy to (M/L)_H > 690h_65 (1200h_65) for Omega_0=0.1 (1.0) and H_0=65h_65 km/s/Mpc. We would expect to detect a galaxy somewhere in the field because of the very strong Mg II absorption lines at z=1.1262 in the Q 1634+267 A spectrum, but the HST H-band, I-band (F785LP) and V-band (F555W) images require that any associated galaxy be very under-luminous less than 0.1 L^*_H (1.0 L^*_I) if it lies within less than 40 h^{-1} (100 h^{-1}) kpc from Q 1634+267 A,B. While the large image separation (3.86 arcsec) and the lack of a lens galaxy strongly favor interpreting Q 1634+267A,B as a binary quasar system, the spectral similarity remains a puzzle. We estimate that at most 0.06% of randomly selected quasar pairs would have spectra as similar to each other as the spectra of Q 1634+267 A and B. Moreover, spectral similarities observed for the 14 quasar pairs are significantly greater than would be expected for an equivalent sample of randomly selected field quasars. Depending on how strictly we define similarity, we estimate that only 0.01--3% of randomly drawn samples of 14 quasar pairs would have as many similar pairs as the observational sample.

We determine 37 differential extinctions in 23 gravitational lens galaxies over the range 0 < z_l < 1. Only 7 of the 23 systems have spectral differences consistent with no differential extinction. The median differential extinction for the optically-selected (radio-selected) subsample is E(B-V)=0.04 (0.06) mag. The extinction is patchy and shows no correlation with impact parameter. The median total extinction of the bluest images is E(B-V)=0.08 mag, although the total extinction distribution is dominated by the uncertainties in the intrinsic colors of quasars. The directly measured extinction distributions are consistent with the mean extinction estimated by comparing the statistics of quasar and radio lens surveys, thereby confirming the need for extinction corrections when using the statistics of lensed quasars to estimate the cosmological model. A disjoint subsample of two face-on, radio-selected spiral lenses shows both high differential and total extinctions, but standard dust-to-gas ratios combined with the observed molecular gas column densities overpredict the amount of extinction by factors of 2-5. For several systems we can estimate the extinction law, ranging from R_V=1.5+/-0.2 for a z_l=0.96 elliptical, to R_V=7.2+/-0.1 for a z_l=0.68 spiral. For the four radio lenses where we can construct non-parametric extinction curves we find no evidence for gray dust over the IR-UV wavelength range. The dust can be used to estimate lens redshifts with reasonable accuracy, although we sometimes find two degenerate redshift solutions.

We have obtained and modeled new NICMOS images of the lens system MG1131+0456, which show that its lens galaxy is an H=18.6 mag, transparent, early-type galaxy at a redshift of about z_l = 0.85; it has a major axis effective radius R_e=0.68+/-0.05 arcsec, projected axis ratio b/a=0.77+/-0.02, and major axis PA=60+/-2 degrees. The lens is the brightest member of a group of seven galaxies with similar R-I and I-H colors, and the two closest group members produce sufficient tidal perturbations to explain the ring morphology. The host galaxy of the MG1131+0456 source is a z_s > 2 ERO (``extremely red object'') which is lensed into optical and infrared rings of dramatically different morphologies. These differences imply a strongly wavelength-dependent source morphology that could be explained by embedding the host in a larger, dusty disk. At 1.6 micron (H), the ring is spectacularly luminous, with a total observed flux of H=17.4 mag and a de-magnified flux of 19.3 mag, corresponding to a 1-2L_* galaxy at the probable source redshift of z_s > 2. Thus, it is primarily the stellar emission of the radio source host galaxy that produces the overall colors of two of the reddest radio lenses, MG1131+0456 and B~1938+666, aided by the suppression of optical AGN emission by dust in the source galaxy. The dusty lens hypothesis -- that many massive early-type galaxies with 0.2 < z_l < 1.0 have large, uniform dust opacities -- is ruled out.

Hubble Space Telescope observations of the gravitational lens PG 1115+080 in the infrared show the known z =0.310 lens galaxy and reveal the z = 1.722 quasar host galaxy. The main lens galaxy G is a nearly circular (ellipticity < 0.07) elliptical galaxy with a de Vaucouleurs profile and an effective radius of R_e = 0.59 +/- 0.06 arcsec (1.7 +/- 0.2 h^{-1} kpc for Omega = 1 and h = H_0/100 km/s/Mpc). G is part of a group of galaxies that is a required component of all successful lens models. The new quasar and lens positions (3 milliarcsecond errors) yield constraints for these models that are statistically degenerate, but several conclusions are firmly established. (1) The principal lens galaxy is an elliptical galaxy with normal structural properties, lying close to the fundamental plane for its redshift. (2) The potential of the main lens galaxy is nearly round, even when not constrained by the small ellipticity of the light of this galaxy. (3) All models involving two mass distributions place the group component near the luminosity-weighted centroid of the brightest nearby group members. (4) All models predict a time delay ratio r_{ABC} = 1.3. (5) Our lens models predict H_0 = 44 +/- 4 km/s/Mpc if the lens galaxy contains dark matter and has a flat rotation curve, and H_0 = 65 +/- 5 km/s/Mpc if it has a constant mass-to-light ratio. (6) Any dark halo of the main lens galaxy must be truncated near 1.5 arcsec (4 h^{-1} kpc) before the inferred Ho rises above 60 km/s/Mpc. (7) The quasar host galaxy is lensed into an Einstein ring connecting the four quasar images, whose shape is reproduced by the models. Improved NICMOS imaging of the ring could be used to break the degeneracy of the lens models.

We present a high-resolution spectrum of the high redshift, lensed quasar Q1208+1101, obtained with the echellette spectrograph on the Multiple Mirror Telescope. We examine the new and published spectra and provide an updated list of high-confidence metal-line absorption systems at z=1.1349, 2.8626, 2.9118, 2.9136, 2.9149. Combining this with a simple model of the gravitational lens system allows us to constrain the possible lens redshifts. The high-redshift (z > 2.5) and low-redshift (z < 0.4) candidates can be ruled out with high confidence. The current spectra effectively probe about 40% of the redshift range in which the lens is expected. In that range, there is only one known metal-line absorption system, an MgII absorber at z=1.1349. We consider the possibility that this system is the lensing galaxy and discuss the implied parameters of the galaxy.

We combine photometry and lens modeling to study the properties of 17 gravitational lens galaxies between z=0.1 and 1. Most of the lens galaxies are passively evolving early-type galaxies, with a few spirals. The colors, scale lengths, and ellipticities of lens galaxies are similar to those of the general population of early-type galaxies, although there may be a deficit of apparently round lens galaxies produced by the inclination dependence of lensing cross sections. The projected mass distributions are aligned with the projected light distributions to <~ 10 deg, except in the presence of a strong external tidal perturbation, suggesting that dark matter halos have orbits that are significantly modified by interactions with the baryonic component and are not far out of alignment with the stars. Lens galaxies obey image separation/lens luminosity correlations analogous to the Faber-Jackson and Tully-Fisher relations, which are consistent with standard dark matter lens models. The lens galaxy mass-to-light ratios decrease with redshift as d(log M/L_B)/dz = -0.3\pm0.1 (-0.5\pm0.1) for Omega_0=1 (0.1), thus providing direct evidence of passive evolution for a sample of early-type galaxies in low-density environments. The evolution-corrected mass-to-light ratios are generally larger than predicted by constant M/L dynamical models, although there is significant scatter; with improved photometry, lens galaxy mass-to-light ratios would better distinguish between constant M/L and dark matter models. These conclusions are limited primarily by the quality of lens galaxy photometry.

Time delay measurements have recently been reported for the lensed quasar PG 1115+080. These measurements can be used to derive Ho, but only if we can constrain the lensing potential. We have applied a recently developed deconvolution technique to analyse sub-arcsecond I band images of PG 1115+080, obtained at the Nordic Optical Telescope (NOT) and the Canada France Hawaii Telescope (CFHT). The high performance of the deconvolution code allows us to derive precise positions and magnitudes for the four lensed images of the quasar, as well as for the lensing galaxy. The new measurement of the galaxy position improves its precision by a factor of 3 and thus strengthens the constraints on the lensing potential. With the new data, a range of models incorporating some of the plausible systematic uncertainties yields Ho = 53 (+10/-7) km/s/mpc.

We study gravitational lensing by spiral galaxies, using realistic models consisting of halo, disk, and bulge components combined to produce a flat rotation curve. Proper dynamical normalization of the models is critical because a disk requires less mass than a spherical halo to produce the same rotation curve---a face-on Mestel disk has a lensing cross section only 41% as large as a singular isothermal sphere with the same rotation curve. The cross section is sensitive to inclination and dominated by edge-on galaxies, which produce lenses with an unobserved 2-image geometry and a smaller number of standard 5-image lenses. Unless the disk is unreasonably massive, disk+halo models averaged over inclination predict \lesssim 10% more lenses than pure halo models. Finally, models with an exponential disk and a central bulge are sensitive to the properties of the bulge. In particular, an exponential disk model normalized to our Galaxy cannot produce multiple images without a bulge, and including a bulge reduces the net flattening of edge-on galaxies. The dependence of the lensing properties on the masses and shapes of the halo, disk, and bulge means that a sample of spiral galaxy lenses would provide useful constraints on galactic structure.

For the quadruple gravitational lens PG 1115+080, we combine recent measurements of the time delays with new lens models to determine the Hubble constant H_0. We explore the effects of systematic uncertainties in the lens models on the estimates of H_0, and we discuss how the uncertainties can be reduced by future observations. We find that the lens cannot be fit by an isolated lens galaxy, but that it can be well fit by including a perturbation from the nearby group of galaxies. To understand the full range of systematic uncertainties it is crucial to use an ellipsoidal galaxy and to let the group position vary. In this case, the existing constraints cannot break degeneracies in the models with respect to the profiles of the galaxy and group and to the position of the group. Combining the known time delays with a range of lens models incorporating most of the plausible systematic effects yields H_0 = 51_{-13}^{+14} km s^{-1} Mpc^{-1}. The constraints on the lens models, and hence on H_0, can be improved by reducing the standard errors in the lens galaxy position from 50 mas to \sim10 mas, reducing the uncertainties in the time delays to \sim0.5 days, and constraining the lens mass distribution using HST photometry and the fundamental plane. In particular, the time delay ratio r_{ABC} = \Delta\tau_{AC} / \Delta\tau_{BA} may provide the best constraint on the mass profile of the galaxy.

Gravitational lenses are a unique new constraint on the structure of galaxies. We review the evidence that most lenses are early-type galaxies, the optical properties of the lens galaxies, the evidence against constant M/L models, recent work on the axis ratios of the mass distribution, and the role stellar dynamics plays in gravitational lensing.

Galaxies modeled as singular isothermal ellipsoids with an axis ratio distribution similar to the observed axis ratio distribution of E and S0 galaxies are statistically consistent with both the observed numbers of two-image and four-image lenses and the inferred ellipticities of individual lenses. However, no four-image lens is well fit by the model (typical $\chi^2/N_{dof} \sim 20$), the axis ratio of the model can be significantly different from that of the observed lens galaxy, and the major axes of the model and the galaxy may be slightly misaligned. We found that models with a second, independent, external shear axis could fit the data well (typical $\chi^2/N_{dof} \sim 1$), while adding the same number of extra parameters to the radial mass distribution does not produce such a dramatic improvement in the fit. An independent shear axis can be produced by misalignments between the luminous galaxy and its dark matter halo, or by external shear perturbations due to galaxies and clusters correlated with the primary lens or along the line of sight. We estimate that the external shear perturbations have no significant effect on the expected numbers of two-image and four-image lenses, but that they can be important perturbations in individual lens models. However, the amplitudes of the external shears required to produce the good fits are larger than our estimates for typical external shear perturbations (10-15% shear instead of 1-3% shear) suggesting that the origin of the extra angular structure must be intrinsic to the primary lens galaxy in most cases.