-
Exploiting the diversity of modeling methods to probe systematic biases in strong lensing analyses
Authors:
A. Galan,
G. Vernardos,
Q. Minor,
D. Sluse,
L. Van de Vyvere,
M. Gomer
Abstract:
In the past decade, the diversity of strong lens modeling methods has exploded, from being purely analytical to pixelated and non-parametric, or based on deep learning. We embrace this diversity by selecting different software packages and use them to blindly model independently simulated Hubble Space Telescope imaging data. To overcome the difficulties arising from using different codes and conve…
▽ More
In the past decade, the diversity of strong lens modeling methods has exploded, from being purely analytical to pixelated and non-parametric, or based on deep learning. We embrace this diversity by selecting different software packages and use them to blindly model independently simulated Hubble Space Telescope imaging data. To overcome the difficulties arising from using different codes and conventions, we use the COde-independent Organized LEns STandard (COOLEST) to store, compare and release all models in a self-consistent and human-readable manner. From an ensemble of six modeling methods, we study the recovery of the lens potential parameters and properties of the reconstructed source. In particular, we simulate and infer parameters of an elliptical power-law mass distribution with external shear for the lens while each modeling method reconstructs the source differently. We find that overall, both lens and source properties are recovered reasonably well, but systematic biases arise in all methods. Interestingly, we do not observe that a single method is significantly more accurate than others, and the amount of bias largely depends on the specific lens or source property of interest. By combining posterior distributions from individual methods using equal weights, the maximal systematic biases on lens model parameters inferred from individual models are reduced by a factor of 5.4 on average. We investigate a selection of modeling effects that partly explain the observed biases, such as the cuspy nature of the background source and the accuracy of the point spread function. This work introduces, for the first time, a generic framework to compare and ease the combination of models obtained from different codes and methods, which will be key to retain accuracy in future strong lensing analyses.
△ Less
Submitted 12 June, 2024;
originally announced June 2024.
-
TDCOSMO. XVI. Measurement of the Hubble Constant from the Lensed Quasar WGD$\,$2038$-$4008
Authors:
Kenneth C. Wong,
Frédéric Dux,
Anowar J. Shajib,
Sherry H. Suyu,
Martin Millon,
Pritom Mozumdar,
Patrick R. Wells,
Adriano Agnello,
Simon Birrer,
Elizabeth J. Buckley-Geer,
Frédéric Courbin,
Christopher D. Fassnacht,
Joshua Frieman,
Aymeric Galan,
Huan Lin,
Philip J. Marshall,
Jason Poh,
Stefan Schuldt,
Dominique Sluse,
Tommaso Treu
Abstract:
Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant ($H_{0}$). The TDCOSMO collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038-4008 using new time-delay measurements and previous mass models by TDCOSMO. This is…
▽ More
Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant ($H_{0}$). The TDCOSMO collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038-4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain $D_{Δt} = 1.68^{+0.40}_{-0.38}$ Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat $Λ$CDM cosmology, we constrain the Hubble constant to be $H_{0} = 65^{+23}_{-14}\, \rm km\ s^{-1}\,Mpc^{-1}$. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory's Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.
△ Less
Submitted 5 July, 2024; v1 submitted 4 June, 2024;
originally announced June 2024.
-
Strong-lensing and kinematic analysis of CASSOWARY 31: can strong lensing constrain the masses of multi-plane lenses?
Authors:
H. Wang,
R. Canameras,
S. H. Suyu,
A. Galan,
C. Grillo,
G. B. Caminha,
L. Christensen
Abstract:
We present a mass measurement for the secondary lens along the line of sight (LoS) in the multi-plane strong lens modeling of the group-scale lens CASSOWARY 31 (CSWA 31). The secondary lens at redshift $z = 1.49$ is a spiral galaxy well aligned along the LoS with the main lens at $z = 0.683$. Using the MUSE integral-field spectroscopy of this spiral galaxy, we measure its rotation velocities and d…
▽ More
We present a mass measurement for the secondary lens along the line of sight (LoS) in the multi-plane strong lens modeling of the group-scale lens CASSOWARY 31 (CSWA 31). The secondary lens at redshift $z = 1.49$ is a spiral galaxy well aligned along the LoS with the main lens at $z = 0.683$. Using the MUSE integral-field spectroscopy of this spiral galaxy, we measure its rotation velocities and determine the mass from the gas kinematics. We compare the mass estimation of the secondary lens from the lensing models to the mass measurement from kinematics, finding that the predictions from strong lensing tend to be higher. By introducing an additional lens plane at $z = 1.36$ for an overdensity known to be present, we find a mass of $\simeq 10^{10}$ M$_\odot$ enclosed within 3.3 kpc from the centroid of the spiral galaxy, approaching the estimate from kinematics. This shows that secondary-lens mass measurements from multiple-plane modeling are affected by systematic uncertainties from the degeneracies between lens planes and the complex LoS structure. Conducting a detailed analysis of the LoS structures is therefore essential to improve the mass measurement of the secondary lens.
△ Less
Submitted 21 July, 2024; v1 submitted 19 April, 2024;
originally announced April 2024.
-
El Gordo needs El Anzuelo: Probing the structure of cluster members with multi-band extended arcs in JWST data
Authors:
A. Galan,
G. B. Caminha,
J. Knollmüller,
J. Roth,
S. H. Suyu
Abstract:
Gravitational lensing by galaxy clusters involves hundreds of galaxies over a large redshift range and increases the likelihood of rare phenomena (supernovae, microlensing, dark substructures, etc.). Characterizing the mass and light distributions of foreground and background objects often requires a combination of high-resolution data and advanced modeling techniques. We present the detailed anal…
▽ More
Gravitational lensing by galaxy clusters involves hundreds of galaxies over a large redshift range and increases the likelihood of rare phenomena (supernovae, microlensing, dark substructures, etc.). Characterizing the mass and light distributions of foreground and background objects often requires a combination of high-resolution data and advanced modeling techniques. We present the detailed analysis of El Anzuelo, a prominent quintuply imaged dusty star forming galaxy ($z_{\rm s}=2.29$), mainly lensed by three members of the massive galaxy cluster ACT-CL$\,$J0102$-$4915, also known as El Gordo ($z_{\rm d}=0.87$). We leverage JWST/NIRCam data containing previously unseen lensing features using a Bayesian, multi-wavelength, differentiable and GPU-accelerated modeling framework that combines Herculens (lens modeling) and NIFTy (field model and inference) software packages. For one of the deflectors, we complement lensing constraints with stellar kinematics measured from VLT/MUSE data. In our lens model, we explicitly include the mass distribution of the cluster, locally corrected by a constant shear field. We find that the two main deflectors (L1 and L2) have logarithmic mass density slopes steeper than isothermal, with $γ_{\rm L1} = 2.23\pm0.05$ and $γ_{\rm L2} = 2.21\pm0.04$. We argue that such steep density profiles can arise due to tidally truncated mass distributions, which we probe thanks to the cluster lensing boost and the strong asymmetry of the lensing configuration. Moreover, our three-dimensional source model captures most of the surface brightness of the lensed galaxy, revealing a clump of at most $400$ parsecs at the source redshift, visible at wavelengths $λ_{\rm rest}\gtrsim0.6$ $μ$m. Finally, we caution on using point-like features within extended arcs to constrain galaxy-scale lens models before securing them with extended arc modeling.
△ Less
Submitted 28 February, 2024;
originally announced February 2024.
-
Nine lensed quasars and quasar pairs discovered through spatially-extended variability in Pan-STARRS
Authors:
Frédéric Dux,
Cameron Lemon,
Frédéric Courbin,
Favio Neira,
Timo Anguita,
Aymeric Galan,
Sam Kim,
Maren Hempel,
Angela Hempel,
Régis Lachaume
Abstract:
We present the proof-of-concept of a method to find strongly lensed quasars using their spatially-extended photometric variability through difference imaging in cadenced imaging survey data. We apply the method to Pan-STARRS, starting with an initial selection of 14 107 Gaia multiplets with quasar-like infrared colours from WISE. We identify 229 candidates showing notable spatially-extended variab…
▽ More
We present the proof-of-concept of a method to find strongly lensed quasars using their spatially-extended photometric variability through difference imaging in cadenced imaging survey data. We apply the method to Pan-STARRS, starting with an initial selection of 14 107 Gaia multiplets with quasar-like infrared colours from WISE. We identify 229 candidates showing notable spatially-extended variability during the Pan-STARRS survey period. These include 20 known lenses, alongside an additional 12 promising candidates for which we obtain long-slit spectroscopy follow-up. This process results in the confirmation of four doubly lensed quasars, four unclassified quasar pairs and one projected quasar pair. Only three are pairs of stars or quasar+star projections, the false positive rate is thereby 25%. The lenses have separations between 0.81" and 1.24" and source redshifts between z = 1.47 and z = 2.46. Three of the unclassified quasar pairs are promising dual quasars candidates with separations ranging from 6.6 to 9.3 kpc. We expect that this technique will be a particularly efficient way to select lensed variables in the upcoming Rubin-LSST, which will be crucial given the expected limitations for spectroscopic follow-up
△ Less
Submitted 25 July, 2023;
originally announced July 2023.
-
Accelerating galaxy dynamical modeling using a neural network for joint lensing and kinematics analyses
Authors:
Matthew R. Gomer,
Sebastian Ertl,
Luca Biggio,
Han Wang,
Aymeric Galan,
Lyne Van de Vyvere,
Dominique Sluse,
Georgios Vernardos,
Sherry H. Suyu
Abstract:
Strong gravitational lensing is a powerful tool to provide constraints on galaxy mass distributions and cosmological parameters, such as the Hubble constant, $H_0$. Nevertheless, inference of such parameters from images of lensing systems is not trivial as parameter degeneracies can limit the precision in the measured lens mass and cosmological results. External information on the mass of the lens…
▽ More
Strong gravitational lensing is a powerful tool to provide constraints on galaxy mass distributions and cosmological parameters, such as the Hubble constant, $H_0$. Nevertheless, inference of such parameters from images of lensing systems is not trivial as parameter degeneracies can limit the precision in the measured lens mass and cosmological results. External information on the mass of the lens, in the form of kinematic measurements, is needed to ensure a precise and unbiased inference. Traditionally, such kinematic information has been included in the inference after the image modeling, using spherical Jeans approximations to match the measured velocity dispersion integrated within an aperture. However, as spatially resolved kinematic measurements become available via IFU data, more sophisticated dynamical modeling is necessary. Such kinematic modeling is expensive, and constitutes a computational bottleneck which we aim to overcome with our Stellar Kinematics Neural Network (SKiNN). SKiNN emulates axisymmetric modeling using a neural network, quickly synthesizing from a given mass model a kinematic map which can be compared to the observations to evaluate a likelihood. With a joint lensing plus kinematic framework, this likelihood constrains the mass model at the same time as the imaging data. We show that SKiNN's emulation of a kinematic map is accurate to considerably better precision than can be measured (better than $1\%$ in almost all cases). Using SKiNN speeds up the likelihood evaluation by a factor of $\sim 200$. This speedup makes dynamical modeling economical, and enables lens modelers to make effective use of modern data quality in the JWST era.
△ Less
Submitted 19 July, 2023;
originally announced July 2023.
-
Strong gravitational lensing by AGNs as a probe of the quasar-host relations in the distant Universe
Authors:
Martin Millon,
Frédéric Courbin,
Aymeric Galan,
Dominique Sluse,
Xuheng Ding,
Malte Tewes,
S. G. Djorgovski
Abstract:
The tight correlations found between masses of supermassive black holes (SMBHs) and the luminosities, total stellar masses, and velocity dispersions of their host galaxies are often interpreted as a sign of their co-evolution. Studying these correlations across redshift provides a powerful insight into the evolutionary path followed by the quasar and its host galaxy. While the mass of the black ho…
▽ More
The tight correlations found between masses of supermassive black holes (SMBHs) and the luminosities, total stellar masses, and velocity dispersions of their host galaxies are often interpreted as a sign of their co-evolution. Studying these correlations across redshift provides a powerful insight into the evolutionary path followed by the quasar and its host galaxy. While the mass of the black hole is accessible from single-epoch spectra, measuring the mass of its host galaxy is challenging as the active nucleus largely overshines its host. Here, we present a technique to probe quasar-host relations beyond the local universe with strong gravitational lensing, hence overcoming the use of stellar population models or velocity dispersion measurements, both prone to degeneracies. We study in detail one of the three known cases of strong lensing by a quasar to accurately measure the mass of its host and to infer a total lensing mass of $\log_{10}(M_{\rm Tot, h}/M_{\odot}) = 10.27^{+0.06}_{-0.07}~$ within the Einstein radius of 1.2 kpc. The lensing measurement is more precise than any other alternative technique and compatible with the local $M_{BH}$-$M_{\star, h}$ scaling relation. The sample of such quasar-galaxy or quasar-quasar lensing systems should reach a few hundreds with Euclid and Rubin-LSST, thus enabling the application of such a method with statistically significant sample sizes.
△ Less
Submitted 24 July, 2024; v1 submitted 1 June, 2023;
originally announced June 2023.
-
SKA Science Data Challenge 2: analysis and results
Authors:
P. Hartley,
A. Bonaldi,
R. Braun,
J. N. H. S. Aditya,
S. Aicardi,
L. Alegre,
A. Chakraborty,
X. Chen,
S. Choudhuri,
A. O. Clarke,
J. Coles,
J. S. Collinson,
D. Cornu,
L. Darriba,
M. Delli Veneri,
J. Forbrich,
B. Fraga,
A. Galan,
J. Garrido,
F. Gubanov,
H. Håkansson,
M. J. Hardcastle,
C. Heneka,
D. Herranz,
K. M. Hess
, et al. (83 additional authors not shown)
Abstract:
The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed t…
▽ More
The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets.
△ Less
Submitted 14 March, 2023;
originally announced March 2023.
-
Modeling lens potentials with continuous neural fields in galaxy-scale strong lenses
Authors:
Luca Biggio,
Georgios Vernardos,
Aymeric Galan,
Austin Peel
Abstract:
Strong gravitational lensing is a unique observational tool for studying the dark and luminous mass distribution both within and between galaxies. Given the presence of substructures, current strong lensing observations demand more complex mass models than smooth analytical profiles, such as power-law ellipsoids. In this work, we introduce a continuous neural field to predict the lensing potential…
▽ More
Strong gravitational lensing is a unique observational tool for studying the dark and luminous mass distribution both within and between galaxies. Given the presence of substructures, current strong lensing observations demand more complex mass models than smooth analytical profiles, such as power-law ellipsoids. In this work, we introduce a continuous neural field to predict the lensing potential at any position throughout the image plane, allowing for a nearly model-independent description of the lensing mass. We apply our method on simulated Hubble Space Telescope imaging data containing different types of perturbations to a smooth mass distribution: a localized dark subhalo, a population of subhalos, and an external shear perturbation. Assuming knowledge of the source surface brightness, we use the continuous neural field to model either the perturbations alone or the full lensing potential. In both cases, the resulting model is able to fit the imaging data, and we are able to accurately recover the properties of both the smooth potential and of the perturbations. Unlike many other deep learning methods, ours explicitly retains lensing physics (i.e., the lens equation) and introduces high flexibility in the model only where required, namely, in the lens potential. Moreover, the neural network does not require pre-training on large sets of labelled data and predicts the potential from the single observed lensing image. Our model is implemented in the fully differentiable lens modeling code Herculens.
△ Less
Submitted 17 October, 2022;
originally announced October 2022.
-
Using wavelets to capture deviations from smoothness in galaxy-scale strong lenses
Authors:
Aymeric Galan,
Georgios Vernardos,
Austin Peel,
Frédéric Courbin,
Jean-Luc Starck
Abstract:
Modeling the mass distribution of galaxy-scale strong gravitational lenses is a task of increasing difficulty. The high-resolution and depth of imaging data now available render simple analytical forms ineffective at capturing lens structures spanning a large range in spatial scale, mass scale, and morphology. In this work, we address the problem with a novel multiscale method based on wavelets. W…
▽ More
Modeling the mass distribution of galaxy-scale strong gravitational lenses is a task of increasing difficulty. The high-resolution and depth of imaging data now available render simple analytical forms ineffective at capturing lens structures spanning a large range in spatial scale, mass scale, and morphology. In this work, we address the problem with a novel multiscale method based on wavelets. We tested our method on simulated Hubble Space Telescope (HST) imaging data of strong lenses containing the following different types of mass substructures making them deviate from smooth models: (1) a localized small dark matter subhalo, (2) a Gaussian random field (GRF) that mimics a nonlocalized population of subhalos along the line of sight, and (3) galaxy-scale multipoles that break elliptical symmetry. We show that wavelets are able to recover all of these structures accurately. This is made technically possible by using gradient-informed optimization based on automatic differentiation over thousands of parameters, which also allow us to sample the posterior distributions of all model parameters simultaneously. By construction, our method merges the two main modeling paradigms - analytical and pixelated - with machine-learning optimization techniques into a single modular framework. It is also well-suited for the fast modeling of large samples of lenses. All methods presented here are publicly available in our new Herculens package.
△ Less
Submitted 21 November, 2022; v1 submitted 12 July, 2022;
originally announced July 2022.
-
Gravitationally lensed quasars in Gaia -- IV. 150 new lenses, quasar pairs, and projected quasars
Authors:
C. Lemon,
T. Anguita,
M. Auger,
F. Courbin,
A. Galan,
R. McMahon,
F. Neira,
M. Oguri,
P. Schechter,
A. Shajib,
T. Treu
Abstract:
We report the spectroscopic follow-up of 175 lensed quasar candidates selected using Gaia Data Release 2 observations following Lemon et al. 2019. Systems include 86 confirmed lensed quasars and a further 17 likely lensed quasars based on imaging and/or similar spectra. We also confirm 11 projected quasar pairs and 11 physical quasar pairs, while 25 systems are left as unclassified quasar pairs --…
▽ More
We report the spectroscopic follow-up of 175 lensed quasar candidates selected using Gaia Data Release 2 observations following Lemon et al. 2019. Systems include 86 confirmed lensed quasars and a further 17 likely lensed quasars based on imaging and/or similar spectra. We also confirm 11 projected quasar pairs and 11 physical quasar pairs, while 25 systems are left as unclassified quasar pairs -- pairs of quasars at the same redshift, which could be either distinct quasars or potential lensed quasars. Especially interesting objects include 8 quadruply imaged quasars of which two have BAL sources, an apparent triple, and a doubly lensed LoBaL quasar. The source redshifts and image separations of these new lenses range between 0.65 - 3.59 and 0.78 - 6.23 arcseconds respectively. We compare the known population of lensed quasars to an updated mock catalogue at image separations between 1 and 4 arcseconds, showing a very good match at z<1.5. At z>1.5, only 47% of the predicted number are known, with 56% of these missing lenses at image separations below 1.5 arcseconds. The missing higher-redshift, small-separation systems will have fainter lensing galaxies, and are partially explained by the unclassified quasar pairs and likely lenses presented in this work, which require deeper imaging. Of the 11 new reported projected quasar pairs, 5 have impact parameters below 10 kpc, almost tripling the number of such systems, which can probe the innermost regions of quasar host galaxies through absorption studies. We also report four new lensed galaxies discovered through our searches, with source redshifts ranging from 0.62 to 2.79.
△ Less
Submitted 15 June, 2022;
originally announced June 2022.
-
Lightweight HI source finding for next generation radio surveys
Authors:
Emma Tolley,
Damien Korber,
Aymeric Galan,
Austin Peel,
Mark T. Sargent,
Jean-Paul Kneib,
Frederic Courbin,
Jean-Luc Starck
Abstract:
Future deep HI surveys will be essential for understanding the nature of galaxies and the content of the Universe. However, the large volume of these data will require distributed and automated processing techniques. We introduce LiSA, a set of python modules for the denoising, detection and characterization of HI sources in 3D spectral data. LiSA was developed and tested on the Square Kilometer A…
▽ More
Future deep HI surveys will be essential for understanding the nature of galaxies and the content of the Universe. However, the large volume of these data will require distributed and automated processing techniques. We introduce LiSA, a set of python modules for the denoising, detection and characterization of HI sources in 3D spectral data. LiSA was developed and tested on the Square Kilometer Array Science Data Challenge 2 dataset, and contains modules and pipelines for easy domain decomposition and parallel execution. LiSA contains algorithms for 2D-1D wavelet denoising using the starlet transform and flexible source finding using null-hypothesis testing. These algorithms are lightweight and portable, needing only a few user-defined parameters reflecting the resolution of the data. LiSA also includes two convolutional neural networks developed to analyse data cubes which separate HI sources from artifacts and predict the HI source properties. All of these components are designed to be as modular as possible, allowing users to mix and match different components to create their ideal pipeline. We demonstrate the performance of the different components of LiSA on the SDC2 dataset, which is able to find 95% of HI sources with SNR > 3 and accurately predict their properties.
△ Less
Submitted 20 April, 2022;
originally announced April 2022.
-
TDCOSMO. IX. Systematic comparison between lens modelling software programs: time delay prediction for WGD 2038-4008
Authors:
A. J. Shajib,
K. C. Wong,
S. Birrer,
S. H. Suyu,
T. Treu,
E. Buckley-Geer,
H. Lin,
C. E. Rusu,
J. Poh,
A. Palmese,
A. Agnello,
M. W. Auger-Williams,
A. Galan,
S. Schuldt,
D. Sluse,
F. Courbin,
J. Frieman,
M. Millon
Abstract:
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 procedur…
▽ More
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 $\sim$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 $\sim$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.
△ Less
Submitted 21 September, 2022; v1 submitted 22 February, 2022;
originally announced February 2022.
-
TDCOSMO. VII. Boxyness/discyness in lensing galaxies : Detectability and impact on $H_0$
Authors:
Lyne Van de Vyvere,
Matthew R. Gomer,
Dominique Sluse,
Dandan Xu,
Simon Birrer,
Aymeric Galan,
Georgios Vernardos
Abstract:
In the context of gravitational lensing, the density profile of lensing galaxies is often considered to be perfectly elliptical. Potential angular structures are generally ignored, except to explain flux ratios anomalies. Surprisingly, the impact of azimuthal structures on extended images of the source has not been characterized, nor its impact on the H0 inference. We address this task by creating…
▽ More
In the context of gravitational lensing, the density profile of lensing galaxies is often considered to be perfectly elliptical. Potential angular structures are generally ignored, except to explain flux ratios anomalies. Surprisingly, the impact of azimuthal structures on extended images of the source has not been characterized, nor its impact on the H0 inference. We address this task by creating mock images of a point source embedded in an extended source, lensed by an elliptical galaxy on which multipolar components are added to emulate boxy/discy isodensity contours. Modeling such images with a density profile free of angular structure allow us to explore the detectability of image deformation induced by the multipoles in the residual frame. Multipole deformations are almost always detectable for our highest signal-to-noise mock data. However the detectability depends on the lens ellipticity and Einstein radius, on the S/N of the data, and on the specific lens modeling strategy. Multipoles also introduce small changes to the time delays. We therefore quantify how undetected multipoles would impact H0 inference. When no multipoles are detected in the residuals, the impact on H0 for a given lens is in general less than a few km/s/Mpc, but in the worst case scenario, combining low S/N in the ring and large intrinsic boxyness/discyness, the bias on H0 can reach 10-12 km/s/Mpc. If we now look at the inference on H0 from a population of lensing galaxies, having a distribution of multipoles representative of what is found in the light-profile of elliptical galaxies, we then find a systematic bias on H0 < 1%. The comparison of our mock systems to the state-of-the-art time delay lens sample studied by the H0LiCOW and TDCOSMO collaborations, indicates that multipoles are currently unlikely to be a source of substantial systematic bias on the inferred value of H0 from time-delay lenses.
△ Less
Submitted 3 March, 2022; v1 submitted 7 December, 2021;
originally announced December 2021.
-
Strong lens systems search in the Dark Energy Survey using Convolutional Neural Networks
Authors:
K. Rojas,
E. Savary,
B. Clément,
M. Maus,
F. Courbin,
C. Lemon,
J. H. H. Chan,
G. Vernardos,
R. Joseph,
R. Cañameras,
A. Galan
Abstract:
We performed a search for strong lens galaxy-scale systems in the first data release of the Dark Energy Survey (DES), from a color-selected parent sample of 18~745~029 Luminous Red Galaxies (LRGs). Our search was based on a Convolutional Neural Network (CNN) to grade our LRG selection with values between 0 (non-lens) and 1 (lens). Our training set was data-driven, i.e. using lensed sources taken f…
▽ More
We performed a search for strong lens galaxy-scale systems in the first data release of the Dark Energy Survey (DES), from a color-selected parent sample of 18~745~029 Luminous Red Galaxies (LRGs). Our search was based on a Convolutional Neural Network (CNN) to grade our LRG selection with values between 0 (non-lens) and 1 (lens). Our training set was data-driven, i.e. using lensed sources taken from HST COSMOS images and where the light distribution of the lens plane was taken directly from DES images of our LRGs. A total of 76~582 cutouts obtained a score above 0.9. These were visually inspected and resulted in two catalogs. The first one contains 405 lens candidates, where 90 present clear lensing features and counterparts, while the others 315 require more evidence, such as higher resolution images or spectra to be conclusive. A total of 186 candidates were totally new identified in this search. The second catalog includes 539 ring galaxy candidates that will be useful to train CNNs against this type of false positives. For the 90 best lens candidates we carried out color-based deblending of the lens and source light without fitting any analytical profile to the data. The method turned out to be very efficient in the deblending, even for very compact objects and for objects with very complex morphology. Finally, from the 90 best lens candidates we selected 52 systems having one single deflector, to test an automated modeling pipeline which successfully modeled 79\% of the sample within an acceptable amount of computing time.
△ Less
Submitted 31 August, 2021;
originally announced September 2021.
-
lenstronomy II: A gravitational lensing software ecosystem
Authors:
Simon Birrer,
Anowar J. Shajib,
Daniel Gilman,
Aymeric Galan,
Jelle Aalbers,
Martin Millon,
Robert Morgan,
Giulia Pagano,
Ji Won Park,
Luca Teodori,
Nicolas Tessore,
Madison Ueland,
Lyne Van de Vyvere,
Sebastian Wagner-Carena,
Ewoud Wempe,
Lilan Yang,
Xuheng Ding,
Thomas Schmidt,
Dominique Sluse,
Ming Zhang,
Adam Amara
Abstract:
lenstronomy is an Astropy-affiliated Python package for gravitational lensing simulations and analyses. lenstronomy was introduced by Birrer and Amara (2018) and is based on the linear basis set approach by Birrer et a. (2015). The user and developer base of lenstronomy has substantially grown since then, and the software has become an integral part of a wide range of recent analyses, such as meas…
▽ More
lenstronomy is an Astropy-affiliated Python package for gravitational lensing simulations and analyses. lenstronomy was introduced by Birrer and Amara (2018) and is based on the linear basis set approach by Birrer et a. (2015). The user and developer base of lenstronomy has substantially grown since then, and the software has become an integral part of a wide range of recent analyses, such as measuring the Hubble constant with time-delay strong lensing or constraining the nature of dark matter from resolved and unresolved small scale lensing distortion statistics. The modular design has allowed the community to incorporate innovative new methods, as well as to develop enhanced software and wrappers with more specific aims on top of the lenstronomy API. Through community engagement and involvement, lenstronomy has become a foundation of an ecosystem of affiliated packages extending the original scope of the software and proving its robustness and applicability at the forefront of the strong gravitational lensing community in an open source and reproducible manner.
△ Less
Submitted 10 June, 2021;
originally announced June 2021.
-
SLITronomy: towards a fully wavelet-based strong lensing inversion technique
Authors:
A. Galan,
A. Peel,
R. Joseph,
F. Courbin,
J. -L. Starck
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 a…
▽ More
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.
△ Less
Submitted 4 December, 2020;
originally announced December 2020.
-
TDCOSMO IV: Hierarchical time-delay cosmography -- joint inference of the Hubble constant and galaxy density profiles
Authors:
S. Birrer,
A. J. Shajib,
A. Galan,
M. Millon,
T. Treu,
A. Agnello,
M. Auger,
G. C. -F. Chen,
L. Christensen,
T. Collett,
F. Courbin,
C. D. Fassnacht,
L. V. E. Koopmans,
P. J. Marshall,
J. -W. Park,
C. E. Rusu,
D. Sluse,
C. Spiniello,
S. H. Suyu,
S. Wagner-Carena,
K. C. Wong,
M. Barnabè,
A. S. Bolton,
O. Czoske,
X. Ding
, et al. (2 additional authors not shown)
Abstract:
The H0LiCOW collaboration inferred via gravitational lensing time delays a Hubble constant $H_0=73.3^{+1.7}_{-1.8}$ km s$^{-1}{\rm Mpc}^{-1}$, describing deflector mass density profiles by either a power-law or stars 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 $H_0$. We qu…
▽ More
The H0LiCOW collaboration inferred via gravitational lensing time delays a Hubble constant $H_0=73.3^{+1.7}_{-1.8}$ km s$^{-1}{\rm Mpc}^{-1}$, describing deflector mass density profiles by either a power-law or stars 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 $H_0$. We quantify any potential effect of the MST with flexible mass models that are maximally degenerate with H0. Our calculation is based on a new hierarchical approach in which the MST is only constrained by stellar kinematics. The approach is validated on hydrodynamically simulated lenses. We apply the method to the TDCOSMO sample of 7 lenses (6 from H0LiCOW) and measure $H_0=74.5^{+5.6}_{-6.1}$ km s$^{-1}{\rm Mpc}^{-1}$. In order to further constrain the deflector mass profiles, we then add imaging and spectroscopy for 33 strong gravitational lenses from the SLACS sample. For 9 of the SLAC lenses we use resolved kinematics to constrain the stellar anisotropy. From the joint analysis of the TDCOSMO+SLACS sample, we measure $H_0=67.4^{+4.1}_{-3.2}$ km s$^{-1}{\rm Mpc}^{-1}$, assuming 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. While our new analysis does not statistically invalidate the mass profile assumptions by H0LiCOW, and thus their $H_0$ measurement relying on those, it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on $H_0$ derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data, chiefly spatially resolved kinematics of lens galaxies.
△ Less
Submitted 19 December, 2020; v1 submitted 6 July, 2020;
originally announced July 2020.
-
Time Delay Lens Modelling Challenge
Authors:
X. Ding,
T. Treu,
S. Birrer,
G. C. -F. Chen,
J. Coles,
P. Denzel,
M. Frigo A. Galan,
P. J. Marshall,
M. Millon,
A. More,
A. J. Shajib,
D. Sluse,
H. Tak,
D. Xu,
M. W. Auger,
V. Bonvin,
H. Chand,
F. Courbin,
G. Despali,
C. D. Fassnacht,
D. Gilman,
S. Hilbert,
S. R. Kumar,
Y. -Y. Lin,
J. W. Park
, et al. (4 additional authors not shown)
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 $H_0$. However, published state-of-the-art analyses require of order 1 year 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…
▽ More
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 $H_0$. However, published state-of-the-art analyses require of order 1 year 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 datasets. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20\%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ($ |A| < 2\%$) and precision ($< 6\%$ per system), even in the presence of 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.
△ Less
Submitted 17 February, 2021; v1 submitted 15 June, 2020;
originally announced June 2020.
-
Testing the Evolution of the Correlations between Supermassive Black Holes and their Host Galaxies using Eight Strongly Lensed Quasars
Authors:
Xuheng Ding,
Tommaso Treu,
Simon Birrer,
Adriano Agnello,
Dominique Sluse,
Chris Fassnacht,
Matthew W. Auger,
Kenneth C. Wong,
Sherry H. Suyu,
Takahiro Morishita,
Cristian E. Rusu,
Aymeric Galan
Abstract:
One of the main challenges in using high redshift active galactic nuclei to study the correlations between the mass of the supermassive Black Hole (MBH) and the properties of their active host galaxies is instrumental resolution. Strong lensing magnification effectively increases instrumental resolution and thus helps to address this challenge. In this work, we study eight strongly lensed active g…
▽ More
One of the main challenges in using high redshift active galactic nuclei to study the correlations between the mass of the supermassive Black Hole (MBH) and the properties of their active host galaxies is instrumental resolution. Strong lensing magnification effectively increases instrumental resolution and thus helps to address this challenge. In this work, we study eight strongly lensed active galactic nuclei (AGN) with deep Hubble Space Telescope imaging, using the lens modelling code Lenstronomy to reconstruct the image of the source. Using the reconstructed brightness of the host galaxy, we infer the host galaxy stellar mass based on stellar population models. MBH are estimated from broad emission lines using standard methods. Our results are in good agreement with recent work based on non-lensed AGN, demonstrating the potential of using strongly lensed AGNs to extend the study of the correlations to higher redshifts. At the moment, the sample size of lensed AGN is small and thus they provide mostly a consistency check on systematic errors related to resolution for the non-lensed AGN. However, the number of known lensed AGN is expected to increase dramatically in the next few years, through dedicated searches in ground and space based wide field surveys, and they may become a key diagnostic of black hole and galaxy co-evolution.
△ Less
Submitted 13 November, 2020; v1 submitted 27 May, 2020;
originally announced May 2020.
-
COSMOGRAIL XIX: Time delays in 18 strongly lensed quasars from 15 years of optical monitoring
Authors:
M. Millon,
F. Courbin,
V. Bonvin,
E. Paic,
G. Meylan,
M. Tewes,
D. Sluse,
P. Magain,
J. H. H. Chan,
A. Galan,
R. Joseph,
C. Lemon,
O. Tihhonova,
R. I. Anderson,
M. Marmier,
B. Chazelas,
M. Lendl,
A. H. M. J. Triaud,
A. Wyttenbach
Abstract:
We present the results of 15 years of monitoring lensed quasars, which was conducted by the COSMOGRAIL programme at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed systems are presented for the first time. We complement our data set with other monitoring data available in the literature to measure the time delays in 18 systems, among which nine reach a relative p…
▽ More
We present the results of 15 years of monitoring lensed quasars, which was conducted by the COSMOGRAIL programme at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed systems are presented for the first time. We complement our data set with other monitoring data available in the literature to measure the time delays in 18 systems, among which nine reach a relative precision better than 15% for at least one time delay. To achieve this, we developed an automated version of the curve-shifting toolbox PyCS to ensure robust estimation of the time delay in the presence of microlensing, while accounting for the errors due to the imperfect representation of microlensing. We also re-analysed the previously published time delays of RX J1131$-$1231 and HE 0435$-$1223, by adding six and two new seasons of monitoring, respectively, and confirming the previous time-delay measurements. When the time delay measurement is possible, we corrected the light curves of the lensed images from their time delay and present the difference curves to highlight the microlensing signal contained in the data. To date, this is the largest sample of decade-long lens monitoring data, which is useful to measure $H_0$ and the size of quasar accretion discs with microlensing as well as to study quasar variability.
△ Less
Submitted 25 August, 2020; v1 submitted 13 February, 2020;
originally announced February 2020.
-
TDCOSMO. I. An exploration of systematic uncertainties in the inference of $H_0$ from time-delay cosmography
Authors:
M. Millon,
A. Galan,
F. Courbin,
T. Treu,
S. H. Suyu,
X. Ding,
S. Birrer,
G. C. -F. Chen,
A. J. Shajib,
D. Sluse,
K. C. Wong,
A. Agnello,
M. W. Auger,
E. J. Buckley-Geer,
J. H. H. Chan,
T. Collett,
C. D. Fassnacht,
S. Hilbert,
L. V. E. Koopmans,
V. Motta,
S. Mukherjee,
C. E. Rusu,
A. Sonnenfeld,
C. Spiniello,
L. Van de Vyvere
Abstract:
Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES…
▽ More
Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.
△ Less
Submitted 22 July, 2020; v1 submitted 17 December, 2019;
originally announced December 2019.
-
STRIDES: a 3.9 per cent measurement of the Hubble constant from the strong lens system DES J0408-5354
Authors:
A. J. Shajib,
S. Birrer,
T. Treu,
A. Agnello,
E. J. Buckley-Geer,
J. H. H. Chan,
L. Christensen,
C. Lemon,
H. Lin,
M. Millon,
J. Poh,
C. E. Rusu,
D. Sluse,
C. Spiniello,
G. C. -F. Chen,
T. Collett,
F. Courbin,
C. D. Fassnacht,
J. Frieman,
A. Galan,
D. Gilman,
A. More,
T. Anguita,
M. W. Auger,
V. Bonvin
, et al. (66 additional authors not shown)
Abstract:
We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for th…
▽ More
We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for this lens system using three band Hubble Space Telescope imaging. We combine the measured time delays, line-of-sight central velocity dispersion of the deflector, and statistically constrained external convergence with our lens models to estimate two cosmological distances. We measure the "effective" time-delay distance corresponding to the redshifts of the deflector and the lensed quasar $D_{Δt}^{\rm eff}=3382^{+146}_{-115}$ Mpc and the angular diameter distance to the deflector $D_{\rm d}=1711^{+376}_{-280}$ Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant $H_0 = 74.2^{+2.7}_{-3.0}$ km s$^{-1}$ Mpc$^{-1}$ assuming a flat $Λ$CDM cosmology and a uniform prior for $Ω_{\rm m}$ as $Ω_{\rm m} \sim \mathcal{U}(0.05, 0.5)$. This measurement gives the most precise constraint on $H_0$ to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analyzed by the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of $H_0$ based on the local distance ladder, reinforcing the tension with the inference from early Universe probes, for example, with 2.2$σ$ discrepancy from the cosmic microwave background measurement.
△ Less
Submitted 14 April, 2020; v1 submitted 14 October, 2019;
originally announced October 2019.