The abundance of dark matter haloes down to Earth mass

Haonan Zheng, Sownak Bose, Carlos S. Frenk, Liang Gao, Adrian Jenkins, Shihong Liao, Yizhou Liu, Jie Wang

Published 2023-10-24Version 1

We use the Voids-within-Voids-within-Voids (VVV) simulations, a suite of successive nested N-body simulations with extremely high resolution (denoted, from low to high resolution, by L0 to L7), to test the Press-Schechter (PS), Sheth-Tormen (ST), and extended Press-Schechter (EPS) formulae for the halo abundance over the entire mass range, from mini-haloes of $10^{-6}\ \mathrm{M_\odot}$, to rich cluster haloes of $10^{15}\ \mathrm{M_\odot}$, at different redshifts, from $z=30$ to the present. We find that at $z=0$ and $z=2$, ST gives the best prediction for the results of L0, which has the mean cosmic density (i.e. overdensity $\delta=0$), at $10^{11-15} ~\mathrm{M_\odot}$. The higher resolution levels (L1-L7) are biased regions at various underdensities ($\delta<-0.6$). The EPS formalism takes this into account since it gives the mass function of a region conditioned, in this case, on having a given underdensity. The EPS gives good predictions for these higher levels, with deviations $\lesssim 20\%$, at $10^{-6-12.5} ~\mathrm{M_\odot}$. At $z \sim 7-15$, the ST predictions for L0 and the EPS predictions for L1-L7 show somewhat larger deviations from the simulation results. However, at even higher redshifts, $z \sim 30$, the EPS prediction fits the simulations well again. We further confirm our results by picking more subvolumes from the full L0 simulation, finding that our conclusions depend only weakly on the size and overdensity of the chosen region. Since at mean density the EPS reduces to the PS mass function, its good agreement with the higher-level simulations implies that the PS (or, even better, the ST) formula gives an accurate description of the total halo mass function in representative regions of the universe.