-
Find the haystacks, then look for needles: The rate of strongly lensed transients in galaxy-galaxy strong gravitational lenses
Authors:
Ana Sainz de Murieta,
Thomas E. Collett,
Mark R. Magee,
Justin D. R. Pierel,
Wolfgang J. R. Enzi,
Martine Lokken,
Alex Gagliano,
Dan Ryczanowski
Abstract:
The time delay between appearances of multiple images of a gravitationally lensed supernova (glSN) is sensitive to the Hubble constant, $H_0$. As well as time delays, a lensed host galaxy is needed to enable precise inference of $H_0$. In this work we investigate the connection between discoverable lensed transients and their host galaxies. We find that LSST will discover 88 glSNe per year, of whi…
▽ More
The time delay between appearances of multiple images of a gravitationally lensed supernova (glSN) is sensitive to the Hubble constant, $H_0$. As well as time delays, a lensed host galaxy is needed to enable precise inference of $H_0$. In this work we investigate the connection between discoverable lensed transients and their host galaxies. We find that LSST will discover 88 glSNe per year, of which $54\%$ will also have a strongly lensed host. The rates can change by approximately 30 percent uncertainty depending primarily on the choice of unlensed SN population and uncertainties in the redshift evolution of the deflector population, but the fraction of glSNe with a lensed host is consistently around a half. LSST will discover 20 glSNe per year in systems that could plausibly have been identified by Euclid as galaxy-galaxy lenses before the discovery of the glSN. Such systems have preferentially longer time delays and therefore are well suited for cosmography. We define a golden sample of glSNe Ia with time delays over 10 days, image separations greater than 0.8 arcseconds, and a multiply imaged host. For this golden sample, we find $91\%$ occur in systems that should already be discoverable as galaxy-galaxy lenses in Euclid. For cosmology with glSNe, monitoring Euclid lenses is a plausible alternative to searching the entire LSST alert stream.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
How to Break the Mass Sheet Degeneracy with the Lightcurves of Microlensed Type Ia Supernovae
Authors:
Luke Weisenbach,
Thomas Collett,
Ana Sainz de Murieta,
Coleman Krawczyk,
Georgios Vernardos,
Wolfgang Enzi,
Andrew Lundgren
Abstract:
The standardizable nature of gravitationally lensed Type Ia supernovae (glSNe Ia) makes them an attractive target for time delay cosmography, since a source with known luminosity breaks the mass sheet degeneracy. It is known that microlensing by stars in the lensing galaxy can add significant stochastic uncertainty to the unlensed luminosity which is often much larger than the intrinsic scatter of…
▽ More
The standardizable nature of gravitationally lensed Type Ia supernovae (glSNe Ia) makes them an attractive target for time delay cosmography, since a source with known luminosity breaks the mass sheet degeneracy. It is known that microlensing by stars in the lensing galaxy can add significant stochastic uncertainty to the unlensed luminosity which is often much larger than the intrinsic scatter of the Ia population. In this work, we show how the temporal microlensing variations as the supernova disc expands can be used to improve the standardization of glSNe Ia. We find that SNe are standardizable if they do not cross caustics as they expand. We estimate that this will be the case for $\approx$6 doubly imaged systems and $\approx$0.3 quadruply imaged systems per year in LSST. At the end of the ten year LSST survey, these systems should enable us to test for systematics in $H_0$ due to the mass sheet degeneracy at the $1.00^{+0.07}_{-0.06}$\% level, or $1.8\pm0.2$\% if we can only extract time delays from the third of systems with counter images brighter than $i=24$ mag.
△ Less
Submitted 5 March, 2024;
originally announced March 2024.
-
Gravitational imaging through a triple source plane lens: revisiting the $Λ$CDM-defying dark subhalo in SDSSJ0946+1006
Authors:
Daniel J. Ballard,
Wolfgang J. R. Enzi,
Thomas E. Collett,
Hannah C. Turner,
Russell J. Smith
Abstract:
The $Λ$CDM paradigm successfully explains the large-scale structure of the Universe, but is less well constrained on sub-galactic scales. Gravitational lens modelling has been used to measure the imprints of dark substructures on lensed arcs, testing the small-scale predictions of $Λ$CDM. However, the methods required for these tests are subject to degeneracies among the lens mass model and the so…
▽ More
The $Λ$CDM paradigm successfully explains the large-scale structure of the Universe, but is less well constrained on sub-galactic scales. Gravitational lens modelling has been used to measure the imprints of dark substructures on lensed arcs, testing the small-scale predictions of $Λ$CDM. However, the methods required for these tests are subject to degeneracies among the lens mass model and the source light profile. We present a case study of the unique compound gravitational lens SDSSJ0946+1006, wherein a dark, massive substructure has been detected, whose reported high concentration would be unlikely in a $Λ$CDM universe. For the first time, we model the first two background sources in both I- and U-band HST imaging, as well as VLT-MUSE emission line data for the most distant source. We recover a lensing perturber at a $5.9σ$ confidence level with mass $\log_{10}(M_\mathrm{sub}/M_{\odot})=9.2^{+0.4}_{-0.1}$ and concentration $\log_{10}c=2.4^{+0.5}_{-0.3}$. The concentration is more consistent with CDM subhalos than previously reported, and the mass is compatible with that of a dwarf satellite galaxy whose flux is undetectable in the data at the location of the perturber. A wandering black hole with mass $\log_{10}(M_\mathrm{BH}/M_{\odot})=8.9^{+0.2}_{-0.1}$ is a viable alternative model. We systematically investigate alternative assumptions about the complexity of the mass distribution and source reconstruction; in all cases the subhalo is detected at around the $\geq5σ$ level. However, the detection significance can be altered substantially (up to $11.3σ$) by alternative choices for the source regularisation scheme.
△ Less
Submitted 28 February, 2024; v1 submitted 8 September, 2023;
originally announced September 2023.
-
Lensed type Ia supernovae in light of SN Zwicky and iPTF16geu
Authors:
Ana Sainz de Murieta,
Thomas E. Collett,
Mark R. Magee,
Luke Weisenbach,
Coleman M. Krawczyk,
Wolfgang Enzi
Abstract:
Strong gravitationally lensed supernovae (glSNe) are a powerful probe to obtain a measure of the expansion rate of the Universe, but they are also extremely rare. To date, only two glSNe with multiple images strongly lensed by galaxies have been found, but their short time delays make them unsuitable for cosmography. Here, we simulate a realistic catalogue of lensed supernovae and study the charac…
▽ More
Strong gravitationally lensed supernovae (glSNe) are a powerful probe to obtain a measure of the expansion rate of the Universe, but they are also extremely rare. To date, only two glSNe with multiple images strongly lensed by galaxies have been found, but their short time delays make them unsuitable for cosmography. Here, we simulate a realistic catalogue of lensed supernovae and study the characteristics of the population of detectable systems for different surveys. Compared to previous studies, our simulations also account for the effect of microlensing and its impact on the glSNe yields. We show that the properties of glSNe in shallow surveys (such as the Zwicky Transient Facility; ZTF) are determined by the need for large magnifications, which favours systems of four images with short time delays and low image separations. This picture is consistent with the properties of iPTF16geu and SN Zwicky, but is not representative of the population found in deeper surveys, which are limited by the volume of the Universe that is strongly lensed. For deeper surveys, such as the Legacy Survey of Space and Time (LSST), glSNe show longer time delays and greater angular separations, and the inclusion of microlensing results in 8$\%$ of glSNe becoming demagnified under the detection threshold. In the 10 years of the survey LSST should be able to find $\approx$ 180 systems, of which 70 will be suited for cosmography enabling a $\approx$ 1.2$\%$ precision $H_0$ measurement with LSST glSNe.
△ Less
Submitted 4 October, 2023; v1 submitted 24 July, 2023;
originally announced July 2023.
-
Cosmology from Large Populations of Galaxy-Galaxy Strong Gravitational Lenses
Authors:
Tian Li,
Thomas E. Collett,
Coleman M. Krawczyk,
Wolfgang Enzi
Abstract:
We present a forecast analysis on the feasibility of measuring the cosmological parameters with a large number of galaxy-galaxy scale strong gravitational lensing systems. Future wide area surveys are expected to discover and measure the properties of more than 10 000 strong lensing systems. We develop a hierarchical model that can simultaneously constrain the lens population and cosmological para…
▽ More
We present a forecast analysis on the feasibility of measuring the cosmological parameters with a large number of galaxy-galaxy scale strong gravitational lensing systems. Future wide area surveys are expected to discover and measure the properties of more than 10 000 strong lensing systems. We develop a hierarchical model that can simultaneously constrain the lens population and cosmological parameters by combining Einstein radius measurements with stellar dynamical mass estimates for every lens. Marginalizing over the lens density profiles and stellar orbital anisotropies, we find that $w$ can be constrained to a precision of $0.11$ with 10 000 galaxy-galaxy lens systems, which would be better than any existing single-probe constraint. We test our method on 161 existing lenses, finding $w=-0.96\pm0.46$. We also show how to mitigate against the potential systematic of redshift evolution in the mean lens density profile of the population.
△ Less
Submitted 18 July, 2023;
originally announced July 2023.
-
A search for gravitationally lensed supernovae within the Zwicky Transient Facility public survey
Authors:
M. R. Magee,
A. Sainz de Murieta,
T. E. Collett,
W. Enzi
Abstract:
Strong gravitational lensing of supernovae is exceedingly rare. To date, only a handful of lensed supernovae are known. Despite this, lensed supernovae have emerged as a promising method for measuring the current expansion rate of the Universe and breaking the Hubble tension. We present an extensive search for gravitationally lensed supernovae within the Zwicky Transient Facility (ZTF) public surv…
▽ More
Strong gravitational lensing of supernovae is exceedingly rare. To date, only a handful of lensed supernovae are known. Despite this, lensed supernovae have emerged as a promising method for measuring the current expansion rate of the Universe and breaking the Hubble tension. We present an extensive search for gravitationally lensed supernovae within the Zwicky Transient Facility (ZTF) public survey, covering 15,215 transients with good light curves discovered during four years of observations. We crossmatch a catalogue of known and candidate lens galaxies with our sample and find three coincident sources, which were due to chance alignment. To search for supernovae magnified by unknown lenses, we test multiple methods suggested in the literature for the first time on real data. This includes selecting objects with extremely red colours, those that appear inconsistent with the host galaxy redshift, and those with bright absolute magnitudes inferred from the host galaxy redshift. We find a few hundred candidates, most of which are due to contamination from activate galactic nuclei, bogus detections, or unlensed supernovae. The false positive rate from these methods presents significant challenges for future surveys. In total, 132 unique transients were identified across all of our selection methods that required detailed manual rejection, which would be infeasible for larger samples. Overall, we do not find any compelling candidates for lensed supernovae, which is broadly consistent with previous estimates for the rate of lensed supernovae from the ZTF public survey alone and the number expected to pass the selection cuts we apply.
△ Less
Submitted 31 May, 2023; v1 submitted 27 March, 2023;
originally announced March 2023.
-
Detecting low-mass haloes with strong gravitational lensing I: the effect of data quality and lensing configuration
Authors:
Giulia Despali,
Simona Vegetti,
Simon D. M. White,
Devon M. Powell,
Hannah R. Stacey,
Christopher D. Fassnacht,
Francesca Rizzo,
Wolfgang Enzi
Abstract:
This paper aims to quantify how the lowest halo mass that can be detected with galaxy-galaxy strong gravitational lensing depends on the quality of the observations and the characteristics of the observed lens systems. Using simulated data, we measure the lowest detectable NFW mass at each location of the lens plane, in the form of detailed \emph{sensitivity maps}. In summary, we find that: (i) th…
▽ More
This paper aims to quantify how the lowest halo mass that can be detected with galaxy-galaxy strong gravitational lensing depends on the quality of the observations and the characteristics of the observed lens systems. Using simulated data, we measure the lowest detectable NFW mass at each location of the lens plane, in the form of detailed \emph{sensitivity maps}. In summary, we find that: (i) the lowest detectable mass $M_{\rm low}$ decreases linearly as the signal-to-noise ratio (SNR) increases and the sensitive area is larger when we decrease the noise; (ii) a moderate increase in angular resolution (0.07" vs 0.09") and pixel scale (0.01" vs 0.04") improves the sensitivity by on average 0.25 dex in halo mass, with more significant improvement around the most sensitive regions; (iii) the sensitivity to low-mass objects is largest for bright and complex lensed galaxies located inside the caustic curves and lensed into larger Einstein rings (i.e $r_{E}\geq1.0"$). We find that for the sensitive mock images considered in this work, the minimum mass that we can detect at the redshift of the lens lies between $1.5\times10^{8}$ and $3\times10^{9}M_{\odot}$. We derive analytic relations between $M_{\rm low}$, the SNR and resolution and discuss the impact of the lensing configuration and source structure. Our results start to fill the gap between approximate predictions and real data and demonstrate the challenging nature of calculating precise forecasts for gravitational imaging. In light of our findings, we discuss possible strategies for designing strong lensing surveys and the prospects for HST, Keck, ALMA, Euclid and other future observations.
△ Less
Submitted 4 December, 2021; v1 submitted 16 November, 2021;
originally announced November 2021.
-
Joint constraints on thermal relic dark matter from strong gravitational lensing, the Lyman-$α$ forest, and Milky Way satellites
Authors:
Wolfgang Enzi,
Riccardo Murgia,
Oliver Newton,
Simona Vegetti,
Carlos Frenk,
Matteo Viel,
Marius Cautun,
Christopher D. Fassnacht,
Matt Auger,
Giulia Despali,
John McKean,
Léon V. E. Koopmans,
Mark Lovell
Abstract:
We derive joint constraints on the warm dark matter (WDM) half-mode scale by combining the analyses of a selection of astrophysical probes: strong gravitational lensing with extended sources, the Lyman-$α$ forest, and the number of luminous satellites in the Milky Way. We derive an upper limit of $λ_{\rm hm}=0.089{\rm~Mpc~h^{-1} }$ at the 95 per cent confidence level, which we show to be stable fo…
▽ More
We derive joint constraints on the warm dark matter (WDM) half-mode scale by combining the analyses of a selection of astrophysical probes: strong gravitational lensing with extended sources, the Lyman-$α$ forest, and the number of luminous satellites in the Milky Way. We derive an upper limit of $λ_{\rm hm}=0.089{\rm~Mpc~h^{-1} }$ at the 95 per cent confidence level, which we show to be stable for a broad range of prior choices. Assuming a Planck cosmology and that WDM particles are thermal relics, this corresponds to an upper limit on the half-mode mass of $M_{\rm hm }< 3 \times 10^{7} {\rm~M_{\odot}~h^{-1}}$, and a lower limit on the particle mass of $m_{\rm th }> 6.048 {\rm~keV}$, both at the 95 per cent confidence level. We find that models with $λ_{\rm hm}> 0.223 {\rm~Mpc~h^{-1} }$ (corresponding to $m_{\rm th }> 2.552 {\rm~keV}$ and $M_{\rm hm }< 4.8 \times 10^{8} {\rm~M_{\odot}~h^{-1}}$) are ruled out with respect to the maximum likelihood model by a factor $\leq 1/20$. For lepton asymmetries $L_6>10$, we rule out the $7.1 {\rm~keV}$ sterile neutrino dark matter model, which presents a possible explanation to the unidentified $3.55 {\rm~keV}$ line in the Milky Way and clusters of galaxies. The inferred 95 percentiles suggest that we further rule out the ETHOS-4 model of self-interacting DM. Our results highlight the importance of extending the current constraints to lower half-mode scales. We address important sources of systematic errors and provide prospects for how the constraints of these probes can be improved upon in the future.
△ Less
Submitted 22 June, 2021; v1 submitted 26 October, 2020;
originally announced October 2020.
-
Systematic errors in strong gravitational lensing reconstructions, a numerical simulation perspective
Authors:
W. Enzi,
S. Vegetti,
G. Despali,
J. -W. Hsueh,
R. B. Metcalf
Abstract:
We present the analysis of a sample of twenty-four SLACS-like galaxy-galaxy strong gravitational lens systems with a background source and deflectors from the Illustris-1 simulation. We study the degeneracy between the complex mass distribution of the lenses, substructures, the surface brightness distribution of the sources, and the time delays. Using a novel inference framework based on Approxima…
▽ More
We present the analysis of a sample of twenty-four SLACS-like galaxy-galaxy strong gravitational lens systems with a background source and deflectors from the Illustris-1 simulation. We study the degeneracy between the complex mass distribution of the lenses, substructures, the surface brightness distribution of the sources, and the time delays. Using a novel inference framework based on Approximate Bayesian Computation, we find that for all the considered lens systems, an elliptical and cored power-law mass density distribution provides a good fit to the data. However, the presence of cores in the simulated lenses affects most reconstructions in the form of a Source Position Transformation. The latter leads to a systematic underestimation of the source sizes by 50 per cent on average, and a fractional error in $H_{0}$ of around $25_{-19}^{+37}$ per cent. The analysis of a control sample of twenty-four lens systems, for which we have perfect knowledge about the shape of the lensing potential, leads to a fractional error on $H_{0}$ of $12_{-3}^{+6}$ per cent. We find no degeneracy between complexity in the lensing potential and the inferred amount of substructures. We recover an average total projected mass fraction in substructures of $f_{\rm sub}<1.7-2.0\times10^{-3}$ at the 68 per cent confidence level in agreement with zero and the fact that all substructures had been removed from the simulation. Our work highlights the need for higher-resolution simulations to quantify the lensing effect of more realistic galactic potentials better, and that additional observational constraint may be required to break existing degeneracies.
△ Less
Submitted 4 May, 2020; v1 submitted 6 November, 2019;
originally announced November 2019.
-
SHARP -- VII. New constraints on the dark matter free-streaming properties and substructure abundance from gravitationally lensed quasars
Authors:
Jen-Wei Hsueh,
Wolfgang Enzi,
Simona Vegetti,
Matt Auger,
Christopher D. Fassnacht,
Giulia Despali,
Leon V. E. Koopmans,
John P. McKean
Abstract:
We present an analysis of seven strongly gravitationally lensed quasars and the corresponding constraints on the properties of dark matter. Our results are derived by modelling the lensed image positions and flux-ratios using a combination of smooth macro models and a population of low-mass haloes within the mass range 10^6 to 10^9 Msun. Our lens models explicitly include higher-order complexity i…
▽ More
We present an analysis of seven strongly gravitationally lensed quasars and the corresponding constraints on the properties of dark matter. Our results are derived by modelling the lensed image positions and flux-ratios using a combination of smooth macro models and a population of low-mass haloes within the mass range 10^6 to 10^9 Msun. Our lens models explicitly include higher-order complexity in the form of stellar discs and luminous satellites, as well as low-mass haloes located along the observed lines of sight for the first time. Assuming a Cold Dark Matter (CDM) cosmology, we infer an average total mass fraction in substructure of f_sub = 0.012^{+0.007}_{-0.004} (68 per cent confidence limits), which is in agreement with the predictions from CDM hydrodynamical simulations to within 1 sigma. This result is closer to the predictions than those from previous studies that did not include line-of-sight haloes. Under the assumption of a thermal relic dark matter model, we derive a lower limit on the particle relic mass of m th > 5.58 keV (95 per cent confidence limits), which is consistent with a value of m_th > 5.3 keV from the recent analysis of the Ly-alpha forest. We also identify two main sources of possible systematic errors and conclude that deeper investigations in the complex structure of lens galaxies as well as the size of the background sources should be a priority for this field.
△ Less
Submitted 19 December, 2019; v1 submitted 10 May, 2019;
originally announced May 2019.
-
Primordial power spectrum and cosmology from black-box galaxy surveys
Authors:
Florent Leclercq,
Wolfgang Enzi,
Jens Jasche,
Alan Heavens
Abstract:
We propose a new, likelihood-free approach to inferring the primordial matter power spectrum and cosmological parameters from arbitrarily complex forward models of galaxy surveys where all relevant statistics can be determined from numerical simulations, i.e. black-boxes. Our approach, which we call simulator expansion for likelihood-free inference (SELFI), builds upon approximate Bayesian computa…
▽ More
We propose a new, likelihood-free approach to inferring the primordial matter power spectrum and cosmological parameters from arbitrarily complex forward models of galaxy surveys where all relevant statistics can be determined from numerical simulations, i.e. black-boxes. Our approach, which we call simulator expansion for likelihood-free inference (SELFI), builds upon approximate Bayesian computation using a novel effective likelihood, and upon the linearisation of black-box models around an expansion point. Consequently, we obtain simple "filter equations" for an effective posterior of the primordial power spectrum, and a straightforward scheme for cosmological parameter inference. We demonstrate that the workload is computationally tractable, fixed a priori, and perfectly parallel. As a proof of concept, we apply our framework to a realistic synthetic galaxy survey, with a data model accounting for physical structure formation and incomplete and noisy galaxy observations. In doing so, we show that the use of non-linear numerical models allows the galaxy power spectrum to be safely fitted up to at least $k_\mathrm{max} = 0.5$ $h$/Mpc, outperforming state-of-the-art backward-modelling techniques by a factor of $\sim 5$ in the number of modes used. The result is an unbiased inference of the primordial matter power spectrum across the entire range of scales considered, including a high-fidelity reconstruction of baryon acoustic oscillations. It translates into an unbiased and robust inference of cosmological parameters. Our results pave the path towards easy applications of likelihood-free simulation-based inference in cosmology. We have made our code pySELFI and our data products publicly available at http://pyselfi.florent-leclercq.eu.
△ Less
Submitted 8 October, 2019; v1 submitted 26 February, 2019;
originally announced February 2019.
-
Constraining sterile neutrino cosmologies with strong gravitational lensing observations at redshift z~0.2
Authors:
S. Vegetti,
G. Despali,
M. R. Lovell,
W. Enzi
Abstract:
We use the observed amount of subhaloes and line-of-sight dark matter haloes in a sample of 11 gravitational lens systems from the Sloan Lens ACS Survey to constrain the free-streaming properties of the dark matter particles. In particular, we combine the detection of a small-mass dark matter halo by Vegetti et al. 2010 with the non-detections by Vegetti et al. 2014 and compare the derived subhalo…
▽ More
We use the observed amount of subhaloes and line-of-sight dark matter haloes in a sample of 11 gravitational lens systems from the Sloan Lens ACS Survey to constrain the free-streaming properties of the dark matter particles. In particular, we combine the detection of a small-mass dark matter halo by Vegetti et al. 2010 with the non-detections by Vegetti et al. 2014 and compare the derived subhalo and halo mass functions with expectations from cold dark matter (CDM) and resonantly produced sterile neutrino models. We constrain the half-mode mass, i.e. the mass scale at which the linear matter power spectrum is reduced by 50 per cent relatively to the CDM model, to be $\log M_{\rm{hm}} \left[M_\odot\right] < 12.0$ (equivalent thermal relic mass $m_{\rm th} > 0.3$ keV) at the 2$σ$ level. This excludes sterile neutrino models with neutrino masses $m_{\rm s} < 0.8$ keV at any value of $L_{\rm 6}$. Our constraints are weaker than currently provided by the number of Milky Way satellites, observations of the 3.5 keV X-ray line, and the Lyman $α$ forest. However, they are more robust than the former as they are less affected by baryonic processes. Moreover, unlike the latter, they are not affected by assumptions on the thermal histories for the intergalactic medium. Gravitational lens systems with higher data quality and higher source and lens redshift are required to obtain tighter constraints.
△ Less
Submitted 9 October, 2018; v1 submitted 4 January, 2018;
originally announced January 2018.