Search

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., Comment: 9 pages, 4 figures, submitted to ApJ Letters

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

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., Comment: 19 pages, 7 figures, 4 tables. Accepted for publication in MNRAS. Light curves and simulated data can be downloaded at http://www.h0licow.org

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., Comment: Version revised to address referee's comments, submitted to MNRAS, 21 pages (incl. Appendix). Data associated to the paper available from the H0LICOW website www.h0licow.org

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., Comment: 16 pages, 3 figures, submitted to MNRAS and revised based on referee's comments

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

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

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

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

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., Comment: 10 pages, 6 figures, accepted in A&A

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, Comment: 11 pages, 8 figures, published in Astronomy & Astrophysics

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., Comment: Replaced with accepted version to ApJ

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.

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., Comment: 8 pages, 7 figures, published in Astronomy & Astrophysics

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)., Comment: 6 pages, 5 figures, accepted for publication in ApJL

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., Comment: White paper submitted to SNOWMASS2013

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., Comment: 10 pages, 10 figures, 6 tables

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., Comment: 6 pages, Proceedings of the Seyfert 2012 conference

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)., Comment: 23 pages, 7 figures; submitted to ApJS

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., Comment: 21 pages, 12 figures, minor revisions based on referee's comments, accepted for publication in ApJ

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., Comment: 11 pages, 9 figures, minor additions to the text only, techniques and results remain unchanged, A&A in press

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., Comment: 15 pages, 11 figures, stylistic improvements only, published in A&A

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)., Comment: 13 pages, 12 figures, Accepted for publication in MNRAS

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., Comment: Accepted for publication in Astronomy and Astrophysics. Main data set available via the German virtual observatory http://dc.g-vo.org/mlqso/q/web/form and soon via CDS. Additional material available on request

What are the frequency, shape, kinematics, and luminosity of Ly\alpha\ envelopes surrounding radio-quiet quasars at high redshift, and is the luminosity of these envelopes related to that of the quasar or not? As a first step towards answering these questions, we have searched for Ly\alpha\ envelopes around six radio-quiet quasars at z~4.5, using deep spectra taken with the FORS2 spectrograph attached to the UT1 of the Very Large Telescope (VLT). Using the multi-slit mode allows us to observe several point spread function stars simultaneously with the quasar, and to remove the point-like emission from the quasar, unveiling the faint underlying Ly\alpha\ envelope with unprecedented depth. An envelope is detected around four of the six quasars, which suggests that these envelopes are very frequent. Their diameter varies in the range 26

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., Comment: 8 pages, 6 figures. Metric values updated according to the final GREAT10 analysis software (Kitching et al. 2012, MNRAS 423, 3163-3208), no qualitative changes. Associated code available at http://lastro.epfl.ch/megalut

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

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., Comment: 4 pages, 1 figure, communique of workshop on the Hubble constant at KIPAC during February 6-8 2012

(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., Comment: Accepted for publication in Astronomy and Astrophysics abridged abstract

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., Comment: 9 pages, 8 figures, accepted for publication in A&A

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., Comment: 116 pages, with executive summary and table of contents

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., Comment: Accepted for publication in A&A, 6 pages, 4 figures

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., Comment: Accepted for publication in A&A. 8 pages, 5 figures