Perturbation-theory informed integrators for cosmological simulations

Abstract

Large-scale cosmological simulations are an indispensable tool for modern cosmology. To enable model-space exploration, fast and accurate predictions are critical. In this paper, we show that the performance of such simulations can be further improved with time-stepping schemes that use input from cosmological perturbation theory. Specifically, we introduce a class of time-stepping schemes derived by matching the particle trajectories in a single leapfrog/Verlet drift-kick-drift step to those predicted by Lagrangian perturbation theory (LPT). As a corollary, these schemes exactly yield the analytic Zel'dovich solution in 1D in the pre-shell-crossing regime (i.e. before particle trajectories cross). One representative of this class is the popular FastPM scheme by Feng et al. 2016, which we take as our baseline. We then construct more powerful LPT-inspired integrators and show that they outperform FastPM and standard integrators in fast simulations in two and three dimensions with $\mathcal{O}(1-100)$ timesteps, requiring less steps to accurately reproduce the power spectrum and bispectrum of the density field. Furthermore, we demonstrate analytically and numerically that, for any integrator, convergence is limited in the post-shell-crossing regime (to order 3/2 for planar wave collapse), owing to the lacking regularity of the acceleration field, which makes the use of high-order integrators in this regime futile. Also, we study the impact of the timestep spacing and of a decaying mode present in the initial conditions. Importantly, we find that symplecticity of the integrator plays a minor role for fast approximate simulations with a small number of timesteps.

40 + 11 pages, 12 + 3 figures. Accepted by JCP. Added 3D convergence study, added kick-drift-kick version for some integrators

References in corpus (19)

• Array Programming with NumPy
• The Magnus expansion and some of its applications
• Transients from Initial Conditions in Cosmological Simulations
• The CAMELS project: Cosmology and Astrophysics with MachinE Learning Simulations
• A critical look at cosmological perturbation theory techniques
• AbacusSummit: A Massive Set of High-Accuracy, High-Resolution $N$-Body Simulations
• Large-scale dark matter simulations
• The $\text{Abacus}$ Cosmological $N$-body Code
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• Linear perturbative theory of the discrete cosmological N-body problem
• Towards quantitative control on discreteness error in the non-linear regime of cosmological N body simulations
• FlowPM: Distributed TensorFlow Implementation of the FastPM Cosmological N-body Solver
• Phase-space structure of protohalos: Vlasov versus Particle-Mesh
• Cosmological perturbations for two cold fluids in $\Lambda$CDM
• Shell-crossing in a $\Lambda$CDM Universe
• Cosmological Vlasov-Poisson equations for dark matter: Recent developments and connections to selected plasma problems
• Gravitational instability via the Schrodinger equation
• Analytical growth functions for cosmic structures in a $\Lambda$CDM Universe
• Herglotz Action for Homogeneous Cosmology

Cited by in corpus (3)

• BullFrog: Multi-step perturbation theory as a time integrator for cosmological simulations
• COmoving Computer Acceleration (COCA): $N$-body simulations in an emulated frame of reference
• Testing the Dark Universe through the Layzer-Irvine Equation