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 DM : N gas = Ω cdm /Ω 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 DM : N 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 DM : N 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 ) ∝ 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 multicomponent uniform and isotropic distributions. Our code based on gadget-2 is available at https://github.com/liaoshong/gadget-2glass .