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

This paper presents the PNG-UNITSIMS 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(NL)(local)=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.24x10(9) h(-1) M-circle dot, 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-UNITSIMS, we measure the PNG response parameter, p, as a function of the halo mass. halos with masses between 1x10(12) and 5x10(13) h(-1) M-circle dot are well described by the universality relation, given by p=1. For halos with masses between 2x10(10) and 1x10(12) h(-1) M-circle dot we find that p<1, at a significance between 1.5 and 3.1 sigma. Combining all the halos between 2x10(10) and 5x10(13) h(-1) M-circle dot, we find p consistent with a value of 0.955 +/- 0.013, which is 3 sigma 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-UNITSIMS 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 sigma biases on f(NL).