{"id":"e424b2d2-5dcd-4060-8f46-d4948f0dc843","arxiv_id":"2501.16292","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"The first measurement of line of sight shear in 23 real strong lenses gives a typical magnitude of about 0.056, larger than expectations from N-body simulations.","lead":"Astronomers modelled 23 strong gravitational lenses and, for the first time, extracted the line of sight shear, the small extra bending of light by matter between us and the source galaxy. The typical shear is about 0.056, larger than cosmic web simulations predict, and the comparison shows that the main uncertainty is the lens galaxy's own mass model.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mass-model complexity beyond the EPL is the load-bearing untested assumption; the m=4 octupole test is inconclusive, so the 0.056 mean may include deflector systematics masquerading as LOS shear.","rationale":"The reader's weakest assumption is exactly the condition I judge most load-bearing: unmodelled main-deflector complexity biasing γ_LOS. The paper's own framing makes the main-deflector model the dominant systematic, and the validation in Hogg et al. (2023) only used mocks whose complexity matched the fitted model, so it does not test the mismatch regime. The octupole experiment is the only direct probe in this work, and its outcome is mixed: one lens is better fit by the octupole with reduced shear, while other apparent improvements come from inflated uncertainties rather than better fits; meanwhile, the recovered |a4| values are often larger than observed isophotal octupoles, suggesting the model is absorbing something else. A controlled mock with known input shear and an ellipticity gradient would settle whether the EPL + minimal-shear model recovers γ_LOS or biases it. Secondary statistical issues (posterior mean/median of |γ| is positively biased for low-SNR lenses; abstract vs. main text uncertainty 0.013 vs 0.016; Gaussian tension statistic applied to non-Gaussian posteriors) affect the numeric significance but not the qualitative finding, and they are already noted by the reader. Since the reader's verdict is already CONDITIONAL and my concern is the same one, no verdict change is needed.","tokens_in":24859,"tokens_out":12660,"duration_ms":124588,"concrete_test":"Using the released dolphin/lenstronomy pipeline and public SLACS inputs, simulate the 23 systems with identical noise, PSF, and source models; inject a known LOS shear of γ_LOS = 0.03 at a random orientation and a deflector with an ellipticity gradient (axis ratio varying from q=0.9 at 0.2 θ_E to q=0.7 at 2 θ_E). Fit with the default EPL + minimal-shear model. If recovered γ_LOS differs from the injected value by more than about 0.01, the unmodelled-complexity concern lands.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline measurement, |γ_LOS| mean = 0.056, is only interpretable as line-of-sight shear if the EPL main deflector plus (optionally) an m=4 multipole captures the deflector's mass distribution well enough that residual complexity does not leak into the five shear parameters. This condition is not established. Section 3.2 calls the EPL 'a relatively simplistic description', and Section 5 concludes 'It is the modelling of the main deflector to which attention must now be devoted.' The only complexity probed is the m=4 octupole (Section 4.3), yet the octupole test is inconclusive: |a4| is recovered at or above 2% in a large fraction of lenses, well above the ~1% observed in early-type galaxy isophotes, and for SDSSJ1204+0358 the octupole model finds a different solution with |a4| consistent with zero, suggesting model instability rather than a clean test. Real deflectors may contain m=1 and m=3 multipoles, ellipticity gradients, isophotal twists, or substructure; the paper explicitly defers these to future work. If such features are present, they can be absorbed into γ_LOS and ω_LOS, and the reported mean would not be a LOS shear measurement.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents the first application of the 'minimal' line-of-sight (LOS) shear model of Fleury et al. (2021) to real strong-lens imaging, modelling 23 SLACS lenses with the dolphin/lenstronomy pipeline using an elliptical power-law (EPL) main deflector, two Sérsic profiles for deflector light, and a Sérsic plus shapelets model for source light. The authors report per-lens LOS shear magnitudes, a sample mean |γ_LOS| = 0.056 (quoted as ±0.016 in §4.2 and ±0.013 in the abstract), compare these values with RayGalGroup N-body weak-lensing expectations, and examine two systematic variations: setting the foreground shear γ_od to zero and adding an m=4 octupole to the mass profile. They conclude that the minimal model should be used, that the octupole does not generally reduce the inferred shear, and that future effort should focus on main-deflector modelling.","tokens_in":25210,"tokens_out":5421,"duration_ms":48201,"significance":"If the measurement is robust, this is a genuinely new probe: strong-lens image shapes would constrain a shear quantity that is theoretically free of the ellipticity degeneracy that plagues conventional external shear, and the comparison with N-body predictions would establish a new test of large-scale structure. The paper is exemplary in reproducibility: the pipeline, inputs, MCMC chains, and figure scripts are publicly released, and the comparison to RayGalGroup simulations is an independent, falsifiable check. The authors are also candid about the main limitation — the simplicity of the EPL deflector — and their negative octupole result is a useful data point for the field. However, the headline value is only interpretable as a LOS shear if unmodelled deflector complexity does not leak into the five shear parameters; this condition is not yet established by the tests in the paper.","major_comments":[{"comment":"The abstract states a mean LOS shear of 0.056 ± 0.013, while Section 4.2 states 0.056 ± 0.016; these are different statements about the same quantity and must be reconciled. Please specify how the mean and its uncertainty are computed (e.g., sample mean of per-lens posterior medians, standard error, or inverse-variance weighted mean including asymmetric errors) and correct the inconsistent value.","section":"Abstract and §4.2"},{"comment":"The Einstein radius θ_E is fixed to the best-fit value from Shajib et al. (2021) and is not sampled, so the reported per-lens shear uncertainties do not include uncertainty in θ_E. Because θ_E is tightly constrained, the effect may be small, but this is an assumption rather than a demonstrated result; please quantify the effect by re-fitting at least a subset with θ_E free or by marginalizing over the Shajib et al. posterior.","section":"§3.2"},{"comment":"The octupole test is inconclusive as a control for deflector complexity. The recovered |a_4| values are larger than the ~1% level observed in early-type isophotes (35% of lenses at or above 2%, with a further 45% above 5%), which is itself a red flag that the m=4 term is absorbing something else; for SDSSJ1204+0358 the octupole model finds a different solution with |a_4| consistent with zero but χ² increases from 0.68 to 1.08, indicating model instability rather than a clean null test. The m=4 multipole is also only one of several plausible complexities (m=1, m=3, ellipticity gradients, twists, substructure), so the conclusion that 'shear does not compensate for unmodelled octupoles' does not establish that the reported mean is free of deflector systematics. Please add recovery tests on mock images with realistic deflector complexity (not just EPL+octupole) and report the induced bias in γ_LOS.","section":"§4.3"},{"comment":"The comparison with RayGalGroup simulations uses the mean expected |γ_LOS| for each (z_d, z_s) pair, but the distribution of |γ_LOS| is positive-definite and likely non-Gaussian, so the tension statistic in Eq. (14) may overstate the significance of a 7.9σ outlier; please report the full simulated distribution of |γ_LOS| and the percentile of each measurement. In addition, the simulation expectation does not condition on the SLACS selection function (e.g., group environments), which is relevant to the claim that measured values are larger than expected.","section":"§4.2 and Fig. 6"}],"minor_comments":[{"comment":"The per-lens |γ_LOS| values in Table 1 disagree with those printed in the corresponding figure panels for several lenses: for example, SDSSJ0029-0055 is 0.029+0.009−0.013 in Table 1 but 0.034+0.012−0.013 in Fig. 1, SDSSJ0330-0020 is 0.049+0.013−0.012 in the table but 0.059+0.012−0.014 in Fig. 1, and SDSSJ0728+3835 is 0.042+0.017−0.013 in the table but 0.028+0.013−0.011 in Fig. 1. These inconsistencies must be resolved or explained.","section":"Table 1 and Figs. 1–4"},{"comment":"The text refers to 'SDSSJ1204+0385' in three places; the correct name is SDSSJ1204+0358.","section":"§4.3"},{"comment":"For SDSSJ1204+0358 the paper reports '|a_4| = −0.02 +0.07 −0.04', which is impossible for a magnitude; please clarify whether a_4 is a signed amplitude and whether the reported quantity is a_4 or |a_4|.","section":"§4.3"},{"comment":"The footnote acknowledges that ordinary least squares is inaccurate when the independent variable has significant uncertainties, but the linear fit in Fig. 10 is still presented without an errors-in-variables treatment; please either use a suitable regression method or label the fit as illustrative only.","section":"Footnote 4"}],"recommendation":"major_revision","confidential_remarks":"This is a well-executed and candid paper, and I would be happy to see it published after revision. The main risk is that the headline 'first measurement of LOS shear' may be over-interpreted: the authors' own conclusion about main-deflector modelling is the right caution, and the octupole test does not yet close the gap. I would ask the editor to insist on the abstract/section consistency and a clear error budget, and to request at least one robustness test (e.g., mock recovery with realistic deflector complexity) before acceptance. The comparison with simulations is valuable, but the tension statistic needs to be based on the full simulated distribution rather than a mean and Gaussian width."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is the first application of the minimal LOS shear model to real strong-lens images, and it is a solid piece of work. The formalism comes from Fleury et al. (2021) and the mock validation from Hogg et al. (2023), but actually carrying this through 23 SLACS lenses with public data, code, and MCMC chains is a genuine step forward. The octupole and no-foreground-shear comparisons are useful, even where their conclusions are negative, and the writing is refreshingly honest about what is and is not established.\n\nThe headline number, mean |γ_LOS| = 0.056, looks believable as a model-dependent measurement. The fits are careful, the reduced chi-squared values are reasonable, and the comparison with RayGalGroup simulations is a sensible first check. But I would not yet call this a clean line-of-sight shear measurement. The paper's own conclusion says it best: attention must now go to the main deflector. The EPL plus optional octupole is a restricted family, and m=1, m=3, ellipticity gradients, isophotal twists, and substructure are all deferred. The octupole test is not a decisive null: many lenses have |a4| above the isophotal expectation, and SDSSJ1204+0358 shows clear model instability. So the 0.056 mean should be read as \"shear inferred under EPL-based models,\" not as a proven cosmological signal.\n\nThere are also a few technical soft spots that a referee should press on. The abstract quotes 0.056 ± 0.013 while Section 4.2 says ±0.016; that inconsistency in the headline number should be fixed. The Einstein radius is fixed from Shajib et al. (2021) without propagating its uncertainty, which is probably minor but should be justified. The posterior mean of |γ_LOS| is a positively biased estimator for small true shear, and this could inflate the sample mean; the paper does not quantify the effect. The Gaussian tension statistic is applied to posteriors that are visibly non-Gaussian and truncated at zero, so the 7.9σ outlier is likely overstated. None of these invalidates the central measurement, but they matter for future cosmological use.\n\nThe self-citation point is a non-issue here: the minimal model is this group's construction, and the mock validation was done properly in earlier papers. The data and code release is exemplary.\n\nThis paper deserves a serious referee. It is a legitimate first measurement with real data, and the community needs this kind of systematic test. A referee should ask for the uncertainty inconsistency to be resolved, a quantification of the |γ| bias, and a more careful treatment of mass-model complexity, but the work should not be desk-rejected.","headline":"First real-data measurement of the minimal-model LOS shear, done carefully and honestly, but the main-deflector model is the acknowledged unknown and the headline uncertainty has a fixable inconsistency.","tokens_in":25779,"tokens_out":1674,"would_cite":true,"duration_ms":18064,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.62.Sb"],"model":"deepseek-v4-flash","headline":"The paper reports the first measurement of line-of-sight shear from real strong-lens images, with a mean magnitude of 0.056 across 23 SLACS lenses.","keywords":["strong gravitational lensing","line-of-sight shear","minimal shear model","SLACS","elliptical power-law profile","octupole multipole","post-Born corrections","cosmic shear"],"falsifier":"Measure the fourth-order isophotal distortions of the 23 lens galaxies from their own images: the fits recover octupole strengths of 2% or more in a third of the sample, whereas observational studies of early-type galaxies find roughly 1% distortions, and a mismatch would show the mass model is absorbing non-octupolar complexity rather than real boxiness. A complementary check is to compare each measured $\\gamma_{\\rm LOS}$ against shear reconstructed from background galaxy shapes along the same sightline, since a systematic excess in the lens-based values would signal leakage from the deflector model.","tokens_in":24618,"feed_emoji":"🔭","tokens_out":12317,"duration_ms":94834,"temperature":0.7,"pith_summary":"This paper attempts the first measurement of line-of-sight (LOS) shear in real strong gravitational lens images, using the 'minimal model' designed to be free from the degeneracy between shear and the lens galaxy's mass ellipticity. Fitting 23 lenses from the SLACS catalogue with an automated modelling pipeline, it reports a mean LOS shear magnitude of $0.056 \\pm 0.013$, a value larger than what N-body simulations of the cosmic web predict for the same lines of sight. The paper also shows that dropping the foreground shear — the post-Born correction to the deflector's potential — biases the inferred shear, and that adding an octupole term to the lens mass does not systematically bring measurements into agreement with simulations. If the result holds, strong-lens images carry a measurable weak-lensing signal that could complement cosmic shear as a cosmological probe, with the modelling of the lens galaxy itself now the dominant systematic.","feed_headline":"First measurement of line-of-sight shear in strong lenses","feed_subtitle":"Fitting 23 SLACS lenses gives mean shear 0.056, exceeding cosmic-web simulation predictions.","key_machinery":"The central object is the minimal LOS shear model. It begins with three separate tidal shears acting on the light path — between observer and source, deflector and source, and observer and deflector — and multiplies the lens equation by $A_{\\rm od} A_{\\rm ds}^{-1}$, a source-position transformation that absorbs the individually unmeasurable shears into a single LOS amplification matrix. What survives in the model is the combination $\\gamma_{\\rm LOS} \\equiv \\gamma_{\\rm od} + \\gamma_{\\rm os} - \\gamma_{\\rm ds}$, argued to be free of degeneracy with the lens mass ellipticity at leading order, alongside a foreground shear $\\gamma_{\\rm od}$ acting through the effective deflector potential $\\psi_{\\rm eff}(\\theta) = \\psi(A_{\\rm od}\\theta)$ and a LOS rotation $\\omega_{\\rm LOS}$. The pipeline fits each lens with an elliptical power-law mass profile, two Sérsic profiles for the deflector light, and a Sérsic-plus-shapelet source model, sampling the posterior with MCMC and reading $\\gamma_{\\rm LOS}$ off the recovered shear components.","core_discovery":"The central claim is that the line-of-sight shear of the minimal model, $\\gamma_{\\rm LOS} \\equiv \\gamma_{\\rm od} + \\gamma_{\\rm os} - \\gamma_{\\rm ds}$, can be recovered from real strong-lens images, and the paper reports the first such measurement: across 23 SLACS lenses the mean magnitude is $0.056 \\pm 0.013$. The individual measurements are said to be consistent with the external shears fitted to the same lenses in earlier work, but the sample mean sits well above the expectations drawn from N-body simulation light cones, with several lenses in tension at more than $3\\sigma$. The paper further claims that fixing the foreground shear to zero (the no-foreground-shear model) propagates degeneracies and shifts the inferred LOS shear, so the full minimal model must be kept for shear science, and that including an octupole ($m=4$) term in the lens mass profile rarely lowers the shear — only SDSSJ1531-0105 is better fit by the octupole model, which reduces its shear by just over $1\\sigma$. Its stated conclusion is that the line-of-sight contribution is adequately captured by the model, and that attention must now shift to the mass model of the main deflector.","pith_inferences":["A direct cleanliness test that the paper does not run: compare each measured $\\gamma_{\\rm LOS}$ with shear reconstructed from galaxy shapes (cosmic shear) in the same sky patch; a systematic excess in the lens-based values would show that part of the 'LOS' signal is still mass-model leakage.","The recovered octupole strengths could be checked against independent isophotal measurements of the same galaxies; a mismatch would indicate that the mass model is absorbing something other than true boxiness or diskyness, with the same leakage possibly contaminating $\\gamma_{\\rm LOS}$.","Extending the minimal-model pipeline to the full 50-lens SLACS sample, or beyond, would turn the mean $|\\gamma_{\\rm LOS}|$ into a distribution-wide comparison with simulations and would clarify whether the high-shear tail (for example SDSSJ1112+0826 at about 0.23) is physical or pathological.","If later, more flexible mass models bring the mean shear down toward the roughly 0.01–0.02 level predicted by simulations, the excess reported here would itself become a diagnostic of how much complexity the elliptical power-law fits miss."],"forward_implications":["If the measurement is right, strong-lens images can be used to measure weak-lensing shear along individual lines of sight, and the mean $|\\gamma_{\\rm LOS}| \\approx 0.056$ becomes the first observational anchor for the signal.","The measured shear being larger than cosmic-web simulations predict independently confirms, for the LOS formalism, the known discrepancy between strong-lens external shear and weak-lensing expectations.","Pipelines that aim to measure shear must keep the foreground shear parameters: fixing $\\gamma_{\\rm od} = 0$ trades a small speed gain for a biased LOS shear and artificially tightens the recovered lens ellipticity.","Because octupole flexibility does not systematically reduce the shear, the explanation for large shear values must lie elsewhere in the mass model, so adding boxiness and diskyness is not a cure.","The result lays the groundwork for using strong-lens shear as a cosmological probe, but only after the mass-model systematics are controlled."],"supporting_citations":[{"why":"Constructs the minimal LOS shear model and proves the source-position transformation that isolates the degeneracy-free combination $\\gamma_{\\rm LOS}$.","marker":"Fleury et al. (2021)"},{"why":"Demonstrates with mock strong-lens images that $\\gamma_{\\rm LOS}$ is recoverable and non-degenerate with lens mass parameters, motivating the real-data attempt.","marker":"Hogg et al. (2023)"},{"why":"First analysis to extract shear from a strong-lens image, whose observed degeneracies between shear and lens ellipticity motivate the minimal-model approach.","marker":"Birrer et al. (2017)"},{"why":"Supplies the 23-lens SLACS sample, the previously measured external shears used for comparison, and the automated modelling pipeline.","marker":"Shajib et al. (2021)"},{"why":"Documents external shears that are suspiciously large relative to weak-lensing expectations and proposes octupolar mass complexity as the culprit, the hypothesis tested here.","marker":"Etherington et al. (2024)"},{"why":"Provides the elliptical power-law mass profile used as the baseline lens model.","marker":"Tessore and Metcalf (2015)"},{"why":"RayGalGroup N-body simulations whose light cones yield the expected LOS shear used as the benchmark for the measurements.","marker":"Breton et al. (2019)"},{"why":"Extends the RayGalGroup suite whose weak-lensing maps, including post-Born effects, set the simulation expectations for $|\\gamma_{\\rm LOS}|$.","marker":"Rasera et al. (2022)"}],"fun_headline_variants":["First LOS shear measured: 0.056, above cosmic-web predictions","23 SLACS lenses yield LOS shear 0.056, exceeding simulations","First LOS shear measurement: 0.056, tension with N-body","LOS shear from 23 lenses: 0.056, above cosmic-web simulations","First measured LOS shear: 0.056, exceeds cosmic-web forecasts"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the smooth elliptical power-law mass profile — and, in the extended model, a single octupole term — describes each lens galaxy well enough that any unmodelled complexity in the deflector, such as ellipticity gradients, isophotal twists, other multipoles, or substructure, does not mimic shear and leak into the measured $\\gamma_{\\rm LOS}$.","fun_headline_variants_meta":{"raw":{"variants":["First LOS shear measured: 0.056, above cosmic-web predictions","23 SLACS lenses yield LOS shear 0.056, exceeding simulations","First LOS shear measurement: 0.056, tension with N-body","LOS shear from 23 lenses: 0.056, above cosmic-web simulations","First measured LOS shear: 0.056, exceeds cosmic-web forecasts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00126,"raw_usage":{"total_tokens":5198,"prompt_tokens":1023,"completion_tokens":4175,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":639,"completion_tokens_details":{"reasoning_tokens":4077}},"tokens_in":639,"tokens_out":4175,"duration_ms":26844,"temperature":1.0,"reasoning_tokens":4077,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T13:34:50.765999+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the fourth-order isophotal distortions of the 23 lens galaxies from their own images: the fits recover octupole strengths of 2% or more in a third of the sample, whereas observational studies of early-type galaxies find roughly 1% distortions, and a mismatch would show the mass model is absorbing non-octupolar complexity rather than real boxiness. A complementary check is to compare each measured $\\gamma_{\\rm LOS}$ against shear reconstructed from background galaxy shapes along the same sightline, since a systematic excess in the lens-based values would signal leakage from the deflector model.","supporting_citations":[],"review_version":1}