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

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

We present a tool to generate mock quasar microlensing light curves and sample them according to any observing strategy. An updated treatment of the fixed and random velocity components of observer, lens, and source is used, together with a proper alignment with the external shear defining the magnification map caustic orientation. Our tool produces quantitative results on high magnification events and caustic crossings, which we use to study three lensed quasars known to display microlensing, viz. RX J1131-1231, HE 0230-2130, and Q 2237+0305, as they would be monitored by The Rubin Observatory Legacy Survey of Space and Time (LSST). We conclude that depending on the location on the sky, the lens and source redshift, and the caustic network density, the microlensing variability may deviate significantly than the expected $\sim$20-year average time scale (Mosquera & Kochanek 2011, arXiv:1104.2356). We estimate that $\sim300$ high magnification events with $\Delta$mag$>1$ mag could potentially be observed by LSST each year. The duration of the majority of high magnification events is between 10 and 100 days, requiring a very high cadence to capture and resolve them. Uniform LSST observing strategies perform the best in recovering microlensing high magnification events. Our web tool can be extended to any instrument and observing strategy, and is freely available as a service at http://gerlumph.swin.edu.au/tools/lsst_generator/, along with all the related code.

We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for this lens system using three band Hubble Space Telescope imaging. We combine the measured time delays, line-of-sight central velocity dispersion of the deflector, and statistically constrained external convergence with our lens models to estimate two cosmological distances. We measure the "effective" time-delay distance corresponding to the redshifts of the deflector and the lensed quasar $D_{\Delta t}^{\rm eff}=3382^{+146}_{-115}$ Mpc and the angular diameter distance to the deflector $D_{\rm d}=1711^{+376}_{-280}$ Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant $H_0 = 74.2^{+2.7}_{-3.0}$ km s$^{-1}$ Mpc$^{-1}$ assuming a flat $\Lambda$CDM cosmology and a uniform prior for $\Omega_{\rm m}$ as $\Omega_{\rm m} \sim \mathcal{U}(0.05, 0.5)$. This measurement gives the most precise constraint on $H_0$ to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analyzed by the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of $H_0$ based on the local distance ladder, reinforcing the tension with the inference from early Universe probes, for example, with 2.2$\sigma$ discrepancy from the cosmic microwave background measurement.

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

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

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

The Hubble Frontier Fields offer an exceptionally deep window into the high-redshift universe, covering a substantially larger area than the Hubble Ultra-Deep field at low magnification and probing 1--2 mags deeper in exceptional high-magnification regions. We aim to leverage Atacama Large Millimetre Array (ALMA) band 6 ($\approx$263\,GHz) mosaics in the central portions of five Frontier Fields to characterize the infrared (IR) properties of $1582$ ultraviolet (UV)-selected Lyman-Break Galaxies (LBGs) at redshifts of $z {\sim}$2--8. We investigated individual and stacked fluxes and IR excess (IRX) values of the LBG sample as functions of stellar mass ($\mathrm{M}_{\bigstar}$), redshift, UV luminosity and slope $\beta$, and lensing magnification. Two (2) LBG candidates were individually detected above a significance of $4.1{-}\sigma$, while stacked samples of the remaining LBG candidates yielded no significant detections. We investigated our detections and upper limits in the context of the IRX-$\mathrm{M}_{\bigstar}$ and IRX-$\beta$ relations, probing at least one dex lower in stellar mass than past studies have done. Our upper limits exclude substantial portions of parameter space and they are sufficiently deep in a handful of cases to create mild tension with the typically assumed attenuation and consensus relations. We observe a clear and smooth trend between $\mathrm{M}_{\bigstar}$ and $\beta$, which extends to low masses and blue (low) $\beta$ values, consistent with expectations from previous works.

Aims: We present and study spatially resolved imaging obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) of multiple $^{12}$CO($J=$6$-$5, 8$-$7 and 9$-$8) and two H$_2$O(2$_{02}-$1$_{11}$ and 2$_{11}-$2$_{02}$) emission lines and cold dust continuum toward the gravitationally lensed dusty star forming galaxy SPT0346-52 at z=$5.656$. Methods: Using a visibility-domain source-plane reconstruction we probe the structure and dynamics of the different components of the interstellar medium (ISM) in this galaxy down to scales of 1 kpc in the source plane. Results: Measurements of the intrinsic sizes of the different CO emission lines indicate that the higher J transitions trace more compact regions in the galaxy. Similarly, we find smaller dust continuum intrinsic sizes with decreasing wavelength, based on observations at rest-frame 130, 300 and 450$\mu$m. The source shows significant velocity structure, and clear asymmetry where an elongated structure is observed in the source plane with significant variations in their reconstructed sizes. This could be attributed to a compact merger or turbulent disk rotation. The differences in velocity structure through the different line tracers, however, hint at the former scenario in agreement with previous [CII] line imaging results. Measurements of the CO line ratios and magnifications yield significant variations as a function of velocity, suggesting that modeling of the ISM using integrated values could be misinterpreted. Modeling of the ISM in SPT0346-52 based on delensed fluxes indicate a highly dense and warm medium, qualitatively similar to that observed in high redshift quasar hosts.

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

Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide information from distinct emission regions in the quasar. Through monitoring of these strong gravitational lenses, a full tomographic view can emerge with Astronomical-Unit scale resolution. Work to date has demonstrated the potential of this technique in about a dozen systems. In the 2020s there will be orders of magnitude more systems to work with. Monitoring of lens systems for caustic-crossing events to enable triggering of multi-platform, multi-wavelength observations in the 2020s will fulfill the potential of quasar microlensing as a unique and comprehensive probe of active black hole structure and dynamics.

Optical and X-ray observations of strongly gravitationally lensed quasars (especially when four separate images of the quasar are produced) determine not only the amount of matter in the lensing galaxy but also how much is in a smooth component and how much is composed of compact masses (e.g., stars, stellar remnants, primordial black holes, CDM sub-halos, and planets). Future optical surveys will discover hundreds to thousands of quadruply lensed quasars, and sensitive X-ray observations will unambiguously determine the ratio of smooth to clumpy matter at specific locations in the lensing galaxies and calibrate the stellar mass fundamental plane, providing a determination of the stellar $M/L$. A modest observing program with a sensitive, sub-arcsecond X-ray imager, combined with the planned optical observations, can make those determinations for a large number (hundreds) of the lensing galaxies, which will span a redshift range of $\sim$$0.25<z<1.5$

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.

Gravitationally lensed quasars are valuable, but extremely rare, probes of observational cosmology and extragalactic astrophysics. Progress in these fields has been limited just by the paucity of systems with good ancillary data. Here we present a first spectroscopic confirmation of lenses discovered in the Southern Sky from the DES and KiDS-DR3 footprints. % optical ground based surveys in the Southern Hemisphere. We have targeted 7 high-graded candidates, selected with new techniques, with NTT-EFOSC2, and confirmed 5 of them. We provide source spectroscopic redshifts, image separations, $gri$ photometry and first lens model parameters. The success rate of ~70% confirms our forecasts, based on the comparison between the number of candidate doubles and quadruplets in our searches over a ~5000 sq.deg footprint and theoretical predictions.

[abridged] Characterizing the number counts of faint, dusty star-forming galaxies is currently a challenge even for deep, high-resolution observations in the FIR-to-mm regime. They are predicted to account for approximately half of the total extragalactic background light at those wavelengths. Searching for dusty star-forming galaxies behind massive galaxy clusters benefits from strong lensing, enhancing their measured emission while increasing spatial resolution. Derived number counts depend, however, on mass reconstruction models that properly constrain these clusters. We estimate the 1.1 mm number counts along the line of sight of three galaxy clusters, i.e. Abell 2744, MACSJ0416.1-2403 and MACSJ1149.5+2223, which are part of the ALMA Frontier Fields Survey. We perform detailed simulations to correct these counts for lensing effects. We use several publicly available lensing models for the galaxy clusters to derive the intrinsic flux densities of our sources. We perform Monte Carlo simulations of the number counts for a detailed treatment of the uncertainties in the magnifications and adopted source redshifts. We find an overall agreement among the number counts derived for the different lens models, despite their systematic variations regarding source magnifications and effective areas. Our number counts span ~2.5 dex in demagnified flux density, from several mJy down to tens of uJy. Our number counts are consistent with recent estimates from deep ALMA observations at a 3$\sigma$ level. Below $\approx$ 0.1 mJy, however, our cumulative counts are lower by $\approx$ 1 dex, suggesting a flattening in the number counts. In our deepest ALMA mosaic, we estimate number counts for intrinsic flux densities $\approx$ 4 times fainter than the rms level. This highlights the potential of probing the sub-10 uJy population in larger samples of galaxy cluster fields with deeper ALMA observations.

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

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

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

Gravitationally lensed quasars are powerful and versatile astrophysical tools, but they are challengingly rare. In particular, only ~25 well-characterized quadruple systems are known to date. To refine the target catalogue for the forthcoming Taipan Galaxy Survey, the images of a large number of sources are being visually inspected in order to identify objects that are confused by a foreground star or galaxies that have a distinct multi-component structure. An unexpected by-product of this work has been the serendipitous discovery of about a dozen galaxies that appear to be lensing quasars, i.e. pairs or quartets of foreground stellar objects in close proximity to the target source. Here we report two diamond-shaped systems. Follow-up spectroscopy with the IMACS instrument on the 6.5m Magellan Baade telescope confirms one of these as a z = 1.975 quasar quadruply lensed by a double galaxy at z = 0.293. Photometry from publicly available survey images supports the conclusion that the other system is a highly sheared quadruply-imaged quasar. In starting with objects thought to be galaxies, our lens finding technique complements the conventional approach of first identifying sources with quasar-like colours and subsequently finding evidence of lensing.

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

Gravitational microlensing is a powerful tool allowing one to probe the structure of quasars on sub-parsec scale. We report recent results, focusing on the broad absorption and emission line regions. In particular microlensing reveals the intrinsic absorption hidden in the P Cygni-type line profiles observed in the broad absorption line quasar H1413+117, as well as the existence of an extended continuum source. In addition, polarization microlensing provides constraints on the scattering region. In the quasar Q2237+030, microlensing differently distorts the H$\alpha$ and CIV broad emission line profiles, indicating that the low- and high-ionization broad emission lines must originate from regions with distinct kinematical properties. We also present simulations of the effect of microlensing on line profiles considering simple but representative models of the broad emission line region. Comparison of observations to simulations allows us to conclude that the H$\alpha$ emitting region in Q2237+030 is best represented by a Keplerian disk.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

[abridged] The Hubble and Spitzer Space Telescope surveys of the Frontier Fields (FF) provide extremely deep images around six massive, strong-lensing clusters of galaxies. The ALMA FF survey aims to cover the same fields at 1.1mm, with maps reaching (unlensed) sensitivities of $<$70$\mu$Jy, in order to explore the properties of background dusty star-forming galaxies. We report on the multi-wavelength photometric analysis of all 12 significantly detected ($>$5$\sigma$) sources in the first three FF clusters observed by ALMA, based on data from Hubble and Spitzer, the VLT and Herschel. We measure the total photometry in all available bands and determine the photometric redshifts and the physical properties of the counterparts via SED-fitting. In particular, we carefully estimate the FIR photometry using 1.1mm priors to limit the misidentification of blended FIR counterparts, which strongly affect some flux estimates in previous FIR catalogs. We identify robust near-infrared (NIR) counterparts for all 11 sources with K$_s$ detection, the majority of which are quite red, with eight having $F814W-K_s\gtrsim 4$ and five having $F160W-[4.5]\gtrsim3$. From the FIR point of view, all our objects have $z_{phot}$$\sim$1--3, whereas based on the optical SED one object prefers a high-$z$ solution ($z\geq\ $7). Five objects among our sample have spectroscopic redshifts from the GLASS survey for which we can reproduce their SEDs with existing templates. This verification confirms the validity of our photometric redshift methodology. The mean redshift of our sample is $z_{phot}$=1.99$\pm$0.27. All 1.1mm selected objects are massive (10.0$<\log[M_{\star}(M_{\odot})]<$ 11.5), with high star formation rates ($<\log[SFR(M_{\odot}/yr)]> \approx$1.6) and high dust contents (8.1 $<\log[M_{dust} (M_{\odot})]<$8.8), consistent with previous ALMA surveys.

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

We present an extensive spectroscopic follow-up campaign of 29 strong lensing (SL) selected galaxy clusters discovered primarily in the Second Red-Sequence Cluster Survey (RCS-2). Our spectroscopic analysis yields redshifts for 52 gravitational arcs present in the core of our galaxy clusters, which correspond to 35 distinct background sources that are clearly distorted by the gravitational potential of these clusters. These lensed galaxies span a wide redshift range of $0.8 \le z \le 2.9$, with a median redshift of $z_s = 1.8 \pm 0.1 $. We also measure reliable redshifts for 1004 cluster members, allowing us to obtain robust velocity dispersion measurements for 23 of these clusters, which we then use to determine their dynamical masses by using a simulation-based $\sigma_{DM} - M_{200}$ scaling relation. The redshift and mass ranges covered by our SL sample are $0.22 \le z \le 1.01$ and $5 \times10^{13} \le M_{200}/h^{-1}_{70}M_{\odot} \le 1.9\times10^{15}$, respectively. We analyze and quantify some possible effects that might bias our mass estimates, such as the presence of substructure, the region where cluster members are selected for spectroscopic follow-up, the final number of confirmed members, and line-of-sight effects. We find that 10 clusters of our sample with $N_{mem} \gtrsim 20$ show signs of dynamical substructure. However, the velocity data of only one system is inconsistent with a uni-modal distribution. We therefore assume that the substructures are only marginal and not of comparable size to the clusters themselves. Consequently, our velocity dispersion and mass estimates can be used as priors for SL mass reconstruction studies and also represent an important step toward a better understanding of the properties of the SL galaxy cluster population.

Dusty star-forming galaxies are among the most prodigious systems at high redshift (z>1), characterized by high star formation rates and huge dust reservoirs. The bright end of this population has been well characterized in recent years, but considerable uncertainties remain for fainter dusty star-forming galaxies, which are responsible for the bulk of star formation at high redshift and thus play a key role in galaxy growth and evolution. In this first paper of our series, we describe our methods for finding high redshift faint dusty galaxies using millimeter observations with ALMA. We obtained ALMA 1.1 mm mosaic images for three strong-lensing galaxy clusters from the Frontier Fields survey. The 2'x2' mosaics overlap with the deep HST WFC3/IR footprints and encompass the high magnification regions of each cluster. The combination of extremely high ALMA sensitivity and the magnification power of these clusters allows us to systematically probe the sub-mJy population of dusty star-forming galaxies over a large surveyed area. We present a description of the reduction and analysis of the ALMA continuum observations for the galaxy clusters Abell 2744 (z=0.308), MACSJ0416.1-2403 (z=0.396) and MACSJ1149.5+2223 (z=0.543), for which we reach observed rms sensitivities of 55, 59 and 71 $\mu$Jy beam$^{-1}$ respectively. We detect 12 dusty star-forming galaxies at S/N>5.0 across the three clusters, all of them presenting coincidence with NIR detected counterparts in the HST images. None of the sources fall close to the lensing caustics, thus they are not strongly lensed. The observed 1.1 mm flux densities for the total sample of galaxies range from 0.41 to 2.82 mJy, with observed effective radii spanning <0.05" to 0.37"$\pm$0.21". The lensing-corrected sizes of the detected sources appear to be in the same range as those measured in brighter samples, albeit with possibly larger dispersion.

We present V and R photometry of the gravitationally lensed quasars WFI2033-4723 and HE0047-1756. The data were taken by the MiNDSTEp collaboration with the 1.54 m Danish telescope at the ESO La Silla observatory from 2008 to 2012. Differential photometry has been carried out using the image subtraction method as implemented in the HOTPAnTS package, additionally using GALFIT for quasar photometry. The quasar WFI2033-4723 showed brightness variations of order 0.5 mag in V and R during the campaign. The two lensed components of quasar HE0047-1756 varied by 0.2-0.3 mag within five years. We provide, for the first time, an estimate of the time delay of component B with respect to A of $\Delta t= 7.6\pm1.8$ days for this object. We also find evidence for a secular evolution of the magnitude difference between components A and B in both filters, which we explain as due to a long-duration microlensing event. Finally we find that both quasars WFI2033-4723 and HE0047-1756 become bluer when brighter, which is consistent with previous studies.

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

We present an analysis of a deep (1$\sigma$=13 $\mu$Jy) cosmological 1.2-mm continuum map based on ASPECS, the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field. In the 1 arcmin$^2$ covered by ASPECS we detect nine sources at $>3.5\sigma$ significance at 1.2-mm. Our ALMA--selected sample has a median redshift of $z=1.6\pm0.4$, with only one galaxy detected at z$>$2 within the survey area. This value is significantly lower than that found in millimeter samples selected at a higher flux density cut-off and similar frequencies. Most galaxies have specific star formation rates similar to that of main sequence galaxies at the same epoch, and we find median values of stellar mass and star formation rates of $4.0\times10^{10}\ M_\odot$ and $\sim40~M_\odot$ yr$^{-1}$, respectively. Using the dust emission as a tracer for the ISM mass, we derive depletion times that are typically longer than 300 Myr, and we find molecular gas fractions ranging from $\sim$0.1 to 1.0. As noted by previous studies, these values are lower than using CO--based ISM estimates by a factor $\sim$2. The 1\,mm number counts (corrected for fidelity and completeness) are in agreement with previous studies that were typically restricted to brighter sources. With our individual detections only, we recover $55\pm4\%$ of the extragalactic background light (EBL) at 1.2 mm measured by the Planck satellite, and we recover $80\pm7\%$ of this EBL if we include the bright end of the number counts and additional detections from stacking. The stacked contribution is dominated by galaxies at $z\sim1-2$, with stellar masses of (1-3)$\times$10$^{10}$ M$_\odot$. For the first time, we are able to characterize the population of galaxies that dominate the EBL at 1.2 mm.

We present a search for [CII] line and dust continuum emission from optical dropout galaxies at $z>6$ using ASPECS, our ALMA Spectroscopic Survey in the Hubble Ultra-Deep Field (UDF). Our observations, which cover the frequency range $212-272$ GHz, encompass approximately the range $6<z<8$ for [CII] line emission and reach a limiting luminosity of L$_{\rm [CII]}\sim$(1.6-2.5)$\times$10$^{8}$ L$_{\odot}$. We identify fourteen [CII] line emitting candidates in this redshift range with significances $>$4.5 $\sigma$, two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60\% of our candidates are expected to be spurious. For one of our blindly selected [CII] line candidates, we tentatively detect the CO(6-5) line in our parallel 3-mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all [CII] candidates results in a tentative detection with $S_{1.2mm}=14\pm5\mu$Jy. This implies a dust-obscured star formation rate (SFR) of $(3\pm1)$ M$_\odot$ yr$^{-1}$. We find that the two highest--SFR objects have candidate [CII] lines with luminosities that are consistent with the low-redshift $L_{\rm [CII]}$ vs. SFR relation. The other candidates have significantly higher [CII] luminosities than expected from their UV--based SFR. At the current sensitivity it is unclear whether the majority of these sources are intrinsically bright [CII] emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for [CII] emitters at $6<z<8$ that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $\mu$Jy beam$^{-1}$ and 12.7 $\mu$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

We investigate the kinematics and ionization structure of the broad emission line region of the gravitationally lensed quasar QSO2237+0305 (the Einstein cross) using differential microlensing in the high- and low-ionization broad emission lines. We combine visible and near-infrared spectra of the four images of the lensed quasar and detect a large-amplitude microlensing effect distorting the high-ionization CIV and low-ionization H$\alpha$ line profiles in image A. While microlensing only magnifies the red wing of the Balmer line, it symmetrically magnifies the wings of the CIV emission line. Given that the same microlensing pattern magnifies both the high- and low-ionization broad emission line regions, these dissimilar distortions of the line profiles suggest that the high- and low-ionization regions are governed by different kinematics. Since this quasar is likely viewed at intermediate inclination, we argue that the differential magnification of the blue and red wings of H$\alpha$ favors a flattened, virialized, low-ionization region whereas the symmetric microlensing effect measured in CIV can be reproduced by an emission line formed in a polar wind, without the need of fine-tuned caustic configurations.

We have obtained spectropolarimetric observations of the four images of the gravitationally lensed broad absorption line quasar H1413+117. The polarization of the microlensed image D is significantly different both in the continuum and in the broad lines from the polarization of image A which is essentially unaffected by microlensing. The observations suggest that the continuum is scattered off two regions spatially separated and producing roughly perpendicular polarizations. These results are compatible with a model in which the microlensed polarized continuum comes from a compact region located in the equatorial plane close to the accretion disk and the non-microlensed continuum from an extended region located along the polar axis.

Testing the standard Shakura-Sunyaev model of accretion is a challenging task because the central region of quasars where accretion takes place is unresolved with telescopes. The analysis of microlensing in gravitationally lensed quasars is one of the few techniques that can test this model, yielding to the measurement of the size and of temperature profile of the accretion disc. We present spectroscopic observations of the gravitationally lensed broad absorption line quasar H1413+117, which reveal partial microlensing of the continuum emission that appears to originate from two separated regions: a microlensed region, corresponding to the compact accretion disc; and a non-microlensed region, more extended and contributing to at least 30\% of the total UV-continuum flux. Because this extended continuum is occulted by the broad absorption line clouds, it is not associated with the host galaxy, but rather with light scattered in the neighbourhood of the central engine. We measure the amplitude of microlensing of the compact continuum over the rest-frame wavelength range 1000-7000 \AA. Following a Bayesian scheme, we confront our measurements to microlensing simulations of an accretion disc with a temperature varying as $T \propto R^{-1/\nu}$. We find a most likely source half-light radius of $R_{1/2} = 0.61 \times 10^{16}\,$cm (i.e., 0.002\,pc) at 0.18\,$\mu$m, and a most-likely index of $\nu=0.4$. The standard disc ($\nu=4/3$) model is not ruled out by our data, and is found within the 95\% confidence interval associated with our measurements. We demonstrate that, for H1413+117, the existence of an extended continuum in addition to the disc emission only has a small impact on the inferred disc parameters, and is unlikely to solve the tension between the microlensing source size and standard disc sizes, as previously reported in the literature.

We report the discovery of a large timescale candidate microlensing event of a bulge stellar source based on near-infrared observations with the VISTA Variables in the Via Lactea Survey (VVV). The new microlensing event is projected only 3.5 arcmin away from the center of the globular cluster NGC 6553. The source appears to be a bulge giant star with magnitude Ks = 13.52, based on the position in the color-magnitude diagram. The foreground lens may be located in the globular cluster, which has well-known parameters such as distance and proper motions. If the lens is a cluster member, we can directly estimate its mass simply following Paczynski et al. (1996) which is a modified version of the more general case due to Refsdal. In that case, the lens would be a massive stellar remnant, with M = 1.5-3.5 Msun. If the blending fraction of the microlensing event appears to be small, and this lens would represent a good isolated black hole (BH) candidate, that would be the oldest BH known. Alternative explanations (with a larger blending fraction) also point to a massive stellar remnant if the lens is located in the Galactic disk and does not belong to the globular cluster.

Aims. Construction of a new quasar candidate catalog from the Red-Sequence Cluster Survey 2 (RCS-2), identified solely from photometric information using an automated algorithm suitable for large surveys. The algorithm performance is tested using a well-defined SDSS spectroscopic sample of quasars and stars. Methods. The Random Forest algorithm constructs the catalog from RCS-2 point sources using SDSS spectroscopically-confirmed stars and quasars. The algorithm identifies putative quasars from broadband magnitudes (g, r, i, z) and colours. Exploiting NUV GALEX measurements for a subset of the objects, we refine the classifier by adding new information. An additional subset of the data with WISE W1 and W2 bands is also studied. Results. Upon analyzing 542,897 RCS-2 point sources, the algorithm identified 21,501 quasar candidates, with a training-set-derived precision (the fraction of true positives within the group assigned quasar status) of 89.5% and recall (the fraction of true positives relative to all sources that actually are quasars) of 88.4%. These performance metrics improve for the GALEX subset; 6,530 quasar candidates are identified from 16,898 sources, with a precision and recall respectively of 97.0% and 97.5%. Algorithm performance is further improved when WISE data are included, with precision and recall increasing to 99.3% and 99.1% respectively for 21,834 quasar candidates from 242,902 sources. We compile our final catalog (38,257) by merging these samples and removing duplicates. An observational follow up of 17 bright (r < 19) candidates with long-slit spectroscopy at DuPont telescope (LCO) yields 14 confirmed quasars. Conclusions. The results signal encouraging progress in the classification of point sources with Random Forest algorithms to search for quasars within current and future large-area photometric surveys.

Using infrared spectra of the z = 1.693 quadruply lensed quasar HE0435-1223 acquired in 2009 with the spectrograph SINFONI at the ESO Very Large Telescope, we have detected a clear microlensing effect in images A and D. While microlensing affects the blue and red wings of the H{\alpha} line profile in image D very differently, it de-magnifies the line core in image A. The combination of these different effects sets constraints on the line-emitting region; these constraints suggest that a rotating ring is at the origin of the H{\alpha} line. Visible spectra obtained in 2004 and 2012 indicate that the MgII line profile is microlensed in the same way as the H{\alpha} line. Our results therefore favour flattened geometries for the low-ionization line-emitting region, for example, a Keplerian disk. Biconical models cannot be ruled out but require more fine-tuning. Flux ratios between the different images are also derived and confirm flux anomalies with respect to estimates from lens models with smooth mass distributions.

Multiply-imaged quasars and AGNs observed in the mid-infrared (MIR) range are commonly assumed to be unaffected by the microlensing produced by the stars in their lensing galaxy. In this paper, we investigate the validity domain of this assumption. Indeed, that premise disregards microlensing of the accretion disc in the MIR range, and does not account for recent progress in our knowledge of the dusty torus. To simulate microlensing, we first built a simplified image of the quasar composed of an accretion disc, and of a larger ring-like torus. The mock quasars are then microlensed using an inverse ray-shooting code. We simulated the wavelength and size dependence of microlensing for different lensed image types and fraction of compact objects projected in the lens. This allows us to derive magnification probabilities as a function of wavelength, as well as to calculate the microlensing-induced deformation of the spectral energy distribution of the lensed images. We find that microlensing variations as large as 0.1 mag are very common at 11 microns (observer-frame). The main signal comes from microlensing of the accretion disc, which may be significant even when the fraction of flux from the disc is as small as 5 % of the total flux. We also show that the torus of sources with Lbol <~ 10^45 erg/s is expected to be noticeably microlensed. Microlensing may thus be used to get insight into the rest near-infrared inner structure of AGNs. Finally, we investigate whether microlensing in the mid-infrared can alter the so-called Rcusp relation that links the fluxes of the lensed images triplet produced when the source lies close to a cusp macro-caustic. This relation is commonly used to identify massive (dark-matter) substructures in lensing galaxies. We find that significant deviations from Rcusp may be expected, which means that microlensing can explain part of the flux ratio problem.

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.

With the aim of characterizing the flux and color variations of the multiple components of the gravitationally lensed quasar UM673 as a function of time, we have performed multi-epoch and multi-band photometric observations with the Danish 1.54m telescope at the La Silla Observatory. The observations were carried out in the VRi spectral bands during four seasons (2008--2011). We reduced the data using the PSF (Point Spread Function) photometric technique as well as aperture photometry. Our results show for the brightest lensed component some significant decrease in flux between the first two seasons (+0.09/+0.11/+0.05 mag) and a subsequent increase during the following ones (-0.11/-0.11/-0.10 mag) in the V/R/i spectral bands, respectively. Comparing our results with previous studies, we find smaller color variations between these seasons as compared with previous ones. We also separate the contribution of the lensing galaxy from that of the fainter and close lensed component.

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dln\Sigma/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

The microlensing event OGLE-2008-BLG-510 is characterised by an evident asymmetric shape of the peak, promptly detected by the ARTEMiS system in real time. The skewness of the light curve appears to be compatible both with binary-lens and binary-source models, including the possibility that the lens system consists of an M dwarf orbited by a brown dwarf. The detection of this microlensing anomaly and our analysis demonstrates that: 1) automated real-time detection of weak microlensing anomalies with immediate feedback is feasible, efficient, and sensitive, 2) rather common weak features intrinsically come with ambiguities that are not easily resolved from photometric light curves, 3) a modelling approach that finds all features of parameter space rather than just the `favourite model' is required, and 4) the data quality is most crucial, where systematics can be confused with real features, in particular small higher-order effects such as orbital motion signatures. It moreover becomes apparent that events with weak signatures are a silver mine for statistical studies, although not easy to exploit. Clues about the apparent paucity of both brown-dwarf companions and binary-source microlensing events might hide here.

We present the first galaxy scale lens catalog from the second Red-Sequence Cluster Survey. The catalog contains 60 lensing system candidates comprised of Luminous Red Galaxy (LRG) lenses at 0.2 < z < 0.5 surrounded by blue arcs or apparent multiple images of background sources. The catalog is a valuable complement to previous galaxy-galaxy lens catalogs as it samples an intermediate lens redshift range and is composed of bright sources and lenses that allow easy follow-up for detailed analysis. Mass and mass-to-light ratio estimates reveal that the lens galaxies are massive (<M>~5.5x10e11 M_sun/h) and rich in dark matter (<M/L>~14 M_sun/L_sun,B*h). Even though a slight increasing trend in the mass-to-light ratio is observed from z=0.2 to z=0.5, current redshift and light profile measurements do not allow stringent constraints on the mass-to-light ratio evolution of LRGs.

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with {\Lambda}CDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent {\Lambda}CDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

We aim to use microlensing taking place in the lensed quasar Q2237+0305 to study the structure of the broad line region and measure the size of the region emitting the CIV and CIII] lines. Methods: Based on 39 spectrophotometric monitoring data points obtained between Oct. 2004 and Dec. 2007, we derived lightcurves for the CIV and CIII] emission lines. We used three different techniques to analyse the microlensing signal. Different components of the lines (narrow, broad and very broad) are identified and studied. We built a library of simulated microlensing lightcurves that reproduce the signal observed in the continuum and in the lines provided only the source size is changed. A Bayesian analysis scheme is then developed to derive the size of the various components of the BLR. Results: 1. The half-light radius of the region emitting the CIV line is found to be R_CIV ~ 66^{+110}_{-46} lt-days = 0.06$^{+0.09}_{-0.04}$ pc = 1.7$^{+2.8}_{-1.1}$\,10$^{17}$ cm (at 68.3% CI). Similar values are obtained for CIII]. Relative sizes of the carbon-line and V-band continuum emitting-regions are also derived with median values of R(line)/R(cont) in the range 4 to 29, depending of the FWHM of the line component. 2. The size of the CIV emitting region agrees with the Radius-Luminosity relationship derived from reverberation mapping. Using the virial theorem we derive the mass of the black hole in Q2237+0305 to be M_BH ~ 10^{8.3+/-0.3} M_sun. 3. We find that the CIV and CIII] lines are produced in at least 2 spatially distinct regions, the most compact one giving rise to the broadest component of the line. The broad and narrow line profiles are slightly different for CIV and CIII]. 4. Our analysis suggests a different structure for the CIV and FeII+III emitting regions, with the latter produced in the inner part of the BLR or in a less extended emitting region than CIV.

aims: We present VRi photometric observations of the quadruply imaged quasar HE 0435-1223, carried out with the Danish 1.54m telescope at the La Silla Observatory. Our aim was to monitor and study the magnitudes and colors of each lensed component as a function of time. methods: We monitored the object during two seasons (2008 and 2009) in the VRi spectral bands, and reduced the data with two independent techniques: difference imaging and PSF (Point Spread Function) fitting.results: Between these two seasons, our results show an evident decrease in flux by ~0.2-0.4 magnitudes of the four lensed components in the three filters. We also found a significant increase (~0.05-0.015) in their V-R and R-i color indices. conclusions: These flux and color variations are very likely caused by intrinsic variations of the quasar between the observed epochs. Microlensing effects probably also affect the brightest "A" lensed component.

Among the 100 strong lens candidates found in the COSMOS field, 20 with redshifts in the range [0.34,1.13], feature multiple images of background sources. Using the multi-wavelength coverage of the field and its spectroscopic follow-up, we characterize the evolution with redshift of the environment and of the dark-matter (DM) fraction of the lens galaxies. We present new redshift of the strong lens candidates. The lens environment is characterized by the projected 10 closest galaxies around each lens and by the number of galaxies with a projected distance less than 1Mpc at the lens galaxy redshift. In both cases, we perform similar measurements on a control sample of twin non-lens early type galaxies (ETGs). In addition, we identify group members and field galaxies in the X-ray and optical catalogs of galaxy groups. From those catalogs, we measure the external shear contribution at the lens galaxy positions. The systems are then modeled using a SIE plus the external shear due to the groups. We observe that the average stellar mass of lens galaxies increases with z and that the environment of lens galaxies is compatible with that of the twins. During the lens modeling, we notice that, when let free, the external shear points in a direction which is the mean direction of the external shear due to groups and of the closest galaxy to the lens. We notice that the DM fraction of the lens galaxies within the Einstein radius decreases as the redshift increases. Given these, we conclude that, while the environment of lens galaxies is compatible with that of non-lens ETGS, their mass properties evolves significantly with redshift: it is still not clear whether this advocates in favor of a stronger lensing bias toward massive objects at high redshift or is simply representative of the high proportion of massive and high stellar density galaxies at high redshift.

On the basis of 16 years of spectroscopic observations of the four components of the gravitationally lensed broad absorption line (BAL) quasar H1413+117, covering the ultraviolet to visible rest-frame spectral range, we analyze the spectral differences observed in the P Cygni-type line profiles and have used the microlensing effect to derive new clues to the BAL profile formation. We confirm that the spectral differences observed in component D can be attributed to a microlensing effect lasting at least a decade. We show that microlensing magnifies the continuum source in image D, leaving the emission line region essentially unaffected. We interpret the differences seen in the absorption profiles of component D as the result of an emission line superimposed onto a nearly black absorption profile. We also find that the continuum source and a part of the broad emission line region are likely de-magnified in component C, while components A and B are not affected by microlensing. We show that microlensing of the continuum source in component D has a chromatic dependence compatible with the thermal continuum emission of a standard Shakura-Sunyaev accretion disk. Using a simple decomposition method to separate the part of the line profiles affected by microlensing and coming from a compact region from the part unaffected by this effect and coming from a larger region, we disentangle the true absorption line profiles from the true emission line profiles. The extracted emission line profiles appear double-peaked, suggesting that the emission is occulted by a strong absorber, narrower in velocity than the full absorption profile, and emitting little by itself. We propose that the outflow around H1413+117 is constituted by a high-velocity polar flow and a denser, lower velocity disk seen nearly edge-on.

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

Gravitational microlensing is not only a successful tool for discovering distant exoplanets, but it also enables characterization of the lens and source stars involved in the lensing event. In high magnification events, the lens caustic may cross over the source disk, which allows a determination of the angular size of the source and additionally a measurement of its limb darkening. When such extended-source effects appear close to maximum magnification, the resulting light curve differs from the characteristic Paczynski point-source curve. The exact shape of the light curve close to the peak depends on the limb darkening of the source. Dense photometric coverage permits measurement of the respective limb-darkening coefficients. In the case of microlensing event OGLE 2008-BLG-290, the K giant source star reached a peak magnification of about 100. Thirteen different telescopes have covered this event in eight different photometric bands. Subsequent light-curve analysis yielded measurements of linear limb-darkening coefficients of the source in six photometric bands. The best-measured coefficients lead to an estimate of the source effective temperature of about 4700 +100-200 K. However, the photometric estimate from colour-magnitude diagrams favours a cooler temperature of 4200 +-100 K. As the limb-darkening measurements, at least in the CTIO/SMARTS2 V and I bands, are among the most accurate obtained, the above disagreement needs to be understood. A solution is proposed, which may apply to previous events where such a discrepancy also appeared.

[abridged] We present VLT/FORS1 spectroscopic follow-up observations and HST/WFPC2 imaging of the system COSMOS 5921+0638, which exhibits quadruply lensed images and a perfect Einstein ring. We investigate the nature of COSMOS 5921+0638 by studying its photometric, spectroscopic and physical properties. We completed both an environmental analysis and detailed analytical and grid-based mass modeling of the system. We measured the redshifts of the lensing galaxy in COSMOS 5921+0638 (z_l=0.551+/-0.001) and 9 additional galaxies in the field (5 of them at z~0.35). The redshift of the lensed source was inferred by identifying a candidate Lya line at z_s=3.14+/-0.05. The mass modeling reveals the requirement of a small external shear (gamma=0.038), which is slightly larger than the lensing contribution expected by galaxy groups along the line-of-sight (kappa_groups~0.01 and gamma_groups~0.005). The estimated time-delays between the different images are of the order of hours to half a week and the total magnification of the background source is mu~150. The measured mass-to-light ratio of the lensing galaxy within the Einstein ring is M/L_B~8.5+/-1.6. Our analysis indicates that the ring and point-like structures in COSMOS 5921+0638 consist of a lensed high redshift galaxy hosting a low luminosity AGN (LLAGN). Flux ratio anomalies observed in the lensed AGN images are probably due to microlensing by stars in the lensing galaxy and/or a combination of static phenomena. Because of its short time-delays and the possibility of microlensing, COSMOS 5921+0638 is a promising laboratory for future studies of LLAGNs.

We present and analyse light curves of four transits of the Southern hemisphere extrasolar planetary system WASP-4, obtained with a telescope defocussed so the radius of each point spread function was 17 arcsec (44 pixels). This approach minimises both random and systematic errors, allowing us to achieve scatters of between 0.60 and 0.88 mmag per observation over complete transit events. The light curves are augmented by published observations and analysed using the JKTEBOP code. The results of this process are combined with theoretical stellar model predictions to derive the physical properties of the WASP-4 system. We find that the mass and radius of the planet are M_b = 1.289 {+0.090 -0.090} {+0.039 -0.000} MJup and R_b = 1.371 {+0.032 -0.035} {+0.021 -0.000} RJup, respectively (statistical and systematic uncertainties). These quantities give a surface gravity and density of g_b = 17.03 +0.97 -0.54 m/s2 and rho_b = 0.500 {+0.032 -0.021} {+0.000 -0.008} rhoJup, and fit the trends for short-period extrasolar planets to have relatively high masses and surface gravities. WASP-4 is now one of the best-quantified transiting extrasolar planetary systems, and significant further progress requires improvements to our understanding of the physical properties of low-mass stars.

We present high-precision photometry of two transit events of the extrasolar planetary system WASP-5, obtained with the Danish 1.54m telescope at ESO La Silla. In order to minimise both random and flat-fielding errors, we defocussed the telescope so its point spread function approximated an annulus of diameter 40 pixels (16 arcsec). Data reduction was undertaken using standard aperture photometry plus an algorithm for optimally combining the ensemble of comparison stars. The resulting light curves have point-to-point scatters of 0.50 mmag for the first transit and 0.59 mmag for the second. We construct detailed signal to noise calculations for defocussed photometry, and apply them to our observations. We model the light curves with the JKTEBOP code and combine the results with tabulated predictions from theoretical stellar evolutionary models to derive the physical properties of the WASP-5 system. We find that the planet has a mass of M_b = 1.637 +/- 0.075 +/- 0.033 Mjup, a radius of R_b = 1.171 +/- 0.056 +/- 0.012 Rjup, a large surface gravity of g_b = 29.6 +/- 2.8 m/s2 and a density of rho_b = 1.02 +/- 0.14 +/- 0.01 rhojup (statistical and systematic uncertainties). The planet's high equilibrium temperature of T_eq = 1732 +/- 80 K makes it a good candidate for detecting secondary eclipses.

We outline a method for fitting binary-lens caustic-crossing microlensing events based on the alternative model parameterisation proposed and detailed in Cassan (2008). As an illustration of our methodology, we present an analysis of OGLE-2007-BLG-472, a double-peaked Galactic microlensing event with a source crossing the whole caustic structure in less than three days. In order to identify all possible models we conduct an extensive search of the parameter space, followed by a refinement of the parameters with a Markov Chain-Monte Carlo algorithm. We find a number of low-chi2 regions in the parameter space, which lead to several distinct competitive best models. We examine the parameters for each of them, and estimate their physical properties. We find that our fitting strategy locates several minima that are difficult to find with other modelling strategies and is therefore a more appropriate method to fit this type of events.

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 present the continuation of our long-term spectroscopic monitoring of the gravitationally lensed quasar QSO 2237+0305. We investigate the chromatic variations observed in the UV/optical continuum of both quasar images A and B, and compare them with numerical simulations to infer the energy profile of the quasar accretion disk. Our procedure combines the microlensing ray-shooting technique with Bayesian analysis, and derives probability distributions for the source sizes as a function of wavelength. We find that the effective caustic crossing timescale is 4.0+/-1.0 months. Using a robust prior on the effective transverse velocity, we find that the source responsible for the UV/optical continuum has an energy profile well reproduced by a power-law R lambda^{zeta} with zeta=1.2+/-0.3, where R is the source size responsible for the emission at wavelength lambda. This is the first accurate, model-independent determination of the energy profile of a quasar accretion disk on such small scales.

CONTEXT: Gravitationally lensed quasars constitute an independent tool to derive H0 through time-delays; they offer as well the opportunity to study the mass distribution and interstellar medium of their lensing galaxies and, through microlensing they also allow one to study details of the emitting source. AIMS: For such studies, one needs to have an excellent knowledge of the close environment of the lensed images in order to model the lensing potential: this means observational data over a large field-of-view and spectroscopy at high spatial resolution. METHODS: We present VIMOS integral field observations around four lensed quasars: HE 0230-2130, RX J0911.4+0551, H 1413+117 and B 1359+154. Using the low, medium and high resolution modes, we study the quasar images and the quasar environments, as well as provide a detailed report of the data reduction. RESULTS: Comparison between the quasar spectra of the different images reveals differences for HE 0230-2130, RX J0911.4+0551 and H 1413+117: flux ratios between the images of the same quasar are different when measured in the emission lines and in the continuum. We have also measured the redshifts of galaxies in the neighborhood of HE 0230-2130 and RX J0911.4+0551 which possibly contribute to the total lensing potential. CONCLUSIONS: A careful analysis reveals that microlensing is the most natural explanation for the (de)magnification of the continuum emitting region of the background sources. In HE 0230-2130, image D is likely to be affected by microlensing magnification; in RX J0911.4+0551, images A1 and A3 are likely to be modified by microlensing de-magnification and in H 1413+117, at least image D is affected by microlensing.

We use the high magnification event seen in the 1999 OGLE campaign light curve of image C of the quadruply imaged gravitational lens Q2237+0305 to study the structure of the quasar engine. We have obtained g'- and r'-band photometry at the Apache Point Observatory 3.5m telescope where we find that the event has a smaller amplitude in the r'-band than in the g'- and OGLE V-bands. By comparing the light curves with microlensing simulations we obtain constraints on the sizes of the quasar regions contributing to the g'- and r'-band flux. Assuming that most of the surface mass density in the central kiloparsec of the lensing galaxy is due to stars and by modeling the source with a Gaussian profile, we obtain for the Gaussian width 1.20 x 10^15 sqrt(M/0.1M_sun)cm < sigma_g' < 7.96 x 10^15 sqrt(M/0.1Msun) cm, where M is the mean microlensing mass, and a ratio sigma_r'/sigma_g'=1.25^{+0.45}_{-0.15}. With the limits on the velocity of the lensing galaxy from Gil-Merino et al. (2005) as our only prior, we obtain 0.60 x 10^15 sqrt(M/0.1Msun) cm < sigma_g' < 1.57 x 10^15 sqrt(M/0.1Msun) cm and a ratio sigma_r'/sigma_g'=1.45^{+0.90}_{-0.25} (all values at 68 percent confidence). Additionally, from our microlensing simulations we find that, during the chromatic microlensing event observed, the continuum emitting region of the quasar crossed a caustic at >72 percent confidence.