• We present results of a systematic search for gravitationally lensed quasars in Pan-STARRS1. Our final sample of candidates comprise of 46 new doubles, 2 new quads, and 6 known lenses, which we discovered independently. In the absence of spectroscopy to verify the lensing nature of the candidates, the main sources of contaminants are likely to be quasar pairs and quasar$-$star associations. Amongst the independently discovered quads is 2M1134$-$2103, for which we obtained spectroscopy for the first time, finding a redshift of 2.77 for the quasar. There is evidence for microlensing in at least one image. We perform detailed mass modeling of this system using archival imaging data, and find that the unusually large shear responsible for the diamond-like configuration can mainly be attributed to a faint companion $\sim4''$ away, and to a galaxy group $\sim30''$ away. We also set limits of $z\sim0.5-1.5$ on the redshift of the lensing galaxy, based on its brightness, the image separation of the lensed images, and an analysis of the observed photometric flux ratios.
  • 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.
  • Strong gravitational lenses with measured time delays between the multiple images allow a direct measurement of the time-delay distance to the lens, and thus a measure of cosmological parameters, particularly the Hubble constant, $H_{0}$. We present a blind lens model analysis of the quadruply-imaged quasar lens HE 0435-1223 using deep Hubble Space Telescope imaging, updated time-delay measurements from the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL), a measurement of the velocity dispersion of the lens galaxy based on Keck data, and a characterization of the mass distribution along the line of sight. HE 0435-1223 is the third lens analyzed as a part of the $H_{0}$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) project. We account for various sources of systematic uncertainty, including the detailed treatment of nearby perturbers, the parameterization of the galaxy light and mass profile, and the regions used for lens modeling. We constrain the effective time-delay distance to be $D_{\Delta t} = 2612_{-191}^{+208}~\mathrm{Mpc}$, a precision of 7.6%. From HE 0435-1223 alone, we infer a Hubble constant of $H_{0} = 73.1_{-6.0}^{+5.7}~\mathrm{km~s^{-1}~Mpc^{-1}}$ assuming a flat $\Lambda$CDM cosmology. The cosmographic inference based on the three lenses analyzed by H0LiCOW to date is presented in a companion paper (H0LiCOW Paper V).
  • The empirical correlation between the mass of a super-massive black hole (MBH) and its host galaxy properties is widely considered to be evidence of their co-evolution. A powerful way to test the co-evolution scenario and learn about the feedback processes linking galaxies and nuclear activity is to measure these correlations as a function of redshift. Unfortunately, currently MBH can only be estimated in active galaxies at cosmological distances. At these distances, bright active galactic nuclei (AGN) can outshine the host galaxy, making it extremely difficult to measure the host's luminosity. Strongly lensed AGNs provide in principle a great opportunity to improve the sensitivity and accuracy of the host galaxy luminosity measurements as the host galaxy is magnified and more easily separated from the point source, provided the lens model is sufficiently accurate. In order to measure the MBH-L correlation with strong lensing, it is necessary to ensure that the lens modelling is accurate, and that the host galaxy luminosity can be recovered to at least a precision and accuracy better than that of the typical MBH measurement. We carry out extensive and realistic simulations of deep Hubble Space Telescope observations of lensed AGNs obtained by our collaboration. We show that the host galaxy luminosity can be recovered with better accuracy and precision than the typical uncertainty on MBH(~ 0.5 dex) for hosts as faint as 2-4 magnitudes dimmer than the AGN itself. Our simulations will be used to estimate bias and uncertainties on the actual measurements to be presented in a future paper.
  • Based on spectroscopy and multiband wide-field observations of the gravitationally lensed quasar HE 0435-1223, we determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for this system. We measure the under/overdensity of the line of sight towards the lens system and compare it to the average line of sight throughout the universe, determined by using the CFHTLenS as a control field. Aiming to constrain $\kappa_\mathrm{ext}$ as tightly as possible, we determine under/overdensities using various combinations of relevant informative weighing schemes for the galaxy counts, such as projected distance to the lens, redshift, and stellar mass. We then convert the measured under/overdensities into a $\kappa_\mathrm{ext}$ distribution, using ray-tracing through the Millennium Simulation. We explore several limiting magnitudes and apertures, and account for systematic and statistical uncertainties relevant to the quality of the observational data, which we further test through simulations. Our most robust estimate of $\kappa_\mathrm{ext}$ has a median value $\kappa^\mathrm{med}_\mathrm{ext} = 0.004$ and a standard deviation of $\sigma_\kappa = 0.025$. The measured $\sigma_\kappa$ corresponds to $2.5\%$ uncertainty on the time delay distance, and hence the Hubble constant $H_0$ inference from this system. The median $\kappa^\mathrm{med}_\mathrm{ext}$ value is robust to $\sim0.005$ (i.e. $\sim0.5\%$ on $H_0$) regardless of the adopted aperture radius, limiting magnitude and weighting scheme, as long as the latter incorporates galaxy number counts, the projected distance to the main lens, and a prior on the external shear obtained from mass modeling. The availability of a well-constrained $\kappa_\mathrm{ext}$ makes \hequad\ a valuable system for measuring cosmological parameters using strong gravitational lens time delays.
  • We present the results of an imaging observation campaign conducted with the Subaru Telescope adaptive optics system (IRCS+AO188) on 28 gravitationally lensed quasars (23 doubles, 1 quad, and 1 possible triple, and 3 candidates) from the SDSS Quasar Lens Search. We develop a novel modelling technique that fits analytical and hybrid point spread functions (PSFs), while simultaneously measuring the relative astrometry, photometry, as well as the lens galaxy morphology. We account for systematics by simulating the observed systems using separately observed PSF stars. The measured relative astrometry is comparable with that typically achieved with the Hubble Space Telescope, even after marginalizing over the PSF uncertainty. We model for the first time the quasar host galaxies in 5 systems, without a-priory knowledge of the PSF, and show that their luminosities follow the known correlation with the mass of the supermassive black hole. For each system, we obtain mass models far more accurate than those previously published from low-resolution data, and we show that in our sample of lensing galaxies the observed light profile is more elliptical than the mass, for ellipticity > 0.25. We also identify eight doubles for which the sources of external and internal shear are more reliably separated, and should therefore be prioritized in monitoring campaigns aimed at measuring time-delays in order to infer the Hubble constant.
  • We report the discovery of 13 confirmed two-image quasar lenses from a systematic search for gravitationally lensed quasars in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). We adopted a methodology similar to that used in the SDSS Quasar Lens Search (SQLS). In addition to the confirmed lenses, we report 11 quasar pairs with small angular separations ($\lesssim$2") confirmed from our spectroscopy, which are either projected pairs, physical binaries, or possibly quasar lens systems whose lens galaxies have not yet been detected. The newly discovered quasar lens system, SDSS J1452+4224 at zs$\approx$4.8 is one of the highest redshift multiply imaged quasars found to date. Furthermore, we have over 50 good lens candidates yet to be followed up. Owing to the heterogeneous selection of BOSS quasars, the lens sample presented here does not have a well-defined selection function.
  • We present spectroscopic confirmation of two new lensed quasars via data obtained at the 6.5m Magellan/Baade Telescope. The lens candidates have been selected from the Dark Energy Survey (DES) and WISE based on their multi-band photometry and extended morphology in DES images. Images of DES J0115-5244 show two blue point sources at either side of a red galaxy. Our long-slit data confirm that both point sources are images of the same quasar at $z_{s}=1.64.$ The Einstein Radius estimated from the DES images is $0.51$". DES J2200+0110 is in the area of overlap between DES and the Sloan Digital Sky Survey (SDSS). Two blue components are visible in the DES and SDSS images. The SDSS fiber spectrum shows a quasar component at $z_{s}=2.38$ and absorption compatible with Mg II and Fe II at $z_{l}=0.799$, which we tentatively associate with the foreground lens galaxy. The long-slit Magellan spectra show that the blue components are resolved images of the same quasar. The Einstein Radius is $0.68$" corresponding to an enclosed mass of $1.6\times10^{11}\,M_{\odot}.$ Three other candidates were observed and rejected, two being low-redshift pairs of starburst galaxies, and one being a quasar behind a blue star. These first confirmation results provide an important empirical validation of the data-mining and model-based selection that is being applied to the entire DES dataset.
  • We present the result of Subaru Telescope multi-band adaptive optics observations of the complex gravitationally lensed quasar SDSS J1405+0959, which is produced by two lensing galaxies. These observations reveal dramatically enhanced morphological detail, leading to the discovery of an additional object 0. 26'' from the secondary lensing galaxy, as well as three collinear clumps located in between the two lensing galaxies. The new object is likely to be the third quasar image, although the possibility that it is a galaxy cannot be entirely excluded. If confirmed via future observations, it would be the first three image lensed quasar produced by two galaxy lenses. In either case, we show based on gravitational lensing models and photometric redshift that the collinear clumps represent merging images of a portion of the quasar host galaxy, with a magnification factor of 15 - 20, depending on the model.
  • We derive the average mass profile of elliptical galaxies from the ensemble of 161 strong gravitational lens systems selected from several surveys, assuming that the mass profile scales with the stellar mass and effective radius of each lensing galaxy. The total mass profile is well fitted by a power-law \rho(r) \propto r^\gamma with best-fit slope \gamma = -2.11+/-0.05. The decomposition of the total mass profile into stellar and dark matter distributions is difficult due to a fundamental degeneracy between the stellar initial mass function (IMF) and the dark matter fraction f_DM. We demonstrate that this IMF-f_DM degeneracy can be broken by adding direct stellar mass fraction measurements by quasar microlensing observations. Our best-fit model prefers the Salpeter IMF over the Chabrier IMF, and a smaller central dark matter fraction than that predicted by adiabatic contraction models.
  • We report the discovery of a pair of quasars at $z=1.487$, with a separation of $8\farcs585\pm0\farcs002$. Subaru Telescope infrared imaging reveals the presence of an elliptical and a disk-like galaxy located almost symmetrically between the quasars, creating a cross-like configuration. Based on absorption lines in the quasar spectra and the colors of the galaxies, we estimate that both galaxies are located at redshift $z=0.899$. This, as well as the similarity of the quasar spectra, suggests that the system is a single quasar multiply imaged by a galaxy group or cluster acting as a gravitational lens, although the possibility of a binary quasar cannot be fully excluded. We show that the gravitational lensing hypothesis implies these galaxies are not isolated, but must be embedded in a dark matter halo of virial mass $\sim 4 \times 10^{14}\ h_{70}^{-1}\ {M}_\odot$ assuming an NFW model with a concentration parameter of $c_{vir}=6$, or a singular isothermal sphere profile with a velocity dispersion of $\sim 670$ km s$^{-1}$. We place constraints on the location of the dark matter halo, as well as the velocity dispersions of the galaxies. In addition, we discuss the influence of differential reddening, microlensing and intrinsic variability on the quasar spectra and broadband photometry.
  • We present the final statistical sample of lensed quasars from the Sloan Digital Sky Survey (SDSS) Quasar Lens Search (SQLS). The well-defined statistical lens sample consists of 26 lensed quasars brighter than i=19.1 and in the redshift range of 0.6<z<2.2 selected from 50,836 spectroscopically confirmed quasars in the SDSS Data Release 7 (DR7), where we restrict the image separation range to 1"<\theta<20" and the i-band magnitude differences in two image lenses to be smaller than 1.25 mag. The SDSS DR7 quasar catalog also contains 36 additional lenses identified with various techniques. In addition to these lensed quasars, we have identified 81 pairs of quasars from follow-up spectroscopy, 26 of which are physically associated binary quasars. The statistical lens sample covers a wide range of image separations, redshifts, and magnitudes, and therefore is suitable for systematic studies of cosmological parameters and surveys of the structure and evolution of galaxies and quasars.
  • We present a statistical analysis of the final lens sample from the Sloan Digital Sky Survey Quasar Lens Search (SQLS). The number distribution of a complete subsample of 19 lensed quasars selected from 50,836 source quasars is compared with theoretical expectations, with particular attention to the selection function. Assuming that the velocity function of galaxies does not evolve with redshift, the SQLS sample constrains the cosmological constant to \Omega_\Lambda=0.79^{+0.06}_{-0.07}(stat.)^{+0.06}_{-0.06}(syst.) for a flat universe. The dark energy equation of state is found to be consistent with w=-1 when the SQLS is combined with constraints from baryon acoustic oscillation (BAO) measurements or results from the Wilkinson Microwave Anisotropy Probe (WMAP). We also obtain simultaneous constraints on cosmological parameters and redshift evolution of the galaxy velocity function, finding no evidence for redshift evolution at z<1 in any combinations of constraints. For instance, number density evolution quantified as \nu_n=d\ln\phi_*/d\ln(1+z) and the velocity dispersion evolution \nu_\sigma=d\ln\sigma_*/d\ln(1+z) are constrained to \nu_n=1.06^{+1.36}_{-1.39}(stat.)^{+0.33}_{-0.64}(syst.) and \nu_\sigma=-0.05^{+0.19}_{-0.16}(stat.)^{+0.03}_{-0.03}(syst.) respectively when the SQLS result is combined with BAO and WMAP for flat models with a cosmological constant. We find that a significant amount of dark energy is preferred even after fully marginalizing over the galaxy evolution parameters. Thus the statistics of lensed quasars robustly confirm the accelerated cosmic expansion.
  • The quasar SDSS J133401.39+331534.3 at z = 2.426 is found to be a two-image gravitationally lensed quasar with the image separation of 0.833. The object is first identified as a lensed quasar candidate in the Sloan Digital Sky Survey Quasar Lens Search, and then confirmed as a lensed system from follow-up observations at the Subaru and University of Hawaii 2.2-meter telescopes. We estimate the redshift of the lensing galaxy to be 0.557 based on absorption lines in the quasar spectra as well as the color of the galaxy. In particular, we observe the system with the Subaru Telescope AO188 adaptive optics with laser guide star, in order to derive accurate astrometry, which well demonstrates the usefulness of the laser guide star adaptive optics imaging for studying strong lens systems. Our mass modeling with improved astrometry implies that a nearby bright galaxy $\sim 4"$ apart from the lensing galaxy is likely to affect the lens potential.