An alternative method to generate pre-initial conditions for cosmological $N$-body simulations

Abstract

Currently, grid and glass methods are the two most popular choices to generate uniform particle distributions (i.e., pre-initial conditions) for cosmological $N$-body simulations. In this article, we introduce an alternative method called the capacity constrained Voronoi tessellation (CCVT), which originates from computer graphics. As a geometrical equilibrium state, a CCVT particle configuration satisfies two constraints: (i) the volume of the Voronoi cell associated with each particle is equal; (ii) every particle is in the center-of-mass position of its Voronoi cell. We show that the CCVT configuration is uniform and isotropic, follows perfectly the minimal power spectrum, $P(k)\propto k^4$, and is quite stable under gravitational interactions. It is natural to incorporate periodic boundary conditions during CCVT making, therefore, we can obtain a larger CCVT by tiling with a small periodic CCVT. When applying the CCVT pre-initial condition to cosmological $N$-body simulations, we show that it plays as good as grid and glass schemes. The CCVT method will be helpful in studying the numerical convergence of pre-initial conditions in cosmological simulations. It can also be used to set up pre-initial conditions in smoothed-particle hydrodynamics simulations.

11 pages, 10 figures, minor revision on figures to match the version accepted by MNRAS, a code to generate CCVT pre-initial conditions is publicly available at https://github.com/liaoshong/ccvt-preic

References in corpus (9)

• Transients from Initial Conditions in Cosmological Simulations
• Dark matter and cosmic structure
• ZOBOV: a parameter-free void-finding algorithm
• Bayesian Strong Gravitational-Lens Modeling on Adaptive Grids: Objective Detection of Mass Substructure in Galaxies
• Discreteness effects in simulations of Hot/Warm dark matter
• VOBOZ: An Almost-Parameter-Free Halo-Finding Algorithm
• The Voronoi Tessellation cluster finder in 2+1 dimensions
• Towards quantitative control on discreteness error in the non-linear regime of cosmological N body simulations
• An alternative to grids and glasses: Quaquaversal pre-initial conditions for N-body simulations

Cited by in corpus (19)

• Large-scale dark matter simulations
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• A Fully Self-Consistent Cosmological Simulation Pipeline for Interacting Dark Energy Models
• Dark Matter Halos in Interacting Dark Energy Models: Formation History, Density Profile, Spin and Shape
• More accurate simulations with separate initial conditions for baryons and dark matter
• On the Cosmic Web Elongation in Fuzzy Dark Matter Cosmologies: Effects on Density Profiles, Shapes and Alignments of Halos
• Constraining interacting dark energy models with the halo concentration-mass relation
• Testing a quintessence model with Yukawa interaction from cosmological observations and N-body simulations
• A study on Cubic Galileon Gravity Using N-body Simulations
• N-body simulations, halo mass functions, and halo density profile in $f(T)$ gravity
• Improved two-point correlation function estimates using glass-like distributions as a reference sample
• Numerical convergence of pre-initial conditions on dark matter halo properties
• On the statistical theory of self-gravitating collisionless dark matter flow
• Void Lensing in Cubic Galileon Gravity
• Impacts of pre-initial conditions on anisotropic separate universe simulations: a boosted tidal response in the epoch of reionization
• Probing the Type 3 interacting dark-energy model using matter pairwise velocity
• Prediction of Individual Halo Concentrations Across Cosmic Time Using Neural Networks
• Halo spin and orientation in Interacting Dark Matter Dark Energy Cosmology
• Particle loads for cosmological simulations with equal-mass dark matter and baryonic particles