Characterising the environment and the line of sight of gravitationally lensed quasars using weak lensing technique

Estimates of the Hubble constant from different cosmological probes produce discrepant results in the framework of the standard cosmological model. In order to assess whether this tension is coming from unaccounted systematic errors, or is hinting towards an extended cosmological model, independent and precise techniques to determine the Hubble constant are required.

The time-delay cosmography is an emerging probe, which relies on strongly lensed variable sources. It enables to measure the Hubble constant with just a few percent precision. Current efforts are capitalised by the H0LiCOW collaboration, which uses lensed quasars as variable sources. Three analysis steps contribute to the error budget of the final estimation: measurements of the time delays of the lensed images, modeling of the gravitational potential of the lens, characterisation of the environment and the line of sight of the lensed system.

H0LiCOW benefits from two different methods to estimate the contribution of the environment and the line of sight: the number counts technique and the weak lensing technique. The number counts technique is a statistical approach, which uses weighted galaxy number density in the field of the lensed systems coupled with dark matter simulations. The weak lensing technique uses the weak gravitational lensing phenomenon, encompassing the fact that the shapes of the background galaxies are unbiased tracers of matter. These two methods are complimentary and almost independent, which allows for an efficient cross-check of the systematic errors.

During this thesis I developed the pipeline for the weak lensing technique, which provides quantitative characterisation of the environment and the line of sight of gravitationally lensed quasars. I worked with ground- and space-based data, and utilised state-of-art techniques, including inpainting and wavelet transform, to carry out the analysis. I applied the developed pipeline on two lensed systems, HE 0435-1223 and B1608+656. I determined that while the contribution of the structures along the line of sight is negligible for the first, it is quite significant for the second system, requiring a thorough analysis to extract the cosmological information. These results are in perfect agreement with the estimates from number counts technique. By verifying possible systematics, this gives an important support for time-delay cosmography in face of the tension in the estimates of the Hubble constant.

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses X. Modeling based on high-precision astrometry of a sample of 25 lensed quasars: consequences for ellipticity, shear, and astrometric anomalies

Frédéric Courbin, Georges Meylan

Gravitationally lensed quasars can be used as powerful cosmological and astrophysical probes. We can (i) infer the Hubble constant H-0 based on the so-called time-delay technique, (ii) unveil substructures along the line-of-sight toward distant galaxies, and (iii) compare the shape and the slope of baryons and dark matter distributions in the inner regions of galaxies. To reach these goals, we need high-accuracy astrometry of the quasar images relative to the lensing galaxy and morphology measurements of the lens. In this work, we first present new astrometry for 11 lenses with measured time delays, namely, JVAS B0218+357, SBS 0909+532, RX J0911.4+0551, FBQS J0951+2635, HE 1104-1805, PG 1115+080, JVAS B1422+231, SBS 1520+530, CLASS B1600+434, CLASS B1608+656, and HE 2149-2745. These measurements proceed from the use of the Magain-Courbin-Sohy (MCS) deconvolution algorithm applied in an iterative way (ISMCS) to near-IR HST images. We obtain a typical 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 lensing systems, mostly from the COSMOGRAIL set of targets, to present new mass models of these lenses. The modeling of these 25 gravitational lenses led to the following results: 1) in four double-image quasars (HE0047-1746, J1226-006, SBS 1520+530, and HE 2149-2745), 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 problem) is overcome by explicitly including the nearest visible galaxy/satellite in the lens model. However, one anomalous system (RXS J1131-1231) does not show any luminous perturber in its vicinity, and three others (WFI 2026-4536, WFI 2033-4723, and B2045+265) have problematic modeling. These four systems are the best candidates for a pertubation by a dark matter substructure along the line-of-sight; 3) we revisit the correlation between the position angle (PA) and ellipticity of the light and of the mass distribution in lensing galaxies. As in previous studies, we find a significant correlation between the PA of the light and of the mass distributions. However, in contrast with these same 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 similar to 5 kpc of lensing galaxies. This offers a new test for high-resolution hydrodynamical simulations.

2012

Gravitational lenses

Christel Vuissoz

The Hubble constant H0 is one of the most important parameters in cosmology, as it encodes the age of the Universe and is necessary for any distance determination at a cosmological scale. It is, however, only poorly constrained by traditional methods. The current favored value, H0 = 72±8 km s-1 Mpc-1, is provided by the HST Hubble constant Key Project (Freedman et al. 2001), which combines several Cepheid-calibrated distance indicators. This roughly 10% error nevertheless denotes only the statistical uncertainty in the determination of H0, while the possible systematical errors in the first step of the distance ladder (the distance to the Large Magellanic Cloud) may be of the same order of magnitude. Time delays between gravitationally lensed images of distant quasars can yield a more precise measurement of the Hubble constant, on a truly cosmic scale, and independently of any local distance calibrator. At the beginning of this thesis, time delays had been measured in only ten lensed systems, nine of which gave H0 estimates. However before 2004, no concerted and long term action has succeeded to apply the time delay method at a level of precision really competitive with other techniques. The major difficulties arise from the modeling of the lens mass distribution, and from the uncertainty on the time delay measurement itself, which was typically of about 10% in past monitoring programs. COSMOGRAIL (COSmological MOnitoring of GRAvItational Lenses) is an international collaboration initiated in April 2004 at the Laboratory of Astrophysics of EPFL, and which aims at measuring precise time delays for most known lensed quasars, in order to determine the Hubble constant down to an uncertainty of a few percent. This thesis took place at the beginning of COSMOGRAIL and consisted in setting up this large photometric monitoring. It addressed both issues of carrying out accurate photometry of faint blended sources and of obtaining well sampled light curves, in order to measure precise time delays. As part of the COSMOGRAIL project, I have been managing the monitoring of over twenty gravitationally lensed quasars with the 1-2m telescopes involved both in the Northern and Southern hemispheres, and organizing the data. The first crucial work of this thesis was then to develop an automated reduction pipeline able to produce an homogeneous data set from images acquired with very different telescopes. This pipeline was also needed to perform aperture photometry of all lensed quasars, in order to study their variability and define the monitoring priorities. The powerful MCS deconvolution algorithm (Magain, Courbin, & Sohy 1998) was greatly used in this work and allowed to highly improve the image resolution, with the aim of obtaining accurate photometric measurements of the individual quasar lensed images. I have finally tested and improved three different numerical techniques to determine time delays between the quasar components from their light curves. In this thesis, time delays have been determined in four systems. The first one was measured in the doubly imaged quasar SDSS J1650+4251, after two years of monitoring with the 1.5m telescope of Maidanak Observatory, in Uzbekistan. The quadruply lensed system RXS J1131–1231 was then studied and three time delays determined from 3-year observations with the Swiss Euler 1.2m telescope located at La Silla, in Chile. The photometric monitoring of the quadruple WFI J2033–4723 was also carried out with the Euler telescope, and data were then merged with those obtained by a second monitoring group, with the SMARTS 1.3m telescope at the Cerro Tololo Interamerican Observatory (CTIO), also located in Chile. Two time delays were measured in this system, after three years of observations, the close pair A1 – A2 remaining unresolved. Three time delays were determined in the quadruply imaged quasar HE0435–1223, after four years of optical monitoring with Euler, Mercator and Maidanak telescopes, to which photometric measurements by SMARTS 1.3m telescope were added. Euler and SMARTS merged data for the doubly imaged quasar QJ0158–4325 were also analysed, the size of the source accretion disk was measured, but we failed to determine a time delay due to the high amplitude of the microlensing variability in this system. The accuracies on time delay measurements reached in this thesis are of the order of 3-4% and show a clear improvement from the typical 10% uncertainties of past monitoring programs. These results were finally converted into estimates of the Hubble constant following different models of the lensing mass potential. The H0 mean value obtained when considering the individual determinations from twelve gravitationally lensed quasars with known time delays is H0 = 60 ± 7 km s-1 Mpc-1. This result is consistent with the current favored value, and above all promising, as including additional systems to this ensemble will surely provide tighter bounds on H0. In conclusion, the increasing number of time delay measurements and improvements in lens modeling should reduce the errors on the Hubble constant estimate provided by gravitational lensing. Conversely, the determination of more time delays should put further constraints on lens galaxy density profiles when using a prior on H0 from other studies.

EPFL2008

Astrophysical applications of gravitationally lensed quasars

Alexander Eigenbrod

Gravitational lensing describes how light is deflected as it passes in the vicinity of a mass distribution. The amplitude of the deflection is proportional to the mass of the deflector, called "gravitational lens", and is generally weak, even for large masses. The faintness of this phenomenon explains why gravitational lensing remained essentially unobserved until the late 1970s (only gravitational lensing by the Sun has been observed during the solar eclipse of 1919). Before that time, gravitational lensing was considered merely as a theoretical curiosity. However, the situation dramatically changed with the discovery of the first extragalactic gravitational lens in 1979. Since then, together with the technological progress of astronomical instruments, gravitational lensing has turned from a curiosity into a powerful tool to address important astrophysical and cosmological questions. The present thesis focuses on applications related to gravitationally lensed quasars. Quasars are active galactic nuclei, where matter is heated up as it spirals down onto the central supermassive black hole. When a galaxy is located on the line of sight to a distant quasar, it acts as a gravitational lens and produces multiple images of this background source. The light of the quasar follows different paths for each of its images. Thus, variations of the intrinsic quasar luminosity are observed at different times in each image. The time delays between the images can be used to determine the Hubble constant H0, because they are inversely proportional to H0. This constant describes the current expansion rate of the Universe, and is one of the fundamental parameters of cosmological models. Many efforts have been spent over the years to determine H0, but its value is still poorly constrained. Gravitational lensing has the potential to noticeably decrease the uncertainty of H0. In practice, this requires regular and long-term monitoring of lensed quasars. We have run a series of numerical simulations to both optimize the available telescope time, and measure the time delays with an accuracy of a few percent. The results of these simulations are presented in the form of compact plots to be used to optimize the observational strategy of present and future monitoring programs. Once the time delays are measured, one can infer estimates of H0, provided several other observational constraints are available. A key element to accurately convert time delays into H0 is the redshift of the lensing galaxy. These redshift measurements are difficult because lensing galaxies are generally hidden in the glare of the much brighter quasar images. As a consequence, lens redshifts are often poorly constrained or even completely unknown. We have acquired spectroscopic data of sixteen lensing galaxies with the Very Large Telescope located in Chile. In combination with a powerful deconvolution algorithm, we determine the redshift of these sixteen lensing galaxies, which represents about 25% of all currently known quasar lensing galaxies. These results are useful for both H0 determinations and statistical studies of gravitational lenses, which can be used to provide new constraints on cosmological parameters. While the first part of this thesis focuses on the acquisition of observational constraints for the lens models, the main part consists in using the phenomenon of microlensing to determine the energy profile (or spatial structure) of quasar accretion disks. Microlensing is produced by the stars located in the lensing galaxy. These stars act as secondary lenses, and are called microlenses. Since the stars are moving in the galaxy, they induce flux and color variations in the images of the lensed quasar. These effects can be used as a natural telescope to probe the still mysterious inner structures of quasars with a spatial resolution about ten thousand times better than the capacities of current astronomical instruments, including the Very Large Telescope Interferometer. We present a three-year long spectrophotometric monitoring of the lensed quasar QSO 2237+0305, also known as the Einstein Cross, conducted at the Very Large Telescope. This monitoring reveals significant microlensing-induced variations in the spectra of the quasar images. In a subsequent analysis, we find that the source responsible for the optical and ultraviolet continuum has an energy profile well reproduced by a power-law R α λζ with ζ = 1.2 ± 0.3, where R is the size of the source emitting at wavelength λ. This agrees with the predictions of the standard thin accretion disk model and is, so far, the most accurate determination of a quasar energy profile. As a complement to our microlensing study, we have obtained high spectral and spatial resolution observations of the lensing galaxy of QSO 2237+0305. Our spectroscopic data are acquired with the SINFONI, FLAMES, and FORS2 spectrographs of the Very Large Telescope. We describe the reduction of these data, and provide the currently best and most complete determination of the kinematics of a gravitational lens. The comparison of our data with previously published dynamical models suggests that those may have overestimated the mass of the galaxy bulge. Thus, new and more sophisticated models are required. These models, combined with gravitational lensing, will provide two independent constraints on the mass distribution. This will allow to better determine the quantity and distribution of dark matter in this lensing galaxy, and especially in its extended halo.

EPFL2008

Frédéric Courbin, Georges Meylan

2012

Astrophysical applications of gravitationally lensed quasars

Alexander Eigenbrod

EPFL2008

Gravitational lenses

Christel Vuissoz

EPFL2008