Engineering Papers⌕ Search

Engineering topics

Birrer, S.

Publications and source records attributed to Birrer, S..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Strong Lensing by Galaxies

Strong gravitational lensing at the galaxy scale is a valuable tool for various applications in astrophysics and cosmology. Some of the primary uses of galaxy-scale lensing are to study elliptical galaxies’ mass structure and evolution, constrain the stellar initial mass function, and measure cosmological parameters. Since the discovery of the first galaxy-scale lens in the 1980s, this field has made significant advancements in data quality and modeling techniques. In this review, we describe the most common methods for modeling lensing observables, especially imaging data, as they are the most accessible and informative source of lensing observables. We then summarize the primary findings from the literature on the astrophysical and cosmological applications of galaxy-scale lenses. We also discuss the current limitations of the data and methodologies and provide an outlook on the expected improvements in both areas in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

JWST lensed quasar dark matter survey – II. Strongest gravitational lensing limit on the dark matter free streaming length to date

This is the second in a series of papers in which we use JWST Mid Infrared Instrument multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of nine systems. The warm dust region is compact and sensitive to perturbations by populations of haloes down to masses $\sim 10^6$ M$_{\odot }$. Using these warm dust flux-ratio measurements in combination with five previous narrow-line flux-ratio measurements, we constrain the halo mass function. In our model, we allow for complex deflector macromodels with flexible third- and fourth-order multipole deviations from ellipticity, and we introduce an improved model of the tidal evolution of subhaloes. We constrain a WDM model and find an upper limit on the half-mode mass of $10^{7.6}\, {\rm M}_\odot$ at posterior odds of 10:1. This corresponds to a lower limit on a thermally produced dark matter particle mass of 6.1 keV. This is the strongest gravitational lensing constraint to date, and comparable to those from independent probes such as the Ly $\alpha$ forest and Milky Way satellite galaxies.

dark matter↗
DOE OSTI2024

JWST lensed quasar dark matter survey – I. Description and first results

ABSTRACT The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar’s light is emitted from a physical region of the quasar with an angular scale of milliarcseconds or larger; this minimizes microlensing effects by stars within the deflector. The warm dust region of quasars fits this criterion, as it has parsec-size physical scales and dominates the spectral energy distribution of quasars at wavelengths greater than 10 μm. The JWST Mid-Infrared Instrument is adept at detecting redshifted light in this wavelength range, offering both the spatial resolution and sensitivity required for accurate gravitational lensing flux ratio measurements. Here, we introduce our survey designed to measure the warm dust flux ratios of 31 lensed quasars. We discuss the flux-ratio measurement technique and present results for the first target, DES J0405-3308. We find that we can measure the quasar warm dust flux ratios with 3 per cent precision. Our simulations suggest that this precision makes it feasible to detect the presence of 107 M⊙ dark matter haloes at cosmological distances. Such haloes are expected to be completely dark in cold dark matter models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

LensWatch. I. Resolved HST Observations and Constraints on the Strongly Lensed Type Ia Supernova 2022qmx (“SN Zwicky”)

Supernovae (SNe) that have been multiply imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply imaged Type Ia SN 2022qmx (aka, "SN Zwicky") at z = 0.3544. SN Zwicky was discovered by the Zwicky Transient Facility in spatially unresolved data. Here we present follow-up Hubble Space Telescope observations of SN Zwicky, the first from the multicycle "LensWatch (www.lenswatch.org)" program. We measure photometry for each of the four images of SN Zwicky, which are resolved in three WFC3/UVIS filters (F475W, F625W, and F814W) but unresolved with WFC3/IR F160W, and present an analysis of the lensing system using a variety of independent lens modeling methods. We find consistency between lens-model-predicted time delays (≲1 day), and delays estimated with the single epoch of Hubble Space Telescope colors (≲3.5 days), including the uncertainty from chromatic microlensing (~1–1.5 days). Our lens models converge to an Einstein radius of ${\theta }_{{\rm{E}}}=({0.168}_{-0.005}^{+0.009})^{\prime\prime} $, the smallest yet seen in a lensed SN system. The "standard candle" nature of SN Zwicky provides magnification estimates independent of the lens modeling that are brighter than predicted by $\sim {1.7}_{-0.6}^{+0.8}$ mag and $\sim {0.9}_{-0.6}^{+0.8}$ mag for two of the four images, suggesting significant microlensing and/or additional substructure beyond the flexibility of our image-position mass models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Dark Energy Survey Year 3 results: Constraints on extensions to Λ CDM with weak lensing and galaxy clustering

We constrain six possible extensions to the Λ cold dark matter (CDM) model using measurements from the Dark Energy Survey’s first three years of observations, alone and in combination with external cosmological probes. The DES data are the two-point correlation functions of weak gravitational lensing, galaxy clustering, and their cross-correlation. We use simulated data vectors and blind analyses of real data to validate the robustness of our results to astrophysical and modeling systematic errors. In many cases, constraining power is limited by the absence of theoretical predictions beyond the linear regime that are reliable at our required precision. The Λ CDM extensions are dark energy with a time-dependent equation of state, nonzero spatial curvature, additional relativistic degrees of freedom, sterile neutrinos with eV-scale mass, modifications of gravitational physics, and a binned σ 8 ( z ) model which serves as a phenomenological probe of structure growth. For the time-varying dark energy equation of state evaluated at the pivot redshift we find ( w p , w a ) = ( - 0.9 9 - 0.17 + 0.28 , - 0.9 ± 1.2 ) at 68% confidence with z p = 0.24 from the DES measurements alone, and ( w p , w a ) = ( - 1.0 3 - 0.03 + 0.04 , - 0. 4 - 0.3 + 0.4 ) with z p = 0.21 for the combination of all data considered. Curvature constraints of Ω k = 0.0009 ± 0.0017 and effective relativistic species N eff = 3.1 0 - 0.16 + 0.15 are dominated by external data, though adding DES information to external low-redshift probes tightens the Ω k constraints that can be made without cosmic microwave background observables by 20%. For massive sterile neutrinos, DES combined with external data improves the upper bound on the mass m eff by a factor of 3 compared to previous analyses, giving 95% limits of ( Δ N eff , m eff ) ≤ ( 0.28 , 0.20 eV ) when using priors matching a comparable Planck analysis. For modified gravity, we constrain changes to the lensing and Poisson equations controlled by functions Σ ( k , z ) = Σ 0 Ω Λ ( z ) / Ω Λ , 0 and μ ( k , z ) = μ 0 Ω Λ ( z ) / Ω Λ , 0 , respectively, to Σ 0 = 0. 6 - 0.5 + 0.4 from DES alone and ( Σ 0 , μ 0 ) = ( 0.04 ± 0.05 , 0.0 8 - 0.19 + 0.21 ) for the combination of all data, both at 68% confidence. Overall, we find no significant evidence for physics beyond Λ CDM .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

TDCOSMO. X. Automated modeling of nine strongly lensed quasars and comparison between lens-modeling software

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 with GLEE, 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 the Hubble Space 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 software Lenstronomy revealed 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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Early Results from GLASS-JWST. XI. Stellar Masses and Mass-to-light Ratio of z > 7 Galaxies

We exploit James Webb Space Telescope (JWST) NIRCam observations from the GLASS-JWST-Early Release Science program to investigate galaxy stellar masses at z > 7. We first show that JWST observations reduce the uncertainties on the stellar mass by a factor of at least 5–10, when compared with the highest-quality data sets available to date. We then study the UV mass-to-light ratio, finding that galaxies exhibit a a two orders of magnitude range of M / L UV values for a given luminosity, indicative of a broad variety of physical conditions and star formation histories. As a consequence, previous estimates of the cosmic stellar-mass density—based on an average correlation between UV luminosity and stellar mass—can be biased by as much as a factor of ~6. Our first exploration demonstrates that JWST represents a new era in our understanding of stellar masses at z > 7 and, therefore, of the growth of galaxies prior to cosmic reionization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Early Results from GLASS-JWST. VII. Evidence for Lensed, Gravitationally Bound Protoglobular Clusters at z = 4 in the Hubble Frontier Field A2744

We investigate the blue and optical rest-frame sizes (λ ≃ 2300–4000 Å) of three compact star-forming regions in a galaxy at z = 4 strongly lensed (×30, ×45, and ×100) by the Hubble Frontier Field galaxy cluster A2744 using GLASS-ERS James Webb Space Telescope (JWST)/NIRISS imaging at 1.15 μm, 1.50 μm, and 2.0 μm with a point-spread function ≲0$_{.}^{"}$1. In particular, the Balmer break is probed in detail for all multiply imaged sources of the system. With ages of a few tens of Myr, stellar masses in the range (0.7–4.0) ×10 6 M⊙ and optical/ultraviolet effective radii spanning the interval 3 < R eff < 20 pc, such objects are currently the highest-redshift (spectroscopically confirmed) gravitationally bound young massive star clusters (YMCs), with stellar mass surface densities resembling those of local globular clusters. Optical (4000 Å, JWST-based) and ultraviolet (1600 Å, Hubble Space Telescope–based) sizes are fully compatible. The contribution to the ultraviolet underlying continuum emission (1600 Å) is ~30%, which decreases by a factor of 2 in the optical for two of the YMCs (~4000 Å rest-frame), reflecting the young ages (<30 Myr) inferred from the spectral energy distribution fitting and supported by the presence of high-ionization lines secured with the Very Large Telescope/MUSE. Such bursty forming regions enhance the specific star formation rate of the galaxy, which is ≃10 Gyr –1 . This galaxy would be among the extreme analogs observed in the local universe having a high star formation rate surface density and a high occurrence of massive stellar clusters in formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

TDCOSMO. IX. Systematic comparison between lens modelling software programs: time delay prediction for WGD 2038-4008

The importance of alternative methods for measuring the Hubble constant, such as time-delay cosmography, is highlighted by the recent Hubble tension. It is paramount to thoroughly investigate and rule out systematic biases in all measurement methods before we can accept new physics as the source of this tension. In this study, we perform a check for systematic biases in the lens modelling procedure of time-delay cosmography by comparing independent and blind time-delay predictions of the system WGD 2038-4008 from two teams using two different software programs: GLEE and LENSTRONOMY. The predicted time delays from the two teams incorporate the stellar kinematics of the deflector and the external convergence from line-of-sight structures. The un-blinded time-delay predictions from the two teams agree within 1.2σ, implying that once the time delay is measured the inferred Hubble constant will also be mutually consistent. However, there is a ~4σ discrepancy between the power-law model slope and external shear, which is a significant discrepancy at the level of lens models before the stellar kinematics and the external convergence are incorporated. We identify the difference in the reconstructed point spread function (PSF) to be the source of this discrepancy. When the same reconstructed PSF was used by both teams, we achieved excellent agreement, within ~0.6σ<, indicating that potential systematics stemming from source reconstruction algorithms and investigator choices are well under control. We recommend that future studies supersample the PSF as needed and marginalize over multiple algorithms or realizations for the PSF reconstruction to mitigate the systematics associated with the PSF. A future study will measure the time delays of the system WGD 2038-4008 and infer the Hubble constant based on our mass models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

STRIDES: automated uniform models for 30 quadruply imaged quasars

ABSTRACT 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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Early Results from GLASS-JWST. V: The First Rest-frame Optical Size–Luminosity Relation of Galaxies at z > 7

We present the first rest-frame optical size–luminosity relation of galaxies at z > 7, using the NIRCam imaging data obtained by the GLASS James Webb Space Telescope Early Release Science (GLASS-JWST-ERS) program, providing the deepest extragalactic data of the ERS campaign. Our sample consists of 19 photometrically selected bright galaxies with m F444W ≤ 27.8 at 7 < z < 9 and m F444W < 28.2 at z ~ 9–15. We measure the size of the galaxies in five bands, from rest-frame optical (~4800 Å) to the UV (~1600 Å) based on the Sérsic model, and analyse the size–luminosity relation as a function of wavelength. Remarkably, the data quality of the NIRCam imaging is sufficient to probe the half-light radius r e down to ~100 pc at z > 7. Given the limited sample size and magnitude range, we first fix the slope to that observed for larger samples in rest-frame UV using Hubble Space Telescope samples. The median size r 0 at the reference luminosity M = –21 decreases slightly from rest-frame optical (600 ± 80 pc) to UV (450 ± 130 pc). We then refit the size–luminosity relation allowing the slope to vary. The slope is consistent with β ~ 0.2 for all bands except F150W, where we find a marginally steeper slope of β = 0.53 ± 0.15. The steep UV slope is mainly driven by the smallest and faintest galaxies. If confirmed by larger samples, it implies that the UV size–luminosity relation breaks toward the faint end, as suggested by lensing studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Finding quadruply imaged quasars with machine learning – I. Methods

ABSTRACT Strongly lensed quadruply imaged quasars (quads) are extraordinary objects. They are very rare in the sky and yet they provide unique information about a wide range of topics, including the expansion history and the composition of the Universe, the distribution of stars and dark matter in galaxies, the host galaxies of quasars, and the stellar initial mass function. Finding them in astronomical images is a classic ‘needle in a haystack’ problem, as they are outnumbered by other (contaminant) sources by many orders of magnitude. To solve this problem, we develop state-of-the-art deep learning methods and train them on realistic simulated quads based on real images of galaxies taken from the Dark Energy Survey, with realistic source and deflector models, including the chromatic effects of microlensing. The performance of the best methods on a mixture of simulated and real objects is excellent, yielding area under the receiver operating curve in the range of 0.86–0.89. Recall is close to 100 per cent down to total magnitude i ∼ 21 indicating high completeness, while precision declines from 85 per cent to 70 per cent in the range i ∼ 17–21. The methods are extremely fast: training on 2 million samples takes 20 h on a GPU machine, and 108 multiband cut-outs can be evaluated per GPU-hour. The speed and performance of the method pave the way to apply it to large samples of astronomical sources, bypassing the need for photometric pre-selection that is likely to be a major cause of incompleteness in current samples of known quads.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

DeepZipper: A Novel Deep-learning Architecture for Lensed Supernovae Identification

Large-scale astronomical surveys have the potential to capture data on large numbers of strongly gravitationally lensed supernovae (LSNe). To facilitate timely analysis and spectroscopic follow-up before the supernova fades, an LSN needs to be identified soon after it begins. To quickly identify LSNe in optical survey data sets, we designed ZipperNet, a multibranch deep neural network that combines convolutional layers (traditionally used for images) with long short-term memory layers (traditionally used for time series). We tested ZipperNet on the task of classifying objects from four categories—no lens, galaxy-galaxy lens, lensed Type-Ia supernova, lensed core-collapse supernova—within high-fidelity simulations of three cosmic survey data sets: the Dark Energy Survey, Rubin Observatory’s Legacy Survey of Space and Time (LSST), and a Dark Energy Spectroscopic Instrument (DESI) imaging survey. Among our results, we find that for the LSST-like data set, ZipperNet classifies LSNe with a receiver operating characteristic area under the curve of 0.97, predicts the spectroscopic type of the lensed supernovae with 79% accuracy, and demonstrates similarly high performance for LSNe 1–2 epochs after first detection. We anticipate that a model like ZipperNet, which simultaneously incorporates spatial and temporal information, can play a significant role in the rapid identification of lensed transient systems in cosmic survey experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Dark Energy Survey Year 3 results: Cosmological constraints from galaxy clustering and weak lensing

We present the first cosmology results from large-scale structure using the full 5000 deg 2 of imaging data from the Dark Energy Survey (DES) Data Release 1. We perform an analysis of large-scale structure combining three two-point correlation functions ( 3 × 2 pt ): (i) cosmic shear using 100 million source galaxies, (ii) galaxy clustering, and (iii) the cross-correlation of source galaxy shear with lens galaxy positions, galaxy–galaxy lensing. To achieve the cosmological precision enabled by these measurements has required updates to nearly every part of the analysis from DES Year 1, including the use of two independent galaxy clustering samples, modeling advances, and several novel improvements in the calibration of gravitational shear and photometric redshift inference. The analysis was performed under strict conditions to mitigate confirmation or observer bias; we describe specific changes made to the lens galaxy sample following unblinding of the results and tests of the robustness of our results to this decision. We model the data within the flat Λ CDM and w CDM cosmological models, marginalizing over 25 nuisance parameters. We find consistent cosmological results between the three two-point correlation functions; their combination yields clustering amplitude S 8 = 0.77 6 - 0.017 + 0.017 and matter density Ω m = 0.33 9 - 0.031 + 0.032 in Λ CDM , mean with 68% confidence limits; S 8 = 0.77 5 - 0.024 + 0.026 , Ω m = 0.35 2 - 0.041 + 0.035 , and dark energy equation-of-state parameter w = - 0.9 8 - 0.20 + 0.32 in w CDM . These constraints correspond to an improvement in signal-to-noise of the DES Year 3 3 × 2 pt data relative to DES Year 1 by a factor of 2.1, about 20% more than expected from the increase in observing area alone. This combination of DES data is consistent with the prediction of the model favored by the Planck 2018 cosmic microwave background (CMB) primary anisotropy data, which is quantified with a probability-to-exceed p = 0.13 –0.48. We find better agreement between DES 3 × 2 pt and Planck than in DES Y1, despite the significantly improved precision of both. When combining DES 3 × 2 pt data with available baryon acoustic oscillation, redshift-space distortion, and type Ia supernovae data, we find p = 0.34 . Combining all of these datasets with Planck CMB lensing yields joint parameter constraints of S 8 = 0.81 2 - 0.008 + 0.008 , Ω m = 0.30 6 - 0.005 + 0.004 , h = 0.68 0 - 0.003 + 0.004 , and ∑ m ν < 0.13 eV (95% C.L.) in Λ CDM ; S 8 = 0.81 2 - 0.008 + 0.008 , Ω m = 0.30 2 - 0.006 + 0.006 , h = 0.68 7 - 0.007 + 0.006 , and w = - 1.03 1 - 0.027 + 0.030 in w CDM .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Time delay lens modelling challenge

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.

Ding, X.↗
DOE OSTI2020

TDCOSMO - IV. Hierarchical time-delay cosmography – joint inference of the Hubble constant and galaxy density profiles

The H0LiCOW collaboration inferred via strong gravitational lensing time delays a Hubble constant value of H0 = 73.3−1.8+1.7 km s−1 Mpc−1, describing deflector mass density profiles by either a power-law or stars (constant mass-to-light ratio) plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in H0. We quantify any potential effect of the MST with a flexible family of mass models, which directly encodes it, and they are hence maximally degenerate with H0. Our calculation is based on a new hierarchical Bayesian approach in which the MST is only constrained by stellar kinematics. The approach is validated on mock lenses, which are generated from hydrodynamic simulations. We first applied the inference to the TDCOSMO sample of seven lenses, six of which are from H0LiCOW, and measured H0 = 74.5−6.1+5.6 km s−1 Mpc−1. Secondly, in order to further constrain the deflector mass density profiles, we added imaging and spectroscopy for a set of 33 strong gravitational lenses from the Sloan Lens ACS (SLACS) sample. For nine of the 33 SLAC lenses, we used resolved kinematics to constrain the stellar anisotropy. From the joint hierarchical analysis of the TDCOSMO+SLACS sample, we measured H0 = 67.4−3.2+4.1 km s−1 Mpc−1. This measurement assumes that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. Without relying on the form of the mass density profile used by H0LiCOW, we achieve a ∼5% measurement of H0. While our new hierarchical analysis does not statistically invalidate the mass profile assumptions by H0LiCOW – and thus the H0 measurement relying on them – it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on H0 derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data.Key words: gravitational lensing: strong / galaxies: general / galaxies: kinematics and dynamics / distance scale / cosmological parameters / cosmology: observations⋆ The full analysis is available at https://github.com/TDCOSMO/hierarchy_analysis_2020_public.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

STRIDES: Spectroscopic and photometric characterization of the environment and effects of mass along the line of sight to the gravitational lenses DES J0408–5354 and WGD 2038–4008

In time-delay cosmography, three of the key ingredients are (1) determining the velocity dispersion of the lensing galaxy, (2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and (3) estimating the external convergence κ_ext from less massive structures that are not included in the mass model. We present results on all three of these ingredients for two time-delay lensed quad quasar systems, DES J0408–5354 and WGD 2038–4008 . We use the Gemini, Magellan, and VLT telescopes to obtain spectra to both measure the stellar velocity dispersions of the main lensing galaxies and to identify the line-of-sight galaxies in these systems. Next, we identify 10 groups in DES J0408–5354 and two groups in WGD 2038–4008 using a group-finding algorithm. We then identify the most significant galaxy and galaxy-group perturbers using the ‘flexion shift’ criterion. We determine the probability distribution function of the external convergence κ_ext for both of these systems based on our spectroscopy and on the DES-only multiband wide-field observations. Using weighted galaxy counts, calibrated based on the Millennium Simulation, we find that DES J0408–5354 is located in a significantly underdense environment, leading to a tight (width |$\sim 3{{\ \rm per\ cent}}$|⁠), negative-value κ_ext distribution. On the other hand, WGD 2038–4008 is located in an environment of close to unit density, and its low source redshift results in a much tighter κ_ext of |$\sim 1{{\ \rm per\ cent}}$|⁠, as long as no external shear constraints are imposed.

WGD 2038-4008↗