-
On the ellipticity parameterization for an NFW profile: an overlooked angular structure in strong lens modeling
Authors:
Matthew R. Gomer,
Dominique Sluse,
Lyne Van de Vyvere,
Simon Birrer,
Anowar J. Shajib,
Frederic Courbin
Abstract:
Galaxy-scale gravitational lenses are often modeled with two-component mass profiles where one component represents the stellar mass and the second is an NFW profile representing the dark matter. Outside of the spherical case, the NFW profile is costly to implement, and so it is approximated via two different methods; ellipticity can be introduced via the lensing potential (NFWp) or via the mass b…
▽ More
Galaxy-scale gravitational lenses are often modeled with two-component mass profiles where one component represents the stellar mass and the second is an NFW profile representing the dark matter. Outside of the spherical case, the NFW profile is costly to implement, and so it is approximated via two different methods; ellipticity can be introduced via the lensing potential (NFWp) or via the mass by approximating the NFW profile as a sum of analytical profiles (NFWm). While the NFWp method has been the default for lensing applications, it gives a different prescription of the azimuthal structure, which we show introduces ubiquitous gradients in ellipticity and boxiness in the mass distribution rather than having a constant elliptical shape. Because unmodeled azimuthal structure has been shown to be able to bias lens model results, we explore the degree to which this introduced azimuthal structure can affect the model accuracy. We construct input profiles using composite models using both the NFWp and NFWm methods and fit these mocks with a power-law elliptical mass distribution (PEMD) model with external shear. As a measure of the accuracy of the recovered lensing potential, we calculate the value of the Hubble parameter $H_0$ one would determine from the lensing fit. We find that the fits to the NFWp input return $H_0$ values which are systematically biased by about $3\%$ lower than the NFWm counterparts. We explore whether such an effect is attributable to the mass sheet transformation (MST) by using an MST-independent quantity, $ξ_2$. We show that, as expected, the NFWm mocks are degenerate with PEMD through an MST. For the NFWp, an additional bias is found beyond the MST due to azimuthal structures {\it exterior to the Einstein radius}. We recommend modelers use an NFWm prescription in the future, such that azimuthal structure can be introduced explicitly rather than implicitly.
△ Less
Submitted 4 October, 2023;
originally announced October 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.
-
TDCOSMO VIII: A key test of systematics in the hierarchical method of time-delay cosmography
Authors:
Matthew R. Gomer,
Dominique Sluse,
Lyne Van de Vyvere,
Simon Birrer,
Frederic Courbin
Abstract:
The largest source of systematic errors in the time-delay cosmography method likely arises from the lens model mass distribution, where an inaccurate choice of model could in principle bias the value of $H_0$. A Bayesian hierarchical framework has been proposed which combines lens systems with kinematic data, constraining the mass profile shape at a population level. The framework has been previou…
▽ More
The largest source of systematic errors in the time-delay cosmography method likely arises from the lens model mass distribution, where an inaccurate choice of model could in principle bias the value of $H_0$. A Bayesian hierarchical framework has been proposed which combines lens systems with kinematic data, constraining the mass profile shape at a population level. The framework has been previously validated on a small sample of lensing galaxies drawn from hydro-simulations. The goal of this work is to expand the validation to a more general set of lenses consistent with observed systems, as well as confirm the capacity of the method to combine two lens populations: one which has time delay information and one which lacks time delays and has systematically different image radii. For this purpose, we generate samples of analytic lens mass distributions made of baryons+dark matter and fit the subsequent mock images with standard power-law models. Corresponding kinematics data are also emulated. The hierarchical framework applied to an ensemble of time-delay lenses allows us to correct the $H_0$ bias associated with model choice, finding $H_0$ within $1.5σ$ of the fiducial value. We then combine this set with a sample of corresponding lens systems which have no time delays and have a source at lower $z$, resulting in a systematically smaller image radius relative to their effective radius. The hierarchical framework successfully accounts for this effect, recovering a value of $H_0$ which is both more precise ($σ\sim2\%$) and more accurate ($0.7\%$ median offset) than the time-delay set alone. This result confirms that non-time-delay lenses can nonetheless contribute valuable constraining power to the determination of $H_0$ via their kinematic constraints, assuming they come from the same global population as the time-delay set.
△ Less
Submitted 28 September, 2022; v1 submitted 5 September, 2022;
originally announced September 2022.
-
Consequences of the lack of azimuthal freedom in the modeling of lensing galaxies
Authors:
Lyne Van de Vyvere,
Dominique Sluse,
Matthew R. Gomer,
Sampath Mukherjee
Abstract:
Massive elliptical galaxies can display structures that deviate from a pure elliptical shape, such as a twist of the principal axis or variations in the axis ratio with galactocentric distance. Although satisfactory lens modeling is generally achieved without accounting for these azimuthal structures, the question about their impact on inferred lens parameters remains, in particular, on time delay…
▽ More
Massive elliptical galaxies can display structures that deviate from a pure elliptical shape, such as a twist of the principal axis or variations in the axis ratio with galactocentric distance. Although satisfactory lens modeling is generally achieved without accounting for these azimuthal structures, the question about their impact on inferred lens parameters remains, in particular, on time delays as they are used in time-delay cosmography. This paper aims at characterizing these effects and quantifying their impact considering realistic amplitudes of the variations. We achieved this goal by creating mock lensing galaxies with morphologies based on two data sets: observational data of local elliptical galaxies, and hydrodynamical simulations of elliptical galaxies at a typical lens redshift. We then simulated images of the lensing systems with space-based data quality and modeled them in a standard way to assess the impact of a lack of azimuthal freedom in the lens model. We find that twists in lensing galaxies are easily absorbed in homoeidal lens models by a change in orientation of the lens up to 10° with respect to the reference orientation at the Einstein radius, and of the shear by up to 20° with respect to the input shear orientation. The ellipticity gradients, on the other hand, can introduce a substantial amount of shear that may impact the radial mass model and consequently bias $H_0$, up to 10 km/s/Mpc. However, we find that light is a good tracer of azimuthal structures, meaning that direct imaging should be capable of diagnosing their presence. This in turn implies that such a large bias is unlikely to be unaccounted for in standard modeling practices. Furthermore, the overall impact of twists and ellipticity gradients averages out at a population level. For the galaxy populations we considered, the cosmological inference remains unbiased.
△ Less
Submitted 31 May, 2022;
originally announced June 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.
-
Galaxy-lens determination of $H_0$: the effect of the ellipse+shear modeling assumption
Authors:
Matthew R. Gomer,
Liliya L. R. Williams
Abstract:
Galaxy lenses are frequently modeled as an elliptical mass distribution with external shear and isothermal spheres to account for secondary and line-of-sight galaxies. There is statistical evidence that some fraction of observed quads are inconsistent with these assumptions, and require a dipole-like contribution to the mass with respect to the light. Simplifying assumptions about the shape of mas…
▽ More
Galaxy lenses are frequently modeled as an elliptical mass distribution with external shear and isothermal spheres to account for secondary and line-of-sight galaxies. There is statistical evidence that some fraction of observed quads are inconsistent with these assumptions, and require a dipole-like contribution to the mass with respect to the light. Simplifying assumptions about the shape of mass distributions can lead to the incorrect recovery of parameters such as $H_0$. We create several tests of synthetic quad populations with different deviations from an elliptical shape, then fit them with an ellipse+shear model, and measure the recovered values of $H_0$. Kinematic constraints are not included. We perform two types of fittings -- one with a single point source and one with an array of sources emulating an extended source. We carry out two model-free comparisons between our mock quads and the observed population. One result of these comparisons is a statistical inconsistency not yet mentioned in the literature: the image distance ratios with respect to the lens center of observed quads appear to span a much wider range than those of synthetic or simulated quads. Bearing this discrepancy in mind, our mock populations can result in biases on $H_0$ $\sim10\%$.
△ Less
Submitted 30 March, 2021;
originally announced March 2021.
-
Galaxy-lens determination of $H_0$: constraining density slope in the context of the mass sheet degeneracy
Authors:
Matthew R. Gomer,
Liliya L. R. Williams
Abstract:
Gravitational lensing offers a competitive method to measure $H_0$ with the goal of 1% precision. A major obstacle comes in the form of lensing degeneracies, such as the mass sheet degeneracy (MSD), which make it possible for a family of density profiles to reproduce the same lensing observables but return different values of $H_0$. The modeling process artificially selects one choice from this fa…
▽ More
Gravitational lensing offers a competitive method to measure $H_0$ with the goal of 1% precision. A major obstacle comes in the form of lensing degeneracies, such as the mass sheet degeneracy (MSD), which make it possible for a family of density profiles to reproduce the same lensing observables but return different values of $H_0$. The modeling process artificially selects one choice from this family, potentially biasing $H_0$. The effect is more pronounced when the profile of a given lens is not perfectly described by the lens model, which will always be the case to some extent. To explore this, we quantify the bias and spread in $H_0$ by creating quads from two-component mass models and fitting them with a power-law ellipse+shear model. We find that the bias does not correspond to the estimate one would calculate by transforming the profile into a power law near the image radius. We also emulate the effect of including stellar kinematics by performing fits where the slope is constrained to the true value. Informing the fit using the true value near the image radius can introduce substantial bias (0-23% depending on the model). We confirm using Jeans arguments that kinematic constraints can result in a biased value of $H_0$, though the degree of bias depends on the region kinematic modeling probes in specific lenses. We conclude that lensing degeneracies manifest through commonplace modeling approaches in a more complicated way than is assumed in the literature. If stellar kinematics incorrectly break the MSD, their inclusion may introduce more bias than their omission.
△ Less
Submitted 3 November, 2020; v1 submitted 19 July, 2019;
originally announced July 2019.
-
The impact of $Λ$CDM substructure and baryon-dark matter transition on the image positions of quad galaxy lenses
Authors:
Matthew R. Gomer,
Liliya L. R. Williams
Abstract:
The positions of multiple images in galaxy lenses are related to the galaxy mass distribution. Smooth elliptical mass profiles were previously shown to be inadequate in reproducing the quad population. In this paper, we explore the deviations from such smooth elliptical mass distributions. Unlike most other work, we use a model-free approach based on the relative polar image angles of quads, and t…
▽ More
The positions of multiple images in galaxy lenses are related to the galaxy mass distribution. Smooth elliptical mass profiles were previously shown to be inadequate in reproducing the quad population. In this paper, we explore the deviations from such smooth elliptical mass distributions. Unlike most other work, we use a model-free approach based on the relative polar image angles of quads, and their position in 3D space with respect to the Fundamental Surface of Quads. The FSQ is defined by quads produced by elliptical lenses. We have generated thousands of quads from synthetic populations of lenses with substructure consistent with $Λ$CDM simulations, and found that such perturbations are not sufficient to match the observed distribution of quads relative to the FSQ. The result is unchanged even when subhalo masses are increased by a factor of ten, and the most optimistic lensing selection bias is applied. We then produce quads from galaxies created using two components, representing baryons and dark matter. The transition from the mass being dominated by baryons in inner radii to being dominated by dark matter in outer radii can carry with it asymmetries, which would affect relative image angles. We run preliminary experiments using lenses with two elliptical mass components with nonidentical axis ratios and position angles, perturbations from ellipticity in the form of nonzero Fourier coefficients $a_4$ and $a_6$, and artificially offset ellipse centers as a proxy for asymmetry at image radii. We show that combination of these effects is a promising way of accounting for quad population properties. We conclude that the quad population provides a unique and sensitive tool for constraining detailed mass distribution in the centers of galaxies.
△ Less
Submitted 21 February, 2018; v1 submitted 20 December, 2017;
originally announced December 2017.