Pith. sign in

REVIEW 3 major objections 5 minor 1 references

This paper claims that the total mass-density slope of massive early-type lens galaxies evolves with both redshift and surface mass density, measured for the first time from a sample of lensed quasars.

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 · deepseek-v4-flash

2026-08-02 23:29 UTC pith:XJ7PIUFI

load-bearing objection Useful new data and a genuine first application, but the headline γs>0 result is not robust because the covariate and the outcome come from the same measured velocity dispersion. the 3 major comments →

arxiv 2602.13578 v2 pith:XJ7PIUFI submitted 2026-02-14 astro-ph.GA

New Dynamical Measurements from a Lensed Quasar Sample: Joint Analysis Constrains the Mass Profile Evolution of Lens Galaxies

Ziyu Guo , Yun Chen , Yiping Shu , Jiaze Gao , Hui Li , Zizhao He , Jun Wang This is my paper
classification astro-ph.GA
keywords strong gravitational lensinglensed quasarsgalaxy mass density profilesstellar velocity dispersionearly-type galaxiesgalaxy evolutionjoint lensing dynamicspower-law slope
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.

The paper tries to establish that the internal mass structure of massive early-type galaxies changes with cosmic time and with galaxy density. By combining strong-lensing geometry with stellar kinematics for 24 lensed quasars, it measures the power-law slope of the total mass profile and finds it decreases with redshift and increases with surface mass density. This matters because it tests galaxy formation models and supports the inside-out growth picture, while validating lensed quasars as an independent probe of galaxy structure.

Core claim

Modeling each lens galaxy's total mass profile as a power law, ρ ∝ r^(−γ), the authors parameterize the slope as γ = γ0 + γz·z_l + γs·log Σ̃, where z_l is the lens redshift and Σ̃ is a normalized surface mass density. From joint lensing–dynamics analysis of 24 systems, they infer γ0 = 1.62(+0.11/−0.12), γz = −0.35(+0.08/−0.09), and γs = +0.37(+0.08/−0.07). The negative γz says that at fixed density, higher-redshift galaxies have shallower (less centrally concentrated) mass profiles; the positive γs says that at fixed redshift, denser galaxies have steeper profiles. The paper argues these trends are statistically significant and consistent with earlier galaxy–galaxy lensing results, and it pr

What carries the argument

The central identity is the joint lensing–dynamics mass equality, M_grl = M_dyn, which equates the projected mass inside the Einstein radius from lens geometry with the dynamical mass from stellar kinematics. The analysis combines power-law density and luminosity profiles, the spherical Jeans equation, and an aperture correction to predict the velocity dispersion within a standardized aperture. The other key piece is a new iterative spectral-fitting procedure that masks contaminating quasar emission lines, allowing velocity dispersions to be measured for lens galaxies whose spectra are otherwise overwhelmed by the background quasar.

Load-bearing premise

The 24 lens systems with measurable velocity dispersions are treated as representative of the 106-system parent population, even though selection depends on the availability and quality of archival spectra; if those systems are biased toward brighter or more massive lenses, the inferred redshift and density trends could be wrong.

What would settle it

Measure stellar velocity dispersions for all 106 parent-sample lenses from new deep spectroscopy and rerun the joint analysis; if the recovered γz and γs shift outside the quoted error bars, the reported trends are an artifact of which systems happened to have archival spectra.

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

If this is right

  • If the measured trends are correct, massive early-type galaxies were systematically less centrally concentrated at higher redshift, supporting inside-out assembly.
  • Denser galaxies having steeper profiles implies a tight coupling between stellar surface density and total mass distribution, useful for galaxy formation models.
  • Lensed quasars become a validated, independent population for joint lensing–dynamics studies, complementing galaxy–galaxy lenses.
  • The measured γ(z, Σ) relation can serve as a prior for time-delay cosmography, reducing a key systematic in Hubble constant measurements.
  • The spectroscopic technique for handling quasar contamination extends directly to the much larger lensed-quasar samples expected from upcoming wide-field surveys.

Where Pith is reading between the lines

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

  • The authors do not model the selection of the 24 systems; because only systems with usable spectra enter the analysis, the reported γz and γs could be biased if measurable lenses are brighter or more massive. A complete kinematic census of the parent sample would settle this.
  • If the γs > 0 relation holds generally, it suggests a unified scaling between galaxy mass profile and stellar surface density that could be applied to non-lensed galaxies as a dynamical mass estimator.
  • The steeper γz compared with galaxy–galaxy compilations may reflect the higher redshift range of lensed quasars, but it could also indicate sample selection; targeted spectroscopy of the remaining parent lenses would test this directly.
  • The parameterized γ(z, Σ) provides a quantitative prediction that cosmological simulations of massive early-type galaxies should reproduce if the trends are physical.

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

3 major / 5 minor

Summary. The paper compiles a parent sample of 106 galaxy-scale strong gravitational lenses with quasar sources and selects a 24-system subset for joint lensing and stellar-kinematics analysis. Newly presented data include 11 stellar velocity dispersions measured from archival BOSS/DESI spectra with an iterative pPXF-based masking procedure, and 4 effective radii from DESI-LS imaging. Assuming a power-law total mass profile, the authors parameterize the slope as γ = γ0 + γz z_l + γs log Σ̃ and report γ0=1.62^{+0.11}_{-0.12}, γz=-0.35^{+0.08}_{-0.09}, γs=0.37^{+0.08}_{-0.07} (68%), interpreting γz<0 and γs>0 as robust evidence for redshift and surface-density trends consistent with earlier galaxy-galaxy strong-lensing studies.

Significance. If the central result holds, the paper would provide a valuable independent confirmation of mass-profile evolution trends using lensed quasars, a sample that is complementary to galaxy-galaxy lenses and that extends to higher lens redshifts. The new spectroscopic measurements for 11 lens galaxies and effective radii for 4 systems are useful observational contributions, and the joint lensing-dynamics machinery is applied in a relatively clean Bayesian framework. However, the headline claims of 'robustly demonstrate' are currently undermined by a self-referential covariate in the γ parameterization and by an incompletely characterized selection from the parent sample. The paper's contribution is therefore promising but not yet established.

major comments (3)
  1. [§4.1, Eqs. (10)–(11); §3.5, Eq. (9)] The surface-density covariate Σ̃ in Eq. (11) is constructed from σ_e2, the aperture-corrected observed velocity dispersion defined in Eq. (8). This same σ_obs,e2 is the dependent variable in the Gaussian likelihood, Eq. (9). Through Eq. (10), the model slope γ_i becomes a function of the noisy observation. A positive measurement error in σ_obs inflates Σ̃, increases γ_i, and through Eq. (6) raises the predicted σ_th, moving it toward the observed value; a negative error acts oppositely. The posterior can therefore absorb measurement noise, producing a spuriously positive γ_s. The reported uncertainty on γ_s (~0.08) reflects only statistical scatter under the assumed model, not this feedback. This is an endogenous-regressor problem. I request a mock-injection test: generate synthetic samples with known γ_s and realistic σ_obs errors, then check recovery; and/or a reanalysis where Σ̃ is ba
  2. [§2.1–2.2] The sample selection from 106 to 24 systems is based on data availability, but the selection process is not fully described and is not modeled. Section 2.1 states that 13 systems had all required values and that 11 new σ_obs and 4 new θ_eff measurements expanded the sample to 24. Section 2.2 then reports that a cross-match with BOSS/DESI yielded 70 matches, three spectra were discarded as lacking reliable absorption features, and 'this effort produced new velocity dispersion measurements for 11 systems'. This leaves about 56 matches unused without explanation. If the 11 were selected by S/N, visual quality, or other criteria, the analysis sample is not a random subset of the parent population. Differences in redshift, Einstein radius, velocity dispersion, or effective radius between the included and excluded systems could bias both γ_z and γ_s. At minimum, compare the distributions of th
  3. [§4.1 and §3.4–3.5] The statistical model is incompletely specified. The parameter vector is stated as Θ = {γ0, γz, γs, Ωm}, but the likelihood in Eq. (9) uses a total uncertainty Δσ_tot that 'incorporates the measurement error of σ_ap, the uncertainty propagated from the aperture correction (η), and a systematic term'. The prior on η is given in §3.4, but it is not stated whether η is sampled, analytically marginalized, or included as an added variance term. The 'systematic term' is never defined, and no prior or fixed value is provided for it. This prevents exact reproduction of the likelihood and leaves the amplitude of the systematic uncertainty unquantified. Please spell out the full generative model, including how η and σ_sys enter, and the priors used.
minor comments (5)
  1. [Section 2.2] Please clarify the count discrepancy: 70 cross-matches, 3 discarded, but only 11 velocity dispersions reported. The selection among the remaining spectra is a key part of sample construction and should be explicit.
  2. [Table 1 / Appendix A] There are naming inconsistencies: 'SDSSJ1640+1932' in Table 1 appears as 'SDSSJ1640-1932' in Appendix A; 'WFI2033-4723' in Table 1 appears as 'WFJ2033-4723' in Appendix A. Please unify.
  3. [§3.4 / Table 1] For slit-based observations, the table lists a single θ_ap; please specify whether an equivalent circular aperture radius was adopted for the correction in Eq. (8), and how this was derived.
  4. [Figure 5] The caption states that dashed lines show 'the mean value'. For skewed posteriors, the posterior median or mode would be more standard; please clarify which statistic is plotted.
  5. [Abstract / §4.2] The phrase 'robustly demonstrate' is stronger than the current evidence warrants given the issues above. I suggest tempering the wording unless the mock-injection test and selection checks support the claim.

Circularity Check

1 steps flagged

γ_s>0 is partly an artifact of using the observed velocity dispersion as both the dependent variable and the surface-density covariate (Eq. 11 vs Eq. 9).

specific steps
  1. self definitional [Section 4.1, Eq. (11), with Eq. (8), Eq. (9), and Eq. (6)]
    "Σ̃ = (σe2/100 km s−1)^2 / (Reff/10 h−1 kpc), where σe2 is the aperture-corrected velocity dispersion taken from Equation (8). ... χ2 = Σ_i ((σth,i − σobs e2,i)/Δσtot,i)^2."

    The covariate defining γ_i (Eq. 10–11) is built from the same observed σ_e2 that is the dependent variable in the Gaussian likelihood (Eq. 9). Since the model prediction σ_th (Eq. 6) depends on γ_i, it is therefore a function of the observed σ_e2. A positive measurement error in σ_e2 inflates Σ̃, raises γ_i, and moves σ_th toward the observed value; a negative error does the opposite. The posterior can absorb this noise by preferring a spuriously positive γ_s. The paper reports γ_s = 0.37 ± 0.08 and declares 'γ_s > 0' robust, but the analysis regresses σ on a function of itself and includes no mock-injection or covariance test to show the endogeneity is negligible. Thus the headline surface-density dependence reduces partly by construction.

full rationale

The central claim that γ increases with surface mass density (γ_s > 0) is compromised by a genuine self-referential step: the paper defines the surface-density covariate using the aperture-corrected velocity dispersion (Eq. 11), and the same aperture-corrected velocity dispersion is the quantity predicted in the likelihood (Eq. 9). This is the textbook endogenous-regressor problem, and it directly targets the paper's primary new quantitative result. The redshift trend γ_z < 0 is less affected because z_l is an independent covariate, and the new velocity-dispersion measurements themselves are a useful observational contribution. The self-citations to Chen et al. (2019) for the δ and β priors are not load-bearing for the headline claim, and the parameterization following Chen et al. (2019) is a standard external framework, not a uniqueness theorem. The selection-effect concern is real but secondary. Overall, the surface-density dependence is not independently established by the quoted equations; a mock-injection or explicit covariance treatment is needed before 'robustly demonstrate γ_s > 0' is supportable. This warrants a partial-circularity score of 6 rather than a higher score because the redshift evolution and the data acquisition retain independent content.

Axiom & Free-Parameter Ledger

8 free parameters · 8 axioms · 0 invented entities

No new physical entities are introduced. The analysis uses standard lensing and dynamical equations, with fitted slope coefficients and external priors on nuisance parameters. The main free parameters are the three slope coefficients; δ, β, Ωm, and η are sampled/marginalized with priors, and an unquantified systematic term is included in the uncertainty.

free parameters (8)
  • γ0 = 1.62^{+0.11}_{-0.12}
    Intercept of the power-law slope relation (Eq 10), fitted to the 24-lens sample.
  • γz = -0.35^{+0.08}_{-0.09}
    Redshift-evolution coefficient in Eq 10.
  • γs = 0.37^{+0.08}_{-0.07}
    Surface-density coefficient in Eq 10.
  • δ (stellar luminosity slope) = 2.173±0.085 (prior)
    Global nuisance parameter, marginalized with a Gaussian prior from Chen et al. (2019).
  • β (orbital anisotropy) = 0.18±0.13 (prior)
    Global nuisance parameter, marginalized with a Gaussian prior from Chen et al. (2019).
  • Ωm = 0.2975±0.0086 (prior)
    Matter density parameter, sampled with a tight Gaussian prior from DESI BAO; enters distance calculations.
  • η (aperture-correction exponent) = -0.066±0.035 (prior)
    Adopted from Cappellari et al. (2006)/Chen et al. (2019) for Eq 8.
  • σsys (systematic uncertainty term)
    Added in quadrature to Δσtot in Eq 9, but its value and derivation are not specified in the text.
axioms (8)
  • domain assumption Lens galaxies are spherically symmetric and in dynamical equilibrium (Section 3.3).
    Used to derive the Jeans equation and equate lensing and dynamical mass.
  • domain assumption Total mass density and stellar luminosity density follow power laws, ρ∝r^{-γ}, ν∝r^{-δ} (Eq 4).
    Underlies the dynamical mass and velocity-dispersion predictions.
  • domain assumption Stellar orbital anisotropy β is constant with radius (Section 3.3(iii)).
    Simplifies the Jeans solution; uncertainty is marginalized with a prior.
  • standard math General-relativistic equivalence of lensing mass and dynamical mass within the Einstein radius (Eq 2).
    Standard result in lensing; not an ad hoc assumption.
  • domain assumption Flat ΛCDM expansion with Ωm prior from DESI BAO (Eqs 12-14).
    Cosmological distances depend on Ωm; the prior is external and tight.
  • standard math The velocity-dispersion prediction from the spherical Jeans equation (Eq 6).
    Standard derivation; assumed to hold for the lens galaxies.
  • domain assumption Aperture correction of observed velocity dispersion with exponent η from Cappellari et al. (2006) (Eq 8).
    Adopted prior; applies to all galaxies regardless of slit geometry.
  • domain assumption The Gaussian priors on δ and β from Chen et al. (2019) are applicable to this sample.
    These priors originate from earlier work by overlapping authors; their applicability to a lensed-quasar sample is assumed.

pith-pipeline@v1.3.0-alltime-deepseek · 20345 in / 17335 out tokens · 149881 ms · 2026-08-02T23:29:32.194830+00:00 · methodology

0 comments
read the original abstract

We present a systematic study of the internal mass structure of early-type galaxies (ETGs) based on 106 galaxy-scale strong gravitational lenses with background quasars, all having spectroscopic redshifts. From this parent sample, we select 24 systems with high-quality ancillary data for joint analysis of strong lensing geometry and stellar kinematics. A key contribution is the derivation of new single-aperture stellar velocity dispersions for 11 lens galaxies via an iterative spectroscopic fitting procedure that mitigates quasar contamination, providing previously unavailable data. We model the total mass-density profile as a power law, $\rho \propto r^{-\gamma}$, and parameterise its logarithmic slope as $\gamma = \gamma_0 + \gamma_z \cdot z_l + \gamma_s \cdot \log \tilde{\Sigma}$, where $z_l$ is the lens redshift and $\tilde{\Sigma}$ the surface mass density. Within a flat $\Lambda$CDM framework and using DESI BAO measurements as a prior, we constrain the parameters via Monte Carlo nested sampling to $\gamma_0 = 1.62^{+0.11}_{-0.12}$, $\gamma_z = -0.35^{+0.08}_{-0.09}$, and $\gamma_s = 0.37^{+0.08}_{-0.07}$ ($68\%$ confidence intervals). Our results robustly demonstrate that $\gamma$ increases with surface mass density ($\gamma_s > 0$) and decreases with redshift ($\gamma_z < 0$). This implies that, at fixed redshift, galaxies with denser stellar cores have steeper mass profiles, while at fixed density, profiles become shallower at higher redshifts. By successfully applying the joint lensing--dynamics method to a substantial, independently acquired sample of lensed quasars, this work provides crucial validation of structural trends previously observed in galaxy--galaxy lensing systems, reinforcing the established evolutionary picture for massive ETGs and establishing lensed quasars as a potent probe of galaxy structure.

Figures

Figures reproduced from arXiv: 2602.13578 by Hui Li, Jiaze Gao, Jun Wang, Yiping Shu, Yun Chen, Ziyu Guo, Zizhao He.

Figure 1
Figure 1. Figure 1: Redshift distributions of the lens and source samples. Left panel: distribution of lens galaxy redshifts (𝑧𝑙). Right panel: distribution of background source (quasar) redshifts (𝑧𝑠). The black bars represent the Parent Sample of 106 systems, and the black shaded region represents the Analysis Sample of 24 systems used for the joint analysis [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Stellar velocity dispersion measurements for the four lensing galaxies, derived from SDSS-III BOSS spectra. Each system is displayed in two panels. Upper panel: The black line shows the extracted and smoothed one-dimensional spectrum of the galaxy, with the wavelength range adjusted to emphasise absorption features. The red line indicates the best-fitting model obtained from the pPXF software applied to th… view at source ↗
Figure 3
Figure 3. Figure 3: Stellar velocity dispersion measurements for the seven lensing galaxies, derived from DESI DR1 spectra. The layout follows that of [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Half-light radius fitting results. Each row corresponds to a lens system. Columns show (from left to right): (1) the original DESI-LS image; (2) the model for the lensed quasar images; (3) the quasar-light-subtracted data used for fitting the galaxy light; (4) the best-fitting de Vaucouleurs galaxy light model; and (5) the normalised residual map. The derived effective radius is annotated. on cosmological … view at source ↗
Figure 5
Figure 5. Figure 5: Posterior probability distributions for the mass density slope parameters 𝛾0, 𝛾𝑧 , and 𝛾𝑠. Diagonal panels show the one-dimensional marginalised posterior distributions, with dashed lines indicating the 68% credible intervals and the mean value. Off-diagonal panels show the two-dimensional joint posterior distributions, with contours marking the 68% and 95% credible regions. The constraints, derived from t… 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

1 extracted references

  1. [1]

    H., 2018a, MNRAS, 474, 3391 Agnello A., et al., 2018b, MNRAS, 475, 2086 Agnello A., Lin H., Kuropatkin N., Buckley-Geer E., Anguita T., Schechter P

    Agnello A., Spiniello C., 2019, MNRAS, 489, 2525 Agnello A., et al., 2015, MNRAS, 454, 1260 Agnello A., Grillo C., Jones T., Treu T., Bonamigo M., Suyu S. H., 2018a, MNRAS, 474, 3391 Agnello A., et al., 2018b, MNRAS, 475, 2086 Agnello A., Lin H., Kuropatkin N., Buckley-Geer E., Anguita T., Schechter P. L., et al., 2018c, MNRAS, 479, 4345 Anguita T., Faure...

    arXiv 2019