PNG-UNITsims: Halo clustering response to primordial non-Gaussianities as a function of mass

This paper presents the PNG-UNIT SIMS suite, which includes the largest full N -body simulation to date with local primordial non-Gaussianities (local PNG), the PNG-UNIT. The amplitude of the PNGs is given by f local NL =100. The simulation follows the evolution of 4096 3 particles in a periodic box with L box = 1 h −1 Gpc, resulting in a mass resolution of m p = 1.24 × 10 9 h −1 M ⊙ , enough to finely resolve the galaxies targeted by stage-IV spectroscopic surveys. The PNG-UNIT has fixed initial conditions with phases also matching the pre-existing UNIT simulation with Gaussian initial conditions. The fixed and matched initial conditions reduce the simulation uncertainty significantly. In this first study of the PNG-UNIT SIMS , we measure the PNG response parameter, p , as a function of the halo mass. halos with masses between 1 × 10 12 and 5 × 10 13 h −1 M ⊙ are well described by the universality relation, given by p = 1. For halos with masses between 2 × 10 10 and 1 × 10 12 h −1 M ⊙ we find that p < 1, at a significance between 1.5 and 3.1 σ . Combining all the halos between 2 × 10 10 and 5 × 10 13 h −1 M ⊙ , we find p consistent with a value of 0.955 ± 0.013, which is 3 σ away from the universality relation. We demonstrate that these findings are robust to mass resolution, scale cuts and uncertainty estimation. We also compare our measurements to separate universe simulations, finding that the PNG-UNIT SIMS constraints outperform the former for the setup considered. Using a prior on p as tight as the one reported here for DESI-like forecast can result in f NL constraints comparable to fixing p . At the same time, fixing p to a wrong value ( p = 1) may result in up to 2 σ biases on f NL .