Recognition: no theorem link
Supernova 2025wny: High-angular resolution Keck/NIRC2 observations and preliminary lens modeling
Pith reviewed 2026-05-13 20:39 UTC · model grok-4.3
The pith
High-resolution Keck imaging yields precise astrometry for the five images of lensed supernova 2025wny and consistent mass models from two independent codes.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
High-angular resolution Kp-band adaptive optics imaging from Keck/NIRC2 determines the precise astrometric positions of the five multiple images of SN 2025wny. Lens mass models using a singular isothermal ellipsoid, a singular isothermal sphere, and external shear in two independent codes reproduce the positions with sub-milli-arcsecond residuals. The inferred masses are M1 = 4.44^{+0.06}_{-0.05} x 10^{11} solar masses and M2 = 0.96^{+0.02}_{-0.02} x 10^{11} solar masses, with the primary lens velocity dispersion of 277.4^{+0.9}_{-0.7} km/s consistent with the DESI value of 298 +/- 37 km/s.
What carries the argument
Singular isothermal ellipsoid or sphere with external shear, constrained by sub-arcsecond astrometric positions from adaptive optics imaging and fitted independently in two lens modeling codes.
If this is right
- The system is suitable for time-delay cosmography studies.
- The inferred masses and velocity dispersion are consistent with prior LBT-based modeling using the same codes.
- One image shows an anomalous flux excess by a factor of 2-3 not accounted for by the smooth mass models.
Where Pith is reading between the lines
- Precise time-delay measurements in this system could supply an independent route to cosmological parameters.
- The flux anomaly may indicate microlensing or substructure that future monitoring could isolate.
- The dual-code, high-resolution approach demonstrated here can be applied to additional galaxy-scale lensed transients.
Load-bearing premise
The lens mass distribution can be adequately described by a singular isothermal ellipsoid or sphere plus external shear.
What would settle it
A direct dynamical or weak-lensing mass measurement that deviates by more than the reported uncertainties from the modeled values, or image-position residuals that remain larger than sub-milliarcseconds after the same modeling procedure.
Figures
read the original abstract
Multiply imaged, gravitationally lensed supernovae are rare but powerful tools for providing independent measurements on cosmological parameters. Supernova (SN) 2025wny ("SN Winny") is the first gravitationally-lensed Type I superluminous supernova and the first lensed supernova in a galaxy-scale system that is suitable for time-delay cosmography studies. In this work, we present high-resolution $K_p$-band adaptive optics imaging of SN Winny obtained with the near-infrared camera (NIRC2) on the W. M. Keck II telescope. With exquisite image quality (FWHM$\approx0.^{\prime\prime}065$) we determine and make use of the precise astrometric positions of the five multiple images as constraints for our lens mass models. With lenstronomy and Glee, we parameterize the total mass of the system with a singular isothermal ellipsoid, a singular isothermal sphere, and external shear. The two independent models are in excellent agreement and reproduce the observed image positions with sub-milli-arcsecond residuals. The inferred projected total masses enclosed within the Einstein radii of the primary and secondary lens galaxies are M$_1$ = 4.44$^{+0.06}_{-0.05}\times10^{11} M_\odot$ and M$_2$ = 0.96$^{+0.02}_{-0.02}\times10^{11} M_\odot$, respectively. Likewise, the inferred effective velocity dispersion of the primary lens is $\sigma_{1} = $ 277.4$^{+0.9}_{-0.7}$ km/s, consistent with the independent spectroscopic measurement made by DESI of $\sigma_{\star,1} = $ 298$\,\pm\,37$ km/s. Our modeling results are also consistent with previous results for the same system with data from the Large Binocular Telescope (LBT), using the same lens modeling codes. We also corroborate their finding that the SN multiple image A has an anomalous excess of flux by a factor of ~2-3 beyond what our smooth mass models predict.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports high-resolution Keck/NIRC2 adaptive-optics imaging of the gravitationally lensed Type I superluminous supernova SN 2025wny, providing precise astrometric positions for its five multiple images. These positions are used to constrain parametric lens models (singular isothermal ellipsoid and singular isothermal sphere plus external shear) implemented independently in lenstronomy and Glee. The models agree to sub-milliarcsecond residuals and yield projected enclosed masses M1 = 4.44^{+0.06}_{-0.05}×10^{11} M_⊙ and M2 = 0.96^{+0.02}_{-0.02}×10^{11} M_⊙ within the Einstein radii together with an effective velocity dispersion σ1 = 277.4^{+0.9}_{-0.7} km s^{-1} that is consistent with the DESI spectroscopic value. The results are also stated to be consistent with earlier LBT-based modeling of the same system.
Significance. If the modeling assumptions hold, the work supplies the first galaxy-scale lensed superluminous supernova suitable for time-delay cosmography, with mass and velocity-dispersion measurements that can serve as anchors for future cosmological analyses once time delays are measured. The use of two independent codes, sub-mas residuals, and cross-checks with spectroscopy and prior data constitute concrete strengths that increase the reliability of the reported quantities.
major comments (3)
- [Modeling section] Modeling section (lens parameterization): The analysis restricts the total mass distribution to singular isothermal ellipsoid and sphere profiles (γ = 2 exactly). No power-law or composite profiles with free density slope are explored. Because the enclosed mass within the Einstein radius scales directly with the slope for a fixed image configuration, even modest deviations (γ = 1.8 or 2.2) would rescale the reported M1 and M2 while still reproducing the observed positions to sub-mas precision.
- [Results] Results and error analysis: The abstract and summary provide no quantitative error budget separating statistical uncertainties from systematics arising from profile choice, external shear parameterization, or the unmodeled flux anomaly in image A. The quoted ±0.06/0.05 and ±0.02 uncertainties on M1 and M2 therefore cannot be assessed for robustness against the weakest assumption identified in the stress test.
- [Results] Consistency check with spectroscopy: The reported agreement between the modeled σ1 = 277.4 km s^{-1} and the DESI value 298 ± 37 km s^{-1} lies well within the large spectroscopic uncertainty. This comparison therefore supplies only weak validation of the isothermal assumption and does not mitigate the need for explicit tests with non-isothermal profiles.
minor comments (1)
- [Discussion] The flux excess in image A (factor 2–3) is acknowledged but not folded into the mass model; a brief quantitative statement on whether this excess could be absorbed by small adjustments to the shear or by adding a substructure component would improve completeness.
Simulated Author's Rebuttal
We thank the referee for their thorough and constructive review of our manuscript. We address each major comment below and have revised the manuscript to strengthen the analysis and clarify the limitations of the current modeling.
read point-by-point responses
-
Referee: [Modeling section] Modeling section (lens parameterization): The analysis restricts the total mass distribution to singular isothermal ellipsoid and sphere profiles (γ = 2 exactly). No power-law or composite profiles with free density slope are explored. Because the enclosed mass within the Einstein radius scales directly with the slope for a fixed image configuration, even modest deviations (γ = 1.8 or 2.2) would rescale the reported M1 and M2 while still reproducing the observed positions to sub-mas precision.
Authors: We agree that the isothermal assumption (γ = 2) is a limitation and that allowing a free power-law slope would better quantify systematic uncertainties on the enclosed masses. In the revised manuscript we will add power-law models (with free γ) implemented in both lenstronomy and Glee, using the same image positions as constraints. This will demonstrate the sensitivity of M1 and M2 to the density slope and provide a more robust error budget. We note that the isothermal profiles were chosen for consistency with prior LBT modeling and because they are standard for preliminary time-delay cosmography work, but we accept the referee’s point that explicit tests with non-isothermal profiles are needed. revision: yes
-
Referee: [Results] Results and error analysis: The abstract and summary provide no quantitative error budget separating statistical uncertainties from systematics arising from profile choice, external shear parameterization, or the unmodeled flux anomaly in image A. The quoted ±0.06/0.05 and ±0.02 uncertainties on M1 and M2 therefore cannot be assessed for robustness against the weakest assumption identified in the stress test.
Authors: We acknowledge that the current error bars reflect only the statistical uncertainties from the MCMC sampling under the isothermal assumption. In the revised manuscript we will expand the results section to present a quantitative error budget that separates statistical and systematic contributions. This will include the impact of profile choice (from the new power-law runs), variations in external shear parameterization, and a discussion of the flux anomaly in image A. The revised text will clearly state that the quoted uncertainties are statistical only and will report the additional systematic range obtained from the expanded model suite. revision: yes
-
Referee: [Results] Consistency check with spectroscopy: The reported agreement between the modeled σ1 = 277.4 km s^{-1} and the DESI value 298 ± 37 km s^{-1} lies well within the large spectroscopic uncertainty. This comparison therefore supplies only weak validation of the isothermal assumption and does not mitigate the need for explicit tests with non-isothermal profiles.
Authors: We agree that the large spectroscopic uncertainty makes the comparison a relatively weak validation of the isothermal assumption. In the revised manuscript we will clarify this point in the text, noting that the agreement is consistent but does not strongly constrain the density slope. The addition of power-law models will allow us to compare the inferred velocity dispersions across different slopes and thereby provide a stronger internal test of the modeling assumptions. revision: yes
Circularity Check
Minor self-citation to prior LBT modeling; central mass inference from image-position fits remains independent
full rationale
The derivation fits standard SIE/SIS+shear profiles in lenstronomy and Glee directly to the five observed image positions, yielding enclosed masses and velocity dispersion as free parameters of the fit. Image-position residuals are minimized by construction, but the reported masses (M1, M2) and σ1 are outputs, not inputs. The only self-reference is a consistency statement with prior LBT results using the same codes; this is not load-bearing for the Keck-derived values. No self-definitional loops, fitted inputs renamed as predictions, or ansatz smuggling appear in the provided text. The isothermal assumption is an explicit modeling choice whose bias risk is noted separately from circularity.
Axiom & Free-Parameter Ledger
free parameters (3)
- Einstein radius of primary lens
- External shear components
- Velocity dispersion scaling
axioms (2)
- domain assumption Lens mass distribution follows a singular isothermal ellipsoid or sphere profile
- domain assumption External shear is spatially uniform
Forward citations
Cited by 1 Pith paper
-
The Quadruply Lensed Supernova SN 2025wny: Implications for LSST
A protocol is outlined for maximizing LSST quadruply lensed supernova science through pre-listed hosts, real-time modeling, and immediate follow-up to capture all images and time delays.
Reference graph
Works this paper leans on
-
[1]
The Open Journal of Astrophysics , keywords =
Abe, K. T., Oguri, M., Birrer, S., et al. 2025, The Open Journal of Astrophysics, 8, 8, doi: 10.33232/001c.128482
-
[2]
2024, MNRAS, 531, 3509, doi: 10.1093/mnras/stae1356
Arendse, N., Dhawan, S., Sagu´ es Carracedo, A., et al. 2024, MNRAS, 531, 3509, doi: 10.1093/mnras/stae1356
-
[3]
Aryan, A., Chen, T. W., Lee, Y. H., et al. 2025, Transient Name Server AstroNote, 330, 1
work page 2025
-
[4]
2026, arXiv e-prints, arXiv:2601.18787, doi: 10.48550/arXiv.2601.18787
Baltasar, S., Ratier-Werbin, N., Huang, X., et al. 2026, arXiv e-prints, arXiv:2601.18787, doi: 10.48550/arXiv.2601.18787
-
[5]
Bellm, E. C., Kulkarni, S. R., Graham, M. J., et al. 2019, PASP, 131, 018002, doi: 10.1088/1538-3873/aaecbe
-
[6]
2018, Physics of the Dark Universe, 22, 189, doi: 10.1016/j.dark.2018.11.002
Birrer, S., & Amara, A. 2018, Physics of the Dark Universe, 22, 189, doi: 10.1016/j.dark.2018.11.002
-
[7]
Arendse, N. 2025, Philosophical Transactions of the Royal Society of London Series A, 383, 20240130, doi: 10.1098/rsta.2024.0130
-
[8]
Birrer, S., Shajib, A. J., Galan, A., et al. 2020, arXiv e-prints, arXiv:2007.02941. https://arxiv.org/abs/2007.02941
-
[9]
2024, SSRv, 220, 48, doi: 10.1007/s11214-024-01079-w
Birrer, S., Millon, M., Sluse, D., et al. 2024, SSRv, 220, 48, doi: 10.1007/s11214-024-01079-w
-
[10]
Bolton, A. S., Burles, S., Koopmans, L. V. E., et al. 2008, ApJ, 682, 964, doi: 10.1086/589327
-
[11]
Burke, W. L. 1981, ApJL, 244, L1, doi: 10.1086/183466 Ca˜ nameras, R., Schuldt, S., Suyu, S. H., et al. 2020, arXiv e-prints, arXiv:2004.13048. https://arxiv.org/abs/2004.13048
-
[12]
Chen, G. C.-F., Suyu, S. H., Wong, K. C., et al. 2016, MNRAS, 462, 3457, doi: 10.1093/mnras/stw991
-
[13]
Chen, J., Rozo, E., Dalal, N., & Taylor, J. E. 2007, ApJ, 659, 52, doi: 10.1086/512002
-
[14]
Dobler, G., & Keeton, C. R. 2006, MNRAS, 365, 1243, doi: 10.1111/j.1365-2966.2005.09809.x
-
[15]
2025, A&A, 697, A139, doi: 10.1051/0004-6361/202553807
Dux, F., Millon, M., Galan, A., et al. 2025, A&A, 697, A139, doi: 10.1051/0004-6361/202553807
-
[16]
Ecker, L. R., Schweinfurth, A. G., Saglia, R., et al. 2026, arXiv e-prints, arXiv:2602.16620, doi: 10.48550/arXiv.2602.16620
-
[17]
Goobar, A., Amanullah, R., Kulkarni, S. R., et al. 2017, Science, 356, 291, doi: 10.1126/science.aal2729
-
[18]
2022, arXiv e-prints, arXiv:2211.00656
Goobar, A., Johansson, J., Schulze, S., et al. 2022, arXiv e-prints, arXiv:2211.00656. https://arxiv.org/abs/2211.00656
-
[19]
Grayling, M., Thorp, S., Mandel, K. S., et al. 2026, MNRAS, doi: 10.1093/mnras/stag340
-
[20]
2008, A&A, 477, 397, doi: 10.1051/0004-6361:20077534
Grillo, C., Lombardi, M., & Bertin, G. 2008, A&A, 477, 397, doi: 10.1051/0004-6361:20077534
-
[21]
Grillo, C., Pagano, L., Rosati, P., & Suyu, S. H. 2024, A&A, 684, L23, doi: 10.1051/0004-6361/202449278 15
-
[22]
Huang, X., Storfer, C., Ravi, V., et al. 2020, The Astrophysical Journal, 894, 78, doi: 10.3847/1538-4357/ab7ffb
-
[23]
2021, ApJ, 909, 27, doi: 10.3847/1538-4357/abd62b
Huang, X., Storfer, C., Gu, A., et al. 2021, ApJ, 909, 27, doi: 10.3847/1538-4357/abd62b
-
[24]
Inoue, K. T. 2016, MNRAS, 461, 164, doi: 10.1093/mnras/stw1270
-
[25]
Jacobs, C., Collett, T., Glazebrook, K., et al. 2019, ApJS, 243, 17, doi: 10.3847/1538-4365/ab26b6
-
[26]
Johansson, J., Perley, D. A., Goobar, A., et al. 2025, ApJL, 995, L17, doi: 10.3847/2041-8213/ae1d61
-
[27]
Kelly, P. L., Rodney, S. A., Treu, T., et al. 2015, Science, 347, 1123, doi: 10.1126/science.aaa3350
-
[28]
L., Rodney, S., Treu, T., et al
Kelly, P. L., Rodney, S., Treu, T., et al. 2023, Science, 380, abh1322, doi: 10.1126/science.abh1322
-
[29]
2020, A&A, 640, A105, doi: 10.1051/0004-6361/202037740
Millon, M., Courbin, F., Bonvin, V., et al. 2020, A&A, 640, A105, doi: 10.1051/0004-6361/202037740
-
[30]
2022, emolter/nirc2 reduce: pre-alpha, 0.0.1, Zenodo, doi: 10.5281/zenodo.6584662
Molter, N. 2022, emolter/nirc2 reduce: pre-alpha, 0.0.1, Zenodo, doi: 10.5281/zenodo.6584662
-
[31]
Pascale, M., Frye, B. L., Pierel, J. D. R., et al. 2025, ApJ, 979, 13, doi: 10.3847/1538-4357/ad9928
-
[32]
Pierel, J. D. R., Frye, B. L., Pascale, M., et al. 2024, ApJ, 967, 50, doi: 10.3847/1538-4357/ad3c43
-
[33]
2022, PhRvD, 106, 023520, doi: 10.1103/PhysRevD.106.023520
Qi, J.-Z., Cui, Y., Hu, W.-H., et al. 2022, PhRvD, 106, 023520, doi: 10.1103/PhysRevD.106.023520
-
[34]
Rodney, S. A., Brammer, G. B., Pierel, J. D. R., et al. 2021, Nature Astronomy, 5, 1118, doi: 10.1038/s41550-021-01450-9 Sainz de Murieta, A., Collett, T. E., Magee, M. R., et al. 2024, MNRAS, 535, 2523, doi: 10.1093/mnras/stae2486
-
[35]
Schechter, P. L., & Wambsganss, J. 2002, ApJ, 580, 685, doi: 10.1086/343856
-
[36]
Service, M., Lu, J. R., Campbell, R., et al. 2016, PASP, 128, 095004, doi: 10.1088/1538-3873/128/967/095004
-
[37]
Shajib, A. J., Mozumdar, P., Chen, G. C.-F., et al. 2023, A&A, 673, A9, doi: 10.1051/0004-6361/202345878
-
[38]
2024, ApJS, 274, 16, doi: 10.3847/1538-4365/ad527e
Storfer, C., Huang, X., Gu, A., et al. 2024, ApJS, 274, 16, doi: 10.3847/1538-4365/ad527e
-
[39]
Storfer, C. J., Magnier, E. A., Huang, X., et al. 2026, ApJ, 999, 1, doi: 10.3847/1538-4357/ae381c
-
[40]
Suyu, S. H., & Halkola, A. 2010, A&A, 524, A94, doi: 10.1051/0004-6361/201015481
-
[41]
Suyu, S. H., Hensel, S. W., McKean, J. P., et al. 2012, ApJ, 750, 10, doi: 10.1088/0004-637X/750/1/10
-
[42]
Takahashi, R., & Inoue, K. T. 2014, MNRAS, 440, 870, doi: 10.1093/mnras/stu328
-
[43]
Taubenberger, S., Acebron, A., Ca˜ nameras, R., et al. 2025, arXiv e-prints, arXiv:2510.21694, doi: 10.48550/arXiv.2510.21694 TDCosmo Collaboration, Birrer, S., Buckley-Geer, E. J., et al. 2025, A&A, 704, A63, doi: 10.1051/0004-6361/202555801
-
[44]
Wong, K. C., Suyu, S. H., Chen, G. C.-F., et al. 2020, MNRAS, 498, 1420, doi: 10.1093/mnras/stz3094
-
[45]
Yelda, S., Lu, J. R., Ghez, A. M., et al. 2010, ApJ, 725, 331, doi: 10.1088/0004-637X/725/1/331
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.