Abstract:Strongly lensed quasars can be used to constrain cosmological parameters through time-delay cosmography. Models of the lens masses are a necessary component of this analysis. To enable time-delay cosmography from a sample of $\mathcal{O}(10^3)$ lenses, which will soon become available from surveys like the Rubin Observatory's Legacy Survey of Space and Time (LSST) and the Euclid Wide Survey, we require fast and standardizable modeling techniques. To address this need, we apply neural posterior estimation (NPE) for modeling galaxy-scale strongly lensed quasars from the Strong Lensing Insights into the Dark Energy Survey (STRIDES) sample. NPE brings two advantages: speed and the ability to implicitly marginalize over nuisance parameters. We extend this method by employing sequential NPE to increase precision of mass model posteriors. We then fold individual lens models into a hierarchical Bayesian inference to recover the population distribution of lens mass parameters, accounting for out-of-distribution shift. After verifying our method using simulated analogs of the STRIDES lens sample, we apply our method to 14 Hubble Space Telescope single-filter observations. We find the population mean of the power-law elliptical mass distribution slope, $\gamma_{\text{lens}}$, to be $\mathcal{M}_{\gamma_{\text{lens}}}=2.13 \pm 0.06$. Our result represents the first population-level constraint for these systems. This population-level inference from fully automated modeling is an important stepping stone towards cosmological inference with large samples of strongly lensed quasars.

What problem does this paper attempt to address?

### What problems does this paper attempt to solve? This paper aims to solve the problems of large - scale modeling and parameter estimation in strong gravitational lensing systems (especially strong - gravitational - lensing quasars) to support the application of time - delay cosmology. Specifically, researchers hope to achieve fast and automated modeling of large - scale strong - gravitational - lensing samples through Neural Posterior Estimation (NPE) and Sequential Neural Posterior Estimation (SNPE) techniques. These problems include: 1. **Increasing the modeling speed**: Traditional methods require a large amount of computing resources and time when dealing with a single lens (about several thousand CPU hours per lens), and upcoming large - scale survey projects such as LSST will provide tens of thousands of strong - gravitational - lensing samples. Therefore, faster modeling methods must be developed. 2. **Reducing human intervention**: Traditional methods rely on a large amount of manual adjustment and verification work, which is unrealistic for dealing with large - scale data sets. An automated modeling process can significantly reduce the need for human intervention. 3. **Implicitly marginalizing nuisance parameters**: NPE and SNPE can implicitly handle nuisance parameters in the model, thereby simplifying the model construction process and improving the accuracy of parameter estimation. 4. **Population - level inference**: By combining Bayesian hierarchical inference, researchers hope to infer the distribution of lens mass parameters at the population level from multiple lenses, which will help to more accurately constrain cosmological parameters (such as the Hubble constant $ H_0 $ and the dark - energy equation - of - state parameters $ w_0, w_a $). 5. **Verifying the effectiveness of new techniques**: Researchers also hope to verify the effectiveness and reliability of NPE and SNPE techniques through simulated data and actual observational data (such as single - band observations by the Hubble Space Telescope). In summary, the goal of this paper is to develop an efficient, fully - automated lens - modeling pipeline so that accurate measurements of cosmological parameters can be made using large - scale strong - gravitational - lensing samples in the future.

Lens Modeling of STRIDES Strongly Lensed Quasars using Neural Posterior Estimation

Neural network prediction of model parameters for strong lensing samples from Hyper Suprime-Cam Survey

Modeling lens potentials with continuous neural fields in galaxy-scale strong lenses

Domain-Adaptive Neural Posterior Estimation for Strong Gravitational Lens Analysis

A model for galaxy-galaxy strong lensing statistics in surveys

COOL-LAMPS. VII. Quantifying Strong-lens Scaling Relations with 177 Cluster-scale Gravitational Lenses in DECaLS

Unveiling lens light complexity with a novel multi-Gaussian expansion approach for strong gravitational lensing

A model for galaxy–galaxy strong lensing statistics in surveys

Quantifying Environmental and Line-of-Sight Effects in Models of Strong Gravitational Lens Systems

A Strong Gravitational Lens Is Worth a Thousand Dark Matter Halos: Inference on Small-Scale Structure Using Sequential Methods

Elucidating Galaxy Population Properties Using a Model-Free Analysis of Quadruply Imaged Quasar Lenses From Large Surveys

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses VIII. Deconvolution of high resolution near-IR images and simple mass models for 7 gravitationally lensed quasars

Cosmology using Strong Gravitational Lensing

Gravitational lens modelling in a citizen science context

Exploiting the diversity of modeling methods to probe systematic biases in strong lensing analyses

Measuring the Substructure Mass Power Spectrum of 23 SLACS Strong Galaxy-Galaxy Lenses with Convolutional Neural Networks

A Bayesian Approach to Strong Lens Finding in the Era of Wide-area Surveys

Mass Reconstruction of Galaxy-scale Strong Gravitational Lenses Using a Broken Power-law Model

Accelerating galaxy dynamical modeling using a neural network for joint lensing and kinematics analyses