Predictions for local PNG bias in the galaxy power spectrum and bispectrum and the consequences for f NL constraints

We use hydrodynamical separate universe simulations with the IllustrisTNG model to predict the local primordial non-Gaussianity (PNG) bias parameters b _ϕ and b _ϕδ, which enter at leading order in the galaxy power spectrum and bispectrum. This is the first time that b _ϕδ is measured from either gravity-only or galaxy formation simulations. For dark matter halos, the popular assumption of universality overpredicts the b _ϕδ(b ₁) relation in the range 1 ≲ b ₁ ≲ 3 by up to Δ b _ϕδ ~ 3 (b ₁ is the linear density bias). The adequacy of the universality relation is worse for the simulated galaxies, with the relations b _ϕ(b ₁) and b _ϕδ(b ₁) being generically redshift-dependent and very sensitive to how galaxies are selected (we test total, stellar and black hole mass, black hole mass accretion rate and color). The uncertainties on b _ϕ and b _ϕδ have a direct, often overlooked impact on the constraints of the local PNG parameter f _NL, which we study and discuss. For a survey with V = 100 Gpc³/h³ at z=1, uncertainties Δ b _ϕ ≲ 1 and Δ b _ϕδ ≲ 5 around values close to the fiducial can yield relatively unbiased constraints on f _NL using power spectrum and bispectrum data. We also show why priors on galaxy bias are useful even in analyses that fit for products f _NL b _ϕ and f _NL b _ϕδ. The strategies we discuss to deal with galaxy bias uncertainties can be straightforwardly implemented in existing f _NL constraint analyses (we provide fits for some of the bias relations). Our results motivate more works with galaxy formation simulations to refine our understanding of b _ϕ and b _ϕδ towards improved constraints on f _NL.