The H0LiCOW collaboration inferred via gravitational lensing time delays a Hubble constant $H_0=73.3^{+1.7}_{-1.8}$ km s$^{-1}{\rm Mpc}^{-1}$, describing deflector mass density profiles by either a power-law or stars plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in $H_0$. We quantify any potential effect of the MST with flexible mass models that are maximally degenerate with H0. Our calculation is based on a new hierarchical approach in which the MST is only constrained by stellar kinematics. The approach is validated on hydrodynamically simulated lenses. We apply the method to the TDCOSMO sample of 7 lenses (6 from H0LiCOW) and measure $H_0=74.5^{+5.6}_{-6.1}$ km s$^{-1}{\rm Mpc}^{-1}$. In order to further constrain the deflector mass profiles, we then add imaging and spectroscopy for 33 strong gravitational lenses from the SLACS sample. For 9 of the SLAC lenses we use resolved kinematics to constrain the stellar anisotropy. From the joint analysis of the TDCOSMO+SLACS sample, we measure $H_0=67.4^{+4.1}_{-3.2}$ km s$^{-1}{\rm Mpc}^{-1}$, assuming that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. While our new analysis does not statistically invalidate the mass profile assumptions by H0LiCOW, and thus their $H_0$ measurement relying on those, it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on $H_0$ derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data, chiefly spatially resolved kinematics of lens galaxies.

MAAT is proposed as a visitor mirror-slicer optical system that will allow the OSIRIS spectrograph on the 10.4-m Gran telescopio CANARIAS (GTC) the capability to perform Integral Field Spectroscopy (IFS) over a seeing-limited FoV 14.20''x10'' with a slice width of 0.303''. MAAT@GTC will enhance the resolution power of OSIRIS by 1.6 times as compared to its 0.6'' wide long-slit. All the eleven OSIRIS grisms and volume-phase holographic gratings will be available to provide broad spectral coverage with moderate resolution (R=600 up to 4100) in the 3600 - 10000 {\AA} wavelength range. MAAT unique observing capabilities will broaden its use to the needs of the GTC community to unveil the nature of most striking phenomena in the universe well beyond time-domain astronomy. The GTC equipped with OSIRIS+MAAT will also play a fundamental role in synergy with other facilities, some of them operating on the northern ORM at La Palma. This White Paper presents the different aspects of MAAT@GTC - including scientific and technical specifications, outstanding science cases, and an outline of the instrument concept.

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

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.

The MATLAS deep imaging survey has uncovered a plethora of dwarf galaxies in the low density environment it has mapped. A fraction of them are unusually extended and have a low-surface brightness. Among these so-called ultra-diffuse galaxies, a few seem to host an excess of globular clusters. With the integral-field unit spectrograph MUSE we have observed one of these galaxies - MATLAS J15052031+0148447 (MATLAS-2019) - located towards the nearby group NGC 5846 and measured its systemic velocity,age, and metallicity, and that of its globular clusters candidates. For the stellar body of MATLAS-2019 we derive a metallicity of -1.33+0.19-0.01 dex and an age of 11.2+1.8-0.8 Gyr. For some of the individual GCs and the stacked GC population, we derive consistent ages and metallicities. From the 11 confirmed globular clusters and using a Markov Chain Monte Carlo approach we derived a dynamical mass-to-light ratio of 4.2+8.6-3.4M/L. This is at the lower end of the luminosity-mass scaling relation defined by the Local Group dwarf galaxies. Furthermore, we couldn't confirm nor reject the possibility of a rotational component of the GC system. If present, this would further modify the inferred mass. Follow-up observations of the globular cluster population and of the stellar body of the galaxy are needed to assess whether this galaxy is lacking dark matter like it was suggested for the pair of dwarf galaxies in the field of NGC 1052, or if this is a miss-interpretation arising from systematic uncertainties of the method commonly used for these systems and the large uncertainties of the individual globular cluster velocities.

One of the main challenges in using high redshift active galactic nuclei to study the correlations between the mass of the supermassive Black Hole (MBH) and the properties of their active host galaxies is instrumental resolution. Strong lensing magnification effectively increases instrumental resolution and thus helps to address this challenge. In this work, we study eight strongly lensed active galactic nuclei (AGN) with deep Hubble Space Telescope imaging, using the lens modelling code Lenstronomy to reconstruct the image of the source. Using the reconstructed brightness of the host galaxy, we infer the host galaxy stellar mass based on stellar population models. MBH are estimated from broad emission lines using standard methods. Our results are in good agreement with recent work based on non-lensed AGN, providing additional evidence that the correlation evolves over cosmic time. At the moment, the sample size of lensed AGN is small and thus they provide mostly a consistency check on systematic errors related to resolution for the non-lensed AGN. However, the number of known lensed AGN is expected to increase dramatically in the next few years, through dedicated searches in ground and space based wide field surveys, and they may become a key diagnostic of black hole and galaxy co-evolution.

Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Calibrating the Baryon Acoustic Oscillations (BAO) and Type Ia supernovae with combinations of the distance indicators, we obtain a joint and self-consistent measurement of $H_0$ and $r_d$ at low redshift, independent of cosmological models and CMB inference. In an attempt to alleviate the tension between late-time and CMB-based measurements, we consider four extensions of the standard $\Lambda$CDM model. The sound horizon from our different measurements is $r_d=(137\pm3^{stat.}\pm2^{syst.})$~Mpc. Depending on the adopted distance indicators, the $combined$ tension in $H_0$ and $r_d$ ranges between 2.3 and 5.1 $\sigma$. We find that modifications of $\Lambda$CDM that change the physics after recombination fail to solve the problem, for the reason that they only resolve the tension in $H_0$, while the tension in $r_d$ remains unchanged. Pre-recombination extensions (with early dark energy or the effective number of neutrinos $\rm{N}_{\rm{eff}}=3.24 \pm 0.16$) are allowed by the data, unless the calibration from Cepheids is included. Results from time-delay lenses are consistent with those from distance-ladder calibrations and point to a discrepancy between absolute distance scales measured from the CMB (assuming the standard cosmological model) and late-time observations. New proposals to resolve this tension should be examined with respect to reconciling not only the Hubble constant but also the sound horizon derived from the CMB and other cosmological probes.

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.

The matter sound horizon can be inferred from the cosmic microwave background within the Standard Model. Independent direct measurements of the sound horizon are then a probe of possible deviations from the Standard Model. We aim at measuring the sound horizon $r_s$ from low-redshift indicators, which are completely independent of CMB inference. We used the measured product $H(z)r_s$ from baryon acoustic oscillations (BAO) together with supernovae~\textsc{I}a to constrain $H(z)/H_{0}$ and time-delay lenses analysed by the H0LiCOW collaboration to anchor cosmological distances ($\propto H_{0}^{-1}$). {Additionally, we investigated the influence of adding a sample of quasars with higher redshift with standardisable UV-Xray luminosity distances. We adopted polynomial expansions in $H(z)$ or in comoving distances} so that our inference was completely independent of any cosmological model on which the expansion history might be based. Our measurements are independent of Cepheids and systematics from peculiar motions {to within percent-level accuracy.} The inferred sound horizon $r_s$ varies between $(133 \pm 8)$~Mpc and $(138 \pm 5)$~Mpc across different models. The discrepancy with CMB measurements is robust against model choice. Statistical uncertainties are comparable to systematics. The combination of time-delay lenses, supernovae, and BAO yields a distance ladder that is independent of cosmology (and of Cepheid calibration) and a measurement of $r_s $ that is independent of the CMB. These cosmographic measurements are then a competitive test of the Standard Model, regardless of the hypotheses on which the cosmology is based.

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 the photometric properties of 2210 newly identified dwarf galaxy candidates in the MATLAS fields. The Mass Assembly of early Type gaLAxies with their fine Structures (MATLAS) deep imaging survey mapped $\sim$142 deg$^2$ of the sky around nearby isolated early type galaxies using MegaCam on the Canada-France-Hawaii Telescope, reaching surface brightnesses of $\sim$ 28.5 - 29 in the g-band. The dwarf candidates were identified through a direct visual inspection of the images and by visually cleaning a sample selected using a partially automated approach, and were morphologically classified at the time of identification. Approximately 75% of our candidates are dEs, indicating that a large number of early type dwarfs also populate low density environments, and 23.2% are nucleated. Distances were determined for 13.5% of our sample using pre-existing $z_{spec}$ measurements and HI detections. We confirm the dwarf nature for 99% of this sub-sample based on a magnitude cut $M_g$ = -18. Additionally, most of these ($\sim$90%) have relative velocities suggesting that they form a satellite population around nearby massive galaxies rather than an isolated field sample. Assuming that the candidates over the whole survey are satellites of the nearby galaxies, we demonstrate that the MATLAS dwarfs follow the same scaling relations as dwarfs in the Local Group as well as the Virgo and Fornax clusters. We also find that the nucleated fraction increases with $M_g$, and find evidence of a morphology-density relation for dwarfs around isolated massive galaxies.

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.

We describe the observation and confirmation of \nbconfirmtext\ new strong gravitational lenses discovered in Year 1 data from the Dark Energy Survey (DES). We created candidate lists based on a) galaxy group and cluster samples and b) photometrically selected galaxy samples. We selected 46 candidates through visual inspection and then used the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South telescope to acquire spectroscopic follow-up of 21 of these candidates. Through analysis of this spectroscopic follow-up data, we confirmed nine new lensing systems and rejected 2 candidates, but the analysis was inconclusive on 10 candidates. For each of the confirmed systems, we report measured spectroscopic properties, estimated \einsteinradiussub, and estimated enclosed masses. The sources that we targeted have an i-band surface brightness range of iSB ~ 22 - 24 mag arcsec^2 and a spectroscopic redshift range of zspec ~0.8 - 2.6. The lens galaxies have a photometric redshift range of zlens ~ 0.3 - 0.7. The lensing systems range in image-lens separation 2 - 9 arcsec and in enclosed mass 10^12 - 10^13 Msol.

We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $\Lambda$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1\sigma$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3\sigma$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $\Lambda$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $\Lambda$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints.

We present an improved determination of the total mass distribution of three CLASH/HFF massive clusters, MACS J1206.2-0847 (z=0.44), MACS J0416.1-2403 (z=0.40), Abell S1063 (z=0.35). We specifically reconstruct the sub-halo mass component with robust stellar kinematics information of cluster galaxies, in combination with precise strong lensing models based on large samples of spectroscopically identified multiple images. We use VLT/MUSE integral-field spectroscopy in the cluster cores to measure the stellar velocity dispersion, $\sigma$, of 40-60 member galaxies per cluster, covering 4-5 magnitudes to $m_{F160W}\simeq 21.5$. We verify the robustness and quantify the accuracy of the velocity dispersion measurements with extensive spectral simulations. With these data, we determine the normalization and slope of the galaxy $L\mbox{-}\sigma$ Faber-Jackson relation in each cluster and use these parameters as a prior for the scaling relations of the sub-halo population in the mass distribution modeling. When compared to our previous lens models, the inclusion of member galaxies' kinematics provides a similar precision in reproducing the positions of the multiple images. However, the inherent degeneracy between the central effective velocity dispersion, $\sigma_0$, and truncation radius, $r_{cut}$, of sub-halos is strongly reduced, thus significantly alleviating possible systematics in the measurements of sub-halo masses. The three independent determinations of the $\sigma_0\mbox{-}r_{cut}$ scaling relation in each cluster are found to be fully consistent, enabling a statistical determination of sub-halo sizes as a function of $\sigma_0$, or halo masses. We derive galaxy central velocity dispersion functions of the three clusters and found them in agreement with each other. Sub-halo mass functions determined with this method can be compared with those obtained from cosmological simulations.

Galaxies and galaxy groups located along the line of sight towards gravitationally lensed quasars produce high-order perturbations of the gravitational potential at the lens position. When these perturbation are too large, they can induce a systematic error on $H_0$ of a few-percent if the lens system is used for cosmological inference and the perturbers are not explicitly accounted for in the lens model. In this work, we present a detailed characterization of the environment of the lens system WFI2033-4723 ($z_{\rm src} = 1.662$, $z_{\rm lens}$ = 0.6575), one of the core targets of the H0LICOW project for which we present cosmological inferences in a companion paper (Rusu et al. 2019). We use the Gemini and ESO-Very Large telescopes to measure the spectroscopic redshifts of the brightest galaxies towards the lens, and use the ESO-MUSE integral field spectrograph to measure the velocity-dispersion of the lens ($\sigma_{\rm {los}}= 250^{+15}_{-21}$ km/s) and of several nearby galaxies. In addition, we measure photometric redshifts and stellar masses of all galaxies down to $i < 23$ mag, mainly based on Dark Energy Survey imaging (DR1). Our new catalog, complemented with literature data, more than doubles the number of known galaxy spectroscopic redshifts in the direct vicinity of the lens, expanding to 116 (64) the number of spectroscopic redshifts for galaxies separated by less than 3 arcmin (2 arcmin) from the lens. Using the flexion-shift as a measure of the amplitude of the gravitational perturbation, we identify 2 galaxy groups and 3 galaxies that require specific attention in the lens models. The ESO MUSE data enable us to measure the velocity-dispersions of three of these galaxies. These results are essential for the cosmological inference analysis presented in Rusu et al. (2019).

We present the first public data release of the VISTA EXtension to Auxiliary Surveys (VEXAS), comprising of 9 cross-matched multi-wavelength photometric catalogs where each object has a match in at least two surveys. We aim at a spatial coverage as uniform as possible in the multi-wavelength sky, with the purpose of providing the astronomical community with reference magnitudes and colours for various scientific uses, including: object classification (e.g. quasars, galaxies, and stars; high-z galaxies, white dwarfs, etc.); photometric redshifts of large galaxy samples; searches of exotic objects such as, for example, extremely red objects and lensed quasars. We have cross-matched the wide-field VISTA catalogs (the VISTA Hemisphere Survey and the VISTA Kilo Degree Infrared Galaxy Survey) with the AllWISE mid-infrared Survey, requiring that a match exists within 10 arcsec. We have further matched this table with X-Ray and radio data (ROSAT, XMM, SUMSS). We also performed a second cross-match between VISTA and AllWISE, with a smaller matching radius (3"), including WISE magnitudes. We have then cross-matched this resulting table ($\approx138\times10^6$ objects) with three photometric wide-sky optical deep surveys (DES, SkyMapper, PanSTARRS). We finally include matches to objects with spectroscopic follow-up by the SDSS and 6dFGS. To demonstrate the power of all-sky multi-wavelength cross-match tables, we show two examples of scientific applications of VEXAS, in particular using the publicly released tables to discover strong gravitational lenses (beyond the reach of previous searches), and to build a statistically large sample of extremely red objects. The VEXAS catalog is currently the widest and deepest, public, optical-to-IR photometric and spectroscopic database in the Southern Hemisphere.

We revisit the claimed tension, or lack thereof, of measured values of the Hubble-Lema\^{i}tre parameter $H_0$ from Cosmic Microwave Background (CMB) data and low-redshift indicators. Baryon Acoustic Oscillations (BAO) rely on the scale of the sound horizon at recombination $r_s$ to convert angular measurements into angular-diameter distances, so fixing $r_s$ from CMB measurements already constrains $H_0.$ If departures from concordance cosmology are to be constrained, truly independent measurements of $H_0$ are needed. We use the angular-diameter distances to three time-delay lenses from the H0LiCOW collaboration to calibrate the distance ladder, combine them with relative distances from Supernovae Ia and BAO, leaving $r_s$ completely free, and provide the inferred coefficients ($q_{0},j_{0},s_{0}$) in the polynomial expansion of H(z). We obtain $H_{0}r_{s}=(9895\pm161)$km/s and $H_0=(72\pm7)$km/s/Mpc. Combined with $H_0$ from H0LiCOW, then $r_s=(137\pm4.5)$Mpc is consistent with previous work and systematically lower than the CMB-inferred value. Our results are independent of the adopted cosmology, and removing Supernovae with z<0.1 has a negligible effect.

The KiDS Strongly lensed QUAsar Detection project (KiDS-SQuaD) aims at finding as many previously undiscovered gravitational lensed quasars as possible in the Kilo Degree Survey. This is the second paper of this series where we present a new, automatic object classification method based on machine learning technique. The main goal of this paper is to build a catalogue of bright extragalactic objects (galaxies and quasars), from the KiDS Data Release 4, with a minimum stellar contamination, preserving the completeness as much as possible, to then apply morphological methods to select reliable gravitationally lensed (GL) quasar candidates. After testing some of the most used machine learning algorithms, decision trees based classifiers, we decided to use CatBoost, that was specifically trained with the aim of creating a sample of extragalactic sources as clean as possible from stars. We discuss the input data, define the training sample for the classifier, give quantitative estimates of its performances, and finally describe the validation results with Gaia DR2, AllWISE, and GAMA catalogues. We have built and make available to the scientific community the KiDS Bright EXtraGalactic Objects catalogue (KiDS-BEXGO), specifically created to find gravitational lenses. This is made of $\approx6$ millions of sources classified as quasars ($\approx 200\,000$) and galaxies ($\approx 5.7$M), up to $r<22^m$. From this catalog we selected 'Multiplets': close pairs of quasars or galaxies surrounded by at least one quasar, presenting the 12 most reliable gravitationally lensed quasar candidates, to demonstrate the potential of the catalogue, which will be further explored in a forthcoming paper. We compared our search to the previous one, presented in the first paper from this series, showing that employing a machine learning method decreases the stars-contaminators within the GL candidates.

The recently discovered discrepancy between galaxy mass measurements from weak lensing and predictions from abundance matching questions our understanding of cosmology, or of the galaxy-halo connection, or of both. We re-examined this tension by considering, as models, different cosmological simulations in the Illustris suite. We produced excess profiles $R\Delta\Sigma$ from subhalo snapshots at different redshifts in Illustris-1 and IllustrisTNG (TNG100 and TNG300) simulations, enabling a direct comparison with weak-lensing measurements. We separate the individual contributions of stars, dark matter and gas within $\approx1$ Mpc (comoving length), beyond which correlated two-halo terms dominate. The mismatch between measurements and predictions is more severe than in previous studies: $R\Delta\Sigma$ profiles from IllustrisTNG are $\approx2$ times higher than the measured ones. Contrary to abundance matching results, the mismatch is mostly unchanged with increasing redshifts. The contribution of gas to the $R\Delta\Sigma$ profiles is $5-10\%$ over the scales dominated by one-halo terms. Different procedures to link stellar and halo masses (abundance matching, cosmological simulations) are still significantly discrepant with weak lensing measurements, but their trends are different. Therefore, the change in cosmological parameters advocated through abundance-matching arguments may not resolve this tension. Also, current criteria to select isolated massive galaxies in simulations are susceptible to resolution issues and may not correspond to observational criteria. The (currently subdominant) contribution of gas is non-negligible, and even if the major discrepancy within stellar and halo masses is resolved, it will be an appreciable source of systematics in the LSST era, when uncertainties on the $R\Delta\Sigma$ profiles are expected to be $\approx10$ times smaller.

We present a blind time-delay strong lensing (TDSL) cosmographic analysis of the doubly imaged quasar SDSS 1206+4332. We combine the relative time delay between the quasar images, Hubble Space Telescope imaging, the Keck stellar velocity dispersion of the lensing galaxy, and wide-field photometric and spectroscopic data of the field to constrain two angular diameter distance relations. The combined analysis is performed by forward modelling the individual data sets through a Bayesian hierarchical framework, and it is kept blind until the very end to prevent experimenter bias. After unblinding, the inferred distances imply a Hubble constant $H_0 = 68.8^{+5.4}_{-5.1}$ kms$^{-1}$Mpc$^{-1}$, assuming a flat Lambda cold dark matter cosmology with uniform prior on $\Omega_{\rm m}$ in [0.05, 0.5]. The precision of our cosmographic measurement with the doubly imaged quasar SDSS 1206+4332 is comparable with those of quadruply imaged quasars and opens the path to perform on selected doubles the same analysis as anticipated for quads. Our analysis is based on a completely independent lensing code than our previous three H0LiCOW systems and the new measurement is fully consistent with those. We provide the analysis scripts paired with the publicly available software to facilitate independent analysis. The consistency between blind measurements with independent codes provides an important sanity check on lens modelling systematics. By combining the likelihoods of the four systems under the same prior, we obtain $H_0 = 72.5^{+2.1}_{-2.3}$kms$^{-1}$Mpc$^{-1}$. This measurement is independent of the distance ladder and other cosmological probes.

The so-called ultra-diffuse galaxy NGC~1052-DF2 was announced to be a galaxy lacking dark matter based on a spectroscopic study of its constituent globular clusters. Here we present the first spectroscopic analysis of the stellar body of this galaxy using the MUSE integral-field spectrograph at the (ESO) Very Large Telescope. The MUSE datacube simultaneously provides DF2's stellar velocity field and systemic velocities for seven globular clusters (GCs). We further discovered three planetary nebulae (PNe) that are likely part of this galaxy. While five of the clusters had velocities measured in the literature, we were able to confirm the membership of two more candidates through precise radial velocity measurements, which increases the measured specific frequency of GCs in DF2. The mean velocity of the diffuse stellar body, 1792.9$^{-1.8}_{+1.4}$~\kms, is consistent with the mean globular cluster velocity. We detect a weak but significant velocity gradient within the stellar body, with a kinematic axis close to the photometric major-axis, making it a prolate-like rotator. We estimate a velocity dispersion from the clusters and PNe of $\sigma_{\mathrm{int}}=10.6^{+3.9}_{-2.3}$~\kms. The velocity dispersion $\sigma_{\rm{DF2}\star}$(\re) for the stellar body within one effective radius is $10.8^{-4.0}_{+3.2}$~\kms. Considering various sources of systemic uncertainties this central value varies between 5 and 13~\kms, and we conservatively report a 95\% confidence upper limit to the dispersion within one \re\ of 21~\kms. We provide updated mass estimates based on these dispersions corresponding to the different distances to NGC~1052-DF2 that have been reported in the recent literature.

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.

NGC 1052-DF2, an ultra diffuse galaxy (UDG), has been the subject of intense debate. Its alleged absence of dark matter, and the brightness and number excess of its globular clusters (GCs) at an initially assumed distance of 20Mpc, suggested a new formation channel for UDGs. We present the first systematic spectroscopic analysis of both the stellar body and the GCs (six of which were previously known, and one newly confirmed member) of this galaxy using MUSE@VLT. Even though NGC 1052-DF2 does not show any spatially extended emission lines we report the discovery of three planetary nebulae (PNe). We conduct full spectral fitting on the UDG and the stacked spectra of all GCs. The UDG's stellar population is old, 8.9$\pm$1.5 Gyr, metal-poor, with [M/H] = $-$1.07$\pm$0.12 with little or no $\alpha$-enrichment. The stacked spectrum of all GCs indicates a similar age of 8.9$\pm$1.8 Gyr, but lower metallicity, with [M/H] = $-$1.63$\pm$0.09, and similarly low $\alpha$-enrichment. There is no evidence for a variation of age and metallicity in the GC population with the available spectra. The significantly more metal-rich stellar body with respect to its associated GCs, the age of the population, its metallicity and alpha enrichment, are all in line with other dwarf galaxies. NGC 1052-DF2 thus falls on the same empirical mass-metallicity relation as other dwarfs, for the full distance range assumed in the literature. We find that both debated distance estimates (13 and 20 Mpc) are similarly likely, given the three discovered PNe.

Quadruply lensed quasars are extremely rare objects, but incredibly powerful cosmological tools. Only few dozen are known in the whole sky. Here we present the spectroscopic confirmation of two new quadruplets WG0214-2105 and WG2100-4452 discovered by Agnello & Spiniello (2018) within the Dark Energy Survey (DES) public footprints. We have conducted spectroscopic follow-up of these systems with the Southern African Large Telescope as part of a program that aims at confirming the largest possible number of optically selected strong gravitational lensing systems in the Equatorial and Southern Hemisphere. For both systems, we present the spectra for the sources and deflectors that allowed us to estimate the source redshifts and unambiguously confirm their lensing nature. For the brighter deflector (WG2100-4452), we measure the stellar velocity dispersion from the spectrum. We also obtain photometry for both lenses, directly from DES multi-band images, isolating the lens galaxies from the quasar images. One of the quadruplets, WG0214-2105, was also observed by Pan-STARRS, allowing us to estimate the apparent brightness of each quasar image at two different epochs, and thus to find evidence for flux variability. This result could suggest a microlensing event for the faintest components, although intrinsic variability cannot be excluded with only two epochs. Finally, we present simple lens models for both quadruplets, obtaining Einstein radii, SIE velocity dispersions, ellipticities, and position angles of the lens systems, as well as time delay predictions assuming a concordance cosmological model.

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.

The virial masses of ultra-diffuse galaxies (UDGs) have been estimated using the kinematics and abundance of their globular cluster populations, leading to disparate results. Some studies conclude that UDGs reside in massive dark matter halos while others, controversially, argue for the existence of UDGs with no dark matter at all. Here we show that these results arise because the uncertainties of these mass estimates have been substantially underestimated. Indeed, applying the same procedure to the well-studied Fornax dwarf spheroidal would conclude that it has an "overmassive" dark halo or, alternatively, that it lacks dark matter. We corroborate our argument with self-consistent mocks of tracers in cosmological halos, showing that masses from samples with $5 < N < 10$ tracers (assuming no measurement errors) are uncertain by at least an order of magnitude. Finally, we estimate masses of UDGs with HST imaging in Coma and show that their recent mass measurements (with adequate uncertainties) are in agreement with that of other dwarfs, such as Fornax.. We also provide bias and scatter factors for a range of sample sizes and measurement errors, of wider applicability.

We report the discovery of a candidate to quadrupole gravitationally lensed system KiDS0239-3211 based on the public data release 3 of the KiDS survey and machine learning techniques.

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.

New methods have been recently developed to search for strong gravitational lenses, in particular lensed quasars, in wide-field imaging surveys. Here, we compare the performance of three different, morphology- and photometry- based methods to find lens candidates over the Kilo-Degree Survey (KiDS) DR3 footprint (440 deg$^2$). The three methods are: i) a multiplet detection in KiDS-DR3 and/or Gaia-DR1, ii) direct modeling of KiDS cutouts and iii) positional offsets between different surveys (KiDS-vs-Gaia, Gaia-vs-2MASS), with purpose-built astrometric recalibrations. The first benchmark for the methods has been set by the recovery of known lenses. We are able to recover seven out of ten known lenses and pairs of quasars observed in the KiDS DR3 footprint, or eight out of ten with improved selection criteria and looser colour pre-selection. This success rate reflects the combination of all methods together, which, taken individually, performed significantly worse (four lenses each). One movelty of our analysis is that the comparison of the performances of the different methods has revealed the pros and cons of the approaches and, most of all, the complementarities. We finally provide a list of high-grade candidates found by one or more methods, awaiting spectroscopic follow-up for confirmation. Of these, KiDS 1042+0023 is to our knowledge the first confirmed lensed quasar from KiDS, exhibiting two quasar spectra at the same source redshift at either sides of a red galaxy, with uniform flux-ratio $f\approx1.25$ over the wavelength range $0.45\mu\mathrm{m}<\lambda<0.75\mu\mathrm{m}.$

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 have scanned 5000 deg$^{2}$ of Southern Sky to search for strongly lensed quasars with five methods, all source-oriented, but based on different assumptions and selection criteria. We analyse morphological searches based on Gaia multiplet detection and chromatic offsets, fibre-spectroscopic preselection, and X-ray and radio preselection. The performance and complementarity of the methods are evaluated on a common sample of known lenses in the Dark Energy Survey public DR1 footprint. We recovered in total 13 known lenses, of which 8 quadruplets. The method that found the largest number of known lenses is the one based on morphological and colour selection of objects from the WISE and Gaia-DR2 Surveys. We finally present a list of high-grade candidates from each method, to facilitate follow-up spectroscopic campaigns, including two previously unknown quadruplets: WG210014.9-445206.4 and WG021416.37-210535.3.

The Southern Hemisphere has just recently begun to be charted by wide-field surveys, with a sufficient depth and image quality to enable the discovery of strongly lensed quasars. The quadruply imaged quasar WG0214-2105 (r.a.=02:14:16.37, dec.=-21:05:35.3) is a previously unknown lens, with `blue' mid-IR colors and high UV deficit, found in the intersection of three survey footprints: the Dark Energy Survey public DR1 (DES, Abbott et al. 2018), The VST-ATLAS (Shanks et al. 2015) and Pan-STARRS (Chambers et al. 2016). Its discovery relied on high spatial resolution from the Gaia mission (Lindegren et al. 2016) and mid-IR color preselection in the WISE catalog (Wright et al. 2010).

M49 (NGC 4472) is the dominant galaxy in subcluster B of the Virgo Cluster, and a benchmark for studying the build-up of the extended halos of brightest group galaxies in the outskirts of galaxy clusters. We investigate the kinematics in the outer halo of M49, look for substructures, and describe the transition to the surrounding intra-group light. As kinematic tracers we use planetary nebulae (PNe), combining kinematics from the extended Planetary Nebula Spectrograph (PN.S) early-type galaxy survey with our recent deep photometric sample. We study the position-velocity-plane for bright and faint PN populations out to 95 kpc radius, and employ a multi-Gaussian model for the velocity distribution to identify stellar populations with distinct kinematics and histories. We report the detection of stellar-kinematic substructure associated with the interaction of M49 with the dwarf irregular galaxy VCC 1249. We find two kinematically distinct PN populations associated with the main M49 halo and the extended intra-group light (IGL). These have velocity dispersions $\sigma_{halo}\simeq 170$ km/s and $\sigma_{IGL} \simeq400$ km/s at 10-80 kpc radii. The overall luminosity profile and velocity dispersion at $\sim80$ kpc are consistent with a flat circular velocity curve extrapolated from X-ray observations. The dispersion of the PNe associated with the IGL joins onto that of the satellite galaxies in subcluster B at $\sim100$ kpc radius. This is the first time that the transition from halo to IGL is observed based on the velocities of individual stars. Therefore the halo of M49, consisting of at least three distinct components, has undergone an extended accretion history within its parent group potential. The blue colours of the IGL component are consistent with a population of stars formed in low-mass galaxies at redshift $\sim0.5$ that has since evolved passively, as suggested by other data.

Strongly gravitational lensed quasars can be used to measure the so-called time-delay distance $D_{\Delta t}$, and thus the Hubble constant $H_0$ and other cosmological parameters. Stellar kinematics of the deflector galaxy play an essential role in this measurement by: (i) helping break the mass-sheet degeneracy; (ii) determining in principle the angular diameter distance $D_{\rm d}$ to the deflector and thus further improving the cosmological constraints. In this paper we simulate observations of lensed quasars with integral field spectrographs and show that spatially resolved kinematics of the deflector enable further progress by helping break the mass-anisotropy degeneracy. Furthermore, we use our simulations to obtain realistic error estimates with current/upcoming instruments like OSIRIS on Keck and NIRSPEC on the James Webb Space Telescope for both distances (typically $\sim6$ per cent on $D_{\Delta t}$ and $\sim10$ per cent on $D_{\rm d}$). We use the error estimates to compute cosmological forecasts for the sample of nine lenses that currently have well measured time delays and deep Hubble Space Telescope images and for a sample of 40 lenses that is projected to be available in a few years through follow-up of candidates found in ongoing wide field surveys. We find that $H_0$ can be measured with 2 per cent (1 per cent) precision from nine (40) lenses in a flat $\Lambda$cold dark matter cosmology. We study several other cosmological models beyond the flat $\Lambda$cold dark matter model and find that time-delay lenses with spatially resolved kinematics can greatly improve the precision of the cosmological parameters measured by cosmic microwave background data.

We use data from the HST Coma Cluster Treasury program to assess the richness of the Globular Cluster Systems (GCSs) of 54 Coma ultra-diffuse galaxies (UDGs), 18 of which have a half-light radius exceeding 1.5 kpc. We use a hierarchical Bayesian method tested on a large number of mock datasets to account consistently for the high and spatially varying background counts in Coma. These include both background galaxies and intra-cluster GCs (ICGCs), which are disentangled from the population of member GCs in a probabilistic fashion. We find no candidate for a GCS as rich as that of the Milky Way, our sample has GCSs typical of dwarf galaxies. For the standard relation between GCS richness and halo mass, 33 galaxies have a virial mass $M_{vir}\leq10^{11}M_\odot$ at 90% probability. Only three have $M_{vir}>10^{11}M_\odot$ with the same confidence. The mean colour and spread in colour of the UDG GCs are indistinguishable from those of the abundant population of ICGCs. The majority of UDGs in our sample are consistent with the relation between stellar mass and GC richness of 'normal' dwarf galaxies. Nine systems, however, display GCSs that are richer by a factor of 3 or more (at 90% probability). Six of these have sizes $\lesssim1.4$ kpc. Our results imply that the physical mechanisms responsible for the extended size of the UDGs and for the enhanced GC richness of some cluster dwarfs are at most weakly correlated.

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.

Strongly lensed active galactic nuclei (AGN) provide a unique opportunity to make progress in the study of the evolution of the correlation between the mass of supermassive black holes ($\mathcal M_{BH}$) and their host galaxy luminosity ($L_{host}$). We demonstrate the power of lensing by analyzing two systems for which state-of-the-art lens modelling techniques have been applied to Hubble Space Telescope imaging data. We use i) the reconstructed images to infer the total and bulge luminosity of the host and ii) published broad-line spectroscopy to estimate $\mathcal M_{BH}$ using the so-called virial method. We then enlarge our sample with new calibration of previously published measurements to study the evolution of the correlation out to z~4.5. Consistent with previous work, we find that without taking into account passive luminosity evolution, the data points lie on the local relation. Once passive luminosity evolution is taken into account, we find that BHs in the more distant Universe reside in less luminous galaxies than today. Fitting this offset as $\mathcal M_{BH}$/$L_{host}$ $\propto$ (1+z)$^{\gamma}$, and taking into account selection effects, we obtain $\gamma$ = 0.6 $\pm$ 0.1 and 0.8$\pm$ 0.1 for the case of $\mathcal M_{BH}$-$L_{bulge}$ and $\mathcal M_{BH}$-$L_{total}$, respectively. To test for systematic uncertainties and selection effects we also consider a reduced sample that is homogeneous in data quality. We find consistent results but with considerably larger uncertainty due to the more limited sample size and redshift coverage ($\gamma$ = 0.7 $\pm$ 0.4 and 0.2$\pm$ 0.5 for $\mathcal M_{BH}$-$L_{bulge}$ and $\mathcal M_{BH}$-$L_{total}$, respectively), highlighting the need to gather more high-quality data for high-redshift lensed quasar hosts. Our result is consistent with a scenario where the growth of the black hole predates that of the host galaxy.

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.

I describe two novel techniques originally devised to select strongly lensed quasar candidates in wide-field surveys. The first relies on outlier selection in optical and mid-infrared magnitude space; the second combines mid-infrared colour selection with GAIA spatial resolution, to identify multiplets of objects with quasar-like colours. Both methods have already been applied successfully to the SDSS, ATLAS and DES footprints: besides recovering known lenses from previous searches, they have led to new discoveries, including quadruply lensed quasars, which are rare within the rare-object class of quasar lenses. As a serendipitous by-product, at least four candidate Galactic streams in the South have been identified among foreground contaminants. There is considerable scope for tailoring the WISE-GAIA multiplet search to stellar-like objects, instead of quasar-like, and to automatically detect Galactic streams.

We report the discovery of the quadruply lensed quasar J1433+6007, mined in the SDSS DR12 photometric catalogues using a novel outlier-selection technique, without prior spectroscopic or UV excess information. Discovery data obtained at the Nordic Optical telescope (NOT, La Palma) show nearly identical quasar spectra at $z_s=2.74$ and four quasar images in a fold configuration, one of which sits on a blue arc. The deflector redshift is $z_{l}=0.407,$ from Keck-ESI spectra. We describe the selection procedure, discovery and follow-up, image positions and $BVRi$ magnitudes, and first results and forecasts from simple lens models.

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 report the discovery and spectroscopic confirmation of the quad-like lensed quasar system DES J0408-5354 found in the Dark Energy Survey (DES) Year 1 (Y1) data. This system was discovered during a search for DES Y1 strong lensing systems using a method that identified candidates as red galaxies with multiple blue neighbors. DES J0408-5354 consists of a central red galaxy surrounded by three bright (i < 20) blue objects and a fourth red object. Subsequent spectroscopic observations using the Gemini South telescope confirmed that the three blue objects are indeed the lensed images of a quasar with redshift z = 2.375, and that the central red object is an early-type lensing galaxy with redshift z = 0.597. DES J0408-5354 is the first quad lensed quasar system to be found in DES and begins to demonstrate the potential of DES to discover and dramatically increase the sample size of these very rare objects.

Strong gravitational lens systems with time delays between the multiple images allow measurements of time-delay distances, which are primarily sensitive to the Hubble constant that is key to probing dark energy, neutrino physics, and the spatial curvature of the Universe, as well as discovering new physics. We present H0LiCOW ($H_0$ Lenses in COSMOGRAIL's Wellspring), a program that aims to measure $H_0$ with $<3.5\%$ uncertainty from five lens systems (B1608+656, RXJ1131-1231, HE0435-1223, WFI2033-4723 and HE1104-1805). We have been acquiring (1) time delays through COSMOGRAIL and Very Large Array monitoring, (2) high-resolution Hubble Space Telescope imaging for the lens mass modeling, (3) wide-field imaging and spectroscopy to characterize the lens environment, and (4) moderate-resolution spectroscopy to obtain the stellar velocity dispersion of the lenses for mass modeling. In cosmological models with one-parameter extension to flat $\Lambda$CDM, we expect to measure $H_0$ to $<3.5\%$ in most models, spatial curvature $\Omega_{\rm k}$ to 0.004, $w$ to 0.14, and the effective number of neutrino species to 0.2 (1$\sigma$ uncertainties) when combined with current CMB experiments. These are, respectively, a factor of $\sim15$, $\sim2$, and $\sim1.5$ tighter than CMB alone. Our data set will further enable us to study the stellar initial mass function of the lens galaxies, and the co-evolution of supermassive black holes and their host galaxies. This program will provide a foundation for extracting cosmological distances from the hundreds of time-delay lenses that are expected to be discovered in current and future surveys.

Wide-field photometric surveys enable searches of rare yet interesting objects, such as strongly lensed quasars or quasars with a bright host galaxy. Past searches for lensed quasars based on their optical and near infrared properties have relied on photometric cuts and spectroscopic pre-selection (as in the Sloan Quasar Lens Search), or neural networks applied to photometric samples. These methods rely on cuts in morphology and colours, with the risk of losing many interesting objects due to scatter in their population properties, restrictive training sets, systematic uncertainties in catalog-based magnitudes, and survey-to-survey photometric variations. Here, we explore the performance of a Gaussian Mixture Model to separate point-like quasars, quasars with an extended host, and strongly lensed quasars using griz psf and model magnitudes and WISE W1, W2. The choice of optical magnitudes is due to their presence in all current and upcoming releases of wide-field surveys, whereas UV information is not always available. We then assess the contamination from blue galaxies and the role of additional features such as W3 magnitudes or psf-model terms as morphological information. As a demonstration, we conduct a search in a random 10% of the SDSS footprint, and we provide the catalog of the 43 SDSS object with the highest `lens' score in our selection that survive visual inspection, and are spectroscopically confirmed to host active nuclei. We inspect archival data and find images of 5/43 objects in the Hubble Legacy Archive, including 2 known lenses. The code and materials are available to facilitate follow-up.

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.

The arrival times, positions, and fluxes of multiple images in strong lens systems can be used to infer the presence of dark subhalos in the deflector, and thus test predictions of cold dark matter models. However, gravitational lensing does not distinguish between perturbations to a smooth gravitational potential arising from baryonic and non-baryonic mass. In this work, we quantify the extent to which the stellar mass distribution of a deflector can reproduce flux ratio and astrometric anomalies typically associated with the presence of a dark matter subhalo. Using Hubble Space Telescope images of nearby galaxies, we simulate strong lens systems with real distributions of stellar mass as they would be observed at redshift $z_d=0.5$. We add a dark matter halo and external shear to account for the smooth dark matter field, omitting dark substructure, and use a Monte Carlo procedure to characterize the distributions of image positions, time delays, and flux ratios for a compact background source of diameter 5 pc. By convolving high-resolution images of real galaxies with a Gaussian PSF, we simulate the most detailed smooth potential one could construct given high quality data, and find scatter in flux ratios of $\approx 10\%$, which we interpret as a typical deviation from a smooth potential caused by large and small scale structure in the lensing galaxy. We demonstrate that the flux ratio anomalies arising from galaxy-scale baryonic structure can be minimized by selecting the most massive and round deflectors, and by simultaneously modeling flux ratio and astrometric data.

We report the observation and confirmation of the first group- and cluster-scale strong gravitational lensing systems found in Dark Energy Survey (DES) data. Through visual inspection of data from the Science Verification (SV) season, we identified 53 candidate systems. We then obtained spectroscopic follow-up of 21 candidates using the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South telescope and the Inamori-Magellan Areal Camera and Spectrograph (IMACS) at the Magellan/Baade telescope. With this follow-up, we confirmed six candidates as gravitational lenses: Three of the systems are newly discovered, and the remaining three were previously known. Of the 21 observed candidates, the remaining 15 were either not detected in spectroscopic observations, were observed and did not exhibit continuum emission (or spectral features), or were ruled out as lensing systems. The confirmed sample consists of one group-scale and five galaxy cluster-scale lenses. The lensed sources range in redshift z ~ 0.80-3.2, and in i-band surface brightness i_{SB} ~ 23-25 mag/sq.-arcsec. (2" aperture). For each of the six systems, we estimate the Einstein radius and the enclosed mass, which have ranges ~ 5.0 - 8.6" and ~ 7.5 x 10^{12} - 6.4 x 10^{13} solar masses, respectively.

Lens time delays are a powerful probe of cosmology, provided that the gravitational potential of the main deflector can be modeled with sufficient precision. Recent work has shown that this can be achieved by detailed modeling of the host galaxies of lensed quasars, which appear as "Einstein Rings" in high resolution images. We carry out a systematic exploration of the high resolution imaging required to exploit the thousands of lensed quasars that will be discovered by current and upcoming surveys with the next decade. Specifically, we simulate realistic lens systems as imaged by the Hubble Space Telescope (HST), James Webb Space Telescope (JWST), and ground based adaptive optics images taken with Keck or the Thirty Meter Telescope (TMT). We compare the performance of these pointed observations with that of images taken by the Euclid (VIS), Wide-Field Infrared Survey Telescope (WFIRST) and Large Synoptic Survey Telescope (LSST) surveys. We use as our metric the precision with which the slope $\gamma'$ of the total mass density profile $\rho_{tot}\propto r^{-\gamma'}$ for the main deflector can be measured. Ideally, we require that the statistical error on $\gamma'$ be less than 0.02, such that it is subdominant to other sources of random and systematic uncertainties. We find that survey data will likely have sufficient depth and resolution to meet the target only for the brighter gravitational lens systems, comparable to those discovered by the SDSS survey. For fainter systems, that will be discovered by current and future surveys, targeted follow-up will be required. However, the exposure time required with upcoming facilitites such as JWST, the Keck Next Generation Adaptive Optics System, and TMT, will only be of order a few minutes per system, thus making the follow-up of hundreds of systems a practical and efficient cosmological probe.

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 spectroscopy and laser guide star adaptive optics (LGSAO) images of the doubly imaged lensed quasar SDSS J1206+4332. We revise the deflector redshift proposed previously to $z_{d}=0.745,$ and measure for the first time its velocity dispersion $\sigma=(290\pm30)$ km/s. The LGSAO data show the lensed quasar host galaxy stretching over the astroid caustic thus forming an extra pair of merging images, which was previously thought to be an unrelated galaxy in seeing limited data. Owing to the peculiar geometry, the lens acts as a natural coronagraph on the broad-line region of the quasar so that only narrow [O III] emission is found in the fold arc. We use the data to reconstruct the source structure and deflector potential, including nearby perturbers. We reconstruct the point-spread function (PSF) from the quasar images themselves, since no additional point source is present in the field of view. From gravitational lensing and stellar dynamics, we find the slope of the total mass density profile to be $\gamma^{\prime}=-\log\rho/\log r =1.93\pm0.09.$ We discuss the potential of SDSS J1206+4332 for measuring time delay distance (and thus H$_0$ and other cosmological parameters), or as a standard ruler, in combination with the time delay published by the COSMOGRAIL collaboration. We conclude that this system is very promising for cosmography. However, in order to achieve competitive precision and accuracy, an independent characterization of the PSF is needed. Spatially resolved kinematics of the deflector would reduce the uncertainties further. Both are within the reach of current observational facilities.

Gravitationally lensed (GL) quasars are brighter than their unlensed counterparts and produce images with distinctive morphological signatures. Past searches and target selection algorithms, in particular the Sloan Quasar Lens Search (SQLS), have relied on basic morphological criteria, which were applied to samples of bright, spectroscopically confirmed quasars. The SQLS techniques are not sufficient for searching into new surveys (e.g. DES, PS1, LSST), because spectroscopic information is not readily available and the large data volume requires higher purity in target/candidate selection. We carry out a systematic exploration of machine learning techniques and demonstrate that a two step strategy can be highly effective. In the first step we use catalog-level information ($griz$+WISE magnitudes, second moments) to preselect targets, using artificial neural networks. The accepted targets are then inspected with pixel-by-pixel pattern recognition algorithms (Gradient-Boosted Trees), to form a final set of candidates. The results from this procedure can be used to further refine the simpler SQLS algorithms, with a twofold (or threefold) gain in purity and the same (or $80\%$) completeness at target-selection stage, or a purity of $70\%$ and a completeness of $60\%$ after the candidate-selection step. Simpler photometric searches in $griz$+WISE based on colour cuts would provide samples with $7\%$ purity or less. Our technique is extremely fast, as a list of candidates can be obtained from a stage III experiment (e.g. DES catalog/database) in {a few} CPU hours. The techniques are easily extendable to Stage IV experiments like LSST with the addition of time domain information.

We study dynamical models for elliptical galaxies, deriving the projected kinematic profiles in a form that is valid for general surface-brightness profiles and (spherical) total mass profiles, without the need for any explicit deprojection. We then show that an almost flat rotation curve, combined with modest velocity anisotropy, is already sufficient to recover the kinematic profiles of nearby ellipticals. As an application, we provide two different sets of mass estimators for elliptical galaxies, based on either the velocity dispersion at a specific location near the effective radius, or the aperture-averaged velocity dispersion. In the large aperture (virial) limit, mass estimators are naturally independent of anisotropy. The spherical mass enclosed within the effective radius $R_{\rm e}$ can be estimated as $2.4 R_{\rm e} \langle \sigma^{2}_{\rm p} \rangle/ G$, where $\langle \sigma^2_{\rm p} \rangle$ is the average of the squared velocity dispersion over a finite aperture. This formula does not depend on assumptions such as mass-follows-light, and is a compromise between the cases of small and large apertures sizes. Its general agreement with results from other methods in the literature makes it a reliable means to infer masses in the absence of detailed kinematic information. If on the other hand the velocity dispersion profile is available, tight mass estimates can be found that are independent of the mass-model and anisotropy profile (within $\approx$ 10% accuracy). Explicit formulae are given for small anisotropy, large radii and/or power-law total densities. Motivated by recent observational claims, we also discuss the issue of weak homology of elliptical galaxies, emphasizing the interplay between morphology and orbital structure.

We study the Globular Cluster (GC) system of the nearby elliptical galaxy M87 using the newly available dataset with accurate kinematics (Strader et al.2011). We find evidence for three distinct sub-populations of GCs in terms of colours, kinematics and radial profiles. We show that a decomposition into three populations (blue, intermediate and red GCs) is statistically preferred to one with two or four populations, and relate them to the stellar profile. We exploit the sub-populations to derive dynamical constraints on the mass and Dark Matter (DM) content of M87 out to $\sim100$ kpc. We use a class of global mass-estimators (from Paper I), obtaining mass measurements at different locations. M87's DM fraction changes from $\approx$0.2 at the starlight's effective radius (6 kpc) to $\approx$0.95 at the distance probed by the most extended, blue GCs (135 kpc). We supplement this with \textit{virial decompositions}, exploiting the dynamical model to produce a separation into multiple components. These yield the luminous mass as $5.5^{+1.5}_{-2.0}\times 10^{11}M_\odot$ and the DM within 135 kpc as $8.0^{+1.0}_{-4.0}\times 10^{12}M_\odot.$ The inner DM density behaves as $\rho \sim r^{-\gamma}$ with $\gamma\approx 1.6$. This is steeper than the cosmologically preferred cusp $\rho \sim r^{-1},$ providing evidence of DM contraction. Finally, we combine the GC separation into three sub-populations and the Jeans equations, obtaining information on the orbits of the GC system. The centrally concentrated red GCs exhibit tangential anisotropy, consistent with radial-orbit depletion by tidal shredding. The most extended blue GCs have an isotropic velocity dispersion tensor in the central parts, which becomes more tangential moving outwards, consistent with adiabatic contraction of the DM halo.

We exploit the detection of three distinct stellar subpopulations in the red giant branch of the Fornax dwarf Spheroidal to probe its density distribution. This allows us to resolve directly the evolution with radius of the dark matter mass profile. We find that a cored dark matter halo provides a perfect fit to the data, being consistent with all three stellar populations well within 1-sigma, and for the first time we are able to put constraints on the core size of such a halo. With respect to previous work, we do not strengthen the statistical exclusion of a dark matter cusp in Fornax, but we find that Navarro-Frenk-White haloes would be required to have unrealistically large scale radii in order to be compatible with the data, hence low values of the concentration parameter. We are then forced to conclude that the Fornax dwarf Spheroidal sits within a dark matter halo having a constant density core, with a core size of between 0.6 and 1.8 kpc.

We present a ready-to-use method to constrain the density distribution in early-type galaxy lenses. Assuming a power-law density profile, then joint use of the virial theorem and the lens equation yields simple formulae for the power-law index (or logarithmic density gradient). Any dependence on orbital anisotropy can be tightly constrained or even erased completely. Our results rely just on surface brightnesses and line-of-sight kinematics, making deprojection unnecessary. We revisit three systems that have already been examined in the literature (the Cosmic Horseshoe, the Jackpot and B1608+656) and provide our estimates. Finally, we show that the method yields a good approximation for the density profile even when the true profile is a broken power-law, albeit with a mild bias towards isothermality.

The projected virial theorem is applied to the case of multiple stellar populations in the nearby dwarf spheroidal galaxies. As each population must reside in the same gravitational potential, this provides strong constraints on the nature of the dark matter halo. We derive necessary conditions for two populations with Plummer or exponential surface brightnesses to reside in a cusped Navarro-Frenk-White (NFW) halo. We apply our methods to the Sculptor dwarf spheroidal, and show that there is no NFW halo compatible with the energetics of the two populations. The dark halo must possess a core radius of ~ 120 pc for the virial solutions for the two populations to be consistent. This conclusion remains true, even if the effects of flattening or self-gravity of the stellar populations are included.

The virial theorem prescribes the ratio of the globally-averaged equatorial to vertical velocity dispersion of a tracer population in spherical and flattened dark haloes. This gives sequences of physical models in the plane of global anisotropy and flattening. The tracer may have any density, though there are particularly simple results for power-laws and exponentials. We prove the flattening theorem: for a spheroidally stratified tracer density with axis ratio q in a dark density potential with axis ratio g, the ratio of globally averaged equatorial to vertical velocity dispersion depends only on q/g. As the stellar halo density and velocity dispersion of the Milky Way are accessible to observations, this provides a new method for measuring the flattening of the dark matter. If the kinematics of the local halo subdwarfs are representative, then the Milky Way's dark halo is oblate with a flattening in the potential of g ~ 0.85, corresponding to a flattening in the dark matter density of ~ 0.7. The fractional pressure excess for power-law populations is roughly proportional to both the ellipticity and the fall-off exponent. Given the same pressure excess, if the density profile of one stellar population declines more quickly than that of another, then it must be rounder. This implies that the dual halo structure claimed by Carollo et al. (2007) for the Galaxy, a flatter inner halo and a rounder outer halo, is inconsistent with the virial theorem. For the thick disc, we provide formulae for the virial sequences of double-exponential discs in logarithmic and Navarro-Frenk-White (NFW) haloes. There are good matches to the observational data on the flattening and anisotropy of the thick disc if the thin disc is exponential with a short scalelength ~ 2.6 kpc and normalisation of 56 solar masses per square parsec, together with a logarithmic dark halo.