skip to main content

Title: STRIDES: automated uniform models for 30 quadruply imaged quasars

Gravitational time delays provide a powerful one-step measurement of H0, independent of all other probes. One key ingredient in time-delay cosmography are high-accuracy lens models. Those are currently expensive to obtain, both, in terms of computing and investigator time (105–106 CPU hours and ∼0.5–1 yr, respectively). Major improvements in modelling speed are therefore necessary to exploit the large number of lenses that are forecast to be discovered over the current decade. In order to bypass this roadblock, we develop an automated modelling pipeline and apply it to a sample of 31 lens systems, observed by the Hubble Space Telescope in multiple bands. Our automated pipeline can derive models for 30/31 lenses with few hours of human time and <100 CPU hours of computing time for a typical system. For each lens, we provide measurements of key parameters and predictions of magnification as well as time delays for the multiple images. We characterize the cosmography-readiness of our models using the stability of differences in the Fermat potential (proportional to time delay) with respect to modelling choices. We find that for 10/30 lenses, our models are cosmography or nearly cosmography grade (<3 per cent and 3–5 per cent variations). For 6/30 lenses, the models are close to cosmography grade (5–10 per cent). These results utilize informative priors and will need to be confirmed by further analysis. However, they are also likely to improve by extending the pipeline modelling sequence and options. In conclusion, we show that uniform cosmography grade modelling of large strong lens samples is within reach.

more » « less
Award ID(s):
1906976 1907396
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Corporate Creator(s):
Publisher / Repository:
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
1260 to 1300
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. When strong gravitational lenses are to be used as an astrophysical or cosmological probe, models of their mass distributions are often needed. We present a new, time-efficient automation code for the uniform modeling of strongly lensed quasars withGLEE, a lens-modeling software for multiband data. By using the observed positions of the lensed quasars and the spatially extended surface brightness distribution of the host galaxy of the lensed quasar, we obtain a model of the mass distribution of the lens galaxy. We applied this uniform modeling pipeline to a sample of nine strongly lensed quasars for which images were obtained with the Wide Field Camera 3 of theHubbleSpace Telescope. The models show well-reconstructed light components and a good alignment between mass and light centroids in most cases. We find that the automated modeling code significantly reduces the input time during the modeling process for the user. The time for preparing the required input files is reduced by a factor of 3 from ~3 h to about one hour. The active input time during the modeling process for the user is reduced by a factor of 10 from ~ 10 h to about one hour per lens system. This automated uniform modeling pipeline can efficiently produce uniform models of extensive lens-system samples that can be used for further cosmological analysis. A blind test that compared our results with those of an independent automated modeling pipeline based on the modeling softwareLenstronomyrevealed important lessons. Quantities such as Einstein radius, astrometry, mass flattening, and position angle are generally robustly determined. Other quantities, such as the radial slope of the mass density profile and predicted time delays, depend crucially on the quality of the data and on the accuracy with which the point spread function is reconstructed. Better data and/or a more detailed analysis are necessary to elevate our automated models to cosmography grade. Nevertheless, our pipeline enables the quick selection of lenses for follow-up and further modeling, which significantly speeds up the construction of cosmography-grade models. This important step forward will help us to take advantage of the increase in the number of lenses that is expected in the coming decade, which is an increase of several orders of magnitude.

    more » « less
  2. null (Ed.)
    ABSTRACT In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant H0. However, published state-of-the-art analyses require of order 1 yr of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge, we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated data sets. The results in Rungs 1 and 2 show that methods that use only the point source positions tend to have lower precision ($10\!-\!20{{\ \rm per\ cent}}$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic data sets can recover H0 within the target accuracy (|A| < 2 per cent) and precision (<6 per cent per system), even in the presence of a poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges. 
    more » « less
  3. Abstract We present a measurement of the Hubble constant (H0) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat ΛCDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}~\mathrm{km~s^{-1}~Mpc^{-1}}$, a $2.4{{\ \rm per\ cent}}$ precision measurement, in agreement with local measurements of H0 from type Ia supernovae calibrated by the distance ladder, but in 3.1σ tension with Planck observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in 5.3σ tension with Planck CMB determinations of H0 in flat ΛCDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat ΛCDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from Planck, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with Planck. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our H0 inference to cosmological model assumptions. For six different cosmological models, our combined inference on H0 ranges from ∼73–78 km s−1 Mpc−1, which is consistent with the local distance ladder constraints. 
    more » « less
  4. null (Ed.)
    Strong lensing time delays can measure the Hubble constant H 0 independently of any other probe. Assuming commonly used forms for the radial mass density profile of the lenses, a 2% precision has been achieved with seven Time-Delay Cosmography (TDCOSMO) lenses, in tension with the H 0 from the cosmic microwave background. However, without assumptions on the radial mass density profile – and relying exclusively on stellar kinematics to break the mass-sheet degeneracy – the precision drops to 8% with the current data obtained using the seven TDCOSMO lenses, which is insufficient to resolve the H 0 tension. With the addition of external information from 33 Sloan Lens ACS (SLACS) lenses, the precision improves to 5% if the deflectors of TDCOSMO and SLACS lenses are drawn from the same population. We investigate the prospect of improving the precision of time-delay cosmography without relying on mass profile assumptions to break the mass-sheet degeneracy. Our forecasts are based on a previously published hierarchical framework. With existing samples and technology, 3.3% precision on H 0 can be reached by adding spatially resolved kinematics of the seven TDCOSMO lenses. The precision improves to 2.5% with the further addition of kinematics for 50 nontime-delay lenses from SLACS and the Strong Lensing Legacy Survey. Expanding the samples to 40 time-delay and 200 nontime-delay lenses will improve the precision to 1.5% and 1.2%, respectively. Time-delay cosmography can reach sufficient precision to resolve the Hubble tension at 3–5 σ , without assumptions on the radial mass profile of lens galaxies. By obtaining this precision with and without external datasets, we will test the consistency of the samples and enable further improvements based on even larger future samples of time-delay and nontime-delay lenses (e.g., from the Rubin , Euclid , and Roman Observatories). 
    more » « less

    Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilizing automated techniques and human visual classification in conjunction with realistic mock images produced using the NewHorizon cosmological simulation, we investigate the nature, frequency, and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-yr depth of the Legacy Survey of Space and Time (30–31 mag arcsec−2), falling to 60 per cent assuming a shallower final depth of 29.5 mag arcsec−2. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M⋆ ∼ 1011.5 M⊙). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterization of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimized, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z < 0.2).

    more » « less