Pith. sign in

REVIEW 2 major objections 5 minor 2 cited by

One formula gives exact N-point functions of non-Gaussian fields, even when the local mapping is non-analytic.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-03 03:03 UTC pith:NR2BRN74

load-bearing objection Useful resummation for locally non-Gaussian fields, but the 'exact' N≥3 claim outruns the convergence proof—treat the expansion as asymptotic beyond two points. the 2 major comments →

arxiv 2602.10151 v2 pith:NR2BRN74 submitted 2026-02-09 astro-ph.CO cond-mat.stat-mechmath-phmath.MP

A non-perturbative framework for N-point functions of locally non-Gaussian fields

Hardi Veerm\"ae This is my paper
classification astro-ph.CO cond-mat.stat-mechmath-phmath.MP
keywords locally non-Gaussian fieldsKibble–Slepian formulaN-point correlation functionsnon-analytic mappingexponential-tail modelprimordial black holesscalar-induced gravitational wavespower spectrum IR tail
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper claims that every N-point function of a locally non-Gaussian field ζ(x)=F(ζ_G(x))—including non-analytic mappings—can be written exactly as a sum over powers of the Gaussian correlation function, with all non-Gaussian information packed into one-point coefficients C_s. The derivation uses the Kibble–Slepian decomposition of the multivariate Gaussian density, replacing the usual Taylor expansion of F with a semi-perturbative series in the Gaussian correlation function. For a model with exponential tails, the paper derives an exact strongly non-Gaussian two-point function and shows that the normalized correlation saturates while the power spectrum develops a k^3 infrared tail. The result matters because it extends analytic control to strongly non-Gaussian regimes relevant to primordial black holes and scalar-induced gravitational waves.

Core claim

The paper establishes the exact expansion (3.9): any N-point function of a locally non-Gaussian field is a sum over multiplicity matrices ν_ij of products (ξ_ij^{ν_ij}/ν_ij!) times s_i! C_{s_i}, where s_i is the total multiplicity at point i and C_s is a one-point Gaussian expectation value of the rescaled Hermite polynomial ar H_s with F. This remains meaningful when F(ζ_G) is not analytic, because C_s is defined as an integral rather than through a Taylor expansion. In the exponential-tail model, the paper evaluates the coefficients and derives the closed-form two-point function G_2 = ξ_0/(8β²) arcsin²(ξ_G/ξ_0) in the strong non-Gaussianity limit, from which the saturation of the normalize

What carries the argument

The Kibble–Slepian decomposition: an identity that expands the n-dimensional Gaussian probability density in products of one-dimensional Hermite-polynomial densities, with powers of the correlation matrix serving as expansion coefficients. Applied to the local field ζ=F(ζ_G), it replaces the multidimensional path integral with sums of products of ξ_ij times resummed one-point coefficients C_s. These coefficients absorb all information about F and can be computed by a one-dimensional Gaussian integral even when F is non-analytic.

Load-bearing premise

The load-bearing premise is that term-by-term integration of the Kibble–Slepian series against arbitrary non-analytic F is legitimate and that the resulting power series in the Gaussian correlation function represents the exact N-point function; the paper gives no convergence proof and notes the series is poorly behaved near full correlation.

What would settle it

Evaluate G_2(ξ_G/ξ_0) for the exponential-tail model at an intermediate β (say β√ξ_0 ≈ 1) by direct numerical double integration of the bivariate Gaussian average, and compare against truncated versions of Eq. (3.35) for |ξ_G/ξ_0| between 0.9 and 1; if the truncated series does not approach the exact integral as more terms are included, the claimed exactness of the expansion is falsified. Alternatively, test Eq. (4.9) directly at large β by comparing the closed-form arcsin² expression to Monte Carlo integration of the original bivariate average.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • For any locally non-Gaussian field, the N-point function is determined entirely by the Gaussian correlation function and the coefficients C_s, so the mapping G_n can be precomputed independently of the shape of the power spectrum.
  • The formalism extends to non-analytic F, so strongly non-Gaussian models with exponential tails, logarithms, or piecewise continuations can be treated without a Taylor expansion.
  • For the exponential-tail model, the exact strong-NG two-point function gives a normalized correlation ξ(x)/ξ(0) = (4/π²) arcsin²(ξ_G(x)/ξ_0), independent of β, and a power spectrum with a k^3 infrared tail.
  • The series in C_s captures nonlinearities where the conventional F_NL,n expansion fails; the paper estimates that perturbativity breaks down when ξ_0 f_NL² ≲ O(10^-2).

Where Pith is reading between the lines

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

  • If the exact expansion holds for non-analytic F, the standard practice of expanding curvature perturbations as a polynomial in the Gaussian field may be unnecessary for a whole class of models; this could change how strongly non-Gaussian SIGW and PBH constraints are computed.
  • The k^3 infrared tail is a generic prediction of local non-Gaussianity whenever G_2 is unrelated to ξ_G; this gives a target observable feature that could be searched for in induced gravitational-wave backgrounds.
  • The paper's method still leaves n≥3 practical: the parameter space of covariances grows rapidly, so fully non-perturbative bispectra and trispectra may need interpolation or compressed representations; the paper notes this as an obstacle.
  • One could test the formalism on a known analytic example (e.g., F quadratic) to benchmark truncation errors, then extend to F with a branch point; the paper does not perform such an explicit convergence study.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 5 minor

Summary. The paper develops a formalism for N-point functions of locally non-Gaussian fields ζ(x)=F(ζ_G(x)) without Taylor-expanding F. Using the Kibble–Slepian decomposition of the multivariate Gaussian distribution, Eq. (3.9) expresses the N-point function as a power series in the Gaussian correlation ξ_ij, with all F-dependence encoded in one-point coefficients C_s defined by Eq. (3.10). The authors give diagrammatic Feynman rules, explicit formulas for bispectra and trispectra, and apply the method to an exponential-tail model, deriving an exact strongly non-Gaussian two-point function, Eq. (4.9), and studying the resulting deformation of the power spectrum.

Significance. If the convergence question is settled, this is a valuable framework: it factorizes the model-dependent and spectrum-dependent parts of the computation, extends beyond analytic mappings F, and provides non-perturbative benchmarks for PBH and SIGW applications. The central identity is derived from a rigorous mathematical formula; the coefficients C_s are computed from F rather than fitted to the target correlation function; numerical checks for N=2 are reported; and the limitations of the exponential-tail model are stated explicitly. The main unresolved issue is whether the multivariate series (3.9) is an exact convergent representation for N≥3 and non-analytic F; the paper currently provides no proof of this.

major comments (2)
  1. [Sec. 3.1, Eq. (3.9)] The derivation of Eq. (3.9) integrates the Kibble–Slepian series term by term. For N=2, Eq. (3.35) is absolutely convergent for |ψ|≤1 because ∑ n!C_n²=⟨F²⟩<∞. For N≥3 no analogous bound is given: the terms contain products ∏ s_i!C_{s_i}; writing C_s=a_s/√(s!) with ∑a_s²<∞, the factors √(s_i)!|a_{s_i}| can grow. A saddle-point estimate for the exponential-tail model gives dominant balanced diagrams scaling as (√(2ψ))^S S^{-3/4}, suggesting divergence when ψ>1/2. Thus exactness of (3.9) for non-analytic F is not established; the numerical checks in Sec. 4 cover only the two-point function. Please either prove convergence under stated assumptions or explicitly qualify (3.9) as asymptotic/formal for N≥3 and provide practical convergence criteria.
  2. [Sec. 4, footnote 11 / Eq. (4.9)] The exact strong-NG result Eq. (4.9) is derived in the footnote only as a 'rough outline' via derivatives of the hypergeometric function 2F1, with intermediate steps omitted. Since this is a headline exact analytic result, please provide a complete derivation, or alternatively state it as a conjecture supported by numerical evidence.
minor comments (5)
  1. [Sec. 3.1, near Eq. (3.14)] The text refers to 'defining C_n in Eq. (3.35)', but Eq. (3.35) is the two-point expansion; the definition appears in Eq. (3.10). Cross-reference error.
  2. [Eq. (3.23)] The notation P^{*ν_ij}_G(q_ij) ≡ (2π²/q³) P^{*ν_ij}_G(q_ij) is self-referential; the convolution power should be denoted with a different symbol or explicitly defined.
  3. [Eq. (2.12)] Extra closing parenthesis in 'F(ζ_G)) = C^{-1}(C_G(ζ_G))'.
  4. [Sec. 3.5 heading] Typo: '4P: T rispectra' should be 'Trispectra'.
  5. [Sec. 4, Eq. (4.6)] The text says the series (4.6) is asymptotic; it would help to also state clearly which of the subsequent series in the paper are exact and which are formal/asymptotic, especially in relation to (3.9) for N≥3.

Circularity Check

0 steps flagged

No significant circularity: the N-point expansion and the exponential-tail two-point result are derived from the model F and the Kibble–Slepian identity, not fitted to the target correlations.

full rationale

The paper's central object is the Kibble–Slepian decomposition (3.9), whose coefficients C_s are defined as one-dimensional Gaussian averages of the local mapping F (Eq. 3.10). These coefficients are computed from the model itself, not adjusted to reproduce any N-point function. The two-point function in the strongly non-Gaussian exponential-tail model, G2 = (ξ0/8β²) arcsin²(ξ_G/ξ0), is obtained by evaluating the double integral (3.28) for F(ζ) ∝ −ln|ζ|, a legitimate independent derivation. The subsequent extraction of C̄_n from this G2 (Eq. 4.10) is an inverse consistency step, not a fit masquerading as a prediction: the arcsin² result already stands on its own from the integral identity. The only self-citation, Ref. [33], is used to compare perturbativity bounds and is not load-bearing for the main derivation. The concern that the multivariate series (3.9) may not converge for N≥3 when F is non-analytic is a mathematical correctness issue about the interchange of expectation and series, not a circularity: it does not make the claimed result equivalent to its inputs by construction. Overall, no fitted parameter is renamed as a prediction, no load-bearing premise is justified solely by self-citation, and no derived quantity is defined in terms of the target outcome.

Axiom & Free-Parameter Ledger

3 free parameters · 5 axioms · 0 invented entities

The expansion uses only standard mathematics plus the assumed local mapping F; no new physical entities are introduced. The only ad hoc input is the illustrative model (4.1) and the template spectra.

free parameters (3)
  • β (or β̄ = β√ξ0)
    Model parameter controlling the strength of non-Gaussianity in the toy model (4.1); not fitted to data in this paper.
  • µ in BPL1 template = 0.15
    Shape parameter of the illustrative Gaussian power spectrum template in Eq. (4.12); chosen by hand, not fitted.
  • Normalization 2.86 in BPL2 = 2.86
    Peak normalization of the second template in Eq. (4.12); chosen for illustration, not fitted.
axioms (5)
  • standard math Kibble-Slepian formula (Eq. 3.6) for the multivariate Gaussian density with unit diagonal correlation matrix.
    Used as the foundation of the expansion in Sec. 3.1; cited to Refs. [21,22].
  • standard math The N-point function is defined by the multidimensional Gaussian average (2.14) over the auxiliary Gaussian field.
    Standard path-integral representation of Gaussian statistics; used throughout.
  • domain assumption The series (3.9) can be integrated term-by-term against non-analytic F(ζ_i).
    Not proven; needed to convert the exact density expansion into the N-point function expansion. The paper does not discuss convergence conditions.
  • domain assumption The field statistics are homogeneous and isotropic, so ξ_G depends only on |x1-x2|.
    Stated in Sec. 2.1, Eq. (2.7); restricts the framework to statistically homogeneous fields.
  • ad hoc to paper The absolute-value continuation of the log model (4.1) approximates the non-Gaussianity in USR/curvaton/phase-transition scenarios.
    Introduced following Ref. [20]; the paper explicitly notes lattice simulations of USR suggest a different continuation beyond 2βζ_G≳1 (Sec. 4).

pith-pipeline@v1.3.0-alltime-deepseek · 17962 in / 18110 out tokens · 145800 ms · 2026-08-03T03:03:52.798835+00:00 · methodology

0 comments
read the original abstract

We present a non-perturbative approach to correlation functions and polyspectra of locally non-Gaussian fields and develop a semi-perturbative framework that does not rely on a local expansion. This enables the computation of $N$-point functions of non-Gaussian fields even when the mapping between the auxiliary Gaussian field $\zeta_{\rm G}$ and the non-Gaussian field $\zeta({\bf x}) = F(\zeta_{\rm G}({\bf x}))$ is non-analytic. As an example, we consider non-Gaussian fields with exponentially tailed distributions, which can arise, for instance, in ultra-slow-roll models of inflation, and derive some exact analytic results in the strongly non-Gaussian regime.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Memoirs of the curvaton: non-perturbative non-Gaussianity and supermassive primordial black holes

    astro-ph.CO 2026-06 unverdicted novelty 6.0

    Curvaton self-interactions in non-quadratic potentials produce a local non-Gaussian map that enables supermassive primordial black hole formation at peak amplitudes of order 10^{-5} while remaining consistent with μ-d...

  2. Probing non-Gaussianity during reheating with SIGW in the LISA band

    astro-ph.CO 2026-04 unverdicted novelty 4.0

    Non-standard reheating imprints detectable features on SIGW spectra via non-Gaussianity, with dynamics that can suppress or boost the signal amplitude for LISA.

Reference graph

Works this paper leans on

30 extracted references · 1 canonical work pages · cited by 2 Pith papers

  1. [4]

    Tomita,Non-Linear Theory of Gravitational Instability in the Expanding Universe,Prog

    K. Tomita,Non-Linear Theory of Gravitational Instability in the Expanding Universe,Prog. Theor. Phys.37(1967) 831

    1967

  2. [5]

    Matarrese, O

    S. Matarrese, O. Pantano and D. Saez,A General relativistic approach to the nonlinear evolution of collisionless matter,Phys. Rev. D47(1993) 1311

    1993

  3. [6]

    Bruni, S

    M. Bruni, S. Matarrese, S. Mollerach and S. Sonego,Perturbations of space-time: Gauge transformations and gauge invariance at second order and beyond,Class. Quant. Grav.14 (1997) 2585 [gr-qc/9609040]

    Pith/arXiv arXiv 1997

  4. [7]

    Matarrese, S

    S. Matarrese, S. Mollerach and M. Bruni,Second order perturbations of the Einstein-de Sitter universe,Phys. Rev. D58(1998) 043504 [astro-ph/9707278]

    Pith/arXiv arXiv 1998

  5. [8]

    Ananda, C

    K.N. Ananda, C. Clarkson and D. Wands,The Cosmological gravitational wave background from primordial density perturbations,Phys. Rev. D75(2007) 123518 [gr-qc/0612013]

    Pith/arXiv arXiv 2007

  6. [9]

    Baumann, P.J

    D. Baumann, P.J. Steinhardt, K. Takahashi and K. Ichiki,Gravitational Wave Spectrum Induced by Primordial Scalar Perturbations,Phys. Rev. D76(2007) 084019 [hep-th/0703290]

    Pith/arXiv arXiv 2007

  7. [10]

    Dom` enech,Scalar Induced Gravitational Waves Review,Universe7(2021) 398 [2109.01398]

    G. Dom` enech,Scalar Induced Gravitational Waves Review,Universe7(2021) 398 [2109.01398]

    Pith/arXiv arXiv 2021

  8. [11]

    Zel’dovich and I.D

    Y.B. Zel’dovich and I.D. Novikov,The Hypothesis of Cores Retarded during Expansion and the Hot Cosmological Model,Sov. Astron.10(1967) 602

    1967

  9. [12]

    Hawking,Gravitationally collapsed objects of very low mass,Mon

    S. Hawking,Gravitationally collapsed objects of very low mass,Mon. Not. Roy. Astron. Soc. 152(1971) 75

    1971

  10. [13]

    Carr and S.W

    B.J. Carr and S.W. Hawking,Black holes in the early Universe,Mon. Not. Roy. Astron. Soc. 168(1974) 399

    1974

  11. [14]

    Carr,The Primordial black hole mass spectrum,Astrophys

    B.J. Carr,The Primordial black hole mass spectrum,Astrophys. J.201(1975) 1

    1975

  12. [15]

    B. Carr, F. Kuhnel and M. Sandstad,Primordial Black Holes as Dark Matter,Phys. Rev. D94 (2016) 083504 [1607.06077]

    Pith/arXiv arXiv 2016

  13. [16]

    Carr and F

    B. Carr and F. Kuhnel,Primordial Black Holes as Dark Matter: Recent Developments,Ann. Rev. Nucl. Part. Sci.70(2020) 355 [2006.02838]

    Pith/arXiv arXiv 2020

  14. [17]

    Carr, A.J

    B. Carr, A.J. Iovino, G. Perna, V. Vaskonen and H. Veerm¨ ae,Primordial black holes: constraints, potential evidence and prospects,2601.06024

    arXiv

  15. [18]

    Caravano, G

    A. Caravano, G. Franciolini and S. Renaux-Petel,Ultraslow-roll inflation on the lattice: Backreaction and nonlinear effects,Phys. Rev. D111(2025) 063518 [2410.23942]

    Pith/arXiv arXiv 2025

  16. [19]

    Caravano, G

    A. Caravano, G. Franciolini and S. Renaux-Petel,Ultraslow-roll inflation on the lattice. II. Nonperturbative curvature perturbation,Phys. Rev. D112(2025) 083508 [2506.11795]

    arXiv 2025

  17. [20]

    Iovino, S

    A.J. Iovino, S. Matarrese, G. Perna, A. Ricciardone and A. Riotto,How well do we know the scalar-induced gravitational waves?,Phys. Lett. B872(2026) 140039 [2412.06764]. – 21 –

    Pith/arXiv arXiv 2026

  18. [21]

    Kibble,An extension of a theorem of mehler’s on hermite polynomials,Mathematical Proceedings of the Cambridge Philosophical Society41(1945) 12–15

    W.F. Kibble,An extension of a theorem of mehler’s on hermite polynomials,Mathematical Proceedings of the Cambridge Philosophical Society41(1945) 12–15

    1945

  19. [22]

    D. Slepian,On the symmetrized kronecker power of a matrix and extensions of mehler’s formula for hermite polynomials,SIAM Journal on Mathematical Analysis3(1972) 606 [https://doi.org/10.1137/0503060]

    work page doi:10.1137/0503060 1972

  20. [23]

    Unal,Imprints of Primordial Non-Gaussianity on Gravitational Wave Spectrum,Phys

    C. Unal,Imprints of Primordial Non-Gaussianity on Gravitational Wave Spectrum,Phys. Rev. D99(2019) 041301 [1811.09151]

    Pith/arXiv arXiv 2019

  21. [24]

    Adshead, K.D

    P. Adshead, K.D. Lozanov and Z.J. Weiner,Non-Gaussianity and the induced gravitational wave background,JCAP10(2021) 080 [2105.01659]

    Pith/arXiv arXiv 2021

  22. [25]

    Perna, C

    G. Perna, C. Testini, A. Ricciardone and S. Matarrese,Fully non-Gaussian Scalar-Induced Gravitational Waves,JCAP05(2024) 086 [2403.06962]

    Pith/arXiv arXiv 2024

  23. [26]

    J.-P. Li, S. Wang, Z.-C. Zhao and K. Kohri,Isotropy, anisotropies and non-Gaussianity in the scalar-induced gravitational-wave background: diagrammatic approach for primordial non-Gaussianity up to arbitrary order,2505.16820

    Pith/arXiv arXiv

  24. [27]

    Vaskonen and H

    V. Vaskonen and H. Veerm¨ ae,Did NANOGrav see a signal from primordial black hole formation?,Phys. Rev. Lett.126(2021) 051303 [2009.07832]

    Pith/arXiv arXiv 2021

  25. [28]

    Franciolini, A

    G. Franciolini, A. Iovino, Junior., V. Vaskonen and H. Veermae,Recent Gravitational Wave Observation by Pulsar Timing Arrays and Primordial Black Holes: The Importance of Non-Gaussianities,Phys. Rev. Lett.131(2023) 201401 [2306.17149]

    Pith/arXiv arXiv 2023

  26. [29]

    Atal and C

    V. Atal and C. Germani,The role of non-gaussianities in Primordial Black Hole formation, Phys. Dark Univ.24(2019) 100275 [1811.07857]

    Pith/arXiv arXiv 2019

  27. [30]

    Tomberg,Stochastic constant-roll inflation and primordial black holes,Phys

    E. Tomberg,Stochastic constant-roll inflation and primordial black holes,Phys. Rev. D108 (2023) 043502 [2304.10903]

    Pith/arXiv arXiv 2023

  28. [31]

    Sasaki, J

    M. Sasaki, J. Valiviita and D. Wands,Non-Gaussianity of the primordial perturbation in the curvaton model,Phys. Rev. D74(2006) 103003 [astro-ph/0607627]

    Pith/arXiv arXiv 2006

  29. [32]

    Lewicki, P

    M. Lewicki, P. Toczek and V. Vaskonen,Black Holes and Gravitational Waves from Slow First-Order Phase Transitions,Phys. Rev. Lett.133(2024) 221003 [2402.04158]

    Pith/arXiv arXiv 2024

  30. [33]

    Iovino, Junior., G

    A. Iovino, Junior., G. Perna and H. Veerm¨ ae,The impact of non-Gaussianity when searching for Primordial Black Holes with LISA,2512.13648. – 22 –

    arXiv