pith. sign in

arxiv: 2605.16186 · v1 · pith:XRHHJOZTnew · submitted 2026-05-15 · 🌌 astro-ph.CO

Analytical method for computing the covariance matrix of cosmic shear two-point correlation function

Kosuke Nagura , Ryo Terasawa , Taisei Terawaki , Masahiro Takada This is my paper

Pith reviewed 2026-05-20 16:06 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords cosmic shearcovariance matrixtwo-point correlation functionweak lensingsurvey geometryharmonic spaceLegendre transformationGaussian covariance
0
0 comments X

The pith

Residual inaccuracies in the off-diagonal elements of harmonic-space covariance propagate through Legendre transforms and spoil real-space covariance estimates for cosmic shear under survey masks.

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

The paper tests analytical routes to the Gaussian covariance of the cosmic shear two-point correlation function that incorporate finite survey geometry. It applies the improved Narrow Kernel Approximation to the harmonic-space covariance and projects the result to real space with a Legendre transformation. Direct comparison with Gaussian simulations that include the HSC Year 3 mask shows that the resulting real-space covariance does not match the measured one. The mismatch is traced to small errors in the off-diagonal terms of the harmonic covariance, which are amplified by the transformation. A weighted quartic-counts estimator matches the simulations more closely, indicating that faithful modeling of those off-diagonal terms is required for trustworthy real-space covariance in the presence of survey windows.

Core claim

The covariance of the real-space two-point correlation function obtained from the improved Narrow Kernel Approximation does not reproduce the covariance measured directly from Gaussian simulations that include the survey mask. Although the approximation reproduces the diagonal structure of the harmonic-space covariance accurately, residual inaccuracies in the off-diagonal components propagate through the Legendre transformation and significantly affect the real-space result. In contrast, the weighted quartic-counts method shows better agreement with the simulations.

What carries the argument

Improved Narrow Kernel Approximation (iNKA) for the harmonic-space covariance, followed by projection via the Legendre transformation to real space.

If this is right

  • Reliable analytical covariance for real-space cosmic shear two-point functions requires accurate modeling of off-diagonal terms in harmonic space when survey windows are present.
  • The weighted quartic-counts method provides a better match to simulation covariances than the iNKA-plus-Legendre route under the same mask.
  • Simple f_sky approximations are likely to inherit similar limitations from incomplete treatment of survey geometry.
  • Covariance estimates used in cosmological analyses of current and future weak-lensing surveys must incorporate faithful off-diagonal structure to avoid biased parameter constraints.

Where Pith is reading between the lines

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

  • Surveys with more complex or irregular masks may amplify the propagation of off-diagonal errors, making direct simulation-based covariances preferable until analytical models improve.
  • Extending the comparison to include non-Gaussian covariance contributions could test whether the Gaussian off-diagonal mismatch remains the leading error source at smaller scales.
  • The result suggests that hybrid approaches combining analytical harmonic-space modeling with targeted simulation corrections for off-diagonals could be developed for upcoming wide-area surveys.

Load-bearing premise

The discrepancies between the iNKA-derived covariance and the simulation covariance are caused primarily by residual inaccuracies in the off-diagonal components of the harmonic-space covariance rather than by details of the Legendre transformation, the specific survey mask, or non-Gaussian contributions.

What would settle it

Directly comparing the off-diagonal elements of the harmonic-space covariance computed by iNKA against the same elements measured from the Gaussian simulations would show whether those off-diagonal inaccuracies are the dominant source of the mismatch.

read the original abstract

Accurate estimation of the covariance matrix of cosmic shear statistics is essential for cosmological analyses using current and upcoming wide-area weak lensing surveys. In this work, we investigate analytical methods for computing the Gaussian covariance matrix of the cosmic shear two-point correlation function (2PCF), taking into account the effects of finite survey geometry. We compute the covariance of 2PCF based on the improved Narrow Kernel Approximation (iNKA), with a projection using the Legendre transformation. We also consider other analytical covariance estimators, the $f_{\mathrm{sky}}$ approximation and the weighted quartic-counts method. We evaluate the accuracy of those analytical methods using the convergence fields with the HSC Year 3 survey mask as a test case. We find that the covariance of the 2PCF obtained by using the iNKA does not reproduce the covariance measured directly from Gaussian simulations. Although the iNKA accurately models the diagonal structure of the harmonic-space covariance, residual inaccuracies in the off-diagonal components propagate through the Legendre transformation and significantly affect the real-space covariance. In contrast, the weighted quartic-counts method shows better agreement with the simulations. Our results demonstrate that accurate modeling of the off-diagonal structure of the harmonic-space covariance is crucial for obtaining reliable covariance estimates of real-space weak lensing statistics in the presence of survey window effects.

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 / 1 minor

Summary. The paper develops analytical methods for the Gaussian covariance matrix of the cosmic shear two-point correlation function (2PCF) that incorporate finite survey geometry effects. It computes the covariance via the improved Narrow Kernel Approximation (iNKA) followed by a Legendre transformation, and compares this to the f_sky approximation and the weighted quartic-counts estimator. Accuracy is assessed by comparing the resulting real-space 2PCF covariances against direct measurements from Gaussian convergence simulations that employ the HSC Year 3 survey mask. The central conclusion is that iNKA fails to reproduce the simulation covariances because, while it captures the diagonal of the harmonic-space covariance, residual inaccuracies in the off-diagonals propagate through the Legendre transform; the weighted quartic-counts method agrees better with the simulations, demonstrating that accurate off-diagonal modeling in harmonic space is essential for reliable real-space covariance estimates.

Significance. If the reported discrepancies and their attribution hold, the work is significant for weak-lensing cosmology: it identifies a concrete limitation of a widely used analytical covariance estimator under realistic survey masks and supplies an external benchmark from independent Gaussian simulations. This directly informs covariance modeling choices for ongoing and future surveys (HSC, LSST, Euclid) where accurate covariance matrices are required for cosmological parameter inference.

major comments (2)
  1. [Abstract and results section] Abstract and results section: the statement that 'the covariance of the 2PCF obtained by using the iNKA does not reproduce the covariance measured directly from Gaussian simulations' is presented without quantitative metrics (e.g., Frobenius norm of the difference matrix, element-wise fractional residuals, or reduced chi-squared values with error bars). This leaves the central claim only partially supported.
  2. The attribution of the observed discrepancies to 'residual inaccuracies in the off-diagonal components' of the harmonic-space covariance (Abstract) is load-bearing for the conclusion that off-diagonal modeling is crucial. However, the manuscript does not provide an explicit side-by-side comparison of the full harmonic-space covariance matrices (iNKA versus simulation) that isolates off-diagonal residuals from possible differences in mask convolution, binning, or the numerical implementation of the Legendre projection itself.
minor comments (1)
  1. Specify the precise binning scheme and multipole range used when constructing the harmonic-space covariance before the Legendre transform.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and insightful comments on our manuscript. We appreciate the opportunity to clarify and strengthen our presentation of the results. Below, we address each major comment in detail.

read point-by-point responses

  1. Referee: [Abstract and results section] Abstract and results section: the statement that 'the covariance of the 2PCF obtained by using the iNKA does not reproduce the covariance measured directly from Gaussian simulations' is presented without quantitative metrics (e.g., Frobenius norm of the difference matrix, element-wise fractional residuals, or reduced chi-squared values with error bars). This leaves the central claim only partially supported.

    Authors: We agree with the referee that including quantitative metrics would provide stronger support for our central claim. In the revised manuscript, we have added quantitative comparisons, including the Frobenius norm of the difference between the analytical and simulated covariance matrices, as well as average element-wise fractional residuals. We have also computed a reduced chi-squared statistic for the agreement between the matrices, accounting for the simulation variance where possible. These additions are now included in the results section and referenced in the abstract. revision: yes

  2. Referee: [—] The attribution of the observed discrepancies to 'residual inaccuracies in the off-diagonal components' of the harmonic-space covariance (Abstract) is load-bearing for the conclusion that off-diagonal modeling is crucial. However, the manuscript does not provide an explicit side-by-side comparison of the full harmonic-space covariance matrices (iNKA versus simulation) that isolates off-diagonal residuals from possible differences in mask convolution, binning, or the numerical implementation of the Legendre projection itself.

    Authors: We acknowledge that a direct visual or quantitative comparison of the full harmonic-space covariance matrices would help to more clearly isolate the off-diagonal effects. In our original analysis, we verified that the iNKA reproduces the diagonal elements of the harmonic-space covariance accurately (see Section 3.2 and associated figures), and the discrepancies appear only after applying the Legendre transformation to obtain the real-space 2PCF covariance. To address the referee's point, we have included in the revision a new figure comparing the harmonic-space covariance matrices from iNKA and the simulations, with a focus on the off-diagonal elements. We used identical mask convolution, binning, and projection methods for both to minimize other sources of difference. This supports our attribution while providing the requested explicit comparison. revision: yes

Circularity Check

0 steps flagged

No significant circularity; external Gaussian simulation benchmarks used for validation

full rationale

The paper derives analytical covariance estimators (iNKA with Legendre projection, f_sky, weighted quartic-counts) and directly compares their outputs to covariance matrices measured from independent Gaussian simulations using the HSC Y3 mask. This supplies an external benchmark rather than reducing any result to a fitted parameter or self-referential equation. The iNKA method is referenced from prior literature, but the central accuracy claims rest on the simulation comparisons, not on self-citation chains. No self-definitional, fitted-input-as-prediction, or uniqueness-imported patterns appear in the derivation or evaluation steps.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the Gaussian random field assumption for the convergence maps used in simulations and on the validity of the Narrow Kernel Approximation framework itself.

axioms (1)
  • domain assumption The convergence field can be treated as a Gaussian random field when computing the covariance of the two-point correlation function.
    The paper evaluates analytical methods against covariance measured directly from Gaussian simulations.

pith-pipeline@v0.9.0 · 5779 in / 1374 out tokens · 77847 ms · 2026-05-20T16:06:02.735901+00:00 · methodology

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

30 extracted references · 30 canonical work pages · 15 internal anchors

  1. [2]

    X. Li, T. Zhang, S. Sugiyama, R. Dalal, R. Terasawa, M. M. Rau et al.,Hyper Suprime-Cam Year 3 results: Cosmology from cosmic shear two-point correlation functions,Phys. Rev. D108 (2023) 123518 [2304.00702]

    work page arXiv 2023

  2. [3]

    Dalal et al., Hyper Suprime-Cam Year 3 results: Cos- mology from cosmic shear power spectra, Phys

    R. Dalal, X. Li, A. Nicola, J. Zuntz, M. A. Strauss, S. Sugiyama et al.,Hyper Suprime-Cam Year 3 results: Cosmology from cosmic shear power spectra,Phys. Rev. D108(2023) 123519 [2304.00701]

    work page arXiv 2023

  3. [4]

    Terasawa, X

    R. Terasawa, X. Li, M. Takada, T. Nishimichi, S. Tanaka, S. Sugiyama et al.,Exploring the baryonic effect signature in the Hyper Suprime-Cam Year 3 cosmic shear two-point correlations on small scales: The S8 tension remains present,Phys. Rev. D111(2025) 063509 [2403.20323]

    work page arXiv 2025

  4. [5]

    The contribution from small scales on two-point shear analysis: comparison between power spectrum and correlation function

    J. Ferri, E. G. M. Ferreira and R. Terasawa,The contribution from small scales on two-point shear analysis: comparison between power spectrum and correlation function,arXiv e-prints (2025) arXiv:2512.17022 [2512.17022]

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  5. [6]

    Takahashi, T

    R. Takahashi, T. Hamana, M. Shirasaki, T. Namikawa, T. Nishimichi, K. Osato et al.,Full-sky gravitational lensing simulation for large-area galaxy surveys and cosmic microwave background experiments,The Astrophysical Journal850(2017) 24

    work page 2017

  6. [7]

    Power Spectrum Super-Sample Covariance

    M. Takada and W. Hu,Power spectrum super-sample covariance,Phys. Rev. D87(2013) 123504 [1302.6994]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  7. [8]

    Terasawa, R

    R. Terasawa, R. Takahashi, T. Nishimichi and M. Takada,Separate universe approach to evaluate nonlinear matter power spectrum for nonflatΛCDM model,Phys. Rev. D106(2022) 083504 [2205.10339]

    work page arXiv 2022

  8. [9]

    Terasawa, R

    R. Terasawa, R. Takahashi, T. Nishimichi and M. Takada,Flat to nonflat: Calculating nonlinear power spectra of biased tracers for a nonflatΛCDM model,Phys. Rev. D109(2024) 063504 [2310.13330]

    work page arXiv 2024

  9. [10]

    M. Sato, M. Takada, T. Hamana and T. Matsubara,Simulations of Wide-field Weak-lensing Surveys. II. Covariance Matrix of Real-space Correlation Functions,Astrophys. J.734(2011) 76 [1009.2558]

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  10. [11]

    CosmoLike - Cosmological Likelihood Analyses for Photometric Galaxy Surveys

    E. Krause and T. Eifler,cosmolike - cosmological likelihood analyses for photometric galaxy surveys,Mon. Not. Roy. Astron. Soc.470(2017) 2100 [1601.05779]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  11. [12]

    Reischke, S

    R. Reischke, S. Unruh, M. Asgari, A. Dvornik, H. Hildebrandt, B. Joachimi et al., KiDS-Legacy: Covariance validation and the unified ONECOVARIANCE framework for projected large-scale structure observables,Astronomy & Astrophysics699(2025) A124 [2410.06962]. – 10 –

    work page arXiv 2025

  12. [13]

    Mock galaxy shape catalogs in the Subaru Hyper Suprime-Cam Survey

    M. Shirasaki, T. Hamana, M. Takada, R. Takahashi and H. Miyatake,Mock galaxy shape catalogues in the Subaru Hyper Suprime-Cam Survey,Mon. Not. Roy. Astron. Soc.486(2019) 52 [1901.09488]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  13. [14]

    Friedrich et al.,Dark energy survey year 3 results: Covariance modelling and its impact on parameter estimation and quality of fit,2012.08568

    O. Friedrich et al.,Dark energy survey year 3 results: Covariance modelling and its impact on parameter estimation and quality of fit,2012.08568

    work page arXiv 2012

  14. [15]

    Performance of internal Covariance Estimators for Cosmic Shear Correlation Functions

    O. Friedrich, S. Seitz, T. F. Eifler and D. Gruen,Performance of internal covariance estimators for cosmic shear correlation functions,Mon. Not. Roy. Astron. Soc.456(2016) 2662 [1508.00895]

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  15. [16]

    Analysis of two-point statistics of cosmic shear: III. Covariances of shear measures made easy

    B. Joachimi, P. Schneider and T. Eifler,Analysis of two-point statistics of cosmic shear. III. Covariances of shear measures made easy,Astronomy & Astrophysics477(2008) 43 [0708.0387]

    work page internal anchor Pith review Pith/arXiv arXiv 2008

  16. [17]

    Cosmological Simulations for Combined-Probe Analyses: Covariance and Neighbour-Exclusion Bias

    J. Harnois-Déraps, A. Amon, A. Choi, V. Demchenko, C. Heymans, A. Kannawadi et al., Cosmological simulations for combined-probe analyses: covariance and neighbour-exclusion bias,Mon. Not. Roy. Astron. Soc.481(2018) 1337 [1805.04511]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  17. [18]

    X. Fang, T. Eifler and E. Krause,2D-FFTLog: efficient computation of real-space covariance matrices for galaxy clustering and weak lensing,Mon. Not. Roy. Astron. Soc.497(2020) 2699 [2004.04833]

    work page arXiv 2020

  18. [19]

    Truttero, J

    O. Truttero, J. Zuntz, A. Pourtsidou and N. Robertson,Baryon-free S 8 tension with Stage IV cosmic shear surveys,The Open Journal of Astrophysics8(2025) 19 [2410.18191]

    work page arXiv 2025

  19. [20]

    DeRose, N

    J. DeRose, N. Weaverdyck, M. White, S.-F. Chen, D. Schlegel and A. Slosar,Steeling Weak Lensing Source Galaxy Samples against Systematics using Wide Field Spectroscopy,arXiv e-prints(2026) arXiv:2603.10113 [2603.10113]

    work page arXiv 2026

  20. [21]

    Myths and Truths Concerning Estimation of Power Spectra

    G. Efstathiou,Myths and truths concerning estimation of power spectra: the case for a hybrid estimator,Mon. Not. Roy. Astron. Soc.349(2004) 603 [astro-ph/0307515]

    work page internal anchor Pith review Pith/arXiv arXiv 2004

  21. [22]

    García-García, D

    C. García-García, D. Alonso and E. Bellini,Disconnected pseudo-Cl covariances for projected large-scale structure data,JCAP2019(2019) 043 [1906.11765]

    work page arXiv 2019

  22. [23]

    Nicola, C

    A. Nicola, C. García-García, D. Alonso, J. Dunkley, P. G. Ferreira, A. Slosar et al.,Cosmic shear power spectra in practice,JCAP2021(2021) 067 [2010.09717]

    work page arXiv 2021

  23. [24]

    M. A. Troxel, E. Krause, C. Chang, T. F. Eifler, O. Friedrich, D. Gruen et al.,Survey geometry and the internal consistency of recent cosmic shear measurements,Mon. Not. Roy. Astron. Soc. 479(2018) 4998 [1804.10663]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  24. [25]

    Analysis of two-point statistics of cosmic shear: I. Estimators and covariances

    P. Schneider, L. van Waerbeke, M. Kilbinger and Y. Mellier,Analysis of two-point statistics of cosmic shear. I. Estimators and covariances,Astronomy & Astrophysics396(2002) 1 [astro-ph/0206182]

    work page internal anchor Pith review Pith/arXiv arXiv 2002

  25. [26]

    Shear Power Spectrum Reconstruction using Pseudo-Spectrum Method

    C. Hikage, M. Takada, T. Hamana and D. Spergel,Shear power spectrum reconstruction using the pseudo-spectrum method,Mon. Not. Roy. Astron. Soc.412(2011) 65 [1004.3542]

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  26. [27]

    A unified pseudo-$C_\ell$ framework

    D. Alonso, J. Sanchez, A. Slosar and LSST Dark Energy Science Collaboration,A unified pseudo-Cℓ framework,Mon. Not. Roy. Astron. Soc.484(2019) 4127 [1809.09603]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  27. [28]

    N. E. Chisari, D. Alonso, E. Krause, C. D. Leonard, P. Bull, J. Neveu et al.,Core Cosmology Library: Precision Cosmological Predictions for LSST,Astrophys. J. Suppl.242(2019) 2 [1812.05995]

    work page arXiv 2019

  28. [29]

    K. M. Górski, E. Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke et al., HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere,Astrophys. J.622(2005) 759 [astro-ph/0409513]

    work page internal anchor Pith review Pith/arXiv arXiv 2005

  29. [30]

    X. Li, H. Miyatake, W. Luo, S. More, M. Oguri, T. Hamana et al.,The three-year shear catalog – 11 – of the Subaru Hyper Suprime-Cam SSP Survey,Publ. Astron. Soc. Japan74(2022) 421 [2107.00136]

    work page arXiv 2022

  30. [31]

    The Skewness of the Aperture Mass Statistic

    M. Jarvis, G. Bernstein and B. Jain,The skewness of the aperture mass statistic, Mon. Not. Roy. Astron. Soc.352(2004) 338 [astro-ph/0307393]. – 12 –

    work page internal anchor Pith review Pith/arXiv arXiv 2004