Pith. sign in

REVIEW 2 major objections 5 minor 133 references

Third- and fourth-order multipoles plus radial iso-density twists can explain B1422+231 flux ratios without dark-matter clumps.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · grok-4.5

2026-07-14 04:45 UTC pith:66QBFR4L

load-bearing objection Clean existence proof that TNG multipoles + radial twists can fit B1422 flux ratios at 2% precision, but the successful macros sit in extreme tails and the sampling of radial variations is coarse. the 2 major comments →

arxiv 2607.11559 v1 pith:66QBFR4L submitted 2026-07-13 astro-ph.GA astro-ph.CO

Can third- and fourth-order multipoles plus radial variation of iso-density ellipses explain the observed flux ratios in B1422+231? YES, and a lesson learned from a TNG100 lensing galaxy sample

Ruizhe Feng , Dandan Xu , Dominique Sluse , Giulia Despali , Anowar Shajib , Cai-Na Hao This is my paper
classification astro-ph.GA astro-ph.CO
keywords strong gravitational lensingflux-ratio anomaliesmultipole perturbationsiso-density twistsB1422+231TNG100macro-model degeneracy
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Flux-ratio anomalies in quadruply imaged quasars have long been read as signs of clumpy dark-matter substructure. This paper tests whether ordinary, non-clumpy deviations from perfect ellipses—third- and fourth-order multipoles together with radial changes in axis ratio and position angle—can do the same job for the classic system B1422+231. The authors pull realistic versions of these perturbations from a sample of strong-lensing galaxies in the TNG100 simulation, then add them to standard smooth models (SIE+γ and EPL+γ) while fitting the observed image positions and fluxes at several precision levels. They show that, once photometric precision reaches 2 percent, the smooth models fail, yet the same models plus the four extracted perturbation types recover the data without any subhalos. The work also demonstrates that model flexibility and hidden degeneracies matter: an EPL macro-model is far more forgiving than an isothermal one, and treating multipoles as free global parameters can artificially break degeneracies that real galaxies keep coupled.

Core claim

Typical macroscopic, non-clumpy density perturbations—m3 and m4 multipoles plus radial variations of iso-density ellipses—extracted from TNG100 strong lenses can rescue an EPL+γ macro-model and fully reproduce the image positions and anomalous flux ratios of B1422+231 at 2 percent photometric precision, without any clumpy mass components.

What carries the argument

Simulation-based “coarse-grain” sampling of the high-dimensional perturbation space: multipole amplitudes and radial Δq, Δϕq profiles are extracted from 536 TNG100 projections and grafted onto SIE/EPL+γ macros via multi-slice and multipole lens models, then tested for successful 3-σ recovery of the observed configuration.

Load-bearing premise

The 536 simulated galaxy projections sample the high-dimensional space of radial ellipticity and twist variations densely enough that a non-detection can be trusted as a physical impossibility rather than incomplete coverage.

What would settle it

A larger simulation sample or a direct multi-slice reconstruction of B1422 that still cannot recover the 2 percent flux ratios within 3σ after all four perturbation types are freely varied would falsify the claim that these macroscopic features suffice.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • At current 2 percent flux-ratio precision, B1422-like anomalies need not imply subhalos once realistic multipoles and iso-density twists are allowed.
  • SIE+γ should be avoided as the default macro-model for flux-ratio studies; free radial slope already absorbs much of the anomaly.
  • Global m3/m4 multipoles alone can artificially break degeneracies that real galaxies keep coupled to radial shape variation, biasing dark-matter inferences.
  • Success-rate differences between multipole-only and twist-inclusive models largely reflect sampling density, not relative physical importance.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same machinery applied to a statistical sample of quads could quantify how often macroscopic perturbations alone erase the need for substructure constraints.
  • If future high-resolution imaging of B1422 reveals strong radial isophote twists, the paper’s successful multi-slice solutions become the preferred null hypothesis before subhalo searches.
  • The demonstrated macro-model degeneracy implies that many existing subhalo mass-function limits may need re-derivation with free radial slopes and full multipole-plus-twist freedom.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 5 minor

Summary. The paper asks whether non-clumpy macroscopic perturbations (global m3 and m4 multipoles plus radial variations of iso-density ellipticity and position angle) extracted from a TNG100 strong-lensing sample can reproduce the image positions and flux ratios of B1422+231 without clumpy substructure. Perturbations are measured via SPH surface-density maps and photutils isophote fits, then added to SIE+γ and EPL+γ macros (global multipoles via EPL_MULTIPOLE_M3M4; radial Δq/Δφq via a 1000-slice ElliSLICE construction). Macro parameters are re-optimized against MERLIN positions (and fluxes) under two astrometric and three photometric precision levels. At σp=10 mas both macros fit positions alone; at 2 mas astrometric anomalies appear and are partially rescued by the perturbations. With positions+fluxes, SIE+γ fails even at 10% photometry, while EPL+γ alone succeeds at 10% and 5% but fails at 2%; adding all four perturbation types yields a small number of successful EPL+γ analogues (Table 4, last row). The authors emphasize model flexibility, degeneracy, and the risk of artificially breaking degeneracies by adding only global multipoles.

Significance. If the existence result holds, it supplies a concrete, simulation-calibrated demonstration that typical non-clumpy azimuthal structure can produce the classic B1422 flux-ratio anomaly at the precision of current radio/mid-IR data, without dark-matter subhalos. The public code, transparent extraction pipeline, and systematic comparison of SIE versus EPL under controlled precision levels are genuine strengths. The caution against treating global m3/m4 as independent of radial ellipticity twists is timely for JWST and VLBI flux-ratio analyses that already include multipoles when constraining subhalo populations and dark-matter models. The work is an existence proof rather than a population inference, but that is still a useful contribution to the modeling literature.

major comments (2)
  1. Table 4 (2%-EPL row) and Fig. 12 / §5: the 10 successful EPL+γ+m3+m4+PAv+Qv cases (and the multipole-only successes) recover macros deep in the tails of both the adopted priors and the TNG/SL2S population (s~2.35, q~0.5, γ~0.27). The paper itself flags physical plausibility. Because the macro is freely re-optimized for every extracted perturbation, the “rescue” may be achieved only by driving the smooth component into a regime that a more realistic multi-component or non-power-law macro would not require. A quantitative comparison of these best-fit macros against the TNG sample’s own radial slopes, ellipticities and shear proxies (or an explicit statement that the YES claim is purely an existence result under the adopted EPL family) is needed to keep the central claim proportionate.
  2. §2.3–2.4 and §5 (discussion of Table 4 success rates): multipoles are implemented as globally averaged, radially constant m3/m4, while Δq/Δφq are fully radial. The paper notes that low success rates for PAv/Qv combinations may simply reflect incomplete sampling of a high-dimensional space (536 projections). This asymmetry is load-bearing for the claim that “all four types” can rescue the model: the few successes may be rare TNG projections whose large |am/a| and |Δq| happen to compensate an already extreme macro. Either a controlled test with radially varying multipoles, or a clearer quantification of how densely the radial-variation space is sampled, is required before the relative importance of multipoles versus ellipse twists can be interpreted.
minor comments (5)
  1. §2.5 / Table 2: the Gaussian prior on q and φq combines Impey et al. light measurements with the TNG mass–light offset; the absolute cut q∈[0.4,1] is stated but the impact of that hard boundary on the extreme-q solutions of Fig. 12 is not shown.
  2. Fig. 10 versus Fig. 12: the striking difference in recovered macros when fitting positions alone versus positions+fluxes is important; a short quantitative statement of how often the position-only solutions already lie near the flux-ratio solutions would help the reader.
  3. §3.1 and Fig. 5: the comparison of a4/a to Hao et al. (2006a) is useful; note explicitly that the annulus used here (0.8–1.2 Rein) differs from theirs, as already mentioned in the footnote.
  4. Typos / notation: “as as the center” (§2.5); occasional switches between φq and φa; ensure consistent use of “principal ellipse” versus “globally averaged” multipoles.
  5. The public repository is a clear plus; a short README note on which exact TNG snapshot and subhalo IDs produce the 10 successful full-perturbation models would aid reproducibility.

Circularity Check

1 steps flagged

No load-bearing circularity: TNG-extracted perturbations are independent of B1422 data; macros are re-optimized against external observations, yielding an existence result rather than a forced prediction.

specific steps
  1. other [Sec. 2.5 + Sec. 3.1 / Fig. 4]
    "we consider the Gaussian priors on q and ϕq to be composed of two parts: one is from the light distribution of the lensing galaxy in C. D. Impey et al. (1996)... the other comes from the simulated misalignment between the mass and light elliptical shapes, which is given by qmass-qlight=0.08±0.06 and ϕq,mass-ϕq,light=0°±3.5° (see Section 3.1 and Fig. 4)."

    The identical TNG100 sample that later supplies the m3/m4/Δq/Δϕq perturbations is also used to construct the mass–light prior. This is a minor self-reference, not a definitional loop: the prior only regularizes the macro search and is not required for the existence of successful fits (which in any case recover macros deep in the prior tails).

full rationale

The derivation chain is: (1) select TNG100 projections by independent photometric/kinematic cuts resembling SL2S; (2) extract m3/m4 amplitudes and radial Δq/Δϕq via isophote fitting (photutils/Ellipse) on the particle density maps; (3) for each fixed extracted perturbation, re-optimize the free macro parameters (θE, q, ϕq, s, γ, ϕγ, source position) of SIE+γ or EPL+γ against the external MERLIN positions and flux ratios of B1422, accepting only those that land inside 3σ. Success is therefore an existence statement over a coarse-grained sample of realistic non-clumpy features, not a prediction that reduces to the fit by construction. The sole mild self-reference is that the same TNG sample also supplies the mass–light misalignment used to set the Gaussian prior on q and ϕq (Sec. 2.5 + Fig. 4); that prior is not load-bearing for the YES claim (the successful high-precision models sit far in the tails of the prior anyway). No equation equates the target flux ratios to a quantity defined by the fit itself, no uniqueness theorem is imported from overlapping authors, and no ansatz is smuggled via self-citation. Extreme recovered macros (s~2.35, q~0.5, γ~0.27) raise a physical-plausibility question, not a circularity one.

Axiom & Free-Parameter Ledger

4 free parameters · 4 axioms · 0 invented entities

The claim rests on standard strong-lensing formalism, the assumption that TNG100 galaxies supply realistic macroscopic perturbations, and a set of conventional macro-model parameters that are fitted rather than predicted. No new physical entities are introduced; the multi-slice technique is taken from Schramm (1994). The free parameters are the usual lens-model degrees of freedom plus a few numerical choices (number of slices, averaging annulus) that do not drive the existence result.

free parameters (4)
  • EPL/SIE macro parameters (θE, q, ϕq, s, γ, ϕγ, source position)
    Fitted by χ² minimization to the observed image positions (and optionally flux ratios) for each perturbation realization; priors listed in Table 2.
  • Gaussian prior means/widths on q and ϕq
    Constructed from Impey et al. light photometry plus the mass–light offset measured from the same TNG sample (Section 2.5); the prior is informative and affects the EPL solutions.
  • Number of elliptical slices (1000) and radial range (0.001–5.5 arcsec)
    Numerical discretization choice for the multi-slice implementation of radial Δq/Δϕq; stated to be dense enough to recover the unperturbed profile.
  • Averaging annulus 0.8–1.2 Rein for global multipoles and principal ellipse
    Choice of radial window used to extract a3/a, a4/a, q and ϕq from each TNG map; results are shown to be stable under a wider annulus but remain a free methodological choice.
axioms (4)
  • standard math Standard thin-lens equation and magnification under the point-source approximation
    Used throughout Sections 2.4–2.5 to compute image positions and flux ratios.
  • domain assumption TNG100 projected mass maps (after SPH smoothing) supply statistically realistic macroscopic multipoles and radial ellipticity twists for real strong lenses
    Core premise of the entire experiment (Sections 2.2–2.3); sample selection is matched to SL2S-like galaxies but remains a simulation-to-reality assumption.
  • ad hoc to paper Global (radially constant) m3 and m4 multipoles plus multi-slice Δq/Δϕq adequately capture the relevant non-elliptical complexity
    Explicit modeling choice (Section 2.4); radial variation of the multipoles themselves is deliberately omitted.
  • domain assumption Lens center coincides with the simulated galaxy center and is fixed at (0,0)
    Stated in Section 2.5; removes two free parameters but assumes perfect centering.

pith-pipeline@v1.1.0-grok45 · 36400 in / 3124 out tokens · 42161 ms · 2026-07-14T04:45:55.481850+00:00 · methodology

0 comments
read the original abstract

Flux ratio anomalies in multiply-imaged quasar lenses are a long-standing issue. Using a classical system B1422+231 as a case study, we investigate how typical non-clumpy perturbations beyond elliptical shapes -- multipoles $m_3, m_4$ and radial variations in $q, \phi_q$ -- can account for the observed image positions and flux ratios under different observational precisions. We extract these perturbations from a pre-selected strong-lensing galaxy sample from the TNG100 simulation. Smooth macroscopic models (SIE+$\gamma$, EPL+$\gamma$) are then fitted to the observed image positions alone and to both positions and flux ratios, with and without including the extracted perturbations. With astrometric uncertainty of $\sigma_{p}=10$ mas, both macro-models alone can already successfully fit image positions within $3\sigma_{p}$. At $\sigma_{p}=2$ mas, however, 'astrometric anomalies' appear if smooth macro-models alone are adopted. In this case, adding the extracted perturbations can explain the anomalous image positions. When both positions and flux ratios are adopted, the SIE+$\gamma$ model family already shows 'flux ratio anomalies' at photometric uncertainty $\sigma_{f} \le 10\%$ (keeping $\sigma_{p}=10$ mas). When EPL+$\gamma$ is used, the smooth model alone can simultaneously fit both positions and flux ratios with $\sigma_{f}=10\%, 5\%$, but not with $\sigma_{f}=2\%$, where 'flux ratio anomalies' appear. Adding all four types of extracted perturbations can rescue the macro-models and explain the observed anomalous flux ratios. We present important lessons learned regarding model flexibility and degeneracy.

Figures

Figures reproduced from arXiv: 2607.11559 by Anowar Shajib, Cai-Na Hao, Dandan Xu, Dominique Sluse, Giulia Despali, Ruizhe Feng.

Figure 1
Figure 1. Figure 1: Image configuration of B1422+231 as in [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Properties of our lensing galaxy sample at z = 0.3 (snapshot 078) from TNG100-1 simulation. Left: distribution of velocity dispersions and stellar masses of our sample, with results of SL2S (A. Sonnenfeld et al. 2013a,b) and C. N. Hao et al. (2006a) overlapped. The stellar masses of SL2S are estimated with two different initial mass functions (IMF). Note that the aperture of velocity dispersion in C. N. Ha… view at source ↗
Figure 3
Figure 3. Figure 3: An example of isophote-fitting results of simu￾lation galaxies (subhalo ID: 208612, projected along z-axis). Background grey-scale map: r-band light distribution; red lines: iso-density contours of total mass density distributions; black lines: elliptical iso-density contours fitted to the mass map. example of visualization, [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Difference in axis ratio q and position angle ϕq between mass and light. Left: principal values of axis ratio q for mass and light maps of our galaxy sample, with both histograms added on the top panel. The magenta dashed line means y = x. Right: difference of averaged ϕq and q (within two annuli) between mass and light maps for our sample. Dashed lines indicate 16%, 50%, and 84% quartiles. 0 50 PDF 0.04 0… view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of a4/a between the mass and light distributions, for our full galaxy sample (black points and histograms) and a subsample with velocity dispersion σv ⩾ 200 km/s (blue points and histograms). The magenta dashed line means y = x, and the red histogram is the dis￾tribution from C. N. Hao et al. (2006a). the observation. In addition, at q ≲ 0.55 (a minor frac￾tion at the tail of q), the simulated g… view at source ↗
Figure 6
Figure 6. Figure 6: Relations between multipole strengths from total mass maps and velocity dispersion. Left: a3/a VS σv; right: a4/a VS σv. For comparison, the multipole strengths are averaged within 2 different annuli (as the right panel of [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Relation between m4 strength a4/a and axis ratio q. Blue points and histograms: our sample; red points and histograms: C. N. Hao et al. (2006a); grey points and shaded histogram: random sample generated by prior in M. S. H. Oh et al. (2026). See also [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Relations between ∆q, ∆ϕq (difference of q and ϕq between 2Rein and 0.5Rein) and velocity dispersion σv. Left: difference in q; right: difference in ϕq. The red dashed line indicates zero difference. 2~6 kpc 6~10 kpc 10~14 kpc 14~18 kpc 18~22 kpc 22~26 kpc 26~30 kpc bins 0.2 0.0 0.2 0.4 qrig ht qleft 2~6 kpc 6~10 kpc 10~14 kpc 14~18 kpc 18~22 kpc 22~26 kpc 26~30 kpc bins 20 10 0 10 20 q, rig ht q, left (d … view at source ↗
Figure 9
Figure 9. Figure 9: Box plot of variations in q and ϕq (value at the right side minus that at the left side) in different radial bins with a width of 4 kpc. Left: difference in q; right: difference in ϕq. For each bin, the red line means median, the box spans from Q1 (first quartile) to Q3 (third quartile), and the whiskers extend from the box to the smallest and largest values within 1.5 × IQR (interquartile range). The dash… view at source ↗
Figure 10
Figure 10. Figure 10: Distributions of best (successful) macro-model parameters, s, q, ϕq, γ, ϕγ − ϕq, for all seven kinds of perturbations each added to EPL+γ to fit to image positions only (assuming σp = 2 mas, essentially the last row in [PITH_FULL_IMAGE:figures/full_fig_p013_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: The left panel presents the distributions of the radial difference ∆q2−6kpc ≡ q6kpc − q2kpc for both the successful and failed cases when adopting EPL+γ in the presence of perturbation Qv alone (assuming σf = 10%, σp = 10 mas), and in the right panel the distributions for ∆ϕq, 2−6kpc ≡ ϕq, 6kpc − ϕq, 2kpc when adding perturbation PAv alone. The samples essentially corresponds to the fourth row in [PITH_F… view at source ↗
Figure 12
Figure 12. Figure 12: Distributions of the best (successful) macro-model parameters, s, q, ϕq, γ, ϕγ −ϕq, for three kinds of perturbations (m3, m4, m3+m4), each added to EPL+γ to fit both image positions and flux ratios (assuming σf = 2% and σp = 10 mas, essentially the last row in [PITH_FULL_IMAGE:figures/full_fig_p016_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: The critical curves and caustics together with the image and source positions predicted by one of the best– fit models in EPL+γ+m3+m4+PAv+Qv (assuming σf = 2% and σp = 10 mas). The best model parameters are also marked out by the vertical lines in [PITH_FULL_IMAGE:figures/full_fig_p017_13.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

133 extracted references · 32 canonical work pages · 1 internal anchor

  1. [1]

    C., Nightingale, J., He, Q., et al

    Amorisco, N. C., Nightingale, J., He, Q., et al. 2022, MNRAS, 510, 2464, doi: 10.1093/mnras/stab3527

    doi:10.1093/mnras/stab3527 2022

  2. [2]

    2023, Nature Astronomy, 7, 736, doi: 10.1038/s41550-023-01943-9 Astropy Collaboration, et al

    Amruth, A., Broadhurst, T., Lim, J., et al. 2023, Nature Astronomy, 7, 736, doi: 10.1038/s41550-023-01943-9 Astropy Collaboration, et al. 2013, A&A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f Astropy Collaboration, et al. 2022, ApJ, 935, 167, doi: 10.3847/1538-4357/ac7c74

    doi:10.1038/s41550-023-01943-9 2023

  3. [3]

    W., Treu, T., Bolton, A

    Auger, M. W., Treu, T., Bolton, A. S., et al. 2010, ApJ, 724, 511, doi: 10.1088/0004-637X/724/1/511

    doi:10.1088/0004-637x/724/1/511 2010

  4. [4]

    1988, A&AS, 74, 385

    Bender, R., Doebereiner, S., & Moellenhoff, C. 1988, A&AS, 74, 385

    1988

  5. [5]

    1989, A&A, 217, 35

    Madejsky, R. 1989, A&A, 217, 35

    1989

  6. [6]

    1999, AJ, 118, 14, doi: 10.1086/300949

    Bernstein, G., & Fischer, P. 1999, AJ, 118, 14, doi: 10.1086/300949

    work page doi:10.1086/300949 1999

  7. [7]

    A., & Rhee, G

    Bernstein, G., Fischer, P., Tyson, J. A., & Rhee, G. 1997, ApJL, 483, L79, doi: 10.1086/310747

    work page doi:10.1086/310747 1997

  8. [8]

    2021, ApJ, 919, 38, doi: 10.3847/1538-4357/ac1108

    Birrer, S. 2021, ApJ, 919, 38, doi: 10.3847/1538-4357/ac1108

    doi:10.3847/1538-4357/ac1108 2021

  9. [9]

    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

    doi:10.1016/j.dark.2018.11.002 2018

  10. [10]

    J., Gilman, D., et al

    Birrer, S., Shajib, A. J., Gilman, D., et al. 2021, Journal of Open Source Software, 6, 3283, doi: 10.21105/joss.03283

    doi:10.21105/joss.03283 2021

  11. [11]

    1986, ApJ, 310, 568, doi: 10.1086/164709 Bradaˇ c, M., Schneider, P., Lombardi, M., et al

    Blandford, R., & Narayan, R. 1986, ApJ, 310, 568, doi: 10.1086/164709 Bradaˇ c, M., Schneider, P., Lombardi, M., et al. 2004, A&A, 423, 797, doi: 10.1051/0004-6361:20040168 Bradaˇ c, M., Schneider, P., Steinmetz, M., et al. 2002, A&A, 388, 373, doi: 10.1051/0004-6361:20020559

    doi:10.1086/164709 1986

  12. [12]

    2025, astropy/photutils: 2.3.0, 2.3.0 Zenodo, doi: 10.5281/zenodo.17129028

    Bradley, L., Sip˝ ocz, B., Robitaille, T., et al. 2025, astropy/photutils: 2.3.0, 2.3.0 Zenodo, doi: 10.5281/zenodo.17129028

    doi:10.5281/zenodo.17129028 2025

  13. [13]

    R., Bolton, A

    Brownstein, J. R., Bolton, A. S., Schlegel, D. J., et al. 2012, ApJ, 744, 41, doi: 10.1088/0004-637X/744/1/41

    doi:10.1088/0004-637x/744/1/41 2012

  14. [14]

    2003, MNRAS, 344, 1000, doi: 10.1046/j.1365-8711.2003.06897.x

    Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000, doi: 10.1046/j.1365-8711.2003.06897.x

    doi:10.1046/j.1365-8711.2003.06897.x 2003

  15. [15]

    W., et al

    Cao, X., Li, R., Nightingale, J. W., et al. 2022, Research in Astronomy and Astrophysics, 22, 025014, doi: 10.1088/1674-4527/ac3f2b

    work page doi:10.1088/1674-4527/ac3f2b 2022

  16. [16]

    2006, MNRAS, 366, 1126, doi: 10.1111/j.1365-2966.2005.09981.x 19

    Cappellari, M., Bacon, R., Bureau, M., et al. 2006, MNRAS, 366, 1126, doi: 10.1111/j.1365-2966.2005.09981.x 19

    doi:10.1111/j.1365-2966.2005.09981.x 2006

  17. [17]

    2020, PhRvL, 125, 111102, doi: 10.1103/PhysRevLett.125.111102

    Broadhurst, T. 2020, PhRvL, 125, 111102, doi: 10.1103/PhysRevLett.125.111102

    doi:10.1103/physrevlett.125.111102 2020

  18. [18]

    V., & Keeton, C

    Chen, J., Kravtsov, A. V., & Keeton, C. R. 2003, ApJ, 592, 24, doi: 10.1086/375639

    work page doi:10.1086/375639 2003

  19. [19]

    2002, ApJ, 565, 17, doi: 10.1086/324493

    Chiba, M. 2002, ApJ, 565, 17, doi: 10.1086/324493

    work page doi:10.1086/324493 2002

  20. [20]

    Inoue, K. T. 2005, ApJ, 627, 53, doi: 10.1086/430403

    doi:10.1086/430403 2005

  21. [21]

    2024, MNRAS, 531, 3431, doi: 10.1093/mnras/stae1228

    Vegetti, S. 2024, MNRAS, 531, 3431, doi: 10.1093/mnras/stae1228

    doi:10.1093/mnras/stae1228 2024

  22. [22]

    B., & Keeton, C

    Congdon, A. B., & Keeton, C. R. 2005, MNRAS, 364, 1459, doi: 10.1111/j.1365-2966.2005.09699.x

    doi:10.1111/j.1365-2966.2005.09699.x 2005

  23. [23]

    Dalal, N., & Kochanek, C. S. 2002, ApJ, 572, 25, doi: 10.1086/340303

    doi:10.1086/340303 2002

  24. [24]

    Despali, G., Vegetti, S., White, S. D. M., Giocoli, C., & van den Bosch, F. C. 2018, MNRAS, 475, 5424, doi: 10.1093/mnras/sty159

    doi:10.1093/mnras/sty159 2018

  25. [25]

    2023, MNRAS, 522, 5434, doi: 10.1093/mnras/stad1313

    Dike, V., Gilman, D., & Treu, T. 2023, MNRAS, 522, 5434, doi: 10.1093/mnras/stad1313

    work page doi:10.1093/mnras/stad1313 2023

  26. [26]

    Dobler, G., & Keeton, C. R. 2006, MNRAS, 365, 1243, doi: 10.1111/j.1365-2966.2005.09809.x

    doi:10.1111/j.1365-2966.2005.09809.x 2006

  27. [27]

    W., Massey, R., et al

    Etherington, A., Nightingale, J. W., Massey, R., et al. 2024, MNRAS, 531, 3684, doi: 10.1093/mnras/stae1375

    doi:10.1093/mnras/stae1375 2024

  28. [29]

    W., & Witt, H

    Evans, N. W., & Witt, H. J. 2003b, MNRAS, 345, 1351, doi: 10.1046/j.1365-2966.2003.07057.x

    doi:10.1046/j.1365-2966.2003.07057.x 2003

  29. [30]

    Fadely, R., & Keeton, C. R. 2012, MNRAS, 419, 936, doi: 10.1111/j.1365-2966.2011.19729.x

    doi:10.1111/j.1365-2966.2011.19729.x 2012

  30. [31]

    R., Nakajima, R., & Bernstein, G

    Fadely, R., Keeton, C. R., Nakajima, R., & Bernstein, G. M. 2010, ApJ, 711, 246, doi: 10.1088/0004-637X/711/1/246

    work page doi:10.1088/0004-637x/711/1/246 2010

  31. [32]

    1989, A&AS, 79, 291

    Fasano, G., & Bonoli, C. 1989, A&AS, 79, 291

    1989

  32. [33]

    Gallazzi, A., Charlot, S., Brinchmann, J., White, S. D. M., & Tremonti, C. A. 2005, MNRAS, 362, 41, doi: 10.1111/j.1365-2966.2005.09321.x

    doi:10.1111/j.1365-2966.2005.09321.x 2005

  33. [34]

    Nierenberg, A. M. 2017, MNRAS, 467, 3970, doi: 10.1093/mnras/stx158

    work page doi:10.1093/mnras/stx158 2017

  34. [35]

    Gilman, D., Birrer, S., Nierenberg, A., & Oh, M. S. H. 2024, MNRAS, 533, 1687, doi: 10.1093/mnras/stae1810

    doi:10.1093/mnras/stae1810 2024

  35. [36]

    2020a, MNRAS, 491, 6077, doi: 10.1093/mnras/stz3480

    Gilman, D., Birrer, S., Nierenberg, A., et al. 2020a, MNRAS, 491, 6077, doi: 10.1093/mnras/stz3480

    doi:10.1093/mnras/stz3480

  36. [37]

    2018, MNRAS, 481, 819, doi: 10.1093/mnras/sty2261

    Nierenberg, A. 2018, MNRAS, 481, 819, doi: 10.1093/mnras/sty2261

    work page doi:10.1093/mnras/sty2261 2018

  37. [38]

    2019, MNRAS, 487, 5721, doi: 10.1093/mnras/stz1593

    Gilman, D., Birrer, S., Treu, T., Nierenberg, A., & Benson, A. 2019, MNRAS, 487, 5721, doi: 10.1093/mnras/stz1593

    doi:10.1093/mnras/stz1593 2019

  38. [39]

    2021, MNRAS, 507, 2432, doi: 10.1093/mnras/stab2335

    Gilman, D., Bovy, J., Treu, T., et al. 2021, MNRAS, 507, 2432, doi: 10.1093/mnras/stab2335

    doi:10.1093/mnras/stab2335 2021

  39. [40]

    2020b, MNRAS, 492, L12, doi: 10.1093/mnrasl/slz173

    Gilman, D., Du, X., Benson, A., et al. 2020b, MNRAS, 492, L12, doi: 10.1093/mnrasl/slz173

    doi:10.1093/mnrasl/slz173

  40. [41]

    2023, PhRvD, 107, 103008, doi: 10.1103/PhysRevD.107.103008

    Gilman, D., Zhong, Y.-M., & Bovy, J. 2023, PhRvD, 107, 103008, doi: 10.1103/PhysRevD.107.103008

    doi:10.1103/physrevd.107.103008 2023

  41. [42]

    M., Treu, T., et al

    Gilman, D., Nierenberg, A. M., Treu, T., et al. 2025, arXiv e-prints, arXiv:2511.07513, doi: 10.48550/arXiv.2511.07513

    Pith/arXiv arXiv doi:10.48550/arxiv.2511.07513 2025

  42. [43]

    N., Mao, S., Deng, Z

    Hao, C. N., Mao, S., Deng, Z. G., Xia, X. Y., & Wu, H. 2006a, MNRAS, 370, 1339, doi: 10.1111/j.1365-2966.2006.10545.x

    doi:10.1111/j.1365-2966.2006.10545.x 2006

  43. [44]

    N., Mao, S., Deng, Z

    Hao, C. N., Mao, S., Deng, Z. G., Xia, X. Y., & Wu, H. 2006b, MNRAS, 373, 1264, doi: 10.1111/j.1365-2966.2006.11137.x

    doi:10.1111/j.1365-2966.2006.11137.x 2006

  44. [45]

    R., Millman, K

    Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357, doi: 10.1038/s41586-020-2649-2

    doi:10.1038/s41586-020-2649-2 2020

  45. [46]

    2019, MNRAS, 491, 4247–4253, doi: 10.1093/mnras/stz3305

    Harvey, D., Valkenburg, W., Tamone, A., et al. 2019, MNRAS, 491, 4247–4253, doi: 10.1093/mnras/stz3305

    work page doi:10.1093/mnras/stz3305 2019

  46. [47]

    2023, MNRAS, 518, 220, doi: 10.1093/mnras/stac2779

    He, Q., Nightingale, J., Robertson, A., et al. 2023, MNRAS, 518, 220, doi: 10.1093/mnras/stac2779

    doi:10.1093/mnras/stac2779 2023

  47. [48]

    W., & Blandford, R

    Hogg, D. W., & Blandford, R. D. 1994, MNRAS, 268, 889, doi: 10.1093/mnras/268.4.889

    work page doi:10.1093/mnras/268.4.889 1994

  48. [49]

    2026, Flux-ratio anomalies in cusp quasars reveal dark matter beyond CDM, doi: 10.48550/arXiv.2601.16818

    Hou, S., Xiang, S., Sming Tsai, Y.-L., et al. 2026, Flux-ratio anomalies in cusp quasars reveal dark matter beyond CDM, doi: 10.48550/arXiv.2601.16818

    doi:10.48550/arxiv.2601.16818 2026

  49. [50]

    2018, MNRAS, 475, 2438, doi: 10.1093/mnras/stx3320

    Hsueh, J.-W., Despali, G., Vegetti, S., et al. 2018, MNRAS, 475, 2438, doi: 10.1093/mnras/stx3320

    work page doi:10.1093/mnras/stx3320 2018

  50. [51]

    W., Enzi, W., Vegetti, S., et al

    Hsueh, J. W., Enzi, W., Vegetti, S., et al. 2020, MNRAS, 492, 3047, doi: 10.1093/mnras/stz3177

    doi:10.1093/mnras/stz3177 2020

  51. [52]

    D., Vegetti, S., et al

    Hsueh, J.-W., Fassnacht, C. D., Vegetti, S., et al. 2016, MNRAS, 463, L51

    2016

  52. [53]

    W., Oldham, L., Spingola, C., et al

    Hsueh, J. W., Oldham, L., Spingola, C., et al. 2017, MNRAS, 469, 3713, doi: 10.1093/mnras/stx1082

    work page doi:10.1093/mnras/stx1082 2017

  53. [54]

    Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90, doi: 10.1109/MCSE.2007.55

    doi:10.1109/mcse.2007.55 2007

  54. [55]

    D., Foltz, C

    Impey, C. D., Foltz, C. B., Petry, C. E., Browne, I. W. A., & Patnaik, A. R. 1996, ApJL, 462, L53, doi: 10.1086/310035

    work page doi:10.1086/310035 1996

  55. [56]

    Inoue, K. T. 2016, MNRAS, 461, 164, doi: 10.1093/mnras/stw1270

    doi:10.1093/mnras/stw1270 2016

  56. [57]

    T., & Takahashi, R

    Inoue, K. T., & Takahashi, R. 2012, MNRAS, 426, 2978, doi: 10.1111/j.1365-2966.2012.21915.x

    doi:10.1111/j.1365-2966.2012.21915.x 2012

  57. [58]

    T., Takahashi, R., Takahashi, T., & Ishiyama, T

    Inoue, K. T., Takahashi, R., Takahashi, T., & Ishiyama, T. 2015, MNRAS, 448, 2704, doi: 10.1093/mnras/stv194

    work page doi:10.1093/mnras/stv194 2015

  58. [59]

    Jackson, N., Xanthopoulos, E., & Browne, I. W. A. 2000, MNRAS, 311, 389, doi: 10.1046/j.1365-8711.2000.03115.x

    doi:10.1046/j.1365-8711.2000.03115.x 2000

  59. [60]

    Jedrzejewski, R. I. 1987, MNRAS, 226, 747, doi: 10.1093/mnras/226.4.747 20

    doi:10.1093/mnras/226.4.747 1987

  60. [61]

    T., & Takahashi, T

    Kamada, A., Inoue, K. T., & Takahashi, T. 2016, PhRvD, 94, 023522, doi: 10.1103/PhysRevD.94.023522

    work page doi:10.1103/physrevd.94.023522 2016

  61. [62]

    E., Nierenberg, A

    Keeley, R. E., Nierenberg, A. M., Gilman, D., et al. 2023, MNRAS, 524, 6159, doi: 10.1093/mnras/stad2251

    doi:10.1093/mnras/stad2251 2023

  62. [63]

    E., Nierenberg, A

    Keeley, R. E., Nierenberg, A. M., Gilman, D., et al. 2024, MNRAS, 535, 1652, doi: 10.1093/mnras/stae2458

    doi:10.1093/mnras/stae2458 2024

  63. [64]

    E., Nierenberg, A

    Keeley, R. E., Nierenberg, A. M., Gilman, D., et al. 2025, arXiv e-prints, arXiv:2511.07765, doi: 10.48550/arXiv.2511.07765

    doi:10.48550/arxiv.2511.07765 2025

  64. [65]

    Keeton, C. R. 2001, arXiv e-prints, astro, doi: 10.48550/arXiv.astro-ph/0102340

    Pith/arXiv arXiv doi:10.48550/arxiv.astro-ph/0102340 2001

  65. [66]

    R., Gaudi, B

    Keeton, C. R., Gaudi, B. S., & Petters, A. O. 2003, ApJ, 598, 138, doi: 10.1086/378934

    doi:10.1086/378934 2003

  66. [67]

    R., Gaudi, B

    Keeton, C. R., Gaudi, B. S., & Petters, A. O. 2005, ApJ, 635, 35, doi: 10.1086/497324

    work page doi:10.1086/497324 2005

  67. [68]

    R., Kochanek, C

    Keeton, C. R., Kochanek, C. S., & Seljak, U. 1997, ApJ, 482, 604, doi: 10.1086/304172

    doi:10.1086/304172 1997

  68. [69]

    R., Falco, E

    Keeton, C. R., Falco, E. E., Impey, C. D., et al. 2000, ApJ, 542, 74, doi: 10.1086/309517

    work page doi:10.1086/309517 2000

  69. [70]

    Kochanek, C. S. 2021, MNRAS, 501, 5021, doi: 10.1093/mnras/staa4033

    work page doi:10.1093/mnras/staa4033 2021

  70. [71]

    S., & Dalal, N

    Kochanek, C. S., & Dalal, N. 2004, ApJ, 610, 69, doi: 10.1086/421436

    doi:10.1086/421436 2004

  71. [72]

    Koopmans, L. V. E., Biggs, A., Blandford, R. D., et al. 2003, ApJ, 595, 712, doi: 10.1086/377434

    work page doi:10.1086/377434 2003

  72. [73]

    1994, A&A, 286, 357, doi: 10.48550/arXiv.astro-ph/9311011

    Kormann, R., Schneider, P., & Bartelmann, M. 1994, A&A, 286, 357, doi: 10.48550/arXiv.astro-ph/9311011

    Pith/arXiv arXiv doi:10.48550/arxiv.astro-ph/9311011 1994

  73. [74]

    B., Cornell, M

    Kormendy, J., Fisher, D. B., Cornell, M. E., & Bender, R. 2009, ApJS, 182, 216, doi: 10.1088/0067-0049/182/1/216 Krajnovi´ c, D., Alatalo, K., Blitz, L., et al. 2013, MNRAS, 432, 1768, doi: 10.1093/mnras/sts315

    doi:10.1088/0067-0049/182/1/216 2009

  74. [75]

    K., Pitrou, A., & Seibert, S

    Lam, S. K., Pitrou, A., & Seibert, S. 2015, in Proceedings of the Second Workshop on the LLVM Compiler Infrastructure in HPC, LLVM ’15 (New York, NY, USA: Association for Computing Machinery), doi: 10.1145/2833157.2833162

    doi:10.1145/2833157.2833162 2015

  75. [76]

    Lauer, T. R. 1985, MNRAS, 216, 429, doi: 10.1093/mnras/216.2.429

    work page doi:10.1093/mnras/216.2.429 1985

  76. [77]

    R., Neugebauer, G., Weir, N., Matthews, K., & Patnaik, A

    Lawrence, C. R., Neugebauer, G., Weir, N., Matthews, K., & Patnaik, A. R. 1992, MNRAS, 259, 5P, doi: 10.1093/mnras/259.1.5P

    work page doi:10.1093/mnras/259.1.5p 1992

  77. [78]

    T., Futamase, T., & Akhlaghi, M

    Lefor, A. T., Futamase, T., & Akhlaghi, M. 2013, NewAR, 57, 1, doi: 10.1016/j.newar.2013.05.001

    work page doi:10.1016/j.newar.2013.05.001 2013

  78. [79]

    2009, MNRAS, 397, 341, doi: 10.1111/j.1365-2966.2009.14912.x

    Liesenborgs, J., de Rijcke, S., Dejonghe, H., & Bekaert, P. 2009, MNRAS, 397, 341, doi: 10.1111/j.1365-2966.2009.14912.x

    doi:10.1111/j.1365-2966.2009.14912.x 2009

  79. [80]

    2014, MNRAS, 437, 2461, doi: 10.1093/mnras/stt2057

    Lubini, M., Sereno, M., Coles, J., Jetzer, P., & Saha, P. 2014, MNRAS, 437, 2461, doi: 10.1093/mnras/stt2057

    work page doi:10.1093/mnras/stt2057 2014

  80. [81]

    1998, MNRAS, 295, 587

    Mao, S., & Schneider, P. 1998, MNRAS, 295, 587

    1998

Showing first 80 references.