{"id":"8cc864c9-a02f-4715-8e56-cf6ca6eaf364","arxiv_id":"2411.10594","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A new cosmic shear estimator shears the weight function instead of the image and derives an analytic noise-bias correction, but the noise-bias correction for the response is not yet complete.","lead":"This paper introduces a new way to measure cosmic shear, the subtle distortion of galaxy shapes by dark matter, that does not require assumptions about galaxy shapes and includes an analytical noise correction. It could help future surveys like LSST reach the precision needed for dark energy measurements.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's own Section 4 admits the noise-bias correction for the response R is still work in progress; without it the headline 'unbiased' is not established.","rationale":"Good-faith reading: this is a methods paper proposing a clever coordinate-shear estimator. The formal equivalence behind Eqs. 5-6 is defensible in the continuous limit, and the noise-free results are encouraging. However, the paper itself contains an explicit limitation statement in Section 4: the noise-bias correction for the distorted moments M^S, needed to correct R, is still work in progress. Because the estimator's output is ⟨R⟩^{-1}⟨e⟩, a biased R means a biased shear regardless of how well e is corrected. The reader's weakest_assumption focused on the coordinate-shear equivalence; I think that equivalence is the more likely-to-be-formally-ok part, while the admitted missing R correction is the more immediate, textually explicit gap. Hence partial agreement with the reader. The recommended verdict remains CONDITIONAL: the claim of unbiasedness should be withheld until the R correction is derived and validated by the end-to-end noisy test described above.","tokens_in":3939,"tokens_out":4979,"duration_ms":49580,"concrete_test":"Run a dedicated noisy end-to-end simulation with the current pipeline: fixed known shear g_true, Galsim galaxy profiles (elliptical and COSMOS), PSF profiles as in Section 3, and Gaussian noise at the σ_noise levels of Fig. 2 (e.g., 3 to 5 in the paper's units). For each noise level, generate many realizations with random intrinsic ellipticities, estimate M and the four distorted M^S, apply the existing M noise correction, form R via Eq. 8 and ⟨g⟩ via Eq. 9 without any R correction. Measure the multiplicative bias per component, m_i = ⟨g_i⟩/g_true,i − 1, and its statistical error. If |m_i| exceeds 10^{-3} at the target noise, the uncorrected R bias invalidates the 'unbiased' claim as stated; if it does not, quantify that the omitted term is subdominant and report it.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the estimator is unbiased. But the shear estimate is ⟨g⟩ = ⟨R⟩^{-1}⟨e⟩ (Eq. 9), so any noise bias in R translates directly into multiplicative bias. Section 4 derives and tests an analytical noise correction for the second moments M (hence for e, Fig. 2), but explicitly states: 'The same kind of noise bias analytical formulas needs to be calculated for the distorted M^S, in order to correct R, but this is still a work in progress.' No noisy end-to-end test of ⟨g⟩ is shown. The noise-free tests of Section 3 and the e-only correction of Fig. 2 therefore do not support the Conclusion's 'unbiased' claim. This is not a hypothetical failure mode: the estimator is nonlinear in the image through the measured centroid x0, and the distorted moments M^S used in the derivatives in Eq. 8 are evaluated at sheared coordinates on the same noisy image, so the noise enters R in a correlated way that cannot be assumed to cancel. A secondary issue is the practical implementation of Eqs. 5-6 (sampled PSF, Fourier division, interpolation on distorted grids), whose δ' cross-term is mentioned but not shown to be subtracted in the reported tests; however, the decisive, paper-admitted gap is the missing R noise-bias correction.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a cosmic shear estimator built from weighted second moments, where the shear response matrix R is computed by applying shear to the coordinate system of the weight kernel rather than to the galaxy image. The authors derive an analytical noise-bias correction for the second moments and test the estimator on noise-free simulations with Gaussian and COSMOS galaxies, reporting multiplicative biases below 1e-3. They also show that the analytical correction improves noisy ellipticity estimates. The Conclusion claims an unbiased self-calibrated estimator, but the paper explicitly states that the analogous noise-bias correction for the sheared moments entering R is still work in progress, and no noisy end-to-end shear test is presented.","tokens_in":4396,"tokens_out":3029,"duration_ms":31791,"significance":"If the full correction were completed and validated, the method would be attractive because it avoids galaxy-profile assumptions, does not shear the image, and offers an analytic route to noise-bias correction. The self-calibration scheme via R is a standard and conceptually sound approach. The paper is honest about its current limitations, and the noise-free residuals are encouraging, but the central claim of an unbiased estimator on noisy images is not yet supported by the evidence presented.","major_comments":[{"comment":"The paper states: 'The same kind of noise bias analytical formulas needs to be calculated for the distorted M^S, in order to correct R, but this is still a work in progress.' Because the estimator in Eq. (9) is ⟨g⟩ = ⟨R⟩^{-1}⟨e⟩, any noise bias in R translates directly into a multiplicative bias on g. No noisy end-to-end test of ⟨g⟩ is shown. Therefore the Conclusion's claim of an 'unbiased self-calibrated shear estimator' is not established by the present results. This is a load-bearing gap, not a cosmetic one.","section":"Section 4, last paragraph"},{"comment":"The residuals in Figure 1 are presented without error bars or scatter estimates. With only 40 random shear values and 20 pairs of intrinsic ellipticities, the claim that the bias is below the 10^-3 requirement needs an uncertainty quantification. As written, the plots show small residuals but do not demonstrate statistically that the bias is below 10^-3; a single realization with 40 points has large sampling variance.","section":"Section 3, Figure 1"},{"comment":"The equivalence between shearing the coordinate system of F and shearing the galaxy image relies on the Fourier-space division by the PSF and interpolation on distorted grids. The paper mentions a cross-effect term δ' in Eq. (7) that must be subtracted, but it never states explicitly whether this correction was applied in the noise-free tests of Section 3. If it was not, the reported residuals could be hiding a known systematic; if it was, the implementation should be described. This ambiguity affects the interpretation of the validation results.","section":"Section 2, Eqs. (5) and (6)"},{"comment":"The analytical noise-bias correction is derived for the second moments M and validated only for the ellipticity e. The response matrix R, however, depends on derivatives ∂M/∂g_i computed from sheared moments M^S, whose noise properties are correlated with the centroid and with the unsheared moments. The paper does not address whether the same correction formulas apply to these derivatives. This reinforces the concern that the missing R correction is not a minor extension but a nontrivial step.","section":"Section 4, Eqs. after Fig. 2"}],"minor_comments":[{"comment":"The author name 'Enya V an den Abeele' contains an erroneous space; it should read 'Enya Van den Abeele'.","section":"Author list and title"},{"comment":"The phrase 'under our of 10^-3 upper limit' is a typo; it should read 'under our 10^-3 upper limit'.","section":"Section 3"},{"comment":"The x-axis label 'noise' is vague; please specify whether it is σ_noise, the noise variance, or something else, and give units.","section":"Figure 2"},{"comment":"The symbols s and ϵ are used in the expansion but their definitions appear only later in the text; please define them near the equation.","section":"Eq. (7)"},{"comment":"Reference 3 is incomplete (only a DOI is given); please provide authors, title, and publication details.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is refreshingly transparent about its limitations, but the Conclusion overstates what has been demonstrated. The missing noise-bias correction for R is the decisive point; until it is supplied or the claim is explicitly softened, the paper cannot support the headline 'unbiased'. The lack of error bars in the noise-free tests is also a basic statistical omission that should be addressed in a revision. Given the work-in-progress nature, a major revision is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The one thing to know: this paper has a genuinely new idea. Instead of shearing the galaxy image like Metacalibration does, it shears the weight/PSF kernel F in the estimator (Eqs. 4-6). That is a real variant with practical appeal—no image distortion, no correlated noise from resampling the galaxy, and a path to handling undersampled images. The noise-free tests on Gaussian and COSMOS galaxies with different PSFs look encouraging, and the analytic noise-bias expansion for the second moments is a solid piece of work that matches simulations in Fig. 2. The paper is also honest: it explicitly says the response R needs its own noise correction and that this is work in progress.\n\nThe soft spots are in proportion. The load-bearing one is the missing R correction. The estimator is ⟨g⟩ = ⟨R⟩^{-1}⟨e⟩, and R is a nonlinear function of the noisy image—through the centroid and through the sheared coordinates—so noise in R does not cancel by symmetry. The paper's own Section 4 admits this, and no end-to-end noisy test of ⟨g⟩ is shown. Without it, the word \"unbiased\" in the Conclusion is not established. A secondary issue is that the noise-free residuals have no error bars; the \"below 1e-3\" claim is qualitative, though the scatter in the figures suggests it is probably fine. The δ' cross-term between shear and sampling is mentioned but not shown to be subtracted in the reported tests, so the practical implementation (Fourier division, interpolation on distorted grids) is under-tested. These are all fixable, and the paper already identifies the main one.\n\nThis is a short work-in-progress paper, and it reads like one. The derivations are sketched, but the core idea is clear and the math is coherent. I do not see a circular step: R is measured from the estimator under artificial shears, which is standard self-calibration, and the analytic correction is tested against simulations, not fit to them. The citation pattern is fine—Metacalibration and GalSim are the relevant references.\n\nWho gets value: weak-lensing method people working on shear calibration for LSST-era surveys. It deserves a serious referee, not a desk reject. The referee should demand the R noise-bias correction and a statistically quantified end-to-end noisy test before the \"unbiased\" claim is credible. If that lands, this could be a useful contribution. I would bring it to a reading group and would engage with it as a referee.","headline":"A genuinely new shear estimator that shears the weight kernel rather than the image, with an honest but load-bearing gap: the noise-bias correction for the response R is still missing.","tokens_in":4719,"tokens_out":1507,"would_cite":false,"duration_ms":16579,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper develops a cosmic shear estimator that measures galaxy shapes from second moments and calibrates its response by shearing the weight-function coordinate system, achieving multiplicative bias below 10^-3 without assuming galaxy…","keywords":["cosmic shear","weak lensing","shear estimation","second moments","multiplicative bias","noise bias","self-calibration","galaxy images"],"falsifier":"Run the estimator on noise-free simulated images with a compact, high-frequency galaxy profile and a Moffat PSF at various pixel scales; if the recovered shear differs from the input by more than roughly $10^{-3}\\,g$ for shear values up to 0.03, the coordinate-shear equivalence is violated. Alternatively, compare the response matrix $R$ computed by the coordinate-shear method with $R$ computed by actually shearing the simulated image using the same pipeline; any statistically significant difference is a direct failure of Eqs. (5-6).","tokens_in":3756,"feed_emoji":"🔭","tokens_out":6928,"duration_ms":64396,"temperature":0.7,"pith_summary":"This paper develops a cosmic shear estimator that measures galaxy shapes from second moments and calibrates itself by shearing the coordinate system of the weight function rather than the galaxy image itself. The estimator requires no assumption about galaxy profiles and no suite of simulated galaxy images for calibration. On noise-free images of elliptical and COSMOS-like galaxies the method recovers input shears with multiplicative bias below $10^{-3}$, the level projected for future surveys. The paper also derives an analytic correction for the leading noise bias, caused by the estimator's nonlinear dependence on the measured centroid and pixel noise, and validates it against noisy simulations. If the self-calibration holds, shear calibration would no longer depend on the realism of galaxy simulations.","feed_headline":"Cosmic shear estimator hits 1e-3 bias without image warping","feed_subtitle":"By shearing coordinates instead of galaxy pixels and correcting noise analytically, the estimator meets survey-grade requirements.","key_machinery":"The load-bearing object is the PSF-convolved weight kernel $F(X) = ([XX^T W] \\ast \\psi)(X)$. Instead of shearing the galaxy image, the method shears the coordinates of $F$ by $S^{-1}$, divides by the PSF in Fourier space to obtain $G(S,X)$, and uses the moments $\\int G(S,X)I(X)\\,d^2X$ with $S$ chosen as $\\pm\\epsilon$ in each shear component. The derivatives of the ellipticity components with respect to $\\epsilon$ define the self-calibration response $R$, and the ratio $\\langle R\\rangle^{-1}\\langle e\\rangle$ cancels shape noise and PSF and pixel effects to first order. A second piece of machinery is the analytic noise-bias expansion: because the measured centroid is re-injected into the weight function, the moment estimator is nonlinear in the pixel noise, and the bias is computed from second derivatives of the moments with respect to image values and centroid position.","core_discovery":"The central claim is that one can measure gravitational shear without ever distorting the observed image and without knowing the galaxy's light profile. The ellipticity $e$ is computed from weighted second moments of the image; a response matrix $R$ is obtained by numerically applying small shears to the coordinate system of the PSF-convolved weight kernel $F$, recovering the image by Fourier division by the PSF, and taking finite differences of the resulting moments. The estimator $\\langle \\mathbf{g}\\rangle = \\langle R\\rangle^{-1}\\langle e\\rangle$ is then unbiased by construction with respect to the galaxy profile, up to a high-order cross-term between shear and pixel sampling that is subtracted analytically. In the noise-free simulations the multiplicative bias is below $10^{-3}$. The paper further computes the second-order noise bias analytically and shows it removes most of the ellipticity bias in noisy images; correcting the response matrix $R$ for noise is left as work in progress.","pith_inferences":["A testable extension is to compare the response $R$ computed by shearing coordinates with $R$ obtained by shearing the actual image in a Metacalibration-style pipeline for the same galaxies; agreement to the claimed tolerance would confirm the equivalence assumption in a regime the paper has not yet fully tested.","The analytic noise correction could be generalised to the moments used in the response matrix, which the paper leaves for future work; if successful, the method would no longer need noisy simulations for any bias term.","The coordinate-shear trick might extend to higher-order shape measurements such as flexion, where response matrices are harder to simulate, since the same Fourier-division trick applies.","Because the method avoids galaxy-profile assumptions, it could be combined with machine-learning shape classifiers as a cross-check, using the classifier only for selection, not calibration."],"forward_implications":["Multiplicative bias can be kept below $10^{-3}$ without generating a large library of galaxy shape simulations for calibration.","The same estimator can be applied to undersampled images, because the shear is applied to the more extended kernel $F$ rather than to the compact galaxy image.","The analytical noise-bias formula allows fast corrections on individual exposures, removing the need to simulate noise to calibrate noise bias.","Since no image distortion is applied, correlated noise from resampling the galaxy image is avoided.","If the noise-bias correction is extended from ellipticity to the response matrix $R$, the full shear estimator would be self-contained and applicable to survey data."],"supporting_citations":[{"why":"LSST Dark Energy Science Collaboration forecast establishing the precision target and motivating the need for shear calibration.","marker":"1"},{"why":"Huterer et al. bias parameterization used to define multiplicative and additive bias $m$ and $c$.","marker":"2"},{"why":"Cropper et al. set the $10^{-3}$ multiplicative bias requirement.","marker":"3"},{"why":"Sheldon and Huff (Metacalibration) provide the method this paper contrasts by not distorting the image.","marker":"4"},{"why":"Rowe et al. (GalSim) provide the simulation package used for all performance tests.","marker":"5"}],"fun_headline_variants":["Coordinate shear, no image warp: unbiased cosmic shear","No shape assumptions: cosmic shear bias below 1e-3","Shear coordinates, not pixels, for unbiased cosmic shear","Analytic noise bias correction for unbiased cosmic shear","Cosmic shear without warping: bias below 1e-3"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central bet is that shearing the coordinate system of the PSF-convolved weight kernel is exactly equivalent to shearing the galaxy image itself, once the PSF is divided out in Fourier space and the distorted pixel grid is interpolated.","fun_headline_variants_meta":{"raw":{"variants":["Coordinate shear, no image warp: unbiased cosmic shear","No shape assumptions: cosmic shear bias below 1e-3","Shear coordinates, not pixels, for unbiased cosmic shear","Analytic noise bias correction for unbiased cosmic shear","Cosmic shear without warping: bias below 1e-3"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001211,"raw_usage":{"total_tokens":4952,"prompt_tokens":880,"completion_tokens":4072,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":3989}},"tokens_in":496,"tokens_out":4072,"duration_ms":30959,"temperature":1.0,"reasoning_tokens":3989,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:31:42.682810+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the estimator on noise-free simulated images with a compact, high-frequency galaxy profile and a Moffat PSF at various pixel scales; if the recovered shear differs from the input by more than roughly $10^{-3}\\,g$ for shear values up to 0.03, the coordinate-shear equivalence is violated. Alternatively, compare the response matrix $R$ computed by the coordinate-shear method with $R$ computed by actually shearing the simulated image using the same pipeline; any statistically significant difference is a direct failure of Eqs. (5-6).","supporting_citations":[{"cited_title":"Cropper et al , DOI 10.1093, 2013","cited_arxiv_id":null,"evidence_quote":"Cropper et al. set the $10^{-3}$ multiplicative bias requirement."},{"cited_title":"Huterer et al , arXiv 0506030, 2006","cited_arxiv_id":null,"evidence_quote":"Huterer et al. bias parameterization used to define multiplicative and additive bias $m$ and $c$."}],"review_version":1}