Phase-space simulations of prompt cusps: simulating the formation of the first haloes without artificial fragmentation

Abstract: The first generation of haloes forms from the collapse of the smallest peaks in the initial density field. $N$-body simulations of this process suggest a prompt formation of a steep power-law cusp, but these calculations are plagued by numerical artifacts which casts some doubt on this result. Here, we develop new simulation methods based on the dark matter phase-space sheet approach and present results which are entirely free of artificial clumps. We find that a cusp with density $\rho\propto r^{-1.5}$ is indeed formed promptly, subsequently accreting a more extended halo and participating in the hierarchical growth of later halo generations. However, our simulations also suggest that the presence of artificial clumps just before peak collapse can significantly shallow the inner profiles of the cusps. We use $N$-body simulations with controlled amounts of small-scale power to place a conservative upper limit on the scales affected by artificial clumps. Finally, we used these results to simulate the collapse of the first generation of peaks of various types and in different cosmologies, finding prompt cusps to form in all cases. We conclude that prompt cusps are a generic feature of the collapse of peaks on the free-streaming scale of the initial density field, and their structure can safely be studied using $N$-body simulations provided care is taken to excise the region potentially affected by artificial clumps.