-
General Multipoles and Their Implications for Dark Matter Inference
Authors:
Jacob S. Cohen,
Christopher D. Fassnacht,
Conor M. O'Riordan,
Simona Vegetti
Abstract:
The flux ratios of strongly lensed quasars have previously been used to infer the properties of dark matter. In these analyses it is crucial to separate the effect of the main lensing galaxy and the low-mass dark matter halo population. In this work, we investigate flux-ratio perturbations resulting from general third- and fourth-order multipole perturbations to the main lensing galaxy's mass prof…
▽ More
The flux ratios of strongly lensed quasars have previously been used to infer the properties of dark matter. In these analyses it is crucial to separate the effect of the main lensing galaxy and the low-mass dark matter halo population. In this work, we investigate flux-ratio perturbations resulting from general third- and fourth-order multipole perturbations to the main lensing galaxy's mass profile. We simulate four lens systems, each with a different lensing configuration, without multipoles. The simulated flux ratios are perturbed by 10-40 per cent by a population of low-mass haloes consistent with CDM and, in one case, also a satellite galaxy. This level of perturbation is comparable to the magnitude of flux-ratio anomalies in real data that has been previously analyzed. We then attempt to fit the simulated systems using multipoles instead of low-mass haloes. We find that multipoles with amplitudes of 0.01 or less can produce flux-ratio perturbations in excess of 40 per cent. In all cases, third- or fourth-order multipoles can individually reduce the magnitude of, if not eliminate, flux-ratio anomalies. When both multipole orders are jointly included, all simulated flux ratios can be fit to within the observational uncertainty. Our results indicate that low-mass haloes and multipoles are highly degenerate when modelling quadruply-imaged quasars based just on image positions and flux ratios. In the presence of this degeneracy, flux-ratio anomalies in lensed quasars alone cannot be used to place strong constraints on the properties of dark matter without additional information that can inform our priors.
△ Less
Submitted 15 March, 2024; v1 submitted 13 March, 2024;
originally announced March 2024.
-
Complex angular structure of three elliptical galaxies from high-resolution ALMA observations of strong gravitational lenses
Authors:
H. R. Stacey,
D. M. Powell,
S. Vegetti,
J. P. McKean,
C. D. Fassnacht,
D. Wen,
C. M. O'Riordan
Abstract:
The large-scale mass distributions of galaxy-scale strong lenses have long been assumed to be well-described by a singular ellipsoidal power-law density profile with external shear. However, the inflexibility of this model could lead to systematic errors in astrophysical parameters inferred with gravitational lensing observables. Here, we present observations with the Atacama Large (sub-)Millimetr…
▽ More
The large-scale mass distributions of galaxy-scale strong lenses have long been assumed to be well-described by a singular ellipsoidal power-law density profile with external shear. However, the inflexibility of this model could lead to systematic errors in astrophysical parameters inferred with gravitational lensing observables. Here, we present observations with the Atacama Large (sub-)Millimetre Array (ALMA) of three strongly lensed dusty star-forming galaxies at $\simeq30$ mas angular resolution and investigate the sensitivity of these data to angular structure in the lensing galaxies. We jointly infer the lensing mass distribution and the full surface brightness of the lensed sources with multipole expansions of the power-law density profile up to fourth order using a technique developed for interferometric data. All three data sets strongly favour third and fourth-order multipole amplitudes of $\approx1$ percent of the convergence. While the infrared stellar isophotes and isodensity shapes agree for one lens system, for the other two the isophotes disagree to varying extents, suggesting contributions to the angular structure from dark matter intrinsic or extrinsic to the lensing galaxy.
△ Less
Submitted 27 May, 2024; v1 submitted 7 March, 2024;
originally announced March 2024.
-
Angular complexity in strong lens substructure detection
Authors:
Conor M. O'Riordan,
Simona Vegetti
Abstract:
Strong gravitational lensing can be used to find otherwise invisible dark matter subhaloes. In such an analysis, the lens galaxy mass model is a significant source of systematic uncertainty. In this paper we analyse the effect of angular complexity in the lens model. We use multipole perturbations which introduce low-order deviations from pure ellipticity in the isodensity contours, keeping the ra…
▽ More
Strong gravitational lensing can be used to find otherwise invisible dark matter subhaloes. In such an analysis, the lens galaxy mass model is a significant source of systematic uncertainty. In this paper we analyse the effect of angular complexity in the lens model. We use multipole perturbations which introduce low-order deviations from pure ellipticity in the isodensity contours, keeping the radial density profile fixed. We find that, in HST-like data, multipole perturbations consistent with those seen in galaxy isophotes are very effective at causing false positive substructure detections. We show that the effectiveness of this degeneracy depends on the deviation from a pure ellipse and the lensing configuration. We find that, when multipoles of one per cent are allowed in the lens model, the area in the observation where a subhalo could be detected drops by a factor of three. Sensitivity away from the lensed images is mostly lost. However, the mass limit of detectable objects on or close to the lensed images does not change. We do not expect the addition of multipole perturbations to lens models to have a significant effect on the ability of strong lensing to constrain the underlying dark matter model. However, given the high rate of false positive detections, angular complexity beyond the elliptical power-law should be included for such studies to be reliable. We discuss implications for previous detections and future work.
△ Less
Submitted 12 January, 2024; v1 submitted 16 October, 2023;
originally announced October 2023.
-
Euclid: Identification of asteroid streaks in simulated images using deep learning
Authors:
M. Pöntinen,
M. Granvik,
A. A. Nucita,
L. Conversi,
B. Altieri,
B. Carry,
C. M. O'Riordan,
D. Scott,
N. Aghanim,
A. Amara,
L. Amendola,
N. Auricchio,
M. Baldi,
D. Bonino,
E. Branchini,
M. Brescia,
S. Camera,
V. Capobianco,
C. Carbone,
J. Carretero,
M. Castellano,
S. Cavuoti,
A. Cimatti,
R. Cledassou,
G. Congedo
, et al. (92 additional authors not shown)
Abstract:
Up to 150000 asteroids will be visible in the images of the ESA Euclid space telescope, and the instruments of Euclid offer multiband visual to near-infrared photometry and slitless spectra of these objects. Most asteroids will appear as streaks in the images. Due to the large number of images and asteroids, automated detection methods are needed. A non-machine-learning approach based on the Strea…
▽ More
Up to 150000 asteroids will be visible in the images of the ESA Euclid space telescope, and the instruments of Euclid offer multiband visual to near-infrared photometry and slitless spectra of these objects. Most asteroids will appear as streaks in the images. Due to the large number of images and asteroids, automated detection methods are needed. A non-machine-learning approach based on the StreakDet software was previously tested, but the results were not optimal for short and/or faint streaks. We set out to improve the capability to detect asteroid streaks in Euclid images by using deep learning.
We built, trained, and tested a three-step machine-learning pipeline with simulated Euclid images. First, a convolutional neural network (CNN) detected streaks and their coordinates in full images, aiming to maximize the completeness (recall) of detections. Then, a recurrent neural network (RNN) merged snippets of long streaks detected in several parts by the CNN. Lastly, gradient-boosted trees (XGBoost) linked detected streaks between different Euclid exposures to reduce the number of false positives and improve the purity (precision) of the sample.
The deep-learning pipeline surpasses the completeness and reaches a similar level of purity of a non-machine-learning pipeline based on the StreakDet software. Additionally, the deep-learning pipeline can detect asteroids 0.25-0.5 magnitudes fainter than StreakDet. The deep-learning pipeline could result in a 50% increase in the number of detected asteroids compared to the StreakDet software. There is still scope for further refinement, particularly in improving the accuracy of streak coordinates and enhancing the completeness of the final stage of the pipeline, which involves linking detections across multiple exposures.
△ Less
Submitted 5 October, 2023;
originally announced October 2023.
-
Selection functions of strong lens finding neural networks
Authors:
A. Herle,
C. M. O'Riordan,
S. Vegetti
Abstract:
Convolution Neural Networks trained for the task of lens finding with similar architecture and training data as is commonly found in the literature are biased classifiers. An understanding of the selection function of lens finding neural networks will be key to fully realising the potential of the large samples of strong gravitational lens systems that will be found in upcoming wide-field surveys.…
▽ More
Convolution Neural Networks trained for the task of lens finding with similar architecture and training data as is commonly found in the literature are biased classifiers. An understanding of the selection function of lens finding neural networks will be key to fully realising the potential of the large samples of strong gravitational lens systems that will be found in upcoming wide-field surveys. We use three training datasets, representative of those used to train galaxy-galaxy and galaxy-quasar lens finding neural networks. The networks preferentially select systems with larger Einstein radii and larger sources with more concentrated source-light distributions. Increasing the detection significance threshold to 12$σ$ from 8$σ$ results in 50 per cent of the selected strong lens systems having Einstein radii $θ_\mathrm{E}$ $\ge$ 1.04 arcsec from $θ_\mathrm{E}$ $\ge$ 0.879 arcsec, source radii $R_S$ $\ge$ 0.194 arcsec from $R_S$ $\ge$ 0.178 arcsec and source Sérsic indices $n_{\mathrm{Sc}}^{\mathrm{S}}$ $\ge$ 2.62 from $n_{\mathrm{Sc}}^{\mathrm{S}}$ $\ge$ 2.55. The model trained to find lensed quasars shows a stronger preference for higher lens ellipticities than those trained to find lensed galaxies. The selection function is independent of the slope of the power-law of the mass profiles, hence measurements of this quantity will be unaffected. The lens finder selection function reinforces that of the lensing cross-section, and thus we expect our findings to be a general result for all galaxy-galaxy and galaxy-quasar lens finding neural networks.
△ Less
Submitted 19 July, 2023;
originally announced July 2023.
-
Strong gravitational lensing as a probe of dark matter
Authors:
S. Vegetti,
S. Birrer,
G. Despali,
C. D. Fassnacht,
D. Gilman,
Y. Hezaveh,
L. Perreault Levasseur,
J. P. McKean,
D. M. Powell,
C. M. O'Riordan,
G. Vernardos
Abstract:
Dark matter structures within strong gravitational lens galaxies and along their line of sight leave a gravitational imprint on the multiple images of lensed sources. Strong gravitational lensing provides, therefore, a key test of different dark matter models in a way that is independent of the baryonic content of matter structures on subgalactic scales. In this chapter, we describe how galaxy-sca…
▽ More
Dark matter structures within strong gravitational lens galaxies and along their line of sight leave a gravitational imprint on the multiple images of lensed sources. Strong gravitational lensing provides, therefore, a key test of different dark matter models in a way that is independent of the baryonic content of matter structures on subgalactic scales. In this chapter, we describe how galaxy-scale strong gravitational lensing observations are sensitive to the physical nature of dark matter. We provide a historical perspective of the field, and review its current status. We discuss the challenges and advances in terms of data, treatment of systematic errors and theoretical predictions, that will enable one to deliver a stringent and robust test of different dark matter models in the near future. With the advent of the next generation of sky surveys, the number of known strong gravitational lens systems is expected to increase by several orders of magnitude. Coupled with high-resolution follow-up observations, these data will provide a key opportunity to constrain the properties of dark matter with strong gravitational lensing.
△ Less
Submitted 20 June, 2023;
originally announced June 2023.
-
Abell 1201: Detection of an Ultramassive Black Hole in a Strong Gravitational Lens
Authors:
James. W. Nightingale,
Russell J. Smith,
Qiuhan He,
Conor M. O'Riordan,
Jacob A. Kegerreis,
Aristeidis Amvrosiadis,
Alastair C. Edge,
Amy Etherington,
Richard G. Hayes,
Ash Kelly,
John R. Lucey,
Richard J. Massey Richard J. Massey
Abstract:
Supermassive black holes (SMBHs) are a key catalyst of galaxy formation and evolution, leading to an observed correlation between SMBH mass $M_{\rm BH}$ and host galaxy velocity dispersion $σ_{\rm e}$. Outside the local Universe, measurements of $M_{\rm BH}$ are usually only possible for SMBHs in an active state: limiting sample size and introducing selection biases. Gravitational lensing makes it…
▽ More
Supermassive black holes (SMBHs) are a key catalyst of galaxy formation and evolution, leading to an observed correlation between SMBH mass $M_{\rm BH}$ and host galaxy velocity dispersion $σ_{\rm e}$. Outside the local Universe, measurements of $M_{\rm BH}$ are usually only possible for SMBHs in an active state: limiting sample size and introducing selection biases. Gravitational lensing makes it possible to measure the mass of non-active SMBHs. We present models of the $z=0.169$ galaxy-scale strong lens Abell~1201. A cD galaxy in a galaxy cluster, it has sufficient `external shear' that a magnified image of a $z = 0.451$ background galaxy is projected just $\sim 1$ kpc from the galaxy centre. Using multi-band Hubble Space Telescope imaging and the lens modeling software $\texttt{PyAutoLens}$ we reconstruct the distribution of mass along this line of sight. Bayesian model comparison favours a point mass with $M_{\rm BH} = 3.27 \pm 2.12\times10^{10}\,$M$_{\rm \odot}$ (3$σ$ confidence limit); an ultramassive black hole. One model gives a comparable Bayesian evidence without a SMBH, however we argue this model is nonphysical given its base assumptions. This model still provides an upper limit of $M_{\rm BH} \leq 5.3 \times 10^{10}\,$M$_{\rm \odot}$, because a SMBH above this mass deforms the lensed image $\sim 1$ kpc from Abell 1201's centre. This builds on previous work using central images to place upper limits on $M_{\rm BH}$, but is the first to also place a lower limit and without a central image being observed. The success of this method suggests that surveys during the next decade could measure thousands more SMBH masses, and any redshift evolution of the $M_{\rm BH}$--$σ_{\rm e}$ relation. Results are available at https://github.com/Jammy2211/autolens_abell_1201.
△ Less
Submitted 27 March, 2023;
originally announced March 2023.
-
Sensitivity of strong lensing observations to dark matter substructure: a case study with Euclid
Authors:
Conor M. O'Riordan,
Giulia Despali,
Simona Vegetti,
Mark R. Lovell,
Ángeles Moliné
Abstract:
We introduce a machine learning method for estimating the sensitivity of strong lens observations to dark matter subhaloes in the lens. Our training data include elliptical power-law lenses, Hubble Deep Field sources, external shear, and noise and PSF for the Euclid VIS instrument. We set the concentration of the subhaloes using a $v_\mathrm{max}$-$r_\mathrm{max}$ relation. We then estimate the da…
▽ More
We introduce a machine learning method for estimating the sensitivity of strong lens observations to dark matter subhaloes in the lens. Our training data include elliptical power-law lenses, Hubble Deep Field sources, external shear, and noise and PSF for the Euclid VIS instrument. We set the concentration of the subhaloes using a $v_\mathrm{max}$-$r_\mathrm{max}$ relation. We then estimate the dark matter subhalo sensitivity in $16{,}000$ simulated strong lens observations with depth and resolution resembling Euclid VIS images. We find that, with a $3σ$ detection threshold, $2.35$ per cent of pixels inside twice the Einstein radius are sensitive to subhaloes with a mass $M_\mathrm{max}\leq 10^{10}M_\odot$, $0.03$ per cent are sensitive to $M_\mathrm{max}\leq 10^{9}M_\odot$, and, the limit of sensitivity is found to be $M_\mathrm{max}=10^{8.8\pm0.2}M_\odot$. Using our sensitivity maps and assuming CDM, we estimate that Euclid-like lenses will yield $1.43^{+0.14}_{-0.11}[f_\mathrm{sub}^{-1}]$ detectable subhaloes per lens in the entire sample, but this increases to $35.6^{+0.9}_{-0.9}[f_\mathrm{sub}^{-1}]$ per lens in the most sensitive lenses. Estimates are given in units of the inverse of the substructure mass fraction $f_\mathrm{sub}^{-1}$. Assuming $f_\mathrm{sub}=0.01$, one in every $70$ lenses in general should yield a detection, or one in every $\sim$ three lenses in the most sensitive sample. From $170,000$ new strong lenses detected by Euclid, we expect $\sim 2500$ new subhalo detections. We find that the expected number of detectable subhaloes in warm dark matter models only changes relative to cold dark matter for models which have already been ruled out, i.e., those with half-mode masses $M_\mathrm{hm}>10^8M_\odot$.
△ Less
Submitted 27 February, 2023; v1 submitted 28 November, 2022;
originally announced November 2022.
-
Realistic galaxy images and improved robustness in machine learning tasks from generative modelling
Authors:
Benjamin J. Holzschuh,
Conor M. O'Riordan,
Simona Vegetti,
Vicente Rodriguez-Gomez,
Nils Thuerey
Abstract:
We examine the capability of generative models to produce realistic galaxy images. We show that mixing generated data with the original data improves the robustness in downstream machine learning tasks. We focus on three different data sets; analytical Sérsic profiles, real galaxies from the COSMOS survey, and galaxy images produced with the SKIRT code, from the IllustrisTNG simulation. We quantif…
▽ More
We examine the capability of generative models to produce realistic galaxy images. We show that mixing generated data with the original data improves the robustness in downstream machine learning tasks. We focus on three different data sets; analytical Sérsic profiles, real galaxies from the COSMOS survey, and galaxy images produced with the SKIRT code, from the IllustrisTNG simulation. We quantify the performance of each generative model using the Wasserstein distance between the distributions of morphological properties (e.g. the Gini-coefficient, the asymmetry, and ellipticity), the surface brightness distribution on various scales (as encoded by the power-spectrum), the bulge statistic and the colour for the generated and source data sets. With an average Wasserstein distance (Fréchet Inception Distance) of $7.19 \times 10^{-2}\, (0.55)$, $5.98 \times 10^{-2}\, (1.45)$ and $5.08 \times 10^{-2}\, (7.76)$ for the Sérsic, COSMOS and SKIRT data set, respectively, our best models convincingly reproduce even the most complicated galaxy properties and create images that are visually indistinguishable from the source data. We demonstrate that by supplementing the training data set with generated data, it is possible to significantly improve the robustness against domain-shifts and out-of-distribution data. In particular, we train a convolutional neural network to denoise a data set of mock observations. By mixing generated images into the original training data, we obtain an improvement of $11$ and $45$ per cent in the model performance regarding domain-shifts in the physical pixel size and background noise level, respectively.
△ Less
Submitted 22 March, 2022;
originally announced March 2022.
-
The inner density profile of an elliptical galaxy at z=1.15 from gravitational lensing
Authors:
H. R. Stacey,
C. M. O'Riordan,
S. Vegetti,
D. M. Powell,
M. W. Auger,
G. Despali
Abstract:
The density profiles of lensing galaxies are typically parameterised by singular power-law models with a logarithmic slope close to isothermal ($ζ=2$). This is sufficient to fit the lensed emission near the Einstein radius but may not be sufficient when extrapolated to smaller or larger radii if the large-scale density profile is more complex. Here, we consider a broken power-law model for the den…
▽ More
The density profiles of lensing galaxies are typically parameterised by singular power-law models with a logarithmic slope close to isothermal ($ζ=2$). This is sufficient to fit the lensed emission near the Einstein radius but may not be sufficient when extrapolated to smaller or larger radii if the large-scale density profile is more complex. Here, we consider a broken power-law model for the density profile of an elliptical galaxy at $z=1.15$ using observations with the Atacama Large (sub-)Millimetre Array of the strong gravitational lens system SPT0532$-$50. This is the first application of such a model to real data. We find the lensed emission is best fit by a density profile that is sub-isothermal ($ζ= 1.87^{+0.02}_{-0.03}$) near the Einstein radius and steepens to super-isothermal ($ζ= 2.14^{+0.03}_{-0.02}$) at around half the Einstein radius, demonstrating that the lensing data probes the mass distribution inside the region probed by the lensed images. Assuming that a broken power-law is the underlying truth, we find that a single power-law would result in a $10\pm1$ percent underestimate of the Hubble constant from time-delay cosmography. Our results suggest that a broken power-law could be useful for precision lens modelling and probing the structural evolution of elliptical galaxies.
△ Less
Submitted 12 October, 2021; v1 submitted 22 September, 2021;
originally announced September 2021.
-
Galaxy mass profiles from strong lensing III: The two-dimensional broken power-law model
Authors:
Conor M. O'Riordan,
Stephen J. Warren,
Daniel J. Mortlock
Abstract:
When modelling strong gravitational lenses, i.e., where there are multiple images of the same source, the most widely used parameterisation for the mass profile in the lens galaxy is the singular power-law model $ρ(r)\propto r^{-γ}$. This model may be insufficiently flexible for very accurate work, for example measuring the Hubble constant based on time delays between multiple images. Here we deri…
▽ More
When modelling strong gravitational lenses, i.e., where there are multiple images of the same source, the most widely used parameterisation for the mass profile in the lens galaxy is the singular power-law model $ρ(r)\propto r^{-γ}$. This model may be insufficiently flexible for very accurate work, for example measuring the Hubble constant based on time delays between multiple images. Here we derive the lensing properties - deflection angle, shear, and magnification - of a more adaptable model where the projected mass surface density is parameterised as a continuous two-dimensional broken power-law (2DBPL). This elliptical 2DBPL model is characterised by power-law slopes $t_1$, $t_2$ either side of the break radius $θ_\mathrm{B}$. The key to the 2DBPL model is the derivation of the lensing properties of the truncated power law (TPL) model, where the surface density is a power law out to the truncation radius $θ_\mathrm{T}$ and zero beyond. This TPL model is also useful by itself. We create mock observations of lensing by a TPL profile where the images form outside the truncation radius, so there is no mass in the annulus covered by the images. We then show that the slope of the profile interior to the images may be accurately recovered for lenses of moderate ellipticity. This demonstrates that the widely-held notion that lensing measures the slope of the mass profile in the annulus of the images, and is insensitive to the mass distribution at radii interior to the images, is incorrect.
△ Less
Submitted 2 December, 2020; v1 submitted 22 July, 2020;
originally announced July 2020.
-
Galaxy mass profiles from strong lensing II: The elliptical power-law model
Authors:
Conor M. O'Riordan,
Stephen J. Warren,
Daniel J. Mortlock
Abstract:
We present a systematic analysis of the constraints $σ_γ$ on the mass profile slope $γ$ obtainable when fitting a singular power-law ellipsoid model to a typical strong lensing observation of an extended source. These results extend our previous analysis of circular systems, Paper I. We draw our results from 676 mock observations covering a range of image configurations, each created with a fixed…
▽ More
We present a systematic analysis of the constraints $σ_γ$ on the mass profile slope $γ$ obtainable when fitting a singular power-law ellipsoid model to a typical strong lensing observation of an extended source. These results extend our previous analysis of circular systems, Paper I. We draw our results from 676 mock observations covering a range of image configurations, each created with a fixed signal to noise ratio $S=100$ in the images. We analyse the results using a combination of theory and a simplified model which identifies the contribution to the constraints of the individual fluxes and positions in each of the lensed images. The main results are: 1. For any lens ellipticity, the constraints $σ_γ$ for two image systems are well described by the results of Paper I, transformed to elliptical coordinates; 2. We derive an analytical expression for $σ_γ$ for systems with the source aligned with the axis of the lens; 3. For both two-image systems and aligned systems, $σ_γ$ is limited by the flux uncertainties; 4. The constraints for off-axis four-image systems are a factor of two to eight better, depending on source size, than for two-image systems, and improve with increasing lens ellipticity. We show that the constraints on $γ$ in these systems derive from the complementary positional information of the images alone, without using flux. The complementarity improves as the offset of the source from the axis increases, such that the best constraints $σ_γ<0.01$, for $S=100$, occur when the source approaches the caustic.
△ Less
Submitted 11 June, 2020; v1 submitted 28 April, 2020;
originally announced April 2020.
-
Galaxy mass profiles from strong lensing I: The circular power-law model
Authors:
Conor M. O'Riordan,
Stephen J. Warren,
Daniel J. Mortlock
Abstract:
In this series of papers we develop a formalism for constraining mass profiles in strong gravitational lenses with extended images, using fluxes in addition to positional information. We start in this paper with a circular power-law profile and show that the slope $γ$ is uniquely determined by only two observables: the flux ratio $f_1/f_2$ and the image position ratio $θ_1/θ_2$ of the two images.…
▽ More
In this series of papers we develop a formalism for constraining mass profiles in strong gravitational lenses with extended images, using fluxes in addition to positional information. We start in this paper with a circular power-law profile and show that the slope $γ$ is uniquely determined by only two observables: the flux ratio $f_1/f_2$ and the image position ratio $θ_1/θ_2$ of the two images. We derive an analytic expression relating these two observables to the slope, a result which does not depend on the Einstein angle or the structure or brightness of the source. We then find an expression for the uncertainty on the slope $σ_γ$ that depends only on the position ratio $θ_1/θ_2$ and the total S/N in the images. For example, in a system with position ratio $θ_1/θ_2=0.5$, S/N $=100$ and $γ=2$ we find that $γ$ is constrained to a precision of $\pm0.03$. We then test these results against a series of mock observations. We invert the images and fit an 11 parameter model, including ellipticity and position angle for both lens and source and measure the uncertainty on $γ$. We find agreement with the theoretical estimate for all mock observations. In future papers we will examine the radial range of the galaxy over which the constraint on the slope applies, and extend the analysis to elliptical lenses.
△ Less
Submitted 7 June, 2019; v1 submitted 17 April, 2019;
originally announced April 2019.