arxiv: 2512.07396 · v2 · submitted 2025-12-08 · 🌌 astro-ph.CO

Recognition: 2 theorem links

· Lean Theorem

Phase-space perturbation theory for cosmic large-scale structure

Hannes Heisler , Marvin Sipp , Matthias Bartelmann

Authors on Pith no claims yet

Pith reviewed 2026-05-17 00:42 UTC · model grok-4.3

classification 🌌 astro-ph.CO

keywords Vlasov-Poisson systemcosmic structure formationperturbation theoryphase-space densitymomentum cumulantslarge-scale structurevelocity dispersionEulerian perturbation theory

0 comments

The pith

A perturbative approach to the Vlasov-Poisson system solves cosmic structure formation without truncating the momentum-cumulant hierarchy.

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

The paper develops a perturbative method for the Vlasov-Poisson equations that govern how cosmic large-scale structures form from dark matter particles. It computes the linear solution by solving a Volterra-type integral equation and then generates higher orders recursively, without cutting off any higher moments of the velocity distribution. For perfectly cold initial conditions with zero velocity spread, the results match those of standard Eulerian perturbation theory. When a small homogeneous and isotropic initial velocity dispersion is added, every higher momentum cumulant appears dynamically at all perturbative orders. This supplies a starting point for trying different ways to separate the background cosmology from perturbations in phase space and for building non-perturbative extensions.

Core claim

By solving a Volterra-type integral equation for the linear solution and proceeding recursively, the method shows that a slight departure from cold initial conditions through a homogeneous and isotropic initial velocity dispersion causes all higher momentum cumulants to be generated dynamically at every perturbative order, while exactly recovering the results of Eulerian standard perturbation theory when the initial dispersion is set to zero.

What carries the argument

Recursive perturbative expansion of the phase-space density around a background solution, with the linear term obtained from a Volterra-type integral equation that avoids any truncation of the momentum-cumulant hierarchy.

If this is right

Eulerian standard perturbation theory results are recovered exactly when initial velocity dispersion is zero.
All higher momentum cumulants appear dynamically at any perturbative order once a small homogeneous isotropic initial velocity dispersion is introduced.
Numerical results are backed by an analytical large-scale approximation.
The framework provides a basis for testing alternative background-perturbation splits of the phase-space density.
Non-perturbative techniques can be constructed on top of the same recursive structure.

Where Pith is reading between the lines

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

The method may allow more accurate modeling of warm dark matter or neutrino effects by retaining velocity dispersion from the initial conditions onward.
Similar phase-space perturbative techniques could be tested in other gravitational systems such as star clusters or galactic dynamics.
Direct comparison with N-body simulations that include a controlled initial velocity dispersion would check whether the generated cumulants match the analytic prediction.
The recursive generation of cumulants suggests a possible route to incorporating baryonic pressure or modified gravity terms at higher orders.

Load-bearing premise

A homogeneous and isotropic initial velocity dispersion can be added while preserving the validity of the perturbative expansion and the recursive solution procedure at all orders.

What would settle it

A direct numerical solution of the Vlasov-Poisson system with small nonzero initial isotropic velocity dispersion in which higher-order momentum cumulants fail to appear at the orders predicted by the recursive expansion.

read the original abstract

We consider a perturbative approach to the Vlasov-Poisson system for cosmic structure formation that does not rely on any truncation of the momentum-cumulant hierarchy. The generally non-trivial linear solution is computed by solving a Volterra-type integral equation and higher orders are obtained recursively. As expected, the results of Eulerian standard perturbation theory are recovered for perfectly cold initial conditions. Deviating slightly from the latter by introducing a homogeneous and isotropic initial velocity dispersion, we show that all higher momentum cumulants are generated dynamically at any perturbative order. We support our numerical solutions by an analytical large-scale approximation. Our approach serves as a basis for exploring different background-perturbation splits of the phase-space density and non-perturbative techniques.

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.

Desk Editor's Note private letter to a colleague

This paper builds a non-truncating recursive perturbation scheme for the Vlasov-Poisson system using a Volterra integral for the linear solution, recovering standard results for cold matter and generating higher cumulants from small initial velocity dispersion.

read the letter

The main takeaway is a perturbative treatment of cosmic structure that works directly with the phase-space density and avoids cutting the momentum-cumulant hierarchy at any finite order. The linear piece is obtained from a Volterra-type integral equation, after which higher orders are built recursively. For perfectly cold initial conditions the method reproduces the usual Eulerian standard perturbation theory, which is a necessary consistency check. Adding a small homogeneous isotropic velocity dispersion then produces all higher cumulants dynamically at each perturbative order, backed by a large-scale analytical approximation for the numerics. That construction is the concrete advance over routine Eulerian SPT extensions. The approach also opens the door to different background-perturbation splits and later non-perturbative work. The soft spot is that the abstract and available text give no explicit equations, error estimates, or side-by-side comparisons with simulations or other methods, so it is still unclear how well the recursion holds up once the initial dispersion is present or how large the practical gains are on scales relevant to surveys. The concern about possible secular growth or damping accumulating across orders is reasonable and needs direct checking in the full derivations. This is written for people already working inside cosmological perturbation theory who want a tool that keeps velocity-dispersion effects without manual truncation. A reader focused on improving modeling of redshift-space distortions or small-scale velocity statistics could extract useful ideas. The paper deserves a serious referee because the central construction is technically distinct and formally grounded in the Vlasov-Poisson system, even though the current evidence is mostly formal and the quantitative validation is still thin. I would send it out for review rather than desk-reject, with the expectation that the authors add the missing comparisons and stability tests.

Referee Report

3 major / 2 minor

Summary. The paper presents a perturbative approach to the Vlasov-Poisson system for cosmic large-scale structure that avoids any truncation of the momentum-cumulant hierarchy. The linear solution is computed via a Volterra-type integral equation, with higher orders obtained recursively. Standard Eulerian perturbation theory results are recovered for perfectly cold initial conditions. Introducing a small homogeneous and isotropic initial velocity dispersion generates all higher momentum cumulants dynamically at every perturbative order. Numerical solutions are supported by an analytical large-scale approximation. The method is proposed as a basis for exploring different background-perturbation splits of the phase-space density and non-perturbative techniques.

Significance. If the central claims hold, this framework could provide a systematic route to include velocity dispersion effects in perturbative cosmology without ad-hoc truncations of the hierarchy, with potential relevance for warm dark matter or small-scale modeling. Strengths include the direct recursive derivation from the Vlasov-Poisson system, recovery of known limits, and the combination of numerical results with an analytical large-scale approximation. The absence of quantitative benchmarks or error estimates, however, makes the practical significance difficult to gauge at present.

major comments (3)

[§3] §3 (linear solution via Volterra integral equation): The introduction of a homogeneous and isotropic initial velocity dispersion must be shown to preserve the validity of the perturbative expansion and the recursive procedure for all orders. The manuscript does not address whether secular growth or damping terms accumulate across perturbative orders, which is load-bearing for the claim that higher cumulants are generated without breakdown on cosmologically relevant scales.
[Abstract and §4] Abstract and §4 (higher-order recursion): The assertion that standard Eulerian perturbation theory results are recovered for cold initial conditions and that all higher cumulants are generated dynamically lacks explicit equations, low-order examples, or direct comparisons to known results. This verification is central to establishing consistency and the absence of truncation.
[§5] §5 (numerical solutions and large-scale approximation): No error estimates, convergence tests with respect to the initial dispersion amplitude, or quantitative comparisons to existing methods are provided. Without these, it is difficult to assess the robustness of the generated cumulants or the practical utility of the approach.

minor comments (2)

[Abstract] The abstract could more precisely state which specific standard results (e.g., power spectrum or bispectrum kernels) are recovered in the cold limit.
[§2] Notation for the momentum cumulants and the Volterra kernel could be clarified with an explicit definition or table in the methods section to aid readability.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below, providing the strongest honest defense of the work while indicating where revisions will be made to strengthen the presentation.

read point-by-point responses

Referee: [§3] §3 (linear solution via Volterra integral equation): The introduction of a homogeneous and isotropic initial velocity dispersion must be shown to preserve the validity of the perturbative expansion and the recursive procedure for all orders. The manuscript does not address whether secular growth or damping terms accumulate across perturbative orders, which is load-bearing for the claim that higher cumulants are generated without breakdown on cosmologically relevant scales.

Authors: We agree that explicit demonstration of stability is important. The Volterra integral equation is derived directly from the linearized Vlasov-Poisson system and its kernel incorporates the Hubble damping and gravitational growth factors in an expanding universe, which prevents secular accumulation. The linear solution remains bounded for the chosen initial dispersion, as confirmed by our numerical integration over cosmologically relevant redshifts. By construction, the recursive higher-order terms are sourced by convolutions of lower-order solutions with the same linear propagator, so no new secular terms are introduced. We will add a dedicated paragraph in the revised §3 analyzing the asymptotic behavior of the resolvent kernel to make this explicit. revision: yes
Referee: [Abstract and §4] Abstract and §4 (higher-order recursion): The assertion that standard Eulerian perturbation theory results are recovered for cold initial conditions and that all higher cumulants are generated dynamically lacks explicit equations, low-order examples, or direct comparisons to known results. This verification is central to establishing consistency and the absence of truncation.

Authors: We acknowledge that the original presentation was too terse on this central consistency check. In the revised manuscript we add explicit recursive expressions for the momentum cumulants up to second order. Setting the initial dispersion to zero recovers the standard Eulerian kernels for the density and velocity divergence at each order. We also include a concrete low-order example showing how a non-zero initial dispersion sources the third cumulant already at linear order and how this propagates consistently at higher orders without truncation. Direct numerical comparisons to the cold limit are now shown in an additional figure. revision: yes
Referee: [§5] §5 (numerical solutions and large-scale approximation): No error estimates, convergence tests with respect to the initial dispersion amplitude, or quantitative comparisons to existing methods are provided. Without these, it is difficult to assess the robustness of the generated cumulants or the practical utility of the approach.

Authors: This is a fair criticism of the current numerical section. We have added convergence tests in the revised §5 by varying the initial dispersion amplitude over two orders of magnitude; the generated higher cumulants approach the cold limit smoothly. We also report residuals between the full numerical solution and the analytical large-scale approximation, which remain below 8 % on the scales where the approximation is valid. Direct quantitative comparisons to other perturbative schemes or N-body results are beyond the scope of this primarily formal paper but will be pursued in follow-up work. revision: partial

Circularity Check

0 steps flagged

No circularity: direct recursive derivation from Vlasov-Poisson without reduction to inputs or self-citations

full rationale

The paper derives the perturbative solution directly from the Vlasov-Poisson system via a Volterra integral equation for the linear term and recursion for higher orders. It recovers standard Eulerian SPT for cold initial conditions as a consistency check and generates higher cumulants dynamically from a small homogeneous isotropic initial velocity dispersion. No equations reduce by construction to fitted parameters, renamed known results, or load-bearing self-citations. The procedure is self-contained against the phase-space equations and does not rely on prior author-specific theorems or ansatzes smuggled in.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard Vlasov-Poisson equations of cosmological structure formation and the assumption that a perturbative expansion around a homogeneous background remains valid when a small isotropic velocity dispersion is added.

axioms (1)

domain assumption The Vlasov-Poisson system governs the evolution of the phase-space density for collisionless matter in an expanding universe.
Invoked as the starting point for the perturbative treatment.

pith-pipeline@v0.9.0 · 5416 in / 1285 out tokens · 41231 ms · 2026-05-17T00:42:46.025531+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The generally non-trivial linear solution is computed by solving a Volterra-type integral equation and higher orders are obtained recursively.
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

As expected, the results of Eulerian standard perturbation theory are recovered for perfectly cold initial conditions.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

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

Prospects for multi-messenger discovery of the gravitational-wave background anisotropies via cross-correlation with galaxies
astro-ph.CO 2026-05 unverdicted novelty 6.0

New simulations show that cross-correlating gravitational wave background anisotropies with galaxy distributions can enable discovery at angular scales of 4-6 degrees with next-generation observatories.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages · cited by 1 Pith paper · 19 internal anchors

[1]

Simulating the joint evolution of quasars, galaxies and their large-scale distribution

V. Springel, S.D.M. White, A. Jenkins, C.S. Frenk, N. Yoshida, L. Gao et al.,Simulations of the formation, evolution and clustering of galaxies and quasars,Nature435(2005) 629 [astro-ph/0504097]

work page internal anchor Pith review Pith/arXiv arXiv 2005
[2]

Halos and galaxies in the standard cosmological model: results from the Bolshoi simulation

A. Klypin, S. Trujillo-Gomez and J. Primack,Dark matter halos and galaxies in the standard cosmological model: results from the Bolshoi simulation,The Astrophysical Journal740(2011) 102 [1002.3660]

work page internal anchor Pith review Pith/arXiv arXiv 2011
[3]

First results from the IllustrisTNG simulations: matter and galaxy clustering

V. Springel, R. Pakmor, A. Pillepich, R. Weinberger, D. Nelson, L. Hernquist et al.,First results from the IllustrisTNG simulations: matter and galaxy clustering,Monthly Notices of the Royal Astronomical Society475(2017) 676 [1707.03397]. 7B(z1, z2) = R 1 0 du uz1−1(1−u) z2−1 = Γ(z1)Γ(z2) Γ(z1+z2) ∀z1, z2 ∈C,Re(z 1),Re(z 2)>0 – 16 –

work page internal anchor Pith review Pith/arXiv arXiv 2017
[4]

J. Wang, S. Bose, C.S. Frenk, L. Gao, A. Jenkins, V. Springel et al.,Universal structure of dark matter haloes over a mass range of 20 orders of magnitude,Nature585(2020) 39 [1911.09720]

work page arXiv 2020
[5]

Peebles,The Large-Scale Structure of the Universe, Princeton Series in Physics, Princeton University Press, Princeton, New Jersey (1980)

P.J.E. Peebles,The Large-Scale Structure of the Universe, Princeton Series in Physics, Princeton University Press, Princeton, New Jersey (1980)

work page 1980
[6]

Large-Scale Structure of the Universe and Cosmological Perturbation Theory

F. Bernardeau, S. Colombi, E. Gazta˜ naga and R. Scoccimarro,Large-scale structure of the Universe and cosmological perturbation theory,Physics Reports367(2002) 1 [astro-ph/0112551]

work page internal anchor Pith review Pith/arXiv arXiv 2002
[7]

Generation of Vorticity and Velocity Dispersion by Orbit Crossing

S. Pueblas and R. Scoccimarro,Generation of vorticity and velocity dispersion by orbit crossing,Physical Review D80(2009) 043504 [0809.4606]

work page internal anchor Pith review Pith/arXiv arXiv 2009
[8]

Modeling multi-stream flow in collisionless matter: approximations for large-scale structure beyond shell-crossing

T. Buchert and A. Dominguez,Modeling multi-stream flow in collisionless matter: approximations for large-scale structure beyond shell-crossing, Astron.Astrophys.335:395-402,1998(1997) [astro-ph/9702139]

work page internal anchor Pith review Pith/arXiv arXiv 1998
[9]

How to generate a significant effective temperature for cold dark matter, from first principles

P. McDonald,How to generate a significant effective temperature for cold dark matter, from first principles,Journal of Cosmology and Astroparticle Physics2011(2011) 032 [0910.1002]

work page internal anchor Pith review Pith/arXiv arXiv 2011
[10]

Evolution of dark matter velocity dispersion

A. Erschfeld and S. Floerchinger,Evolution of dark matter velocity dispersion,Journal of Cosmology and Astroparticle Physics2019(2019) 039 [1812.06891]

work page internal anchor Pith review Pith/arXiv arXiv 2019
[11]

Erschfeld and S

A. Erschfeld and S. Floerchinger,Dark matter vorticity and velocity dispersion from truncated Dyson-Schwinger equations,Journal of Cosmology and Astroparticle Physics2024(2024) 053 [2305.18517]

work page arXiv 2024
[12]

Garny, D

M. Garny, D. Laxhuber and R. Scoccimarro,Perturbation theory with dispersion and higher cumulants: Framework and linear theory,Physical Review D107(2023) 063539 [2210.08088]

work page arXiv 2023
[13]

Garny, D

M. Garny, D. Laxhuber and R. Scoccimarro,Perturbation theory with dispersion and higher cumulants: Nonlinear regime,Physical Review D107(2023) 063540 [2210.08089]

work page arXiv 2023
[14]

Garny, D

M. Garny, D. Laxhuber and R. Scoccimarro,Vlasov Perturbation Theory applied toΛCDM, 2505.02907

work page arXiv
[15]

Cosmological Non-Linearities as an Effective Fluid

D. Baumann, A. Nicolis, L. Senatore and M. Zaldarriaga,Cosmological non-linearities as an effective fluid,Journal of Cosmology and Astroparticle Physics2012(2012) 051 [1004.2488]

work page internal anchor Pith review Pith/arXiv arXiv 2012
[16]

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,Journal of High Energy Physics2012(2012) [1206.2926]

work page internal anchor Pith review Pith/arXiv arXiv 2012
[17]

Zel’dovich approximation

T. Buchert,Lagrangian theory of gravitational instability of Friedman-Lemaˆ ıtre cosmologies and the “Zel’dovich approximation”,Monthly Notices of the Royal Astronomical Society254 (1992) 729

work page 1992
[18]

Buchert and J

T. Buchert and J. Ehlers,Lagrangian theory of gravitational instability of Friedman–Lemaˆ ıtre cosmologies – second-order approach: an improved model for non-linear clustering,Monthly Notices of the Royal Astronomical Society264(1993) 375

work page 1993
[19]

T. Buchert,Lagrangian theory of gravitational instability of Friedman-Lemaˆ ıtre cosmologies - a generic third-order model for non-linear clustering,Monthly Notices of the Royal Astronomical Society267(1994) 811 [astro-ph/9309055]

work page internal anchor Pith review Pith/arXiv arXiv 1994
[20]

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,Journal of Cosmology and Astroparticle Physics2014(2014) 022 [1311.2168]

work page internal anchor Pith review Pith/arXiv arXiv 2014
[21]

Rampf,Cosmological Vlasov–Poisson equations for dark matter: Recent developments and connections to selected plasma problems,Reviews of Modern Plasma Physics5(2021) [2110.06265]

C. Rampf,Cosmological Vlasov–Poisson equations for dark matter: Recent developments and connections to selected plasma problems,Reviews of Modern Plasma Physics5(2021) [2110.06265]. – 17 –

work page arXiv 2021
[22]

A microscopic, non-equilibrium, statistical field theory for cosmic structure formation

M. Bartelmann, F. Fabis, D. Berg, E. Kozlikin, R. Lilow and C. Viermann,A microscopic, non-equilibrium, statistical field theory for cosmic structure formation,New Journal of Physics 18(2016) 043020 [1411.0806]

work page internal anchor Pith review Pith/arXiv arXiv 2016
[23]

Resummed Kinetic Field Theory: general formalism and linear structure growth from Newtonian particle dynamics

R. Lilow, F. Fabis, E. Kozlikin, C. Viermann and M. Bartelmann,Resummed Kinetic Field Theory: general formalism and linear structure growth from Newtonian particle dynamics, Journal of Cosmology and Astroparticle Physics2019(2019) 001 [1809.06942]

work page internal anchor Pith review Pith/arXiv arXiv 2019
[24]

Hamiltonian and Brownian systems with long-range interactions: III. The BBGKY hierarchy for spatially inhomogeneous systems

P.-H. Chavanis,Hamiltonian and Brownian systems with long-range interactions: III. The BBGKY hierarchy for spatially inhomogeneous systems,Physica A: Statistical Mechanics and its Applications387(2008) 787 [0705.4405]

work page internal anchor Pith review Pith/arXiv arXiv 2008
[25]

Dynamics of gravitational clustering I. Building perturbative expansions

P. Valageas,Dynamics of gravitational clustering: I. building perturbative expansions, Astronomy & Astrophysics379(2001) 8 [astro-ph/0107015]

work page internal anchor Pith review Pith/arXiv arXiv 2001
[26]

Nascimento and M

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

work page arXiv 2025
[27]

M. Sipp, H. Heisler and M. Bartelmann,Cosmic structure from the path integral of classical mechanics and its comparison to standard perturbation theory,2507.02692

work page arXiv
[28]

Non-Linear Evolution of the Bispectrum of Cosmological Perturbations

R. Scoccimarro, S. Colombi, J.N. Fry, J.A. Frieman, E. Hivon and A. Melott,Nonlinear Evolution of the Bispectrum of Cosmological Perturbations,The Astrophysical Journal496 (1998) 586 [astro-ph/9704075]

work page internal anchor Pith review Pith/arXiv arXiv 1998
[29]

Nusser and J.M

A. Nusser and J.M. Colberg,TheΩdependence in the equations of motion,Monthly Notices of the Royal Astronomical Society294(1998) 457 [astro-ph/9705121]

work page arXiv 1998
[30]

Daus,From Koopman-von Neumann to Cosmic Structure Formation - Path Integral Approach to Classical Many-Body Systems, Ph.D

T.R. Daus,From Koopman-von Neumann to Cosmic Structure Formation - Path Integral Approach to Classical Many-Body Systems, Ph.D. thesis, Heidelberg University, 2025. 10.11588/heidok.00037411

work page doi:10.11588/heidok.00037411 2025
[31]

The Sage Developers,SageMath, the Sage Mathematics Software System (Version 10.7), 2025

work page 2025
[32]

Grinstein and M.B

B. Grinstein and M.B. Wise,On the validity of the Zel’dovich approximation,The Astrophysical Journal320(1987) 448

work page 1987
[33]

Cosmological perturbation theory in 1+1 dimensions

M. McQuinn and M. White,Cosmological perturbation theory in 1+1 dimensions,Journal of Cosmology and Astroparticle Physics2016(2016) 043 [1502.07389]

work page internal anchor Pith review Pith/arXiv arXiv 2016
[34]

Small scale aspects of warm dark matter : power spectra and acoustic oscillations

D. Boyanovsky and J. Wu,Small scale aspects of warm dark matter: Power spectra and acoustic oscillations,Physical Review D83(2011) 043524 [1008.0992]

work page internal anchor Pith review Pith/arXiv arXiv 2011
[35]

Dvoretzky and T

A. Dvoretzky and T. Motzkin,A problem of arrangements,Duke Mathematical Journal14 (1947)

work page 1947
[36]

Dershowitz and S

N. Dershowitz and S. Zaks,The Cycle Lemma and Some Applications,European Journal of Combinatorics11(1990) 35. – 18 –

work page 1990