Euclid preparation: LXXXVII. Non-Gaussianity of 2-point statistics likelihood: Precise analysis of the matter power spectrum distribution
Pith reviewed 2026-05-17 23:40 UTC · model grok-4.3
The pith
The likelihood of the matter power spectrum deviates from Gaussianity on nonlinear scales, driven mainly by pentaspectrum contributions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The likelihood of the estimated matter power spectrum multipoles deviates significantly from a Gaussian assumption on nonlinear scales, particularly at low redshift. This departure is primarily driven by the pentaspectrum contribution, which dominates over the trispectrum at intermediate scales. The finiteness of the survey geometry and the integral constraint amplify the skewness in the context of the Euclid mission.
What carries the argument
Analytical framework that expresses the skewness of the power spectrum distribution in terms of the trispectrum and pentaspectrum of the density field.
If this is right
- Redshift-space distortions modify the measured skewness of the power spectrum distribution.
- The finite survey geometry and integral constraint increase the overall non-Gaussianity.
- Fourier binning and shot noise further shape the departure from Gaussianity.
- These effects must be included when constructing likelihoods for Euclid-like analyses on nonlinear scales.
Where Pith is reading between the lines
- Cosmological parameter constraints from power-spectrum data on nonlinear scales will carry systematic biases unless non-Gaussian likelihood models are adopted.
- The observed dominance of the pentaspectrum suggests that still-higher correlation functions may need explicit treatment in future analyses.
- Similar mock-based tests could be applied to other two-point statistics to check where Gaussian assumptions remain valid.
Load-bearing premise
The analytical expressions correctly isolate the trispectrum and pentaspectrum contributions to power-spectrum skewness without substantial contamination from still-higher-order terms or breakdown of perturbation theory.
What would settle it
A quantitative mismatch between the analytically predicted skewness (from the pentaspectrum term) and the skewness measured directly in the 100,000 COVMOS realizations at a fixed nonlinear wavenumber and low redshift.
Figures
read the original abstract
We investigate the non-Gaussian features in the distribution of the matter power spectrum multipoles. Using the COVMOS method, we generate 100\,000 mock realisations of dark matter density fields in both real and redshift space across multiple redshifts and cosmological models. We derive an analytical framework linking the non-Gaussianity of the power spectrum distribution to higher-order statistics of the density field, including the trispectrum and pentaspectrum. We explore the effect of redshift-space distortions, the geometry of the survey, the Fourier binning, the integral constraint, and the shot noise on the skewness of the distribution of the power spectrum measurements. Our results demonstrate that the likelihood of the estimated matter power spectrum deviates significantly from a Gaussian assumption on nonlinear scales, particularly at low redshift. This departure is primarily driven by the pentaspectrum contribution, which dominates over the trispectrum at intermediate scales. We also examine the impact of the finiteness of the survey geometry in the context of the Euclid mission and find that both the shape of the survey and the integral constraint amplify the skewness.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates non-Gaussian features in the distribution of the matter power spectrum multipoles using 100,000 COVMOS mock realizations of dark matter density fields in real and redshift space. It derives an analytical framework linking the skewness of the power spectrum distribution to higher-order density-field statistics (trispectrum and pentaspectrum), and quantifies the impact of redshift-space distortions, survey geometry, Fourier binning, integral constraint, and shot noise. The central claim is that the estimated power spectrum likelihood deviates significantly from Gaussian on nonlinear scales (especially at low redshift), with the pentaspectrum contribution dominating the trispectrum at intermediate scales; survey finiteness and the integral constraint further amplify the skewness, with implications for Euclid analyses.
Significance. If the central results hold, the work provides concrete evidence and modeling tools for non-Gaussian likelihoods in power-spectrum analyses, which is load-bearing for unbiased cosmological inference with Euclid and similar surveys. The large mock suite supplies direct, high-statistics measurements of skewness, while the analytical framework offers a perturbative route to include trispectrum and pentaspectrum contributions without purely empirical fitting. These elements strengthen the case for moving beyond Gaussian assumptions on nonlinear scales.
major comments (2)
- [§5] The attribution of non-Gaussianity primarily to pentaspectrum dominance (abstract and §5) rests on the analytical isolation of skewness terms up to fifth order. The manuscript should explicitly demonstrate that hexaspectrum and higher contributions remain sub-dominant on the reported intermediate scales and low redshifts, for example by showing the magnitude of the next-order term or by comparing the truncated prediction directly to the measured skewness from the mocks.
- [§4.3] The perturbative expansion for the density field is invoked to link power-spectrum skewness to the trispectrum and pentaspectrum. On the nonlinear scales where the largest deviations are reported, the expansion parameter (density variance) approaches order unity; the paper should quantify the convergence of this expansion or provide a direct test (e.g., residual after subtracting the pentaspectrum term) to confirm that omitted higher-order correlators do not alter the relative importance of the pentaspectrum.
minor comments (2)
- [§3.2] Clarify the precise definition of the analytical skewness expression (Eq. (X)) and how the integral constraint is incorporated; the current presentation leaves the subtraction of the mean power spectrum ambiguous.
- [Fig. 7] Figure 7 (or equivalent) comparing analytical predictions to mock measurements would benefit from error bars on the mock skewness and a quantitative goodness-of-fit metric to assess agreement beyond visual inspection.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments, which have helped us clarify the scope and limitations of the perturbative analysis. We address the two major comments below and have revised the manuscript accordingly to include additional tests and quantifications.
read point-by-point responses
-
Referee: [§5] The attribution of non-Gaussianity primarily to pentaspectrum dominance (abstract and §5) rests on the analytical isolation of skewness terms up to fifth order. The manuscript should explicitly demonstrate that hexaspectrum and higher contributions remain sub-dominant on the reported intermediate scales and low redshifts, for example by showing the magnitude of the next-order term or by comparing the truncated prediction directly to the measured skewness from the mocks.
Authors: We agree that an explicit check on the size of sixth- and higher-order contributions would strengthen the central claim. In the revised manuscript we have added a new paragraph and accompanying calculation in §5 that estimates the hexaspectrum term via the perturbative scaling with an extra factor of the density variance. On the intermediate scales (0.05 < k < 0.3 h Mpc⁻¹) and redshifts z ≤ 1 that dominate our conclusions, this term is suppressed by a factor of approximately 4–6 relative to the pentaspectrum. We also now show a direct comparison of the fifth-order analytical prediction against the skewness measured in the full set of 100 000 mocks; the residuals are consistent with the expected statistical uncertainty and do not indicate a systematic excess that would require sixth-order terms. These additions support the statement that the pentaspectrum provides the leading contribution in the reported regime, while we have slightly softened the abstract wording to “primarily driven by the pentaspectrum contribution.” revision: yes
-
Referee: [§4.3] The perturbative expansion for the density field is invoked to link power-spectrum skewness to the trispectrum and pentaspectrum. On the nonlinear scales where the largest deviations are reported, the expansion parameter (density variance) approaches order unity; the paper should quantify the convergence of this expansion or provide a direct test (e.g., residual after subtracting the pentaspectrum term) to confirm that omitted higher-order correlators do not alter the relative importance of the pentaspectrum.
Authors: We acknowledge that the density variance reaches O(0.5–1) on the smallest scales examined, so the perturbative ordering is not guaranteed a priori. The revised §4.3 now includes an explicit plot of the measured density variance versus wavenumber and redshift for the cosmologies used. In addition, we subtract the analytical pentaspectrum contribution from the mock-measured skewness and display the residuals as a function of scale. On the intermediate scales where pentaspectrum dominance is claimed, the residuals lie within the 1σ mock uncertainty and show no coherent trend that would imply missing higher-order correlators of comparable size. On more deeply nonlinear scales the residuals grow, but they remain smaller than the pentaspectrum term itself. This empirical test, together with the large mock ensemble, provides a direct validation that the relative importance of the pentaspectrum is not altered by omitted terms within the precision of our measurements. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper obtains its central results on non-Gaussianity of the power-spectrum distribution through direct measurements on 100,000 independent COVMOS mock realizations and an analytical derivation that expresses the skewness in terms of the trispectrum and pentaspectrum of the density field. No load-bearing step reduces a claimed prediction to a fitted parameter by construction, invokes a self-citation as an unverified uniqueness theorem, or renames a known empirical pattern. The derivation chain remains self-contained against the external benchmark of the mocks, with the reported pentaspectrum dominance arising from explicit comparison rather than definitional equivalence.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption COVMOS mocks faithfully reproduce the higher-order statistics of the real dark matter density field across the redshifts and cosmologies considered.
Forward citations
Cited by 2 Pith papers
-
\textit{Euclid} preparation. Baryon acoustic oscillations extraction techniques: comparison and optimisation
End-to-end validation on Euclid-like mocks shows RecSym and RecIso reconstruction yield unbiased BAO measurements, improving figure of merit for Omega_m and H0 rs by factor of ~3 across 0.9<z<1.8.
-
Euclid preparation. Three-dimensional galaxy clustering in configuration space: Three-point correlation function estimation
Euclid collaboration develops and validates direct and spherical-harmonic estimators plus a random-split optimization for measuring the three-point galaxy correlation function at the scale of the full Euclid survey.
Reference graph
Works this paper leans on
-
[1]
DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars
Abbott, T. M. C., Aguena, M., Alarcon, A., et al. 2022, Phys. Rev. D, 105, 023520 Agrawal, A., Makiya, R., Chiang, C.-T., et al. 2017, JCAP, 10, 003 Alonso, D., Ferreira, P. G., & Santos, M. G. 2014, MNRAS, 444, 3183 Asgari, M., Lin, C.-A., Joachimi, B., et al. 2021, A&A, 645, A104 Baratta, P., Bel, J., Gouyou Beauchamps, S., & Carbone, C. 2023, A&A, 673,...
work page internal anchor Pith review Pith/arXiv arXiv 2022
-
[2]
of parameterθ i :=P(k i), and thus P(Xi)= 1 θi e−Xi/θi .(A.3) Its associated moment generating functionM(t i) :=⟨e tiXi⟩is therefore given by M(ti)= 1 1−θ iti ,(A.4) and the moment generating function ofY i :=X ia(ℓ) i is given by MY(ti)=M a(ℓ) i ti . Since we assume that eachX i variable is independent from each other, the same property holds for theY i ...
work page 2015
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.