Linear perturbative theory of the discrete cosmological N-body problem

Abstract

We present a perturbative treatment of the evolution under their mutual self-gravity of particles displaced off an infinite perfect lattice, both for a static space and for a homogeneously expanding space as in cosmological N-body simulations. The treatment, analogous to that of perturbations to a crystal in solid state physics, can be seen as a discrete (i.e. particle) generalization of the perturbative solution in the Lagrangian formalism of a self-gravitating fluid. Working to linear order, we show explicitly that this fluid evolution is recovered in the limit that the initial perturbations are restricted to modes of wavelength much larger than the lattice spacing. The full spectrum of eigenvalues of the simple cubic lattice contains both oscillatory modes and unstable modes which grow slightly faster than in the fluid limit. A detailed comparison of our perturbative treatment, at linear order, with full numerical simulations is presented, for two very different classes of initial perturbation spectra. We find that the range of validity is similar to that of the perturbative fluid approximation (i.e. up to close to ``shell-crossing''), but that the accuracy in tracing the evolution is superior. The formalism provides a powerful tool to systematically calculate discreteness effects at early times in cosmological N-body simulations.

25 pages, 21 figures

Cited by in corpus (29)

• Large-scale dark matter simulations
• The $\text{Abacus}$ Cosmological $N$-body Code
• Accurate estimators of power spectra in N-body simulations
• Evaluating backreaction with the peak model of structure formation
• Higher-order initial conditions for mixed baryon-CDM simulations
• Towards quantitative control on discreteness error in the non-linear regime of cosmological N body simulations
• Aemulus $\nu$: Precise Predictions for Matter and Biased Tracer Power Spectra in the Presence of Neutrinos
• Quantification of discreteness effects in cosmological N-body simulations: II. Evolution up to shell crossing
• Higher-order initial conditions with massive neutrinos
• The DESI $N$-body Simulation Project I: Testing the Robustness of Simulations for the DESI Dark Time Survey
• Gravitational Dynamics of an Infinite Shuffled Lattice of Particles
• Perturbation theory with dispersion and higher cumulants: non-linear regime
• Gravitational Dynamics of an Infinite Shuffled Lattice: Particle Coarse-grainings, Non-linear Clustering and the Continuum Limit
• Force distribution in a randomly perturbed lattice of identical particles with $1/r^2$ pair interaction
• Gravitational dynamics of an infinite shuffled lattice: early time evolution and universality of non-linear correlations
• Particle linear theory on a self-gravitating perturbed cubic Bravais lattice
• Accuracy of power spectra in dissipationless cosmological simulations
• BullFrog: Multi-step perturbation theory as a time integrator for cosmological simulations
• Statistical physics for cosmic structures
• Numerical convergence of pre-initial conditions on dark matter halo properties
• Perturbation-theory informed integrators for cosmological simulations
• Yukawa vs. Newton: gravitational forces in a cubic cosmological simulation box
• Starting Cosmological Simulations from the Big Bang
• Vlasov limit and discreteness effects in cosmological N-body simulations
• The Hidden Role of Anisotropies in Shaping Structure Formation in Cosmological N-Body Simulations
• On the origin of transient features in cosmological N-Body Simulations
• Effect of the cubic torus topology on cosmological perturbations
• Tidal adaptive softening and artificial fragmentation in cosmological simulations
• Particle loads for cosmological simulations with equal-mass dark matter and baryonic particles