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.

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.

Relic galaxies are thought to be the progenitors of high-redshift red nuggets that for some reason missed the channels of size growth and evolved passively and undisturbed since the first star formation burst (at $z>2$). These local ultracompact old galaxies are unique laboratories for studying the star formation processes at high redshift and thus the early stage of galaxy formation scenarios. Counterintuitively, theoretical and observational studies indicate that relics are more common in denser environments, where merging events predominate. To verify this scenario, we compared the number counts of a sample of ultracompact massive galaxies (UCMGs) selected within the third data release of the Kilo Degree Survey, that is, systems with sizes $R_{\rm e} < 1.5 \, \rm kpc$ and stellar masses $M_{\rm \star} > 8 \times 10^{10}\, \rm M_{\odot}$, with the number counts of galaxies with the same masses but normal sizes in field and cluster environments. Based on their optical and near-infrared colors, these UCMGs are likely to be mainly old, and hence representative of the relic population. We find that both UCMGs and normal-size galaxies are more abundant in clusters and their relative fraction depends only mildly on the global environment, with denser environments penalizing the survival of relics. Hence, UCMGs (and likely relics overall) are not special because of the environment effect on their nurture, but rather they are just a product of the stochasticity of the merging processes regardless of the global environment in which they live.

We report new high-quality galaxy scale strong lens candidates found in the Kilo Degree Survey data release 4 using Machine Learning. We have developed a new Convolutional Neural Network (CNN) classifier to search for gravitational arcs, following the prescription by \cite{2019MNRAS.484.3879P} and using only $r-$band images. We have applied the CNN to two "predictive samples": a Luminous red galaxy (LRG) and a "bright galaxy" (BG) sample ($r<21$). We have found 286 new high probability candidates, 133 from the LRG sample and 153 from the BG sample. We have then ranked these candidates based on a value that combines the CNN likelihood to be a lens and the human score resulting from visual inspection (P-value) and we present here the highest 82 ranked candidates with P-values $\ge 0.5$. All these high-quality candidates have obvious arc or point-like features around the central red defector. Moreover, we define the best 26 objects, all with scores P-values $\ge 0.7$ as a "golden sample" of candidates. This sample is expected to contain very few false positives and thus it is suitable for follow-up observations. The new lens candidates come partially from the the more extended footprint adopted here with respect to the previous analyses, partially from a larger predictive sample (also including the BG sample). These results show that machine learning tools are very promising to find strong lenses in large surveys and more candidates that can be found by enlarging the predictive samples beyond the standard assumption of LRGs. In the future, we plan to apply our CNN to the data from next-generation surveys such as the Large Synoptic Survey Telescope, Euclid, and the Chinese Space Station Optical Survey.

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.

Ultra-compact massive galaxies UCMGs, i.e. galaxies with stellar masses $M_{*} > 8 \times 10^{10} M_{\odot}$ and effective radii $R_{e} < 1.5$ kpc, are very rare systems, in particular at low and intermediate redshifts. Their origin as well as their number density across cosmic time are still under scrutiny, especially because of the paucity of spectroscopically confirmed samples. We have started a systematic census of UCMG candidates within the ESO Kilo Degree Survey, together with a large spectroscopic follow-up campaign to build the largest possible sample of confirmed UCMGs. This is the third paper of the series and the second based on the spectroscopic follow-up program. Here, we present photometrical and structural parameters of 33 new candidates at redshifts $0.15 \lesssim z \lesssim 0.5$ and confirm 19 of them as UCMGs, based on their nominal spectroscopically inferred $M_{*}$ and $R_{e}$. This corresponds to a success rate of $\sim 58\%$, nicely consistent with our previous findings. The addition of these 19 newly confirmed objects, allows us to fully assess the systematics on the system selection, and finally reduce the number density uncertainties. Moreover, putting together the results from our current and past observational campaigns and some literature data, we build the largest sample of UCMGs ever collected, comprising 92 spectroscopically confirmed objects at $0.1 \lesssim z \lesssim 0.5$. This number raises to 116, allowing for a $3\sigma$ tolerance on the $M_{*}$ and $R_{e}$ thresholds for the UCMG definition. For all these galaxies we have estimated the velocity dispersion values at the effective radii which have been used to derive a preliminary mass-velocity dispersion correlation.

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.

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.

We present a catalog of quasars selected from broad-band photometric ugri data of the Kilo-Degree Survey Data Release 3 (KiDS DR3). The QSOs are identified by the random forest (RF) supervised machine learning model, trained on SDSS DR14 spectroscopic data. We first cleaned the input KiDS data from entries with excessively noisy, missing or otherwise problematic measurements. Applying a feature importance analysis, we then tune the algorithm and identify in the KiDS multiband catalog the 17 most useful features for the classification, namely magnitudes, colors, magnitude ratios, and the stellarity index. We used the t-SNE algorithm to map the multi-dimensional photometric data onto 2D planes and compare the coverage of the training and inference sets. We limited the inference set to r<22 to avoid extrapolation beyond the feature space covered by training, as the SDSS spectroscopic sample is considerably shallower than KiDS. This gives 3.4 million objects in the final inference sample, from which the random forest identified 190,000 quasar candidates. Accuracy of 97%, purity of 91%, and completeness of 87%, as derived from a test set extracted from SDSS and not used in the training, are confirmed by comparison with external spectroscopic and photometric QSO catalogs overlapping with the KiDS footprint. The robustness of our results is strengthened by number counts of the quasar candidates in the r band, as well as by their mid-infrared colors available from WISE. An analysis of parallaxes and proper motions of our QSO candidates found also in Gaia DR2 suggests that a probability cut of p(QSO)>0.8 is optimal for purity, whereas p(QSO)>0.7 is preferable for better completeness. Our study presents the first comprehensive quasar selection from deep high-quality KiDS data and will serve as the basis for versatile studies of the QSO population detected by this survey.

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.

We present a new sample of galaxy-scale strong gravitational-lens candidates, selected from 904 square degrees of Data Release 4 of the Kilo-Degree Survey (KiDS), i.e., the "Lenses in the Kilo-Degree Survey" (LinKS) sample. We apply two Convolutional Neural Networks (ConvNets) to $\sim88\,000$ colour-magnitude selected luminous red galaxies yielding a list of 3500 strong-lens candidates. This list is further down-selected via human inspection. The resulting LinKS sample is composed of 1983 rank-ordered targets classified as "potential lens candidates" by at least one inspector. Of these, a high-grade subsample of 89 targets is identified with potential strong lenses by all inspectors. Additionally, we present a collection of another 200 strong lens candidates discovered serendipitously from various previous ConvNet runs. A straightforward application of our procedure to future Euclid or LSST data can select a sample of $\sim3000$ lens candidates with less than 10 per cent expected false positives and requiring minimal human intervention.

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.

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.

We present results from an ongoing investigation using the Kilo Degree Survey (KiDS) on the VLT Survey Telescope (VST) to provide a census of ultra-compact massive galaxies (UCMGs), defined as galaxies with stellar masses $M_{\rm \star} > 8 \times 10^{10} \rm M_{\odot}$ and effective radii $R_{\rm e} < 1.5\,\rm kpc$. UCMGs, which are expected to have undergone very few merger events, provide a unique view on the accretion history of the most massive galaxies in the Universe. Over an effective sky area of nearly 330 square degrees, we select UCMG candidates from KiDS multi-colour images, which provide high quality structural parameters, photometric redshifts and stellar masses. Our sample of $\sim 1000$ photometrically selected UCMGs at $z < 0.5$ represents the largest sample of UCMG candidates assembled to date over the largest sky area. In this paper we present the first effort to obtain their redshifts using different facilities, starting with first results for 28 candidates with redshifts $z < 0.5$, obtained at NTT and TNG telescopes. We confirmed, as bona fide UCMGs, 19 out of the 28 candidates with new redshifts. A further 46 UCMG candidates are confirmed with literature spectroscopic redshifts (35 at $z < 0.5$), bringing the final cumulative sample of spectroscopically-confirmed lower-z UCMGs to 54 galaxies, which is the largest sample at redshifts below $0.5$. We use these spectroscopic redshifts to quantify systematic errors in our photometric selection, and use these to correct our UCMG number counts. We finally compare the results to independent datasets and simulations.

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

We present the largest and most spatially extended planetary nebulae (PNe) catalog ever obtained for the Fornax cluster. We measured velocities of 1452 PNe out to 200 kpc in the cluster core using a counter-dispersed slitless spectroscopic technique with data from FORS2 on the VLT. With such extended spatial coverage, we can study separately the stellar halos of some of the cluster main galaxies and the intracluster light. In this second paper of the Fornax Cluster VLT Spectroscopic Survey (FVSS), we identify and classify the emission-line sources, describe the method to select PNe and calculate their coordinates and velocities from the dispersed slitless images. From the PN 2D velocity map we identify stellar streams that are possibly tracing the gravitational interaction of NGC1399 with NGC1404 and NGC1387. We also present the velocity dispersion profile out to $\sim 200$ kpc radii, which shows signatures of a superposition of the bright central galaxy and the cluster potential, with the latter clearly dominating the regions outside R $\sim 1000$" ($\sim 100$ kpc).

We present the results of a wide spectroscopic survey aimed at detecting extragalactic globular clusters (GCs) in the core of the Fornax cluster. About 4500 low resolution spectra (from 4800 to 10000 \AA ) were observed in 25 VLT/VIMOS masks covering the central 1 deg$^{2}$ around the dominant galaxy NGC 1399 corresponding to $\sim$175 kpc galactocentric radius. We describe the methodology used for data reduction and data analysis. We found a total of 387 unique physical objects (372 GCs and 15 ultra compact dwarfs) in the field covered by our observations. Most of these objects are associated with NGC 1399, with only 10% likely belonging to other giant galaxies. The new VIMOS dataset is complementary to the many GC catalogues already present in the literature and it brings the total number of tracer particles around NGC 1399 to more than 1130 objects. With this comprehensive radial velocity sample we have found that the velocity dispersion of the GC population (equally for red and blue GC populations) shows a relatively sharp increase from low velocity dispersion ($\sim250$-$350$ kms$^{-1}$) to high velocity dispersion ($\sim300$-$400$ kms$^{-1}$) at projected radius of $\approx10$ arcmin ($\sim60$ kpc) from the galaxy centre. This suggests that at a projected radius of $\approx60$ kpc both blue and red GC populations begin to be governed by the dominating Fornax cluster potential, rather than by the central NGC 1399 galaxy potential. This kinematic evidence corroborates similar results found using surface brightness analysis and planetary nebulae kinematics.

Based on MUSE integral-field data we present evidence for a radial variation at the low-mass end of the stellar initial-mass function (IMF) in the central regions of the giant early-type galaxy NGC4486 (M87). We used state-of-the-art stellar population models and the observed strength of various IMF-sensitive absorption-line features to solve for the best low-mass tapered "bimodal" form of the IMF, while accounting also for radial variations in stellar metallicity, the overall $\alpha$-elements abundance, and the abundance of individual elements such as Ti, O, Na and Ca. Our analysis reveals a strong IMF gradient in M87, corresponding to an exceeding fraction of low-mass stars compared to the case of the Milky Way toward the center of M87 that drops to nearly Milky-way levels by 0.4 $R_e$. This IMF gradient is found to correlate well with both the radial profile for stellar metallicity and for $\alpha$-elements abundance but not with stellar velocity dispersion. Such IMF variations correspond to over a factor two increase in stellar mass-to-light M/L ratio compared to the case of a Milky-way like IMF, consistent with other investigations into IMF gradients in early-type galaxies, including recent dynamical constraints on M/L radial variations in M87 by Oldham & Auger. In addition to constraining the IMF in M87 we also looked into the abundance of Sodium, which turned up to be super-Solar over the entire radial range of our MUSE observations and to exhibit a considerable negative gradient. These findings suggest an additional role of metallicity in boosting the Na-yields in the central, metal-rich regions of M87 during its early and brief star-formation history. Our work adds the case of M87 to the few objects that as of today have radial constraints on their IMF or [Na/Fe] abundance, while also illustrating the accuracy that MUSE could bring to this kind of investigations.

We present a new measurement of the Hubble Constant H0 and other cosmological parameters based on the joint analysis of three multiply-imaged quasar systems with measured gravitational time delays. First, we measure the time delay of HE0435-1223 from 13-year light curves obtained as part of the COSMOGRAIL project. Companion papers detail the modeling of the main deflectors and line of sight effects, and how these data are combined to determine the time-delay distance of HE 0435-1223. Crucially, the measurements are carried out blindly with respect to cosmological parameters in order to avoid confirmation bias. We then combine the time-delay distance of HE0435-1223 with previous measurements from systems B1608+656 and RXJ1131-1231 to create a Time Delay Strong Lensing probe (TDSL). In flat $\Lambda$CDM with free matter and energy density, we find $H_0$ = 71.9 +2.4 -3.0 km/s/Mpc and $\Omega_{\Lambda}$ = 0.62 +0.24 -0.35 . This measurement is completely independent of, and in agreement with, the local distance ladder measurements of H0. We explore more general cosmological models combining TDSL with other probes, illustrating its power to break degeneracies inherent to other methods. The TDSL and Planck joint constraints are $H_0$ = 69.2 +1.4 -2.2 km/s/Mpc, $\Omega_{\Lambda}$ = 0.70 +0.01 -0.01 and $\Omega_k$ = 0.003 +0.004 -0.006 in open $\Lambda$CDM and $H_0$ = 79.0 +4.4 -4.2 km/s/Mpc, $\Omega_{de}$ = 0.77 +0.02 -0.03 and $w$ = -1.38 +0.14 -0.16 in flat $w$CDM. Combined with Planck and Baryon Acoustic Oscillation data, when relaxing the constraints on the numbers of relativistic species we find $N_{eff}$ = 3.34 +0.21 -0.21 and when relaxing the total mass of neutrinos we find 0.182 eV. In an open $w$CDM in combination with Planck and CMB lensing we find $H_0$ = 77.9 +5.0 -4.2 km/s/Mpc, $\Omega_{de}$ = 0.77 +0.03 -0.03, $\Omega_k$ = -0.003 +0.004 -0.004 and $w$ = -1.37 +0.18 -0.23.

Galaxies located in the environment or on the line of sight towards gravitational lenses can significantly affect lensing observables, and can lead to systematic errors on the measurement of $H_0$ from the time-delay technique. We present the results of a systematic spectroscopic identification of the galaxies in the field of view of the lensed quasar HE0435-1223, using the W. M. Keck, Gemini and ESO-Very Large telescopes. Our new catalog triples the number of known galaxy redshifts in the vicinity of the lens, expanding to 102 the number of measured redshifts for galaxies separated by less than 3 arcmin from the lens. We complement our catalog with literature data to gather redshifts up to 15 arcmin from the lens, and search for galaxy groups or clusters projected towards HE0435-1223. We confirm that the lens is a member of a small group that includes at least 12 galaxies, and find 8 other group candidates near the line of sight of the lens. The flexion shift, namely the shift of lensed images produced by high order perturbation of the lens potential, is calculated for each galaxy/group and used to identify which objects produce the largest perturbation of the lens potential. This analysis demonstrates that i) at most three of the five brightest galaxies projected within 12 arcsec of the lens need to be explicitly used in the lens models, and ii) the groups can be treated in the lens model as an external tidal field (shear) contribution.

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.

We present the X-shooter Lens Survey (XLENS) data. The main goal of XLENS is to disentangle the stellar and dark matter content of massive early-type galaxies (ETGs), through combined strong gravitational lensing, dynamics and spectroscopic stellar population studies. The sample consists of 11 lens galaxies covering the redshift range from $0.1$ to $0.45$ and having stellar velocity dispersions between $250$ and $380\,\mathrm{km}\,\mathrm{s}^{-1}$. All galaxies have multi-band, high-quality HST imaging. We have obtained long-slit spectra of the lens galaxies with X-shooter on the VLT. We are able to disentangle the dark and luminous mass components by combining lensing and extended kinematics data-sets, and we are also able to precisely constrain stellar mass-to-light ratios and infer the value of the low-mass cut-off of the IMF, by adding spectroscopic stellar population information. Our goal is to correlate these IMF parameters with ETG masses and investigate the relation between baryonic and non-baryonic matter during the mass assembly and structure formation processes. In this paper we provide an overview of the survey, highlighting its scientific motivations, main goals and techniques. We present the current sample, briefly describing the data reduction and analysis process, and we present the first results on spatially resolved kinematics.

We present an integral field study of the internal structure, kinematics and stellar population of the almost edge-on, intermediate luminosity ($L_ {*}$) elliptical galaxy NGC 4697. We build extended 2-dimensional (2D) maps of the stellar kinematics and line-strengths of the galaxy up to $\sim 0.7 $ effective radii (R$_{eff}$) using a mosaic of 8 VIMOS (VIsible Multi-Objects Spectrograph on the VLT) integral-field unit pointings. We find clear evidence for a rotation-supported structure along the major axis from the 2D kinematical maps, confirming the previous classification of this system as a `fast-rotator'. We study the correlations between the third and fourth Gauss-Hermite moments of the line-of-sight velocity distribution (LOSVD) $h_3$ and $h_4$ with the rotation parameter ($V/\sigma$), and compare our findings to hydrodynamical simulations. We find remarkable similarities to predictions from gas-rich mergers. Based on photometry, we perform a bulge/disk decomposition and study the stellar population properties of the two components. The bulge and the disk show different stellar populations, with the stars in the bulge being older (age$_{\rm bulge}=13.5^{+1.4}_{-1.4}$ Gyr, age$_{\rm disk}=10.5^{+1.6}_{-2.0}$Gyr) and more metal-poor ($\mathrm{[M/H]_{bulge}} = -0.17^{+0.12}_{-0.1}$, $\mathrm{[M/H]_{disk}}=-0.03^{+0.02}_{-0.1}$). The evidence of a later-formed, more metal-rich disk embedded in an older, more metal-poor bulge, together with the LOSVD structure, supports a mass assembly scenario dominated by gas-rich minor mergers and possibly with a late gas-rich major merger that left a previously rapidly rotating system unchanged. The bulge and the disk do not show signs of different stellar Initial Mass Function slopes, and both match well with a Milky Way-like IMF.

We conduct a detailed lensing, dynamics and stellar population analysis of nine massive lens early-type galaxies (ETGs) from the X-Shooter Lens Survey (XLENS). Combining gravitational lensing constraints from HST imaging with spatially-resolved kinematics and line-indices constraints from VLT X-Shooter (XSH) spectra, we infer the low-mass slope and the low cut-off mass of the stellar Initial Mass Function (IMF): $x_{250}=2.37^{+0.12}_{-0.12}$ and $M_{{\rm low}, 250}= 0.131^{+0.023}_{-0.026}\, M_{\odot}$, respectively, for a reference point with $\sigma \equiv 250\, {{\rm kms}}^{-1}$ and R$_{{\rm eff}} \equiv 10$ kpc. All the XLENS systems are consistent with an IMF slope steeper than Milky Way-like. We find no significant correlations between IMF slope and any other quantity, except for an anti-correlation between total dynamical mass density and low-mass IMF slope at the 87% CL [$dx/d\log(\rho)$ = $ -0.19^{+0.15}_{-0.15}$]. This anti-correlation is consistent with the low redshift lenses found by Smith et al. (2015) that have high velocity dispersions and high stellar mass densities but surprisingly shallow IMF slopes.

We perform a direct comparison of two state-of-the art single stellar population (SSP) models that have been used to demonstrate the non-universality of the low-mass end of the Initial Mass Function (IMF) slope. The two public versions of the SSP models are restricted to either solar abundance patterns or solar metallicity, too restrictive if one aims to disentangle elemental enhancements, metallicity changes and IMF variations in massive early-type galaxies (ETGs) with star formation histories different from the solar neighborhood. We define response functions (to metallicity and \alpha-abundance) to extend the parameter space of each set of models. We compare these extended models with a sample of Sloan Digital Sky Survey (SDSS) ETGs spectra with varying velocity dispersions. We measure equivalent widths of optical IMF-sensitive stellar features to examine the effect of the underlying model assumptions and ingredients, such as stellar libraries or isochrones, on the inference of the IMF slope down to ~0.1 solar masses. We demonstrate that the steepening of the low-mass end of the Initial Mass Function (IMF) based on a non-degenerate set of spectroscopic optical indicators is robust against the choice of the stellar population model. Although the models agree in a relative sense (i.e. both imply more bottom-heavy IMFs for more massive systems), we find non-negligible differences on the absolute values of the IMF slope inferred at each velocity dispersion by using the two different models. In particular, we find large inconsistency in the quantitative predictions of IMF slope variations and abundance patterns when sodium lines are used. We investigate the possible reasons for these inconsistencies.

The combination of strong gravitational lensing and stellar kinematics provides a powerful and robust method to investigate the mass and dynamical structure of early-type galaxies. We demonstrate this approach by analysing two massive ellipticals from the XLENS Survey for which both high-resolution HST imaging and X-Shooter spectroscopic observations are available. We adopt a flexible axisymmetric two-component mass model for the lens galaxies, consisting of a generalised NFW dark halo and a realistic self-gravitating stellar mass distribution. For both systems, we put constraints on the dark halo inner structure and flattening, and we find that they are dominated by the luminous component within one effective radius. By comparing the tight inferences on the stellar mass from the combined lensing and dynamics analysis with the values obtained from stellar population studies, we conclude that both galaxies are characterised by a Salpeter-like stellar initial mass function.

Under the assumption of a flat Lambda-CDM cosmology, recent data from the Planck satellite point toward a Hubble constant that is in tension with that measured by gravitational lens time delays and by the local distance ladder. Prosaically, this difference could arise from unknown systematic uncertainties in some of the measurements. More interestingly -- if systematics were ruled out -- resolving the tension would require a departure from the flat Lambda-CDM cosmology, introducing for example a modest amount of spatial curvature, or a non-trivial dark energy equation of state. To begin to address these issues, we present here an analysis of the gravitational lens RXJ1131-1231 that is improved in one particular regard: we examine the issue of systematic error introduced by an assumed lens model density profile. We use more flexible gravitational lens models with baryonic and dark matter components, and find that the exquisite Hubble Space Telescope image with thousands of intensity pixels in the Einstein ring and the stellar velocity dispersion of the lens contain sufficient information to constrain these more flexible models. The total uncertainty on the time-delay distance is 6.6% for a single system. We proceed to combine our improved time-delay distance measurements with the WMAP9 and Planck posteriors. In an open Lambda-CDM model, the data for RXJ1131-1231 in combination with Planck favor a flat universe with Omega_k = 0.00+0.01/-0.02 (68% CI). In a flat wCDM model, the combination of RXJ1131-1231 and Planck yields w = -1.52+0.19/-0.20 (68% CI).

We investigate the optical spectral region of spectra of 1000 stars searching for IMF-sensitive features to constrain the low-mass end of the initial mass function (IMF) slope in elliptical galaxies. We use the MILES stellar library in the wavelength range [3500-7500] {\AA} to select indices that are sensitive to cool dwarf stars and that only weakly depend on age and metallicity. In this wavelength range, the response of a change in the effective temperature of the cool red giant population is similar to the response of a change in the number of dwarf stars in the galaxy. We investigate the degeneracy between IMF variation and Teff,RGB and show that it is possible to break it with the new IMF indicators defined here. In particular, we define a CaH1 index around {\lambda}6380 {\AA} that arises purely from cool dwarfs, does not strongly depend on age and is anti-correlated with [a/Fe]. This index allows the determination of the the IMF slope when combined with TiO lines and age- and metallicity-dependent features. The use of several indicators is crucial to break degeneracies between IMF, age, abundance pattern and Teff of the RGB population. We measure line-index strengths of the optical IMF indicators in the SSP models and compare these with index strengths of the same features in stacked SSDSS early-type galaxy spectra with varying velocity dispersions. Using 10 different indicators, we find a clear trend of a steepening IMF with increasing velocity dispersion described by the equation x=(2.3+-0.1)*log{\sigma}200+(2.13+-0.15), where x is the IMF slope and {\sigma}200 is the central stellar velocity dispersion measured in units of 200km/s. We compute the 'IMF mismatch' parameter as the ratio of stellar mass-to-light ratio predicted from the x-{\sigma}200 relation to that inferred from SSP models assuming a Salpeter IMF and find good agreement with independent published results.

We conduct a detailed investigation of the properties of the stellar initial mass function (IMF) in two massive early-type lens galaxies with velocity dispersions of sigma ~245 km/s and sigma ~325 km/s, for which both HST imaging and X-Shooter spectra are available. We compare the inferences obtained from two fully independent methods: (i) a combined gravitational lensing and stellar dynamics (L&D) analysis of the data sets employing self-consistent axisymmetric models, and (ii) a spectroscopic simple stellar population (SSP) analysis of optical line-strength indices, assuming single power-law IMFs. The results from the two approaches are found to be in agreement within the 1-sigma uncertainties. Both galaxies are consistent with having a Salpeter IMF (power-law slope of x = 2.35), which is strongly favoured over a Chabrier IMF (x = 1.8), with probabilities inferred from the joint analysis of 89% and 99%, respectively. Bottom-heavy IMFs significantly steeper than Salpeter (x >= 3.0) are ruled out with decisive evidence (Bayes factor B > 1000) for both galaxies, as they exceed the total mass derived from the L&D constraints. Our analysis allows, for the first time, the inference of the low-mass cut-off of the IMF (M_low). Combining the joint L&D and SSP analyses of both galaxies, we infer an IMF slope of x = 2.22 +/- 0.14, consistent with Salpeter IMF, and a low-mass limit M_low = 0.13 +/- 0.03 M_sun, just above the hydrogen burning limit.

We measure equivalent widths (EW) - focussing on two unique features (NaI and TiO2) of low-mass stars (<0.3M\odot) - for luminous red galaxy spectra from the the Sloan Digital Sky Survey (SDSS) and X-Shooter Lens Survey (XLENS) in order to study the low-mass end of the initial mass function (IMF). We compare these EWs to those derived from simple stellar population models computed with different IMFs, ages, [{\alpha}/Fe], and elemental abundances. We find that models are able to simultaneously reproduce the observed NaD {\lambda}5895 and Na I {\lambda}8190 features for lower-mass (\sim {\sigma}\ast) early-type galaxies (ETGs) but deviate increasingly for more massive ETGs, due do strongly mismatching NaD EWs. The TiO2 {\lambda}6230 and the Na I {\lambda}8190 features together appear to be a powerful IMF diagnostic, with age and metallicity effects orthogonal to the effect of IMF. We find that both features correlate strongly with galaxy velocity dispersion. The XLENS ETG (SDSSJ0912+0029) and an SDSS ETG (SDSSJ0041-0914) appear to require both an extreme dwarf-rich IMF and a high sodium enhancement ([Na/Fe] = +0.4). In addition, lensing constraints on the total mass of the XLENS system within its Einstein radius limit a bottom-heavy IMF with a power-law slope to x \leq 3.0 at the 90% C.L. We conclude that NaI and TiO features, in comparison with state-of-the-art SSP models, suggest a mildly steepening IMF from Salpeter (dn/dm \propto m-x with x = 2.35) to x \approx 3.0 for ETGs in the range {\sigma} = 200 - 335 km s-1.

We present the first results from the X-shooter Lens Survey (XLENS): an analysis of the massive early-type galaxy SDSS J1148+1930 at redshift z=0.444. We combine its extended kinematic profile -derived from spectra obtained with X-shooter on the ESO VLT- with strong gravitational lensing and multi-color information derived from SDSS images. Our main results are (i) The luminosity-weighted stellar velocity dispersion is \sigma *(<Reff) = 351\pm10 km/s, extracted from a rectangular aperture of 1.8"x1.6", more accurate and considerably lower than a previously published value of \sim450km/s. (ii) A single-component (stellar plus dark) mass model of the lens galaxy yields a logarithmic total-density slope of {\gamma}' = 1.72+0.05-0.06. (iii) The projected stellar mass fraction, derived solely from lensing, is f*(<Reff) = 0.19+0.04-0.09 inside the Einstein radius for a Hernquist profile and no anisotropy. The dark-matter fraction inside the effective radius f_DM(<Reff) = 0.60+0.15\pm0.1. (iv) Based on the SDSS colors, we find f*(<Reff)=0.17 \pm 0.06 for a Salpeter IMF and f*(<Reff)=0.07 \pm 0.02 for a Chabrier IMF. The lensing and dynamics constraints on the stellar mass fraction agree well with those independently derived from the SDSS colors for a Salpeter IMF, preferred over Chabrier IMF at variance with standard results for lower mass galaxies. Dwarf-rich IMFs in the lower mass range of 0.1-0.7 solar mass, with {\alpha} >= 3 (with dN/dM \propto M^(-{\alpha})) -such as that recently suggested for massive ETGs- are excluded at the > 90 % C.L. and in some cases violate the total lensing-derived mass limit. We conclude that this very massive ETG is dark-matter dominated inside one effective radius, consistent with the trend recently found from massive SLACS galaxies, with a IMF normalization consistent with Salpeter.