Improving Initial Conditions for Cosmological $N$-Body Simulations

Abstract

In cosmological $N$-body simulations, the representation of dark matter as discrete "macroparticles" suppresses the growth of structure, such that simulations no longer reproduce linear theory on small scales near $k_{\mathrm{Nyquist}}$. Marcos et al. demonstrate that this is due to sparse sampling of modes near $k_{\mathrm{Nyquist}}$ and that the often-assumed continuum growing modes are not proper growing modes of the particle system. We develop initial conditions that respect the particle linear theory growing modes and then rescale the mode amplitudes to account for growth suppression. These ICs also allow us to take advantage of our very accurate $N$-body code Abacus to implement 2LPT in configuration space. The combination of 2LPT and rescaling improves the accuracy of the late-time power spectra, halo mass functions, and halo clustering. In particular, we achieve 1% accuracy in the power spectrum down to $k_{\mathrm{Nyquist}}$, versus $k_{\mathrm{Nyquist}}/4$ without rescaling or $k_{\mathrm{Nyquist}}/13$ without 2LPT, relative to an oversampled reference simulation. We anticipate that our 2LPT will be useful for large simulations where FFTs are expensive and that rescaling will be useful for suites of medium-resolution simulations used in cosmic emulators and galaxy survey mock catalogs. Code to generate initial conditions is available at https://github.com/lgarrison/zeldovich-PLT

22 pages, 14 figures. Minor changes to match MNRAS published version. Code available at https://github.com/lgarrison/zeldovich-PLT

References in corpus (7)

• Astrophysics Source Code Library
• Transients from Initial Conditions in Cosmological Simulations
• SPACE: the SPectroscopic All-sky Cosmic Explorer
• 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
• Quantification of discreteness effects in cosmological N-body simulations: II. Evolution up to shell crossing
• Particle mesh multipole method: An efficient solver for gravitational/electrostatic forces based on multipole method and fast convolution over a uniform mesh

Cited by in corpus (55)

• AbacusSummit: A Massive Set of High-Accuracy, High-Resolution $N$-Body Simulations
• Dark Quest. I. Fast and Accurate Emulation of Halo Clustering Statistics and Its Application to Galaxy Clustering
• Corrfunc --- A Suite of Blazing Fast Correlation Functions on the CPU
• Euclid preparation: II. The EuclidEmulator -- A tool to compute the cosmology dependence of the nonlinear matter power spectrum
• Blinded challenge for precision cosmology with large-scale structure: results from effective field theory for the redshift-space galaxy power spectrum
• Large-scale dark matter simulations
• Euclid preparation: IX. EuclidEmulator2 -- Power spectrum emulation with massive neutrinos and self-consistent dark energy perturbations
• The $\text{Abacus}$ Cosmological $N$-body Code
• Limitations to the "basic'' HOD model and beyond
• The Abacus Cosmos: A Suite of Cosmological N-body Simulations
• A High-Fidelity Realization of the Euclid Code Comparison $N$-body Simulation with Abacus
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• Initial Conditions for Accurate N-Body Simulations of Massive Neutrino Cosmologies
• The cosmology dependence of galaxy clustering and lensing from a hybrid $N$-body-perturbation theory model
• Evidence for galaxy assembly bias in BOSS CMASS redshift-space galaxy correlation function
• Emulating galaxy clustering and galaxy-galaxy lensing into the deeply nonlinear regime: methodology, information, and forecasts
• An iterative reconstruction of cosmological initial density fields
• Exploring the squeezed three-point galaxy correlation function with generalized halo occupation distribution models
• Aemulus $\nu$: Precise Predictions for Matter and Biased Tracer Power Spectra in the Presence of Neutrinos
• Higher-order initial conditions with massive neutrinos
• Cosmology with Stacked Cluster Weak Lensing and Cluster-Galaxy Cross-Correlations
• Moving around the cosmological parameter space: a nonlinear power spectrum reconstruction based on high-resolution cosmic responses
• The DESI $N$-body Simulation Project I: Testing the Robustness of Simulations for the DESI Dark Time Survey
• Accurate predictions from small boxes: variance suppression via the Zel'dovich approximation
• A fast semi-discrete optimal transport algorithm for a unique reconstruction of the early Universe
• Cosmological perturbations for two cold fluids in $\Lambda$CDM
• Some statistical remarks on the Giant GRB Ring
• Can Assembly Bias Explain the Lensing Amplitude of the BOSS CMASS Sample in a Planck Cosmology?
• The DESI $N$-body Simulation Project -- II. Suppressing sample variance with fast simulations
• Quantifying resolution in cosmological N-body simulations using self-similarity
• Accurate initial conditions in mixed Dark Matter--Baryon simulations
• Perturbation theory with dispersion and higher cumulants: non-linear regime
• Unequal time correlators and the Zeldovich approximation
• Application of the iterative reconstruction to simulated galaxy fields
• Good and Proper: Self-similarity of N-body Simulations with Proper Force Softening
• Using galaxy pairs to investigate the three-point correlation function in the squeezed limit
• Accuracy of power spectra in dissipationless cosmological simulations
• Constraining galaxy-halo connection with high-order statistics
• BullFrog: Multi-step perturbation theory as a time integrator for cosmological simulations
• Tests of Acoustic Scale Shifts in Halo-based Mock Galaxy Catalogues
• Numerical convergence of pre-initial conditions on dark matter halo properties
• Isolating the linear signal when making redshift space distortion measurements
• Perturbation-theory informed integrators for cosmological simulations
• Self-Similarity of $k$-Nearest Neighbor Distributions in Scale-Free Simulations
• Starting Cosmological Simulations from the Big Bang
• A fully Lagrangian, non-parametric bias model for dark-matter halos
• Convergence of halo statistics: code comparison between Rockstar and CompaSO using scale-free simulations
• Decorrelating the errors of the galaxy correlation function with compact transformation matrices
• Generating Approximate Halo Catalogs for Blind Challenges in Precision Cosmology
• The Effect of Corner Modes in the Initial Conditions of Cosmological Simulations
• The Hidden Role of Anisotropies in Shaping Structure Formation in Cosmological N-Body Simulations
• Constraining Galaxy-Halo Connection Using Machine Learning
• Non-parametric Lagrangian biasing from the insights of neural nets
• Particle loads for cosmological simulations with equal-mass dark matter and baryonic particles
• Tidal adaptive softening and artificial fragmentation in cosmological simulations