F. Courbin

Results of the GREAT08 Challenge: an image analysis competition for cosmological lensing

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

Image analysis for cosmology: results from the GREAT10 Galaxy Challenge

We present the results from the first public blind point-spread function (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 modeled ellipticity e and size R^2 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 nine teams. The best methods reconstructed the PSF with an accuracy of σ(e) ≈ 2.5 × 10^(–4) and σ(R^2)/R^2 ≈ 7.4 × 10^(–4). 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).

H0LiCOW - V. New COSMOGRAIL time delays of HE 0435-1223: H0 to 3.8 per cent precision from strong lensing in a flat ΛCDM model

We present a new measurement of the Hubble Constant H-0 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 HE 0435-1223 from 13-yr light curves obtained as part of the COSMOGRAIL project. Companion papers detail the modelling 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 HE 0435-1223 with previous measurements from systems B1608+656 and RXJ1131-1231 to create a Time Delay Strong Lensing probe (IDSL). In flat A cold dark matter (ACDM) with free matter and energy density, we find H-0 = 71.9(-3.0)(+2.4) km s(-1) Mpc(-1) and Omega(Lambda) = 0.62(-0.35)(+0.24) This measurement is completely independent of, and in agreement with, the local distance ladder measurements of H-0. We explore more general cosmological models combining TDSL with other probes, illustrating its power to break degeneracies inherent to other methods. The joint constraints from IDSL and Planck are H-0 = 69.2(-2.2)(+1.4) km s(-1) Mpc(-1), Omega(Lambda) = 0.70(-0.01)(+0.01) and Omega(k) = 0.003(-0.006)(+0.004) in open ACDM and H-0 = 79.0(-4.2)(+4.4) km s(-1) Mpc(-1), Omega(de) = 0.77(-0.03)(+0.02) and w = -1.38(-0.16)(+0.14) in flat wCDM. In combination 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) in N-eff Lambda CDM and when relaxing the total mass of neutrinos we find Sigma rn(nu) <= 0.182 eV in m(nu) Lambda CDM. Finally, in an open wCDM in combination with Planck and cosmic microwave background lensing, we find H-0 = 77.9(-4.2)(+5.0) km s(-1) Mpc(-1), Omega(de) = 0.77(-0.03)(+0.03), Omega(k) = -0.003(-0.004)(+0.004) and w = -1.37(-0.23)(+0.18).

High accuracy transit photometry of the planet OGLE-TR-113b with a new deconvolution-based method

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

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses - XIII. Time delays and 9-yr optical monitoring of the lensed quasar RX J1131−1231

We present the results from nine years of optically monitoring the gravitationally lensed z(QSO) = 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 measured 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 measured 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.

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses - IX. Time delays, lens dynamics and baryonic fraction in HE 0435-1223

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 ΔtBC = 7.8 ± 0.8 days, ΔtBD = -6.5 ± 0.7 days and ΔtCD = -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, fb, in the Einstein radius. We measured fb = 0.65-0.10+0.13 if the lensing galaxy has a Salpeter IMF and fb = 0.45-0.07+0.04 if it has a Kroupa IMF. The spectrum also allowed us to estimate the velocity dispersion of the lensing galaxy, σap = 222 ± 34 km s-1. We used fb and σap to constrain an analytical model of the lensing galaxy composed of an Hernquist plus generalized NFW profile. We solved 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 fb and σap, we did not constrain H0 yet with high accuracy, i.e., we found a broad range of models with χ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 ourconstraints with lens reconstruction techniques that make use of the full Einstein ring that is visible in this object. Based on observations made with the 1.2 m Euler Swiss Telescope, the 1.5 m telescope of Maidanak Observatory in Uzbekistan, and with the 1.2 m Mercator Telescope, operated on the island of La Palma by the Flemish Community, at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. The NASA/ESA Hubble Space Telescope data was obtained from the data archive at the Space Telescope Science Institute, which is operated by AURA, the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS-5-26555.Light curves are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/536/A53

Microlensing variability in the gravitationally lensed quasar QSO 2237+0305 = the Einstein Cross. II. Energy profile of the accretion disk

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 proportional to lambda(zeta) with zeta = 1.2 +/- 0.3, where R is the source size responsible for the emission at wavelength.. This is the first accurate, model-independent determination of the energy profile of a quasar accretion disk on such small scales.

COSMOLOGY FROM GRAVITATIONAL LENS TIME DELAYS AND PLANCK DATA

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 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 measurement 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.02)(+0.01) (68% credible interval (CI)). In a flat wCDM model, the combination of RXJ1131-1231 and Planck yields w = -1.52 (+0.19)(-0.20) (68% CI).

Further Evidence That Quasar X-Ray Emitting Regions Are Compact: X-Ray And Optical Microlensing In The Lensed Quasar Q J0158-4325

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 the 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.5)(+0.4), 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.

THE THIRD GRAVITATIONAL LENSING ACCURACY TESTING (GREAT3) CHALLENGE HANDBOOK

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