pith. sign in

arxiv: 2605.28943 · v1 · pith:LF6AGDGNnew · submitted 2026-05-27 · 🌌 astro-ph.CO

Dominated-Convergence Failure in Cosmological Perturbation Theory and a Numerical Foundation for BBGKY+ZA

Svetlin V. Tassev This is my paper

Pith reviewed 2026-06-29 10:09 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords dominated convergencecosmological perturbation theoryEulerian perturbation theoryLagrangian perturbation theoryBBGKY hierarchyZel'dovich approximationlarge-scale structurephase-space cumulants
0
0 comments X

The pith

Perturbative solutions for dark matter overdensity diverge due to an invalid interchange of sums and Fourier integrals that violates the dominated convergence theorem.

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

The paper shows that the divergence of perturbative expansions for the overdensity in cosmological perturbation theory does not arise from the expansion in Lagrangian displacement itself. It arises instead from exchanging an infinite sum over orders with a Fourier integral, which violates the conditions of Lebesgue's dominated convergence theorem. This dominated convergence obstruction explains the limited convergence of Eulerian perturbation theory and restricts Lagrangian perturbation theory to a bounded validity region. The authors therefore propose evolving phase-space cumulants through the BBGKY hierarchy initialized by the Zel'dovich approximation, an approach that avoids the obstruction by never expanding the density in displacement.

Core claim

The perturbative solutions for δ diverge not because of the expansion in s per se, but because of an exchange of an infinite sum with a Fourier integral that violates the conditions of Lebesgue's dominated-convergence theorem. This DC obstruction is one clear reason why the convergence of Eulerian PT is controlled by advection terms beyond the linear δ. The same DC obstruction underlies LPT, whose region of validity is the resummation region of a DC-violating series bounded by shell crossing on one side and severely underdense regions on the other. Effective field theories need to smooth at short scales just to recover from that DC obstruction. An alternative is to evolve phase-space cumulan

What carries the argument

The dominated convergence obstruction that appears when an infinite sum over perturbative orders is interchanged with a Fourier integral during the expansion of the overdensity δ in the Lagrangian displacement s.

If this is right

  • The convergence of Eulerian perturbation theory is controlled by advection terms beyond the linear overdensity.
  • Lagrangian perturbation theory remains valid only inside the resummation region of a DC-violating series, bounded by shell crossing and severely underdense regions.
  • Effective field theories of large-scale structure must introduce short-scale smoothing to recover from the DC obstruction even when non-linearities beyond mass conservation are unimportant.
  • Evolving phase-space cumulants via the BBGKY hierarchy initialized with the Zel'dovich approximation removes the DC obstruction, so any remaining EFT corrections can target physics beyond mass conservation.
  • The Zel'dovich-approximation phase-space two-point function can be integrated numerically and used to construct the higher-order correlators required for closure in the BBGKY approach.

Where Pith is reading between the lines

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

  • If the DC obstruction is the dominant source of divergence, similar sum-integral interchanges may limit perturbative methods in other fluid or gravitational systems that expand densities in displacements.
  • Numerical implementation of the BBGKY+ZA scheme could be tested against N-body simulations to measure how far beyond shell crossing the method remains accurate.
  • The choice of Zel'dovich approximation for closure may need to be compared with other possible closures to determine which best captures the missing physics.
  • The separation of mass-conservation effects from other dynamical effects opens the possibility of hybrid schemes that apply standard PT on large scales and BBGKY+ZA on smaller scales.

Load-bearing premise

The divergence of the perturbative series is caused specifically by the sum-integral interchange violating dominated convergence rather than by other properties of the expansion or the underlying dynamics.

What would settle it

A direct numerical evaluation of the perturbative series for δ performed once with the sum and Fourier integral interchanged and once without the interchange, checking whether only the interchanged version diverges in regimes where the other remains finite.

read the original abstract

A common ingredient in cosmological perturbation theory (PT) is the expansion of the dark matter overdensity $\delta$ in the Lagrangian displacement $s$, which amounts to enforcing mass conservation perturbatively. In Eulerian PT (EPT), that expansion occurs already at the level of the continuity equation; in Lagrangian PT (LPT) it is done in the Poisson equation. We show that the resulting perturbative solutions for $\delta$ can diverge not because of the expansion in $s$ per se, but because of an exchange of an infinite sum with a Fourier integral that violates the conditions of Lebesgue's dominated-convergence (DC) theorem. We show that this DC obstruction (DCO) is one clear reason why the convergence of EPT is controlled by advection terms beyond the linear $\delta$. The same DCO underlies LPT: LPT's region of validity is the resummation region of a DC-violating series, bounded by shell crossing on one side and severely underdense regions on the other. Effective field theories (EFT) of large-scale structure need to smooth at short scales just to recover from that DCO, independent of whether non-linearities beyond mass conservation are important or not. An alternative is to never expand $\delta$ in $s$: instead evolve phase-space cumulants using the BBGKY hierarchy, initialized with the Zel'dovich approximation (ZA). The DCO is then absent by construction, so an EFT of BBGKY can focus on physics beyond mass conservation, which may allow pushing PT beyond shell crossing. The trade-off is the need for a closure relation, for which one can again use the ZA. We provide the building blocks for such a BBGKY+ZA recipe. A bottleneck for implementing it has been the ZA phase-space two-point function $\mathcal{P}$, which we successfully integrate numerically; we then write the higher ZA phase-space correlators needed for closure as products and convolutions of $\mathcal{P}$.

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

3 major / 2 minor

Summary. The paper claims that divergences in perturbative solutions for the dark matter overdensity δ in both Eulerian PT (EPT) and Lagrangian PT (LPT) arise not from the expansion in the Lagrangian displacement s itself, but from an interchange of an infinite sum over perturbative orders with a Fourier integral that violates the hypotheses of Lebesgue's dominated convergence theorem. It identifies this 'DC obstruction' (DCO) as controlling the convergence of EPT via advection terms and bounding the validity region of LPT between shell crossing and underdense regions. The manuscript argues that EFTs must smooth at short scales to recover from the DCO independently of other nonlinearities, and proposes an alternative BBGKY hierarchy initialized with the Zel'dovich approximation (ZA) that avoids the DCO by construction; it supplies the building blocks for such a scheme, including a numerical integration of the ZA phase-space two-point function P and expressions for higher-order ZA correlators as products and convolutions of P.

Significance. If the DCO diagnosis is substantiated with explicit kernels and verification of the DC conditions, the result would supply a precise mathematical reason for the limited radius of convergence in standard PT and a principled motivation for BBGKY-based approaches that can focus on physics beyond mass conservation. The successful numerical integration of P removes a stated bottleneck and constitutes a concrete, reproducible technical contribution.

major comments (3)
  1. [Abstract] Abstract and the paragraph beginning 'A common ingredient...': the central claim that the series diverges specifically because of a DC-violating sum-integral interchange (rather than the radius of the s-expansion or the fluid dynamics) requires an explicit verification that (i) the summed series kept inside the Fourier integral satisfies the hypotheses of Lebesgue DC (or converges pointwise), while (ii) the interchanged form violates them, and (iii) this is the dominant mechanism. The manuscript must exhibit the concrete kernel whose absolute integrability fails and the absent dominating function; without this step the diagnosis remains an assertion.
  2. [§3] §3 (LPT discussion): the statement that 'LPT's region of validity is the resummation region of a DC-violating series' needs a direct comparison showing that the resummed LPT expressions coincide with the DC-obstructed series inside its claimed domain and diverge outside it for the same mathematical reason; a concrete counter-example or theorem application at the level of the displacement field would make the claim load-bearing.
  3. [§5] §5 (BBGKY+ZA closure): the proposed ZA closure for the BBGKY hierarchy is presented as sufficient to push PT beyond shell crossing, but the manuscript does not demonstrate that the resulting equations remain consistent with the Vlasov-Poisson system or that truncation errors do not reintroduce an analogous obstruction; an explicit consistency check or numerical test against N-body data is required before the claim that 'an EFT of BBGKY can focus on physics beyond mass conservation' can be assessed.
minor comments (2)
  1. [§4] Notation for the phase-space two-point function P is introduced without an explicit integral representation or Fourier convention; adding the definition immediately after its first appearance would improve readability.
  2. [§4] The numerical integration procedure for P is described at a high level; a short appendix with the quadrature method, convergence tests, and a plot of the resulting P(k) would make the 'successful integration' claim verifiable.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address each major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract and the paragraph beginning 'A common ingredient...': the central claim that the series diverges specifically because of a DC-violating sum-integral interchange (rather than the radius of the s-expansion or the fluid dynamics) requires an explicit verification that (i) the summed series kept inside the Fourier integral satisfies the hypotheses of Lebesgue DC (or converges pointwise), while (ii) the interchanged form violates them, and (iii) this is the dominant mechanism. The manuscript must exhibit the concrete kernel whose absolute integrability fails and the absent dominating function; without this step the diagnosis remains an assertion.

    Authors: We agree that an explicit verification of the DC conditions strengthens the central claim. Section 2 identifies the advection term as the source of the obstruction and shows that the interchange violates absolute integrability while the un-interchanged series converges pointwise for the linear contribution. To make this fully load-bearing as requested, we will add a dedicated subsection exhibiting the concrete kernel (the Fourier transform of exp(ik·s) - 1) and verifying conditions (i)-(iii) directly. revision: yes

  2. Referee: [§3] §3 (LPT discussion): the statement that 'LPT's region of validity is the resummation region of a DC-violating series' needs a direct comparison showing that the resummed LPT expressions coincide with the DC-obstructed series inside its claimed domain and diverge outside it for the same mathematical reason; a concrete counter-example or theorem application at the level of the displacement field would make the claim load-bearing.

    Authors: We will revise §3 to include the requested direct comparison at the level of the displacement field s. This will apply Lebesgue's theorem explicitly, demonstrate coincidence of the resummed LPT form with the DC-obstructed series inside the validity region, and provide a concrete 1D counter-example showing divergence outside that region for the identical mathematical reason. revision: yes

  3. Referee: [§5] §5 (BBGKY+ZA closure): the proposed ZA closure for the BBGKY hierarchy is presented as sufficient to push PT beyond shell crossing, but the manuscript does not demonstrate that the resulting equations remain consistent with the Vlasov-Poisson system or that truncation errors do not reintroduce an analogous obstruction; an explicit consistency check or numerical test against N-body data is required before the claim that 'an EFT of BBGKY can focus on physics beyond mass conservation' can be assessed.

    Authors: The BBGKY hierarchy is derived directly from the Vlasov-Poisson equation, so the equations themselves remain consistent; the ZA enters only as a closure for higher correlators. The manuscript supplies the numerical building blocks (integration of P and expressions for higher correlators) and notes that the DCO is absent by construction, but does not perform a full numerical test or claim a completed implementation beyond shell crossing. We will revise §5 and the abstract to clarify the scope and state that full consistency validation against N-body data is future work. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation applies external DC theorem to standard PT expressions

full rationale

The paper diagnoses divergence in EPT/LPT as arising from a sum-integral interchange violating Lebesgue dominated convergence, then proposes BBGKY+ZA to avoid the expansion in s that triggers the DCO. This follows directly from applying an external theorem to existing perturbative kernels without any reduction of a claimed result to a fitted parameter, self-citation chain, or definitional equivalence. The BBGKY+ZA construction is motivated by the diagnosis rather than presupposing it, and the numerical integration of ZA phase-space correlators is presented as an independent implementation step. No load-bearing step collapses to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper invokes the standard mathematical statement of Lebesgue's dominated convergence theorem and the existing Zel'dovich approximation; no free parameters are introduced and no new physical entities are postulated.

axioms (1)
  • standard math Lebesgue's dominated convergence theorem governs the legitimacy of interchanging the infinite sum over perturbative orders with the Fourier integral in the expressions for δ
    Invoked in the abstract to diagnose the source of divergence in both EPT and LPT.

pith-pipeline@v0.9.1-grok · 5898 in / 1367 out tokens · 30066 ms · 2026-06-29T10:09:52.513020+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

35 extracted references · 15 canonical work pages · 14 internal anchors

  1. [1]

    Solving Large Scale Structure in Ten Easy Steps with COLA

    S. Tassev, M. Zaldarriaga and D. Eisenstein,Solving Large Scale Structure in Ten Easy Steps with COLA,JCAP06(2013) 036 [1301.0322]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  2. [2]

    sCOLA: The N-body COLA Method Extended to the Spatial Domain

    S. Tassev, D.J. Eisenstein, B.D. Wandelt and M. Zaldarriaga,sCOLA: The N-body COLA Method Extended to the Spatial Domain,1502.07751

    work page internal anchor Pith review Pith/arXiv arXiv

  3. [3]

    Valageas,Using the zeldovich dynamics to test expansion schemes,Astronomy & Astrophysics476(2007) 31–58

    P. Valageas,Using the zeldovich dynamics to test expansion schemes,Astronomy & Astrophysics476(2007) 31–58

    2007

  4. [4]

    Large-Scale Structure of the Universe and Cosmological Perturbation Theory

    F. Bernardeau, S. Colombi, E. Gaztanaga and R. Scoccimarro,Large scale structure of the universe and cosmological perturbation theory,Phys. Rept.367(2002) 1 [astro-ph/0112551]

    work page internal anchor Pith review Pith/arXiv arXiv 2002

  5. [5]

    Lagrangian Dynamics in Non-Flat Universes and Non-Linear Gravitational Evolution

    P. Catelan,Lagrangian dynamics in nonflat universes and nonlinear gravitational evolution, Mon. Not. Roy. Astron. Soc.276(1995) 115 [astro-ph/9406016]

    work page internal anchor Pith review Pith/arXiv arXiv 1995

  6. [6]

    Resumming Cosmological Perturbations via the Lagrangian Picture: One-loop Results in Real Space and in Redshift Space

    T. Matsubara,Resumming Cosmological Perturbations via the Lagrangian Picture: One-loop Results in Real Space and in Redshift Space,Phys. Rev. D77(2008) 063530 [0711.2521]

    work page internal anchor Pith review Pith/arXiv arXiv 2008

  7. [7]

    The Mildly Non-Linear Regime of Structure Formation

    S. Tassev and M. Zaldarriaga,The Mildly Non-Linear Regime of Structure Formation,JCAP 04(2012) 013 [1109.4939]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  8. [8]

    Lagrangian or Eulerian; Real or Fourier? Not All Approaches to Large-Scale Structure Are Created Equal

    S. Tassev,Lagrangian or Eulerian; Real or Fourier? Not All Approaches to Large-Scale Structure Are Created Equal,JCAP06(2014) 008 [1311.4884]

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  9. [9]

    The Lagrangian-space Effective Field Theory of Large Scale Structures

    R.A. Porto, L. Senatore and M. Zaldarriaga,The Lagrangian-space Effective Field Theory of Large Scale Structures,JCAP05(2014) 022 [1311.2168]

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  10. [10]

    Porto, L

    R.A. Porto, L. Senatore and M. Zaldarriaga,The lagrangian-space effective field theory of large scale structures,Journal of Cosmology and Astroparticle Physics2014(2014) 022–022

    2014

  11. [11]

    Cosmological Non-Linearities as an Effective Fluid

    D. Baumann, A. Nicolis, L. Senatore and M. Zaldarriaga,Cosmological Non-Linearities as an Effective Fluid,JCAP07(2012) 051 [1004.2488]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  12. [12]

    The Effective Field Theory of Cosmological Large Scale Structures

    J.J.M. Carrasco, M.P. Hertzberg and L. Senatore,The Effective Field Theory of Cosmological Large Scale Structures,JHEP09(2012) 082 [1206.2926]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  13. [13]

    Will,The confrontation between general relativity and experiment,Living Reviews in Relativity17(2014)

    C.M. Will,The confrontation between general relativity and experiment,Living Reviews in Relativity17(2014)

    2014

  14. [14]

    Carlson, B

    J. Carlson, B. Reid and M. White,Convolution lagrangian perturbation theory for biased tracers,Monthly Notices of the Royal Astronomical Society429(2012) 1674–1685

    2012

  15. [15]

    Peebles,The large-scale structure of the universe, Princeton University Press (1980)

    P.J.E. Peebles,The large-scale structure of the universe, Princeton University Press (1980)

    1980

  16. [16]

    Tassev,The helmholtz hierarchy: phase space statistics of cold dark matter,Journal of Cosmology and Astroparticle Physics2011(2011) 022–022

    S.V. Tassev,The helmholtz hierarchy: phase space statistics of cold dark matter,Journal of Cosmology and Astroparticle Physics2011(2011) 022–022

    2011

  17. [17]

    Perturbative interaction approach to cosmological structure formation

    Y. Ali-Ha¨ ımoud,Perturbative interaction approach to cosmological structure formation, Phys. Rev. D91(2015) 103507 [1502.00580]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  18. [18]

    Nascimento and M

    C. Nascimento and M. Loverde,Cosmological perturbation theory for large scale structure in phase space,JCAP06(2025) 002 [2410.05389]

    work page arXiv 2025

  19. [19]

    Connection between Newtonian simulations and general relativity

    N.E. Chisari and M. Zaldarriaga,Connection between Newtonian simulations and general relativity,Phys. Rev. D83(2011) 123505 [1101.3555]

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  20. [20]

    Separate Universe Simulations

    C. Wagner, F. Schmidt, C.-T. Chiang and E. Komatsu,Separate Universe Simulations,Mon. Not. Roy. Astron. Soc.448(2015) L11 [1409.6294]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  21. [21]

    Lorenz,The problem of deducing the climate from the governing equations,Tellus16 (1964) 1

    E.N. Lorenz,The problem of deducing the climate from the governing equations,Tellus16 (1964) 1

    1964

  22. [22]

    Tao,An Introduction to Measure Theory, Graduate Studies in Mathematics, American Mathematical Society (2011)

    T. Tao,An Introduction to Measure Theory, Graduate Studies in Mathematics, American Mathematical Society (2011). – 72 –

    2011

  23. [23]

    NIST Digital Library of Mathematical Functions

    “NIST Digital Library of Mathematical Functions.”https://dlmf.nist.gov/, Release 1.2.1 of 2024-06-15

    2024

  24. [24]

    Zel’dovich,Gravitational instability: An approximate theory for large density perturbations.,Astronomy & Astrophysics5(1970) 84

    Y.B. Zel’dovich,Gravitational instability: An approximate theory for large density perturbations.,Astronomy & Astrophysics5(1970) 84

    1970

  25. [25]

    Valageas,A new approach to gravitational clustering: A path-integral formalism and large-nexpansions,Astronomy & Astrophysics421(2004) 23–40

    P. Valageas,A new approach to gravitational clustering: A path-integral formalism and large-nexpansions,Astronomy & Astrophysics421(2004) 23–40

    2004

  26. [26]

    Aghanim, Y

    Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi et al., Planck 2018 results. VI. Cosmological parameters,Astronomy & Astrophysics641(2020) A6

    2018

  27. [27]

    Hamilton,Uncorrelated modes of the non-linear power spectrum,Monthly Notices of the Royal Astronomical Society312(2000) 257–284

    A.J.S. Hamilton,Uncorrelated modes of the non-linear power spectrum,Monthly Notices of the Royal Astronomical Society312(2000) 257–284

    2000

  28. [28]

    Cosmological Perturbation Theory Using the FFTLog: Formalism and Connection to QFT Loop Integrals

    M. Simonovi´ c, T. Baldauf, M. Zaldarriaga, J.J. Carrasco and J.A. Kollmeier,Cosmological perturbation theory using the FFTLog: formalism and connection to QFT loop integrals,JCAP 04(2018) 030 [1708.08130]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  29. [29]

    Schneider and M

    P. Schneider and M. Bartelmann,The power spectrum of density fluctuations in the zel’dovich approximation,Monthly Notices of the Royal Astronomical Society273(1995) 475

    1995

  30. [30]

    Konrad and M

    S. Konrad and M. Bartelmann,On the asymptotic behaviour of cosmic density-fluctuation power spectra,Monthly Notices of the Royal Astronomical Society515(2022) 2578

    2022

  31. [31]

    Konrad, Y.B

    S. Konrad, Y.B. Ginat and M. Bartelmann,On the asymptotic behaviour of cosmic density-fluctuation power spectra of cold dark matter,Monthly Notices of the Royal Astronomical Society515(2022) 5823

    2022

  32. [32]

    Bardeen, J.R

    J.M. Bardeen, J.R. Bond, N. Kaiser and A.S. Szalay,The Statistics of Peaks of Gaussian Random Fields,Astrophys. J.304(1986) 15

    1986

  33. [33]

    The FLINT team,FLINT: Fast Library for Number Theory, 2023

    2023

  34. [34]

    Abramowitz and I.A

    M. Abramowitz and I.A. Stegun,Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover, New York (1964). – 73 – 10−4 10−3 10−2 10−1 100 k [h/Mpc] 10−1 100 101 102 103 104 105 106 107 108 109 P(k, √ p ˆk, −√ p ˆk) e−σ2 vp [(Mpc/h)3] PZ (k), D= 0.40 PZ (k), D= 1 PZ (k), D= 1.41 PZ (k), D= 100.00 Pe−σ2 v p|√p=0.2 h/Mpc (a) 10−...

    1964

  35. [35]

    Y×(1−µ lm)

    ratio scaled the same way and evaluated at Z= 1/4, 1/2, 3/4, 1. As per eq. (5.35), all curves should asymptote to 1 in the lin- ear regime, which the figure confirms (for the top-right panel, see below). Columns scan µlm ∈ {−1,−0.7,0,+0.7,+1}. The top-right cell (µ lm = +1) uses adifferentnormaliza- tion than the rest of the top row: dropping the (1−µ lm)...