Search

We report the discovery of the first example of an Einstein zig-zag lens, an extremely rare lensing configuration. In this system, J1721+8842, six images of the same background quasar are formed by two intervening galaxies, one at redshift $z_1 = 0.184$ and a second one at $z_2 = 1.885$. Two out of the six multiple images are deflected in opposite directions as they pass the first lens galaxy on one side, and the second on the other side -- the optical paths forming zig-zags between the two deflectors. In this letter, we demonstrate that J1721+8842, previously thought to be a lensed dual quasar, is in fact a compound lens with the more distant lens galaxy also being distorted as an arc by the foreground galaxy. Evidence supporting this unusual lensing scenario includes: 1- identical light curves in all six lensed quasar images obtained from two years of monitoring at the Nordic Optical Telescope; 2- detection of the additional deflector at redshift $z_2 = 1.885$ in JWST/NIRSpec IFU data; and 3- a multiple-plane lens model reproducing the observed image positions. This unique configuration offers the opportunity to combine two major lensing cosmological probes: time-delay cosmography and dual source-plane lensing since J1721+8842 features multiple lensed sources forming two distinct Einstein radii of different sizes, one of which being a variable quasar. We expect tight constraints on the Hubble constant and the equation of state of dark energy by combining these two probes on the same system. The $z_2 = 1.885$ deflector, a quiescent galaxy, is also the highest-redshift strong galaxy-scale lens with a spectroscopic redshift measurement., Comment: 6 pages, 3 figures

This is the second in a series of papers in which we use JWST MIRI multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of 9 systems. The warm dust region is compact and sensitive to perturbations by populations of halos down to masses $\sim 10^6$ M$_{\odot}$. Using these warm dust flux-ratio measurements in combination with 5 previous narrow-line flux-ratio measurements, we constrain the halo mass function. In our model, we allow for complex deflector macromodels with flexible third and fourth-order multipole deviations from ellipticity, and we introduce an improved model of the tidal evolution of subhalos. We constrain a WDM model and find an upper limit on the half-mode mass of $10^{7.6} M_\odot$ at posterior odds of 10:1. This corresponds to a lower limit on a thermally produced dark matter particle mass of 6.1 keV. This is the strongest gravitational lensing constraint to date, and comparable to those from independent probes such as the Ly$\alpha$ forest and Milky Way satellite galaxies.

The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light is emitted from a physical region of the quasar with an angular scale of milli-arcseconds or larger; this minimizes microlensing effects by stars within the deflector. The warm dust region of quasars fits this criterion, as it has parsec-size physical scales and dominates the spectral energy distribution of quasars at wavelengths greater than 10$\mu$m. The JWST Mid-Infrared Instrument (MIRI) is adept at detecting redshifted light in this wavelength range, offering both the spatial resolution and sensitivity required for accurate gravitational lensing flux ratio measurements. Here, we introduce our survey designed to measure the warm dust flux ratios of 31 lensed quasars. We discuss the flux-ratio measurement technique and present results for the first target, DES J0405-3308. We find that we can measure the quasar warm dust flux ratios with 3% precision. Our simulations suggest that this precision makes it feasible to detect the presence of 10$^7$ M$_\odot$ dark matter halos at cosmological distances. Such halos are expected to be completely dark in Cold Dark Matter models.

We present new broadband X-ray observations of the $z \sim 2.5$ lensed quasar 2MASS J1042+1641, combining $XMM$-$Newton$, $Chandra$ and $NuSTAR$ to provide coverage of the X-ray spectrum over the 0.3$-$40 keV bandpass in the observed frame, corresponding to the $\sim$1$-$140 keV band in the rest-frame of 2MASS J1042+1641. The X-ray data show clear evidence for strong (but still Compton-thin) X-ray absorption, $N_{\rm{H}} \sim 3-4 \times 10^{23}$ cm$^{-2}$, in addition to significant reprocessing by Compton-thick material that must lie away from our line-of-sight to the central X-ray source. We test two different interpretations for the latter: first that the reprocessing occurs in a classic AGN torus, as invoked in unification models, and second that the reprocessing occurs in the accretion disc. Both models can successfully reproduce the observed spectra, and both imply that the source is viewed at moderately low inclinations ($i < 50^{\circ}$) despite the heavy line-of-sight absorption. Combining the X-ray data with infrared data from $WISE$, the results seen from 2MASS J1042+1641 further support the recent suggestion that large X-ray and IR surveys may together be able to identify good lensed quasar candidates in advance of detailed imaging studies., Comment: 13 pages, 7 figures, accepted for publication in MNRAS

We report periodic oscillations in the 15-year long optical light curve of the gravitationally lensed quasar QJ0158-4325. The signal is enhanced during a high magnification microlensing event undergone by the fainter lensed image of the quasar, between 2003 and 2010. We measure a period of $P_{o}=172.6\pm0.9$ days. We explore four scenarios to explain the origin of the periodicity: 1- the high magnification microlensing event is due to a binary star in the lensing galaxy, 2- QJ0158-4325 contains a massive binary black hole system in its final dynamical stage before merging, 3- the quasar accretion disk contains a bright inhomogeneity in Keplerian motion around the black hole, and 4- the accretion disk is in precession. Among these four scenarios, only a supermassive binary black hole can account for both the short observed period and the amplitude of the signal, through the oscillation of the accretion disk towards and away from high-magnification regions of a microlensing caustic. The short measured period implies that the semi-long axis of the orbit is $\sim10^{-3}$pc, and the coalescence timescale is $t_{coal}\sim1000$ years, assuming that the decay of the orbit is solely powered by the emission of gravitational waves. The probability of observing a system so close to coalescence suggests either a much larger population of supermassive black hole binaries than predicted, or, more likely, that some other mechanism significantly increases the coalescence timescale. Three tests of the binary black hole hypothesis include: i) the recurrence of oscillations in photometric monitoring during any future microlensing events in either image, ii) spectroscopic detection of Doppler shifts (up to 0.01$c$), and iii) the detection of gravitational waves through Pulsar Timing Array experiments, such as the SKA, which will have the sensitivity to detect the $\sim$100 nano-hertz emission., Comment: 18 pages, 15 figures, published in A&A

We report the spectroscopic follow-up of 175 lensed quasar candidates selected using Gaia Data Release 2 observations following Lemon et al. 2019. Systems include 86 confirmed lensed quasars and a further 17 likely lensed quasars based on imaging and/or similar spectra. We also confirm 11 projected quasar pairs and 11 physical quasar pairs, while 25 systems are left as unclassified quasar pairs -- pairs of quasars at the same redshift, which could be either distinct quasars or potential lensed quasars. Especially interesting objects include 8 quadruply imaged quasars of which two have BAL sources, an apparent triple, and a doubly lensed LoBaL quasar. The source redshifts and image separations of these new lenses range between 0.65 - 3.59 and 0.78 - 6.23 arcseconds respectively. We compare the known population of lensed quasars to an updated mock catalogue at image separations between 1 and 4 arcseconds, showing a very good match at z<1.5. At z>1.5, only 47% of the predicted number are known, with 56% of these missing lenses at image separations below 1.5 arcseconds. The missing higher-redshift, small-separation systems will have fainter lensing galaxies, and are partially explained by the unclassified quasar pairs and likely lenses presented in this work, which require deeper imaging. Of the 11 new reported projected quasar pairs, 5 have impact parameters below 10 kpc, almost tripling the number of such systems, which can probe the innermost regions of quasar host galaxies through absorption studies. We also report four new lensed galaxies discovered through our searches, with source redshifts ranging from 0.62 to 2.79., Comment: 30 pages, 19 figures, submitted to MNRAS

Gravitational time delays provide a powerful one step measurement of $H_0$, independent of all other probes. One key ingredient in time delay cosmography are high accuracy lens models. Those are currently expensive to obtain, both, in terms of computing and investigator time (10$^{5-6}$ CPU hours and $\sim$ 0.5-1 year, respectively). Major improvements in modeling speed are therefore necessary to exploit the large number of lenses that are forecast to be discovered over the current decade. In order to bypass this roadblock, building on the work by Shajib et al. (2019), we develop an automated modeling pipeline and apply it to a sample of 30 quadruply imaged quasars and one lensed compact galaxy, observed by the Hubble Space Telescope in multiple bands. Our automated pipeline can derive models for 30/31 lenses with few hours of human time and <100 CPU hours of computing time for a typical system. For each lens, we provide measurements of key parameters and predictions of magnification as well as time delays for the multiple images. We characterize the cosmography-readiness of our models using the stability of differences in Fermat potential (proportional to time delay) w.r.t. modeling choices. We find that for 10/30 lenses our models are cosmography or nearly cosmography grade (<3% and 3-5% variations). For 6/30 lenses the models are close to cosmography grade (5-10%). These results are based on informative priors and will need to be confirmed by further analysis. However, they are also likely to improve by extending the pipeline modeling sequence and options. In conclusion, we show that uniform cosmography grade modeling of large strong lens samples is within reach., Comment: 40 pages, 24 figures, 11 tables. Submitted to MNRAS

We present a search for galaxy-scale strong gravitational lenses in the initial 2 500 square degrees of the Canada-France Imaging Survey (CFIS). We designed a convolutional neural network (CNN) committee that we applied to a selection of 2 344 002 exquisite-seeing $r$-band images of color-selected luminous red galaxies (LRGs). Our classification uses a realistic training set where the lensing galaxies and the lensed sources are both taken from real data, namely the CFIS $r$-band images themselves and the Hubble Space Telescope (HST). A total of 9 460 candidates obtain a score above 0.5 with the CNN committee. After a visual inspection of the candidates, we find a total of 133 lens candidates, of which 104 are completely new. The set of false positives mainly contains ring, spiral, and merger galaxies, and to a lesser extent galaxies with nearby companions. We classify 32 of the lens candidates as secure lenses and 101 as maybe lenses. For the 32 highest quality lenses, we also fit a singular isothermal ellipsoid mass profile with external shear along with an elliptical Sersic profile for the lens and source light. This automated modeling step provides distributions of properties for both sources and lenses that have Einstein radii in the range $0.5\arcsec<\theta_E<2.5\arcsec$. Finally, we introduce a new lens and/or source single-band deblending algorithm based on auto-encoder representation of our candidates. This is the first time an end-to-end lens-finding and modeling pipeline is assembled together, in view of future lens searches in a single band, as will be possible with Euclid., Comment: 29 pages, 21 figures, accepted by A&A, in press

We report the discovery of five new doubly-imaged lensed quasars from the first 2500 square degrees of the ongoing Canada-France Imaging Survey (CFIS), which is a component of the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), selected from initial catalogues of either Gaia pairs or MILLIQUAS quasars. We take advantage of the deep, 0.6'' median-seeing $r$-band imaging of CFIS to confirm the presence of multiple point sources with similar colour of $u-r$, via convolution of the Laplacian of the point spread function. Requiring similar-colour point sources with flux ratios less than 2.5 mag in $r$-band, reduces the number of candidates from 256314 to 7815. After visual inspection we obtain 30 high-grade candidates, and prioritise spectroscopic follow-up for those showing signs of a lensing galaxy upon subtraction of the point sources. We obtain long-slit spectra for 18 candidates with ALFOSC on the 2.56-m Nordic Optical Telescope (NOT), confirming five new doubly lensed quasars with $1.21

Strongly lensed quadruply imaged quasars (quads) are extraordinary objects. They are very rare in the sky -- only a few tens are known to date -- and yet they provide unique information about a wide range of topics, including the expansion history and the composition of the Universe, the distribution of stars and dark matter in galaxies, the host galaxies of quasars, and the stellar initial mass function. Finding them in astronomical images is a classic "needle in a haystack" problem, as they are outnumbered by other (contaminant) sources by many orders of magnitude. To solve this problem, we develop state-of-the-art deep learning methods and train them on realistic simulated quads based on real images of galaxies taken from the Dark Energy Survey, with realistic source and deflector models, including the chromatic effects of microlensing. The performance of the best methods on a mixture of simulated and real objects is excellent, yielding area under the receiver operating curve in the range 0.86 to 0.89. Recall is close to 100% down to total magnitude i~21 indicating high completeness, while precision declines from 85% to 70% in the range i~17-21. The methods are extremely fast: training on 2 million samples takes 20 hours on a GPU machine, and 10^8 multi-band cutouts can be evaluated per GPU-hour. The speed and performance of the method pave the way to apply it to large samples of astronomical sources, bypassing the need for photometric pre-selection that is likely to be a major cause of incompleteness in current samples of known quads., Comment: 17 pages, 14 figures, submitted to MNRAS

High-redshift binary quasars provide key insights into mergers and quasar activity, and are useful tools for probing the spatial kinematics and chemistry of galaxies along the line-of-sight. However, only three sub-10-kpc binaries have been confirmed above $z=1$. Gravitational lensing would provide a way to easily resolve such binaries, study them in higher resolution, and provide more sightlines, though the required alignment with a massive foreground galaxy is rare. Through image deconvolution of StanCam Nordic Optical Telescope (NOT) monitoring data, we reveal two further point sources in the known, $z \approx 2.38$, quadruply lensed quasar (quad), J1721+8842. An ALFOSC/NOT long-slit spectrum shows that the brighter of these two sources is a quasar with $z = 2.369 \pm 0.007$ based on the C III] line, while the C III] redshift of the quad is $z = 2.364 \pm 0.003$. Lens modelling using point source positions rules out a single source model, favouring an isothermal lens mass profile with two quasar sources separated by $\sim6.0$ kpc (0.73$^{\prime \prime}$) in projection. Given the resolving ability from lensing and current lensed quasar statistics, this discovery suggests a large population of undiscovered, unlensed sub-10-kpc binaries. We also analyse spectra of two images of the quad, showing narrow Ly$\alpha$ emission within the trough of a proximate damped Ly$\alpha$ absorber (PDLA). An apparent mismatch between the continuum and narrow line flux ratios provides a new potential tool for simultaneously studying microlensing and the quasar host galaxy. Signs of the PDLA are also seen in the second source, however a deeper spectrum is still required to confirm this. Thanks to the multiple lines-of-sight from lensing and two quasar sources, this system offers simultaneous sub-parsec and kpc-scale probes of a PDLA., Comment: 10 pages, 7 figures, submitted to A&A

We performed a search for strong lens galaxy-scale systems in the first data release of the Dark Energy Survey (DES), from a color-selected parent sample of 18~745~029 Luminous Red Galaxies (LRGs). Our search was based on a Convolutional Neural Network (CNN) to grade our LRG selection with values between 0 (non-lens) and 1 (lens). Our training set was data-driven, i.e. using lensed sources taken from HST COSMOS images and where the light distribution of the lens plane was taken directly from DES images of our LRGs. A total of 76~582 cutouts obtained a score above 0.9. These were visually inspected and resulted in two catalogs. The first one contains 405 lens candidates, where 90 present clear lensing features and counterparts, while the others 315 require more evidence, such as higher resolution images or spectra to be conclusive. A total of 186 candidates were totally new identified in this search. The second catalog includes 539 ring galaxy candidates that will be useful to train CNNs against this type of false positives. For the 90 best lens candidates we carried out color-based deblending of the lens and source light without fitting any analytical profile to the data. The method turned out to be very efficient in the deblending, even for very compact objects and for objects with very complex morphology. Finally, from the 90 best lens candidates we selected 52 systems having one single deflector, to test an automated modeling pipeline which successfully modeled 79\% of the sample within an acceptable amount of computing time., Comment: 37 pages, 31 figures, submitted to A&A

We present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the Max Planck Institute for Astrophysics (MPIA) 2.2 m 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 measured time delays between all pairs of multiple images with only one or two seasons of monitoring with the exception of the time delays relative to image D of PSJ 1606-2333. The most precise estimate was obtained for the delay between image A and image B of DES 0407-5006, where $\tau_{AB} = -128.4^{+3.5}_{-3.8}$ d (2.8% precision) including systematics due to extrinsic variability in the light curves. For HE 0047-1756, we combined our high-cadence data with measurements from decade-long light curves from previous COSMOGRAIL campaigns, and reach a precision of 0.9 d 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., Comment: 14 pages, 9 figures, 4 Tables, published in A&A

We present a systematic search for wide-separation (Einstein radius >1.5"), galaxy-scale strong lenses in the 30 000 sq.deg of the Pan-STARRS 3pi survey on the Northern sky. With long time delays of a few days to weeks, such systems are particularly well suited for catching strongly lensed supernovae with spatially-resolved multiple images and open new perspectives on early-phase supernova spectroscopy and cosmography. We produce a set of realistic simulations by painting lensed COSMOS sources on Pan-STARRS image cutouts of lens luminous red galaxies with known redshift and velocity dispersion from SDSS. First of all, we compute the photometry of mock lenses in gri bands and apply a simple catalog-level neural network to identify a sample of 1050207 galaxies with similar colors and magnitudes as the mocks. Secondly, we train a convolutional neural network (CNN) on Pan-STARRS gri image cutouts to classify this sample and obtain sets of 105760 and 12382 lens candidates with scores pCNN>0.5 and >0.9, respectively. Extensive tests show that CNN performances rely heavily on the design of lens simulations and choice of negative examples for training, but little on the network architecture. Finally, we visually inspect all galaxies with pCNN>0.9 to assemble a final set of 330 high-quality newly-discovered lens candidates while recovering 23 published systems. For a subset, SDSS spectroscopy on the lens central regions proves our method correctly identifies lens LRGs at z~0.1-0.7. Five spectra also show robust signatures of high-redshift background sources and Pan-STARRS imaging confirms one of them as a quadruply-imaged red source at z_s = 1.185 strongly lensed by a foreground LRG at z_d = 0.3155. In the future, we expect that the efficient and automated two-step classification method presented in this paper will be applicable to the deeper gri stacks from the LSST with minor adjustments., Comment: 18 pages and 11 figures (plus appendix), version published in A&A

We present the results of 15 years of monitoring lensed quasars, which was conducted by the COSMOGRAIL programme at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed systems are presented for the first time. We complement our data set with other monitoring data available in the literature to measure the time delays in 18 systems, among which nine reach a relative precision better than 15% for at least one time delay. To achieve this, we developed an automated version of the curve-shifting toolbox PyCS to ensure robust estimation of the time delay in the presence of microlensing, while accounting for the errors due to the imperfect representation of microlensing. We also re-analysed the previously published time delays of RX J1131$-$1231 and HE 0435$-$1223, by adding six and two new seasons of monitoring, respectively, and confirming the previous time-delay measurements. When the time delay measurement is possible, we corrected the light curves of the lensed images from their time delay and present the difference curves to highlight the microlensing signal contained in the data. To date, this is the largest sample of decade-long lens monitoring data, which is useful to measure $H_0$ and the size of quasar accretion discs with microlensing as well as to study quasar variability., Comment: Published in A&A, 31 pages, 8 figures, 3 Tables

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., Comment: 22 pages, 16 figures, submitted to MNRAS

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., Comment: 32 pages, 17 figures, 10 tables. Accepted by MNRAS

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., Comment: (MNRAS in press; color scheme of Figures 1-4 improved with respect to v1)

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., Comment: 22 pages, 9 figures, 8 tables, Accepted to MNRAS. This version: updated Table 4 from the erratum

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

We present the first rest-frame UV population study of 17 heavily reddened, high-luminosity (E(B-V)$_{\rm{QSO}}\gtrsim$ 0.5; L$_{\rm{bol}}>$ 10$^{46}$ergs$^{-1}$) broad-line quasars at $1.5 < z < 2.7$. We combine the first year of deep, optical, ground-based observations from the Dark Energy Survey (DES) with the near infrared VISTA Hemisphere Survey (VHS) and UKIDSS Large Area Survey (ULAS) data, from which the reddened quasars were initially identified. We demonstrate that the significant dust reddening towards the quasar in our sample allows host galaxy emission to be detected at the rest-frame UV wavelengths probed by the DES photometry. By exploiting this reddening effect, we disentangle the quasar emission from that of the host galaxy via spectral energy distribution (SED) fitting. We find evidence for a relatively unobscured, star-forming host galaxy in at least ten quasars, with a further three quasars exhibiting emission consistent with either star formation or scattered light. From the rest-frame UV emission, we derive instantaneous, dust-corrected star formation rates (SFRs) in the range 25 < SFR$_{\rm{UV}}$ < 365 M$_{\odot}$yr$^{-1}$, with an average SFR$_{\rm{UV}}$ = 130 $\pm$ 95 M$_{\odot}$yr$^{-1}$. We find a broad correlation between SFR$_{\rm{UV}}$ and the bolometric quasar luminosity. Overall, our results show evidence for coeval star formation and black hole accretion occurring in luminous, reddened quasars at the peak epoch of galaxy formation., Comment: 21 pages, 8 figures, MNRAS Accepted

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., Comment: 12 pages, 8 figures, MNRAS subm. 30/01/2017. Based on the discovery results of H~Lin et al. (2017), ApJL subm. This paper has been approved for submission by DES, via Collaboration-Wide Review