Search

Modeling the mass distribution of galaxy-scale strong gravitational lenses is a task of increasing difficulty. The high-resolution and depth of imaging data now available render simple analytical forms ineffective at capturing lens structures spanning a large range in spatial scale, mass scale, and morphology. In this work, we address the problem with a novel multiscale method based on wavelets. We tested our method on simulated Hubble Space Telescope (HST) imaging data of strong lenses containing the following different types of mass substructures making them deviate from smooth models: (1) a localized small dark matter subhalo, (2) a Gaussian random field (GRF) that mimics a nonlocalized population of subhalos along the line of sight, and (3) galaxy-scale multipoles that break elliptical symmetry. We show that wavelets are able to recover all of these structures accurately. This is made technically possible by using gradient-informed optimization based on automatic differentiation over thousands of parameters, which also allow us to sample the posterior distributions of all model parameters simultaneously. By construction, our method merges the two main modeling paradigms - analytical and pixelated - with machine-learning optimization techniques into a single modular framework. It is also well-suited for the fast modeling of large samples of lenses. All methods presented here are publicly available in our new Herculens package., 24 pages, 12 figures, accepted for publication in A&A

For many years, observations of the Universe suggest a series problems with our theoretical models, particularly its dark energy and dark matter components. Fortunately, the Universe also provides us with a tool to solve these problems, called strong gravitational lensing. This natural phenomenon, observed as multiple images of a distant background source, distorted by the mass of a foreground galaxy, offers a unique opportunity to detect the otherwise invisible total mass of that galaxy. Yet, this detection is only possible if both background and foreground objects are accurately modeled, a task of increasing difficulty because our telescopes reveal more and more of the intrinsic complexity of galaxies. The primary goal of this thesis is to demonstrate how novel and well-motivated techniques can go beyond the current simplifying assumptions to improve the modeling of luminous and dark components of galaxies. Gravitational lens modeling is an under-constrained problem that has several possible solutions. The new techniques I introduce here are based on two key concepts to overcome this difficulty: sparsity and the wavelet transform. Sparsity selects the solution that contains the fewest parameters, namely the least complex one that best fits the observation. The wavelet transform separates the various spatial scales of the solution and enables the reconstruction of the small-scale compact features up to the larger, smoother variations, which are all found in real galaxies. Since the complexity and size of data sets are dramatically increasing, modeling techniques must also be fast and scalable. I address these requirements using differentiable programming to enable unprecedented gains in computation time. The proposed modeling framework allows us to effortlessly combine simple and more complex techniques together, if required by the observations. In this work, I demonstrate that sparsity and wavelets can address the limitations of current methods for modeling the full complexity of gravitational lenses. Compared to the many methods based on smoothness assumptions, I show how multi-scale modeling techniques significantly improve the reconstruction of lensed galaxies at high resolution. Moreover, I demonstrate that those same techniques are well-suited to characterize the invisible mass distributions of galaxies, notably when it deviates from the widely used smooth elliptical profiles. These results offer exciting possibilities to measure better the properties of galaxies via gravitational lensing, including their dark matter content, ultimately improving our understanding of galaxy evolution. Additionally, I take part to the long-standing debate regarding the role of dark energy in the expansion of the Universe using the method of time-delay cosmography. Based on the gravitational lensing of distant quasars, this method plays a central role in this context because it can measure the expansion rate of the Universe (the Hubble constant) independently of all other methods. The results I present are various: searching for systematic errors in past measurements, testing modeling techniques on a blind challenge, and modeling recent Hubble Space Telescope observations of lensed quasars to measure their absolute distance. I also describe the new approach introduced within the TDCOSMO collaboration, based on relaxing most assumptions on the mass distribution of galaxies and replacing those with observations of stellar kinematics.

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.

There is now very strong evidence that our Universe is undergoing an accelerated expansion period as if it were under the influence of a gravitationally repulsive “dark energy” component. Furthermore, most of the mass of the Universe seems to be in the form of non-luminous matter, the so-called “dark matter”. Together, these “dark” components, whose nature remains unknown today, represent around 96 % of the matter-energy budget of the Universe. Unraveling the true nature of the dark energy and dark matter has thus, obviously, become one of the primary goals of present-day cosmology. Weak gravitational lensing, or weak lensing for short, is the effect whereby light emitted by distant galaxies is slightly deflected by the tidal gravitational fields of intervening foreground structures. Because it only relies on the physics of gravity, weak lensing has the unique ability to probe the distribution of mass in a direct and unbiased way. This technique is at present routinely used to study the dark matter, typical applications being the mass reconstruction of galaxy clusters and the study of the properties of dark halos surrounding galaxies. Another and more recent application of weak lensing, on which we focus in this thesis, is the analysis of the cosmological lensing signal induced by large-scale structures, the so-called “cosmic shear”. This signal can be used to measure the growth of structures and the expansion history of the Universe, which makes it particularly relevant to the study of dark energy. Of all weak lensing effects, the cosmic shear is the most subtle and its detection requires the accurate analysis of the shapes of millions of distant, faint galaxies in the near infrared. So far, the main factor limiting cosmic shear measurement accuracy has been the relatively small sky areas covered. Next-generation of wide-field, multicolor surveys will, however, overcome this hurdle by covering a much larger portion of the sky with improved image quality. The resulting statistical errors will then become subdominant compared to systematic errors, the latter becoming instead the main source of uncertainty. In fact, uncovering key properties of dark energy will only be achievable if these systematics are well understood and reduced to the required level. The major sources of uncertainty resides in the shape measurement algorithm used, the convolution of the original image by the instrumental and possibly atmospheric point spread function (PSF), the pixelation effect caused by the integration of light falling on the detector pixels and the degradation caused by various sources of noise. Measuring the Cosmic shear thus entails solving the difficult inverse problem of recovering the shear signal from blurred, pixelated and noisy galaxy images while keeping errors within the limits demanded by future weak lensing surveys. Reaching this goal is not without challenges. In fact, the best available shear measurement methods would need a tenfold improvement in accuracy to match the requirements of a space mission like Euclid from ESA, scheduled at the end of this decade. Significant progress has nevertheless been made in the last few years, with substantial contributions from initiatives such as GREAT (GRavitational lEnsing Accuracy Testing) challenges. The main objective of these open competitions is to foster the development of new and more accurate shear measurement methods. We start this work with a quick overview of modern cosmology: its fundamental tenets, achievements and the challenges it faces today. We then review the theory of weak gravitational lensing and explains how it can make use of cosmic shear observations to place constraints on cosmology. The last part of this thesis focuses on the practical challenges associated with the accurate measurement of the cosmic shear. After a review of the subject we present the main contributions we have brought in this area: the development of the gfit shear measurement method, new algorithms for point spread function (PSF) interpolation and image denoising. The gfit method emerged as one of the top performers in the GREAT10 Galaxy Challenge. It essentially consists in fitting two-dimensional elliptical Sérsic light profiles to observed galaxy image in order to produce estimates for the shear power spectrum. PSF correction is automatic and an efficient shape-preserving denoising algorithm can be optionally applied prior to fitting the data. PSF interpolation is also an important issue in shear measurement because the PSF is only known at star positions while PSF correction has to be performed at any position on the sky. We have developed innovative PSF interpolation algorithms on the occasion of the GREAT10 Star Challenge, a competition dedicated to the PSF interpolation problem. Our participation was very successful since one of our interpolation method won the Star Challenge while the remaining four achieved the next highest scores of the competition. Finally we have participated in the development of a wavelet-based, shape-preserving denoising method particularly well suited to weak lensing analysis.

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.