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.

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

References in corpus (64)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• The cosmological simulation code GADGET-2
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES
• The Magnus expansion and some of its applications
• The frontier of simulation-based inference
• Transients from Initial Conditions in Cosmological Simulations
• Solving Large Scale Structure in Ten Easy Steps with COLA
• Transients from Initial Conditions: A Perturbative Analysis
• The Quijote simulations
• The CAMELS project: Cosmology and Astrophysics with MachinE Learning Simulations
• The Uchuu Simulations: Data Release 1 and Dark Matter Halo Concentrations
• HACC: Simulating Sky Surveys on State-of-the-Art Supercomputing Architectures
• A critical look at cosmological perturbation theory techniques
• AbacusSummit: A Massive Set of High-Accuracy, High-Resolution $N$-Body Simulations
• FastPM: a new scheme for fast simulations of dark matter and halos
• Contact Hamiltonian Mechanics
• Tracing the Dark Matter Sheet in Phase Space
• Large-scale dark matter simulations
• Reconstruction of the early Universe as a convex optimization problem
• Physical Bayesian modelling of the non-linear matter distribution: new insights into the Nearby Universe
• The $\text{Abacus}$ Cosmological $N$-body Code
• Newtonian Cosmology in Lagrangian Formulation: Foundations and Perturbation Theory
• TreePM: A code for Cosmological N-Body Simulations
• ICE-COLA: Towards fast and accurate synthetic galaxy catalogues optimizing a quasi $N$-body method
• A new approach to simulating collisionless dark matter fluids
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• Tree-Particle-Mesh: an adaptive, efficient, and parallel code for collisionless cosmological simulation
• Recursive Solutions of Lagrangian Perturbation Theory
• A Convenient Set of Comoving Cosmological Variables and Their Application
• Fast and accurate mock catalogue generation for low-mass galaxies
• Cosmological perturbation theory in 1+1 dimensions
• Dark matter statistics for large galaxy catalogs: power spectra and covariance matrices
• Time-analyticity of Lagrangian particle trajectories in ideal fluid flow
• Improving Initial Conditions for Cosmological $N$-Body Simulations
• Fourth Order Gradient Symplectic Integrator Methods for Solving the Time-Dependent Schrödinger Equation
• The growth of structure in the intergalactic medium
• The recursion relation in Lagrangian perturbation theory
• Demonstrating Discreteness and Collision Error in Cosmological N-body Simulations of Dark Matter Gravitational Clustering
• Cosmological neutrino simulations at extreme scale
• 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
• Zeldovich flow on cosmic vacuum background: new exact nonlinear analytical solution
• Contact variational integrators
• A fast particle-mesh simulation of non-linear cosmological structure formation with massive neutrinos
• Post-collapse perturbation theory in 1D cosmology -- beyond shell-crossing
• Phase-space structure of protohalos: Vlasov versus Particle-Mesh
• Structure formation beyond shell-crossing: nonperturbative expansions and late-time attractors
• Lagrangian cosmological perturbation theory at shell-crossing
• Cosmological perturbations for two cold fluids in $\Lambda$CDM
• Semiclassical path to cosmic large-scale structure
• Shell-crossing in a $\Lambda$CDM Universe
• HACC Cosmological Simulations: First Data Release
• Unveiling the singular dynamics in the cosmic large-scale structure
• Gravitational instability via the Schrodinger equation
• Scalar Fields and the FLRW Singularity
• Cosmological Vlasov-Poisson equations for dark matter: Recent developments and connections to selected plasma problems
• Numerical integration in celestial mechanics: a case for contact geometry
• Herglotz Action for Homogeneous Cosmology
• Analytical growth functions for cosmic structures in a $\Lambda$CDM Universe
• The Layzer-Irvine Equation for Scalar-Tensor Theories: A Test of Modified Gravity N-body Simulations
• Evolution of density perturbations in a realistic universe

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