The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

Weak lensing, which is the deflection of light by matter along the line of sight, has proven to be an efficient method for constraining models of structure formation and reveal the nature of dark energy. So far, most weak-lensing studies have focused on the shear field that can be measured directly from the ellipticity of background galaxies. However, within the context of forthcoming full-sky weak-lensing surveys such as Euclid, convergence maps (mass maps) offer an important advantage over shear fields in terms of cosmological exploitation. While it carry the same information, the lensing signal is more compressed in the convergence maps than in the shear field. This simplifies otherwise computationally expensive analyses, for instance, non-Gaussianity studies. However, the inversion of the non-local shear field requires accurate control of systematic effects caused by holes in the data field, field borders, shape noise, and the fact that the shear is not a direct observable (reduced shear). We present the two mass-inversion methods that are included in the official Euclid data-processing pipeline: the standard Kaiser & Squires method (KS), and a new mass-inversion method (KS+) that aims to reduce the information loss during the mass inversion. This new method is based on the KS method and includes corrections for mass-mapping systematic effects. The results of the KS+ method are compared to the original implementation of the KS method in its simplest form, using the Euclid Flagship mock galaxy catalogue. In particular, we estimate the quality of the reconstruction by comparing the two-point correlation functions and third- and fourth- order moments obtained from shear and convergence maps, and we analyse each systematic effect independently and simultaneously...

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.

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

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.

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

In time-delay cosmography, three of the key ingredients are 1) determining the velocity dispersion of the lensing galaxy, 2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and 3) estimating the external convergence $\kappa_\mathrm{ext}$ from less massive structures that are not included in the mass model. We present results on all three of these ingredients for two time-delay lensed quasar systems, DES J0408-5354 and WGD 2038-4008. We use the Gemini, Magellan and VLT telescopes to obtain spectra to both measure the stellar velocity dispersions of the main lensing galaxies and to identify the line-of-sight galaxies in these systems. Next, we identify 10 groups in DES J0408-5354 and 2 groups in WGD 2038-4008using a group-finding algorithm. We then identify the most significant galaxy and galaxy-group perturbers using the "flexion shift" criterion. We determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for both of these systems based on our spectroscopy and on the DES-only multiband wide-field observations. Using weighted galaxy counts, calibrated based on the Millennium Simulation, we find that DES J0408-5354 is located in a significantly underdense environment, leading to a tight (width $\sim3\%$), negative-value $\kappa_\mathrm{ext}$ distribution. On the other hand, WGD 2038-4008 is located in an environment of close to unit density, and its low source redshift results in a much tighter $\kappa_\mathrm{ext}$ of $\sim1\%$, as long as no external shear constraints are imposed.

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

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

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

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

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

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.

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

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.

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

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 use the Wilkinson Microwave Anisotropy Probe 7-year data (WMAP7) to further probe point source detection technique in the sky maps of the cosmic microwave background (CMB) radiation. The method by Tegmark et al. for foreground reduced maps and the Kolmogorov parameter as the descriptor are adopted for the analysis of WMAP satellite CMB temperature data. Part of the detected points coincide with point sources already revealed by other methods. However, we have also found 2 source candidates for which still no counterparts are known, and identified 7 point sources listed in Planck Early Release Compact Source Catalogue as high reliability sources.

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.

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.

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.

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.

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.

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 deep Hubble Space Telescope (HST) NICMOS 2 F160W band observations of the central 56*57" (14pc*14.25pc) region around R136 in the starburst cluster 30 Dor (NGC 2070) located in the Large Magellanic Cloud. Our aim is to derive the stellar Initial Mass Function (IMF) down to ~1 Msun in order to test whether the IMF in a massive metal-poor cluster is similar to that observed in nearby young clusters and the field in our Galaxy. We estimate the mean age of the cluster to be 3 Myr by combining our F160W photometry with previously obtained HST WFPC2 optical F555W and F814W band photometry and comparing the stellar locus in the color-magnitude diagram with main sequence and pre-main sequence isochrones. The color-magnitude diagrams show the presence of differential extinction and possibly an age spread of a few megayears. We convert the magnitudes into masses adopting both a single mean age of 3 Myr isochrone and a constant star formation history from 2 to 4 Myr. We derive the IMF after correcting for incompleteness due to crowding. The faintest stars detected have a mass of 0.5 Msun and the data are more than 50% complete outside a radius of 5 pc down to a mass limit of 1.1 Msun for 3 Myr old objects. We find an IMF of dN/dlog(M) M^(-1.20+-0.2) over the mass range 1.1--20 Msun only slightly shallower than a Salpeter IMF. In particular, we find no strong evidence for a flattening of the IMF down to 1.1 Msun at a distance of 5 pc from the center, in contrast to a flattening at 2 Msun at a radius of 2 pc, reported in a previous optical HST study. We examine several possible reasons for the different results. If the IMF determined here applies to the whole cluster, the cluster would be massive enough to remain bound and evolve into a relatively low-mass globular cluster.

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 a sample of $i_{775}$-dropout candidates identified in five Hubble Advanced Camera for Surveys fields centered on Sloan Digital Sky Survey QSOs at redshift $z\sim 6$. Our fields are as deep as the Great Observatory Origins Deep Survey (GOODS) ACS images which are used as a reference field sample. We find them to be overdense in two fields, underdense in two fields, and as dense as the average density of GOODS in one field. The two excess fields show significantly different color distributions from that of GOODS at the 99% confidence level, strengthening the idea that the excess objects are indeed associated with the QSO. The distribution of $i_{775}$-dropout counts in the five fields is broader than that derived from GOODS at the 80% to 96% confidence level, depending on which selection criteria were adopted to identify $i_{775}$-dropouts; its width cannot be explained by cosmic variance alone. Thus, QSOs seem to affect their environments in complex ways. We suggest the picture where the highest redshift QSOs are located in very massive overdensities and are therefore surrounded by an overdensity of lower mass halos. Radiative feedback by the QSO can in some cases prevent halos from becoming galaxies, thereby generating in extreme cases an underdensity of galaxies. The presence of both enhancement and suppression is compatible with the expected differences between lines of sight at the end of reionization as the presence of residual diffuse neutral hydrogen would provide young galaxies with shielding from the radiative effects of the QSO.

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 a high signal-to-noise spectrum of a bright galaxy at z = 4.9 in 14 h of integration on VLT FORS2. This galaxy is extremely bright, i_850 = 23.10 +/- 0.01, and is strongly-lensed by the foreground massive galaxy cluster Abell 1689 (z=0.18). Stellar continuum is seen longward of the Ly-alpha emission line at ~7100 \AA, while intergalactic H I produces strong absorption shortward of Ly-alpha. Two transmission spikes at ~6800 Angstroms (A) and ~7040 A are also visible, along with other structures at shorter wavelengths. Although fainter than a QSO, the absence of a strong central ultraviolet flux source in this star forming galaxy enables a measurement of the H I flux transmission in the intergalactic medium (IGM) in the vicinity of a high redshift object. We find that the effective H I optical depth of the IGM is remarkably high within a large 14 Mpc (physical) region surrounding the galaxy compared to that seen towards QSOs at similar redshifts. Evidently, this high-redshift galaxy is located in a region of space where the amount of H I is much larger than that seen at similar epochs in the diffuse IGM. We argue that observations of high-redshift galaxies like this one provide unique insights on the nascent stages of baryonic large-scale structures that evolve into the filamentary cosmic web of galaxies and clusters of galaxies observed in the present universe.

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 quantify the intrinsic width of the red giant branches of three massive globular clusters in M31 in a search for metallicity spreads within these objects. We present HST/ACS observations of three massive clusters in M31, G78, G213, and G280. A thorough description of the photometry extraction and calibration is presented. After derivation of the color-magnitude diagrams, we quantify the intrinsic width of the red giant branch of each cluster. This width translates into a metallicity dispersion that indicates a complex star formation history for this type of system. For G78, sigma_[Fe/H]}=0.86 \pm 0.37; for G213, 0.89 \pm 0.20; and for G280, 1.03 \pm 0.26. We find that the metallicity dispersion of the clusters does not scale with mean metallicity. We also find no trend with the cluster mass. We discuss some possible formation scenarios that would explain our results.

The relationship between galaxies and supermassive black holes (SMBH) found in their cores plays a key role in the formation and evolution of both of these major constituents of the universe, as well as the evolution of the intergalactic medium. Neither can be fully understood on their own, and studies of galaxy and SMBH co-formation and co-evolution are now among the central topics of research in cosmology. Yet the very origins, and the early growth phases of the SMBH are still not firmly established. We review our current understanding of the relevant processes and their astrophysical and cosmological context, with an emphasis on the observability of the SMBH growth mechanisms at high redshifts, and their leftover progeny at low redshifts.

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.

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 results of the HST Advanced Camera for Surveys spectroscopic ground-based redshift survey in the field of A1689. We measure 98 redshifts, increasing the number of spectroscopically confirmed objects by sixfold. We present two spectra from this catalog of the Sextet Arcs, images which arise from a strongly-lensed Lyman Break Galaxy (LBG) at a redshift of z=3.038. Gravitational lensing by the cluster magnifies its flux by a factor of ~16 and produces six separate images with a total r-band magnitude of r_625=21.7. The two spectra, each of which represents emission from different regions of the LBG, show H I and interstellar metal absorption lines at the systemic redshift. Significant variations are seen in Ly-alpha profile across a single galaxy, ranging from strong absorption to a combination of emission plus absorption. A spectrum of a third image close to the brightest arc shows Ly-alpha emission at the same redshift as the LBG, arising from either another spatially distinct region of the galaxy, or from a companion galaxy close to the LBG. Taken as a group, the Ly-alpha equivalent width in these three spectra decreases with increasing equivalent width of the strongest interstellar absorption lines. We discuss how these variations can be used to understand the physical conditions in the LBG. Intrinsically, this LBG is faint, ~0.1L*, and forming stars at a modest rate, ~4 solar masses per year. We also detect absorption line systems toward the Sextet Arcs at z=2.873 and z=2.534. The latter system is seen across two of our spectra.

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 investigate whether the brightest globular clusters (GCs) in the giant elliptical galaxies are similar to the less luminous GCs like those found in Local Group galaxies, we study the velocity dispersion and structural parameter correlations of a sample of bright GCs in the nearest gE galaxy NGC 5128. UVES echelle spectrograph on the ESO VLT, and EMMI on the ESO NTT were used to obtain high resolution spectra of bright GCs in NGC 5128. The velocity dispersions were obtained for all the targets. The structural parameters were either taken from the existing literature, or derived from our VLT FORS1 images using the ISHAPE software. The velocity dispersion and structural parameter measurements were used to obtain masses and M/L_V ratios of 22 clusters. The masses of the clusters in our sample range from M_vir=10^5-10^7 M_sun and the average M/L_V is 3+/-1. The three GCs harbouring X-ray point sources are the second, third and sixth most massive in our sample. The most massive cluster, HCH99-18, is also the brightest and the largest in size. It has the mass (M_vir=1.4x10^7 M_sun) an order of magnitude larger than the most massive clusters in the Local Group, and a high M/L_V ratio (4.7+/-1.2). We discuss briefly possible formation scenarios for this object. The correlations of structural parameters, velocity dispersion, masses and M/L_V for the bright GCs in NGC 5128 extend the properties established for the most massive Local Group clusters towards those characteristic of dE galaxy nuclei and Ultra Compact Dwarfs (UCDs). The detection of the mass-radius and the mass-M/L_V relations for the GCs with masses greater than ~2x10^6 M_sun provides the missing link between ``normal'' old globular clusters, young massive clusters, and evolved objects like UCDs. (Abridged)

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.

We acquired radial velocities of a significant number of globular clusters (GCs) on wide fields between galaxies in the nearby Fornax cluster of galaxies, in order to derive their velocity dispersion radial profile and to probe the dynamics of the cluster. We used FLAMES on the VLT to obtain accurate velocities for 149 GCs, within a ~500x150 kpc strip centered on NGC 1399, the Fornax central galaxy. These objects are at the very bright tail (M_V < -9.5) of the GC luminosity function, overlapping the so-called ``ultra-compact dwarfs'' magnitude range. Eight of the brightest FLAMES-confirmed members indeed show hints of resolution in the subarcsecond pre-imaging data we used for selecting the ~500 targets for FLAMES spectroscopy. Ignoring the GCs around galaxies by applying 3d_25 diameter masks, we find 61 GCs of 20.0 < V < 22.2 lying in the intra-cluster (IC) medium. The velocity dispersion of the population of ICGCs is 200 km/s at ~150 kpc from the central NGC 1399 and rises to nearly 400 km/s at 200 kpc, a value which compares with the velocity dispersion of the population of dwarf galaxies, thought to be infalling from the surroundings of the cluster.

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.

We have used HST imaging of the central regions (R<100 arcsec, about 5 core radii) of the globular cluster 47 Tucanae to derive proper motions and U- and V-band magnitudes for 14,366 cluster members. We also present a catalogue of astrometry and F475W photometry for nearly 130,000 stars in a rather larger central area. These data are made available in their entirety, in the form of downloadable electronic tables. We use them first to obtain a new estimate for the position of the cluster center and to define the stellar density profile into essentially zero radius. We then search in particular for any very fast-moving stars, such as might be expected to result from very close stellar encounters. Likely fewer than 0.1% (and no more than about 0.3%) of stars have total speeds above the nominal central escape velocity in 47 Tuc, and at lower speeds the velocity distribution is described very well by a regular King model. Considerations of only the proper-motion velocity dispersion then lead to a number of results: (1) Blue stragglers in the core of 47 Tuc have a velocity dispersion lower than that of the cluster giants by a factor of sqrt{2}. (2) The velocity distribution in the cluster center is essentially isotropic, as expected. (3) Using a sample of radial velocities for stars in the core, we estimate the distance to 47 Tuc: D = 4.0 +/- 0.35 kpc. And (4) we infer a 1-sigma upper limit of M<1000-1500 solar masses for any central, intermediate-mass black hole. We can neither confirm nor refute the hypothesis that 47 Tuc might lie on an extension of the M-sigma relation observed for supermassive black holes in galaxy bulges. [Abridged]

We present a high resolution (R ~ 43000) abundance analysis of a total of nine stars in three of the five globular clusters associated with the nearby Fornax dwarf spheroidal galaxy. These three clusters (1, 2 and 3) trace the oldest, most metal-poor stellar populations in Fornax. We determine abundances of O, Mg, Ca, Ti, Cr, Mn, Fe, Ni, Zn, Y, Ba, La, Nd and Eu in most of these stars, and for some stars also Mn and La. We demonstrate that classical indirect methods (isochrone fitting and integrated spectra) of metallicity determination lead to values of [Fe/H] which are 0.3 to 0.5 dex too high, and that this is primarily due to the underlying reference calibration typically used by these studies. We show that Cluster 1, with [Fe /H] = -2.5, now holds the record for the lowest metallicity globular cluster. We also measure an over-abundance of Eu in Cluster 3 stars that has only been previously detected in a subgroup of stars in M15. We find that the Fornax globular cluster properties are a global match to what is found in their Galactic counterparts; including deep mixing abundance patterns in two stars. We conclude that at the epoch of formation of globular clusters both the Milky Way and the Fornax dwarf spheroidal galaxy shared the same initial conditions, presumably pre-enriched by the same processes, with identical nucleosynthesis patterns.

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

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 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 tentative evidence for primordial clustering, manifested as an excess of color-selected objects in the field of the QSO SDSS J1030+0524 at redshift z=6.28. We have selected objects red in i_{775}-z_{850} on the basis of Hubble Space Telescope Advanced Camera for Surveys imaging of a field centered on the QSO. Compared to data at comparable depth obtained by the GOODS survey, we find an excess of objects with (i_{775}-z_{850}) \geq 1.5 in the QSO field. The significance of the detection is estimated to be ~97% on the basis of the counts alone and increases to 99.4% if one takes into account the color distribution. If confirmed this would represent the highest redshift example of galaxy clustering and would have implications on models for the growth of structure. Bias-driven clustering of first luminous objects forming in the highest peaks of the primordial density field is expected in most models of early structure formation. The redshift of one of the candidates has been found to be z=5.970 by our spectroscopy with Keck I/LRIS, confirming the validity of our color selection.

Since its launch in April 1990, the Hubble Space Telescope (HST) has produced an increasing flow of scientific results. The large number of refereed publications based on HST data allows a detailed evaluation of the effectiveness of this observatory and of its scientific programs. This paper presents the results of selected science metrics related to paper counts, citation counts, citation history, high-impact papers, and the most productive programs and most cited papers, through the end of 2003. All these indicators point towards the high-quality scientific impact of HST.

Recently, combining radial velocities from Keck/HIRES echelle spectra with published proper motion membership probabilities, Cote et al (2002) observed a sample of 21 stars, probable members of Palomar 13, a globular cluster in the Galactic halo. Their projected velocity dispersion sigma_p = 2.2 +/-0.4 km/s gives a mass-to-light ratio M/L_V = 40 +24/-17, about one order of magnitude larger than the usual estimate for globular clusters. We present here radial velocities measured from three different CCD frames of commissioning observations obtained with the new ESO/VLT instrument FLAMES (Fibre Large Array Multi Element Spectrograph). From these data, now publicly available, we measure the homogeneous radial velocities of eight probable members of this globular cluster. A new projected velocity dispersion sigma_p = 0.6-0.9 +/-0.3 km/s implies Palomar 13 mass-to-light ratio M/L_V = 3-7, similar to the usual value for globular clusters. We discuss briefly the two most obvious reasons for the previous unusual mass-to-light ratio finding: binaries, now clearly detected, and more homogeneous data from the multi-fibre FLAMES spectrograph.

We discuss the main properties of the Galactic globular cluster (GC) blue straggler stars (BSS), as inferred from our new catalog containing nearly 3000 BSS. The catalog has been extracted from the photometrically homogeneous V vs. (B-V) color-magnitude diagrams (CMD) of 56 GCs, based on WFPC2 images of their central cores. In our analysis we used consistent relative distances based on the same photometry and calibration. The number of BSS has been normalized to obtain relative frequencies (F_{BSS}) and specific densities (N_S) using different stellar populations extracted from the CMD. The cluster F_{BSS} is significantly smaller than the relative frequency of field BSS. We find a significant anti-correlation between the BSS relative frequency in a cluster and its total absolute luminosity (mass). There is no statistically significant trend between the BSS frequency and the expected collision rate. F_{BSS} does not depend on other cluster parameters, apart from a mild dependence on the central density. PCC clusters act like normal clusters as far as the BSS frequency is concerned. We also show that the BSS luminosity function for the most luminous clusters is significantly different, with a brighter peak and extending to brighter luminosities than in the less luminous clusters. These results imply that the efficiency of BSS production mechanisms and their relative importance vary with the cluster mass.

Correlations among the independently measured physical properties of globular clusters (GCs) can provide powerful tests for theoretical models and new insights into their dynamics, formation, and evolution. We review briefly some of the previous work, and present preliminary results from a comparative study of GC correlations in the Local Group galaxies. The results so far indicate that these diverse GC systems follow the same fundamental correlations, suggesting a commonality of formative and evolutionary processes which produce them.

Scarpa et al. (2003) have used their new measurements of the line-of-sight velocity dispersion at radii R ~ 30-45 pc in omega Centauri to argue for modifications to Newtonian dynamics (MOND) at low accelerations in the outer reaches of this globular cluster. Alternatively, any such evidence could be taken as an indication of the presence of an extended dark-matter halo. We show that the data of Scarpa et al. are in fact consistent with earlier observations by Seitzer and by Meylan et al., which are in turn well described at large radii by simple, self-consistent dynamical models. There is no need to invoke MOND or dark matter to explain the current data.

The distant globular cluster Palomar 13 has been found to have a very high mass-to-light ratio and its profile can be well fitted either by a King model with a tail, or with a NFW model. This cluster may be the first case of the many clumps predicted by CDM simulations that would not be disrupted by the galactic halo potential. We make the hypothesis that Pal 13 is made of neutralinos and run the DarkSuspect code to estimate the high-energy photon flux due to the annihilation of neutralinos through various channels in some benchmark scenarios. These low fluxes may be used as targets to be reached in proposals for future ground-based high altitude Cerenkov telescopes.

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.

We present the complete photometric database and the color-magnitude diagrams for 74 Galactic globular clusters observed with the HST/WFPC2 camera in the F439W and F555W bands. A detailed discussion of the various reduction steps is also presented, and of the procedures to transform instrumental magnitudes into both the HST F439W and F555W flight system and the standard Johnson B and V systems. We also describe the artificial star experiments which have been performed to derive the star count completeness in all the relevant branches of the color magnitude diagram. The entire photometric database and the completeness function will be made available on the Web immediately after the publication of the present paper.

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.

We have measured Keck/HIRES radial velocities for 30 candidate red giants in the direction of Palomar 13: an object traditionally cataloged as a compact, low-luminosity globular cluster. From a sample of 21 confirmed members, we find a systemic velocity of 24.1 km/s and a projected, intrinsic velocity dispersion of 2.2 km/s. Although small, this dispersion is several times larger than that expected for a globular cluster of this luminosity and central concentration. Taken at face value, this dispersion implies a mass-to-light ratio of ~ 40 (in solar units) based on the best-fit King-Michie model. The surface density profile of Palomar 13 also appears to be anomalous among Galactic globular clusters -- depending upon the details of background subtraction and model-fitting, Palomar 13 either contains a substantial population of "extra-tidal" stars, or it is far more spatially extended than previously suspected. The full surface density profile is equally well-fit by a King-Michie model having a high concentration and large tidal radius, or by a NFW model. We examine -- and tentatively reject -- a number of possible explanations for the observed characteristics of Palomar 13 (e.g., velocity "jitter" among the red giants, spectroscopic binary stars, non-standard mass functions, modified Newtonian dynamics), and conclude that the two most plausible scenarios are either catastrophic heating during a recent perigalacticon passage, or the presence of a massive dark halo. Thus, the available evidence suggests that Palomar 13 is either a globular cluster which is now in the process of dissolving into the Galactic halo, or a faint, dark-matter-dominated stellar system (ABRIDGED).

Using current observational data and simple dynamical modeling, we demonstrate that Omega Cen is not special among the Galactic globular clusters in its ability to produce and retain the heavy elements dispersed in the AGB phase of stellar evolution. The multiple stellar populations observed in Omega Cen cannot be explained if it had formed as an isolated star cluster. The formation within a progenitor galaxy of the Milky Way is more likely, although the unique properties of Omega Cen still remain a mystery.

We present the first detection of weak gravitational shear at infrared wavelengths, using observations of the lensing cluster Abell 1689, taken with the SofI camera on the ESO-NTT telescope. The imprint of cluster lenses on the shapes of the background galaxy population has previously been harnessed at optical wavelengths, and this gravitational shear signal enables cluster mass distributions to be probed, independent of whether the matter is luminous or dark. At near-infrared wavelengths, the spectrophotometric properties of galaxies facilitate a clean selection of background objects for use in the lensing analysis. A finite-field mass reconstruction and application of the aperture mass (Map) statistic are presented. The probability that the peak of the Map detection S/N~5, arises from a chance alignment of background sources is only ~4.5*10^-7. The velocity dispersion of the best-fit singular isothermal sphere model for the cluster is sigma_1D=1030^{+70}_{-80} km/s, and we find a K-band mass-to-light ratio of ~40 M_solar/L_solar inside a 0.44 Mpc radius.

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.

Mass determinations are difficult to obtain and still frequently characterised by deceptively large uncertainties. We review below the various mass estimators used for star clusters of all ages and luminosities. We highlight a few recent results related to (i) very massive old star clusters, (ii) the differences and similarities between star clusters and cores of dwarf elliptical galaxies, and (iii) the possible strong biases on mass determination induced by tidal effects.

(Abridged version) Mayall II = G1 is one of the brightest globular clusters belonging to M31, the Andromeda galaxy. Our observations with HST/WFPC2 provide data for the (I vs. V-I) and (V vs. V-I) color-magnitude diagrams. From model fitting, we determine a rather high mean metallicity of [Fe/H] = --0.95 +- 0.09, somewhat similar to 47 Tucanae. We find a larger spread in V-I than can be explained by the measurement errors, and we attribute this to an intrinsic metallicity dispersion amongst the stars of G1. So far, only omega Centauri, the giant Galactic globular cluster, has been known to exhibit such an intrinsic metallicity dispersion. Three estimates of the total mass of this globular cluster can be obtained: King mass = 15 x 10^6 with M/Lv ~ 7.5, Virial mass = 7.3 x 10^6 with M/Lv ~ 3.6, and King-Michie mass range from 14 to 17 x 10^6. Although uncertain, all of these mass estimates make G1 more than twice as massive as omega Centauri. Such large masses relate to the metallicity spread whose origin is still unknown (either self-enrichment, an inhomogeneous proto-cluster cloud, or remaining core of a dwarf galaxy). When considering the positions of G1 in the different diagrams defined by Kormendy (1985), G1 always appears on the sequence defined by globular clusters, and definitely away from the other sequences defined by elliptical galaxies, bulges, and dwarf spheroidal galaxies. The same is true for omega Centauri and for the nucleus of the dwarf elliptical NGC 205. This does not prove that all (massive) globular clusters are the remnant cores of nucleated dwarf galaxies.

We report on multi-object spectroscopy of 37 candidate star clusters around the giant early-type radio galaxy NGC 1316 (Fornax A), the brightest galaxy in the Fornax cluster. For the star clusters in this sample, we measure a mean radial velocity of 1698 +/- 46 km/s, a rotation velocity of 175 +/- 70 km/s, and a velocity dispersion of 227 +/- 33 km/s within a galactocentric radius of 24 kpc. Using the projected mass estimator and assuming isotropic orbits, the estimated total mass is (6.6 +/- 1.7) E+11 solar masses within a radius of 24 kpc. The mass is uncertain by about a factor of two, depending on the orbital assumptions. The implied M/L_B ratio is in the range 3-6. Four star clusters in our sample are exceptionally luminous (M_V < -12.3). This means that (1) at least this many clusters in NGC 1316 are up to an order of magnitude more luminous than the most luminous star cluster in our Galaxy or M31, and (2) that the S/N ratio of their spectra allows us to measure line strengths with good accuracy. By comparing the colours and equivalent widths of H-alpha and the Ca II triplet absorption lines for those bright star clusters in our sample with those of single-burst population models, we find that they are coeval with an age of 3.0 +/- 0.5 Gyr. Their metallicities are found to be solar to within 0.15 dex. These results show once again that globular clusters with near-solar metallicity do form during galactic mergers, and, moreover, that they can survive disruption processes taking place during the merger (e.g., dynamical friction, tidal disruption), as well as evaporation. In this respect, NGC 1316 provides a hitherto "missing" evolutionary link between young merger remnants of age ~0.5 Gyr such as NGC 3597, NGC 3921 and NGC 7252 on one side, and older giant ellipticals featuring bimodal colour distributions on the other side.

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

Large-field multi-color images of 20 galactic globular clusters are used to investigate the presence of tidal tails around these stellar systems. Field and cluster stars are sorted with the help of color-magnitude diagrams, and star-count analysis is performed on the selected cluster stars in order to increase the signal-to-noise ratio of their surface density. We study the overdensities of these stars using the wavelet transform of the star counts in order to filter the background density noise and to detect the weak structures, at large scale, formed by the numerous stars previously members of the clusters. We associate these stellar overdensities with the stars evaporated from the clusters because of dynamical relaxation and/or tidal stripping from the clusters by the galactic gravitational field. Most of the globular clusters in our sample display strong evidence of tidal interactions with the galactic plane in the form of large and extended deformations. These tidal tails exhibit projected directions preferentially towards the galactic center. All the clusters observed, which do not suffer from strong observational biases, present such tidal tails, tracing their dynamical evolution (evaporation, tidal shocking, tidal torquing, and bulge shocking) in the Galaxy. The clusters exhibit different regimes of mass loss rate, detected using the radial density slope in the outer parts of the clusters. For NGC 5139 ($\omega$ Centauri), we estimate, taking into account the possible presence of mass segregation in its outer parts, that about 0.6 to 1% of its mass has been lost during the crossing (Abridged)

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.

Star clusters - open and globulars - experience dynamical evolution on time scales shorter than their age. Consequently, open and globular clusters provide us with unique dynamical laboratories for learning about two-body relaxation, mass segregation from equipartition of energy, and core collapse. We review briefly, in the framework of star clusters, some elements related to the theoretical expectation of mass segregation, the results from N-body and other computer simulations, as well as the now substantial clear observational evidence.

Galactic globular clusters are ancient building blocks of our Galaxy. They represent a very interesting family of stellar systems in which some fundamental dynamical processes have been taking place for more than 10 Gyr, but on time scales shorter than the age of the universe. In contrast with galaxies, these star clusters represent unique laboratories for learning about two-body relaxation, mass segregation from equipartition of energy, stellar collisions, stellar mergers, core collapse, and tidal disruption. This review briefly summarizes some of the tremendous developments that have taken place during the last two decades. It ends with some recent results on tidal tails around galactic globular clusters and on a very massive globular cluster in M31.

We present our study of the 2-D structures of the tidal tails associated with 20 galactic globular clusters, structures obtained by using the wavelet transform to detect weak structures at large scale and filter the strong background noise for the low galactic latitude clusters (Leon, Meylan & Combes 2000). We also present N-body simulations of globular clusters, in orbits around the Galaxy, in order to study quantitatively and geometrically the tidal effects they encounter (Combes, Leon & Meylan 2000).

We present N-body simulations of globular clusters, in orbits around the Galaxy, in order to study quantitatively and geometrically the tidal effects they encounter. The clusters are modelised with multi-mass King-Michie models (Michie 1963), including mass segregation at initial conditions. The Galaxy is modelled as realistic as possible, with three components: bulge, disk and dark halo. The main finding is that there exist two giant tidal tails around the globuler cluster in permanence along its orbit, whatever this orbit. The length of these tails is of the order of 5 tidal radii, or greater. The escaped stars are distributed radially as a power law in density, with a slope of -4. The tails present substructures, or clumps, that are the relics of the strongest shocks. Due to the compressive disk-shocking, the clusters display a prolate shape which major axis is precessing around the z axis. The tails are preferentially formed by the lowest mass stars, as expected, so that the tidal truncation increases mass segregation. Internal rotation of the cluster increases the mass loss. The flattening of dark matter cannot influence significantly the dynamics of the clusters. The orientation and the strength of the tidal tails are signatures of the last disk crossing, so that observed tidal tails can constrain strongly the cluster orbit and the galactic model (vertical scale of the disc).

We present observations of the center of the Galactic globular cluster NGC 6273, obtained with the HST/WFPC2 camera as part of the snapshot program GO-7470. A B, V color-magnitude diagram (CMD) for ~28,000 stars is presented and discussed. The most prominent feature of the CMD, identified for the first time in this paper, is the extended horizontal branch blue tail (EBT) with a clear double-peaked distribution and a significant gap. The EBT of NGC 6273 is compared with the EBTs of seven other globular clusters for which we have a CMD in the same photometric system. From this comparison we conclude that all the globular clusters in our sample with an EBT show at least one gap along the HB, which could have similar origins. A comparison with theoretical models suggests that at least some of these gaps may be occuring at a particular value of the stellar mass, common to a number of different clusters. From the CMD of NGC 6273 we obtain a distance modulus (m-M)_V=16.27+-0.20. We also estimate an average reddening E(B-V)=0.47+-0.03, though the CMD is strongly affected by differential reddening, with the relative reddening spanning a 0.2 magnitude interval in E(B-V) in the WFPC2 field. A luminosity function for the evolved stars in NGC 6273 is also presented and compared with the most recent evolutionary models.

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.

We present observations with the Hubble Space Telescope Wide Field Planetary Camera 2 of three fields centered on super-metal-rich globular clusters in the bulge of M31. Our (I,V-I) color-magnitude diagrams reach as faint as I ~ 26.5 mag and clearly reveal the magnitude of the first ascent red giant branch (RGB) tip. We find that the apparent I magnitude of the RGB tip does not become brighter near the center of M31 as concluded by previous investigators. Our observations and artificial star experiments presented in this study strongly support the idea that previous very bright stars were likely the result of spurious detections of blended stars due to crowding in lower resolution images. On the contrary, our observations indicate that, at a mean projected galactocentric distance of 1.1 kpc, the RGB tip is some 1.3 magnitudes fainter than it is at 7 kpc. An analysis of this difference in RGB tip magnitude suggests that the M31 bulge stellar population has a mean metallicity close to that of the Sun.

We investigate the kinematics of a combined sample of 74 globular clusters around NGC 1399. Their high velocity dispersion, increasing with radius, supports their association with the gravitational potential of the galaxy cluster rather than with that of NGC 1399 itself. We find no evidence for rotation in the full sample, although some indication for rotation in the outer regions. The data do not allow us to detect differences between the kinematics of the blue and red sub-populations of globular clusters. A comparison between the globular cluster systems of NGC 1399 and those of NGC 1404 and NGC 1380 indicates that the globular clusters in all three galaxies are likely to have formed via similar mechanisms and at similar epochs. The only property which distinguishes the NGC 1399 globular cluster system from these others is that it is ten times more abundant. We summarize the evidence for associating these excess globulars with the galaxy cluster rather than with NGC 1399 itself, and suggest that the over-abundance can be explained by tidal stripping, at an early epoch, of neighboring galaxies and subsequent accumulation of globulars in the gravitational potential of the galaxy cluster.

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.

We report the discovery of two low redshift HI 21cm absorbers, one at z = 0.2212 towards the z_{em} = 0.630 quasar OI 363 (B0738+313), and the other at z = 0.3127 towards PKS B1127-145 (z_{em} = 1.187). Both were found during a survey of MgII selected systems at redshifts 0.2 < z < 1 using the new UHF-high system at the Westerbork Synthesis Radio Telescope (WSRT). New HST/FOS observations also identify both systems as damped Ly-alpha (DLa) absorbers. By comparing the column density from the DLa line with that from the HI 21cm line, we calculate the spin temperature, and find that T_s is approximately 1000 K for both of these low redshift absorbers. We briefly discuss some implications of this result.

We report on the radial velocity measurements of 18 globular clusters and one dE galaxy in the field of the giant elliptical galaxy NGC 1399, the dominant galaxy of the Fornax cluster. We present also the radial velocity measurements of some candidates for young and super metal-rich globular clusters, which turn out to be foreground or background objects, viz. 28 stars and 7 galaxies. The bonafide globular clusters are selected, on the basis of their magnitudes and colors, to be metal-rich. For this cluster sample we measure a mean radial velocity of v_{rad} = 1353 +/- 79 km/s, and a velocity dispersion of 338 +/- 56 km/s. Using a few different mass estimators, this implies M/L values in the range 50-130 within a radius of 28 kpc, consistent with M/L rising with radius. Our velocity dispersion estimate is intermediate between the value computed from the integrated stellar light at smaller radii, and that computed at large radii from recent X-ray observations.

We have used WFPC2 to construct B, V color-magnitude diagrams of four metal-rich globular clusters, NGC 104 (47 Tuc), NGC 5927, NGC 6388, and NGC 6441. All four clusters have well populated red horizontal branches (RHB), as expected for their metallicity. However, NGC 6388 and 6441 also exhibit a prominent blue HB (BHB) extension, including stars reaching as faint in V as the turnoff luminosity. This discovery demonstrates directly for the first time that a major population of hot HB stars can exist in old, metal-rich systems. This may have important implications for the interpretation of the integrated spectra of elliptical galaxies. The cause of the phenomenon remains uncertain. We examine the possibility that NGC 6388 and 6441 are older than the other clusters, but a simple difference in age may not be sufficient to produce the observed distributions along the HB. The high central densities in NGC 6388 and 6441 suggest that the existence of the blue HB (BHB) tails might be caused by stellar interactions in the dense cores of these clusters, which we calculate to have two of the highest collision rates among globular clusters in the Galaxy. Tidal collisions might act in various ways to enhance loss of envelope mass, and therefore populate the blue side of the HB. However, the relative frequency of tidal collisions does not seem large enough (compared to that of the clusters with pure RHBs) to account for such a drastic difference in HB morphology. While a combination of an age difference and dynamical interactions may help, prima facie the lack of a radial gradient in the BHB/RHB star ratio seems to argue against dynamical effects playing a role.

The stellar dynamics of Omega Centauri are inferred from the radial velocities of 469 stars measured with CORAVEL (Mayor et al. 1997). Rather than fit the data to a family of models, we generate estimates of all dynamical functions nonparametrically, by direct operation on the data. The cluster is assumed to be oblate and edge-on but mass is not assumed to follow light. The mean motions are consistent with axisymmetry but the rotation is not cylindrical. The peak rotational velocity is 7.9 km/s at 11 pc from the center. The apparent rotation of Omega Centauri is attributable in part to its proper motion. We reconstruct the stellar velocity ellipsoid as a function of position, assuming isotropy in the meridional plane. We find no significant evidence for a difference between the velocity dispersions parallel and perpendicular to the meridional plane. The mass distribution inferred from the kinematics is slightly more extended than, though not strongly inconsistent with, the luminosity distribution. We also derive the two-integral distribution function f(E,Lz) implied by the velocity data.

Galactic globular clusters, which are ancient building blocks of our Galaxy, represent a very interesting family of stellar systems in which some fundamental dynamical processes have taken place on time scales shorter than the age of the universe. In contrast with galaxies, these clusters represent unique laboratories for learning about two-body relaxation, mass segregation from equipartition of energy, stellar collisions, stellar mergers, and core collapse. In the present review, we summarize the tremendous developments, as much theoretical as observational, that have taken place during the last two decades, and which have led to a quantum jump in our understanding of these beautiful dynamical systems.

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 have identified 26 new globular cluster candidates in the inner 3 kpc of NGC 5128 (Cen A), the nearest known large galaxy being the probable product of a merger. The clusters are selected on the basis of their structural parameters (observed core diameters and ellipticities), as measured from archival WFPC HST images. IR photometry obtained with IRAC2B at the ESO/MPI 2.2-m telescope is combined with the optical HST photometry. Most of these clusters have normal colors typical of old globular clusters like those found in the Milky Way and M31. We estimate their metal abundances based on the R-K_0 color, confirming the existence of a metallicity gradient in the inner regions of NGC 5128. The presence of metal-rich globular clusters suggests that one of the colliding galaxies was a bulge-dominated galaxy (E or early S). A few clusters have colors and magnitudes similar to intermediate-age clusters containing carbon stars in the Magellanic Clouds. If the intermediate-age clusters were formed during a merger, then this episode must have occurred a few Gyr ago. Alternatively, we are looking at the cluster members of one of the colliding galaxies, which would then have been a late-type disk galaxy.

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.