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

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.

We present the results of 15 years of monitoring of lensed quasars conducted by the COSMOGRAIL program at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed QSOs 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 9 reach a relative precision better than 15% for at least one time delay. To achieve this, we develop 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 imperfect representation of microlensing. We also reanalyze the previously published time delays of RX J1131$-$1231 and HE 0435$-$1223, adding respectively 6 and 2 new seasons of monitoring and confirming the previous time-delay measurements. When the time delay measurement is possible, we correct 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. This is to date the largest sample of decade-long lens monitoring data, useful to measure $H_0$, to measure the size of quasar accretion disks with microlensing, and to study quasar variability.

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

We present two catalogues of active galactic nucleus (AGN) candidates selected from the latest data of two all-sky surveys -- Data Release 2 (DR2) of the \emph{Gaia} mission and the unWISE catalogue of the \emph{Wide-field Infrared Survey Explorer} (\emph{WISE}). We train a random forest classifier to predict the probability of each source in the \emph{Gaia}-unWISE joint sample being an AGN, $P_{\rm RF}$, based on \emph{Gaia} astrometric and photometric measurements and unWISE photometry. The two catalogues, which we designate C75 and R85, are constructed by applying different $P_{\rm RF}$ threshold cuts to achieve an overall completeness of 75\% ($\approx$90\% at \emph{Gaia} $G\leq20$ mag) and reliability of 85\% respectively. The C75 (R85) catalogue contains 2,734,464 (2,182,193) AGN candidates across the effective 36,000 deg$^2$ sky, of which $\approx$0.91 (0.52) million are new discoveries. Photometric redshifts of the AGN candidates are derived by a random forest regressor using \emph{Gaia} and \emph{WISE} magnitudes and colours. The estimated overall photometric redshift accuracy is 0.11. Cross-matching the AGN candidates with a sample of known bright cluster galaxies, we identify a high-probability strongly-lensed AGN candidate system, SDSS\,J1326$+$4806, with a large image separation of 21\farcs06. All the AGN candidates in our catalogues will have $\sim$5-year long light curves from \emph{Gaia} by the end of the mission, and thus will be a great resource for AGN variability studies. Our AGN catalogues will also be helpful in AGN target selections for future spectroscopic surveys, especially ones in the southern hemisphere. The C75 catalogue can be downloaded at https://www.ast.cam.ac.uk/~ypshu/AGN_Catalogues.html.

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

We report the discovery and spectroscopic confirmation of 22 new gravitationally lensed quasars found using $\textit{Gaia}$ data release 2. The selection was made using several techniques: multiple $\textit{Gaia}$ detections around objects in quasar candidate catalogues, modelling of unWISE coadd pixels using $\textit{Gaia}$ astrometry, and $\textit{Gaia}$ detections offset from photometric and spectroscopic galaxies. Spectra of 33 candidates were obtained with the William Herschel Telescope, 22 of which are lensed quasars, 2 highly probably lensed quasars, 5 nearly identical quasar pairs, 1 inconclusive system, and 3 contaminants. Of the 3 confirmed quadruply imaged systems, J2145+6345 is a 2.1 arcsecond separation quad with 4 bright images ($G$=16.86, 17.26, 18.34, 18.56), making it ideal for time delay monitoring. Analysing this new sample alongside known lenses in the Pan-STARRS footprint, and comparing to expected numbers of lenses, we show that, as expected, we are biased towards systems with bright lensing galaxies and low source redshifts. We discuss possible techniques to remove this bias from future searches. A |b|>20 complete sample of lensed quasars detected by $\textit{Gaia}$ and with image separations above 1 arcsecond will provide a valuable statistical sample of around 350 systems. Currently only 96 known lenses satisfy these criteria, yet promisingly, our unWISE modelling technique is able to recover all of these with simple WISE-$\textit{Gaia}$ colour cuts that remove $\sim$80 per cent of previously followed-up contaminants. Finally, we provide an online database of known lenses, quasar pairs, and contaminant systems.

Data mining through large, multi-band ground-based surveys and the recent availability of source catalogues from the {\it Gaia} space mission has led to a surge in the number of known gravitationally lensed quasars, including some unusual systems. One of these newly discovered systems is GraL~J181730853+272940139. This system was reported in Delchambre et al. 2018, and was identified as a quadruply lensed quasar candidate via a machine learning search inside {\it Gaia} Data Release 2, where it consists of three sources with maximum separation $1.8''$. Lemon et al. in prep (private communications) confirmed it spectroscopically as a lensed quasar with source redshift $z=3.07$. A mass model of this system has not been published to date, and the highest resolution imaging data available is from PanSTARRS1 \citep{chambers16}, with seeing $\sim1.0''-1.3''$ in the $grizy$ filters. Here we present the first mass model of J1817+2729, which includes the measured relative astrometry and morphology of the lensing galaxy, using newly acquired Subaru/FOCAS imaging.

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 the discovery, spectroscopic confirmation and preliminary characterisation of 24 gravitationally lensed quasars identified using Gaia observations. Candidates were selected in the Pan-STARRS footprint with quasar-like WISE colours or as photometric quasars from SDSS, requiring either multiple detections in Gaia or a single Gaia detection near a morphological galaxy. The Pan-STARRS grizY images were modelled for the most promising candidates and 60 candidate systems were followed up with the William Herschel Telescope. 13 of the lenses were discovered as Gaia multiples and 10 as single Gaia detections near galaxies. We also discover 1 lens identified through a quasar emission line in an SDSS galaxy spectrum. The lenses have median image separation 2.13 arcsec and the source redshifts range from 1.06 to 3.36. 4 systems are quadruply-imaged and 20 are doubly-imaged. Deep CFHT data reveal an Einstein ring in one double system. We also report 12 quasar pairs, 10 of which have components at the same redshift and require further follow-up to rule out the lensing hypothesis. We compare the properties of these lenses and other known lenses recovered by our search method to a complete sample of simulated lenses to show the lenses we are missing are mainly those with small separations and higher source redshifts. The initial Gaia data release only catalogues all images of ~ 30% of known bright lensed quasars, however the improved completeness of Gaia data release 2 will help find all bright lensed quasars on the sky.

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

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.

We report the discovery, spectroscopic confirmation, and mass modelling of the gravitationally lensed quasar system PS J0630-1201. The lens was discovered by matching a photometric quasar catalogue compiled from Pan-STARRS and WISE photometry to the Gaia DR1 catalogue, exploiting the high spatial resolution of the latter (FWHM $\sim $0.1") to identify the three brightest components of the lens. Follow-up spectroscopic observations with the WHT confirm the multiple objects are quasars at redshift $z_{q}=3.34$. Further follow-up with Keck AO high-resolution imaging reveals that the system is composed of two lensing galaxies and the quasar is lensed into a $\sim$2.8" separation four-image cusp configuration with a fifth image clearly visible, and a 1.0" arc due to the lensed quasar host galaxy. The system is well-modelled with two singular isothermal ellipsoids, reproducing the position of the fifth image. We discuss future prospects for measuring time delays between the images and constraining any offset between mass and light using the faintly detected Einstein arcs associated with the quasar host galaxy.

Gaia's exceptional resolution (FWHM $\sim$ 0.1$^{\prime\prime}$) allows identification and cataloguing of the multiple images of gravitationally lensed quasars. We investigate a sample of 49 known lensed quasars in the SDSS footprint, with image separations less than 2$^{\prime\prime}$, and find that 8 are detected with multiple components in the first Gaia data release. In the case of the 41 single Gaia detections, we generally are able to distinguish these lensed quasars from single quasars when comparing Gaia flux and position measurements to those of Pan-STARRS and SDSS. This is because the multiple images of these lensed quasars are typically blended in ground-based imaging and therefore the total flux and a flux-weighted centroid are measured, which can differ significantly from the fluxes and centroids of the individual components detected by Gaia. We compare the fluxes through an empirical fit of Pan-STARRS griz photometry to the wide optical Gaia bandpass values using a sample of isolated quasars. The positional offsets are calculated from a recalibrated astrometric SDSS catalogue. Applying flux and centroid difference criteria to spectroscopically confirmed quasars, we discover 4 new sub-arcsecond-separation lensed quasar candidates which have two distinct components of similar colour in archival CFHT or HSC data. Our method based on single Gaia detections can be used to identify the $\sim$ 1400 lensed quasars with image separation above 0.5$^{\prime\prime}$, expected to have only one image bright enough to be detected by Gaia.

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

We present the discovery and preliminary characterization of a gravitationally lensed quasar with a source redshift $z_{s}=2.74$ and image separation of $2.9"$ lensed by a foreground $z_{l}=0.40$ elliptical galaxy. Since the images of gravitationally lensed quasars are the superposition of multiple point sources and a foreground lensing galaxy, we have developed a morphology independent multi-wavelength approach to the photometric selection of lensed quasar candidates based on Gaussian Mixture Models (GMM) supervised machine learning. Using this technique and $gi$ multicolour photometric observations from the Dark Energy Survey (DES), near IR $JK$ photometry from the VISTA Hemisphere Survey (VHS) and WISE mid IR photometry, we have identified a candidate system with two catalogue components with $i_{AB}=18.61$ and $i_{AB}=20.44$ comprised of an elliptical galaxy and two blue point sources. Spectroscopic follow-up with NTT and the use of an archival AAT spectrum show that the point sources can be identified as a lensed quasar with an emission line redshift of $z=2.739\pm0.003$ and a foreground early type galaxy with $z=0.400\pm0.002$. We model the system as a single isothermal ellipsoid and find the Einstein radius $\theta_E \sim 1.47"$, enclosed mass $M_{enc} \sim 4 \times 10^{11}$M$_{\odot}$ and a time delay of $\sim$52 days. The relatively wide separation, month scale time delay duration and high redshift make this an ideal system for constraining the expansion rate beyond a redshift of 1.