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

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

We report the quadruple nature of the source WISE 025942.9-163543 as observed in the VST-ATLAS survey. Spectra of the two brightest images show quasar emission lines at z=2.16. The system was discovered by splitting ATLAS "cutouts" of WISE sources with W1-W2 > 0.7, when possible, into three components. Followup Magellan images were used to obtain astrometry and g and i photometry, with i=18.78 and 19.73, respectively, for the brightest and faintest components. Absorption lines are observed at z=0.905 but there is little evidence for a lensing galaxy after PSF fitting and subtraction. Saha and Williams (2003) would classify ATLAS 0259-1635 as a short-axis quad. The larger and smaller diameters are 1.57 and 1.32 arcseconds, respectively. Modeling the lensing galaxy as a singular isothermal sphere with external shear, the largest and smallest predicted magnifications are 10.8 and 6.4 respectively.

We report on discovery results from a quasar lens search in the ATLAS public footprint, extending quasar lens searches to a regime without $u-$band or fiber-spectroscopic information, using a combination of data mining techniques on multi-band catalog magnitudes and image-cutout modelling. Spectroscopic follow-up campaigns, conducted at the 2.6m Nordic Optical Telescope (La Palma) and 3.6m New Technology Telescope (La Silla) in 2016, yielded seven pairs of quasars exhibiting the same lines at the same redshift and monotonic flux-ratios with wavelength (hereafter NIQs, Nearly Identical Quasar pairs). The quasar redshifts range between $\approx1.2$ and $\approx 2.7;$ contaminants are typically pairs of bright blue stars, quasar-star alignments along the line of sight, and narrow-line galaxies at $0.3<z<0.7.$ Magellan data of A0140-1152 (01$^h$40$^m$03.0$^s$-11$^d$52$^m$19.0$^s$, $z_{s}=1.807$) confirm it as a lens with deflector at $z_{l}=0.277$ and Einstein radius $\theta_{\rm E}=(0.73\pm0.02)^\ase$. We show the use of spatial resolution from the Gaia mission to select lenses and list additional systems from a WISE-Gaia-ATLAS search, yielding three additional lenses (02$^h$35$^m$27.4$^s$-24$^d$33$^m$13.2$^s$, 02$^h$59$^m$33.$^s$-23$^d$38$^m$01.8$^s$, 01$^h$46$^m$32.9$^s$-11$^d$33$^m$39.0$^s$). The overall sample consists of 11 lenses/NIQs, plus three lenses known before 2016, over the ATLAS-DR3 footprint ($\approx3500$~deg$^2$). Finally, we discuss future prospects for objective classification of pair/NIQ/contaminant spectra.

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

We have analyzed images from the VST ATLAS survey to identify candidate gravitationally lensed quasar systems in a sample of WISE sources with W1 - W2 > 0.7. Results from followup spectroscopy with the Baade 6.5 m telescope are presented for eight systems. One of these is a quadruply lensed quasar and two are doubly lensed systems. Two are projected superpositions of two quasars at different redshifts. In one system two quasars, though at the same redshift, have very different emission line profiles, and constitute a physical binary. In two systems the component spectra are consistent with the lensing hypothesis, after allowing for micro-lensing. But as no lensing galaxy is detected in these two, we classify them as lensless twins. More extensive observations are needed to establish whether they are in fact lensed quasars or physical binaries.

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

We present an analysis of observations of the doubly-lensed double hotspot in the giant radio galaxy J0816+5003 from MERLIN, MDM, WIYN, WHT, UKIRT and the VLA. The images of the two hotspot components span a factor of two in radius on one side of the lensing galaxy at impact parameters of less than 500pc. Hence we measure the slope of the lensing potential over a large range in radius, made possible by significant improvement in the accuracy of registration of the radio and optical frame and higher resolution imaging data than previously available. We also infer the lens and source redshifts to be 0.332 and > 1 respectively. Purely on the basis of lens modelling, and independently of stellar velocity dispersion measurements, we find the potential to be very close to isothermal.

We use the microlensing variability observed for eleven gravitationally lensed quasars to show that the accretion disk size at a rest-frame wavelength of 2500 Angstroms is related to the black hole mass by log(R_{2500}/cm)=(15.78\pm0.12) + (0.80\pm0.17)\log(M_BH/10^9M_sun). This scaling is consistent with the expectation from thin disk theory (R ~ M_BH^{2/3}), but when interpreted in terms of the standard thin disk model (T ~ R^{-3/4}), it implies that black holes radiate with very low efficiency, log(eta) = -1.77\pm0.29 + log(L/L_E) where eta=L/(Mdot*c^2). Only by making the maximum reasonable shifts in the average inclination, Eddington factors and black hole masses can we raise the efficiency estimate to be marginally consistent with typical efficiency estimates (eta ~ 10%). With one exception, these sizes are larger by a factor of ~4 than the size needed to produce the observed 0.8 micron quasar flux by thermal radiation from a thin disk with the same T ~ R^{-3/4} temperature profile. While scattering a significant fraction of the disk emission on large scales or including a large fraction of contaminating line emission can reduce the size discrepancy, resolving it also appears to require that accretion disks have flatter temperature/surface brightness profiles.

We report on the redshift of the lensing galaxy and of the quasar QJ 0158-4325 and on the lens model of the system. A deep VLT/FORS2 spectrum and HST/NICMOS-F160W images are deconvolved. From the images we derive the light profile of the lensing galaxy and an accurate relative astrometry for the system. In addition we measure the flux ratio between the quasar images in the MgII emission line to constrain the mass model. From the spectrum we measure the redshift of the lensing galaxy (z=0.317+/-0.001) and of the quasar (z=1.294+/-0.008). Using the flux ratio in the lens model allows to discard the SIE as a suitable approximation of the lens potential. On the contrary the truncated-PIEMD gives a good fit to the lens and leads to a time delay of t(A-B)=-14.5+/-0.1 days, with H0=73 km/s/Mpc. Using the flux ratio to constrain the mass model favors the truncated-PIEMD over the SIE, while ignoring this constraint leaves the choice open.

We analyzed the microlensing of the X-ray and optical emission of the lensed quasar PG 1115+080. We find that the effective radius of the X-ray emission is 1.3(+1.1 -0.5) dex smaller than that of the optical emission. Viewed as a thin disk observed at inclination angle i, the optical accretion disk has a scale length, defined by the point where the disk temperature matches the rest frame energy of the monitoring band (kT=hc/lambda_rest with lambda_rest=0.3 micron), of log[(r_{s,opt}/cm)(cos(i) / 0.5)^{1/2}] = 16.6 \pm 0.4. The X-ray emission region (1.4-21.8 keV in the rest frame) has an effective half-light radius of log[r_{1/2,X}/cm] = 15.6 (+0.6-0.9}. Given an estimated black hole mass of 1.2 * 10^9 M_sun, corresponding to a gravitational radius of log[r_g/cm] = 14.3, the X-ray emission is generated near the inner edge of the disk while the optical emission comes from scales slightly larger than those expected for an Eddington-limited thin disk. We find a weak trend supporting models with low stellar mass fractions near the lensed images, in mild contradiction to inferences from the stellar velocity dispersion and the time delays.

We expand our Bayesian Monte Carlo method for analyzing the light curves of gravitationally lensed quasars to simultaneously estimate time delays and quasar structure including their mutual uncertainties. We apply the method to HE1104-1805 and QJ0158-4325, two doubly-imaged quasars with microlensing and intrinsic variability on comparable time scales. For HE1104-1805 the resulting time delay of (Delta t_AB) = t_A - t_B = 162.2 -5.9/+6.3 days and accretion disk size estimate of log(r_s/cm) = 15.7 -0.5/+0.4 at 0.2 micron in the rest frame are consistent with earlier estimates but suggest that existing methods for estimating time delays in the presence of microlensing underestimate the uncertainties. We are unable to measure a time delay for QJ0158-4325, but the accretion disk size is log(r_s/cm) = 14.9 +/- 0.3 at 0.3 micron in the rest frame.

We analyze V, I and H band HST images and two seasons of R-band monitoring data for the gravitationally lensed quasar SDSS0924+0219. We clearly see that image D is a point-source image of the quasar at the center of its host galaxy. We can easily track the host galaxy of the quasar close to image D because microlensing has provided a natural coronograph that suppresses the flux of the quasar image by roughly an order of magnitude. We observe low amplitude, uncorrelated variability between the four quasar images due to microlensing, but no correlated variations that could be used to measure a time delay. Monte Carlo models of the microlensing variability provide estimates of the mean stellar mass in the lens galaxy (0.02 Msun < M < 1.0 Msun), the accretion disk size (the disk temperature is 5 x 10^4 K at 3.0 x 10^14 cm < rs < 1.4 x 10^15 cm), and the black hole mass (2.0 x 10^7 Msun < MBH \eta_{0.1}^{-1/2} (L/LE)^{1/2} < 3.3 x 10^8 Msun), all at 68% confidence. The black hole mass estimate based on microlensing is consistent with an estimate of MBH = 7.3 +- 2.4 x 10^7 Msun from the MgII emission line width. If we extrapolate the best-fitting light curve models into the future, we expect the the flux of images A and B to remain relatively stable and images C and D to brighten. In particular, we estimate that image D has a roughly 12% probability of brightening by a factor of two during the next year and a 45% probability of brightening by an order of magnitude over the next decade.

We use the microlensing variability observed for nine gravitationally lensed quasars to show that the accretion disk size at 2500 Angstroms is related to the black hole mass by log(R_2500/cm) = (15.6+-0.2) + (0.54+-0.28)log(M_BH/10^9M_sun). This scaling is consistent with the expectation from thin disk theory (R ~ M_BH^(2/3)), but it implies that black holes radiate with relatively low efficiency, log(eta) = -1.29+-0.44 + log(L/L_E) where eta=L/(Mdot c^2). These sizes are also larger, by a factor of ~3, than the size needed to produce the observed 0.8 micron quasar flux by thermal radiation from a thin disk with the same T ~ R^(-3/4) temperature profile. More sophisticated disk models are clearly required, particularly as our continuing observations improve the precision of the measurements and yield estimates of the scaling with wavelength and accretion rate.

Based on the microlensing variability of the two-image gravitational lens HE1104-1805 observed between 0.4 and 8 microns, we have measured the size and wavelength-dependent structure of the quasar accretion disk. Modeled as a power law in temperature, T proportional to R^-beta, we measure a B-band (0.13 microns in the rest frame) half-light radius of R_{1/2,B} = 6.7 (+6.2 -3.2) x 10^15 cm (68% CL) and a logarithmic slope of beta=0.61 (+0.21 -0.17) for our standard model with a logarithmic prior on the disk size. Both the scale and the slope are consistent with simple thin disk models where beta=3/4 and R_{1/2,B} = 5.9 x 10^15 cm for a Shakura-Sunyaev disk radiating at the Eddington limit with 10% efficiency. The observed fluxes favor a slightly shallower slope, beta=0.55 (+0.03 -0.02), and a significantly smaller size for beta=3/4.

The mid-IR flux ratios F_A/F_B = 2.84 +/- 0.06 of the two images of the gravitationally lensed quasar HE 1104-1805 show no wavelength dependence to within 3% across 3.6-8.0 um, no time dependence over 6 months and agree with the broad emission line flux ratios. This indicates that the mid-IR emission likely comes from scales large enough to be little affected by microlensing and that there is little differential extinction between the images. We measure a revised time-delay between these two images of 152.2 +2.8-3.0 days from R and V-band data covering 1997 to 2006. This time-delay indicates that the lens has an approximately flat rotation curve over scales of 1-2 R_e. We also observed uncorrelated variations of ~0.05 mag/yr which we attribute to microlensing of the optical emission from the accretion disk. The optical colors have also changed significantly in the sense that image A is now redder than image B, rather than bluer as it was in 1993.

We stacked the X-ray data from the ROSAT All Sky Survey for over 4,000 clusters selected from the 2MASS catalog and divided into five richness classes. We detected excess X-ray emission over background at the center of the stacked images in all five richness bins. The interrelationships between the mass, X-ray temperature and X-ray luminosity of the stacked clusters agree well with those derived from catalogs of X-ray clusters. Poisson variance in the number of galaxies occupying halos of a given mass leads to significant differences between the average richness at fixed mass and the average mass at fixed richness that we can model relatively easily using a simple model of the halo occupation distribution. These statistical effects probably explain recent results in which optically-selected clusters lie on the same X-ray luminosity-temperature relations as local clusters but have lower optical richnesses than observed for local clusters with the same X-ray properties. When we further binned the clusters by redshift, we did not find significant redshift-dependent biases in the sense that the X-ray luminosities for massive clusters of fixed optical richness show little dependence on redshift beyond that expected from the effects of Poisson fluctuations. Our results demonstrate that stacking of RASS data from optically selected clusters can be a powerful test for biases in cluster selection algorithms.

We obtained 98 R-band and 18 B, r', i' images of the optical afterglow of GRB 060526 (z=3.21) with the MDM 1.3m, 2.4m, and the PROMPT telescopes in Cerro Tololo over the 5 nights following the burst trigger. Combining these data with other optical observations reported in GCN and the Swift-XRT observations, we compare the optical and X-ray afterglow light curves of GRB 060526. Both the optical and X-ray afterglow light curves show rich features, such as flares and breaks. The densely sampled optical observations provide very good coverage at T>1.e4 sec. We observed a break at 2.4e5 sec in the optical afterglow light curve. Compared with the X-ray afterglow light curve, the break is consistent with an achromatic break supporting the beaming models of GRBs. However, the pre-break and post-break temporal decay slopes are difficult to explain in simple afterglow models. We estimated a jet angle of \theta_j ~ 7 degrees and a prompt emission size of R_{prompt} ~ 2e14 cm. In addition, we detected several optical flares with amplitudes of \Delta m ~ 0.2, 0.6, and 0.2 mag. The X-ray afterglows detected by Swift have shown complicated decay patterns. Recently, many well-sampled optical afterglows also show decays with flares and multiple breaks. GRB 060526 provides an additional case of such a complex, well observed optical afterglow. The accumulated well-sampled afterglows indicate that most of the optical afterglows are complex.

Microlensing of gravitationally lensed quasars by the stars in the foreground lens galaxy can be used to probe the nature of dark matter, to determine the mean stellar mass in the lens galaxy, and to measure the internal structure of quasar accretion disks. Until recently, little progress has been made toward using microlensing for these purposes because of the difficulty in obtaining the necessary data and the lack of good analysis methods. In the last few years, both problems have been solved. In particular, Bayesian analysis methods provide a general approach to measuring quantities of physical interest and their uncertainties from microlensing light curves. We discuss the data and the analysis methods and show preliminary results for all three astrophysical applications.

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.

We have measured the three time delays for the quadruple gravitational lens RX J1131-1231 using two seasons of monitoring data. The short delays between cusp images are A-B = 11.98 +1.52,-1.27 days and A-C = 9.61 +1.97,-1.57 days. The long A-D delay for the counter image is not as precisely determined because of the season gaps, but the data suggest a delay of -87 +- 8 days. The short delays are difficult to explain using standard isothermal halo models of the lensing potential, which instead prefer A-B and A-C delays of ~1 day for reasonable values of the Hubble constant. Matching the cusp delays is possible by adding a significant (~5 x 10^10 M_sun) amount of matter nearly coincident (~0.05" South-East) with the A image. Adding such a satellite also helps improve the quasar and lens astrometry of the model, reduces the velocity dispersion of the main lens and shifts it closer to the Fundamental Plane. This is suggestive of a satellite galaxy to the primary lens, but its expected luminosity and proximity to both image A and the system's bright Einstein ring make visual identification impossible even with the existing HST data. We also find evidence for significant structure along the line of sight toward the lens. Archival Chandra observations show two nearby regions of extended X-ray emission, each with bolometric X-ray luminosities of 2-3 x 10^43 ergs/s. The brighter region is located 153" from the lens and centered on a z=0.1 foreground cD galaxy, and the fainter and presumably more distant region is 4-5 times closer (in angular separation) to the lens and likely corresponds to the weaker of two galaxy red sequences (which includes the lens galaxy) previously detected at optical wavelengths.

We demonstrate that gravitationally lensed quasars are easily recognized using image subtraction methods as time variable sources that are spatially extended. For Galactic latitudes |b|>20 deg, lensed quasars dominate the population of spatially extended variable sources, although there is some contamination from variable star pairs, variable star-quasar pairs and binary quasars that can be easily controlled using other information in the survey such as the object light curves and colors. This will allow planned large-scale synoptic surveys to find lensed quasars almost down to their detection limits without the need for extensive follow-up observations.

We present Hubble Space Telescope images and 2 years of optical photometry of the quadruple quasar HE0435-1223. The time delays between the intrinsic quasar variations are 14.4+/-0.8 (A-D), 8.0+/-0.8 (A-B) and 2.1+/-0.8 (A-C) days. We also observed non-intrinsic variations of ~0.1 mag/yr that we attribute to microlensing. Instead of the traditional approach of assuming a rotation curve for the lens galaxy and then deriving the Hubble constant (H_0), we assume H_0=(72+/-7) km/s/Mpc and derive constraints on the rotation curve. On the scale over which the lensed images occur (1.2"=5kpc/h=1.5R_e), the lens galaxy must have a rising rotation curve, and it cannot have a constant mass-to-light ratio. These results add to the evidence that the structures of early-type galaxies are heterogeneous.

The redshift of the galaxy lensing HE 0435-1223 is 0.4546 +/- 0.0002, based on observations obtained with the Low Dispersion Survey Spectrograph 2 (LDSS2) on the Magellan Consortium's 6.5 m Clay telescope. HST/ACS observations of the system also reveal a spiral-rich group of 10 galaxies within 40" of the elliptical lensing galaxy. The redshifts for two of these galaxies were measured to be in the foreground (at z=0.419) with respect to the lens, thus at least some of the nearby galaxies are not part of the same physical group as the lensing galaxy. Mass models of the system (assuming same-plane deflectors) that take the local group environment into account do better at explaining the observed emission-line flux ratios (which are presumably unaffected by microlensing) than single halo models, but the match is still not perfect. In particular, component A (a minimum of the light travel time) is observed to be 0.20 mag brighter than predicted and component C (also a minimum image) is observed to be 0.16 mag fainter than predicted. Mass models for the system predict an A-D time delay of either 15.8 or 17.6 days (Ho = 72 km/s/Mpc) depending on the details of the local galaxy environment.

Exploration of the time variability on the sky over a broad range of flux levels and wavelengths is rapidly becoming a new frontier of astronomical research. We describe here briefly the Palomar-QUEST survey being carried out from the Samuel Oschin 48-inch Schmidt telescope at Palomar. The following features make the survey an attractive candidate for studying time variability: anticipated survey area of 12,000 - 15,000 sq. degrees in the drift scan mode, point source depth of 21st mag. in I under good conditions, near simultaneous observations in four filters, and at least four passes per year at each location covered. The survey will yield a large number of transients and highly variable sources in the near future and in that sense is a prototype of LSST and Pan-STARRS. We briefly outline our strategy for searching such objects and the proposed pipeline for detecting transients in real-time.

The Supernova / Acceleration Probe (SNAP) is a proposed space-based experiment designed to study the dark energy and alternative explanations of the acceleration of the Universe's expansion by performing a series of complementary systematics-controlled measurements. We describe a self-consistent reference mission design for building a Type Ia supernova Hubble diagram and for performing a wide-area weak gravitational lensing study. A 2-m wide-field telescope feeds a focal plane consisting of a 0.7 square-degree imager tiled with equal areas of optical CCDs and near infrared sensors, and a high-efficiency low-resolution integral field spectrograph. The SNAP mission will obtain high-signal-to-noise calibrated light-curves and spectra for several thousand supernovae at redshifts between z=0.1 and 1.7. A wide-field survey covering one thousand square degrees resolves ~100 galaxies per square arcminute. If we assume we live in a cosmological-constant-dominated Universe, the matter density, dark energy density, and flatness of space can all be measured with SNAP supernova and weak-lensing measurements to a systematics-limited accuracy of 1%. For a flat universe, the density-to-pressure ratio of dark energy can be similarly measured to 5% for the present value w0 and ~0.1 for the time variation w'. The large survey area, depth, spatial resolution, time-sampling, and nine-band optical to NIR photometry will support additional independent and/or complementary dark-energy measurement approaches as well as a broad range of auxiliary science programs. (Abridged)

The advent of wide-area multicolour synoptic sky surveys is leading to data sets unprecedented in size, complexity and data throughput. VO technology offers a way to exploit these to the full but requires changes in design philosophy. The Palomar-QUEST survey is a major new survey being undertaken by Caltech, Yale, JPL and Indiana University to repeatedly observe 1/3 of the sky (~15000 sq. deg. between -27 < Dec <27 in seven passbands. Utilising the 48-inch Oschin Schmidt Telescope at the Palomar Observatory with the 112-CCD QUEST camera covering the full 4 x 4 sq. deg. field of view, it will generate \~1TB of data per month. In this paper, we review the design of QUEST as a VO resource, a federated data set and an exemplar of VO standards.

We report the discovery of two new gravitationally lensed quasars, WFI J2026-4536 and WFI J2033-4723, at respective source redshifts of z=2.23 and z=1.66. Both systems are quadruply imaged and have similar PG1115-like image configurations. WFI J2026-4536 has a maximum image separation of 1.4", a total brightness of g = 16.5, and a relatively simple lensing environment, while WFI J2033-4723 has a maximum image separation of 2.5", an estimated total brightness of g = 17.9, and a more complicated environment of at least six galaxies within 20". The primary lensing galaxies are detected for both systems after PSF subtraction. Several of the broadband flux ratios for the two lenses show a strong (0.1-0.4 mags) trend with wavelength, suggesting either microlensing or differential extinction through the lensing galaxy. For WFI J2026-4536, the total quasar flux has dimmed by 0.1 mag in the blue but only half as much in the red over three months, suggestive of microlensing-induced variations. For WFI J2033-4723, resolved spectra of some of the quasar components reveal emission line flux ratios that agree better with the macromodel predictions than either the broadband or continuum ratios, also indicative of microlensing. The predicted differential time delays for WFI J2026-4536 are short, ranging from 1-2 weeks for the long delay, but are longer for WFI J2033-4723, ranging from 1-2 months. Both systems hold promise for future monitoring campaigns aimed at microlensing or time delay studies.

We report that the Sloan Digital Sky Survey quasar SDSS J165043.44+425149.3 is gravitationally lensed into two images, based on observations obtained with the WIYN 3.5 m telescope at the Kitt Peak National Observatory. The lensed quasar, at a redshift of z=1.54, appears as two images separated by 1.2" with B-band magnitudes of 17.8 and 20.0. The lensing galaxy is clearly detected in I-band images obtained in 0.3" seeing, after point spread function subtraction of the two quasar images. A strong metal-line absorption system is also identified in the unresolved SDSS spectrum of the double quasar, suggesting a plausible lens redshift of z=0.58. The UBRI flux ratios of the pair vary significantly from 8.5:1 in the blue to 5.4:1 in the red, a difference of 0.5 magnitudes, and may indicate the presence of microlensing in one or both quasar components. The predicted differential time delay between quasar images is on the order of one month, assuming the intervening absorption system is due to the lensing galaxy.

We present results from monitoring of the distant (z = 2.64), gravitationally lensed quasar MG J0414+0534 with the Chandra X-ray Observatory. An Fe Ka line at 6.49 +/- 0.09 keV (rest-frame) with an equivalent width of ~ 190eV consistent with fluorescence from a cold medium is detected at the 99 percent confidence level in the spectrum of the brightest image A. During the last two observations of our monitoring program we detected a five-fold increase of the equivalent width of a narrow Fe Ka line in the spectrum of image B but not in the brighter image A whereas image C is too faint to resolve the line. The continuum emission component of image B did not follow the sudden enhancement of the iron line in the last two observations. We propose that the sudden increase in the iron line strength from ~ 190eV to 900eV only in image B can be explained with a caustic crossing due to microlensing that selectively enhances a strip of the line emission region of the accretion disk. The non-enhancement of the continuum emission in the spectrum of image B suggests that the continuum emission region is concentrated closer to the center of the accretion disk than the iron line emission region and the magnification caustic has not reached close enough to the former region to amplify it. A model of a caustic crossing event predicts discontinuities in the light-curve of the magnification and provides an upper limit of ~ 5 x 10^(-4) pc on the outer radius of the Fe Ka emission region. The non-detection of any relativistic or Doppler shifts of the iron line in the spectrum of image B implies that the magnification caustic for the last two observations was located at a radius greater than ~ 100 gravitational radii.

We report the discovery of a gravitationally lensed quasar resulting from our survey for lenses in the southern sky. Radio images of PMN J1632-0033 with the VLA and ATCA exhibit two compact, flat-spectrum components with separation 1.47" and flux density ratio 13.2. Images with the HST reveal the optical counterparts to the radio components and also the lens galaxy. An optical spectrum of the bright component, obtained with the first Magellan telescope, reveals quasar emission lines at redshift 3.42. Deeper radio images with MERLIN and the VLBA reveal a faint third radio component located near the center of the lens galaxy, which is either a third image of the background quasar or faint emission from the lens galaxy.

The mini Broad Absorption Line (BAL) quasar RX J0911.4+0551 was observed with the Advanced CCD Imaging Spectrometer (ACIS) of the Chandra X-ray Observatory for ~ 29 ks as part of a gravitational lens (GL) survey aimed at measuring time-delays. Timing analysis of the light-curve of the lensed image A2 shows a rapid flux variation with a duration of about 2000s. A Kolmogorov-Smirnov test shows that the probability that a constant-intensity source would produce the observed variability is less than ~ 0.2 percent. We discuss possible origins for the observed short-term X-ray variability. Our gravitational lens models for the RX J0911.4+0551 GL system predict a time-delay of less than a day between images A1 and A2. The rapid variability combined with the predicted short-time delay make RX J0911.4+0551 an ideal system to apply the GL method for estimating the Hubble constant. We describe the prospects of measuring H_0 within single X-ray observations of GL systems with relatively short time delays. Modeling of the spectrum of the mini-BAL quasar RX J0911.4+0551 suggests the presence of an intrinsic absorber. Partial covering models are slightly preferred over models that contain absorption due to intrinsic ionized or neutral gas.

We present results from X-ray observations of the quadruply lensed quasar RX J0911.4+0551 using data obtained with the Advanced CCD Imaging Spectrometer (ACIS) on board the Chandra X-ray Observatory. The 29 ks observation detects a total of ~404 X-ray photons (0.3 to 7.0 keV) from the four images of the lensed quasar. Deconvolution of the aspect corrected data resolves all four lensed images, with relative positions in good agreement with optical measurements. When compared to contemporaneous optical data, one of the lensed images (component A3) is dimmer by a factor of ~6 in X-rays with respect to the 2 brighter images (components A1 and A2). Spectral fitting for the combined images shows significant intrinsic absorption in the soft (0.2 to 2.4 keV) energy band, consistent with the mini-BAL nature of this quasar, while a comparison with ROSAT PSPC observations from 1990 shows a drop of ~6.5 in the total soft bandpass flux. The observations also detect ~157 X-ray photons arising from extended emission of the nearby cluster (peaked ~42" SW of RXJ0911.4+0551) responsible for the large external shear present in the system. The Chandra observation reveals the cluster emission to be complex and non-spherical, and yields a cluster temperature of kT = 2.3^{+1.8}_{-0.8} keV and a 2.0 to 10 keV cluster luminosity within a 1 Mpc radius of L_X = 7.6_{-0.2}^{+0.6} x 10^{43} ergs/s (error bars denote 90% confidence limits). Our mass estimate of the cluster within its virial radius is 2.3^{+1.8}_{-0.7} x 10^{14} solar, and is a factor of 2 smaller than, although consistent with, previous mass estimates based on the observed cluster velocity dispersion.

We present our ground-based optical imaging, spectral analysis, and high resolution radio mapping of the gravitational lens candidate FBQ 1633+3134. This z=1.52, B=17.7 quasar appears double on CCD images with an image separation of 0.66 arcseconds and a flux ratio of ~3:1 across BVRI filters. A single 0.27 mJy radio source is detected at 8.46 GHz, coincident to within an arcsecond of both optical components, but no companion at radio wavelengths is detected down to a flux level of 0.1 mJy (3 sigma). Spectral observations reveal a rich metal-line absorption system consisting of a strong Mg II doublet and associated Fe I and Fe II absorption features, all at an intervening redshift of z=0.684, suggestive of a lensing galaxy. Point spread function subtraction however shows no obvious signs of a third object between the two quasar images, and places a detection limit of I > 23.0 if such an object exists. Although the possibility that FBQ 1633+3134 is a binary quasar cannot be ruled out, the evidence is consistent with it being a single quasar lensed by a faint, metal-rich galaxy.

We report the discovery of a new double-image quasar that was found during a search for gravitational lenses in the southern sky. Radio source PMN J1838-3427 is composed of two flat-spectrum components with separation 1", flux density ratio 14:1 and matching spectral indices, in VLA and VLBA images. Ground-based BRI images show the optical counterpart (total I=18.6) is also double with the same separation and position angle as the radio components. An HST/WFPC2 image reveals the lens galaxy. The optical flux ratio (27:1) is higher than the radio value probably due to differential extinction of the components by the lens galaxy. An optical spectrum of the bright component contains quasar emission lines at z=2.78 and several absorption features, including prominent Ly-alpha absorption. The lens galaxy redshift could not be measured but is estimated to be z=0.36 +/- 0.08. The image configuration is consistent with the simplest plausible models for the lens potential. The flat radio spectrum and observed variability of PMN J1838-3427 suggest the time delay between flux variations of the components is measurable, and could thus provide an independent measurement of H_0.

We report the discovery of the new double quasar CTQ 839. This B = 18.3, radio quiet quasar pair is separated by 2.1" in BRIH filters with magnitude differences of delta m_B = 2.5, delta m_R = delta m_I = 1.9, and delta m_H = 2.3. Spectral observations reveal both components to be z = 2.24 quasars, with relative redshifts that agree at the 100 km/s level, but exhibit pronounced differences in the equivalent widths of related emission features, as well as an enhancement of blue continuum flux in the brighter component longward of the Ly alpha emission feature. In general, similar redshift double quasars can be the result of a physical binary pair, or a single quasar multiply imaged by gravitational lensing. Empirical PSF subtraction of R and H band images of CTQ 839 reveal no indication of a lensing galaxy, and place a detection limit of R = 22.5 and H = 17.4 for a third component in the system. For an Einstein-de Sitter cosmology and SIS model, the R band detection limit constrains the characteristics of any lensing galaxy to z_lens >= 1 with a corresponding luminosity of L >~ 5 L_*, while an analysis based on the redshift probability distribution for the lensing galaxy argues against the existence of a z_lens >~ 1 lens at the 2 sigma level. A similar analysis for a Lambda dominated cosmology, however, does not significantly constrain the existence of any lensing galaxy. The broadband flux differences, spectral dissimilarities, and failure to detect a lensing galaxy make the lensing hypothesis for CTQ 839 unlikely. The similar redshifts of the two components would then argue for a physical quasar binary. At a projected separation of 8.3/h kpc (Omega_matter = 1), CTQ 839 would be the smallest projected separation binary quasar currently known.

We report the discovery of the new gravitational lens system HE~0230$-$2130, a QSO at redshift $z=2.162$ consisting of at least five distinct components. Three of these are clearly lensed images of the QSO, one is most likely the lensing galaxy, while for the fifth component the identity is unclear: It could be a fourth QSO image (if so, then highly reddened), or another intervening galaxy, or a superposition of the two. Differential reddening seems to be important also for the first three QSO images. The surface density of faint galaxies near the QSO appears to be enhanced by a factor of $\ga 2$, indicating the presence of a distant cluster close to the line of sight.

We report the discovery and ground based observations of the new gravitational lens CTQ 414. The source quasar lies at a redshift of z = 1.29 with a B magnitude of 17.6. Ground based optical imaging reveals two point sources separated by 1.2 arcsec with a magnitude difference of roughly 1 mag. Subtraction of two stellar point spread functions from images obtained in subarcsecond seeing consistently leaves behind a faint, residual object. Fits for two point sources plus an extended object places the fainter object collinear with the two brighter components. Subsequent HST/NICMOS observations have confirmed the identification of the fainter object as the lensing galaxy. VLA observations at 8.46 GHz reveal that all components of the lensing system are radio quiet down to the 0.2 mJy flux level.