Time-delay Cosmography with Strongly Lensed Quasars

The LCDM model has emerged as the concordance model of cosmology for its ability to explain a variety of observations, ranging from the anisotropies of the Cosmic Microwave Background (CMB) to the accelerated expansion of the Universe. However, with the ever-increasing precision of the measurements, tensions have recently emerged between the latest CMB observations and other cosmological probes. The most prominent one concerns the expansion rate of the Universe, that is the Hubble constant, H0, with a statistically significant discrepancy between local and early measurements. If not arising from unaccounted systematic errors, this discrepancy is an exciting opportunity to explore missing physics beyond the standard LCDM model.Time-delay cosmography has emerged as a competitive and single-step method to measure locally the Hubble constant, providing an alternative to the cosmic distance ladder. The two techniques do not share any known source of errors, allowing us to obtain truly independent estimates. This technique is based on the strong lensing phenomenon, occurring when the light coming from a background source is split into multiple images by a massive foreground galaxy. In asymmetric configurations, the travel time is then slightly different for each of the lensed images, leading to measurable time delays in the signals received from the background object. The Hubble constant can then be directly inferred from these time delays, provided that an accurate reconstruction of the mass distribution of the lens galaxy is available.The first part of my PhD work focuses on the precise determination of the time delays in strongly lensed quasars, an essential ingredient of the method. The COSMOGRAIL collaboration has been monitoring lensed quasars for two decades to record enough quasar variations that can be unambiguously matched in all light curves. I compiled these data obtained from the Euler 1.2m Swiss telescope in La Silla, Chile, and measured the time delays in 18 systems, more than doubling the current sample of lensed quasars with known delays. This work required the development and automation of curve-shifting algorithms presented in a dedicated chapter. As the sample of known lensed quasars is growing rapidly, I also present a novel monitoring strategy to obtain the time delays efficiently, in only one or two monitoring seasons. Thanks to this long-term monitoring effort, the TDCOSMO collaboration published the most precise determination of the Hubble constant obtained with time-delay cosmography from a sample of seven lensed quasars. We found H0 = 73.7+/- 1.5 km/s/Mpc, at 2% precision. This result confirmed the existing tension, now reaching a high significance when combined with the distance ladder results. "As extraordinary claims require extraordinary evidence", the second part of my work focused on the verification of the assumptions made to obtain this result. In the last chapter of this thesis, I present my work on the search for unaccounted systematic errors in every step of the analysis. I also show the results of our participation to a blind data challenge designed to test lens modelling codes and the assumptions used to reconstruct the mass distribution of the galaxies. I conclude this thesis by presenting the new approach proposed by the TDCOSMO collaboration to relax all assumptions about the density profile of massive elliptical galaxies and adopt a parametrisation entirely constrained by stellar kinematic.

Precise measurements of time delays in gravitationally lensed quasars for competitive and independent determination of the Hubble constant

Vivien François Bonvin

During these last decades, by virtue of observations, the Standard Cosmological Model has emerged, providing a description of the Universe's evolution using a minimal set of independent constraints - the cosmological parameters. Among them is the expansion rate of the Universe, the so-called Hubble constant or H0, first measured by Lemaître in 1927. The century that followed this cornerstone measurement saw numerous attempts to refine the initial value, and for good reason: a precise and independent measurement of H0 will bring strong constraints on the cosmological models. It could notably help the astronomers to better understand the nature of dark energy, thus making it one of the most sought-after prizes in modern cosmology. My work at the Laboratory of Astrophysics of EPFL is embedded in this context. I am part of the COSMOGRAIL and H0LiCOW collaborations, aiming to measure the Hubble constant with the highest level of precision using time-delay cosmography, a method based on the theory of strong gravitational lensing. This effect occurs when an observer looks at a light source located behind a massive foreground galaxy. The mass of the galaxy acts similarly to an optical lens and focuses the light rays emitted by the source. As a consequence, multiple lensed images of the source appear around the lens galaxy. If the luminosity of the source changes over time, the variations will be seen in all the lensed images but with a temporal delay due to the different travel paths of the light rays. By carefully monitoring the luminosity variations of each lensed image, one can precisely measure the temporal delays between them. Combined to high-resolution observations of the foreground galaxy and its surroundings, it is possible to directly measure the Hubble constant upon the sole assumption that the General Relativity is correct. Since more than 13 years, COSMOGRAIL monitors dozens of lensed quasars to produce high-quality light curves and time-delay measurements. During these last four years, I took care of the monitoring schedule, continuous data reduction and time-delay measurements through the development of curve-shifting techniques. I produced light curves and measured time delays on a variety of lenses. After more than a decade of endeavours, COSMOGRAIL and H0LiCOW finally revealed their measurement of the expansion rate of the Universe from a blind analysis of three lensed sources. I had the privilege to be the lead author of the publication presenting our measurement of the Hubble constant, H0=71.9 -3.0+2.4 km/s/Mpc 3.8% precision in the Standard Cosmological Model. Such a precision allows a direct comparison with the results of the distance ladder technique in the local Universe and the Planck satellite Cosmic Microwave Background observations in the distant Universe, both of which being currently in a significant tension of unknown source.

EPFL2017

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

Gravitational lenses

Christel Vuissoz

EPFL2008

Astrophysical applications of gravitationally lensed quasars

Alexander Eigenbrod

EPFL2008

Precise measurements of time delays in gravitationally lensed quasars for competitive and independent determination of the Hubble constant

Vivien François Bonvin

EPFL2017