Abstract
We present a blind time-delay strong lensing (TDSL) cosmographic analysis of the doubly imaged quasar SDSS 1206+4332 . We combine the relative time delay between the quasar images, Hubble ...Space Telescope imaging, the Keck stellar velocity dispersion of the lensing galaxy, and wide-field photometric and spectroscopic data of the field to constrain two angular diameter distance relations. The combined analysis is performed by forward modelling the individual data sets through a Bayesian hierarchical framework, and it is kept blind until the very end to prevent experimenter bias. After unblinding, the inferred distances imply a Hubble constant H0 = 68.8$^{+5.4}_{-5.1}$ km s−1 Mpc−1, assuming a flat Λ cold dark matter cosmology with uniform prior on Ωm in 0.05, 0.5. The precision of our cosmographic measurement with the doubly imaged quasar SDSS 1206+4332 is comparable with those of quadruply imaged quasars and opens the path to perform on selected doubles the same analysis as anticipated for quads. Our analysis is based on a completely independent lensing code than our previous three H0LiCOW systems and the new measurement is fully consistent with those. We provide the analysis scripts paired with the publicly available software to facilitate independent analysis (footnote with link to www.h0licow.org). The consistency between blind measurements with independent codes provides an important sanity check on lens modelling systematics. By combining the likelihoods of the four systems under the same prior, we obtain H0 = 72.5$^{+2.1}_{-2.3}$ km s−1 Mpc−1. This measurement is independent of the distance ladder and other cosmological probes.
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.
We present a new measurement of the Hubble Constant H sub( 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 (TDSL). In flat Lambda cold dark matter ( Lambda CDM) with free matter and energy density, we find H sub( 0) =71.9... km s super( -1) Mpc super( -1) and Omega sub( Lambda )=0.62... This measurement is completely independent of, and in agreement with, the local distance ladder measurements of H sub( 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 TDSL and Planck are H sub( 0) = 69.2... km s super( -1) Mpc super( -1) , Omega sub( Lambda )=0.70... and Omega sub( k)=0.003+0.004-0.006 in open ...CDM and H sub( 0) =79.0... km s super( -1) Mpc super( -1), Omega sub( de)=0.77... and w=-1.38... 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 sub( eff) = 3.34... in N sub( eff) Lambda CDM and when relaxing the total mass of neutrinos we find ...m sub( ...) less than or equal to ...0.182 eV in m sub( nu )...CDM. Finally, in an open wCDM in combination with Planck and cosmic microwave background lensing, we find H sub( 0) =77.9... km s super( -1) Mpc super( -1), Omega sub( de) = 0.77..., Omega sub( k) = -0.003... and w=-1.37... (ProQuest: ... denotes formulae/symbols omitted.)
COSMOGRAIL is a long-term photometric monitoring of gravitationally lensed quasars aimed at implementing Refsdal's time-delay method to measure cosmological parameters, in particular H sub(0). Given ...the long and well sampled light curves of strongly lensed quasars, 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 chi super(2) 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.
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). In all our methods we model the PSF using a single Moffat profile and we interpolate the fitted parameters at a set of required positions. This allowed us to win the Star-challenge of GREAT10, with the B-splines method. However, 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 σ2sys better than the 1 × 10-7 upper bound expected by future space-based surveys, with the local interpolators performing better than the global ones.
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, the COSmological MOnitoring of GRAvItational Lenses, (2) archival Hubble Space Telescope imaging of the lens system, (3) a new velocity-dispersion measurement of the lens galaxy of 323 + or- 20 km ssup -1 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. The joint analysis of the data sets allows a time-delay distance measurement to 6% precision that takes into account all known systematic uncertainties. 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.
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 with bolometric luminosities between 1044.7 − 47.4 erg/s and black hole masses 107.6 − 9.8 M⊙. This sample is composed of 13 pairs of previously unpublished spectra and four pairs of spectra from literature. Our analysis is based on a simple 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 emission lines in our spectra (C III and Mg II) 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.
Abstract
Galaxies located in the environment or along 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 HE 0435−1223 using the W. M. Keck, Gemini and ESO-Very Large telescopes. Our new catalogue triples the number of known galaxy redshifts in the direct 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 catalogue with literature data to gather redshifts up to 15 arcmin from the lens, and search for galaxy groups or clusters projected towards HE 0435−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.
ABSTRACT
We investigate the environment and line of sight (LoS) of the H0LiCOW (H0 Lenses in COSMOGRAIL’s Wellspring) lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to ...perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser–Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and ${\rm {\small {glimpse}}}$, a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data are less sensitive to small-scale structures than the space-based observations. Marginalizing over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, κext, at the position of B1608+656 is $\kappa _{\mathrm{ext}}= 0.11^{+0.06}_{-0.04}$, where the error bars correspond, respectively, to the 16th and 84th quartiles. This result is compatible with previous estimates using the number counts technique, suggesting that B1608+656 resides in an overdense LoS, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favour truncated halo models.