pith. sign in

arxiv: 2605.19615 · v1 · pith:MUJBKDU3new · submitted 2026-05-19 · 🌌 astro-ph.CO

Reassessing the Statistical Necessity of Stellar Velocity Anisotropy in Strong-Lensing Cosmology with Lens-by-Lens Photometric Constraints

Jian Hu , Yi Liu , Jian-Ping Hu , Zhongmu Li This is my paper

Pith reviewed 2026-05-20 02:42 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords strong gravitational lensingstellar velocity anisotropysupernova distancesgalaxy density profilesearly-type galaxiescosmological parametersphotometric constraints
0
0 comments X

The pith

Incorporating lens-by-lens photometric constraints shows stellar velocity anisotropy must be treated as a free parameter in strong-lensing cosmology.

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

The paper tests whether the stellar orbital anisotropy parameter can be fixed or must vary freely when strong-lensing systems are used for cosmological distance measurements. It builds a sample of 107 matched strong-lensing and supernova pairs and adds per-lens luminosity density slopes from imaging data. The models then compare cases where anisotropy is fixed to isotropy or a local value against cases where it is left free. A sympathetic reader cares because fixing the parameter produces biased density profiles and larger scatter that would affect inferred cosmic distances and dark-energy constraints.

Core claim

When individual photometric constraints are incorporated, the stellar orbital anisotropy parameter β_ani is statistically required as a free parameter to prevent severe dynamical modeling biases. Fixing β_ani to isotropy (β=0) or a local prior (β=0.18) is strongly disfavored by BIC differences of 14.2 and 48.9 in the P2 framework and similar values in the P3 framework, while also inflating intrinsic scatter. The same analysis detects negative redshift evolution of the lens density slope (γ_z ≈ -0.42 to -0.46) at 1.5–2.0σ across variants.

What carries the argument

The stellar orbital anisotropy parameter β_ani, which quantifies the difference between radial and tangential stellar velocities and degenerates with the lens density slope until broken by per-lens luminosity density slopes from high-resolution imaging.

If this is right

  • Fixing β_ani to isotropy or a constant local value is disfavored by BIC differences of 13–49 and increases model scatter.
  • Early-type galaxy density profiles show negative redshift evolution (γ_z ≈ -0.42 to -0.46) at modest significance.
  • The requirement for a free β_ani holds across both the P2 redshift-only and P3 redshift-plus-surface-density frameworks.
  • Photometric constraints per lens are essential to separate anisotropy from density-profile effects in cosmological applications.

Where Pith is reading between the lines

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

  • Future large lensing surveys would benefit from routinely obtaining high-resolution imaging for every system to allow per-lens anisotropy constraints.
  • Ignoring the need for free anisotropy could systematically shift cosmological parameters inferred from strong lensing relative to other probes.
  • The detected redshift evolution of density slopes could be tested by applying the same photometric approach to independent lens samples at higher redshifts.

Load-bearing premise

The 5% comoving-distance tolerance used to construct the 107 SGL-SN pairs yields accurate distance ratios that are free of significant selection biases or mismatches between the lensing geometry and the supernova distances.

What would settle it

Repeating the full analysis on a sample built with a tighter distance-matching tolerance or with independent integral-field spectroscopic measurements of anisotropy that remove the BIC penalty for fixed β_ani would falsify the necessity of a free parameter.

Figures

Figures reproduced from arXiv: 2605.19615 by Jian Hu, Jian-Ping Hu, Yi Liu, Zhongmu Li.

Figure 1
Figure 1. Figure 1: Overlaid 1D marginalized and 2D joint posterior distributions for the four kinematic scenarios under the primary P2 framework. The systematic shift of the fixed-β models toward lower γ0 and elevated δint is clearly visible, particularly for the fixed-0.18 case (blue contours) [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Posterior distributions for the primary P2 model with βani free under a flat prior. The posterior bulk lies at negative βani, well away from isotropy under the flat-prior assumption [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Posterior distributions for the P2 model with βani under the local Gaussian prior (0.18 ± 0.13). The posterior shifts to βani ≈ −0.15, illustrating the prior sensitivity [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Posterior distributions for the P2 model with βani fixed to isotropy (β = 0) [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Posterior distributions for the P2 model with βani fixed to 0.18. The sub-isothermal central density slope (γ0 ≈ 1.86) and elevated δint reflect the compensation required to absorb the kinematic mismatch [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Posterior distributions for the complementary P3 model with βani free under a flat prior, confirming the primary P2 conclusions within a more complete structural-evolutionary framework [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Posterior distributions for the P3 model with βani under the local Gaussian prior (0.18 ± 0.13) [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Posterior distributions for the P3 model with βani fixed to isotropy (β = 0) [PITH_FULL_IMAGE:figures/full_fig_p016_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Posterior distributions for the P3 model with βani fixed to 0.18. The sub-isothermal γ0 ≈ 1.70 and the negative γs highlight the parameter tension introduced by imposing the local kinematic prior on high-redshift lenses [PITH_FULL_IMAGE:figures/full_fig_p017_9.png] view at source ↗
read the original abstract

The stellar orbital anisotropy parameter ($\beta_{\rm ani}$) is a persistent systematic uncertainty in galaxy-scale strong gravitational lensing (SGL) cosmology. Typically fixed to isotropy or a local prior, it frequently degenerates with the lens density profile. We demonstrate this apparent redundancy largely arises from incomplete photometric constraints. We cross-matched 130 SGL systems with the Pantheon+ SN~Ia compilation, constructing a strictly matched sample of 107 SGL-SN pairs using a 5\% comoving-distance tolerance. Assuming a flat universe ($\Omega_k = 0$), the distance ratio is derived from apparent magnitude differences between paired SNe~Ia, eliminating $H_0$ and absolute magnitude dependence without fitting explicit dark-energy models. To break the kinematic degeneracy, we incorporate lens-by-lens luminosity density slopes ($\delta_i$) from high-resolution imaging. Adopting the quasi-model-independent P2 redshift-evolutionary framework ($\gamma(z) = \gamma_0 + \gamma_z z$), we find very strong statistical evidence for a free $\beta_{\rm ani}$. Fixing $\beta_{\rm ani}$ to isotropy ($\beta=0$) or a local prior ($\beta=0.18$) is strongly disfavored ($\Delta\bic = 14.2$ and $48.9$) and artificially inflates intrinsic scatter. A complementary P3 framework ($\gamma(z,\tilde{\Sigma}) = \gamma_0 + \gamma_z z + \gamma_s \log_{10}\tilde{\Sigma}$) confirms these penalties ($\Delta\bic = 13.5$ and $49.1$). Across all P2 variants, we consistently detect a negative redshift evolution of the density slope ($\gamma_z \approx -0.42$ to $-0.46$; ${\sim}1.5{-}2.0\sigma$), indicating ETG density profiles become shallower at higher redshifts. We conclude that when individual photometric constraints are incorporated, $\beta_{\rm ani}$ is statistically required as a free parameter to prevent severe dynamical modeling biases.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript claims that cross-matching 130 SGL systems with Pantheon+ SNe Ia under a 5% comoving-distance tolerance yields 107 pairs whose magnitude differences provide distance ratios (assuming flatness) that, when combined with lens-by-lens δ_i photometric constraints, demonstrate β_ani must be treated as free in the P2 (γ(z) = γ0 + γz z) and P3 frameworks; fixing it to 0 or 0.18 produces ΔBIC = 14.2/48.9 (P2) and 13.5/49.1 (P3), inflates scatter, and biases the inferred negative γz ≈ −0.42 to −0.46.

Significance. If the paired distance ratios faithfully represent the lensing geometry without selection bias, the result would establish that individual photometric constraints break the β_ani–γ degeneracy and that prior fixed-β_ani analyses incur systematic errors; the consistent γz detection across variants would also support evolving ETG density profiles at the ~1.5–2σ level.

major comments (2)
  1. [§3] §3 (Sample Construction): The 5% comoving-distance tolerance used to build the 107 SGL-SN pairs permits redshift mismatches Δz ~ 0.05–0.1 at typical lens redshifts; without explicit tests for resulting bias in the derived D ratios or robustness checks against the tolerance, the central ΔBIC values (14.2, 48.9) that declare β_ani statistically required cannot be considered secure.
  2. [§5] §5 (Results and BIC Analysis): The abstract and reported BIC differences omit details on error propagation, covariance between magnitude-derived ratios and δ_i, and sensitivity to the 5% matching; these omissions are load-bearing because they directly affect whether the penalties for fixed β_ani remain significant after realistic uncertainties are included.
minor comments (2)
  1. [Abstract] Abstract: The phrase 'very strong statistical evidence' should be accompanied by the conventional BIC interpretation thresholds used (e.g., ΔBIC > 10).
  2. [§2] Notation: The definition of the normalized surface density Σ̃ in the P3 framework should be stated explicitly when first introduced.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful and constructive review of our manuscript. We have addressed each major comment below with additional analyses and clarifications. Revisions have been made to strengthen the robustness checks and statistical details.

read point-by-point responses

  1. Referee: [§3] §3 (Sample Construction): The 5% comoving-distance tolerance used to build the 107 SGL-SN pairs permits redshift mismatches Δz ~ 0.05–0.1 at typical lens redshifts; without explicit tests for resulting bias in the derived D ratios or robustness checks against the tolerance, the central ΔBIC values (14.2, 48.9) that declare β_ani statistically required cannot be considered secure.

    Authors: We appreciate the referee's concern regarding potential bias from redshift mismatches in the sample construction. In the revised manuscript, we have added explicit robustness tests in a new subsection of §3. We reconstructed the paired samples using alternative tolerances of 3%, 4%, 6%, and 7%, and recomputed the full P2 and P3 analyses. Across these variants, the ΔBIC penalties for fixed β_ani remain significant (minimum ΔBIC > 12 for isotropy and > 40 for the local prior), and the recovered γ_z values stay consistent within 1σ (ranging from -0.40 to -0.48). We also verified that magnitude-derived distance ratios exhibit no systematic trend with Δz, confirming negligible bias in the central results. These checks are documented with tables and figures in the revision. revision: yes

  2. Referee: [§5] §5 (Results and BIC Analysis): The abstract and reported BIC differences omit details on error propagation, covariance between magnitude-derived ratios and δ_i, and sensitivity to the 5% matching; these omissions are load-bearing because they directly affect whether the penalties for fixed β_ani remain significant after realistic uncertainties are included.

    Authors: We agree that greater transparency on the statistical pipeline is warranted. The revised manuscript now includes an expanded description in §5 and the methods section detailing error propagation: distance-ratio uncertainties are obtained via standard Gaussian propagation from Pantheon+ magnitude errors, with the flat-universe assumption removing H0 and M dependence. The joint likelihood explicitly incorporates the covariance matrix between the distance ratios and lens-specific δ_i constraints during MCMC sampling. Sensitivity to the 5% tolerance is addressed via the robustness tests described in our response to the §3 comment, which show the BIC differences and γ_z inference are stable. These additions ensure the statistical significance is presented with full context. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is data-driven and self-contained

full rationale

The paper selects 107 SGL-SN pairs via a 5% comoving-distance tolerance on Pantheon+ data, derives distance ratios from apparent magnitude differences under the explicit flat-universe assumption (Ω_k=0), and then performs BIC model comparisons on the resulting distance-ratio constraints while allowing or fixing β_ani. This is a standard statistical test of parameter necessity given external data; the evolutionary parameters γ(z) and β_ani are fitted quantities, not defined in terms of each other or renamed from the input ratios. No self-definitional loop, fitted-input-as-prediction, or load-bearing self-citation chain appears in the abstract or described chain. The 5% tolerance is an acknowledged selection cut whose validity is an external assumption, not a circular reduction.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The central claim depends on fitted evolutionary parameters in the P2/P3 frameworks and the flat-universe assumption used to convert magnitude differences into distance ratios; no new entities are postulated.

free parameters (3)
  • γ0
    Baseline value of the density slope in the redshift-evolutionary model γ(z) = γ0 + γz z
  • γz
    Coefficient describing redshift evolution of the density slope
  • β_ani
    Stellar orbital anisotropy parameter treated as free in the preferred models
axioms (2)
  • domain assumption The universe is spatially flat (Ω_k = 0)
    Required to derive distance ratios from SN Ia apparent magnitude differences without explicit dark-energy models
  • domain assumption The 5% comoving-distance tolerance produces unbiased SGL-SN pairs
    Used to select the final sample of 107 systems from 130 candidates

pith-pipeline@v0.9.0 · 5928 in / 1519 out tokens · 53731 ms · 2026-05-20T02:42:59.729774+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We cross-matched 130 SGL systems with the Pantheon+ SN Ia compilation, constructing a strictly matched sample of 107 SGL-SN pairs using a 5% comoving-distance tolerance... Adopting the quasi-model-independent P2 redshift-evolutionary framework (γ(z) = γ0 + γz z)

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

46 extracted references · 46 canonical work pages · 2 internal anchors

  1. [1]

    F., Aboubrahim, A., et al

    Abdalla, E., Abell´ an, G. F., Aboubrahim, A., et al. 2022, JHEAp, 34, 49

    work page 2022

  2. [2]

    LSST Science Book, Version 2.0

    Abell, P. A., Allison, J., Anderson, S. F., et al. 2009, arXiv:0912.0201

    work page internal anchor Pith review Pith/arXiv arXiv 2009

  3. [3]

    W., & Evans, N

    Agnello, A., Auger, M. W., & Evans, N. W. 2013, MNRAS, 429, L35

    work page 2013

  4. [4]

    W., Treu, T., Bolton, A

    Auger, M. W., Treu, T., Bolton, A. S., et al. 2010, ApJ, 724, 511

    work page 2010

  5. [5]

    2006, ApJ, 640, 514

    Biesiada, M., Godlowski, W., & Szydowski, M. 2006, ApJ, 640, 514

    work page 2006

  6. [6]

    2010, MNRAS, 406, 1055

    Biesiada, M., Pi´ orkowska, A., & Malec, B. 2010, MNRAS, 406, 1055

    work page 2010

  7. [7]

    1980, MNRAS, 190, 873

    Binney, J. 1980, MNRAS, 190, 873

    work page 1980

  8. [8]

    2016, JCAP, 08, 020

    Birrer, S., Amara, A., & Refsdal, S. 2016, JCAP, 08, 020

    work page 2016

  9. [9]

    J., Galan, A., et al

    Birrer, S., Shajib, A. J., Galan, A., et al. 2020, A&A, 643, A165

    work page 2020

  10. [10]

    S., Brownstein, J

    Bolton, A. S., Brownstein, J. R., Kochanek, C. S., et al. 2012, ApJ, 757, 82

    work page 2012

  11. [11]

    2022, ApJ, 938, 110

    Brout, D., Scolnic, D., Popovic, B., et al. 2022, ApJ, 938, 110

    work page 2022

  12. [12]

    2012, JCAP, 03, 016

    Cao, S., Pan, Y., Biesiada, M., et al. 2012, JCAP, 03, 016

    work page 2012

  13. [13]

    2015, ApJ, 806, 185

    Cao, S., Biesiada, M., Gavazzi, R., et al. 2015, ApJ, 806, 185

    work page 2015

  14. [14]

    2019, MNRAS, 488, 3745

    Chen, Y., Li, R., Shu, Y., & Cao, X. 2019, MNRAS, 488, 3745

    work page 2019

  15. [15]

    Collett, T. E. 2015, ApJ, 811, 20 Di Valentino, E., Mena, O., Pan, S., et al. 2021, CQGra, 38, 153001

    work page 2015

  16. [16]

    E., Gorenstein, M., & Shapiro, I

    Falco, E. E., Gorenstein, M., & Shapiro, I. I. 1985, ApJL, 289, L1

    work page 1985

  17. [17]

    W., Lang, D., & Goodman, J

    Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306

    work page 2013

  18. [18]

    2025, A&A, 694, A196

    Geng, S., Grespan, M., Thuruthipilly, H., et al. 2025, A&A, 694, A196

    work page 2025

  19. [19]

    2026, arXiv:2603.04279

    Geng, S., Cao, S., Biesiada, M., et al. 2026, arXiv:2603.04279

    work page arXiv 2026

  20. [20]

    2019, ApJ, 883, 203 Statistical Necessity ofβ ani in SGL Cosmology7

    Gong, Y., Liu, X., Cao, Y., et al. 2019, ApJ, 883, 203 Statistical Necessity ofβ ani in SGL Cosmology7

    work page 2019

  21. [21]

    2026, arXiv:2603.12790

    He, Z., Deng, L., Chen, Q., et al. 2026, arXiv:2603.12790

    work page arXiv 2026

  22. [22]

    Holanda, R. F. L., Lima, J. A. S., & Ribeiro, M. B. 2010, ApJL, 722, L233

    work page 2010

  23. [23]

    2023, PhRvD, 108, 083024

    Hu, J., et al. 2023, PhRvD, 108, 083024

    work page 2023

  24. [24]

    Koopmans, L. V. E., Treu, T., Bolton, A. S., et al. 2006, ApJ, 649, 599

    work page 2006

  25. [25]

    Koopmans, L. V. E., Bolton, A., Treu, T., et al. 2009, ApJL, 703, L51 Krajnovi´ c, D., Emsellem, E., Cappellari, M., et al. 2011, MNRAS, 414, 2923

    work page 2009

  26. [26]

    Euclid Definition Study Report

    Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, arXiv:1110.3193

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  27. [27]

    E., Krawczyk, C

    Li, T., Collett, T. E., Krawczyk, C. M., & Enzi, W. 2024, MNRAS, 527, 5311

    work page 2024

  28. [28]

    2016, ApJ, 822, 74

    Liao, K., Li, Z., Wang, G.-J., & Fan, X.-L. 2016, ApJ, 822, 74

    work page 2016

  29. [29]

    A., & Bou´ e, G

    Mamon, G. A., & Bou´ e, G. 2010, MNRAS, 401, 2433

    work page 2010

  30. [30]

    2020, A&A, 639, A101

    Millon, M., Galan, A., Courbin, F., et al. 2020, A&A, 639, A101

    work page 2020

  31. [31]

    H., & Ostriker, J

    Naab, T., Johansson, P. H., & Ostriker, J. P. 2009, ApJL, 699, L178

    work page 2009

  32. [32]

    F., Frenk, C

    Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493

    work page 1997

  33. [33]

    2010, ApJ, 725, 2312

    Burkert, A. 2010, ApJ, 725, 2312

    work page 2010

  34. [34]

    2022, NewAR, 95, 101659 Planck Collaboration, Aghanim, N., Akrami, Y., et al

    Perivolaropoulos, L., & Skara, F. 2022, NewAR, 95, 101659 Planck Collaboration, Aghanim, N., Akrami, Y., et al. 2020, A&A, 641, A6

    work page 2022

  35. [35]

    Z., Zhao, J

    Qi, J. Z., Zhao, J. W., Cao, S., et al. 2022, MNRAS, 514, 495

    work page 2022

  36. [36]

    G., Yuan, W., Macri, L

    Riess, A. G., Yuan, W., Macri, L. M., et al. 2022, ApJL, 934, L7

    work page 2022

  37. [37]

    J., Gavazzi, R., Marshall, P

    Ruff, A. J., Gavazzi, R., Marshall, P. J., et al. 2011, ApJ, 727, 96

    work page 2011

  38. [38]

    2013, A&A, 559, A37

    Schneider, P., & Sluse, D. 2013, A&A, 559, A37

    work page 2013

  39. [39]

    2013, ApJ, 777, 98

    Sonnenfeld, A., Treu, T., Gavazzi, R., et al. 2013, ApJ, 777, 98

    work page 2013

  40. [40]

    H., Bonvin, V., Courbin, F., et al

    Suyu, S. H., Bonvin, V., Courbin, F., et al. 2017, MNRAS, 468, 2590

    work page 2017

  41. [41]

    P., Bender, R., et al

    Thomas, J., Saglia, R. P., Bender, R., et al. 2009, ApJ, 691, 770

    work page 2009

  42. [42]

    R., Saglia, R

    Tortora, C., Napolitano, N. R., Saglia, R. P., et al. 2014, MNRAS, 445, 115

    work page 2014

  43. [43]

    Treu, T., & Koopmans, L. V. E. 2002, ApJ, 575, 87

    work page 2002

  44. [44]

    2020, ApJ, 904, 153

    Wang, B.-Z., Li, R., Liao, K., et al. 2020, ApJ, 904, 153

    work page 2020

  45. [45]

    C., Suyu, S

    Wong, K. C., Suyu, S. H., Chen, G. C. F., et al. 2020, MNRAS, 498, 1420

    work page 2020

  46. [46]

    decisive evidence

    Wu, H., Chen, Y., Liu, T., et al. 2025, arXiv:2511.08030 8Hu et al. APPENDIX A.ROBUSTNESS TESTS AND MATHEMATICAL FRAMEWORK A.1. EXPLICIT LIKELIHOOD CONSTRUCTION AND JACOBIAN ERROR PROPAGATION The theoretical velocity dispersion is governed by the structural functionF(γ, δ i, βani). To decouple the total mass density slope (γ) from the lens-by-lens luminos...

    work page arXiv 2025