Particle loads for cosmological simulations with equal-mass dark matter and baryonic particles

Abstract

Traditional cosmological hydrodynamical simulations usually assume equal-numbered but unequal-mass dark matter and baryonic particles, which can lead to spurious collisional heating due to energy equipartition. To avoid such a numerical heating effect, a simulation setup with equal-mass dark matter and baryonic particles, which corresponds to a particle number ratio of $N_{\mathrm{DM}}:N_{\mathrm{gas}}={\Omega }_{\mathrm{cdm}}/{\Omega }_{\mathrm{b}}$, is preferred. However, previous studies have typically used grid-based particle loads to prepare such initial conditions, which can only reach specific values for $N_{\mathrm{DM}}:N_{\mathrm{gas}}$ due to symmetry requirements. In this study, we propose a method based on the glass approach that can generate two-component particle loads with more general $N_{\mathrm{DM}}:N_{\mathrm{gas}}$ ratios. The method simultaneously relaxes two Poisson particle distributions by introducing an additional repulsive force between particles of the same component. We show that the final particle load closely follows the expected minimal power spectrum, $P(k)\propto k^4$, exhibits good homogeneity and isotropy properties, and remains sufficiently stable under gravitational interactions. Both the dark matter and gas components individually also exhibit uniform and isotropic distributions. We apply our method to two-component cosmological simulations and demonstrate that an equal-mass particle setup effectively mitigates the spurious collisional heating that arises in unequal-mass simulations. Our method can be extended to generate multi-component uniform and isotropic distributions. Our code based on Gadget-2 is available at https://github.com/liaoshong/gadget-2glass .

16 pages, 11 figures, accepted for publication in The Astrophysical Journal. The code is publicly available at https://github.com/liaoshong/gadget-2glass

References in corpus (36)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• The cosmological simulation code GADGET-2
• Simulating the joint evolution of quasars, galaxies and their large-scale distribution
• Stable clustering, the halo model and nonlinear cosmological power spectra
• Scaling relations for galaxy clusters in the Millennium-XXL simulation
• Discreteness effects in simulations of Hot/Warm dark matter
• Large-scale dark matter simulations
• A New Parallel N-body Gravity Solver: TPM
• Energy equipartition between stellar and dark matter particles in cosmological simulations results in spurious growth of galaxy sizes
• A comparison of the evolution of the density field in perturbation theory and numerical simulations - II Counts in cells analysis
• Improving Initial Conditions for Cosmological $N$-Body Simulations
• Spurious heating of stellar motions in simulated galactic disks by dark matter halo particles
• Higher-order initial conditions for mixed baryon-CDM simulations
• Linear perturbative theory of the discrete cosmological N-body problem
• How closely do baryons follow dark matter on large scales?
• Towards quantitative control on discreteness error in the non-linear regime of cosmological N body simulations
• The Evolution of Baryon Density Fluctuations in Multi-Component Cosmological Simulations
• Fast Estimation of Gravitational and Primordial Bispectra in Large Scale Structures
• The impact of spurious collisional heating on the morphological evolution of simulated galactic discs
• Effects of the initial conditions on cosmological $N$-body simulations
• Galaxy Formation: Warm Dark Matter, Missing Satellites, and the Angular Momentum Problem
• Quantification of discreteness effects in cosmological N-body simulations: II. Evolution up to shell crossing
• Gravitational evolution of a perturbed lattice and its fluid limit
• More accurate simulations with separate initial conditions for baryons and dark matter
• The segregation of baryons and dark matter during halo assembly
• Accurate initial conditions in mixed Dark Matter--Baryon simulations
• An alternative method to generate pre-initial conditions for cosmological $N$-body simulations
• An alternative to grids and glasses: Quaquaversal pre-initial conditions for N-body simulations
• Modelling the accretion and feedback of supermassive black hole binaries in gas-rich galaxy mergers
• A high-resolution cosmological simulation of a strong gravitational lens
• Numerical convergence of pre-initial conditions on dark matter halo properties
• Numerical Discreteness Errors in Multi-Species Cosmological N-body Simulations
• Impacts of pre-initial conditions on anisotropic separate universe simulations: a boosted tidal response in the epoch of reionization
• N-body simulations of dark matter-baryon interactions
• On the origin of transient features in cosmological N-Body Simulations