Time delay measurements with Broken Power Law model
Pith reviewed 2026-05-16 14:56 UTC · model grok-4.3
The pith
A Broken Power Law lens mass model for WGD 2038-4008 yields H0 of 75.3 km s^{-1} Mpc^{-1} when the internal mass-sheet factor varies, compared to 61.1 from the elliptical power-law model.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
When the internal mass-sheet factor is allowed to vary, the inferred Hubble constant in a flat Lambda-CDM cosmology with fixed Omega_m=0.3 is H0 = 75.3^{+23.1}_{-16.3} km s^{-1} Mpc^{-1} for the BPL model and H0 = 61.1^{+19.2}_{-13.2} km s^{-1} Mpc^{-1} for the EPL model; with the factor fixed to unity the values become 74.2^{+20.3}_{-13.8} and 66.1^{+18.8}_{-12.8} respectively.
What carries the argument
The Broken Power Law (BPL) mass model implemented in Lenstronomy, which parametrizes the radial mass distribution with two power-law indices to capture deviations from a single power-law profile, paired with numerical evaluation of time delays at the image positions.
If this is right
- Both BPL and EPL models reproduce the observed imaging and stellar kinematics of WGD 2038-4008, with a modest statistical preference for BPL.
- Allowing the internal mass-sheet factor to vary amplifies the difference between the two models' H0 posteriors.
- Time-delay cosmography remains sensitive to the choice of lens mass parametrization even after kinematic constraints are included.
Where Pith is reading between the lines
- If future data on more lenses consistently favor BPL over EPL, the current tension between lensing H0 values and early-Universe determinations may shift depending on the model family adopted.
- Cross-checking BPL-derived mass profiles against independent stellar-dynamics or weak-lensing constraints on the same galaxies would test whether the extra freedom reduces systematic error.
- Extending the numerical time-delay module to other flexible profiles could quantify how much of the H0 spread arises from radial-profile assumptions versus other modeling choices.
Load-bearing premise
The Broken Power Law model accurately represents the true radial mass distribution of the lens galaxy and the numerical time-delay computation adds no biases beyond those already present in standard models.
What would settle it
An independent high-resolution measurement of the radial mass-density slope of the lens galaxy in WGD 2038-4008 that lies outside the BPL posterior range would invalidate the higher H0 value obtained with that model.
read the original abstract
One challenge in strong gravitational lensing cosmography is the measurement of time delays between multiple lensed images, which are essential for constraining the Hubble constant (\(H_0\)). In this study, we investigate how assumptions about the lens mass profile affect time-delay measurements in lensing systems. Specifically, we implement a Broken Power Law (BPL) mass model within the \textsc{Lenstronomy} framework (Birrer & Amara 2018), which introduces additional flexibility in the radial mass distribution and can phenomenologically capture deviations from a single power-law profile. This model is combined with a numerical approach to compute time delays at the image positions. We validate the BPL implementation using simulated lens systems and compare the results with those from the elliptical power-law (EPL) model. We then apply both model families to the quadruply imaged quasar WGD~2038--4008. Both models provide good fits to the imaging and kinematic data, with a slight preference for the BPL model. When the internal mass-sheet factor is allowed to vary, the inferred Hubble constant in a flat \(\Lambda\)CDM cosmology with fixed \(\Omega_{\rm m}=0.3\) is \(H_0 = 75.3^{+23.1}_{-16.3} \ \mathrm{km \ s^{-1} \ Mpc^{-1}}\) for the BPL model and \(H_0 = 61.1^{+19.2}_{-13.2} \ \mathrm{km \ s^{-1} \ Mpc^{-1}}\) for the EPL model. For comparison, in the diagnostic case with the internal mass-sheet factor fixed to unity under the same setup, we obtain \(H_0 = 74.2^{+20.3}_{-13.8} \ \mathrm{km \ s^{-1} \ Mpc^{-1}}\) for the BPL model and \(H_0 = 66.1^{+18.8}_{-12.8} \ \mathrm{km \ s^{-1} \ Mpc^{-1}}\) for the EPL model. This highlights how time-delay cosmography remains sensitive to assumptions about the lens mass profile. With current precision, this difference does not favor one cosmological scenario over another, but underscores the importance of flexible mass modeling.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript implements a Broken Power Law (BPL) lens mass model in Lenstronomy, validates it on simulated lenses, and applies both BPL and elliptical power-law (EPL) models to the quadruply imaged quasar WGD 2038-4008. With the internal mass-sheet factor free, it reports H0 = 75.3+23.1−16.3 km s−1 Mpc−1 for BPL versus 61.1+19.2−13.2 km s−1 Mpc−1 for EPL in flat ΛCDM with fixed Ωm=0.3; the values shift modestly when the mass-sheet factor is fixed to unity. The work concludes that time-delay cosmography remains sensitive to radial mass-profile assumptions.
Significance. If the BPL implementation and its time-delay numerics prove unbiased, the result quantifies how additional radial flexibility alters the inferred Hubble constant at the level of the current error bars, reinforcing that mass-model choice is a leading systematic in strong-lensing cosmography. The large reported uncertainties prevent any decisive cosmological statement, but the explicit comparison supplies a concrete diagnostic for future analyses.
major comments (2)
- [Validation on simulations] Validation section: the simulations are stated to show good fits, yet no quantitative recovery statistics (bias, scatter, or coverage) are supplied for time delays or H0 when the internal mass-sheet factor is left free—the exact regime that produces the headline H0 offset between BPL and EPL.
- [Numerical time-delay approach] Time-delay computation: the numerical evaluation of the Fermat potential at image positions for the BPL profile (with its extra break radius and slope parameters) is not demonstrated to be free of profile-dependent truncation or interpolation errors relative to the EPL case, which directly enters the time-delay distance and therefore H0.
minor comments (2)
- [Abstract] The abstract states “a slight preference for the BPL model” without specifying the quantitative metric (Δχ², Bayesian evidence, or information criterion).
- The reported H0 uncertainties are given to one decimal place; the text does not detail how the posterior is marginalized over the additional BPL parameters and the mass-sheet factor.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major point below and have revised the manuscript to incorporate additional quantitative validation and numerical checks as requested.
read point-by-point responses
-
Referee: [Validation on simulations] Validation section: the simulations are stated to show good fits, yet no quantitative recovery statistics (bias, scatter, or coverage) are supplied for time delays or H0 when the internal mass-sheet factor is left free—the exact regime that produces the headline H0 offset between BPL and EPL.
Authors: We agree that quantitative recovery statistics are essential to substantiate the validation, especially for the free mass-sheet case that drives the reported H0 difference. In the revised manuscript we have added a new table (Table 2) reporting bias, scatter, and 68% coverage for recovered time delays and H0 on the simulated sample, separately for BPL and EPL with the internal mass-sheet factor left free. The statistics show that both models recover the input H0 within 1σ on average, with the BPL model exhibiting modestly larger scatter consistent with its additional degrees of freedom. revision: yes
-
Referee: [Numerical time-delay approach] Time-delay computation: the numerical evaluation of the Fermat potential at image positions for the BPL profile (with its extra break radius and slope parameters) is not demonstrated to be free of profile-dependent truncation or interpolation errors relative to the EPL case, which directly enters the time-delay distance and therefore H0.
Authors: We acknowledge that explicit demonstration of numerical robustness for the BPL Fermat-potential evaluation was not provided. The implementation re-uses Lenstronomy’s adaptive quadrature routine without additional truncation or interpolation steps specific to the break radius. In the revised manuscript we have added a dedicated subsection (Section 3.3) that compares numerical time delays against high-resolution reference calculations for both BPL and EPL on the same simulated lenses, confirming fractional differences below 0.05% that do not propagate into H0 at the reported uncertainty level. We have also included convergence tests varying the quadrature tolerance. revision: yes
Circularity Check
No significant circularity; derivation remains self-contained
full rationale
The paper implements the BPL profile inside Lenstronomy, validates the numerical time-delay integrator on independent simulated lenses (recovering known inputs), then fits the same functional form plus external time-delay measurements and imaging/kinematics to WGD 2038-4008. The reported H0 values are downstream outputs of that fit under two different radial profiles; they are not re-injected as inputs, nor is any parameter fitted to a subset and then relabeled a prediction. No self-citation supplies a uniqueness theorem or ansatz that the present work merely renames. The mass-sheet degeneracy is explicitly varied and its effect on H0 is shown, which is the intended scientific comparison rather than a definitional loop.
Axiom & Free-Parameter Ledger
free parameters (2)
- BPL break radius and inner/outer slopes
- internal mass-sheet factor
axioms (2)
- domain assumption The lens mass distribution can be adequately described by a broken power-law profile
- domain assumption Numerical computation of time delays from the lens potential yields unbiased values at image positions
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
BPL density κ(R_el) piecewise with inner slope α_c, outer α, break r_c; numerical ψ via u-integration of φ_r(ξ) (Eqs. 1,11); λ_int rescaling of time delays (Eq. 23)
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Validation of Fermat-potential solver on EPL analytic limit; BIC-weighted posteriors for H0 under free λ_int
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.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.