REVIEW 2 major objections 2 minor 1 cited by
TDCOSMO XXVI: Uniform lens modeling of eight doubly imaged quasars
T0 review · 2 major / 2 minor · reviewed 2026-05-08 · grok-4.3
Pith's one-line read Surface brightness of lensed host arcs sets the precision of mass models in doubly imaged quasars.
desk verdict The paper shows arc brightness sets Fermat potential precision in doubles via a uniform Lenstronomy pipeline on eight systems, with solid literature checks and a conjugate-point test. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Correlation between Fermat potential uncertainty and host-arc surface brightness, measured through full extended-image modeling versus a conjugate-point analysis restricted to quasar positions.
What would settle it
Finding no correlation between measured arc surface brightness and Fermat-potential uncertainty when the same eight systems are re-modeled with an independent code or additional high-resolution data.
Extended reading notes
Core claim
Uniform modeling of the eight systems yields Einstein radii consistent with literature values at 1.5 sigma and image separations matching Gaia DR2 to 3.6 mas rms. Full image reconstruction and a conjugate-point analysis that uses only the quasar positions demonstrate that Fermat-potential precision improves directly with the surface brightness of the spatially extended host arcs. An anti-correlation between mass-parameter hypervolume and arc magnitude further isolates arc brightness as the primary driver of how well the lens mass profile can be recovered in doubly imaged systems.
Load-bearing premise
The tailored Lenstronomy pipeline accurately recovers the true lensing geometry and mass profiles of doubles without large biases from data heterogeneity or unmodeled systematics.
Editorial extensions
If this is right
- Doubly imaged systems can now be ranked for cosmographic usefulness by a directly observable quantity: arc surface brightness.
- The larger population of doubles can be incorporated into hierarchical H0 analyses once arc brightness is used to select or weight targets.
- Uniform pipelines become feasible for the thousands of new lenses expected from LSST, Roman, and Euclid.
- Model precision in doubles is shown to be limited by extended emission rather than by the mere number of point images.
Reading between the lines
- Survey strategies could prioritize follow-up of doubles that already show bright arcs in discovery imaging to maximize cosmographic return.
- The same brightness trend may allow statistical marginalization over lens-model uncertainty when building large H0 samples.
- Applying the conjugate-point test to triples or other configurations could reveal whether arc brightness remains the dominant constraint across image multiplicities.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents the first uniform gravitational lens modeling of eight doubly imaged quasars from HST multi-band data using a tailored open-source Lenstronomy pipeline. It reports 1.5σ average agreement on Einstein radii with literature values, 3.6 mas rms difference on image separations versus Gaia DR2, a strong correlation between Fermat potential precision and the surface brightness of the lensed host arcs, and an anti-correlation with mass-parameter hypervolume. A conjugate-point analysis using only quasar image positions yields substantially broader posteriors, confirming that the arcs provide the dominant constraints on the lens mass profiles. The work is positioned as preparation for a follow-up hierarchical cosmographic analysis to constrain H0.
Significance. If the reported correlation holds, the result is significant for time-delay cosmography because doubly imaged systems are far more abundant than quads yet have been underutilized due to fewer constraints. Demonstrating that arc surface brightness is the primary driver of Fermat-potential precision, backed by internal cross-checks (conjugate-point posteriors and hypervolume anti-correlation) and external consistency (Einstein radii and image positions), supplies a practical selection criterion for the large samples expected from LSST, Roman, and Euclid. The uniform, open-source pipeline is a clear strength that improves reproducibility across the TDCOSMO series.
major comments (2)
- [§4] §4 (correlation analysis): The central claim that arc surface brightness is the 'primary driver' of mass-model precision rests on a reported 'strong correlation,' yet no quantitative statistic (Spearman coefficient, Pearson r with uncertainty, or p-value) or explicit fit to the eight data points is provided; without this, it is impossible to judge whether the trend is statistically robust or dominated by outliers.
- [Conjugate-point analysis] Conjugate-point analysis section: The statement that posteriors are 'substantially broader' when only quasar positions are used is presented without a quantitative metric (e.g., ratio of credible-interval widths or hypervolume ratio) or the exact number and selection of conjugate points; this leaves the confirmation that arcs dominate the constraints qualitative rather than rigorous.
minor comments (2)
- [Abstract] Abstract: The phrase 'a subsequent publication' for the H0 analysis should include the planned paper number or title to aid readers following the TDCOSMO series.
- [Figure captions] Figure captions (e.g., those showing posterior contours): Captions should explicitly state the number of MCMC samples, burn-in length, and whether the displayed contours are 68 % and 95 % credible regions.
Simulated Author's Rebuttal
We thank the referee for the positive assessment of our manuscript and the recommendation for minor revision. We address each major comment below and will incorporate the suggested quantitative improvements.
read point-by-point responses
-
Referee: [§4] §4 (correlation analysis): The central claim that arc surface brightness is the 'primary driver' of mass-model precision rests on a reported 'strong correlation,' yet no quantitative statistic (Spearman coefficient, Pearson r with uncertainty, or p-value) or explicit fit to the eight data points is provided; without this, it is impossible to judge whether the trend is statistically robust or dominated by outliers.
Authors: We agree that quantitative statistics are needed to substantiate the correlation claim. In the revised manuscript we will report the Spearman rank correlation coefficient with its p-value for Fermat potential precision versus arc surface brightness. We will also add a linear fit to the eight data points, including uncertainties on the slope and intercept, to allow evaluation of robustness and outlier influence. These results will be included in §4. revision: yes
-
Referee: [Conjugate-point analysis] Conjugate-point analysis section: The statement that posteriors are 'substantially broader' when only quasar positions are used is presented without a quantitative metric (e.g., ratio of credible-interval widths or hypervolume ratio) or the exact number and selection of conjugate points; this leaves the confirmation that arcs dominate the constraints qualitative rather than rigorous.
Authors: We agree that quantitative metrics would make this section more rigorous. In the revision we will state the exact number of conjugate points (the four quasar image positions) and their selection criteria. We will report the ratios of 68% credible-interval widths for key parameters (Einstein radius, power-law slope) between the full and conjugate-point-only models, and the posterior hypervolume ratio where relevant. These details will be added to the conjugate-point analysis section. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper performs uniform Lenstronomy modeling on HST imaging of eight doubly-imaged quasars and reports an empirical correlation between Fermat-potential precision and host-arc surface brightness. This correlation is obtained directly from the posterior widths of the full-image fits and is cross-checked by a conjugate-point analysis that deliberately omits the arcs and recovers substantially broader posteriors. External validation is provided by Gaia DR2 image separations (3.6 mas rms) and literature Einstein radii (1.5σ average agreement). No self-definitional equations, fitted parameters renamed as predictions, or load-bearing self-citations appear in the derivation chain; the central claim rests on data-driven posteriors and independent external benchmarks rather than on any reduction to the paper's own inputs.
Assumptions & free parameters
free parameters (1)
- lens mass profile parameters
assumptions (2)
- standard math Standard general-relativity light deflection applies to these systems
- domain assumption Chosen mass profiles adequately describe the foreground galaxies
Cite this review
Pith. "Pith review of TDCOSMO XXVI: Uniform lens modeling of eight doubly imaged quasars." pith.science (2026). https://pith.science/paper/OQSZEZHU
@misc{pith2026260424908,
author = {Pith},
title = {Pith review of: TDCOSMO XXVI: Uniform lens modeling of eight doubly imaged quasars},
year = {2026},
howpublished = {\url{https://pith.science/paper/OQSZEZHU}},
note = {Machine review of arXiv:2604.24908}
}
abstract
We present the first uniform gravitational lens modeling analysis of eight doubly imaged quasars from multi-band observations with the Hubble Space Telescope. Previous time-delay cosmography analyses by the TDCOSMO Collaboration have primarily relied on quadruply imaged quasars, while doubly imaged systems, despite being more abundant, remain underutilized due to their fewer geometric constraints. Using an open-source $\texttt{Lenstronomy}$ framework, we reconstruct the lensing systems with a pipeline tailored for doubles. Comparing our results to the literature, the modeled Einstein radii agree at an average of 1.5$\sigma$, which is expected given data and modeling heterogeneity, while modeled image separations differ from Gaia DR2 measurements with an r.m.s of only 3.6 mas. We find a strong correlation between Fermat potential precision and the surface brightness of the spatially extended host arcs, establishing that arc surface brightness is the primary driver of mass model precision in doubly imaged systems. To further quantify the information contributed by the lensed arcs, we performed a conjugate point analysis that uses only the quasar image positions to constrain the lens mass profiles. The resulting posteriors are substantially broader than those from full image modeling, and a strong anti-correlation between mass parameter hypervolume and arc magnitude additionally confirms that arc brightness determines the degree to which the lens mass profile can be constrained in doubles. A hierarchical cosmographic analysis incorporating time-delay measurements and stellar kinematics to infer $\text{H}_0$ will be presented in a subsequent publication. The uniform pipeline and arc surface brightness trends established here will significantly accelerate the construction of time-delay cosmography samples from the large lens populations expected from LSST, Roman, and Euclid.
Forward citations
Cited by 1 Pith paper
-
Disentangling the dark and stellar mass through precise lens modelling of the JWST observation of lensed quasar WFI2033--4723
Composite time-delay lens modelling of JWST WFI2033–4723 yields intermediate stellar IMF normalisation and a steep gNFW inner slope γ_in≃1.3 that is robust to the choice of external cosmology.
Reviewed May 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.