SLITronomy: towards a fully wavelet-based strong lensing inversion technique

Abstract

Strong gravitational lensing provides a wealth of astrophysical information on the baryonic and dark matter content of galaxies. It also serves as a valuable cosmological probe by allowing us to measure the Hubble constant independently of other methods. These applications all require the difficult task of inverting the lens equation and simultaneously reconstructing the mass profile of the lens along with the original light profile of the unlensed source. As there is no reason for either the lens or the source to be simple, we need methods that both invert the lens equation with a large number of degrees of freedom and also enforce a well-controlled regularisation that avoids the appearance of spurious structures. This can be beautifully accomplished by representing signals in wavelet space. Building on the Sparse Lens Inversion Technique (SLIT), in this work we present an improved sparsity-based method that describes lensed sources using wavelets and optimises over the parameters given an analytical lens mass profile. We apply our technique on simulated HST and E-ELT data, as well as on real HST images of lenses from the Sloan Lens ACS (SLACS) sample, assuming a lens model. We show that wavelets allow us to reconstruct lensed sources containing detailed substructures when using both present-day data and high-resolution images from future thirty-meter-class telescopes. Wavelets moreover provide a much more tractable solution in terms of quality and computation time compared to using a source model that combines smooth analytical profiles and shapelets. Requiring very little human interaction, our pixel-based technique fits into the effort to devise automated modelling schemes. It can be incorporated in the standard workflow of sampling analytical lens model parameters. The method, which we call SLITronomy, is freely available as a new plug-in to the modelling software Lenstronomy.

References in corpus (73)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Astropy: A Community Python Package for Astronomy
• The NumPy array: a structure for efficient numerical computation
• Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions
• H0LiCOW XIII. A 2.4% measurement of $H_0$ from lensed quasars: $5.3\sigma$ tension between early and late-Universe probes
• The Sloan Lens ACS Survey. X. Stellar, Dynamical, and Total Mass Correlations of Massive Early-type Galaxies
• The Sloan Lens ACS Survey. I. A Large Spectroscopically Selected Sample of Massive Early-Type Lens Galaxies
• Shapes and Shears, Stars and Smears: Optimal Measurements for Weak Lensing
• The Sloan Lens ACS Survey. IV: the mass density profile of early-type galaxies out to 100 effective radii
• Massive Dark-matter halos and Evolution of Early-type Galaxies to z=1
• H0LiCOW - IX. Cosmographic analysis of the doubly imaged quasar SDSS 1206+4332 and a new measurement of the Hubble constant
• TDCOSMO IV: Hierarchical time-delay cosmography -- joint inference of the Hubble constant and galaxy density profiles
• Two accurate time-delay distances from strong lensing: Implications for cosmology
• Shapelets: I. A Method for Image Analysis
• Detection of lensing substructure using ALMA observations of the dusty galaxy SDP.81
• Lenstronomy: multi-purpose gravitational lens modelling software package
• Detection of a Dark Substructure through Gravitational Imaging
• Bayesian Strong Gravitational-Lens Modeling on Adaptive Grids: Objective Detection of Mass Substructure in Galaxies
• Fast Automated Analysis of Strong Gravitational Lenses with Convolutional Neural Networks
• A Bayesian analysis of regularised source inversions in gravitational lensing
• Mass-sheet degeneracy, power-law models and external convergence: Impact on the determination of the Hubble constant from gravitational lensing
• STRIDES: a 3.9 per cent measurement of the Hubble constant from the strong lens system DES J0408-5354
• Dark Matter Halo Profiles of Massive Clusters: Theory vs. Observations
• Semi-linear gravitational lens inversion
• The BOSS Emission-Line Lens Survey (BELLS). I. A large spectroscopically selected sample of Lens Galaxies at redshift ~ 0.5
• The Halos of Satellite Galaxies: the Companion of the Massive Elliptical Lens SL2S J08544-0121
• Gravitational lens modeling with basis sets
• Uncertainties in Parameters Estimated with Neural Networks: Application to Strong Gravitational Lensing
• AutoLens: Automated Modeling of a Strong Lens's Light, Mass and Source
• Is every strong lens model unhappy in its own way? Uniform modelling of a sample of 13 quadruply+ imaged quasars
• Lensing Degeneracies Revisited
• Low-mass halo perturbations in strong gravitational lenses at redshift z$\sim$0.5 are consistent with CDM
• Nuisance hardened data compression for fast likelihood-free inference
• Lensing substructure quantification in RXJ1131-1231: A 2 keV lower bound on dark matter thermal relic mass
• Dissecting the Gravitational Lens B1608+656. I. Lens Potential Reconstruction
• Decomposition of the visible and dark matter in the Einstein ring 0047-2808 by semi-linear inversion
• Adaptive Semi-linear Inversion of Strong Gravitational Lens Imaging
• Models of the Cosmic Horseshoe Gravitational Lens
• Lens galaxies in the Illustris simulation: power-law models and the bias of the Hubble constant from time-delays
• Galaxy and Mass Assembly (GAMA): The halo mass of galaxy groups from maximum-likelihood weak lensing
• The SL2S Galaxy-scale Gravitational Lens Sample. I. The alignment of mass and light in massive early-type galaxies at z=0.2-0.9
• Substructure in lensing clusters and simulations
• SZ and CMB reconstruction using Generalized Morphological Component Analysis
• Simulations of Strong Gravitational Lensing with Substructure
• Data-Driven Reconstruction of Gravitationally Lensed Galaxies using Recurrent Inference Machines
• Could quasar lensing time delays hint to a core component in halos, instead of H0 tension?
• Robust sparse image reconstruction of radio interferometric observations with purify
• High Resolution Weak Lensing Mass-Mapping Combining Shear and Flexion
• LensMEM: A Gravitational Lens Inversion Algorithm Using the Maximum Entropy Method
• On the Choice of Lens Density Profile in Time Delay Cosmography
• Deep Learning for space-variant deconvolution in galaxy surveys
• Limitations for shapelet-based weak-lensing measurements
• Time Delay Lens Modelling Challenge
• The use of convolutional neural networks for modelling large optically-selected strong galaxy-lens samples
• A novel 3D technique to study the kinematics of lensed galaxies
• Models of the giant quadruple quasar SDSS J1004+4112
• A novel approach to visibility-space modelling of interferometric gravitational lens observations at high angular resolution
• TDCOSMO III: Dark matter substructure meets dark energy -- the effects of (sub)halos on strong-lensing measurements of $H_0$
• A Compact Early-type Galaxy at z = 0.6 Under a Magnifying Lens: Evidence For Inside-out Growth
• Statistical and systematic uncertainties in pixel-based source reconstruction algorithms for gravitational lensing
• Differentiable Strong Lensing: Uniting Gravity and Neural Nets through Differentiable Probabilistic Programming
• An E-ELT Case Study: Colour-Magnitude Diagrams of an Old Galaxy in the Virgo Cluster
• Multi-band morpho-Spectral Component Analysis Deblending Tool (MuSCADeT): Deblending colourful objects
• PySAP: Python Sparse Data Analysis Package for Multidisciplinary Image Processing
• Mapping UV Properties Throughout the Cosmic Horseshoe: Lessons from VLT-MUSE
• Testing the Evolution of the Correlations between Supermassive Black Holes and their Host Galaxies using Eight Strongly Lensed Quasars
• On the use of shapelets in modelling resolved, gravitationally lensed images
• Sparse Lens Inversion Technique (SLIT): lens and source separability from linear inversion of the source reconstruction problem
• Joint Multichannel Deconvolution and Blind Source Separation
• Sparse reconstruction of the merging A520 cluster system
• Precision cosmology with time delay lenses: high resolution imaging requirements
• The Einstein Ring 0047-2808 Revisited: A Bayesian Inversion
• SimCADO - an instrument data simulator package for MICADO at the E-ELT

Cited by in corpus (20)

• lenstronomy II: A gravitational lensing software ecosystem
• Strong lensing time-delay cosmography in the 2020s
• Automated galaxy-galaxy strong lens modelling: no lens left behind
• A lensed radio jet at milli-arcsecond resolution I: Bayesian comparison of parametric lens models
• TDCOSMO. VII. Boxyness/discyness in lensing galaxies : Detectability and impact on $H_0$
• Beyond the bulge-halo conspiracy? Density profiles of Early-type galaxies from extended-source strong lensing
• GIGA-Lens: Fast Bayesian Inference for Strong Gravitational Lens Modeling
• Using wavelets to capture deviations from smoothness in galaxy-scale strong lenses
• Strong Lensing by Galaxies
• Estimating the warm dark matter mass from strong lensing images with truncated marginal neural ratio estimation
• The Very Knotty Lenser: exploring the role of regularization in source and potential reconstructions using Gaussian Process Regression
• Strong-lensing source reconstruction with variationally optimised Gaussian processes
• Pixelated Reconstruction of Foreground Density and Background Surface Brightness in Gravitational Lensing Systems using Recurrent Inference Machines
• Modeling lens potentials with continuous neural fields in galaxy-scale strong lenses
• Exploiting the diversity of modeling methods to probe systematic biases in strong lensing analyses
• El Gordo needs El Anzuelo: Probing the structure of cluster members with multi-band extended arcs in JWST data
• Strong gravitational lensing by AGNs as a probe of the quasar-host relations in the distant Universe
• Lightweight HI source finding for next generation radio surveys
• Measuring the Substructure Mass Power Spectrum of 23 SLACS Strong Galaxy-Galaxy Lenses with Convolutional Neural Networks
• Resolving star-forming clumps in a z $\sim$ 2 lensed galaxy: a pixelated Bayesian approach