A revisited correction to the halo mass function for local-type primordial non-Gaussianity

We investigate the effect of primordial non-Gaussianities on halo number counts using N-body simulations with different values of f NLloc. We show how current theoretical models fail to adequately describe the non-Gaussian mass function of halos identified with different overdensity thresholds, Δb. We explain how these discrepancies are related to a variation in the density profile of dark matter halos, finding that the internal steepness (i.e. the compactness) of halos depends on the value of f NLloc. We then parametrize these deviations in halo number counts with a factor κ(Δb) that modifies the linear density threshold for collapse according to the halo identification threshold used, defined with respect to the Universe background density. We rely on a second-degree polynomial to describe κ and employ a Bayesian analysis to determine the coefficients of this polynomial. In addition, we verify the independence of the latter on the sign and absolute value of f NLloc. Finally, we show how this re-parametrization prevents the extraction of biased constraints on f NLloc, correcting for large systematic errors especially in the case of halos identified with high density thresholds. This improvement is crucial in the perspective of deriving cosmological constraints with the non-Gaussian mass function from real data, as different mass definitions can be employed depending on the properties of the survey.