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.

We present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the MPIA 2.2m telescope at La Silla observatory between October 2016 and February 2020. The lensed quasars HE 0047-1756, WG 0214-2105, DES 0407-5006, 2M 1134-2103, PSJ 1606-2333 and DES 2325-5229 were observed almost daily at high signal-to-noise ratio to obtain high-quality light curves where we can record fast and small-amplitude variations of the quasars. We measure time delays between all pairs of multiple images with only one or two seasons of monitoring at the exception of the time delays relative to image D of PSJ 1606-2333. The most precise estimate is obtained for the delay between image A and image B of DES 0407-5006, with $\tau_{AB} = -128.4^{+3.5}_{-3.8}$ days (2.8% precision) including systematics due to extrinsic variability in the light curves. For HE 0047-1756, we combine our high-cadence data with measurements from decade-long light-curves from previous COSMOGRAIL campaigns, and reach a precision of 0.9 day on the final measurement. The present work demonstrates the feasibility of measuring time delays in lensed quasars in only one or two seasons, provided high signal-to-noise ratio data are obtained at a cadence close to daily.

We present a study of the dependencies of shear and ellipticity bias on simulation (input) and measured (output) parameters, noise, PSF anisotropy, pixel size and the model bias coming from two different and independent shape estimators. We use simulated images from Galsim based on the GREAT3 control-space-constant branch and we measure ellipticity and shear bias from a model-fitting method (gFIT) and a moment-based method (KSB). We show the bias dependencies found on input and output parameters for both methods and we identify the main dependencies and causes. We find consistent results between the two methods (given the precision of the analysis) and important dependencies on orientation and morphology properties such as flux, size and ellipticity. We show cases where shear bias and ellipticity bias behave very different for the two methods due to the different nature of these measurements. We also show that noise and pixelization play an important role on the bias dependences on the output properties. We find a large model bias for galaxies consisting of a bulge and a disk with different ellipticities or orientations. We also see an important coupling between several properties on the bias dependences. Because of this we need to study several properties simultaneously in order to properly understand the nature of shear bias.

In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated datasets presented in paper I. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20\%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ($ |A| < 2\%$) and precision ($< 6\%$ per system), even in the presence of poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix, and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges.

We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data is less sensitive to small-scale structures than the space-based observations. Marginalising over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, $\kappa_{\rm ext}$ at the position of B1608+656 is $\kappa_{\rm ext} = 0.11^{+0.06}_{-0.04}$, where the error bars corresponds respectively to the 16th and 84th quartiles. This result is compatible with previous estimates using the number-counts technique, suggesting that B1608+656 resides in an over-dense line of sight, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favor truncated halo models.

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.

We present the HOLISMOKES project on strong gravitational lensing of supernovae (SNe) as a probe of SN physics and cosmology. We investigate the effects of microlensing on early-phase SN Ia spectra, and find that within 10 rest-frame days after SN explosion, distortions of SN Ia spectra due to microlensing are typically negligible (<1% distortion within the 68% credible region, and ~10% distortion within the 95% credible region). This shows great prospects of using lensed SNe Ia to obtain intrinsic early-phase SN spectra for deciphering SN Ia progenitors. As a demonstration of the usefulness of lensed SNe Ia for cosmology, we simulate a sample of mock lensed SN Ia systems that are expected to have accurate and precise time-delay measurements in the era of the Legacy Survey of Space and Time (LSST). Adopting realistic yet conservative uncertainties on their time-delay distances and lens angular diameter distances (of 6.6% and 5%, respectively), we find that a sample of 20 lensed SNe Ia would allow a constraint on the Hubble constant ($H_0$) with 1.3% uncertainty in the flat $\Lambda$CDM cosmology. We find a similar constraint on $H_0$ in an open $\Lambda$CDM cosmology, while the constraint degrades to 3% in a flat wCDM cosmology. We anticipate lensed SNe to be an independent and powerful probe of SN physics and cosmology in the upcoming LSST era.

We present the results of 15 years of monitoring of lensed quasars conducted by the COSMOGRAIL program at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed QSOs are presented for the first time. We complement our data set with other monitoring data available in the literature to measure the time delays in 18 systems, among which 9 reach a relative precision better than 15% for at least one time delay. To achieve this, we develop an automated version of the curve-shifting toolbox PyCS to ensure robust estimation of the time delay in the presence of microlensing while accounting for the errors due to imperfect representation of microlensing. We also reanalyze the previously published time delays of RX J1131$-$1231 and HE 0435$-$1223, adding respectively 6 and 2 new seasons of monitoring and confirming the previous time-delay measurements. When the time delay measurement is possible, we correct the light curves of the lensed images from their time delay and present the difference curves to highlight the microlensing signal contained in the data. This is to date the largest sample of decade-long lens monitoring data, useful to measure $H_0$, to measure the size of quasar accretion disks with microlensing, and to study quasar variability.

We use thirteen seasons of R-band photometry from the 1.2m Leonard Euler Swiss Telescope at La Silla to examine microlensing variability in the quadruply-imaged lensed quasar WFI 2026-4536. The lightcurves exhibit ${\sim}\,0.2\,\text{mag}$ of uncorrelated variability across all epochs and a prominent single feature of ${\sim}\,0.1\,\text{mag}$ within a single season. We analyze this variability to constrain the size of the quasar's accretion disk. Adopting a nominal inclination of 60$^\text{o}$, we find an accretion disk scale radius of $\log(r_s/\text{cm}) = 15.74^{+0.34}_{-0.29}$ at a rest-frame wavelength of $2043\,\unicode{xC5}$, and we estimate a black hole mass of $\log(M_{\text{BH}}/M_{\odot}) = 9.18^{+0.39}_{-0.34}$, based on the CIV line in VLT spectra. This size measurement is fully consistent with the Quasar Accretion Disk Size - Black Hole Mass relation, providing another system in which the accretion disk is larger than predicted by thin disk theory.

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 lens mass model of the quadruply-imaged gravitationally lensed quasar WFI2033-4723, and perform a blind cosmographical analysis based on this system. Our analysis combines (1) time-delay measurements from 14 years of data obtained by the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL) collaboration, (2) high-resolution $\textit{Hubble Space Telescope}$ imaging, (3) a measurement of the velocity dispersion of the lens galaxy based on ESO-MUSE data, and (4) multi-band, wide-field imaging and spectroscopy characterizing the lens environment. We account for all known sources of systematics, including the influence of nearby perturbers and complex line-of-sight structure, as well as the parametrization of the light and mass profiles of the lensing galaxy. After unblinding, we determine the effective time-delay distance to be $4784_{-248}^{+399}~\mathrm{Mpc}$, an average precision of $6.6\%$. This translates to a Hubble constant $H_{0} = 71.6_{-4.9}^{+3.8}~\mathrm{km~s^{-1}~Mpc^{-1}}$, assuming a flat $\Lambda$CDM cosmology with a uniform prior on $\Omega_\mathrm{m}$ in the range [0.05, 0.5]. This work is part of the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration, and the full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper (H0LiCOW XIII).

Flux ratio anomalies in strong gravitationally lensed quasars constitute a unique way to probe the abundance of non-luminous dark matter haloes, and hence the nature of dark matter. In this paper we identify double imaged quasars as a statistically efficient probe of dark matter, since they are 20 times more abundant than quadruply imaged quasars. Using N-body simulations that include realistic baryonic feedback, we measure the full distribution of flux ratios in doubly imaged quasars for cold (CDM) and warm dark matter (WDM) cosmologies. Through this method, we fold in two key systematics - quasar variability and line-of-sight structures. We find that WDM cosmologies predict a ~6 per cent difference in the cumulative distribution functions of flux ratios relative to CDM, with CDM predicting many more small ratios. Finally, we estimate that ~600 doubly imaged quasars will need to be observed in order to be able to unambiguously discern between CDM and the two WDM models studied here. Such sample sizes will be easily within reach of future large scale surveys such as Euclid. In preparation for this survey data we require discerning the scale of the uncertainties in modelling lens galaxies and their substructure in simulations, plus a strong understanding of the selection function of observed lensed quasars.

With the advent of high-cadence and multi-band photometric monitoring facilities, continuum reverberation mapping is becoming of increasing importance to measure the physical size of quasar accretion disks. The method is based on the measurement of the time it takes for a signal to propagate from the center to the outer parts of the central engine, assuming the continuum light curve at a given wavelength has a time shift of the order of a few days with respect to light curves obtained at shorter wavelengths. We show that with high-quality light curves, this assumption is not valid anymore and that light curves at different wavelengths are not only shifted in time but also distorted: in the context of the lamp-post model and thin-disk geometry, the multi-band light curves are in fact convolved by a transfer function whose size increase with wavelength. We illustrate the effect with simulated light curves in the LSST ugrizy bands and examine the impact on the delay measurements when using three different methods, namely JAVELIN, CREAM, and PyCS. We find that current accretion disk sizes estimated from JAVELIN and PyCS are underestimated by $\sim30\%$ and that unbiased measurement are only obtained with methods that properly take the skewed transfer functions into account, as the CREAM code does. With the LSST-like light curves, we expect to achieve measurement errors below $5\%$ with typical 2-day photometric cadence.

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.

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of earth-massed habitable worlds around the nearest starts, and also into distant Milky way objects up to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry ESA missions: NEAT for M3, micro-NEAT for S1 mission, and Theia for M4 and M5. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this white paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review quickly instrumentation and mission profiles.

The upcoming Large Synoptic Survey Telescope (LSST) will detect many strongly lensed Type Ia supernovae (LSNe Ia) for time-delay cosmography. This will provide an independent and direct way for measuring the Hubble constant $H_0$, which is necessary to address the current $4.4 \sigma$ tension in $H_0$ between the local distance ladder and the early Universe measurements. We present a detailed analysis of different observing strategies for the LSST, and quantify their impact on time-delay measurement between multiple images of LSNe Ia. For this, we produced microlensed mock-LSST light curves for which we estimated the time delay between different images. We find that using only LSST data for time-delay cosmography is not ideal. Instead, we advocate using LSST as a discovery machine for LSNe Ia, enabling time delay measurements from follow-up observations from other instruments in order to increase the number of systems by a factor of 2 to 16 depending on the observing strategy. Furthermore, we find that LSST observing strategies, which provide a good sampling frequency (the mean inter-night gap is around two days) and high cumulative season length (ten seasons with a season length of around 170 days per season), are favored. Rolling cadences subdivide the survey and focus on different parts in different years; these observing strategies trade the number of seasons for better sampling frequency. In our investigation, this leads to half the number of systems in comparison to the best observing strategy. Therefore rolling cadences are disfavored because the gain from the increased sampling frequency cannot compensate for the shortened cumulative season length. We anticipate that the sample of lensed SNe Ia from our preferred LSST cadence strategies with rapid follow-up observations would yield an independent percent-level constraint on $H_0$.

In this work, we provide an efficient and realistic data-driven approach to simulate astronomical images using deep generative models from machine learning. Our solution is based on a variant of the generative adversarial network (GAN) with progressive training methodology and Wasserstein cost function. The proposed solution generates naturalistic images of galaxies that show complex structures and high diversity, which suggests that data-driven simulations using machine learning can replace many of the expensive model-driven methods used in astronomical data processing.

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

Context. The precise determination of the present-day expansion rate of the Universe, expressed through the Hubble constant $H_0$, is one of the most pressing challenges in modern cosmology. Assuming flat $\Lambda$CDM, $H_0$ inference at high redshift using cosmic-microwave-background data from Planck disagrees at the 4.4$\sigma$ level with measurements based on the local distance ladder made up of parallaxes, Cepheids and Type Ia supernovae (SNe Ia), often referred to as "Hubble tension". Independent, cosmological-model-insensitive ways to infer $H_0$ are of critical importance. Aims. We apply an inverse-distance-ladder approach, combining strong-lensing time-delay-distance measurements with SN Ia data. By themselves, SNe Ia are merely good relative distance indicators, but by anchoring them to strong gravitational lenses one can obtain an $H_0$ measurement that is relatively insensitive to other cosmological parameters. Methods. A cosmological parameter estimate is performed for different cosmological background models, both for strong-lensing data alone and for the combined lensing + SNe Ia data sets. Results. The cosmological-model dependence of strong-lensing $H_0$ measurements is significantly mitigated through the inverse distance ladder. In combination with SN Ia data, the inferred $H_0$ consistently lies around 73-74 km s$^{-1}$ Mpc$^{-1}$, regardless of the assumed cosmological background model. Our results agree nicely with those from the local distance ladder, but there is a >2$\sigma$ tension with Planck results, and a ~1.5$\sigma$ discrepancy with results from an inverse distance ladder including Planck, Baryon Acoustic Oscillations and SNe Ia. Future strong-lensing distance measurements will reduce the uncertainties in $H_0$ from our inverse distance ladder.

We present the measurement of the Hubble Constant, $H_0$, with three strong gravitational lens systems. We describe a blind analysis of both PG1115+080 and HE0435-1223 as well as an extension of our previous analysis of RXJ1131-1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP AO effort, with Hubble Space Telescope (HST) imaging, velocity dispersion measurements, and a description of the line-of-sight mass distribution to build an accurate and precise lens mass model. This mass model is then combined with the COSMOGRAIL measured time delays in these systems to determine $H_{0}$. We do both an AO-only and an AO+HST analysis of the systems and find that AO and HST results are consistent. After unblinding, the AO-only analysis gives $H_{0}=82.8^{+9.4}_{-8.3}~\rm km\,s^{-1}\,Mpc^{-1}$ for PG1115+080, $H_{0}=70.1^{+5.3}_{-4.5}~\rm km\,s^{-1}\,Mpc^{-1}$ for HE0435-1223, and $H_{0}=77.0^{+4.0}_{-4.6}~\rm km\,s^{-1}\,Mpc^{-1}$ for RXJ1131-1231. The joint AO-only result for the three lenses is $H_{0}=75.6^{+3.2}_{-3.3}~\rm km\,s^{-1}\,Mpc^{-1}$. The joint result of the AO+HST analysis for the three lenses is $H_{0}=76.8^{+2.6}_{-2.6}~\rm km\,s^{-1}\,Mpc^{-1}$. All of the above results assume a flat $\Lambda$ cold dark matter cosmology with a uniform prior on $\Omega_{\textrm{m}}$ in [0.05, 0.5] and $H_{0}$ in [0, 150] $\rm km\,s^{-1}\,Mpc^{-1}$. This work is a collaboration of the SHARP and H0LiCOW teams, and shows that AO data can be used as the high-resolution imaging component in lens-based measurements of $H_0$. The full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper.

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.

We present new measurements of the time delays of WFI2033-4723. The data sets used in this work include 14 years of data taken at the 1.2m Leonhard Euler Swiss telescope, 13 years of data from the SMARTS 1.3m telescope at Las Campanas Observatory and a single year of high-cadence and high-precision monitoring at the MPIA 2.2m telescope. The time delays measured from these different data sets, all taken in the R-band, are in good agreement with each other and with previous measurements from the literature. Combining all the time-delay estimates from our data sets results in Dt_AB = 36.2-0.8+0.7 days (2.1% precision), Dt_AC = -23.3-1.4+1.2 days (5.6%) and Dt_BC = -59.4-1.3+1.3 days (2.2%). In addition, the close image pair A1-A2 of the lensed quasars can be resolved in the MPIA 2.2m data. We measure a time delay consistent with zero in this pair of images. We also explore the prior distributions of microlensing time-delay potentially affecting the cosmological time-delay measurements of WFI2033-4723. There is however no strong indication in our measurements that microlensing time delay is neither present nor absent. This work is part of a H0LiCOW series focusing on measuring the Hubble constant from WFI2033-4723.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

We present a sample of 12 Quasi-Stellar Objects (QSOs) potentially acting as strong gravitational lenses on background Emission Line Galaxies (ELG) or Lyman-$\alpha$ Emitters (LAEs) selected selected through a systematic search of the 297301 QSOs in the Sloan Digital Sky Survey (SDSS)-III Data Release 12. Candidates are identified by looking for compound spectra, where emission lines at a redshift larger than that of the quasar can be identified in the residuals after a QSO spectral template is subtracted from the observed spectra. The narrow diameter of BOSS fibers (2") then ensures that the object responsible for the additional emission lines must lie close to the line of sight of the QSO and hence provides a high probability of lensing. Among the 12 candidates identified, 9 have definite evidence for the presence of a background ELG identified by at least 4 higher-redshift nebular emission lines. The remaining $3$ probable candidates present a strong asymmetrical emission line attributed to a background Lyman-$\alpha$ emitter (LAE). The QSO-ELG (QSO-LAE) lens candidates have QSO lens redshifts in the range $0.24\lesssim z_{QSO} \lesssim 0.66$ ($0.75 \lesssim z_{QSO} \lesssim 1.23$ ) and background galaxy redshifts in the range $0.48 \lesssim z_{S,ELG} \lesssim 0.94$ ($2.17 \lesssim z_{S,LAE} \lesssim 4.48$). We show that the algorithmic search is complete at >90% for QSO-ELG systems, whereas it falls at 40-60% for QSO-LAE, depending on the redshift of the source. Upon confirmation of the lensing nature of the systems, this sample may quadruple the number of known QSOs acting as strong lenses. We have determined the completeness of our search, which allows future studies to compute lensing probabilities of galaxies by QSOs and differentiate between different QSO models. (Abridged)

It has recently been proposed that gravitationally lensed type-Ia supernovae can provide microlensing-free time-delay measurements provided that the measurement is taken during the achromatic expansion phase of the explosion and that color light curves are used rather than single-band light curves. If verified, this would provide both precise and accurate time-delay measurements, making lensed type-Ia supernovae a new golden standard for time-delay cosmography. However, the 3D geometry of the expanding shell can introduce an additional bias that has not yet been fully explored. In this work, we present and discuss the impact of this effect on time-delay cosmography with lensed supernovae and find that on average it leads to a bias of a few tenths of a day for individual lensed systems. This is negligible in view of the cosmological time delays predicted for typical lensed type-Ia supernovae but not for the specific case of the recently discovered type-Ia supernova iPTF16geu, whose time delays are expected to be smaller than a day.

Cosmic shear is a primary cosmological probe for several present and upcoming surveys investigating dark matter and dark energy, such as Euclid or WFIRST. The probe requires an extremely accurate measurement of the shapes of millions of galaxies based on imaging data. Crucially, the shear measurement must address and compensate for a range of interwoven nuisance effects related to the instrument optics and detector, noise, unknown galaxy morphologies, colors, blending of sources, and selection effects. This paper explores the use of supervised machine learning (ML) as a tool to solve this inverse problem. We present a simple architecture that learns to regress shear point estimates and weights via shallow artificial neural networks. The networks are trained on simulations of the forward observing process, and take combinations of moments of the galaxy images as inputs. A challenging peculiarity of this ML application is the combination of the noisiness of the input features and the requirements on the accuracy of the inverse regression. To address this issue, the proposed training algorithm minimizes bias over multiple realizations of individual source galaxies, reducing the sensitivity to properties of the overall sample of source galaxies. Importantly, an observational selection function of these source galaxies can be straightforwardly taken into account via the weights. We first introduce key aspects of our approach using toy-model simulations, and then demonstrate its potential on images mimicking Euclid data. Finally, we analyze images from the GREAT3 challenge, obtaining competitively low shear biases despite the use of a simple training set. We conclude that the further development of ML approaches is of high interest to meet the stringent requirements on the shear measurement in current and future surveys. A demonstration implementation of our technique is publicly available.

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

We present 13 seasons of $R$-band photometry of the quadruply-lensed quasar WFI 2033-4723 from the 1.3m SMARTS telescope at CTIO and the 1.2m Euler Swiss Telescope at La Silla, in which we detect microlensing variability of $\sim0.2$ mags on a timescale of $\sim$6 years. Using a Bayesian Monte Carlo technique, we analyze the microlensing signal to obtain a measurement of the size of this system's accretion disk of $\log (r_s/{\rm cm}) = 15.86^{+0.25}_{-0.27}$ at $\lambda_{rest} = 2481{\rm \AA}$, assuming a $60^\circ$ inclination angle. We confirm previous measurements of the BC and AB time delays, and we obtain a tentative measurement of the delay between the closely spaced A1 and A2 images of $\Delta t_{A1A2} = t_{A1} - t_{A2} = -3.9^{+3.4}_{-2.2}$ days. We conclude with an update to the Quasar Accretion Disk Size - Black Hole Mass Relation, in which we confirm that the accretion disk size predictions from simple thin disk theory are too small.

Strong gravitational lensing offers a wealth of astrophysical information on the background source it affects, provided the lensed source can be reconstructed as if it was seen in the absence of lensing. In the present work, we illustrate how sparse optimisation can address the problem. As a first step towards a full free-form lens modelling technique, we consider linear inversion of the lensed source under sparse regularisation and joint deblending from the lens light profile. The method is based on morphological component analysis, assuming a known mass model. We show with numerical experiments that representing the lens and source light using an undecimated wavelet basis allows us to reconstruct the source and to separate it from the foreground lens at the same time. Both the source and lens light have a non-analytic form, allowing for the flexibility needed in the inversion to represent arbitrarily small and complex luminous structures in the lens and source. in addition, sparse regularisation avoids over-fitting the data and does not require the use of any adaptive mesh or pixel grid. As a consequence, our reconstructed sources can be represented on a grid of very small pixels. Sparse regularisation in the wavelet domain also allows for automated computation of the regularisation parameter, thus minimising the impact of arbitrary choice of initial parameters. Our inversion technique for a fixed mass distribution can be incorporated in future lens modelling technique iterating over the lens mass parameters. The python package corresponding to the algorithms described in this article can be downloaded via the github platform at https://github.com/herjy/SLIT.

The primary goals of the STRong lensing Insights into the Dark Energy Survey (STRIDES) collaboration are to measure the dark energy equation of state parameter and the free streaming length of dark matter. To this aim, STRIDES is discovering strongly lensed quasars in the imaging data of the Dark Energy Survey and following them up to measure time delays, high resolution imaging, and spectroscopy sufficient to construct accurate lens models. In this paper, we first present forecasts for STRIDES. Then, we describe the STRIDES classification scheme, and give an overview of the Fall 2016 follow-up campaign. We continue by detailing the results of two selection methods, the Outlier Selection Technique and a morphological algorithm, and presenting lens models of a system, which could possibly be a lensed quasar in an unusual configuration. We conclude with the summary statistics of the Fall 2016 campaign. Including searches presented in companion papers (Anguita et al.; Ostrovski et al.), STRIDES followed up 117 targets identifying 7 new strongly lensed systems, and 7 nearly identical quasars (NIQs), which could be confirmed as lenses by the detection of the lens galaxy. 76 candidates were rejected and 27 remain otherwise inconclusive, for a success rate in the range 6-35\%. This rate is comparable to that of previous searches like SQLS even though the parent dataset of STRIDES is purely photometric and our selection of candidates cannot rely on spectroscopic information.

We present a study of cold gas absorption from a damped Lyman-$\alpha$ absorber (DLA) at redshift $z_{\rm abs}=1.946$ towards two lensed images of the quasar J144254.78+405535.5 at redshift $z_{\rm QSO} = 2.590$. The physical separation of the two lines of sight at the absorber redshift is $d_{\rm abs}=0.7$~kpc based on our lens model. We observe absorption lines from neutral carbon and H$_2$ along both lines of sight indicating that cold gas is present on scales larger than $d_{\rm abs}$. We measure column densities of HI to be $\log N(\rm H\,\i) = 20.27\pm0.02$ and $20.34\pm0.05$ and of H$_2$ to be $\log N(\rm H_2) = 19.7\pm0.1$ and $19.9\pm0.2$. The metallicity inferred from sulphur is consistent with Solar metallicity for both sightlines: $[{\rm S/H}]_A = 0.0\pm0.1$ and $[{\rm S/H}]_B = -0.1\pm0.1$. Based on the excitation of low rotational levels of H$_2$, we constrain the temperature of the cold gas phase to be $T=109\pm20$ and $T=89\pm25$ K for the two lines of sight. From the relative excitation of fine-structure levels of CI, we constrain the hydrogen volumetric densities in the range of $40-110$ cm$^{-3}$. Based on the ratio of observed column density and volumetric density, we infer the average individual `cloud' size along the line of sight to be $l\approx0.1$ pc. Using the transverse line-of-sight separation of 0.7 kpc together with the individual cloud size, we are able to put an upper limit to the volume filling factor of cold gas of $f_{\rm vol} < 0.2$ %. Nonetheless, the projected covering fraction of cold gas must be large (close to unity) over scales of a few kpc in order to explain the presence of cold gas in both lines of sight. Compared to the typical extent of DLAs (~10-30 kpc), this is consistent with the relative incidence rate of CI absorbers and DLAs.

Time-delay strong lensing provides a unique way to directly measure the Hubble constant ($H_{0}$). The precision of the $H_{0}$ measurement depends on the uncertainties in the time-delay measurements, the mass distribution of the main deflector(s), and the mass distribution along the line of sight. Tie and Kochanek (2018) have proposed a new microlensing effect on time delays based on differential magnification of the coherent accretion disc variability of the lensed quasar. If real, this effect could significantly broaden the uncertainty on the time delay measurements by up to $30\%$ for lens systems such as PG1115+080, which have relatively short time delays and monitoring over several different epochs. In this paper we develop a new technique that uses the time-delay ratios and simulated microlensing maps within a Bayesian framework in order to limit the allowed combinations of microlensing delays and thus to lessen the uncertainties due to the proposed effect. We show that, under the assumption of Tie and Kochanek (2018), the uncertainty on the time-delay distance ($D_{\Delta t}$, which is proportional to 1/$H_{0}$) of short time-delay ($\sim18$ days) lens, PG1115+080, increases from $\sim7\%$ to $\sim10\%$ by simultaneously fitting the three time-delay measurements from the three different datasets across twenty years, while in the case of long time-delay ($\sim90$ days) lens, the microlensing effect on time delays is negligible as the uncertainty on $D_{\Delta t}$ of RXJ1131-1231 only increases from $\sim2.5\%$ to $\sim2.6\%$.

We review three distance measurement techniques beyond the local universe: (1) gravitational lens time delays, (2) baryon acoustic oscillation (BAO), and (3) HI intensity mapping. We describe the principles and theory behind each method, the ingredients needed for measuring such distances, the current observational results, and future prospects. Time delays from strongly lensed quasars currently provide constraints on $H_0$ with < 4% uncertainty, and with 1% within reach from ongoing surveys and efforts. Recent exciting discoveries of strongly lensed supernovae hold great promise for time-delay cosmography. BAO features have been detected in redshift surveys up to z <~ 0.8 with galaxies and z ~ 2 with Ly-$\alpha$ forest, providing precise distance measurements and $H_0$ with < 2% uncertainty in flat $\Lambda$CDM. Future BAO surveys will probe the distance scale with percent-level precision. HI intensity mapping has great potential to map BAO distances at z ~ 0.8 and beyond with precisions of a few percent. The next years ahead will be exciting as various cosmological probes reach 1% uncertainty in determining $H_0$, to assess the current tension in $H_0$ measurements that could indicate new physics.

We present time-delay estimates for the quadruply imaged quasar PG 1115+080. Our resuls are based on almost daily observations for seven months at the ESO MPIA 2.2m telescope at La Silla Observatory, reaching a signal-to-noise ratio of about 1000 per quasar image. In addition, we re-analyse existing light curves from the literature that we complete with an additional three seasons of monitoring with the Mercator telescope at La Palma Observatory. When exploring the possible source of bias we consider the so-called microlensing time delay, a potential source of systematic error so far never directly accounted for in previous time-delay publications. In fifteen years of data on PG 1115+080, we find no strong evidence of microlensing time delay. Therefore not accounting for this effect, our time-delay estimates on the individual data sets are in good agreement with each other and with the literature. Combining the data sets, we obtain the most precise time-delay estimates to date on PG 1115+080, with Dt(AB) = 8.3+1.5-1.6 days (18.7% precision), Dt(AC) = 9.9+1.1-1.1 days (11.1%) and Dt(BC) = 18.8+1.6-1.6 days (8.5%). Turning these time delays into cosmological constraints is done in a companion paper that makes use of ground-based Adaptive Optics (AO) with the Keck telescope.

We present ALMA 2-mm continuum and CO (2-1) spectral line imaging of the gravitationally lensed z=0.654 star-forming/quasar composite RX J1131-1231 at 240-400 mas angular resolution. The continuum emission is found to be compact and coincident with the optical emission, whereas the molecular gas forms a complete Einstein ring, which shows strong differential magnification. The de-lensed source structure is determined on 400-pc resolution using a visibility-fitting lens modelling technique. The reconstructed molecular gas velocity-field is consistent with a rotating disk with a maximum rotational velocity of 280 km/s. From dynamical model fitting we find an enclosed mass M(r<5 kpc)=(1.46+/-0.31)*10^11 M_sol. The molecular gas distribution is highly structured, with clumps that are co-incident with higher gas velocity dispersion regions 40-50 km/s and with the intensity peaks in the optical emission, which are associated with sites of on-going turbulent star-formation. The peak in the CO (2-1) distribution is not co-incident with the AGN, where there is a paucity of molecular gas emission, possibly due to radiative feedback from the central engine. The intrinsic molecular gas luminosity is L'_CO=(1.2+/-0.3)*10^10 K km/s pc^2 and the inferred gas mass is M(H2)=(8.3+/-3.0)*10^10 M_sol, which given its dynamical mass is consistent with a CO-H2 conversion factor of alpha = 5.5+/-2.0 M_solar(K km/s pc^2)^-1. This suggests that the star-formation efficiency is dependent on the host galaxy morphology as opposed to the nature of the AGN. The far-infrared continuum spectral energy distribution shows evidence for heated dust, equivalent to an obscured star-formation rate of SFR=69^(+41)_(-25)*(7.3/u_IR)M_sol/yr, which demonstrates the composite star-forming/AGN nature of this system. RX J1131-1231

We present a weak gravitational lensing measurement of the external convergence along the line of sight to the quadruply lensed quasar HE$\,$0435$-$1223. Using deep r-band images from Subaru-Suprime-Cam we observe galaxies down to a 3$\sigma$ limiting magnitude of $\sim 26$ mags resulting in a source galaxy density of 14 galaxies / arcmin$^2$ after redshift-based cuts. Using an inpainting technique and Multi-Scale Entropy filtering algorithm, we find that the region in close proximity to the lens has an estimated external convergence of $\kappa=-0.012^{+0.020}_{-0.013}$ and is hence marginally under-dense. We also rule out the presence of any halo with a mass greater than $M_{\rm vir}=1.6\times10^{14}h^{-1}M_\odot$ (68$\%$ confidence limit). Our results, consistent with previous studies of this lens, confirm that the intervening mass along the line of sight to HE$\,$0435$-$1223 does not affect significantly the cosmological results inferred from the time delay measurements of that specific object.

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 a new formulation of a family of proximity operators that generalize the projector step for phase retrieval. These proximity operators for noisy intensity measurements can replace the classical "noise free" projection in any projection-based algorithm. They are derived from a maximum likelihood formulation and admit closed form solutions for both the Gaussian and the Poisson cases. In addition, we extend these proximity operators to undersampled intensity measurements. To assess their performance, these operators are exploited in a classical Gerchberg Saxton algorithm. We present numerical experiments showing that the reconstructed complex amplitudes with these proximity operators perform always better than using the classical intensity projector while their computational overhead is moderate.

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.

Measuring a weak gravitational lensing signal to the level required by the next generation of space-based surveys demands exquisite reconstruction of the point-spread function (PSF). However, unresolved binary stars can significantly distort the PSF shape. In an effort to mitigate this bias, we aim at detecting unresolved binaries in realistic Euclid stellar populations. We tested methods in numerical experiments where (i) the PSF shape is known to Euclid requirements across the field of view, and (ii) the PSF shape is unknown. We drew simulated catalogues of PSF shapes for this proof-of-concept paper. Following the Euclid survey plan, the objects were observed four times. We propose three methods to detect unresolved binary stars. The detection is based on the systematic and correlated biases between exposures of the same object. One method is a simple correlation analysis, while the two others use supervised machine-learning algorithms (random forest and artificial neural network). In both experiments, we demonstrate the ability of our methods to detect unresolved binary stars in simulated catalogues. The performance depends on the level of prior knowledge of the PSF shape and the shape measurement errors. Good detection performances are observed in both experiments. Full complexity, in terms of the images and the survey design, is not included, but key aspects of a more mature pipeline are discussed. Finding unresolved binaries in objects used for PSF reconstruction increases the quality of the PSF determination at arbitrary positions. We show, using different approaches, that we are able to detect at least binary stars that are most damaging for the PSF reconstruction process.

We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The new spectroscopic data enable better subtraction of foreground light, and better identification of multiple background images. The inferred distribution of dark matter is consistent with being centered on the galaxies, as expected by LCDM. Each galaxy's dark matter also appears to be symmetric. Whilst we do not find an offset between mass and light (suggestive of self-interacting dark matter) as previously reported, the numerical simulations that have been performed to calibrate Abell 3827 indicate that offsets and asymmetry are still worth looking for in collisions with particular geometries. Meanwhile, ALMA proves exceptionally useful for strong lens image identifications.

A striking signal of dark matter beyond the standard model is the existence of cores in the centre of galaxy clusters. Recent simulations predict that a Brightest Cluster Galaxy (BCG) inside a cored galaxy cluster will exhibit residual wobbling due to previous major mergers, long after the relaxation of the overall cluster. This phenomenon is absent with standard cold dark matter where a cuspy density profile keeps a BCG tightly bound at the centre. We test this hypothesis using cosmological simulations and deep observations of 10 galaxy clusters acting as strong gravitational lenses. Modelling the BCG wobble as a simple harmonic oscillator, we measure the wobble amplitude, $A_w$, in the BAHAMAS suite of cosmological hydrodynamical simulations, finding an upper limit for the CDM paradigm of $A_w < 2$kpc at the 95% confidence limit. We carry out the same test on the data finding a non-zero amplitude of $A_w=11.82^{+7.3}_{-3.0}$kpc, with the observations dis-favouring $A_w = 0$ at the 3$\sigma$ confidence level. This detection of BCG wobbling is evidence for a dark matter core at the heart of galaxy clusters. It also shows that strong lensing models of clusters cannot assume that the BCG is exactly coincident with the large scale halo. While our small sample of galaxy clusters already indicates a non-zero $A_w$, with larger surveys, e.g. Euclid, we will be able to not only to confirm the effect but also to use it to determine whether or not the wobbling finds its origin in new fundamental physics or astrophysical process.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

Accurate astronomical distance determination is crucial for all fields in astrophysics, from Galactic to cosmological scales. Despite, or perhaps because of, significant efforts to determine accurate distances, using a wide range of methods, tracers, and techniques, an internally consistent astronomical distance framework has not yet been established. We review current efforts to homogenize the Local Group's distance framework, with particular emphasis on the potential of RR Lyrae stars as distance indicators, and attempt to extend this in an internally consistent manner to cosmological distances. Calibration based on Type Ia supernovae and distance determinations based on gravitational lensing represent particularly promising approaches. We provide a positive outlook to improvements to the status quo expected from future surveys, missions, and facilities. Astronomical distance determination has clearly reached maturity and near-consistency.

The radio-quiet quasar Q2059-360 at redshift $z=3.08$ is known to be close to a small Lyman $\alpha$ blob (LAB) and to be absorbed by a proximate damped Ly$\alpha$ (PDLA) system. Here, we present the Multi Unit Spectroscopic Explorer (MUSE) integral field spectroscopy follow-up of this quasi-stellar object (QSO). Our primary goal is to characterize this LAB in detail by mapping it both spatially and spectrally using the Ly$\alpha$ line, and by looking for high-ionization lines to constrain the emission mechanism. Combining the high sensitivity of the MUSE integral field spectrograph mounted on the Yepun telescope at ESO-VLT with the natural coronagraph provided by the PDLA, we map the LAB down to the QSO position, after robust subtraction of QSO light in the spectral domain. In addition to confirming earlier results for the small bright component of the LAB, we unveil a faint filamentary emission protruding to the south over about 80 pkpc (physical kpc); this results in a total size of about 120 pkpc. We derive the velocity field of the LAB (assuming no transfer effects) and map the Ly$\alpha$ line width. Upper limits are set to the flux of the N V $\lambda 1238-1242$, C IV $\lambda 1548-1551$, He II $\lambda 1640$, and C III] $\lambda 1548-1551$ lines. We have discovered two probable Ly$\alpha$ emitters at the same redshift as the LAB and at projected distances of 265 kpc and 207 kpc from the QSO; their Ly$\alpha$ luminosities might well be enhanced by the QSO radiation. We also find an emission line galaxy at $z=0.33$ near the line of sight to the QSO. This LAB shares the same general characteristics as the 17 others surrounding radio-quiet QSOs presented previously. However, there are indications that it may be centered on the PDLA galaxy rather than on the QSO.

Simulations of self-interacting dark matter (SIDM) predict that dark matter should lag behind galaxies during a collision. If the interaction is mediated by a high-mass force carrier, the distribution of dark matter can also develop asymmetric dark matter tails. To search for this asymmetry, we compute the gravitational lensing properties of a mass distribution with a free {\em skewness} parameter. We apply this to the dark matter around the four central galaxies in cluster Abell~3827. In the galaxy whose dark matter peak has previously been found to be offset, we tentatively measure a skewness $s=0.23^{+0.05}_{-0.22}$ in the same direction as the peak offset. Our method may be useful in future gravitational lensing analyses of colliding galaxy clusters and merging galaxies.

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.

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.

Strong gravitational lenses with measured time delays between the multiple images allow a direct measurement of the time-delay distance to the lens, and thus a measure of cosmological parameters, particularly the Hubble constant, $H_{0}$. We present a blind lens model analysis of the quadruply-imaged quasar lens HE 0435-1223 using deep Hubble Space Telescope imaging, updated time-delay measurements from the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL), a measurement of the velocity dispersion of the lens galaxy based on Keck data, and a characterization of the mass distribution along the line of sight. HE 0435-1223 is the third lens analyzed as a part of the $H_{0}$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) project. We account for various sources of systematic uncertainty, including the detailed treatment of nearby perturbers, the parameterization of the galaxy light and mass profile, and the regions used for lens modeling. We constrain the effective time-delay distance to be $D_{\Delta t} = 2612_{-191}^{+208}~\mathrm{Mpc}$, a precision of 7.6%. From HE 0435-1223 alone, we infer a Hubble constant of $H_{0} = 73.1_{-6.0}^{+5.7}~\mathrm{km~s^{-1}~Mpc^{-1}}$ assuming a flat $\Lambda$CDM cosmology. The cosmographic inference based on the three lenses analyzed by H0LiCOW to date is presented in a companion paper (H0LiCOW Paper V).

The empirical correlation between the mass of a super-massive black hole (MBH) and its host galaxy properties is widely considered to be evidence of their co-evolution. A powerful way to test the co-evolution scenario and learn about the feedback processes linking galaxies and nuclear activity is to measure these correlations as a function of redshift. Unfortunately, currently MBH can only be estimated in active galaxies at cosmological distances. At these distances, bright active galactic nuclei (AGN) can outshine the host galaxy, making it extremely difficult to measure the host's luminosity. Strongly lensed AGNs provide in principle a great opportunity to improve the sensitivity and accuracy of the host galaxy luminosity measurements as the host galaxy is magnified and more easily separated from the point source, provided the lens model is sufficiently accurate. In order to measure the MBH-L correlation with strong lensing, it is necessary to ensure that the lens modelling is accurate, and that the host galaxy luminosity can be recovered to at least a precision and accuracy better than that of the typical MBH measurement. We carry out extensive and realistic simulations of deep Hubble Space Telescope observations of lensed AGNs obtained by our collaboration. We show that the host galaxy luminosity can be recovered with better accuracy and precision than the typical uncertainty on MBH(~ 0.5 dex) for hosts as faint as 2-4 magnitudes dimmer than the AGN itself. Our simulations will be used to estimate bias and uncertainties on the actual measurements to be presented in a future paper.

We present a study on galaxy detection and shape classification using topometric clustering algorithms. We first use the DBSCAN algorithm to extract, from CCD frames, groups of adjacent pixels with significant fluxes and we then apply the DENCLUE algorithm to separate the contributions of overlapping sources. The DENCLUE separation is based on the localization of pattern of local maxima, through an iterative algorithm which associates each pixel to the closest local maximum. Our main classification goal is to take apart elliptical from spiral galaxies. We introduce new sets of features derived from the computation of geometrical invariant moments of the pixel group shape and from the statistics of the spatial distribution of the DENCLUE local maxima patterns. Ellipticals are characterized by a single group of local maxima, related to the galaxy core, while spiral galaxies have additional ones related to segments of spiral arms. We use two different supervised ensemble classification algorithms, Random Forest, and Gradient Boosting. Using a sample of ~ 24000 galaxies taken from the Galaxy Zoo 2 main sample with spectroscopic redshifts, and we test our classification against the Galaxy Zoo 2 catalog. We find that features extracted from our pipeline give on average an accuracy of ~ 93%, when testing on a test set with a size of 20% of our full data set, with features deriving from the angular distribution of density attractor ranking at the top of the discrimination power.

We propose a wavefront-based method to estimate the PSF over the whole field of view. This method estimate the aberrations of all the mirrors of the telescope using only field stars. In this proof of concept paper, we described the method and present some qualitative results.

We present the results of a new search for galaxy-scale strong lensing systems in CFHTLS Wide. Our lens-finding technique involves a preselection of potential lens galaxies, applying simple cuts in size and magnitude. We then perform a Principal Component Analysis of the galaxy images, ensuring a clean removal of the light profile. Lensed features are searched for in the residual images using the clustering topometric algorithm DBSCAN. We find 1098 lens candidates that we inspect visually, leading to a cleaned sample of 109 new lens candidates. Using realistic image simulations we estimate the completeness of our sample and show that it is independent of source surface brightness, Einstein ring size (image separation) or lens redshift. We compare the properties of our sample to previous lens searches in CFHTLS. Including the present search, the total number of lenses found in CFHTLS amounts to 678, which corresponds to ~4 lenses per square degree down to i=24.8. This is equivalent to ~ 60.000 lenses in total in a survey as wide as Euclid, but at the CFHTLS resolution and depth.

Information on the spectral types of stars is of great interest in view of the exploitation of space-based imaging surveys. In this article, we investigate the classification of stars into spectral types using only the shape of their diffraction pattern in a single broad-band image. We propose a supervised machine learning approach to this endeavour, based on principal component analysis (PCA) for dimensionality reduction, followed by artificial neural networks (ANNs) estimating the spectral type. Our analysis is performed with image simulations mimicking the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) in the F606W and F814W bands, as well as the Euclid VIS imager. We first demonstrate this classification in a simple context, assuming perfect knowledge of the point spread function (PSF) model and the possibility of accurately generating mock training data for the machine learning. We then analyse its performance in a fully data-driven situation, in which the training would be performed with a limited subset of bright stars from a survey, and an unknown PSF with spatial variations across the detector. We use simulations of main-sequence stars with flat distributions in spectral type and in signal-to-noise ratio, and classify these stars into 13 spectral subclasses, from O5 to M5. Under these conditions, the algorithm achieves a high success rate both for Euclid and HST images, with typical errors of half a spectral class. Although more detailed simulations would be needed to assess the performance of the algorithm on a specific survey, this shows that stellar classification from single-band images is well possible.

We introduce a new algorithm for colour separation and deblending of multi-band astronomical images called MuSCADeT which is based on Morpho-spectral Component Analysis of multi-band images. The MuSCADeT algorithm takes advantage of the sparsity of astronomical objects in morphological dictionaries such as wavelets and their differences in spectral energy distribution (SED) across multi-band observations. This allows us to devise a model independent and automated approach to separate objects with different colours. We show with simulations that we are able to separate highly blended objects and that our algorithm is robust against SED variations of objects across the field of view. To confront our algorithm with real data, we use HST images of the strong lensing galaxy cluster MACS J1149+2223 and we show that MuSCADeT performs better than traditional profile-fitting techniques in deblending the foreground lensing galaxies from background lensed galaxies. Although the main driver for our work is the deblending of strong gravitational lenses, our method is fit to be used for any purpose related to deblending of objects in astronomical images. An example of such an application is the separation of the red and blue stellar populations of a spiral galaxy in the galaxy cluster Abell 2744. We provide a python package along with all simulations and routines used in this paper to contribute to reproducible research efforts. Codes can be found at http://lastro.epfl.ch/page-126973.html.

We present a two-channel deconvolution method that decomposes images into a parametric point-source channel and a pixelized extended-source channel. Based on the central idea of the deconvolution algorithm proposed by Magain, Courbin & Sohy (1998), the method aims at improving the resolution of the data rather than at completely removing the point spread function (PSF). Improvements over the original method include a better regularization of the pixel channel of the image, based on wavelet filtering and multiscale analysis, and a better controlled separation of the point source vs. the extended source. In addition, the method is able to simultaneously deconvolve many individual frames of the same object taken with different instruments under different PSF conditions. For this purpose, we introduce a general geometric transformation between individual images. This transformation allows the combination of the images without having to interpolate them. We illustrate the capability of our algorithm using real and simulated images with complex diffraction-limited PSF.

COSMOGRAIL is a long-term photometric monitoring of gravitationally lensed QSOs aimed at implementing Refsdal's time-delay method to measure cosmological parameters, in particular H0. Given long and well sampled light curves of strongly lensed QSOs, time-delay measurements require numerical techniques whose quality must be assessed. To this end, and also in view of future monitoring programs or surveys such as the LSST, a blind signal processing competition named Time Delay Challenge 1 (TDC1) was held in 2014. The aim of the present paper, which is based on the simulated light curves from the TDC1, is double. First, we test the performance of the time-delay measurement techniques currently used in COSMOGRAIL. Second, we analyse the quantity and quality of the harvest of time delays obtained from the TDC1 simulations. To achieve these goals, we first discover time delays through a careful inspection of the light curves via a dedicated visual interface. Our measurement algorithms can then be applied to the data in an automated way. We show that our techniques have no significant biases, and yield adequate uncertainty estimates resulting in reduced chi2 values between 0.5 and 1.0. We provide estimates for the number and precision of time-delay measurements that can be expected from future time-delay monitoring campaigns as a function of the photometric signal-to-noise ratio and of the true time delay. We make our blind measurements on the TDC1 data publicly available

We present a detailed study of the colours in late-type galaxy discs for ten of the EDisCS galaxy clusters with 0.5 < z < 0.8. Our cluster sample contains 172 spiral galaxies, and our control sample is composed of 96 field disc galaxies. We deconvolve their ground-based V and I images obtained with FORS2 at the VLT with initial spatial resolutions between 0.4 and 0.8 arcsec to achieve a final resolution of 0.1 arcsec with 0.05 arcsec pixels, which is close to the resolution of the ACS at the HST. After removing the central region of each galaxy to avoid pollution by the bulges, we measured the V-I colours of the discs. We find that 50% of cluster spiral galaxies have disc V-I colours redder by more than 1 sigma of the mean colours of their field counterparts. This is well above the 16% expected for a normal distribution centred on the field disc properties. The prominence of galaxies with red discs depends neither on the mass of their parent cluster nor on the distance of the galaxies to the cluster cores. Passive spiral galaxies constitute 20% of our sample. These systems are not abnormally dusty. They are are made of old stars and are located on the cluster red sequences. Another 24% of our sample is composed of galaxies that are still active and star forming, but less so than galaxies with similar morphologies in the field. These galaxies are naturally located in the blue sequence of their parent cluster colour-magnitude diagrams. The reddest of the discs in clusters must have stopped forming stars more than ~5 Gyr ago. Some of them are found among infalling galaxies, suggesting preprocessing. Our results confirm that galaxies are able to continue forming stars for some significant period of time after being accreted into clusters, and suggest that star formation can decline on seemingly long (1 to 5 Gyr) timescales.

The next generation of space-based telescopes used for weak lensing surveys will require exquisite point spread function (PSF) determination. Previously negligible effects may become important in the reconstruction of the PSF, in part because of the improved spatial resolution. In this paper, we show that unresolved multiple star systems can affect the ellipticity and size of the PSF and that this effect is not cancelled even when using many stars in the reconstruction process. We estimate the error in the reconstruction of the PSF due to the binaries in the star sample both analytically and with image simulations for different PSFs and stellar populations. The simulations support our analytical finding that the error on the size of the PSF is a function of the multiple stars distribution and of the intrinsic value of the size of the PSF, i.e. if all stars were single. Similarly, the modification of each of the complex ellipticity components (e1,e2) depends on the distribution of multiple stars and on the intrinsic complex ellipticity. Using image simulations, we also show that the predicted error in the PSF shape is a theoretical limit that can be reached only if large number of stars (up to thousands) are used together to build the PSF at any desired spatial position. For a lower number of stars, the PSF reconstruction is worse. Finally, we compute the effect of binarity for different stellar magnitudes and show that bright stars alter the PSF size and ellipticity more than faint stars. This may affect the design of PSF calibration strategies and the choice of the related calibration fields.

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 analyze the optical, UV, and X-ray microlensing variability of the lensed quasar SDSS J0924+0219 using six epochs of Chandra data in two energy bands (spanning 0.4-8.0 keV, or 1-20 keV in the quasar rest frame), 10 epochs of F275W (rest-frame 1089A) Hubble Space Telescope data, and high-cadence R-band (rest-frame 2770A) monitoring spanning eleven years. Our joint analysis provides robust constraints on the extent of the X-ray continuum emission region and the projected area of the accretion disk. The best-fit half-light radius of the soft X-ray continuum emission region is between 5x10^13 and 10^15 cm, and we find an upper limit of 10^15 cm for the hard X-rays. The best-fit soft-band size is about 13 times smaller than the optical size, and roughly 7 GM_BH/c^2 for a 2.8x10^8 M_sol black hole, similar to the results for other systems. We find that the UV emitting region falls in between the optical and X-ray emitting regions at 10^14 cm < r_1/2,UV < 3x10^15 cm. Finally, the optical size is significantly larger, by 1.5*sigma, than the theoretical thin-disk estimate based on the observed, magnification-corrected I-band flux, suggesting a shallower temperature profile than expected for a standard disk.

We compared the total mass density profiles of three different types of galaxies using weak gravitational lensing: (i) 29 galaxies that host quasars at z~0.32 that are in a post-starburst (PSQ) phase with high star formation indicating recent merger activity, (ii) 22 large elliptical galaxies from the SLACS sample that do not host a quasar at z~0.23, and (iii) 17 galaxies that host moderately luminous quasars at z~0.36 powered by disk instabilities, but with no intense star formation. On an initial test we found no evidence for a connection between the merger state of a galaxy and the profile of the halo, with the PSQ profile comparable to that of the other two samples and consistent with the Leauthaud et al. (2014) study of moderately luminous quasars in COSMOS. Given the compatibility of the two quasar samples, we combined these and found no evidence for any connection between black hole activity and the dark matter halo. All three mass profiles remained compatible with isothermality given the present data.

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the S\'{e}rsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $\chi^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $\chi^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f = $20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.

Weak gravitational lensing has the potential to place tight constraints on the equation of the state of dark energy. However, this will only be possible if shear measurement methods can reach the required level of accuracy. We present a new method to measure the ellipticity of galaxies used in weak lensing surveys. The method makes use of direct deconvolution of the data by the total Point Spread Function (PSF). We adopt a linear algebra formalism that represents the PSF as a Toeplitz matrix. This allows us to solve the convolution equation by applying the Hopfield Neural Network iterative scheme. The ellipticity of galaxies in the deconvolved images are then measured using second order moments of the autocorrelation function of the images. To our knowledge, it is the first time full image deconvolution is used to measure weak lensing shear. We apply our method to the simulated weak lensing data proposed in the GREAT10 challenge and obtain a quality factor of Q=87. This result is obtained after applying image denoising to the data, prior to the deconvolution. The additive and multiplicative biases on the shear power spectrum are then +0.000009 and +0.0357, respectively.

We present a weak lensing mass map covering ~124 square degrees of the Canada-France-Hawaii Telescope Stripe 82 Survey (CS82). We study the statistics of rare peaks in the map, including peak abundance, the peak-peak correlation functions and the tangential-shear profiles around peaks. We find that the abundance of peaks detected in CS82 is consistent with predictions from a Lambda-CDM cosmological model, once noise effects are properly included. The correlation functions of peaks with different signal-to-noise ratio (SNR) are well described by power laws, and there is a clear cross-correlation between the Sloan Digital Sky Survey III/Constant Mass galaxies and high SNR peaks. The tangential-shear profiles around peaks increase with peak SNR. We fit analytical models to the tangential-shear profiles, including a projected singular isothermal sphere (SIS) model and a projected Navarro, Frenk & White (NFW) model, plus a two-halo term. For the high SNR peaks, the SIS model is rejected at ~3-sigma. The NFW model plus a two-halo term gives more acceptable fits to the data. Some peaks match the positions of optically detected clusters, while others are relatively dark. Comparing dark and matched peaks, we find a difference in lensing signal of a factor of 2, suggesting that about half of the dark peaks are false detections.

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 present an algorithm using Principal Component Analysis (PCA) to subtract galaxies from imaging data, and also two algorithms to find strong, galaxy-scale gravitational lenses in the resulting residual image. The combined method is optimized to find full or partial Einstein rings. Starting from a pre-selection of potential massive galaxies, we first perform a PCA to build a set of basis vectors. The galaxy images are reconstructed using the PCA basis and subtracted from the data. We then filter the residual image with two different methods. The first uses a curvelet (curved wavelets) filter of the residual images to enhance any curved/ring feature. The resulting image is transformed in polar coordinates, centered on the lens galaxy center. In these coordinates, a ring is turned into a line, allowing us to detect very faint rings by taking advantage of the integrated signal-to-noise in the ring (a line in polar coordinates). The second way of analysing the PCA-subtracted images identifies structures in the residual images and assesses whether they are lensed images according to their orientation, multiplicity and elongation. We apply the two methods to a sample of simulated Einstein rings, as they would be observed with the ESA Euclid satellite in the VIS band. The polar coordinates transform allows us to reach a completeness of 90% and a purity of 86%, as soon as the signal-to-noise integrated in the ring is higher than 30, and almost independent of the size of the Einstein ring. Finally, we show with real data that our PCA-based galaxy subtraction scheme performs better than traditional subtraction based on model fitting to the data. Our algorithm can be developed and improved further using machine learning and dictionary learning methods, which would extend the capabilities of the method to more complex and diverse galaxy shapes.

The GRavitational lEnsing Accuracy Testing 3 (GREAT3) challenge is the third in a series of image analysis challenges, with a goal of testing and facilitating the development of methods for analyzing astronomical images that will be used to measure weak gravitational lensing. This measurement requires extremely precise estimation of very small galaxy shape distortions, in the presence of far larger intrinsic galaxy shapes and distortions due to the blurring kernel caused by the atmosphere, telescope optics, and instrumental effects. The GREAT3 challenge is posed to the astronomy, machine learning, and statistics communities, and includes tests of three specific effects that are of immediate relevance to upcoming weak lensing surveys, two of which have never been tested in a community challenge before. These effects include realistically complex galaxy models based on high-resolution imaging from space; spatially varying, physically-motivated blurring kernel; and combination of multiple different exposures. To facilitate entry by people new to the field, and for use as a diagnostic tool, the simulation software for the challenge is publicly available, though the exact parameters used for the challenge are blinded. Sample scripts to analyze the challenge data using existing methods will also be provided. See http://great3challenge.info and http://great3.projects.phys.ucl.ac.uk/leaderboard/ for more information.

We present the results of a visual search for galaxy-scale gravitational lenses in nearly 7 square degrees of Hubble Space Telescope (HST) images. The dataset comprises the whole imaging data ever taken with the Advanced Camera for Surveys (ACS) in the filter F814W (I-band) up to August 31st, 2011, i.e. 6.03 square degrees excluding the field of the Cosmic Evolution Survey (COSMOS) which has been the subject of a separate visual search. In addition, we have searched for lenses in the whole Wide Field Camera 3 (WFC3) near-IR imaging dataset in all filters (1.01 square degrees) up to the same date. Our primary goal is to provide a sample of lenses with a broad range of different morphologies and lens-source brightness contrast in order estimate a lower limit to the number of galaxy-scale strong lenses in the future Euclid survey in its VIS band. Our criteria to select lenses are purely morphological as we do not use any colour or redshift information.The final candidate selection is very conservative hence leading to a nearly pure but incomplete sample. We find 49 new lens candidates: 40 in the ACS images and 9 in the WFC3 images. Out of these, 16 candidates are secure lenses owing to their striking morphology, 21 more are very good candidates, and 12 more have morphologies compatible with gravitational lensing but also compatible with other astrophysical objects. It is therefore insensitive to cosmic variance and allows to estimate the number of galaxy-scale strong lenses on the sky for a putative survey depth, which is the main result of the present work. Because of the incompleteness of the sample, the estimated lensing rates should be taken as lower limits. Using these, we anticipate that a 15 000 square degrees space survey such as Euclid will find at least 60 000 galaxy-scale strong lenses down to a limiting AB magnitude of I = 24.5 (10-sigma) or I = 25.8 (3-sigma).

This paper presents optical R-band light curves and the time delay of the doubly imaged gravitationally lensed quasar SDSS J1001+5027 at a redshift of 1.838. We have observed this target for more than six years, between March 2005 and July 2011, using the 1.2-m Mercator Telescope, the 1.5-m telescope of the Maidanak Observatory, and the 2-m Himalayan Chandra Telescope. Our resulting light curves are composed of 443 independent epochs, and show strong intrinsic quasar variability, with an amplitude of the order of 0.2 magnitudes. From this data, we measure the time delay using five different methods, all relying on distinct approaches. One of these techniques is a new development presented in this paper. All our time-delay measurements are perfectly compatible. By combining them, we conclude that image A is leading B by 119.3 +/- 3.3 days (1 sigma, 2.8% uncertainty), including systematic errors. It has been shown recently that such accurate time-delay measurements offer a highly complementary probe of dark energy and spatial curvature, as they independently constrain the Hubble constant. The next mandatory step towards using SDSS J1001+5027 in this context will be the measurement of the velocity dispersion of the lensing galaxy, in combination with deep Hubble Space Telescope imaging.

We present the results from nine years of optically monitoring the gravitationally lensed z=0.658 quasar RX J1131-1231. The R-band light curves of the four individual images of the quasar were obtained using deconvolution photometry for a total of 707 epochs. Several sharp quasar variability features strongly constrain the time delays between the quasar images. Using three different numerical techniques, we measure these delays for all possible pairs of quasar images while always processing the four light curves simultaneously. For all three methods, the delays between the three close images A, B, and C are compatible with being 0, while we measure the delay of image D to be 91 days, with a fractional uncertainty of 1.5% (1 sigma), including systematic errors. Our analysis of random and systematic errors accounts in a realistic way for the observed quasar variability, fluctuating microlensing magnification over a broad range of temporal scales, noise properties, and seasonal gaps. Finally, we find that our time-delay measurement methods yield compatible results when applied to subsets of the data.

Measuring time delays between the multiple images of gravitationally lensed quasars is now recognized as a competitive way to constrain the cosmological parameters, and it is complementary with other cosmological probes. This requires long and well sampled optical light curves of numerous lensed quasars, such as those obtained by the COSMOGRAIL collaboration. High-quality data from our monitoring campaign call for novel numerical techniques to robustly measure the delays, as well as the associated random and systematic uncertainties, even in the presence of microlensing variations. We propose three different point estimators to measure time delays, which are explicitly designed to handle light curves with extrinsic variability. These methods share a common formalism, which enables them to process data from n-image lenses. Since the estimators rely on significantly contrasting ideas, we expect them to be sensitive to different bias sources. For each method and data set, we empirically estimate both the precision and accuracy (bias) of the time delay measurement using simulated light curves with known time delays that closely mimic the observations. Finally, we test the self-consistency of our approach, and we demonstrate that our bias estimation is serviceable. These new methods, including the empirical uncertainty estimator, will represent the standard benchmark for analyzing the COSMOGRAIL light curves.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

Aims. Within the framework of the COSMOGRAIL collaboration we present 7- and 8.5-year-long light curves and time-delay estimates for two gravitationally lensed quasars: SDSS J1206+4332 and HS 2209+1914. Methods. We monitored these doubly lensed quasars in the R-band using four telescopes: the Mercator, Maidanak, Himalayan Chandra, and Euler Telescopes, together spanning a period of 7 to 8.5 observing seasons from mid-2004 to mid-2011. The photometry of the quasar images was obtained through simultaneous deconvolution of these data. The time delays were determined from these resulting light curves using four very different techniques: a dispersion method, a spline fit, a regression difference technique, and a numerical model fit. This minimizes the bias that might be introduced by the use of a single method. Results. The time delay for SDSS J1206+4332 is Delta_t AB = 111.3 +/- 3 days with A leading B, confirming a previously published result within the error bars. For HS 2209+1914 we present a new time delay of Delta_t BA = 20.0 +/- 5 days with B leading A. Conclusions. The combination of data from up to four telescopes have led to well-sampled and nearly 9-season-long light curves, which were necessary to obtain these results, especially for the compact doubly lensed quasar HS 2209+1914.

In this paper we present results from the weak lensing shape measurement GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. This marks an order of magnitude step change in the level of scrutiny employed in weak lensing shape measurement analysis. We provide descriptions of each method tested and include 10 evaluation metrics over 24 simulation branches. GREAT10 was the first shape measurement challenge to include variable fields; both the shear field and the Point Spread Function (PSF) vary across the images in a realistic manner. The variable fields enable a variety of metrics that are inaccessible to constant shear simulations including a direct measure of the impact of shape measurement inaccuracies, and the impact of PSF size and ellipticity, on the shear power spectrum. To assess the impact of shape measurement bias for cosmic shear we present a general pseudo-Cl formalism, that propagates spatially varying systematics in cosmic shear through to power spectrum estimates. We also show how one-point estimators of bias can be extracted from variable shear simulations. The GREAT10 Galaxy Challenge received 95 submissions and saw a factor of 3 improvement in the accuracy achieved by shape measurement methods. The best methods achieve sub-percent average biases. We find a strong dependence in accuracy as a function of signal-to-noise, and indications of a weak dependence on galaxy type and size. Some requirements for the most ambitious cosmic shear experiments are met above a signal-to-noise ratio of 20. These results have the caveat that the simulated PSF was a ground-based PSF. Our results are a snapshot of the accuracy of current shape measurement methods and are a benchmark upon which improvement can continue. This provides a foundation for a better understanding of the strengths and limitations of shape measurement methods.

Strong gravitational lenses with measured time delays between the multiple images and models of the lens mass distribution allow a one-step determination of the time-delay distance, and thus a measure of cosmological parameters. We present a blind analysis of the gravitational lens RXJ1131-1231 incorporating (1) the newly measured time delays from COSMOGRAIL, (2) archival HST imaging of the lens system, (3) a new velocity-dispersion measurement of the lens galaxy of 323+/-20km/s based on Keck spectroscopy, and (4) a characterization of the line-of-sight structures via observations of the lens' environment and ray tracing through the Millennium Simulation. Our blind analysis is designed to prevent experimenter bias. The joint analysis of the data sets allows a time-delay distance measurement to 6% precision that takes into account all known systematic uncertainties. In combination with the WMAP7 data set in flat wCDM cosmology, our unblinded cosmological constraints for RXJ1131-1231 are: H_0=80.0+5.8/-5.7km/s/Mpc, OmegaDE=0.79+/-0.03 and w=-1.25+0.17/-0.21. We find the results to be statistically consistent with those from the analysis of the gravitational lens B1608+656. The joint constraints from the two lenses and WMAP7 are H_0=75.2+4.4/-4.2km/s/Mpc, OmegaDE=0.76+0.02/-0.03 and w=-1.14+0.17/-0.20 in flat wCDM, and H_0=73.1+2.4/-3.6km/s/Mpc, OmegaL=0.75+0.01/-0.02 and OmegaK=0.003+0.005/-0.006 in open LCDM. Time-delay lenses constrain especially tightly the Hubble constant (5.7% and 4.0% respectively in wCDM and open LCDM) and curvature of the universe. They complement well other cosmological probes, and provide an independent check of unknown systematics. Our measurement of the Hubble constant is completely independent of those based on the local distance ladder method, providing an important consistency check of the standard cosmological model and of general relativity.

Microlensing is a powerful technique which can be used to study the continuum and the broad line emitting regions in distant AGNs. After a brief description of the methods and required data, we present recent applications of this technique. We show that microlensing allows one to measure the temperature profile of the accretion disc, estimate the size and study the geometry of the region emitting the broad emission lines.

It is anticipated that the large sky areas covered by planned wide-field weak lensing surveys will reduce statistical errors to such an extent that systematic errors will instead become the dominant source of uncertainty. It is therefore crucial to devise numerical methods to measure galaxy shapes with the least possible systematic errors. We present a simple "forward deconvolution" method, gfit, to measure galaxy shapes given telescope and atmospheric smearings, in the presence of noise. The method consists in fitting a single 2D elliptical S\'ersic profile to the data, convolved with the point spread function. We applied gfit to the data proposed in the GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. In spite of its simplicity, gfit obtained the lowest additive bias {($\sqrt{\mathcal{A}}=0.057\times10^{-4}$)} on the shear power spectrum among twelve different methods and the second lowest multiplicative bias {($\mathcal{M}/2=0.583\times10^{-2}$)}. It remains that gfit is a fitting method and is therefore affected by noise bias. However, the simplicity of the underlying galaxy model combined with the use of an efficient customized minimization algorithm allow very competitive performances, at least on the GREAT10 data, for a relatively low computing time.

Planned wide-field weak lensing surveys are expected to reduce the statistical errors on the shear field to unprecedented levels. In contrast, systematic errors like those induced by the convolution with the point spread function (PSF) will not benefit from that scaling effect and will require very accurate modeling and correction. While numerous methods have been devised to carry out the PSF correction itself, modeling of the PSF shape and its spatial variations across the instrument field of view has, so far, attracted much less attention. This step is nevertheless crucial because the PSF is only known at star positions while the correction has to be performed at any position on the sky. A reliable interpolation scheme is therefore mandatory and a popular approach has been to use low-order bivariate polynomials. In the present paper, we evaluate four other classical spatial interpolation methods based on splines (B-splines), inverse distance weighting (IDW), radial basis functions (RBF) and ordinary Kriging (OK). These methods are tested on the Star-challenge part of the GRavitational lEnsing Accuracy Testing 2010 (GREAT10) simulated data and are compared with the classical polynomial fitting (Polyfit). We also test all our interpolation methods independently of the way the PSF is modeled, by interpolating the GREAT10 star fields themselves (i.e., the PSF parameters are known exactly at star positions). We find in that case RBF to be the clear winner, closely followed by the other local methods, IDW and OK. The global methods, Polyfit and B-splines, are largely behind, especially in fields with (ground-based) turbulent PSFs. In fields with non-turbulent PSFs, all interpolators reach a variance on PSF systematics $\sigma_{sys}^2$ better than the $1\times10^{-7}$ upper bound expected by future space-based surveys, with the local interpolators performing better than the global ones.

We present the results from the first public blind PSF reconstruction challenge, the GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Star Challenge. Reconstruction of a spatially varying PSF, sparsely sampled by stars, at non-star positions is a critical part in the image analysis for weak lensing where inaccuracies in the modelled ellipticity and size-squared can impact the ability to measure the shapes of galaxies. This is of importance because weak lensing is a particularly sensitive probe of dark energy, and can be used to map the mass distribution of large scale structure. Participants in the challenge were presented with 27,500 stars over 1300 images subdivided into 26 sets, where in each set a category change was made in the type or spatial variation of the PSF. Thirty submissions were made by 9 teams. The best methods reconstructed the PSF with an accuracy of ~0.00025 in ellipticity and ~0.00074 in size squared. For a fixed pixel scale narrower PSFs were found to be more difficult to model than larger PSFs, and the PSF reconstruction was severely degraded with the inclusion of an atmospheric turbulence model (although this result is likely to be a strong function of the amplitude of the turbulence power spectrum).

We describe a simple and fast method to correct ellipticity measurements of galaxies from the distortion by the instrumental and atmospheric point spread function (PSF), in view of weak lensing shear measurements. The method performs a classification of galaxies and associated PSFs according to measured shape parameters, and corrects the measured galaxy ellipticites by querying a large lookup table (LUT), built by supervised learning. We have applied this new method to the GREAT10 image analysis challenge, and present in this paper a refined solution that obtains the competitive quality factor of Q = 104, without any shear power spectrum denoising or training. Of particular interest is the efficiency of the method, with a processing time below 3 ms per galaxy on an ordinary CPU.

When an image of a strongly lensed quasar is microlensed, the different components of its spectrum are expected to be differentially magnified owing to the different sizes of the corresponding emitting region. Chromatic changes are expected to be observed in the continuum while the emission lines should be deformed as a function of the size, geometry and kinematics of the regions from which they originate. Microlensing of the emission lines has been reported only in a handful of systems so far. In this paper we search for microlensing deformations of the optical spectra of pairs of images in 17 lensed quasars. This sample is composed of 13 pairs of previously unpublished spectra and four pairs of spectra from literature. Our analysis is based on a spectral decomposition technique which allows us to isolate the microlensed fraction of the flux independently of a detailed modeling of the quasar emission lines. Using this technique, we detect microlensing of the continuum in 85% of the systems. Among them, 80% show microlensing of the broad emission lines. Focusing on the most common lines in our spectra (CIII] and MgII) we detect microlensing of either the blue or the red wing, or of both wings with the same amplitude. This observation implies that the broad line region is not in general spherically symmetric. In addition, the frequent detection of microlensing of the blue and red wings independently but not simultaneously with a different amplitude, does not support existing microlensing simulations of a biconical outflow. Our analysis also provides the intrinsic flux ratio between the lensed images and the magnitude of the microlensing affecting the continuum. These two quantities are particularly relevant for the determination of the fraction of matter in clumpy form in galaxies and for the detection of dark matter substructures via the identification of flux ratio anomalies.

We present four new seasons of optical monitoring data and six epochs of X-ray photometry for the doubly-imaged lensed quasar Q J0158-4325. The high-amplitude, short-period microlensing variability for which this system is known has historically precluded a time delay measurement by conventional methods. We attempt to circumvent this limitation by application of a Monte Carlo microlensing analysis technique, but we are only able to prove that the delay must have the expected sign (image A leads image B). Despite our failure to robustly measure the time delay, we successfully model the microlensing at optical and X-ray wavelengths to find a half light radius for soft X-ray emission log(r_{1/2,X,soft}/cm) = 14.3^{+0.4}_{-0.5}, an upper limit on the half-light radius for hard X-ray emission log(r_{1/2,X,hard}/cm) <= 14.6 and a refined estimate of the inclination-corrected scale radius of the optical R-band (rest frame 3100 Angstrom) continuum emission region of log(r_s/cm) = 15.6+-0.3.

What are the frequency, shape, kinematics, and luminosity of Ly\alpha\ envelopes surrounding radio-quiet quasars at high redshift, and is the luminosity of these envelopes related to that of the quasar or not? As a first step towards answering these questions, we have searched for Ly\alpha\ envelopes around six radio-quiet quasars at z~4.5, using deep spectra taken with the FORS2 spectrograph attached to the UT1 of the Very Large Telescope (VLT). Using the multi-slit mode allows us to observe several point spread function stars simultaneously with the quasar, and to remove the point-like emission from the quasar, unveiling the faint underlying Ly\alpha\ envelope with unprecedented depth. An envelope is detected around four of the six quasars, which suggests that these envelopes are very frequent. Their diameter varies in the range 26<d<64 kpc, their surface brightness in the range 3x10^{-19}<\mu<2x10^{-17} erg/s/cm^2/arcsec^2, and their luminosity in the range 10^{42}<L(Ly\alpha)<10^{44} erg/s. Their shape may be strongly asymmetric. The Ly\alpha\ emission line full width at half maximum (FWHM) is 900<FWHM<2200 km/s and its luminosity correlates with that of the broad line region (BLR) of the quasar, with the notable exception of BR2237-0607, the brightest object in our sample. The same holds for the relation between the envelope Ly\alpha\ luminosity and the ionizing luminosity of the quasar. While the deep slit spectroscopy presented in this paper is very efficient at detecting very faint Ly\alpha\ envelopes, narrow-band imaging is now needed to measure accurately their spatial extent, radial luminosity profile, and total luminosity. These observables are crucial to help us discriminate between the three possible radiation processes responsible for the envelope emission: (i) cold accretion, (ii) fluorescence induced by the quasar, and (iii) scattering of the BLR photons by cool gas.

We report the outcome of a 3-day workshop on the Hubble constant (H_0) that took place during February 6-8 2012 at the Kavli Institute for Particle Astrophysics and Cosmology, on the campus of Stanford University. The participants met to address the following questions. Are there compelling scientific reasons to obtain more precise and more accurate measurements of H_0 than currently available? If there are, how can we achieve this goal? The answers that emerged from the workshop are (1) better measurements of H_0 provide critical independent constraints on dark energy, spatial curvature of the Universe, neutrino physics, and validity of general relativity, (2) a measurement of H_0 to 1% in both precision and accuracy, supported by rigorous error budgets, is within reach for several methods, and (3) multiple paths to independent determinations of H_0 are needed in order to access and control systematics.

(abridged) Gravitationally lensed quasars can be used as powerful cosmological and astrophysical probes. We can (i) infer the Hubble constant based on the time-delay technique, (ii) unveil substructures along the l.o.s. toward distant galaxies, and (iii) compare the shape and the slope of baryons and dark matter distributions in galaxies. To reach these goals, we need high-accuracy astrometry and morphology measurements of the lens. In this work, we first present new astrometry for 11 lenses with measured time delays. Using MCS deconvolution on NIC2 HST images, we reached an astrometric accuracy of about 1-2.5 mas and an accurate shape measurement of the lens galaxy. Second, we combined these measurements with those of 14 other systems to present new mass models of these lenses. This led to the following results: 1) In 4 double-image quasars, we show that the influence of the lens environment on the time delay can easily be quantified and modeled, hence putting these lenses with high priority for time-delay determination. 2) For quadruple-image quasars, the difficulty often encountered in reproducing the image positions to milli-arcsec accuracy (astrometric anomaly) is overcome by explicitly including the nearest visible galaxy in the model. However, one anomalous system (J1131-1231) does not show any luminous perturber in its vicinity, and three others (WFI2026-4536, WFI2033-4723, and B2045+265) have problematic modeling. These 4 systems are the best candidates for a pertubation by a dark matter substructure. 3) We find a significant correlation between the PA of the light and of the mass distributions in lensing galaxies. In contrast with other studies, we find that the ellipticity of the light and of the mass also correlate well, suggesting that the overall spatial distribution of matter is not very different from the baryon distribution in the inner \sim 5 kpc of lensing galaxies.

GRavitational lEnsing Accuracy Testing 2010 (GREAT10) is a public image analysis challenge aimed at the development of algorithms to analyze astronomical images. Specifically, the challenge is to measure varying image distortions in the presence of a variable convolution kernel, pixelization and noise. This is the second in a series of challenges set to the astronomy, computer science and statistics communities, providing a structured environment in which methods can be improved and tested in preparation for planned astronomical surveys. GREAT10 extends upon previous work by introducing variable fields into the challenge. The "Galaxy Challenge" involves the precise measurement of galaxy shape distortions, quantified locally by two parameters called shear, in the presence of a known convolution kernel. Crucially, the convolution kernel and the simulated gravitational lensing shape distortion both now vary as a function of position within the images, as is the case for real data. In addition, we introduce the "Star Challenge" that concerns the reconstruction of a variable convolution kernel, similar to that in a typical astronomical observation. This document details the GREAT10 Challenge for potential participants. Continually updated information is also available from http://www.greatchallenges.info.

We report the discovery of three new cases of QSOs acting as strong gravitational lenses on background emission line galaxies: SDSS J0827+5224 (zQSO = 0.293, zs = 0.412), SDSS J0919+2720 (zQSO = 0.209, zs = 0.558), SDSS J1005+4016 (zQSO = 0.230, zs = 0.441). The selection was carried out using a sample of 22,298 SDSS spectra displaying at least four emission lines at a redshift beyond that of the foreground QSO. The lensing nature is confirmed from Keck imaging and spectroscopy, as well as from HST/WFC3 imaging in the F475W and F814W filters. Two of the QSOs have face-on spiral host galaxies and the third is a QSO+galaxy pair. The velocity dispersion of the host galaxies, inferred from simple lens modeling, is between \sigma_v = 210 and 285 km/s, making these host galaxies comparable in mass with the SLACS sample of early-type strong lenses.

Aims. Ly-alpha emitters (LAEs) can be detected out to very high redshifts during the epoch of reionization. The evolution of the LAE luminosity function with redshift is a direct probe of the Ly-alpha transmission of the intergalactic medium (IGM), and therefore of the IGM neutral-hydrogen fraction. Measuring the Ly-alpha luminosity function (LF) of LAEs at redshift z = 7.7 therefore allows us to constrain the ionizing state of the Universe at this redshift. Methods. We observed three 7.5'x7.5' fields with the HAWK-I instrument at the VLT with a narrow band filter centred at 1.06 $\mu$m and targeting LAEs at redshift z ~ 7.7. The fields were chosen for the availability of multiwavelength data. One field is a galaxy cluster, the Bullet Cluster, which allowed us to use gravitational amplification to probe luminosities that are fainter than in the field. The two other fields are subareas of the GOODS Chandra Deep Field South and CFHTLS-D4 deep field. We selected z=7.7 LAE candidates from a variety of colour criteria, in particular from the absence of detection in the optical bands. Results. We do not find any LAE candidates at z = 7.7 in ~2.4 x 10^4 Mpc^3 down to a narrow band AB magnitude of ~ 26, which allows us to infer robust constraints on the Ly-alpha LAE luminosity function at this redshift. Conclusions. The predicted mean number of objects at z = 6.5, derived from somewhat different LFs of Hu et al. (2010), Ouchi et al. (2010), and Kashikawa et al. (2011) are 2.5, 13.7, and 11.6, respectively. Depending on which of these LFs we refer to, we exclude a scenario with no evolution from z = 6.5 to z = 7.7 at 85% confidence without requiring a strong change in the IGM Ly-alpha transmission, or at 99% confidence with a significant quenching of the IGM Ly-alpha transmission, possibly from a strong increase in the high neutral-hydrogen fraction between these two redshifts.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We present accurate time delays for the quadruply imaged quasar HE 0435-1223. The delays were measured from 575 independent photometric points obtained in the R-band between January 2004 and March 2010. With seven years of data, we clearly show that quasar image A is affected by strong microlensing variations and that the time delays are best expressed relative to quasar image B. We measured Delta_t(BC) = 7.8+/-0.8 days, Delta_t(BD) = -6.5+/-0.7 days and Delta_t_CD = -14.3+/-0.8 days. We spacially deconvolved HST NICMOS2 F160W images to derive accurate astrometry of the quasar images and to infer the light profile of the lensing galaxy. We combined these images with a stellar population fitting of a deep VLT spectrum of the lensing galaxy to estimate the baryonic fraction, $f_b$, in the Einstein radius. We measured f_b = 0.65+0.13-0.10 if the lensing galaxy has a Salpeter IMF and f_b = 0.45+0.04-0.07 if it has a Kroupa IMF. The spectrum also allowed us to estimate the velocity dispersion of the lensing galaxy, sigma_ap = 222+/-34 km/s. We used f_b and sigma_ap to constrain an analytical model of the lensing galaxy composed of an Hernquist plus generalized NFW profile. We solve the Jeans equations numerically for the model and explored the parameter space under the additional requirement that the model must predict the correct astrometry for the quasar images. Given the current error bars on f_b and sigma_ap, we did not constrain H0 yet with high accuracy, i.e., we found a broad range of models with chi^2 < 1. However, narrowing this range is possible, provided a better velocity dispersion measurement becomes available. In addition, increasing the depth of the current HST imaging data of HE 0435-1223 will allow us to combine our constraints with lens reconstruction techniques that make use of the full Einstein ring that is visible in this object.

Anomalous flux ratios between quasar images are suspected to be caused by substructures in lens galaxies. We present new deep and high resolution H and Ks imaging of the strongly lensed quasar SDSS J0924+0219 obtained using the ESO VLT with adaptive optics and the Laser Guide Star system. SDSS J0924+0219 is particularly interesting as the observed flux ratio between the quasar images vastly disagree with the predictions from smooth mass models. With our adaptive optics observations we find a luminous object, Object L, located ~0.3" to the North of the lens galaxy, but we show that it can not be responsible for the anomalous flux ratios. Object L as well as a luminous extension of the lens galaxy to the South are seen in the archival HST/ACS image in the F814W filter. This suggests that Object L is part of a bar in the lens galaxy, as also supported by the presence of a significant disk component in the light profile of the lens galaxy. Finally, we do not find evidence for any other luminous substructure that may explain the quasar images flux ratios. However, owe to the persistence of the flux ratio anomaly over time (~ 7 years) a combination of microlensing and milli-lensing is the favorite explanation for the observations.

Weak gravitational lensing is a very sensitive way of measuring cosmological parameters, including dark energy, and of testing current theories of gravitation. In practice, this requires exquisite measurement of the shapes of billions of galaxies over large areas of the sky, as may be obtained with the EUCLID and WFIRST satellites. For a given survey depth, applying image denoising to the data both improves the accuracy of the shape measurements and increases the number density of galaxies with a measurable shape. We perform simple tests of three different denoising techniques, using synthetic data. We propose a new and simple denoising method, based on wavelet decomposition of the data and a Wiener filtering of the resulting wavelet coefficients. When applied to the GREAT08 challenge dataset, this technique allows us to improve the quality factor of the measurement (Q; GREAT08 definition), by up to a factor of two. We demonstrate that the typical pixel size of the EUCLID optical channel will allow us to use image denoising.

We aim to use microlensing taking place in the lensed quasar Q2237+0305 to study the structure of the broad line region and measure the size of the region emitting the CIV and CIII] lines. Methods: Based on 39 spectrophotometric monitoring data points obtained between Oct. 2004 and Dec. 2007, we derived lightcurves for the CIV and CIII] emission lines. We used three different techniques to analyse the microlensing signal. Different components of the lines (narrow, broad and very broad) are identified and studied. We built a library of simulated microlensing lightcurves that reproduce the signal observed in the continuum and in the lines provided only the source size is changed. A Bayesian analysis scheme is then developed to derive the size of the various components of the BLR. Results: 1. The half-light radius of the region emitting the CIV line is found to be R_CIV ~ 66^{+110}_{-46} lt-days = 0.06$^{+0.09}_{-0.04}$ pc = 1.7$^{+2.8}_{-1.1}$\,10$^{17}$ cm (at 68.3% CI). Similar values are obtained for CIII]. Relative sizes of the carbon-line and V-band continuum emitting-regions are also derived with median values of R(line)/R(cont) in the range 4 to 29, depending of the FWHM of the line component. 2. The size of the CIV emitting region agrees with the Radius-Luminosity relationship derived from reverberation mapping. Using the virial theorem we derive the mass of the black hole in Q2237+0305 to be M_BH ~ 10^{8.3+/-0.3} M_sun. 3. We find that the CIV and CIII] lines are produced in at least 2 spatially distinct regions, the most compact one giving rise to the broadest component of the line. The broad and narrow line profiles are slightly different for CIV and CIII]. 4. Our analysis suggests a different structure for the CIV and FeII+III emitting regions, with the latter produced in the inner part of the BLR or in a less extended emitting region than CIV.

The low-mass end of the stellar Initial Mass Function (IMF) is constrained by focusing on the baryon-dominated central regions of strong lensing galaxies. We study in this letter the Einstein Cross (Q2237+0305), a z=0.04 barred galaxy whose bulge acts as lens on a background quasar. The positions of the four quasar images constrain the surface mass density on the lens plane, whereas the surface brightness (H-band NICMOS/HST imaging) along with deep spectroscopy of the lens (VLT/FORS1) allow us to constrain the stellar mass content, for a range of IMFs. We find that a classical single power law (Salpeter IMF) predicts more stellar mass than the observed lensing estimates. This result is confirmed at the 99% confidence level, and is robust to systematic effects due to the choice of population synthesis models, the presence of dust, or the complex disk/bulge population mix. Our non-parametric methodology is more robust than kinematic estimates, as we do not need to make any assumptions about the dynamical state of the galaxy or its decomposition into bulge and disk. Over a range of low-mass power law slopes (with Salpeter being Gamma=+1.35) we find that at a 90% confidence level, slopes with Gamma>0 are ruled out.

We present the first case of strong gravitational lensing by a QSO : SDSS J0013+1523, at z = 0.120. The discovery is the result of a systematic search for emission lines redshifted behind QSOs, among 22298 spectra of the SDSS data release 7. Apart from the z = 0.120 spectral features of the foreground QSO, the spectrum of SDSS J0013+1523 also displays the OII and Hbeta emission lines and the OIII doublet, all at the same redshift, z = 0.640. Using sharp Keck adaptive optics K-band images obtained using laser guide stars, we unveil two objects within a radius of 2 arcsec from the QSO. Deep Keck optical spectroscopy clearly confirms one of these objects at z = 0.640 and shows traces of the OIII, emission line of the second object, also at z = 0.640. Lens modeling suggests that they represent two images of the same z = 0.640 emission-line galaxy. Our Keck spectra also allow us to measure the redshift of an intervening galaxy at z = 0.394, located 3.2 arcsec away from the line of sight to the QSO. If the z = 0.120 QSO host galaxy is modeled as a singular isothermal sphere, its mass within the Einstein radius is M_E(r < 1 kpc) = 2.16e10 M_Sun and its velocity dispersion is sigma_SIS = 169 km/s. This is about 1 sigma away from the velocity dispersion estimated from the width of the QSO Hbeta emission line, sigma_*(M_BH) = 124 +/- 47 km/s. Deep optical HST imaging will be necessary to constrain the total radial mass profile of the QSO host galaxy using the detailed shape of the lensed source. This first case of a QSO acting as a strong lens on a more distant object opens new directions in the study of QSO host galaxies.

We present a direct detection of the stellar halo of the edge-on S0 galaxy NGC 3957, using ultra-deep VLT/VIMOS V and R images. This is achieved with a sky subtraction strategy based on infrared techniques. These observations allow us to reach unprecedented high signal-to-noise ratios up to 15 kpc away from the galaxy center, rendering photon-noise negligible. The 1 sigma detection limits are R = 30.6 mag/arcsec^2 and V = 31.4 mag/arcsec^2. We conduct a thorough analysis of the possible sources of systematic errors that could affect the data: flat-fielding, differences in CCD responses, scaling of the sky background, the extended halo itself, and PSF wings. We conclude that the V-R colour of the NGC 3957 halo, calculated between 5 and 8 kpc above the disc plane where the systematic errors are modest, is consistent with an old and preferentially metal-poor normal stellar population, like that revealed in nearby galaxy halos from studies of their resolved stellar content. We do not find support for the extremely red colours found in earlier studies of diffuse halo emission, which we suggest might have been due to residual systematic errors.

We present the results of the GREAT08 Challenge, a blind analysis challenge to infer weak gravitational lensing shear distortions from images. The primary goal was to stimulate new ideas by presenting the problem to researchers outside the shear measurement community. Six GREAT08 Team methods were presented at the launch of the Challenge and five additional groups submitted results during the 6 month competition. Participants analyzed 30 million simulated galaxies with a range in signal to noise ratio, point-spread function ellipticity, galaxy size, and galaxy type. The large quantity of simulations allowed shear measurement methods to be assessed at a level of accuracy suitable for currently planned future cosmic shear observations for the first time. Different methods perform well in different parts of simulation parameter space and come close to the target level of accuracy in several of these. A number of fresh ideas have emerged as a result of the Challenge including a re-examination of the process of combining information from different galaxies, which reduces the dependence on realistic galaxy modelling. The image simulations will become increasingly sophisticated in future GREAT challenges, meanwhile the GREAT08 simulations remain as a benchmark for additional developments in shear measurement algorithms.

We present the first VLT near-IR observations of a gravitationally lensed quasar, using adaptive optics and laser guide star. These observations can be considered as a test bench for future systematic observations of lensed quasars with adaptive optics, even when bright natural guide stars are not available in the nearby field. With only 14 minutes of observing time, we derived very accurate astrometry of the quasar images and of the lensing galaxy, with 0.05 \arcsec spatial resolution, comparable to the Hubble Space Telescope (HST). In combination with deep VLT optical spectra of the quasar images, we use our adaptive optics images to constrain simple models for the mass distribution of the lensing galaxy. The latter is almost circular and does not need any strong external shear to fit the data. The time delay predicted for SDSS0806+2006, assuming a singular isothermal ellipsoid model and the concordance cosmology, is Delta t \simeq 50 days. Our optical spectra indicate a flux ratio between the quasar images of A/B=1.3 in the continuum and A/B=2.2 in both the MgII and in the CIII] broad emission lines. This suggests that microlensing affects the continuum emission. However, the constant ratio between the two emission lines indicates that the broad emission line region is not microlensed. Finally, we see no evidence of reddening by dust in the lensing galaxy.

We present the continuation of our long-term spectroscopic monitoring of the gravitationally lensed quasar QSO 2237+0305. We investigate the chromatic variations observed in the UV/optical continuum of both quasar images A and B, and compare them with numerical simulations to infer the energy profile of the quasar accretion disk. Our procedure combines the microlensing ray-shooting technique with Bayesian analysis, and derives probability distributions for the source sizes as a function of wavelength. We find that the effective caustic crossing timescale is 4.0+/-1.0 months. Using a robust prior on the effective transverse velocity, we find that the source responsible for the UV/optical continuum has an energy profile well reproduced by a power-law R lambda^{zeta} with zeta=1.2+/-0.3, where R is the source size responsible for the emission at wavelength lambda. This is the first accurate, model-independent determination of the energy profile of a quasar accretion disk on such small scales.

We present the main results of the first long-term spectrophotometric monitoring of the ``Einstein cross'' Q2237+0305 and of the single-epoch spectra of the lensed quasar J1131-1231. From October 2004 to December 2006, we find that two prominent microlensing events affect images A & B in Q2237+0305 while images C & D remain grossly unaffected by microlensing on a time scale of a few months. Microlensing in A & B goes with chromatic variations of the quasar continuum. We observe stronger micro-amplification in the blue than in the red part of the spectrum, as expected for continuum emission arising from a standard accretion disk. Microlensing induced variations of the CIII] emission are observed both in the integrated line intensity and profile. Finally, we also find that images C & D are about 0.1-0.3 mag redder than images A & B. The spectra of images A-B-C in J1131-1231 reveal that, in April 2003, microlensing was at work in images A and C. We find that microlensing de-amplifies the continuum emission and the Broad Line Region (BLR) in these images. Contrary to the case of Q2237+0305, we do not find evidence for chromatic microlensing of the continuum emission. On the other hand, we observe that the Balmer and MgII broad line profiles are deformed by microlensing. These deformations imply an anti-correlation between the width of the emission line and the size of the corresponding emitting region. Finally, the differential microlensing of the FeII emission suggests that the bulk of FeII is emitted in the outer parts of the BLR while another fraction of FeII is produced in a compact region.

Gravitationally lensed quasars can be used to map the mass distribution in lensing galaxies and to estimate the Hubble constant H0 by measuring the time delays between the quasar images. Here we report the measurement of two independent time delays in the quadruply imaged quasar WFI J2033-4723 (z = 1.66). Our data consist of R-band images obtained with the Swiss 1.2 m EULER telescope located at La Silla and with the 1.3 m SMARTS telescope located at Cerro Tololo. The light curves have 218 independent epochs spanning 3 full years of monitoring between March 2004 and May 2007, with a mean temporal sampling of one observation every 4th day. We measure the time delays using three different techniques, and we obtain Dt(B-A) = 35.5 +- 1.4 days (3.8%) and Dt(B-C) = 62.6 +4.1/-2.3 days (+6.5%/-3.7%), where A is a composite of the close, merging image pair. After correcting for the time delays, we find R-band flux ratios of F_A/F_B = 2.88 +- 0.04, F_A/F_C = 3.38 +- 0.06, and F_A1/F_A2 = 1.37 +- 0.05 with no evidence for microlensing variability over a time scale of three years. However, these flux ratios do not agree with those measured in the quasar emission lines, suggesting that longer term microlensing is present. Our estimate of H0 agrees with the concordance value: non-parametric modeling of the lensing galaxy predicts H0 = 67 +13/-10 km s-1 Mpc-1, while the Single Isothermal Sphere model yields H0 = 63 +7/-3 km s-1 Mpc-1 (68% confidence level). More complex lens models using a composite de Vaucouleurs plus NFW galaxy mass profile show twisting of the mass isocontours in the lensing galaxy, as do the non-parametric models. As all models also require a significant external shear, this suggests that the lens is a member of the group of galaxies seen in field of view of WFI J2033-4723.

We present the results of the first long-term (2.2 years) spectroscopic monitoring of a gravitationally lensed quasar, namely the Einstein Cross Q2237+0305. We spatially deconvolve deep VLT/FORS1 spectra to accurately separate the spectrum of the lensing galaxy from the spectra of the quasar images. Accurate cross-calibration of the observations at 31 epochs from October 2004 to December 2006 is carried out using foreground stars observed simultaneously with the quasar. The quasar spectra are further decomposed into a continuum component and several broad emission lines. We find prominent microlensing events in the quasar images A and B, while images C and D are almost quiescent on a timescale of a few months. The strongest variations are observed in the continuum, and their amplitude is larger in the blue than in the red, consistent with microlensing of an accretion disk. Variations in the intensity and profile of the broad emission lines are also reported, most prominently in the wings of the CIII] and in the center of the CIV emission lines. During a strong microlensing episode observed in quasar image A, the broad component of the CIII] is more magnified than the narrow component. In addition, the emission lines with higher ionization potentials are more magnified than the lines with lower ionization potentials, consistent with the stratification of the broad line region (BLR) infered from reverberation mapping observations.

We conducted an exploratory search for quasars at z~ 6 - 8, using the Early Data Release from United Kingdom Infrared Deep Sky survey (UKIDSS) cross-matched to panoramic optical imagery. High redshift quasar candidates are chosen using multi-color selection in i,z,Y,J,H and K bands. After removal of apparent instrumental artifacts, our candidate list consisted of 34 objects. We further refined this list with deeper imaging in the optical for ten of our candidates. Twenty-five candidates were followed up spectroscopically in the near-infrared and in the optical. We confirmed twenty-five of our spectra as very low-mass main-sequence stars or brown dwarfs, which were indeed expected as the main contaminants of this exploratory search. The lack of quasar detection is not surprising: the estimated probability of finding a single z>6 quasar down to the limit of UKIDSS in the 27.3 square degrees of the EDR is <5%. We find that the most important limiting factor in this work is the depth of the available optical data. Experience gained in this pilot project can help refine high-redshift quasar selection criteria for subsequent UKIDSS data releases.

We report the first results of a spectroscopic search for Lyman alpha, envelopes around three z=4.5 radio-quiet quasars. Our observational strategy uses the FORS2 spectrograph attached to the UT1 of the Very Large Telescope (VLT) in the multi-slit mode. This allows us to observe simultaneously the quasars and several PSF stars. The spectra of the latter are used to remove the point-like quasar from the data, and to unveil the faint underlying Lyman alpha, envelopes associated with the quasars with unprecedented depth. We clearly detect an envelope around two of the three quasars. These envelopes measure respectively 10" and 13" in extent (i.e. 67 kpc and 87 kpc). This is 5 to 10 times larger than predicted by the models of Haiman & Rees (2001) and up to 100 times fainter. Our observations better agree with models involing a clumpy envelope as in Alam & Miralda-Escude (2002) or Chelouche et al. (2008). We find that the brighter quasars also have the brighter envelopes but that the extend of the envelopes does not depend on the quasar luminosity. Although our results are based on only two objects with a detected Lyman alpha, envelope, the quality of the spatial deblending of the spectra lends considerable hope to estimate the luminosity function and surface brightness profiles of high redshift Lyman alpha, envelopes down to F= 2-3 10^{-21} erg/s/cm^2/A. We find that the best strategy to carry out such a project is to obtain both narrow-band images and deep slit-spectra.

The host galaxies of six nearby QSOs are studied on the basis of high resolution HST optical images and spatially resolved VLT slit spectra. The gas ionization and velocity are mapped as a function of the distance to the central QSO. In the majority of the cases, the QSO significantly contributes to the gas ionization in its whole host galaxy, and sometimes even outside. Reflection or scattering of the QSO \ha line from remote regions of the galaxy is detected in several instances. The line shifts show that, in all cases, the matter responsible for the light reflection moves away from the QSO, likely accelerated by its radiation pressure. The two faintest QSOs reside in spirals, with some signs of a past gravitational perturbation. One of the intermediate luminosity QSOs resides in a massive elliptical containing gas ionized (and probably pushed away) by the QSO radiation. The other medium-power object is found in a spiral galaxy displaying complex velocity structure, with the central QSO moving with respect to the bulge, probably as a result of a galactic collision. The two most powerful objects are involved in violent gravitational interactions and one of them has no detected host. These results suggest that (1) large-scale phenomena, such as galactic collisions, are closely related to the triggering and the feeding of the QSO and (2) once ignited, the QSO has significant influence on its large-scale neighborhood (often the whole host and sometimes further away).

We present the results of the first long-term (2.2 years) spectroscopic monitoring of a gravitationally lensed quasar, namely the Einstein Cross Q2237+0305. The goal of this paper is to present the observational facts to be compared in follow-up papers with theoretical models to constrain the inner structure of the source quasar. We spatially deconvolve deep VLT/FORS1 spectra to accurately separate the spectrum of the lensing galaxy from the spectra of the quasar images. Accurate cross-calibration of the 58 observations at 31-epoch from October 2004 to December 2006 is carried out with non-variable foreground stars observed simultaneously with the quasar. The quasar spectra are further decomposed into a continuum component and several broad emission lines to infer the variations of these spectral components. We find prominent microlensing events in the quasar images A and B, while images C and D are almost quiescent on a timescale of a few months. The strongest variations are observed in the continuum of image A. Their amplitude is larger in the blue (0.7 mag) than in the red (0.5 mag), consistent with microlensing of an accretion disk. Variations in the intensity and profile of the broad emission lines are also reported, most prominently in the wings of the CIII] and center of the CIV emission lines. During a strong microlensing episode observed in June 2006 in quasar image A, the broad component of the CIII] is more highly magnified than the narrow component. In addition, the emission lines with higher ionization potentials are more magnified than the lines with lower ionization potentials, consistent with the results obtained with reverberation mapping. Finally, we find that the V-band differential extinction by the lens, between the quasar images, is in the range 0.1-0.3 mag.

Interest in the quasar HE0450-2958 arose following the publication of the non-detection of its expected massive host, leading to various interpretations. This article investigates the gaseous and stellar contents of the system through additional VLT/FORS slit spectra and integral field spectroscopy from VLT/VIMOS. We apply our MCS deconvolution algorithm on slit spectra for the separation of the QSO and diffuse components, and develop a new method to remove the point sources in Integral Field Spectra, allowing extraction of velocity maps, narrow-line images, spatially resolved spectra or ionization diagrams of the surroundings of HE0450-2958. The whole system is embedded in gas, mostly ionized by the QSO radiation field and shocks associated with radio jets. The observed gas and star dynamics are unrelated, revealing a strongly perturbed system. Despite longer spectroscopic observations, the host galaxy remains undetected.

We expand our Bayesian Monte Carlo method for analyzing the light curves of gravitationally lensed quasars to simultaneously estimate time delays and quasar structure including their mutual uncertainties. We apply the method to HE1104-1805 and QJ0158-4325, two doubly-imaged quasars with microlensing and intrinsic variability on comparable time scales. For HE1104-1805 the resulting time delay of (Delta t_AB) = t_A - t_B = 162.2 -5.9/+6.3 days and accretion disk size estimate of log(r_s/cm) = 15.7 -0.5/+0.4 at 0.2 micron in the rest frame are consistent with earlier estimates but suggest that existing methods for estimating time delays in the presence of microlensing underestimate the uncertainties. We are unable to measure a time delay for QJ0158-4325, but the accretion disk size is log(r_s/cm) = 14.9 +/- 0.3 at 0.3 micron in the rest frame.

We present the results of a continuing survey to detect Lya emitting galaxies at redshifts z~9: the ZEN ("z equals nine'') survey. We have obtained deep VLT/ISAAC observations in the narrow J-band filter NB119 directed towards three massive lensing clusters: Abell clusters 1689, 1835, and 114. The foreground clusters provide a magnified view of the distant universe and permit a sensitive test for the presence of very high-redshift galaxies. We search for z~9 Lya emitting galaxies displaying a significant narrow-band excess relative to accompanying J-band observations that remain undetected in HST/ACS optical images of each field. No sources consistent with this criterion are detected above the unlensed 90% point-source flux limit of the narrow-band image, F_NB=3.7e-18 ergs/s/cm2. To date, the total coverage of the ZEN survey has sampled a volume at z~9 of approximately 1700 co-moving Mpc3 to a Lya emission luminosity of 1e43 erg/s. We conclude by considering the prospects for detecting z~9 Lya emitting galaxies in light of both observed galaxy properties at z<7 and simulated populations at z>7.

We present the analysis of a sample of 20 bright low-redshift quasars (M_B<-23 and z < 0.35) observed spectroscopically with the VLT. The FORS1 spectra, obtained in Multi Object Spectroscopy (MOS) mode, allow to observe simultaneously the quasars and several reference stars used to spatially deconvolve the data. Applying the MCS deconvolution method, we are able to separate the individual spectra of the quasar and of the underlying host galaxy. Contrary to some previous claims, we find that luminous quasars are not exclusively hosted by massive ellipticals. Most quasar host galaxies harbour large amounts of gas, irrespective of their morphological type. Moreover, the stellar content of half of the hosts is a young Sc-like population, associated with a rather low metallicity interstellar medium. A significant fraction of the galaxies contain gas ionized at large distances by the quasar radiation. This large distance ionization is always associated with signs of gravitational interactions (as detected from images or disturbed rotation curves). The spectra of the quasars themselves provide evidence that gravitational interactions bring dust and gas in the immediate surrounding of the super massive black hole, allowing to feed it. The quasar activity might thus be triggered (1) in young gas-rich spiral galaxies by local events and (2) in more evolved galaxies by gravitational interactions or collisions. The latter mechanism gives rises to the most powerful quasars. Finally, we derive mass models for the isolated spiral host galaxies and we show that the most reliable estimators of the systemic redshift in the quasar spectrum are the tips of the Ha and Hb lines.

In order to study the host galaxies of type 1 (broad-line) quasars, we present a semi-analytic modelling method to decompose the on-nucleus spectra of quasars into nuclear and host galaxy channels. The method uses the spatial information contained in long-slit or slitlet spectra. A routine determines the best fitting combination of the spatial distribution of the point like nucleus and extended host galaxy. Inputs are a simultaneously observed PSF, and external constraints on galaxy morphology from imaging. We demonstrate the capabilities of the method to two samples of a total of 18 quasars observed with EFOSC at the ESO 3.6m telescope and FORS1 at the ESO VLT. ~50% of the host galaxies with sucessful decomposition show distortions in their rotation curves or peculiar gas velocities above normal maximum velocities for disks. This is consistent with the fraction from optical imaging. All host galaxies have quite young stellar populations, typically 1-2 Gyr. For the disk dominated hosts these are consistent with their inactive counterparts, the luminosity weighted stellar ages are much younger for the bulge dominated hosts, compared to inactive early type galaxies. While this presents further evidence for a connection of galaxy interaction and AGN activity for half of the sample, this is not clear for the other half: These are often undistorted disk dominated host galaxies, and interaction on a smaller level might be detected in deeper high-resolution images or deeper spectroscopic data. The velocity information does not show obvious signs for large scale outflows triggered by AGN feedback - the data is consistent with velocity fields created by galaxy interaction.

A high accuracy photometric reduction method is needed to take full advantage of the potential of the transit method for the detection and characterization of exoplanets, especially in deep crowded fields. In this context, we present DECPHOT, a new deconvolution-based photometry algorithm able to deal with a very high level of crowding and large variations of seeing. It also increases the resolution of astronomical images, an important advantage for the discrimination of false positives in transit photometry.

Aims: The knowledge of the redshift of a lensing galaxy that produces multiple images of a background quasar is essential to any subsequent modeling, whether related to the determination of the Hubble constant H_0 or to the mass profile of the lensing galaxy. We present the results of our ongoing spectroscopic observations of gravitationally lensed quasars in order to measure the redshift of their lensing galaxies. We report on the determination of the lens redshift in seven gravitationally lensed systems. Methods: Our deep VLT/FORS1 spectra are spatially deconvolved in order to separate the spectrum of the lensing galaxies from the glare of the much brighter quasar images. Our observing strategy involves observations in Multi-Object-Spectroscopy (MOS) mode which allows the simultaneous observation of the target and of several crucial PSF and flux calibration stars. The advantage of this method over traditional long-slit observations is that it allows a much more reliable extraction and flux calibration of the spectra. Results: We obtain the first reliable spectra of the lensing galaxies in six lensed quasars: FBQ 0951+2635 (z=0.260), BRI 0952-0115 (z=0.632), HE 2149-2745 (z=0.603), Q 0142-100 (z=0.491), SDSS J0246-0825 (z=0.723), and SDSS J0806+2006 (z=0.573). The last 3 redshifts also correspond to the MgII doublet seen in absorption in the quasar spectra at the lens redshift. Our spectroscopic redshifts of HE 2149-2745 and FBQ 0951+2635 are higher than previously reported, which means that H_0 estimates from these two systems must be revised to higher values. Finally, we reanalyse our spectra of Q 1355-2257 and find MgII in absorption at z=0.702, confirming our previous redshift estimate. The spectra of all lenses are typical of early-type galaxies.

Aims: Our aim is to measure the time delay between the two gravitationally lensed images of the z = 1.547 quasar SDSS J1650+4251, in order to estimate the Hubble constant H_0. Methods: Our measurement is based on R-band light curves with 57 epochs obtained at Maidanak Observatory, in Uzbekistan, from May 2004 to September 2005. The photometry is performed using simultaneous deconvolution of the data, which provides the individual light curves of the otherwise blended quasar images. The time delay is determined from the light curves using two very different numerical techniques, i.e., polynomial fitting and direct cross-correlation. The time delay is converted into H_0 following analytical modeling of the potential well. Results: Our best estimate of the time delay is Dt = 49.5 +/- 1.9 days, i.e., we reach a 3.8% accuracy. The R-band flux ratio between the quasar images, corrected for the time delay and for slow microlensing, is F_A /F_B = 6.2 +/- 5%. Conclusions: The accuracy reached on the time delay allows us to discriminate well between families of lens models. As for most other multiply imaged quasars, only models of the lensing galaxy that have a de Vaucouleurs mass profile plus external shear give a Hubble constant compatible with the current most popular value (H_0 = 72 +/- 8 km s-1 Mpc-1). A more realistic singular isothermal sphere model plus external shear gives H_0 = 51.7 +4.0 -3.0 km s-1 Mpc-1.

We report the discovery of the first known probable case of a physical triple quasar (not a gravitational lens). A previously known double system, QQ 1429-008 at z = 2.076, is shown to contain a third, fainter QSO component at the same redshift within the measurement errors. Deep optical and IR imaging at the Keck and VLT telescopes has failed to reveal a plausible lensing galaxy group or a cluster, and moreover, we are unable to construct any viable lensing model which could lead to the observed distribution of source positions and relative intensities of the three QSO image components. Furthermore, there are hints of differences in broad-band spectral energy distributions of different components, which are more naturally understood if they are physically distinct AGN. Therefore, we conclude that this system is most likely a physical triple quasar, the first such close QSO grouping known at any redshift. The projected component separations in the restframe are ~ 30 - 50 kpc for the standard concordance cosmology, typical of interacting galaxy systems. The existence of this highly unusual system supports the standard picture in which galaxy interactions lead to the onset of QSO activity.

A new method is presented for determining the Point Spread Function (PSF) of images that lack bright and isolated stars. It is based on the same principles as the MCS (Magain, Courbin, Sohy, 1998) image deconvolution algorithm. It uses the information contained in all stellar images to achieve the double task of reconstructing the PSFs for single or multiple exposures of the same field and to extract the photometry of all point sources in the field of view. The use of the full information available allows to construct an accurate PSF. The possibility to simultaneously consider several exposures makes it very well suited to the measurement of the light curves of blended point sources from data that would be very difficult or even impossible to analyse with traditional PSF fitting techniques. The potential of the method for the analysis of ground-based and space-based data is tested on artificial images and illustrated by several examples, including HST/NICMOS images of a lensed quasar and VLT/ISAAC images of a faint blended Mira star in the halo of the giant elliptical galaxy NGC5128 (Cen A).

A high accuracy photometry algorithm is needed to take full advantage of the potential of the transit method for the characterization of exoplanets, especially in deep crowded fields. It has to reduce to the lowest possible level the negative influence of systematic effects on the photometric accuracy. It should also be able to cope with a high level of crowding and with large scale variations of the spatial resolution from one image to another. A recent deconvolution-based photometry algorithm fulfills all these requirements, and it also increases the resolution of astronomical images, which is an important advantage for the detection of blends and the discrimination of false positives in transit photometry. We made some changes to this algorithm in order to optimize it for transit photometry and used it to reduce NTT/SUSI2 observations of two transits of OGLE-TR-113b. This reduction has led to two very high precision transit light curves with a low level of systematic residuals, used together with former photometric and spectroscopic measurements to derive new stellar and planetary parameters in excellent agreement with previous ones, but significantly more precise.

(Abridged) We present our VLT/FORS1 deep spectroscopic observations of the gravitationally lensed quasar SDSS J0924+0219, as well as archival HST/NICMOS and ACS images of the same object. The two-epoch spectra, obtained in the Multi Object Spectroscopy (MOS) mode, allow for very accurate flux calibration, spatial deconvolution of the data, and provide the redshift of the lensing galaxy z=0.394 +/- 0.001. These spectra, taken 15 days apart, show only slight continuum variations, while the broad emission lines display obvious changes in the red wing of the Mg II line, in the Fe II bands, and in the central part of the C III] line. Even though variations in the line profiles are present, we do not see any significant differences between the continuum and emission line flux ratios of images A and B of the quasar. Spatial deconvolution of the HST images reveals a double Einstein ring. One ring corresponds to the lensed quasar host galaxy at z=1.524 and a second bluer one, is the image either of a star-forming region in the host galaxy, or of another unrelated lower redshift object. We find that a broad range of lens models gives a satisfactory fit to the data. However, they predict very different time delays, making SDSS J0924+0219 an object of particular interest for photometric monitoring. In addition, the lens models reconstructed using exclusively the constraints from the Einstein rings, or using exclusively the astrometry of the quasar images, are not compatible. This suggests that substructures play an important role in SDSS J0924+0219.

Aims: To predict time delays for a sample of gravitationally lensed quasars and to evaluate the accuracy that can be realistically achieved on the value of H_0. Methods: We consider 14 lensed quasars that are candidates for time-delay monitoring and model them in detail using pixelized lens models. For each system, we provide a mass map, arrival-time surface and the distribution of predicted time-delays in a concordance cosmology, assuming H_0^{-1}=14 Gyr (H_0=70 in local units). Based on the predicted time-delays and on the observational circumstances, we rate each lens as `excellent' or `good' or `unpromising' for time-delay monitoring. Finally, we analyze simulated time delays for the 11 lens rated excellent or good, and show that H_0 can be recovered to a precision of 5%. Results: In combination with COSMOGRAIL paper I on the temporal sampling of lensed quasar light curves, the present work will help design monitoring campaigns of lensed quasars.

Aims: We measure the redshift of the lensing galaxy in eight gravitationally lensed quasars in view of determining the Hubble parameter H_0 from the time delay method. Methods: Deep VLT/FORS1 spectra of lensed quasars are spatially deconvolved in order to separate the spectrum of the lensing galaxies from the glare of the much brighter quasar images. A new observing strategy is devised. It involves observations in Multi-Object-Spectroscopy (MOS) which allows the simultaneous observation of the target and of several PSF and flux calibration stars. The advantage of this method over traditional long-slit observations is a much more reliable extraction and flux calibration of the spectra. Results: For the first time we measure the redshift of the lensing galaxy in three multiply-imaged quasars: SDSS J1138+0314 (z=0.445), SDSS J1226-0006 (z=0.517), SDSS J1335+0118 (z=0.440), and we give a tentative estimate of the redshift of the lensing galaxy in Q 1355-2257 (z=0.701). We confirm four previously measured redshifts: HE 0047-1756 (z=0.407), HE 0230-2130 (z=0.523), HE 0435-1223 (z=0.454) and WFI J2033-4723 (z=0.661). In addition, we determine the redshift of the second lensing galaxy in HE 0230-2130 (z=0.526). The spectra of all lens galaxies are typical for early-type galaxies, except for the second lensing galaxy in HE 0230-2130 which displays prominent [OII] emission.

We present the ZEN (z equals nine) survey: a deep, narrow J-band search for proto-galactic Lya emission at redshifts z=9. In the first phase of the survey, dubbed ZEN1, we combine an exceptionally deep image of the Hubble Deep Field South, obtained using a narrow band filter centred on the wavelength 1.187 microns, with existing deep, broad band images covering optical to near infrared wavelengths. Candidate z=9 Lya-emitting galaxies display a significant narrow band excess relative to the Js-band that are undetected at optical wavelengths. We detect no sources consistent with this criterion to the 90% point source flux limit of the NB image, F_NB = 3.28e-18 ergs/s/cm2. The survey selection function indicates that we have sampled a volume of approximately 340 h^{-3} Mpc3 to a Lya emission luminosity of 10e43 h^{-2} ergs/s. When compared to the predicted properties of z=9 galaxies based upon no evolution of observed z=6 Lya-emitting galaxies, the `volume shortfall' of the current survey, i.e. the volume required to detect this putative population, is a factor of at least 8 to 10. We also discuss continuing narrow J-band imaging surveys that will reduce the volume shortfall factor to the point where the no-evolution prediction from z=6 is probed in a meaningful manner.

Quasars are thought to be powered by the infall of matter onto a supermassive black hole at the centre of massive galaxies. As the optical luminosity of quasars exceeds that of their host galaxy, disentangling the two components can be difficult. This led in the 1990's to the controversial claim of the discovery of 'naked' quasars. Since then, the connection between quasars and galaxies has been well established. Here we report on the observation of a quasar lying at the edge of a gas cloud, whose size is comparable to that of a small galaxy, but whose spectrum shows no evidence for stars. The gas cloud is excited by the quasar itself. If a host galaxy is present, it is at least six times fainter than would normally be expected for such a bright quasar. The quasar is interacting dynamically with a neighbouring galaxy - which matter might be feeding the black hole.

We analyse the respective benefits and drawbacks of ground-based and space-based transit surveys for extrasolar planets. Based on simple but realistic assumptions about the fraction of lower main sequence stars harboring telluric and giant planets within the outer limit of the habitable zone, we predict the harvests of fictitious surveys with three existing wide field optical and near-IR cameras: the CFHT-Megacam, SUBARU-Suprime and VISTA-IR. An additional promising instrument is considered, VISTA-Vis, currently under development. The results are compared with the harvests predicted under exactly the same assumptions, for the space missions COROT and KEPLER. We show that ground-based wide field surveys may discover more giant planets than space missions. However, space surveys seem to constitute the best strategy to search for telluric planets. In this respect, the KEPLER mission appears 50 times more efficient than any of the ground--based surveys considered here. KEPLER might even discover telluric planets in the habitable zone of their host star.

We present the results of two-band CCD photometric monitoring of the gravitationally lensed quasar Q 0142-100 (UM 673).The data, obtained at ESO-La Silla with the 1.54 m Danish telescope in the Gunn i-band (October 1998 - September 1999) and in the Johnson V-band (October 1998 to December 2001), were analyzed using three different photometric methods. The light-curves obtained with all methods show variations, with a peak-to-peak amplitude of 0.14 magnitude in $V$. Although it was not possible to measure the time delay between the two lensed QSO images, the brighter component displays possible evidence for microlensing: it becomes bluer as it gets brighter, as expected under the assumption of differential magnification of a quasar accretion disk

We report on new high-resolution near-infrared images of the gravitationally lensed radio source PKS1830-211, a quasar at z=2.507. These adaptive optics observations, taken with the Very Large Telescope (VLT), are further improved through image deconvolution. They confirm the presence of a second object along the line of sight to the quasar, in addition to the previously known spiral galaxy. This additional object is clearly extended in our images. However, its faint luminosity does not allow to infer any photometric redshift. If this galaxy is located in the foreground of PKS1830-211, it complicates the modeling of this system and decreases the interest in using PKS1830-211 as a means to determine H0 via the time delay between the two lensed images of the quasar.

We use numerical simulations to test a broad range of plausible observational strategies designed to measure the time delay between the images of gravitationally lensed quasars. Artificial quasar light curves are created along with Monte-Carlo simulations in order to determine the best temporal sampling to adopt when monitoring the photometric variations of systems with time delays between 5 and 120 days, i.e., always shorter than the visibility window across the year. Few and realistic assumptions are necessary on the quasar photometric variations (peak-to-peak amplitude and time-scale of the variations) and on the accuracy of the individual photometric points. The output of the simulations is the (statistical) relative error made on the time delay measurement, as a function of 1- the object visibility over the year, 2- the temporal sampling of the light curves and 3- the time delay. Also investigated is the effect of long term microlensing variations which must be below the 5 % level (either intrinsically or by subtraction) if the goal is to measure time delays with an accuracy of 1-2 %. However, while microlensing increases the random error on the time delay, it does not significantly increase the systematic error, which is always a factor 5 to 10 smaller than the random error. Finally, it is shown that, when the time delay is comparable to the visibility window of the object, a logarithmic sampling can significantly improve the time delay determination. All results are presented in the form of compact plots to be used to optimize the observational strategy of future monitoring programs.

We present a deep, narrow J-band search for proto--galactic Lyman Alpha emission at redshifts z~9. We combine an exceptionally deep image of the Hubble Deep Field South, obtained using a narrow band filter centred on the wavelength 1.187 microns using the VLT/ISAAC facility, with existing deep, broad band images covering optical to near infrared wavelengths. We search for z~9 Lyman Alpha emitting galaxies displaying a significant narrow band excess relative to the Js-band that are undetected at optical wavelengths. We detect no sources consistent with this criterion to the 90% point source flux limit of the NB image, F_{NB} = 3.28 x 10^{-18} ergs/s/cm2. The survey selection function indicates that we have sampled a volume of approximately 340 h^{-3} Mpc^3 to a Lyman Alpha emission luminosity of 10^{43} h^{-2} ergs/s. We conclude by considering the potential implications for the physics of the high--redshift universe.

We present optical R-band light curves of the double gravitationally lensed quasar FBQ 0951+2635 from observations obtained at the Nordic Optical Telescope between March 1999 and June 2001. A time delay of 16 +/- 2 days (1 sigma) is determined from the light curves. New constraints on the lensing geometry are provided by the position and ellipticity of the lensing galaxy. For a (Omega_m, Omega_Lambda) = (0.3, 0.7) cosmology, the time delay yields a Hubble parameter of H_0 = 60^{+9}_{-7} (random, 1 sigma) +/- 2 (systematic) km/s/Mpc for a singular isothermal ellipsoid model and H_0 = 63^{+9}_{-7} (random, 1 sigma) +/- 1 (systematic) km/s/Mpc for a constant mass-to-light ratio model. In both models, the errors are mainly due to the time-delay uncertainties. Non-parametric models yield H_0 = 64^{+9}_{-7} (random, 1 sigma) +/- 14 (systematic) km/s/Mpc.

In this paper we present new ESO/VLT FORS1 and ISAAC images of the fields around eight gravitationally lensed quasars: CTQ414, HE0230-2130, LBQS1009-0252, B1030+074, HE1104-1805, B1359+154, H1413+117 and HE2149-2745. When available and deep enough, HST/WFPC2 data were also used to infer the photometric redshifts of the galaxies around the quasars. The search of galaxy overdensities in space and redshift, as well as a weak-shear analysis and a mass reconstruction are presented in this paper. We find that there are most probably galaxy groups towards CTQ414, HE0230-2130, B1359+154, H1413+117 and HE2149-2745, with a mass ~ 4x10^14 M_sol h^-1. Considering its photometric redshift, the galaxy group discovered in the field around HE1104-1805 is associated with the quasar rather than with the lensing potential.

VLT on-axis optical spectroscopy of the z=0.144 radio-loud quasar HE 1434-1600 is presented. The spatially resolved spectra of the host galaxy are deconvolved and separated from those of the central quasar in order to study the dynamics of the stars and gas as well as the physical conditions of the ISM. We find that the host of HE 1434-1600 is an elliptical galaxy that resides in a group of at least 5 member galaxies, and that most likely experienced a recent collision with its nearest companion. Compared with other quasar host galaxies, HE 1434-1600 has a highly ionized ISM. The ionization state corresponds to that of typical Seyferts, but the ionized regions are not distributed in a homogeneous way around the QSO, and are located preferentially several kiloparsecs away from it. While the stellar absorption lines do not show any significant velocity field, the gas emission lines do. The observed gas velocity field is hard to reconcile with dynamical models involving rotating disk, modified Hubble laws or power laws, that all require extreme central masses (M > 10^9 M_sun) to provide only poor fit to the data. Power law models, which best fit the data, provide a total mass of M(< 10 kpc)=9.2 10^10 M_sun. We conclude that the recent interaction between HE 1434-1600 and its closest companion has strongly affected the gas velocity and ionization state, from the center of the galaxy to its most external parts.

We report the discovery of a new quadruply imaged quasar surrounded by an optical Einstein ring candidate. Spectra of the different components of 1RXS J113155.4-123155 reveal a source at z=0.658. Up to now, this object is the closest known gravitationally lensed quasar. The lensing galaxy is clearly detected. Its redshift is measured to be z=0.295. Additionally, the total V magnitude of the system has varied by 0.3 mag between two epochs separated by 33 weeks. The measured relative astrometry of the lensed images is best fitted with an SIS model plus shear. This modeling suggests very high magnification of the source (up to 50 for the total magnification) and predicts flux ratios between the lensed images significantly different from what is actually observed. This suggests that the lensed images may be affected by a combination of micro or milli-lensing and dust extinction effects.

The determination of the Hubble parameter H_0 is probably one of the most important applications of quasar lensing. The method, based on the measurement of the so-called ``time-delay'' between the lensed images of distant sources, e.g., quasars, and on detailed mass modeling of the potential well responsible for the multiple images, yields an accuracy at least comparable with other techniques and that can be improved further with high precision observations, as can be obtained with intrumentation of constantly increasing quality. The basics of the ``time-delay'' method are described, and the emphasis is put on the observational constraints available to the astrophysicist in order to implement the method and to derive an accurate value for H_0, independent of any standard candle or any strong prior on the other cosmological parameters.

Very deep ESO/VLT FORS1 and ISAAC images, as well as HST NICMOS2 data are used to infer the nature of the quasar pair LBQS 1429-0053 A and B, either a binary quasar or a doubly-imaged lensed quasar. Direct search of a putative lensing galaxy is unsuccessful . Moreover, no galaxy overdensity close to the quasar pair is found. A weak shear analysis of the field also fails at detecting any concentration of dark matter and weakens the hypothesis of a dark lens. The only sign of a possible lens consists in a group of 5 objects at z~1, within a radius of 5'', from the quasar pair. Considering this group as the lensing potential does not allow to reproduce the image position and flux ratio of the quasars. Our deep R-band image shows a blue, extended object at the position of quasar A, which is consistent with either being the lensed quasar A host, or being an intervening galaxy at lower redshift. Unless future very deep optical images demonstrate that this object is actually the lensed host of LBQS 1429-0053, we conclude that there is very little evidence for this quasar being lensed. Therefore, we are led to declare LBQS 1429-0053 A and B a genuine binary quasar.

We present the first result of a comprehensive spectroscopic study of quasar host galaxies. On-axis, spatially resolved spectra of low redshift quasars have been obtained with FORS1, mounted on the 8.2m ESO Very Large Telescope, Antu. The spectra are spatially deconvolved using a spectroscopic version of the ``MCS deconvolution algorithm''. The algorithm decomposes two dimensional spectra into the individual spectra of the central point-like nucleus and of its host galaxy. Applied to HE1503+0228, at z=0.135 (M_B=-23.0), it provides us with the spectrum of the host galaxy between 3600A, and 8500A, (rest-frame), at a mean resolving power of 700. The data allow us to measure several of the important Lick indices. The stellar populations and gas ionization state of the host galaxy of HE1503+0228, are very similar to the ones measured for normal non-AGN galaxies. Dynamical information is also available for the gas and stellar components of the galaxy. Using deconvolution and a deprojection algorithm, velocity curves are derived for emission lines, from the center up to 4arcsec, away from the nucleus of the galaxy. Fitting a simple three-components mass model (point mass, spherical halo of dark matter, disk) to the position-velocity diagram, we infer a mass of M(r<1kpc) = (2.0 +- 0.3) 10^{10} solar mass within the central kiloparsec of the galaxy, and a mass integrated over 10 kpc of M(r<10kpc) = (1.9 +- 0.3) 10^{11} solar mass, with an additional 10% error due to the uncertainty on the inclination of the galaxy. This, in combination with the analysis of the stellar populations indicates that the host galaxy of HE1503+0228, is a normal spiral galaxy.

Massive structures, such as galaxies, act as strong gravitational lenses on background sources. When the background source is a quasar, several lensed images are seen, as magnified or de-magnified versions of the same object. The detailed study of the image configuration and the measurement of ``time-delays'' between the images yield estimates of the Hubble parameter H_0. We describe in a simple way the phenomenon of strong lensing and review recent progress made in the field, including microlensing by stars in the main lensing galaxy.

We present optical R-band light curves of the gravitationally lensed quasar SBS1520+530 derived from data obtained at the Nordic Optical Telescope. A time delay of 130+/-3 days (1 sigma) is determined from the light curves. In addition, spectra of SBS1520+530 obtained at the Keck Observatory are spatially deconvolved in order to extract the spectrum of the faint lensing galaxy, free of any contamination by the light from the bright quasar images. This spectrum indicates a lens redshift z=0.717, in agreement with one of the absorption systems found in the quasar spectra. The best mass model of the system includes a second nearby galaxy and a cluster of galaxies in addition to the main lensing galaxy. Adopting this model and an Omega=0.3, Lambda=0.7 cosmology, our time-delay measurement yields a Hubble constant of H_0=51+/- 9, km/s/Mpc (1 sigma error).

Optical and near-IR Hubble Space Telescope and Gemini-North adaptive optics images, further improved through deconvolution, are used to explore the gravitationally lensed radio source PKS 1830-211. The line of sight to the quasar at z=2.507 appears to be very busy, with the presence, within 0.5 arcsec from the source of: (i) a possible galactic main-sequence star, (ii) a faint red lensing galaxy visible only in H-band and (iii) a new object whose colors and morphology match those of an almost face-on spiral. The V-I color and faint I magnitude of the latter suggest that it is associated with the molecular absorber seen towards PKS 1830-211, at z=0.89 rather than with the z=0.19 HI absorber previously reported in the spectrum of PKS 1830-211. While this discovery might ease the interpretation of the observed absorption lines, it also complicates the modeling of the lensing potential well, hence decreasing the interest in using this system as a mean to measure H_0 through the time delay between the lensed images. This is the first case of a quasar lensed by an almost face-on spiral galaxy.

The gravitational potential responsible for the lensing effect in SBS 1520+530 is studied over length scales from a few arc-seconds to a few arc-minutes. For this purpose, we use sharply deconvolved Hubble Space Telescope images in the optical and near-IR, in combination with ground based optical data obtained over a wider field-of-view. In particular, we have carried out a multi-color analysis in order to identify groups or clusters of galaxies along the line of sight. Photometric redshifts are measured for 139 galaxies unveiling significant excesses of galaxies 1.0 arcmin, NW and 1.7 arcmin, SW of the main lensing galaxy. The photometric redshift inferred both for the main lensing galaxy and for the galaxy concentrations is z=0.9^{+0.10}_{-0.25}. This is in rough agreement with the measured spectroscopic redshift of the main lensing galaxy, z=0.71 (Burud et al. 2002), suggesting that it is part of a larger group or cluster. We investigate the impact of including the galaxy cluster, first on the modelling of the lensing system, and second on the expected time--delay between the two quasar images.

We present optical V and i-band light curves of the gravitationally lensed BAL quasar HE2149-2745. The data, obtained with the 1.5m Danish Telescope (ESO-La Silla) between October 1998 and December 2000, are the first from a long-term project aimed at monitoring selected lensed quasars in the Southern Hemisphere. A time delay of 103+/-12 days is determined from the light curves. In addition, VLT/FORS1 spectra of HE2149-2745 are deconvolved in order to obtain the spectrum of the faint lensing galaxy, free of any contamination by the bright nearby two quasar images. By cross-correlating the spectrum with galaxy-templates we obtain a tentative redshift estimate of z=0.495+/-0.01. Adopting this redshift, a Omega=0.3, Lambda=0.7 cosmology, and a chosen analytical lens model, our time-delay measurement yields a Hubble constant of H_0=66+/-8 km/s/Mpc with an estimated systematic error of +/-3 km/s/Mpc. Using non-parametric models yields H_0=65+/-8 km/s/Mpc and confirms that the lens exhibits a very dense/concentrated mass profile.

We used the direct CCD camera at the Magellan I telescope at Las Campanas Observatory and FORS1 at Antu VLT at ESO Paranal Observatory to image fields centered on the inner and outer optical filaments in the halo of NGC 5128. In the V vs. U-V color-magnitude diagrams we have identified young blue supergiants associated with these line-emitting filaments located between the inner radio lobe and the northern middle lobe. Around the outer filament stars as young as 10 Myr were detected. They are principally aligned with the direction of the radio jet, but a vertical north-east alignment along the edge of the HI cloud is also present. Young stars in the inner filament field are found inside the bright knots of photoionized gas and are strongly aligned in the direction towards the center of the galaxy at the same position angle as the inner radio jet. Fitting the Padova isochrones on UV color-magnitude diagrams we find that the blue stars around the inner filaments have ages similar to the ones around the outer filaments ~10-15 Myr and the same abundance of Z=0.004. The presence of young blue supergiants clearly shows that the bright blue knots in the north-eastern halo of NGC 5128 are associations of young stars with photoionized gas. The temperature of the brightest stars is T \~ 12000-16000 K, insufficient to account alone for the high excitation lines observed in the surrounding ionized gas. Thus the optical emission jet is principally seen due to its alignment with the radio structure of the AGN. The highly collimated star formation is present only in the north eastern halo of the galaxy, suggesting the interaction of the jet with the gas clouds deposited during the last accretion event as the preferred triggering mechanism. (abridged)

A large fraction of the present-day stellar mass was formed between z=0.5 and z~3 and our understanding of the formation mechanisms at work at these epochs requires both high spatial and high spectral resolution: one shall simultaneously} obtain images of objects with typical sizes as small as 1-2kpc(~0''.1), while achieving 20-50 km/s (R >= 5000) spectral resolution. The obvious instrumental solution to adopt in order to tackle the science goal is therefore a combination of multi-object 3D spectrograph with multi-conjugate adaptive optics in large fields. A partial, but still competitive correction shall be prefered, over a much wider field of view. This can be done by estimating the turbulent volume from sets of natural guide stars, by optimizing the correction to several and discrete small areas of few arcsec2 selected in a large field (Nasmyth field of 25 arcmin) and by correcting up to the 6th, and eventually, up to the 60th Zernike modes. Simulations on real extragalactic fields, show that for most sources (>80%), the recovered resolution could reach 0".15-0".25 in the J and H bands. Detection of point-like objects is improved by factors from 3 to >10, when compared with an instrument without adaptive correction. The proposed instrument concept, FALCON, is equiped with deployable mini-integral field units (IFUs), achieving spectral resolutions between R=5000 and 20000. Its multiplex capability, combined with high spatial and spectral resolution characteristics, is a natural ground based complement to the next generation of space telescopes.

We present the results of an extension of the FIRST Bright Quasar Survey (FBQS) to the South Galactic cap, and to a fainter optical magnitude limit. Radio source counterparts with SERC R magnitudes brighter than 18.9 which meet the other FBQS criteria are included. We supplement this list with a modest number of additional objects to test our completeness for quasars with extended radio morphologies. The survey covers 589 square degrees in two equatorial strips in the southern cap. We have obtained spectra for 86% of the 522 candidates, and find 321 radio-selected quasars of which 264 are reported here for the first time. A comparison of this fainter sample with the FBQS sample shows the two to be generally similar. Fourteen new broad absorption line (BAL) quasars are included in this sample. When combined with the previously identified BAL quasars in our earlier papers, we can discern a break in the frequency of BAL quasars with radio loudness, namely that the relative number of high-ionization BAL quasars drops by a factor of four for quasars with a radio-loudness parameter R* > 100.

Using FORS2, mounted on Kueyen (UT2 of the VLT), we have obtained the redshift of the lensing galaxy in the gravitational lens system HE1104-1805. We measure z=0.729, in good agreement with previous estimates based on the time delay and the position of the lens on the fundamental plane. It also coincides with one of the metal line absorption systems that are seen in high resolution spectra of HE1104-1805.

(Abridged) A new technique for the spatial deconvolution of spectra is applied to near-IR (0.95 - 2.50 micron) NTT/SOFI spectra of the lensed, radio-quiet quasar HE 1104-1805. The continuum of the lensing galaxy is revealed between 1.5 and 2.5 micron. It is used in combination with previous optical and IR photometry to infer a plausible redshift in the range 0.8 < z < 1.2. Modeling of the system shows that the lens is probably composed of the red galaxy seen between the quasar images and a more extended component associated with a galaxy cluster with fairly low velocity dispersion (~ 575 km/s). The spectra of the two lensed images of the source show no trace of reddening at the redshift of the lens nor at the redshift of the source. Additionally, the difference between the spectrum of the brightest component a nd that of a scaled version of the faintest component is a featureless continuum. Broad and narrow emission lines, including the FeII features, are perfectly subtracted. The very good quality of our spectrum makes it possible to fit precisely the optical Fe II feature, taking into account the underlying continuum over a wide wavelength range. HE 1104-1805 can be classified as a weak Fe II emitter. Finally, the slope of the continuum in the brightest image is steeper than the continuum in the faintest image and supports the finding by Wisotzki et al. (1993) that the brightest image is microlensed. This is particularly interesting in view of the new source reconstruction methods from multiwavelength photometric monitoring.

We have applied the MCS image deconvolution algorithm (Magain, Courbin & Sohy 1998) to HST/WFPC2 V, I data of three M31 bulge globular clusters (G170, G177, and G198) and control fields near each cluster. All three clusters are clearly detected, with an increase in stellar density with decreasing radius from the cluster centers; this is the first time that stars have been resolved in bulge clusters in the inner regions of another galaxy. From the RGB slopes of the clusters and the difference in I magnitude between the HB and the top of the RGB, we conclude that these three clusters all have roughly solar metallicity, in agreement with earlier integrated-light spectroscopic measurements. Our data support a picture whereby the M31 bulge clusters and field stars were born from the same metal-rich gas, early in the galaxy formation.

We present deep HST imaging, as well as ground-based imaging and spectroscopy, of the optical afterglow associated with the long-duration gamma-ray burst GRB 990712 and its host galaxy. The data were obtained 48--123 days after the burst occurred. The magnitudes of the host (R = 21.9, V = 22.5) and optical afterglow (R = 25.4, V = 25.8, 47.7 days after the burst) favor a scenario where the optical light follows a pure power-law decay with an index of alpha ~ -1.0. We find no evidence for a contribution from a supernova like SN1998bw. This suggests that either there are multiple classes of long-duration gamma-ray bursts, or that the peak luminosity of the supernova was > 1.5 mag fainter than SN1998bw. The HST images and EFOSC2 spectra indicate that the gamma-ray burst was located in a bright, extended feature (possibly a star-forming region) 1.4 kpc from the nucleus of a 0.2 L^*_B galaxy at z = 0.434, possibly a Seyfert 2 galaxy. The late-time afterglow and host galaxy of GRB 990712 bear some resemblance to those of GRB 970508.

A method for spatial deconvolution of spectra is presented. It follows the same fundamental principles as the ``MCS image deconvolution algorithm'' (Magain, Courbin, Sohy, 1998) and uses information contained in the spectrum of a reference Point Spread Function (PSF) to spatially deconvolve spectra of very blended sources. An improved resolution rather than an infinite one is aimed at, overcoming the well known problem of ``deconvolution artefacts''. As in the MCS algorithm, the data are decomposed into a sum of analytical point sources and a numerically deconvolved background, so that the spectrum of extended sources in the immediate vicinity of bright point sources may be accurately extracted and sharpened. The algorithm has been tested on simulated data including seeing variation as a function of wavelength and atmospheric refraction. It is shown that the spectra of severely blended point sources can be resolved while fully preserving the spectrophotometric properties of the data. Extended objects ``hidden'' by bright point sources (up to 4-5 magnitudes brighter) can be accurately recovered as well, provided the data have a sufficiently high total signal-to-noise ratio (200-300 per spectral resolution element). Such spectra are relatively easy to obtain, even down to faint magnitudes, within a few hours of integration time with 10m class telescopes.

We report on the spectroscopic identification and the long awaited redshift measurement of the heavily obscured, gravitationally lensed radio source PKS 1830-211, which was first observed as a radio Einstein ring. The NE component of the doubly imaged core is identified, in our infrared spectrum covering the wavelength range 1.5-2.5 microns, as an impressively reddened quasar at z=2.507. Our redshift measurement, together with the recently measured time delay (Lovell et al.), means that we are a step closer to determining the Hubble constant from this lens. Converting the time delay into the Hubble constant by using existing models leads to high values for the Hubble constant. Since the lensing galaxy lies very close to the center of the lensed ring, improving the error bars on the Hubble constant will require not only a more precise time delay measurement, but also very precise astrometry of the whole system.

High resolution, ground based R and I band observations of the ultra compact dwarf galaxy POX 186 are presented. The data, obtained with the ESO New Technology Telescope (NTT), are analyzed using a new deconvolution algorithm which allows one to resolve the innermost regions of this stellar-like object into three Super-Star Clusters (SSC). Upper limits to both masses (M\sim 10^5 M_{\odot}) and the physical sizes (\le 60pc) of the SSCs are set. In addition, and maybe most importantly, extended light emission underlying the compact star-forming region is clearly detected in both bands. The R-I color rules out nebular H\alpha contamination and is consistent with an old stellar population. This casts doubt on the hypothesis that Blue Compact Dwarf Galaxies (BCDG) are young galaxies.

We detect a distant cluster of galaxies centered on the QSO lens and luminous X-ray source AXJ2019+112, a.k.a. ``The Dark Cluster'' (Hattori et al 1997). Using deep V,I Keck images and wide-field K_s imaging from the NTT, a tight red sequence of galaxies is identified within a radius of 0.2 h^{-1} Mpc of the known z=1.01 elliptical lensing galaxy. The sequence, which includes the central elliptical galaxy, has a slope in good agreement with the model predictions of Kodama et al (1998) for z~1. We estimate the integrated rest-frame luminosity of the cluster to be L_V > 3.2 x 10^{11}h^{-2}L_{\sun} (after accounting for significant extinction at the low latitude of this field), more than an order of magnitude higher than previous estimates. The central region of the cluster is deconvolved using the technique of Magain, Courbin & Sohy (1998), revealing a thick central arc coincident with an extended radio source. All the observed lensing features are readily explained by differential magnification of a radio loud AGN by a shallow elliptical potential. The QSO must lie just outside the diamond caustic, producing two images, and the arc is a highly magnified image formed from a region close to the center of the host galaxy, projecting inside the caustic. The mass--to--light ratio within an aperture of 0.4 h ^{-1} Mpc is M_x/L_V= 224^{+112}_{-78}h(M/L_V)_{\sun}, using the X-ray temperature. The strong lens model yields a compatible value, M/L_V= 372^{+94}_{-94}h(M/L_V)_{\sun}, whereas an independent weak lensing analysis sets an upper limit of M/L_V <520 h(M/L_V)_{\sun}, typical of massive clusters.

The Einstein Cross, Q~2237+0305, has been photometrically observed in four bands on two successive nights at NOT (La Palma, Spain) in October 1995. Three independent algorithms have been used to analyse the data: an automatic image decomposition technique, a CLEAN algorithm and the new MCS deconvolution code. The photometric and astrometric results obtained with the three methods are presented. No photometric variations were found in the four quasar images. Comparison of the photometry from the three techniques shows that both systematic and random errors affect each method. When the seeing is worse than 1.0", the errors from the automatic image decomposition technique and the Clean algorithm tend to be large (0.04-0.1 magnitudes) while the deconvolution code still gives accurate results (1{sigma} error below 0.04) even for frames with seeing as bad as 1.7". Reddening is observed in the quasar images and is found to be compatible with either extinction from the lensing galaxy or colour dependent microlensing. The photometric accuracy depends on the light distribution used to model the lensing galaxy. In particular, using a numerical galaxy model, as done with the MCS algorithm, makes the method less seeing dependent. Another advantage of using a numerical model is that eventual non-homogeneous structures in the galaxy can be modeled. Finally, we propose an observational strategy for a future photometric monitoring of the Einstein Cross.

We report the detection of four images in the recently discovered lensed QSO RX J0911.4+0551. With a maximum angular separation of 3.1", it is the quadruply imaged QSO with the widest known angular separation. Raw and deconvolved data reveal an elongated lens galaxy. The observed reddening in at least two of the four QSO images suggests differential extinction by this lensing galaxy. We show that both an ellipticity of the galaxy (epsilon_{min}=0.075) and an external shear (gamma_{min}=0.15) from a nearby mass has to be included in the lensing potential in order to reproduce the complex geometry observed in RX J0911.4+0551. A possible galaxy cluster is detected about 38", from RX J0911.4+0551 and could contribute to the X-ray emission observed by ROSAT in this field. The color of these galaxies indicates a plausible redshift in the range of 0.6-0.8.

New high quality Keck and ESO images of PKS 1830-211 are presented. Applying a powerful new deconvolution algorithm to these optical and infrared data, both images of the flat spectrum core of the radio source have been identified. An extended source is also detected in the optical images, consistent with the expected location of the lensing galaxy. The source counterparts are very red at I-K~7, suggesting strong Galactic absorption with additional absorption by the lensing galaxy at z=0.885, and consistent with the detection of high redshift molecules in the lens.

We discuss BIMA (Berkeley Illinois Maryland Association) data and present new high quality optical and near-IR Keck images of the bright radio ring PKS 1830-211. Applying a powerful new deconvolution algorithm we have been able to identify both images of the radio source. In addition we recover an extended source in the optical, consistent with the expected location of the lensing galaxy. The source counterparts are very red, I-K=7, suggesting strong Galactic absorption with additional absorption by the lensing galaxy at z=0.885, and consistent with the detection of high redshift molecules in the lens.

A new method for improving the resolution of astronomical images is presented. It is based on the principle that sampled data cannot be fully deconvolved without violating the sampling theorem. Thus, the sampled image should not be deconvolved by the total Point Spread Function, but by a narrower function chosen so that the resolution of the deconvolved image is compatible with the adopted sampling. Our deconvolution method gives results which are, in at least some cases, superior to those of other commonly used techniques: in particular, it does not produce ringing around point sources superimposed on a smooth background. Moreover, it allows to perform accurate astrometry and photometry of crowded fields. These improvements are a consequence of both the correct treatment of sampling and the recognition that the most probable astronomical image is not a flat one. The method is also well adapted to the optimal combination of different images of the same object, as can be obtained, e.g., from infrared observations or via adaptive optics techniques.

We report on deep IR imaging of the double quasar HE 1104-1805. A new image deconvolution technique has been applied to the data in order to optimally combine the numerous frames obtained. The resulting J and K' images allow us to detect and study the lensing galaxy between the two lensed QSO images. The near infrared images not only confirm the lensed nature of this double quasar, but also support the previous redshift estimate of z=1.66 for the lensing galaxy. No obvious overdensity of galaxies is detected in the immediate region surrounding the lens, down to limiting magnitudes of J=22 and K=20. The geometry of the system, together with the time delays expected for this lensed quasar, make HE 1104-1805 a remarkable target for future photometric monitoring programs, for the study of microlensing and for the determination of the cosmological parameters in the IR and optical domains.

We describe a first attempt to apply adaptive optics to the study of resolved stellar populations in galaxies. Advantages over traditional approaches are (i) improved spatial resolution and point-source sensitivity through adaptive optics, and (ii) use of the near-infrared region, where the peak of the spectral energy distribution for old populations is found. Disadvantages are the small area covered and the need for excellent seeing. We made observations with the ADONIS system at the European Southern Observatory of the peculiar elliptical galaxy NGC 5128; the irregular galaxy IC 5152 (a possible outer member of the Local Group); the Sc galaxy NGC 300 (a member of the Sculptor group); and the Sgr window in the bulge of the Milky Way. These different fields give excellent test cases for the potential of adaptive optics. In the first two cases, we failed to obtain photometry of individual stars, which would have required excellent seeing. For NGC 300 we measured magnitudes for nine individual supergiants (H = 18.3 to 20.2), but did not go deep enough to detect the tip of the RGB of an old population. For the Sgr field we produced a infrared luminosity function and colour-magnitude diagram for 70 stars down to about K = 19.5. These are the deepest yet measured for the Galactic bulge, reaching beyond the turn-off.

Time delay measurements have recently been reported for the lensed quasar PG 1115+080. These measurements can be used to derive Ho, but only if we can constrain the lensing potential. We have applied a recently developed deconvolution technique to analyse sub-arcsecond I band images of PG 1115+080, obtained at the Nordic Optical Telescope (NOT) and the Canada France Hawaii Telescope (CFHT). The high performance of the deconvolution code allows us to derive precise positions and magnitudes for the four lensed images of the quasar, as well as for the lensing galaxy. The new measurement of the galaxy position improves its precision by a factor of 3 and thus strengthens the constraints on the lensing potential. With the new data, a range of models incorporating some of the plausible systematic uncertainties yields Ho = 53 (+10/-7) km/s/mpc.

On the basis of sub-arcsecond optical imaging we resolve and study the core components of the LMC OB association LH 39. The central star of the association, the rare transition object R84, is also investigated using CASPEC echelle spectroscopy. A new, powerful image restoration code that conserves the fluxes allows us to obtain the magnitudes and colors of the components. We bring out some 30 stars in a ~16"x16", area centered on R84. At a resolution of 0.19", the closest components to R84 are shown to lie at 1.1" NW and 1.7" NW respectively of the transition star. One of these stars is the reddest star of the field after R84 but turns out to be too faint to correspond to the red M2 supergiant previously reported to contaminate the spectrum of R84. If the late-type spectrum is due to a line-of-sight supergiant with a luminosity comparable to R84, it should lie closer than 0.12" to R84. The transition star shows spectral variability between 1982 and 1991. We also note some slight radial velocity variations of the Of emission lines over timescales of several years. Furthermore, we discuss the apparent absence of O type stars in this association.

We report on imaging and spectroscopic observations of the radio-loud, highly polarized quasar PKS 1610-771 (z = 1.71). Our long-slit spectroscopy of the companion 4.55 arcseconds NW of the quasar confirms the stellar nature of this object, so ruling out the previously suspected gravitationally lensed nature of this system. PKS 1610-771 looks fuzzy on our sub-arcsecond R and I images and appears located in a rich environment of faint galaxies. Possible magnification, without image splitting of the quasar itself, by some of these maybe foreground galaxies cannot be excluded. The continuum fuzz (made of the closest two objects, viz. A and D) is elongated in a direction orthogonal to the E vector of the optical polarization, as in high-redshift radio-galaxies. The spectrum of PKS 1610-771 appears strongly curved, in a convex way, with a maximum of intensity at ~ 7,600 A (2,800 A rest frame), possibly indicating a strong ultraviolet absorption by dust.

We present the results of a 7 year long photometric monitoring of two components (A and B) of UM 425, thought to be images, separated by 6.5", of the same z = 1.47 quasar. These components have been imaged through an R filter in order to obtain their light curves. The photometry was obtained by simultaneously fitting a stellar two-dimensional profile on each component. The brightest image (component A, m(R) = 15.7) shows a slow and smooth increase in brightness of 0.2 magnitude in seven years, while the faintest one (component B, m(R) = 20.1) displays an outburst of 0.4 magnitude which lasts approximately two years. The variation of component B may be interpreted in two ways, assuming UM 425 is gravitationally lensed. If it is due to an intrinsic variation of the quasar, we derive a lower limit of 3 years on the time delay from the fact that it is not observed in component A. On the other hand, if it is a microlensing ``High Amplification Event'', we estimate the size of the source to be ~ 10^{-3} pc, in agreement with standard models of AGNs. These observations are consistent with the gravitational lens interpretation of the object. Furthermore, all the CCD frames obtained under the best seeing conditions have been co-added, in an attempt to detect the deflector. The final R image reveals a rich field of faint galaxies in the magnitude range m(R) 22-24. No obvious deflector, nor any system of arcs or arclets is detected, down to a limiting magnitude of m(R) 24.