Pith. sign in

REVIEW 3 major objections 5 minor 21 references

Analytical Template for the 4-Point Correlation Function Covariance Beyond the Gaussian Random Field ${\rm II}$: 1-Loop Corrections with Third-Order Densities

T0 review · 3 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash

Pith's one-line read The paper claims that the 1-loop 4PCF covariance from third-order densities reduces to three universal contraction classes, each expressible as low-dimensional radial integrals over the linear power spectrum.

desk verdict A careful, useful template for 1-loop 4PCF covariance, with one load-bearing combinatorial claim (Appendix E) that needs proof and a numerical check before it's survey-ready. read the letter →

arxiv 2509.05422 v1 pith:2VPOZBPG submitted 2025-09-05 astro-ph.CO

classification astro-ph.CO
keywords 4-pointcorrelationfunctioncovariancematrixbeyondGaussianthird-orderdensitycontraststandardperturbationtheorygalaxybiasisotropicbasisfunctionsparityviolation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper aims to supply an analytical, non-Gaussian covariance template for the galaxy 4-point correlation function at next-to-leading order in perturbation theory. It shows that every Wick contraction of two third-order density fields falls into exactly three classes—secondary-secondary, primary-secondary, and primary-primary—and gives explicit formulas for each. The angular integrals are performed analytically, leaving only low-dimensional radial integrals over the linear power spectrum and angular-momentum coefficients. If correct, 4PCF error bars used in parity-violation and cosmological tests can be computed without brute-force simulation or high-dimensional Monte Carlo integration.

What carries the argument

The load-bearing object is the generalized third-order density kernel W^(3), which repackages the standard third-order perturbation kernels and the third-order bias terms into one compact form. Around it, the argument organizes all Wick contractions into three loop configurations based on which vertices carry the third-order density. Angular integrals are performed analytically with the isotropic basis functions and rotation-averaging identities, reducing the problem to one-dimensional radial integrals over the linear power spectrum.

What would settle it

Numerically integrate Eq. (E.3) including all partially-connected terms with a realistic linear power spectrum; if any partially-connected integral contributes at the advertised r_c^3/V order, the template misses same-order terms. Alternatively, estimate the covariance from mock catalogs at scales near 100 Mpc/h and compare the diagonal with the template's prediction.

Watch

Extended reading notes

Core claim

The central claim is that the 1-loop (next-to-leading-order) covariance of the 4-point correlation function, built from third-order density fluctuations and third-order galaxy bias, can be written as the sum of three configuration classes. All four physical contributions—δ^(3), δ_lin G^(2), G^(3), and Γ^(3)—are unified through a single generalized kernel W^(3), so a single Wick-contraction procedure covers all of them. The explicit results, Eqs. (3.20), (3.30), and (3.40), express each class through radial integrals g, h, f, and ξ, with all angular structure captured by tabulated coefficients. Starting from the general covariance expression, one applies Wick's theorem, matches the contractio

Load-bearing premise

The argument assumes the covariance can be ordered by powers of r_c^3/V and that all partially-connected Wick contractions cancel against disconnected pieces at 1-loop order, leaving only the fully-connected terms.

Editorial extensions

If this is right

  • A complete analytical 4PCF covariance from third-order densities is obtained simply by matching any Wick contraction to one of the three classes and substituting W^(3).
  • The same template handles δ^(3), δ_lin G^(2), G^(3), and Γ^(3) bias terms, so third-order bias enters the covariance systematically rather than case by case.
  • Parity-even and parity-odd 4PCF error estimates improve over the Gaussian assumption because mode-coupling corrections beyond the Gaussian random field are included.
  • Together with the companion second-order calculation, this completes the perturbative 1-loop 4PCF covariance template up to the connected terms.
  • The framework is structured for direct extension to redshift-space distortions and for validation against simulations and mock catalogs.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same three-class decomposition likely extends to N-point correlation functions with N>4, since the classification depends only on which external vertices carry the third-order density and the remaining steps are permutations of the same Wick contractions.
  • Precomputing the radial integrals g, h, f, and ξ as interpolation tables would make this covariance cheap enough to embed in survey pipelines and parameter scans, which the paper does not explicitly propose.
  • A direct numerical check of Eq. (E.9) against the full connected covariance including partially-connected terms would settle the power-counting premise; if those terms do not fully cancel, additional same-order contributions would enter.
  • Because the kernel is generalized, non-standard third-order kernels from effective models could be slotted in by changing coefficients, making the template usable beyond standard perturbation theory.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. This paper presents an analytical template for the 1-loop, beyond-Gaussian 4-point correlation function (4PCF) covariance, incorporating third-order matter densities and a third-order galaxy bias expansion. The authors introduce a generalized third-order kernel W^(3), apply Wick contractions, and classify all surviving terms into three topologies: secondary-secondary, primary-secondary, and primary-primary. Using the isotropic basis of Cahn & Slepian, they reduce the angular integrals to Gaunt-type coefficients and express the final result as low-dimensional radial integrals involving g, h, f, and xi, together with 72 permutations per topology. The final expressions are Eqs. (3.20), (3.30), and (3.40). The paper defers numerical validation to future work.

Significance. If the claimed completeness and the reduction to three contraction classes are correct, this is a valuable step toward survey-ready non-Gaussian covariance models for the 4PCF, with direct relevance to DESI and BOSS parity-violation analyses. The paper reuses and extends established mathematical machinery, and the use of a generalized kernel W^(3) is elegant: it separates the SPT kernel structure from the angular/radial evaluation, so that different third-order kernels can be inserted without redoing the full calculation. The work is symbolic and parameter-free in the sense that no fitting or tuning is performed; the bias parameters are inputs. However, the central claim that Eqs. (3.20), (3.30), and (3.40) exhaust the 1-loop third-order covariance is not yet fully secured because the cancellation of partially-connected Wick contractions in Appendix E is asserted rather than proved, and several explicit integral expressions contain index and argument inconsistencies that would need correction before implementation.

major comments (3)
  1. [Appendix E, Eqs. (E.5)-(E.9)] The completeness of the three-class template rests on the claim that 'all the partially-connected terms cancel' after subtracting the disconnected estimator. This is demonstrated for one exemplary permutation only; there is no general enumeration of the Wick contraction classes, no systematic power-counting per class, and no explicit mapping of the surviving 216 terms to the three classes used in the main text. Moreover, the retained patterns in Eqs. (3.1), (3.22), and (3.32) include same-tetrahedron contractions (e.g., ⟨δlin(k1)δlin(k2)⟩ in Eq. (3.1) and internal ⟨q_i δlin⟩ pairs in Eq. (3.32)); whether these are partially-connected terms that should cancel or surviving contributions is never demonstrated. This is a load-bearing gap and should be filled with a complete enumeration or a numerical Wick-contraction check.
  2. [Sec. 3.2, Eqs. (3.23), (3.25), (3.27), (3.31)] The primary-secondary configuration has the q3 and k3 attachment variables swapped relative to the Wick pattern. In Eq. (3.22) the contractions are ⟨q3 δlin(k′1)⟩ and ⟨k3 δlin(k′3)⟩, and Eq. (3.23) correspondingly contains e^{-iq3·(r0-r′1-s)} and e^{-ik3·(r3-r′3-s)}. However Eq. (3.25) evaluates I(q3) with radial argument r′3 while the angular basis uses br′1, and Eq. (3.27) evaluates I(k3) with r′1. The radial integral S in Eq. (3.31) then inherits this mismatch (second h uses r′3, final g uses r′1). If implemented literally, the final template in Eq. (3.30) would not correspond to Eq. (3.23). This needs to be corrected or explicitly explained.
  3. [Sec. 3.3, Eqs. (3.39) and (3.41)] In the primary-primary expression, the radial function accompanying the k0 integral is written as h^{[n′3-n]}_{L33,L′3s,L′33,ℓ2}(r0,s,r′0,r) in Eq. (3.39) and in the definition of S in Eq. (3.41). Since the expansion coefficient is C_{L00,L0s,L′00} and the angular integral is L^{j12,j13,j,ℓ2}_{L00,L0s,L′00}, the h indices should be L00, L0s, L′00, not L33, L′3s, L′33. As written, the expression is internally inconsistent and would prevent direct numerical evaluation.
minor comments (5)
  1. [Sec. 3.1-3.3] The '71 perms.' are never written out. The counting footnotes give 4×6×3, 4×3×3×2, and 3×4×6, but the actual permuted contractions are not displayed or enumerated. For reproducibility, an appendix listing the permutations or a computer-readable contraction list would be very helpful.
  2. [Sec. 3.2] The sentence 'We have already solved the integrals over q2, q3, k2, and k3 in §3.1, given by Eqs. (3.9), (3.9), (3.27), (3.27)' has incorrect cross-references; Eq. (3.9) is the q2 integral of the secondary-secondary configuration, not the q2 integral of the primary-secondary configuration, and the repeated Eq. (3.27) should likely be Eqs. (3.26) and (3.27).
  3. [Eq. (3.35)] In the primary-primary configuration, Eq. (3.35) writes I(k3)(r3,−s,−r′3) with g[0]_{L33,L′3s,L′33}(r3,s,r′3) and PL(br3,−bs,−br′3), which is consistent with Eq. (3.32). However the text above says the k1,k2,k3 integrals have 'different angular momenta dependence from those of Eq. (3.27)'; this comparison to Eq. (3.27) is confusing because Eq. (3.27) is in the primary-secondary section, not an earlier expression of the same integral.
  4. [Eq. (3.20)] The factor ordering and the definitions of the sums over 'All' are not spelled out; it would improve clarity to state explicitly that the sums run over all angular momenta, intermediate momenta, and the coefficients appearing in the generalized kernel, possibly with a notation table. The current presentation makes the formula hard to follow.
  5. [Appendix C] The definitions of the Q coefficients in Eqs. (C.8), (C.12), and (C.15) use a large number of magnetic quantum numbers with no explicit statement of which indices are summed and which are constrained by the Gaunt/Wigner symbols. A sentence clarifying the summation conventions would be useful.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the 1-loop 4PCF covariance template is an SPT/Wick expansion expressed through independent mathematical identities; the only flagged issue is an asserted (not fully enumerated) cancellation in Appendix E, which is a completeness risk, not circularity.

full rationale

The paper's derivation chain is not circular. Eq. (2.6) follows by inserting the Eulerian bias expansion and SPT expansion into the covariance definition Eqs. (2.1)-(2.2), and the Wick contractions lead to the three independent classes in Eqs. (3.20), (3.30), and (3.40). The generalized kernel W^(3) in Eq. (A.18) is a bookkeeping device that re-expresses known third-order SPT/Galilean/Gamma kernels; it is not defined in terms of the desired covariance. No parameter is fitted, and no empirical quantity is 'predicted' from itself; bias coefficients are symbolic inputs. The reused prior results ([3] isotropic basis, [13] radial integral identities, [8] covariance formalism) are stated mathematical identities or published formalism, not the target 1-loop third-order 4PCF covariance, so shared authorship with [3] and [13] does not make the argument circular. Explicit flag for the review rule: Appendix E asserts 'we find that all the partially-connected terms cancel' after demonstrating one permutation, without enumerating all Wick contractions or proving the power-counting cancellation generally; this is a completeness/rigor limitation of the claim that Eqs. (3.20)/(3.30)/(3.40) exhaust the 1-loop answer, but it is not a reduction of the result to its inputs by construction. Overall circularity score 1.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The template rests on standard SPT, bias expansion, and angular-momentum recoupling tools, plus a power-counting assumption about which Wick contractions survive. No parameters are fitted to data in this paper; the bias coefficients are symbolic inputs.

free parameters (4)
  • b1
    Linear Eulerian bias coefficient; input to the template, not fitted here.
  • b_deltaG2
    Second-order Galileon bias coefficient retained in Eq. (2.3); value set by external fits.
  • b_G3
    Third-order Galileon bias coefficient in Eq. (2.3); input to the template.
  • b_Gamma3
    Third-order Gamma bias coefficient in Eq. (2.3); input to the template.
assumptions (5)
  • domain assumption SPT expansion of the matter density contrast to third order, Eq. (2.4)
    The whole 1-loop covariance is built on this perturbative expansion; convergence in the mildly nonlinear regime is assumed.
  • standard math Wick's theorem applies to the linear density field
    Used repeatedly in Sec. 3 to evaluate ensemble averages of products of linear fields.
  • domain assumption Power-counting scaling r_c^3/V orders the covariance terms
    Appendix E uses this to justify cancellation of partially-connected terms; if wrong, the connected covariance is incomplete.
  • domain assumption Eulerian bias expansion with only {1, delta, G2, G3, Gamma3} at third order
    Eq. (2.3); higher-derivative and other bias terms are neglected.
  • standard math Completeness of the isotropic basis and associated PWE
    Sec. 3 and Appendix B rely on the basis of Cahn & Slepian (2023); assumed complete for expanding plane waves and products of spherical harmonics.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Analytical Template for the 4-Point Correlation Function Covariance Beyond the Gaussian Random Field ${\rm II}$: 1-Loop Corrections with Third-Order Densities." pith.science (2026). https://pith.science/paper/2VPOZBPG

@misc{pith2026250905422,
  author       = {Pith},
  title        = {Pith review of: Analytical Template for the 4-Point Correlation Function Covariance Beyond the Gaussian Random Field $\rm II$: 1-Loop Corrections with Third-Order Densities},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2VPOZBPG}},
  note         = {Machine review of arXiv:2509.05422}
}
read the original abstract

Analytical templates for the 4-Point Correlation Function (4PCF) covariance matrix have been developed in the past assuming a Gaussian Random Field (GRF). In this work, we present the second part of the beyond GRF calculation of the 4PCF covariance, incorporating 1-loop corrections stemming from the third-order density contrast. Furthermore, we introduce a non-trivial galaxy biasing scheme at third order. To simplify the calculation, we decompose the covariance into three distinct structures and leverage the isotropic basis of Cahn & Slepian (2023). This approach reduces the complexity of the high-dimensional integrals that would be naively involved, enabling the angular parts to be performed and leaving us with low-dimensional radial integrals. This analytical template will provide a more accurate characterization of the statistical errors on the 4PCF, improving our ability to probe both its parity-even and parity-odd modes. This is the second and final paper in a two-part series.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

21 extracted references · 16 canonical work pages

  1. [1]

    Large-scale structure of the universe and cosmological perturbation theory,

    Bernardeau F., et al., "Large-scale structure of the universe and cosmological perturbation theory," Physics Reports, Sept. 2002

  2. [2]

    Non-Gaussian covariance of the matter power spectrum in the effective field theory of large scale structure,

    Bertolini D., et al., "Non-Gaussian covariance of the matter power spectrum in the effective field theory of large scale structure," PRD, June 2016

  3. [3]

    Isotropic N-point basis functions and their properties,

    Cahn R. N., Slepian Z., "Isotropic N-point basis functions and their properties," J. Phys. A Math. Gen., Aug. 2023

  4. [4]

    Test for Cosmological Parity Violation Using the 3D Distribution of Galaxies,

    Cahn R. N., Slepian Z., Hou J., "Test for Cosmological Parity Violation Using the 3D Distribution of Galaxies," PRL, Volume 130, Issue 20, article id.201002, May 2023

  5. [5]

    Triple-spherical Bessel function integrals with exponential and Gaussian damping: towards an analytic N-point correlation function covariance model,

    Chellino J., Slepian Z., "Triple-spherical Bessel function integrals with exponential and Gaussian damping: towards an analytic N-point correlation function covariance model," PRSA, Vol. 479, Issue 2276, Art. id.20230138, Aug. 2023

  6. [6]

    DESI 2024 V: Full-Shape Galaxy Clustering from Galaxies and Quasars,

    DESI Collaboration, et al., "DESI 2024 V: Full-Shape Galaxy Clustering from Galaxies and Quasars," arXiv e-prints 2411.12021, Nov. 2024

  7. [7]

    http://dlmf.nist.gov/

    DLMF, NIST Digital Library of Mathematical Functions. http://dlmf.nist.gov/

  8. [8]

    Analytic Gaussian covariance matrices for galaxy N -point correlation functions,

    Hou J., et al., "Analytic Gaussian covariance matrices for galaxy N -point correlation functions," PRD, Aug. 2022

Show all 21 references
  1. [9]

    Hou J., et al., "Measurement of parity-odd modes in the large-scale 4-point correlation function of Sloan Digital Sky Survey Baryon Oscillation Spectroscopic Survey twelfth data release CMASS and LOWZ galaxies," MNRAS, May 2023

  2. [10]

    Study of the Connected Four-Point Correlation Function of Galaxies from DESI Data Release 1 Luminous Red Galaxy Sample,

    Hou J., et al., "Study of the Connected Four-Point Correlation Function of Galaxies from DESI Data Release 1 Luminous Red Galaxy Sample," arXiv:2508.09070, Aug. 2025

  3. [11]

    A Non-Degenerate Neutrino Mass Signature in the Galaxy Bispectrum,

    Kamalinejad F., Slepian Z., "A Non-Degenerate Neutrino Mass Signature in the Galaxy Bispectrum," arXiv e-prints 2011.00899, Nov. 2020. – 32 –

  4. [12]

    Perturbative approach to covariance matrix of the matter power spectrum,

    Mohammed I., et al., "Perturbative approach to covariance matrix of the matter power spectrum," MNRAS, April 2017

  5. [13]

    A Model for the Redshift-Space Galaxy 4-Point Correlation Function,

    Ortolá Leonard W., Slepian Z., Hou J., "A Model for the Redshift-Space Galaxy 4-Point Correlation Function," arXiv e-prints 2402.15510, Feb. 2024

  6. [14]

    Analytical Template for the 4-Point Correlation Function Covariance Beyond the Gaussian Random Field I: 1-Loop Corrections with Second-Order Densities,

    Ortolá Leonard W., Slepian Z., "Analytical Template for the 4-Point Correlation Function Covariance Beyond the Gaussian Random Field I: 1-Loop Corrections with Second-Order Densities,"In prep

  7. [15]

    A First Detection of the Connected 4-Point Correlation Function of Galaxies Using the BOSS CMASS Sample,

    Philcox, O.H.E., Hou J., Slepian Z., "A First Detection of the Connected 4-Point Correlation Function of Galaxies Using the BOSS CMASS Sample," eprint arXiv:2108.01670, Aug. 2021

  8. [16]

    Probing parity violation with the four-point correlation function of BOSS galaxies,

    Philcox, O.H.E., "Probing parity violation with the four-point correlation function of BOSS galaxies," PRD, Sept. 2022

  9. [17]

    The Bispectrum as a Signature of Gravitational Instability in Redshift Space,

    Scoccimarro R., et al., "The Bispectrum as a Signature of Gravitational Instability in Redshift Space," AJ, June 1999

  10. [18]

    On a generating function for the isotropic basis functions and other connected results,

    Slepian Z., Chellino J., Greco A., "On a generating function for the isotropic basis functions and other connected results," Journal of Physics A: Mathematical and Theoretical, Volume 57, Issue 50, id.505203, 32 pp., Dec. 2024

  11. [19]

    Modeling the large-scale redshift-space 3-point correlation function of galaxies,

    Slepian Z., Eisenstein D. J., "Modeling the large-scale redshift-space 3-point correlation function of galaxies," MNRAS, Volume 469, Issue 2, p.2059-2076, Aug. 2017

  12. [20]

    Measurement of Parity-Violating Modes of the Dark Energy Spectroscopic Instrument (DESI) Year 1 Luminous Red Galaxies’ 4-Point Correlation Function,

    Slepian Z., et al., "Measurement of Parity-Violating Modes of the Dark Energy Spectroscopic Instrument (DESI) Year 1 Luminous Red Galaxies’ 4-Point Correlation Function," arXiv:2508.09133, Aug. 2025

  13. [21]

    Galaxy power spectrum multipoles covariance in perturbation theory,

    Wadekar D., Scoccimarro R., "Galaxy power spectrum multipoles covariance in perturbation theory," PRD, Dec. 2020. – 33 –

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.