Bridging perturbation theory and simulations: initial conditions and fast integrators for cosmological simulations

Oliver Hahn

REVIEW 3 major objections 5 minor 24 references

Bridging perturbation theory and simulations: initial conditions and fast integrators for cosmological simulations

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

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

Pith's one-line read Cosmological N-body simulations are most accurate when started late with high-order Lagrangian perturbation theory and stepped in growth-factor time.

desk verdict A clean, genuinely useful set of lecture notes that consolidates established LPT/IC material; the only real weakness is that the headline late-start z=24/3LPT recommendation is borrowed from one convergence study and may not generalize as broadly as the takeaways suggest. read the letter →

arxiv 2608.03400 v1 pith:JBVIQTT5 submitted 2026-08-04 astro-ph.CO

Oliver Hahn This is my paper · ORCID

classification astro-ph.CO

keywords cosmologicalN-bodysimulationsinitialconditionsLagrangianperturbationtheoryZel'dovichapproximationdiscretenesserrorsGaussianrandomfieldsVlasov-PoissonequationsPT-informedtimeintegrators

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

These lecture notes connect two usually separate pieces of cosmological simulation practice: how initial conditions are drawn and how particles are advanced in time. They make the case that the standard early-start, first-order practice is the wrong end of the error budget: starting a simulation as late as possible from high-order Lagrangian perturbation theory, for example 3LPT at z=24, keeps both perturbation-theory truncation and lattice discreteness errors small. For time stepping, they argue that integrating in the linear growth factor $D_+$ rather than cosmic time, with velocities rescaled accordingly, yields integrators that exactly reproduce the Zel'dovich solution for one-dimensional data in a single step and can match second-order LPT trajectories, converging on large scales with few time steps. A sympathetic reader comes away with a concrete recipe: late starts, high LPT order, and growth-factor-based integrators give more accurate non-linear structure at fixed computational cost.

What carries the argument

The central machinery is the Lagrangian map $X(q,t) = q + \Psi(q,t)$ of the Vlasov-Poisson system, expanded perturbatively in LPT, together with the change of time variable from cosmic time $t$ to the linear growth factor $D_+$. The key identity is the Zel'dovich-consistency condition $\beta = 1 - \alpha$ for the $D_+$-time drift-kick-drift integrator, which says the kick coefficient must be one minus the drift coefficient for the integrator to reproduce inertial Zel'dovich motion exactly; the BULLFROG variant fixes $\alpha$ from the second-order growth factor $E = D^{(2)}$ so that the trajectory matches 2LPT. The same $D_+$ variable is the time coordinate in which Zel'dovich motion is straight and the growth factor used to assemble nLPT displacements and velocities.

What would settle it

Run identical simulations initialized at early times with 1LPT and at late times with 3LPT, using several softenings and both cubic-lattice and glass particle loadings, and compare their z=0 power spectra against a converged high-resolution reference; the paper's claim predicts the late-start run is closer on small scales, and a reversal or disappearance of that gap would falsify the discreteness freeze-in assumption.

Watch

Extended reading notes

Core claim

The paper's central claim is that the accuracy of a cosmological N-body simulation is governed by two competing errors, LPT truncation and particle discreteness, and that both are minimized by initializing from high-order LPT as late as possible, because discreteness errors from the lattice grow during the linear phase and freeze into the final power spectrum. On the time-integration side, the claim is that re-writing the Vlasov-Poisson characteristics in terms of the growth factor $D_+$ and the rescaled velocity $W = V/(a^2 \dot D_+)$ turns the drift-kick-drift step into a family of Zel'dovich-consistent integrators, with coefficient condition $\beta = 1 - \alpha$; one member, the BULLFROG integrator, matches trajectories to 2LPT and is asserted to produce the most accurate non-linear evolution on large scales with few time steps.

Load-bearing premise

The load-bearing premise is that lattice discreteness errors accumulate during the linear phase and freeze into the final power spectrum, with the growth rate computed for a simple cubic lattice and a specific force softening; if that freeze-in is weaker, or different for other particle loadings or softening, the optimal starting redshift and LPT order could shift.

Editorial extensions

If this is right

Late-start high-order LPT initial conditions, such as 3LPT at z=24, should reduce frozen-in discreteness errors in the z=0 and z=1 matter power spectra relative to early-start 1LPT runs at equal cost.
PT-informed $D_+$ integrators such as BULLFROG require substantially fewer time steps for converged large-scale clustering, making high-accuracy simulations cheaper.
Standard cosmic-time leapfrog integrators are not Zel'dovich consistent, meaning they do not exactly reproduce even inertial Zel'dovich motion for one-dimensional initial data in one step.
LPT is convergent only up to the shell-crossing singularity, so the perturbative initialization and the discrete N-body evolution are complementary: the notes' recipe hands off to N-body just before shell-crossing contaminates the expansion.

Reading between the lines

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

A natural extension the notes leave implicit is that the same late-start logic should apply to higher-order statistics such as the bispectrum, since the discreteness freeze-in mechanism is scale- and statistic-generic; this could be tested by comparing bispectra from early- and late-start runs.
A 3LPT-matched $D_+$ integrator is a natural next step; the notes' 2LPT matching suggests it would converge even faster in time steps while further sacrificing canonical phase-space structure.
The discreteness error analysis is carried out for a simple cubic lattice, so both the optimal starting redshift and the optimal LPT order could shift for glass or other particle loadings and for different force softenings; the z=24/3LPT recipe should therefore be re-validated outside the cubic-lattice setting.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

Short version: this is a textbook-style lecture note, not a research paper, and it should be judged as teaching material. It makes no claim to new results—the author is explicit about that—and almost everything substantive traces to prior work. What it does well is assemble those pieces into a single coherent path from Gaussian random fields to Fourier sampling, the Vlasov-Poisson system, LPT, and PT-informed integrators, with clean derivations and a companion notebook that actually implements the steps. The proof of Fourier diagonality is careful, the characteristic derivation of ZA is clear, and the Zel'dovich-consistency theorem for the D-time integrator family is a nice compact piece of exposition. For someone entering the field or teaching it, this is probably a more efficient entry point than the primary literature.

The soft spots are real but minor, and they concern presentation more than physics. The quantitative anchor of the whole practical message—"3LPT at z=24"—comes from [14], a convergence study for a particular lattice and force softening. The notes cite that study rather than reproduce or bound the analysis, and the abstract and takeaway 7 turn it into a general rule. A student could easily take z=24/3LPT as universal, when in fact the optimal starting redshift and LPT order likely depend on particle loading, force softening, and time stepping. That caveat deserves to be in the main text, not only in the reference. Similarly, the claim that BULLFROG produces "the most accurate non-linear evolution on large scales with few time steps" is a comparative statement from specific tests; it is likely true for the tested cases but is stated without the same hedging the authors would give it in a research paper.

The derivations I checked are correct; no internal contradictions or circular fitting. The citation pattern is self-heavy, but those are the relevant papers and they are published and checkable, so this is not a problem. For the venue—SciPost Physics Lecture Notes—this is exactly the kind of material the venue exists for.

Recommendation: send it to peer review, but ask the authors to soften the generality of the late-start/high-order recommendation, attribute it explicitly to the regime tested in [14], and add a sentence on how the optimum could shift for other setups. I would bring it to reading group when we have new students starting simulations.

Referee Report

3 major / 5 minor

Summary. These lecture notes (arXiv:2608.03400) aim to bridge analytic perturbation theory and N-body simulation practice. Section 1 defines Gaussian random fields, proves the diagonality of their Fourier-space covariance, and describes efficient FFT-based sampling on a periodic box. Section 2 introduces the cosmological Vlasov-Poisson system, derives the characteristic equations, develops the Boltzmann hierarchy, obtains the Zel'dovich approximation as early-time asymptotics, and presents Lagrangian perturbation theory up to all orders via recurrence relations. Section 3 discusses N-body initial conditions, standard leapfrog integrators, PT-informed integrators (FAST-PM, BULLFROG) with a Zel'dovich-consistency theorem, and discreteness errors, concluding that simulations are best initialized with high-order LPT at late times (e.g., 3LPT at z=24) and that PT-informed integrators converge on large scales with few time steps. The notes are accompanied by a Jupyter notebook.

Significance. The notes are clearly written and didactically effective: the proof of Fourier diagonality, the LPT determinant expansion, and the Zel'dovich-consistency theorem are presented cleanly, and the accompanying notebook is a concrete strength. If the practical recommendations are accepted, the notes would usefully update standard practice by moving students away from early 1LPT starts toward late high-order LPT starts and PT-informed stepping. However, the two main practical claims are imported from the author's published papers ([14], [20], [22]) rather than derived or tested here, so the notes' value as a self-contained bridge depends on those references; the caveats attached to those claims in the original papers should be preserved.

major comments (3)

[Section 3.4 and Conclusion (takeaway 7)] The statement that 'simulations are most accurately initialized with high order LPT at very late times, e.g. 3LPT at z=24 in their analysis' is the paper's central practical takeaway, but it is a quotation of the conclusion of [14] for a simple cubic lattice, a specific force softening, and a particular time-integrator setup. The notes do not state this scope, so a reader could apply '3LPT at z=24' to glass or grid-shuffled particle loadings, different softening schemes, or other integrators where the optimum may shift. Please add an explicit caveat after Figure 3 and rephrase takeaway 7 to say, for example, 'for the configurations tested in [14]' or 'when using a simple cubic lattice and the force softening adopted in [14]'.
[Section 2.2, Eq. (15)] The formal solution to the Poisson equation is printed as φ(X,t) = κ/a ∇_X^{-2}(J^{-1}/J). Since n = 1/J and the Poisson equation (10b) has source 1-n, the correct expression is φ(X,t) = κ/a ∇_X^{-2}(1/J - 1), equivalently κ/a ∇_X^{-2}((1-J)/J). The later master equation (28) uses the correct 1-1/J form, so this is a typo rather than a propagated error, but in a lecture note it will mislead readers and should be corrected.
[Section 3.3, after Eq. (44)] The sentence 'This integrator produces the most accurate non-linear evolution on large scales with few time steps' is an unqualified superlative. The notes should specify the comparison class (e.g., the integrators discussed in this section), the error metric (e.g., power spectrum or bispectrum residuals), and the step-count regime, or point to the specific figure in [22] that supports the claim. Without this, takeaway 6 overstates the evidence and invites misuse of BULLFROG outside the tested configurations.

minor comments (5)

[Section 2.4, Eq. (24)] From Eq. (23) and the prefactor κ/(a^3 f_g^2 H^2 D) ≍ 3/(2a), one obtains X' + ∇_X φ = -(2a/3)X''; the printed sign is positive. The limit X' → -∇_q φ_0 is unaffected, but the displayed arithmetic should be corrected.
[Section 1.2, proof of Theorem 2] After the substitution r = x-y and X = x, the exponent should be -i[(k-k')·X + k'·r], not -i[(k-k')·X + k·r]; the subsequent delta function makes the final equality correct, but the intermediate expression is misleading.
[Section 3.1, step 3] Equations (36)-(37) use growth factors D^{(1)}, D^{(2)} without specifying their normalization; please state the convention (e.g., D_+(a=1)=1) so that the relation to the CAMB/CLASS output P_m(k)/D_+(z_target)^2 is unambiguous.
[Figure 3, right panel] The right panel contains many overlapping curves distinguished only by linestyle and color; consider enlarging the panel or splitting it by z, since the comparison of start redshifts and LPT orders is a key visual result.
[Sections 2.5 and 3.1] The lectures referenced in the text ('Romain Teyssier's lecture', 'lectures by Cora Uhlemann') do not appear in the reference list; please add full citations or URLs.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the derivations are self-contained and the numerical recommendations rest on cited external experiments, not on refitted inputs.

full rationale

The notes derive the Gaussian random field sampling, the Vlasov-Poisson system, the Zel'dovich approximation, and 1LPT/2LPT/nLPT recurrence relations from stated equations rather than importing the target conclusions. The late-start/high-order-LPT recommendation (Section 3.4) and the Bullfrog accuracy claim (Section 3.3) are explicitly presented as conclusions of the author's prior numerical studies [14] and [22], which are published, externally checkable benchmarks, and the notes do not fit parameters or derive the target from an input. The Zel'dovich-consistency theorem (Section 3.3) is a definitional algebraic property of the D-time drift-kick-drift family, and the condition alpha + beta = 1 is not an input-output tautology. No equation in the notes is equivalent to a target prediction by construction, and no fitted parameter is renamed as a prediction. Self-citations appear, but because the cited works contain independent numerical evidence rather than definitions that assume the notes' conclusions, they do not constitute circularity.

Assumptions & free parameters 0 free parameters · 7 assumptions · 0 invented entities

The notes introduce no free parameters and no new entities. The central recommendations rest on standard math (Fourier analysis, Isserlis theorem, determinant identities) and on domain assumptions about the cosmological fluid and the accuracy of cited numerical error analyses. The most fragile of these is the discreteness freeze-in assumption, imported from [14] and [23].

assumptions (7)

standard math Isserlis/Wick theorem: centered joint Gaussian moments factor into sums over pairings of covariances.
Used in Section 1.1 to characterize GRF moments and justify the vanishing of odd correlators; proved in textbooks.
standard math Standard Fourier transform identities and the diagonalization of homogeneous covariance on a torus.
Used in Section 1.2 to derive the power-spectrum sampling formula Eq. (9).
standard math Determinant expansion for det(I + epsilon A) and von Neumann series for (I + epsilon A)^-1.
Used in Section 2.5 to obtain the LPT expansion of 1 - 1/J.
domain assumption The primordial density field is a homogeneous, isotropic Gaussian random field with power spectrum from Einstein-Boltzmann codes.
Invoked in Section 1 as the starting point for all cosmological simulations, motivated by CMB observations.
domain assumption Cold collisionless matter is governed by the Vlasov-Poisson system with zero initial velocity dispersion before shell crossing.
Used in Section 2 to derive the characteristic equations and the Euler-Poisson hierarchy.
domain assumption LPT converges up to the shell-crossing singularity, so sufficiently high orders describe the displacement field just before collapse.
Used in Section 2.5 and 3.4 to justify late-time starts; convergence is demonstrated numerically for spherical collapse and attributed to [13].
domain assumption The discreteness error of a lattice-based N-body simulation follows the fluid-lattice growth rate of Joyce et al. (2005) and freezes into the power spectrum during non-linear evolution.
Underlies the Section 3.4 recommendation to start as late as possible; taken from [23] and [14] rather than derived in the notes.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Bridging perturbation theory and simulations: initial conditions and fast integrators for cosmological simulations." pith.science (2026). https://pith.science/paper/JBVIQTT5

@misc{pith2026260803400,
  author       = {Pith},
  title        = {Pith review of: Bridging perturbation theory and simulations: initial conditions and fast integrators for cosmological simulations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JBVIQTT5}},
  note         = {Machine review of arXiv:2608.03400}
}

read the original abstract

These lecture notes provide an introduction to the generation of initial conditions for cosmological N-body simulations. Starting from the definition and properties of Gaussian random fields, we discuss their role in cosmology and the efficient generation of such fields using Fourier methods. The Vlasov-Poisson system is introduced as the governing framework for cold collisionless matter, and its solution via characteristics and Lagrangian perturbation theory (LPT) is detailed. We discuss the use of LPT for initializing N-body simulations, emphasizing the importance of high-order LPT and late-time starts to minimize truncation and discreteness errors. Finally, we discuss time integration schemes, including PT-informed integrators, and their role in accurately evolving the system. These notes aim to bridge the gap between theoretical perturbation methods and practical simulation techniques.

Figures

Figures reproduced from arXiv: 2608.03400 by the authors.

**Figure 2.** Left: Norm of the displacement field ∥ψ(n) ∥ as a function of the initial Lagrangian coordinate q up to 4LPT. A darker color corresponds to a larger displacement, each panel is individually normalised. The slowest converging regions are the almost spherically collapsing and expanding regions. (Adapted from [13]) Right: Convergence of the radius ζ(a) of a spherical perturbation in LPT, where ζ(a) = 1 + PN n=1 a nψ(n… view at source ↗

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 4 canonical work pages

[14]

Michaux, O

M. Michaux, O. Hahn, C. Rampf and R. E. Angulo,Accurate initial conditions for cosmo- logical N-body simulations: minimizing truncation and discreteness errors, MNRAS500(1), 663 (2021), doi:10.1093/mnras/staa3149, 2008.09588

arXiv 2021
[20]

List and O

F . List and O. Hahn,Perturbation-theory informed integrators for cosmological simulations, Journal of Computational Physics513, 113201 (2024), doi:10.1016/j.jcp.2024.113201, 2301.09655

arXiv 2024
[22]

Rampf, F

C. Rampf, F . List and O. Hahn,BULLFROG: multi-step perturbation theory as a time integrator for cosmological simulations, J. Cosmology Astropart. Phys.2025(2), 020 (2025), doi:10.1088/1475-7516/2025/02/020, 2409.19049

arXiv 2025
[1]

R. J. Adler,The Geometry of Random Fields, Society for Industrial and Applied Mathemat- ics, doi:10.1137/1.9780898718980 (2010), https://epubs.siam.org/doi/pdf/10.1137/1. 9780898718980. 14 SciPost Physics Lecture Notes Submission

doi:10.1137/1.9780898718980 2010
[2]

P . J. E. Peebles,The Large-Scale Structure of the Universe, Princeton University Press, ISBN 9780691209838 (1980)

work page 1980
[3]

R. E. Angulo and O. Hahn,Large-scale dark matter simulations, Living Reviews in Computational Astrophysics8(1), 1 (2022), doi:10.1007/s41115-021-00013-z, 2112. 05165

doi:10.1007/s41115-021-00013-z 2022
[4]

C. Rampf,Cosmological Vlasov-Poisson equations for dark matter: Recent developments and connections to selected plasma problems, arXiv e-prints arXiv:2110.06265 (2021), doi:10.48550/arXiv.2110.06265, 2110.06265

arXiv doi:10.48550/arxiv.2110.06265 2021
[5]

A. D. Chernin, D. I. Nagirner and S. V . Starikova,Growth rate of cosmological perturbations in standard model: Explicit analytical solution, A&A399, 19 (2003), doi:10.1051/0004- 6361:20021763, astro-ph/0110107

arXiv 2003

Show all 24 references

[6]

Brenier, U

Y . Brenier, U. Frisch, M. Hénon, G. Loeper, S. Matarrese, R. Mohayaee and A. Sobolevski˘i, Reconstruction of the early Universe as a convex optimization problem, MNRAS346(2), 501 (2003), doi:10.1046/j.1365-2966.2003.07106.x, astro-ph/0304214

2003
[7]

Y. B. Zel’dovich,Gravitational instability: An approximate theory for large density pertur- bations., A&A5, 84 (1970)

1970
[8]

Buchert and J

T . Buchert and J. Ehlers,Lagrangian theory of gravitational instability of Friedman- Lemaitre cosmologies – second-order approach: an improved model for non-linear clustering, MNRAS264, 375 (1993), doi:10.1093/mnras/264.2.375

1993 doi
[9]

F . R. Bouchet, S. Colombi, E. Hivon and R. Juszkiewicz,Perturbative Lagrangian approach to gravitational instability., A&A296, 575 (1995), doi:10.48550/arXiv.astro-ph/9406013, astro-ph/9406013

1995 doi arXiv
[10]

Rampf,The recursion relation in Lagrangian perturbation theory, J

C. Rampf,The recursion relation in Lagrangian perturbation theory, J. Cosmology Astropart. Phys.2012(12), 004 (2012), doi:10.1088/1475-7516/2012/12/004, 1205.5274

2012 arXiv
[11]

Zheligovsky and U

V . Zheligovsky and U. Frisch,Time-analyticity of Lagrangian particle trajectories in ideal fluid flow, Journal of Fluid Mechanics749, 404 (2014), doi:10.1017/jfm.2014.221, 1312.6320

2014 arXiv
[12]

Matsubara,Recursive solutions of Lagrangian perturbation theory, Phys

T . Matsubara,Recursive solutions of Lagrangian perturbation theory, Phys. Rev. D92(2), 023534 (2015), doi:10.1103/PhysRevD.92.023534, 1505.01481

2015 arXiv
[13]

Rampf and O

C. Rampf and O. Hahn,Shell-crossing in a ΛCDM Universe, MNRAS501(1), L71 (2021), doi:10.1093/mnrasl/slaa198, 2010.12584

2021 arXiv
[15]

Efstathiou, M

G. Efstathiou, M. Davis, S. D. M. White and C. S. Frenk,Numerical techniques for large cosmological N-body simulations, ApJS57, 241 (1985), doi:10.1086/191003

1985 doi
[16]

Crocce, S

M. Crocce, S. Pueblas and R. Scoccimarro,Transients from initial conditions in cosmologi- cal simulations, MNRAS373(1), 369 (2006), doi:10.1111/j.1365-2966.2006.11040.x, astro-ph/0606505

2006
[17]

Adamek, D

J. Adamek, D. Daverio, R. Durrer and M. Kunz,General relativity and cosmic structure formation, Nature Physics12(4), 346 (2016), doi:10.1038/nphys3673, 1509.01699. 15 SciPost Physics Lecture Notes Submission

2016 arXiv
[18]

Quinn, N

T . Quinn, N. Katz, J. Stadel and G. Lake,Time stepping N-body simulations, arXiv e-prints astro-ph/9710043 (1997), doi:10.48550/arXiv.astro-ph/9710043, astro-ph/9710043

1997 doi arXiv
[19]

Hairer, C

E. Hairer, C. Lubich and G. Wanner,Geometric numerical integration, vol. 31 ofSpringer Series in Computational Mathematics, Springer-Verlag, Berlin, second edn., ISBN 3-540- 30663-3; 978-3-540-30663-4, Structure-preserving algorithms for ordinary differential equations (2006)

2006
[21]

Feng, M.-Y

Y. Feng, M.-Y. Chu, U. Seljak and P . McDonald,FASTPM: a new scheme for fast simulations of dark matter and haloes, MNRAS463(3), 2273 (2016), doi:10.1093/mnras/stw2123, 1603.00476

2016 arXiv
[23]

Joyce, B

M. Joyce, B. Marcos, A. Gabrielli, T . Baertschiger and F . Sylos Labini,Gravitational Evolution of a Perturbed Lattice and its Fluid Limit, Phys. Rev. Lett.95(1), 011304 (2005), doi:10.1103/PhysRevLett.95.011304, astro-ph/0504213

2005 arXiv
[24]

L. H. Garrison, D. J. Eisenstein, D. Ferrer, M. V . Metchnik and P . A. Pinto,Improving initial conditions for cosmological N-body simulations, MNRAS461(4), 4125 (2016), doi:10.1093/mnras/stw1594, 1605.02333. 16

2016 arXiv

Pith tools

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

[1] [14]

Michaux, O

M. Michaux, O. Hahn, C. Rampf and R. E. Angulo,Accurate initial conditions for cosmo- logical N-body simulations: minimizing truncation and discreteness errors, MNRAS500(1), 663 (2021), doi:10.1093/mnras/staa3149, 2008.09588

arXiv 2021

[2] [20]

List and O

F . List and O. Hahn,Perturbation-theory informed integrators for cosmological simulations, Journal of Computational Physics513, 113201 (2024), doi:10.1016/j.jcp.2024.113201, 2301.09655

arXiv 2024

[3] [22]

Rampf, F

C. Rampf, F . List and O. Hahn,BULLFROG: multi-step perturbation theory as a time integrator for cosmological simulations, J. Cosmology Astropart. Phys.2025(2), 020 (2025), doi:10.1088/1475-7516/2025/02/020, 2409.19049

arXiv 2025

[4] [1]

R. J. Adler,The Geometry of Random Fields, Society for Industrial and Applied Mathemat- ics, doi:10.1137/1.9780898718980 (2010), https://epubs.siam.org/doi/pdf/10.1137/1. 9780898718980. 14 SciPost Physics Lecture Notes Submission

doi:10.1137/1.9780898718980 2010

[5] [2]

P . J. E. Peebles,The Large-Scale Structure of the Universe, Princeton University Press, ISBN 9780691209838 (1980)

work page 1980

[6] [3]

R. E. Angulo and O. Hahn,Large-scale dark matter simulations, Living Reviews in Computational Astrophysics8(1), 1 (2022), doi:10.1007/s41115-021-00013-z, 2112. 05165

doi:10.1007/s41115-021-00013-z 2022

[7] [4]

C. Rampf,Cosmological Vlasov-Poisson equations for dark matter: Recent developments and connections to selected plasma problems, arXiv e-prints arXiv:2110.06265 (2021), doi:10.48550/arXiv.2110.06265, 2110.06265

arXiv doi:10.48550/arxiv.2110.06265 2021

[8] [5]

A. D. Chernin, D. I. Nagirner and S. V . Starikova,Growth rate of cosmological perturbations in standard model: Explicit analytical solution, A&A399, 19 (2003), doi:10.1051/0004- 6361:20021763, astro-ph/0110107

arXiv 2003