A new RTU grid method models the lensing source as a Gaussian process on a ray-transformed uniform grid, achieving comparable fits with roughly half the pixels per dimension and higher ELBOs on mock data.
JWST lensed quasar dark matter survey IV: Stringent warm dark matter constraints from the joint reconstruction of extended lensed arcs and quasar flux ratios
7 Pith papers cite this work. Polarity classification is still indexing.
abstract
We present a measurement of the free-streaming length of dark matter (DM) and subhalo abundance around 28 quadruple image strong lenses using observations from JWST MIRI presented in Paper III of this series. We improve on previous inferences on DM properties from lensed quasars by simultaneously reconstructing extended lensed arcs with image positions and relative magnifications (flux ratios). Our forward modeling framework generates full populations of subhalos, line-of-sight halos, and globular clusters, uses an accurate model for subhalo tidal evolution, and accounts for free-streaming effects on halo abundance and concentration. Modeling lensed arcs leads to more-precise model-predicted flux ratios, breaking covariance between subhalo abundance and the free-streaming scale parameterized by the half-mode mass $m_{\rm{hm}}$. Assuming subhalo abundance predicted by the semi-analytic model {\tt{galacticus}} ($N$-body simulations), we infer (Bayes factor of 10:1) $m_{\rm{hm}} < 10^{7.4} \mathrm{M}_{\odot}$ ($m_{\rm{hm}} < 10^{7.2} \mathrm{M}_{\odot}$), a 0.4 dex improvement relative to omitting lensed arcs. These bounds correspond to lower limits on thermal relic DM particle masses of $6.5$ and $7.4$ keV, respectively. Conversely, assuming DM is cold, we infer a projected mass in subhalos ($10^6 < m/M_{\odot}<10^{10.7}$) of $1.7_{-1.2}^{+2.6} \times 10^7 \ \mathrm{M}_{\odot} \ \rm{kpc^{-2}}$ at $95 \%$ confidence. This is consistent with {\tt{galacticus}} predictions ($0.9 \times 10^7 \mathrm{M}_{\odot} \ \rm{kpc^{-2}}$), but in mild tension with recent $N$-body simulations ($0.6 \times 10^7 \mathrm{M}_{\odot} \ \rm{kpc^{-2}}$). Our results are among the strongest bounds on WDM, and the most precise measurement of subhalo abundance around strong lenses. Further improvements will follow from the large sample of lenses to be discovered by Euclid, Rubin, and Roman.
years
2026 7representative citing papers
First narrow-line flux ratios for lensed quasar RXJ1131-1231 measured with JWST/NIRSpec IFS at ~5% precision, detecting cusp anomaly consistent with prior work.
A calibrated fitting function for the halo mass function that unifies predictions across CDM and non-CDM models over 10 orders of magnitude in mass with typical 12% precision after modeling systematics.
Strong lensing of 28 quasars constrains the CDM halo mass function cutoff to m_low < 10^{8.3} solar masses at 10:1 odds.
A subhalo of M200 = 2.78e10 solar masses and concentration 30 is detected at 5.8 sigma in PJ011646 via ALMA imaging and grid-based NFW search after fitting an elliptical power-law plus multipole macromodel.
The two-point correlation function of lensing deflection fields shows sensitivity to variations in SIDM subhalo core-collapse modeling at small scales.
Reports results from two searches for new radio lenses in existing surveys and discusses the completeness of the population usable for constraining dark matter properties via astrometric perturbations.
citing papers explorer
-
Gaussian processes on ray-guided transformed uniform grids for fast, flexible, and auto-differentiable adaptive source reconstruction in lens modelling
A new RTU grid method models the lensing source as a Gaussian process on a ray-transformed uniform grid, achieving comparable fits with roughly half the pixels per dimension and higher ELBOs on mock data.