One never walks alone: the effect of the perturber population on subhalo measurements in strong gravitational lenses

Abstract

Analyses of extended arcs in strong gravitational lensing images to date have constrained the properties of dark matter by measuring the parameters of one or two individual subhalos. However, since such analyses are reliant on likelihood-based methods like Markov-chain Monte Carlo or nested sampling, they require various compromises to the realism of lensing models for the sake of computational tractability, such as ignoring the numerous other subhalos and line-of-sight halos in the system, assuming a particular form for the source model and requiring the noise to have a known likelihood function. Here we show that a simulation-based inference method calledTruncated Marginal Neural Ratio Estimation (TMNRE) makes it possible to relax these requirements by training neural networks to directly compute marginal posteriors for subhalo parameters from lensing images. By performing a set of inference tasks on mock data, we verify the accuracy of TMNRE and show it can compute posteriors for subhalo parameters marginalized over populations of hundreds of substructures, as well as lens and source uncertainties. We also find the \gls*{mlp} Mixer network works far better for such tasks than the convolutional architectures explored in other lensing analyses. Furthermore, we show that since \gls*{tmnre} learns a posterior function it enables direct statistical checks that would be extremely expensive with likelihood-based methods. Our results show that TMNRE is well-suited for analyzing complex lensing data, and that the full subhalo and line-of-sight halo population must be included when measuring the properties of individual dark matter substructures with this technique.

13 pages, 8 figures - v2: version published by MNRAS

References in corpus (21)

• Array Programming with NumPy
• An Image is Worth 16x16 Words: Transformers for Image Recognition at Scale
• MLP-Mixer: An all-MLP Architecture for Vision
• Bayesian Strong Gravitational-Lens Modeling on Adaptive Grids: Objective Detection of Mass Substructure in Galaxies
• Analytic models of plausible gravitational lens potentials
• Flexible statistical inference for mechanistic models of neural dynamics
• Mining for Dark Matter Substructure: Inferring subhalo population properties from strong lenses with machine learning
• Fast and Credible Likelihood-Free Cosmology with Truncated Marginal Neural Ratio Estimation
• On the Power Spectrum of Dark Matter Substructure in Strong Gravitational Lenses
• Line-of-sight effects in strong gravitational lensing
• From Images to Dark Matter: End-To-End Inference of Substructure From Hundreds of Strong Gravitational Lenses
• A Trust Crisis In Simulation-Based Inference? Your Posterior Approximations Can Be Unfaithful
• Quantifying the Line-of-Sight Halo Contribution to the Dark Matter Convergence Power Spectrum from Strong Gravitational Lenses
• GIGA-Lens: Fast Bayesian Inference for Strong Gravitational Lens Modeling
• Using wavelets to capture deviations from smoothness in galaxy-scale strong lenses
• Estimating the warm dark matter mass from strong lensing images with truncated marginal neural ratio estimation
• Extracting the Subhalo Mass Function from Strong Lens Images with Image Segmentation
• Image segmentation for analyzing galaxy-galaxy strong lensing systems
• Strong-lensing source reconstruction with variationally optimised Gaussian processes
• Inferring subhalo effective density slopes from strong lensing observations with neural likelihood-ratio estimation
• Scanning For Dark Matter Subhalos in Hubble Space Telescope Imaging of 54 Strong Lenses

Cited by in corpus (11)

• Strong Lensing by Galaxies
• Scalable inference with Autoregressive Neural Ratio Estimation
• Simulation-Based Inference of the sky-averaged 21-cm signal from CD-EoR with REACH
• Fast likelihood-free inference in the LSS Stage IV era
• SIDE-real: Supernova Ia Dust Extinction with truncated marginal neural ratio estimation applied to real data
• Robustness of Neural Ratio and Posterior Estimators to Distributional Shifts for Population-Level Dark Matter Analysis in Strong Gravitational Lensing
• Detecting the hidden population of low-mass haloes in strong lenses
• Observable Optimization for Precision Theory: Machine Learning Energy Correlators
• Deep inference of simulated strong lenses in ground-based surveys
• LSST Strong Lensing Systems Dark Matter Sensitivity Analysis with Neural Ratio Estimators
• Quantifying Weighted Morphological Content of Large-Scale Structures via Simulation-Based Inference