Pith. sign in

REVIEW 4 major objections 3 minor 28 references

Model-independent and model-based local lensing properties of B0128+437 from resolved quasar images

T0 review · 4 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper argues that B0128+437's resolved milliarcsecond quasar subcomponents rule out smooth elliptical lens models and require free-form mass reconstructions.

desk verdict Careful, useful application of an existing method to a hard lens; the quantitative results are conditional on a subcomponent matching that is narrower than it should be, but the qualitative conclusion survives. read the letter →

arxiv 1909.01349 v1 pith:X6X5TAV5 submitted 2019-09-03 astro-ph.GA astro-ph.COastro-ph.IM

Jenny Wagner , Liliya L. R. Williams This is my paper · ORCID
classification astro-ph.GAastro-ph.COastro-ph.IM
keywords gravitationallensing:strongdarkmattermethods:analyticalgalaxies:individual:B0128+437quasars:generalmodel-independentlenspropertiesfree-formmassreconstructionVLBIsubcomponents
topics Dark Matter
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper sets out to characterize the galaxy-scale gravitational lens B0128+437 using the milliarcsecond radio subcomponents of its four quasar images. It argues that a model-independent extraction of local lens properties, specifically ratios of scaled mass densities (convergences) and reduced shears at the image positions, together with free-form pixel reconstructions, yields a consistent picture, while smooth elliptical parametric models fail once all subcomponents are included. The upshot is that the lens galaxy's mass distribution is asymmetric on small scales, with possible mass-density gradients on milliarcsecond scales that smooth models cannot capture.

What carries the argument

The load-bearing identity is the mass-sheet-invariant relation between multiple images: $f_{ij} = (1-\kappa_i)/(1-\kappa_j)$ and $g_i = \gamma_i/(1-\kappa_i)$, obtained by linearly mapping three reference points (subcomponent positions, or the endpoints of fitted Gaussian axes) from a reference image to the others. These are the leading-order local lens properties that every valid lens model must share. The comparison machinery has three legs: parametric potentials (softened power-law and boxy power-law potentials with external shear) fitted by a parametric modelling code; a free-form pixel mass reconstruction that averages many regularised solutions; and a model-free check based on the relative polar angles of the quads against the Fundamental Surface of Quads. The argument works by showing that the three legs agree on image scale while only the free-form and local descriptions capture the subcomponent scale.

What would settle it

A smooth, elliptically symmetric lens model (or a smooth power-law potential with external shear) that reproduces all 12 subcomponent positions to within about 0.1 milliarcseconds would falsify the paper's central claim; alternatively, multi-band observations that establish a different physical matching of the subcomponents would shift the derived convergence ratios and reduced shears, as the paper's own comparison of alternate matchings shows.

Watch

Extended reading notes

Core claim

B0128+437 is a quadruple-image quasar lens whose three images resolve into three VLBI subcomponents each, while the fourth image is scatter-broadened. The paper's central claim is that no smooth, elliptically symmetric mass model reproduces all subcomponent positions at sub-milliarcsecond precision, because the underlying mass distribution is asymmetric. Using a model-independent mapping of subcomponents, the paper obtains mass-sheet-invariant ratios $f_{ij}=(1-\kappa_i)/(1-\kappa_j)$ and reduced shears $g_i=\gamma_i/(1-\kappa_i)$ at image and subcomponent scales. These agree within 1-$\sigma$ with free-form pixel reconstructions at image scale, but only about 40 percent of the subcomponent-scale properties overlap in confidence bounds, so mass-density gradients on milliarcsecond scales cannot be excluded. The conclusion is that elliptical symmetry is too simplistic for this lenticular or late-type galaxy and that high-resolution observations demand flexible free-form models.

Load-bearing premise

The load-bearing assumption is that the subcomponents labelled 1, 2, and 3 in images A, C, D, and B correspond to the same source components; if dust changes their brightness ordering, a different physical matching could be correct and every derived local lens property would shift.

Editorial extensions

If this is right

  • Within 1-$\sigma$, the model-independent image-scale local lens properties agree with the free-form reconstruction values, so the fast local mapping can stand in for expensive global reconstructions on image scales.
  • Only 40 percent of the subcomponent-scale local lens properties overlap in confidence bounds, so milliarcsecond-scale mass-density gradients remain plausible and smooth models cannot be validated at that scale.
  • Parametric models fit only when a single subcomponent per image is used; adding all subcomponents degrades the fit, marking the lens mass distribution as asymmetric beyond elliptical symmetry.
  • The model-free polar-angle test places B0128 at effective ellipticity and shear around 0.25, and such a large shear value likely absorbs non-elliptical complexity rather than indicating a dominant nearby perturbing galaxy.
  • Together the methods give a consistent reconstruction of B0128 at current precision including the subcomponent structure, without invoking dark matter substructure as the primary cause of the flux anomalies.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If this pattern holds for other VLBI-resolved quads, survey pipelines that fit smooth elliptical or power-law lens models will systematically miss or misattribute milliarcsecond-scale mass complexity; quads that deviate from the Fundamental Surface of Quads are cheap screening targets for such asymmetry.
  • The 40 percent overlap statistic is a coarse diagnostic; a direct next step is to compute the same local lens properties with more than three reference points per image once deeper VLBI imaging resolves additional knots, which should shrink the confidence bounds and either confirm or dissolve the gradient signal.
  • The paper's interpretation of large external shear as a placeholder for internal small-scale asymmetry is testable by comparing shear estimates at different radii in other resolved lens systems: if the shear angle and magnitude change with scale, smooth-model shear is absorbing substructure rather than measuring a distant perturbing galaxy.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 3 minor

Summary. This paper applies the model-independent local lens-property formalism of Wagner & Tessore (2018) to the quad lens B0128+437 using VLBI-resolved subcomponents in images A, B, C, and D. It derives convergence ratios f_j and reduced shear components g_j on image and subcomponent scales, constructs PixeLens free-form and Lensmodel parametric reconstructions, and performs a Fundamental Surface of Quads (FSQ) symmetry analysis. The authors conclude that elliptically symmetric models are too simplistic to characterise the asymmetric mass distribution of B0128, and that milli-arcsecond-scale mass-density gradients cannot be excluded from the limited overlap of subcomponent-scale local lens properties.

Significance. The paper has clear strengths: the analysis is blinded through an independent mediator, alternative subcomponent matchings are systematically explored in Table 4 and Appendix B, observed flux ratios are held back as a consistency check rather than used as constraints, and the model-independent code is publicly available. If the conclusions stand, the work provides an efficient route to local lens properties on sub-component scales and a concrete example where free-form reconstructions are needed. The evidential weight is nevertheless moderate, because the central conclusion is conditional on the adopted subcomponent labelling and on several assumed positional uncertainties.

major comments (4)
  1. [Abstract and Section 8.1] The abstract states that 'only 40% of the small-scale subcomponent local lens properties overlap within the 1-sigma confidence bounds', while Section 8.1 states that the subcomponent-level local lens properties 'only overlap in 50% of the cases within their 1-sigma confidence bounds' and that the same degree of agreement is found against the Lensmodel reconstruction. This quantitative statement is used to support the milli-arcsecond gradient conclusion, so the manuscript must specify which comparison the 40% refers to and make the abstract consistent with the body.
  2. [Section 4.2.1, Table 4, Appendix B] The choice of the Biggs et al. (2004) labelling as the only viable matching is made because alternative matchings fail to produce the assumed relative parities of J_C and J_D. The parity pattern itself (A and C same parity, D opposite) is stated as an assumption in Section 3, not derived from the VLBI morphology or from an independent parity measurement. Since the PixeLens, Lensmodel, and FSQ analyses all inherit this labelling, the central conclusion that elliptically symmetric models are too simplistic is not model-independent with respect to the matching. Please test parity-consistent alternative matchings, including permutations that move the inner subcomponent 2, and show whether any such matching can be accommodated by an elliptical potential with external shear.
  3. [Section 5.2 and Table 6] Section 5.2 says that the radio data provide an astrometric precision of 0.01 milliarcsecond, while Section 3 and Table 2 state an assumed uncertainty of 0.1 mas for subcomponent positions in images A, C, and D and 1 mas for image B. The Lensmodel rms values in Table 6 (0.0005-0.0028 arcsec) are judged against this assumed uncertainty. The factor-of-ten discrepancy changes the assessment of whether simple parametric models are ruled out, and the positional error budget used for the Lensmodel comparison must be corrected and justified.
  4. [Section 3, Table 2, and Appendix C] The three subcomponent positions of image B are read off Figure 6 of Biggs et al. (2004) with an assumed 1 mas uncertainty, and the four-image results in Table 5 as well as the ABCD PixeLens and Lensmodel comparisons depend on these values. The 3-mas robustness test in Appendix C changes the most likely values of f_B, g_B, and J_B substantially; for example, f_B changes from -16.75 in the fiducial case to -2.50 or -0.72 depending on configuration. The paper should state explicitly whether these shifts are within the quoted confidence bounds and whether the parity-related conclusions for image B are robust to the assumed B-position uncertainty.
minor comments (3)
  1. [Section 8.2] Section 8.2 contains an unremoved author note: 'does that sound better? for me, it's strange to produce something within some error' after the statement about Lensmodel being unable to reproduce images within the astrometric precision. This note must be removed and the sentence rephrased as a finished claim.
  2. [Sections 5.2 and 8.2] There are typographical errors: 'astromentric' in Section 5.2 and 'illustrateed' in Section 8.2 should be corrected to 'astrometric' and 'illustrated'.
  3. [Table 7 and Section 6.2] The PixeLens rms uncertainties in Table 7 are computed as the dispersion between 20 sets of 10 models, but the text does not state whether the quoted local lens properties are means over the 200 accepted models or over the 20 averaged sets; please clarify this in the table caption or text.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the local lens properties are measured from image positions, flux ratios are withheld as a consistency check, and the central claim is supported by independent parametric/free-form reconstructions.

full rationale

The model-independent f_ij and g_i values in Section 4 are derived from the observed subcomponent positions via the linear-mapping method of Wagner & Tessore (2018), not fitted to the paper's conclusions. The observed flux ratios are explicitly not used as constraints and are compared only afterwards as a consistency check, which is the opposite of fitted-input-called-prediction. The Lensmodel, PixeLens, and relative-polar-angle analyses are independently computed: Lensmodel's failure to reproduce all 12 subcomponents is a direct model-fit output, PixeLens reconstructs convergence maps from the same positions under its own stated regularisation constraints, and the FSQ comparison uses observed polar angles against an externally established surface. The only selection step is the subcomponent-labelling choice in Section 4.2.1: alternative swaps of labels 1 and 3 are rejected because they do not yield the assumed relative parities of images A, C, and D. That parity pattern is stated in Section 3 as an external lensing-theory input, not as an output of the analysis, so rejecting matchings that violate it is a consistency requirement rather than a circular reduction. The normalizations f_A = J_A = 1 are definitional but carry no evidential weight. Self-citations to Wagner & Tessore (2018), Wagner et al. (2018), Saha & Williams (2004), and Woldesenbet & Williams (2012) supply algorithms and comparison tools; the central astrophysical claim about asymmetry does not reduce to those citations because the model failures and reconstructed maps are presented as independent, falsifiable results. Thus no circular step is present.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The analysis introduces no new physical entities; it combines standard lensing theory, published observations, and the authors' earlier model-independent formalism. The free parameters are the fitted Lensmodel parameters, PixeLens regularization choices, and the assumed uncertainty for image B positions.

free parameters (3)
  • Lensmodel potential parameters (normalization, core radius, axis ratio, power-law slope, ellipticity, position angle… = e.g., ellipticity and external shear near 0.22-0.26 in fitted models (Section 7.2)
    These parameters are fitted by Lensmodel to the image positions; the fitted values are used to argue that large external shear is needed, but the exact values are not central to the qualitative conclusion.
  • PixeLens pixel scale and regularization constraints = pixel edge 33 mas; mass gradient within 45 degrees of radial; no pixel mass exceeds twice the average of neighbors
    Hand-chosen regularization choices that shape the free-form mass distribution and prevent extraction of subcomponent-scale properties.
  • Assumed positional uncertainty for image B subcomponents = 1 mas, with a 3 mas variant tested in Appendix C
    Because the B subcomponent positions are read off Figure 6 of Biggs et al. (2004) rather than tabulated, the ABCD analysis depends on this assumed input; the ACD analysis does not.
assumptions (5)
  • domain assumption The standard single-plane gravitational lensing formalism with convergence, shear, and deflection potential applies to B0128.
    Invoked throughout Section 4.1 and used to define f_ij and g_i.
  • domain assumption Lens redshift z_l = 1.145 and source redshift z_s = 3.124.
    Section 3 adopts the best-fit redshifts; the critical density entering the convergence ratios depends on these values.
  • domain assumption The subcomponents in images A, B, C, D are matched according to Biggs et al. (2004), selected because this labelling yields the correct relative parities of the magnifications.
    Sections 3 and 4.2.1; alternative matchings change the derived local lens properties substantially (Table B.1).
  • domain assumption Convergence and shear are approximately constant over the area spanned by the reference points in each image or Gaussian subcomponent.
    Section 4.1 states this assumption; it justifies the linear mapping between images.
  • domain assumption Image B is scatter-broadened and its subcomponent positions carry an assumed 1 mas uncertainty.
    Section 3 and Table 2 notes; the impact of a larger uncertainty is tested in Appendix C.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Model-independent and model-based local lensing properties of B0128+437 from resolved quasar images." pith.science (2026). https://pith.science/paper/X6X5TAV5

@misc{pith2026190901349,
  author       = {Pith},
  title        = {Pith review of: Model-independent and model-based local lensing properties of B0128+437 from resolved quasar images},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/X6X5TAV5}},
  note         = {Machine review of arXiv:1909.01349}
}
abstract

The galaxy-scale gravitational lens B0128+437 generates a quadrupole-image configuration of a background quasar that shows milli-arcsecond-scale subcomponents in the multiple images observed with VLBI. As this multiple-image configuration including the subcomponents has eluded a parametric lens-model characterisation so far, we determine local lens properties at the positions of the multiple images with our model-independent approach. Using PixeLens, we also succeed in setting up a global free-form mass density reconstruction including all subcomponents as constraints. We compare the model-independent local lens properties with those obtained by PixeLens and those obtained by the parametric modelling algorithm Lensmodel. A comparison of all three approaches and a model-free analysis based on the relative polar angles of the multiple images corroborate the hypothesis that elliptically symmetric models are too simplistic to characterise the asymmetric mass density distribution of this lenticular or late-type galaxy. In addition, the model-independent approach efficiently determines local lens properties on the scale of the quasar subcomponents, which are computationally intensive to obtain by free-form model-based approaches. As only 40% of the small-scale subcomponent local lens properties overlap within the 1-$\sigma$ confidence bounds, mass density gradients on milli-arcsecond scales cannot be excluded. Hence, aiming at a global reconstruction of the deflecting mass density distribution, increasingly detailed observations require flexible free-form models that allow for density fluctuations on milli-arcsecond scale to replace parametric ones, especially for asymmetric lenses or lenses with localised inhomogeneities like B0128.

Figures

Figures reproduced from arXiv: 1909.01349 by the authors.

Figure 2
Figure 2. Like the lens models by Norbury (2002), the lens mod [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 1
Figure 1. Left: HST I-band observation from Norbury (2002); Right: MERLIN 5GHz observation by Phillips et al. (2000). [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. VLBA 8.4 GHz details of all four multiple images from Biggs et al. (2004). [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figures from the paper (3 more)
Figure 3
Figure 3. Figure 3: Visualisation of image-scale matching: using the posi [PITH_FULL_IMAGE:figures/full_fig_p006_3.png]
Figure 4
Figure 4. Figure 4: PixeLens reconstructed maps of lensing convergence. Left: using all 12 subcomponents of B0128; Right: using only the 9 subcomponents of images A, C, and D. The subcomponents used as constraints in each case are indicated with magenta circles. The grey iso-convergence c…
Figure 5
Figure 5. Figure 5: The result of subtracting mass density distribution of con [PITH_FULL_IMAGE:figures/full_fig_p010_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

28 extracted references · 20 canonical work pages

  1. [1]

    D., Browne, I

    Biggs, A. D., Browne, I. W. A., Jackson, N. J., et al. 2004, MNRAS, 350, 949

    work page 2004

  2. [2]

    S., Burles, S., Koopmans, L

    Bolton, A. S., Burles, S., Koopmans, L. V . E., et al. 2008, ApJ, 682, 964

    2008

  3. [3]

    D., Kochanek, C

    Cohn, J. D., Kochanek, C. S., McLeod, B. A., & Keeton, C. R. 2001, ApJ, 554, 1216

    work page 2001

  4. [4]

    2019, MNRAS, 1618

    Gilman, D., Birrer, S., Treu, T., Nierenberg, A., & Benson, A. 2019, MNRAS, 1618

    work page 2019

  5. [5]

    Gomer, M. R. & Williams, L. L. R. 2018, MNRAS, 475, 1987

    work page 2018

  6. [6]

    2018, Monthly Notices of the Royal Astronomical Society, 475, 2438

    Hsueh, J.-W., Despali, G., Vegetti, S., et al. 2018, Monthly Notices of the Royal Astronomical Society, 475, 2438

    work page 2018

  7. [7]

    2004, gravlens 1.06 Software for Gravitational Lensing, 9th edn., uRL: http://www.physics.rutgers.edu/ keeton/gravlens/manual.pdf

    Keeton, C. 2004, gravlens 1.06 Software for Gravitational Lensing, 9th edn., uRL: http://www.physics.rutgers.edu/ keeton/gravlens/manual.pdf

    work page 2004

  8. [8]

    Keeton, C. R. 2001, ArXiv Astrophysics e-prints [astro-ph/0102341]

    arXiv 2001

Show all 28 references
  1. [9]

    J., Auger, M

    Lagattuta, D. J., Auger, M. W., & Fassnacht, C. D. 2010, ApJ, 716, L185

    2010

  2. [10]

    P., Koopmans, L

    McKean, J. P., Koopmans, L. V . E., Browne, I. W. A., et al. 2004, MNRAS, 350, 167

    2004

  3. [11]

    W., Massey, R

    Nightingale, J. W., Massey, R. J., Harvey, D. R., et al. 2019, arXiv e-prints [arXiv:1901.07801]

    2019 arXiv

  4. [12]

    Norbury, M. A. 2002, PhD thesis, The University of Manchester (United King- dom

    2002

  5. [13]

    O., Levine, H., & Wambsganss, J

    Petters, A. O., Levine, H., & Wambsganss, J. 2001, Singularity Theory and Grav- itational Lensing, Progress in Mathematical Physics, V olume 21 (Birkhäuser)

    2001

  6. [14]

    M., Norbury, M

    Phillips, P. M., Norbury, M. A., Koopmans, L. V . E., et al. 2000, MNRAS, 319, L7

    2000

  7. [15]

    S., Norbury, M., et al

    Rusin, D., Kochanek, C. S., Norbury, M., et al. 2001, ApJ, 557, 594

    2001

  8. [16]

    & Williams, L

    Saha, P. & Williams, L. L. R. 2004, AJ, 127, 2604

    2004

  9. [17]

    Schneider, P., Ehlers, J., & Falco, E. E. 1992, Gravitational Lenses, Astronomy and Astrophysics Library (New York: Springer)

    1992

  10. [18]

    2012, A&A, 538, A99

    Sluse, D., Chantry, V ., Magain, P., Courbin, F., & Meylan, G. 2012, A&A, 538, A99

    2012

  11. [19]

    H., Hensel, S

    Suyu, S. H., Hensel, S. W., McKean, J. P., et al. 2012, ApJ, 750, 10

    2012

  12. [20]

    2017, A&A, 597, L1

    Tessore, N. 2017, A&A, 597, L1

    2017

  13. [21]

    2018, A&A, 620, A86

    Wagner, J. 2018, A&A, 620, A86

    2018

  14. [22]

    2019, MNRAS, 487, 4492

    Wagner, J. 2019, MNRAS, 487, 4492

    2019

  15. [23]

    2018, A&A, 612, A17

    Wagner, J., Liesenborgs, J., & Tessore, N. 2018, A&A, 612, A17

    2018

  16. [24]

    & Tessore, N

    Wagner, J. & Tessore, N. 2018, A&A, 613, A6

    2018

  17. [25]

    Williams, L. L. R. & Liesenborgs, J. 2019, MNRAS, 482, 5666

    2019

  18. [26]

    Woldesenbet, A. G. & Williams, L. L. R. 2012, MNRAS, 420, 2944

    2012

  19. [27]

    Woldesenbet, A. G. & Williams, L. L. R. 2015, MNRAS, 454, 862

    2015

  20. [28]

    2015, MNRAS, 447, 3189 Appendix A: Derivation of the reference point positions on the subcomponent-scale Appendix A.1: Derivation of coordinate positions Fig

    Xu, D., Sluse, D., Gao, L., et al. 2015, MNRAS, 447, 3189 Appendix A: Derivation of the reference point positions on the subcomponent-scale Appendix A.1: Derivation of coordinate positions Fig. A.1: Determining the reference points from an elliptical Gaussian fitted to the subc...

    2015

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.