Measurements of the non-Gaussianity of the primordial density field have the power to considerably improve our understanding of the physics of inflation. Indeed, if we can increase the precision of current measurements by an order of magnitude, a null-detection would rule out many classes of scenarios for generating primordial fluctuations. Large-scale galaxy redshift surveys represent experiments that hold the promise to realise this goal. Thus, we model the galaxy bispectrum and forecast the accuracy with which it will probe the parameter $f_{\rm NL}$, which represents the degree of primordial local-type non Gaussianity. Specifically, we address the problem of modelling redshift space distortions (RSD) in the tree-level galaxy bispectrum including $f_{\rm NL}$. We find novel contributions associated with RSD, with the characteristic large scale amplification induced by local-type non-Gaussianity. These RSD effects must be properly accounted for in order to obtain un-biased measurements of $f_{\rm NL}$ from the galaxy bispectrum. We propose an analytic template for the monopole which can be used to fit against data on large scales, extending models used in the recent measurements. Finally, we perform idealised forecasts on $\sigma_{f_{\rm NL}}$ -- the accuracy of the determination of local non-linear parameter $f_{\rm NL}$ -- from measurements of the galaxy bispectrum. Our findings suggest that current surveys can in principle provide $f_{\rm NL}$ constraints competitive with Planck, and future surveys could improve them further.

Primordial non-Gaussianity can lead to a scale-dependent bias in the density of collapsed halos relative to the underlying matter density. The galaxy power spectrum already provides constraints on local-type primordial non-Gaussianity complementary those from the cosmic microwave background (CMB), while the bispectrum contains additional shape information and has the potential to outperform CMB constraints in future. We develop the bias model for the halo density contrast in the presence of local-type primordial non-Gaussianity, deriving a bivariate expansion up to second order in terms of the local linear matter density contrast and the local gravitational potential in Lagrangian coordinates. Non-linear evolution of the matter density introduces a non-local tidal term in the halo model. Furthermore, the presence of local-type non-Gaussianity in the Lagrangian frame leads to a novel non-local convective term in the Eulerian frame, that is proportional to the displacement field when going beyond the spherical collapse approximation. We use an extended Press-Schechter approach to evaluate the halo mass function and thus the halo bispectrum. We show that including these non-local terms in the halo bispectra can lead to corrections of up to $25\%$ for some configurations, on large scales or at high redshift.

The statistics of large-scale structure in the Universe can be used to probe non-Gaussianity of the primordial density field, complementary to existing constraints from the cosmic microwave background. In particular, the scale dependence of halo bias, which affects the halo distribution at large scales, represents a promising tool for analyzing primordial non-Gaussianity of local form. Future observations, for example, may be able to constrain the trispectrum parameter gNL that is difficult to study and constrain using the CMB alone. We investigate how galaxy and matter bispectra can distinguish between the two non-Gaussian parameters fNL and gNL, whose effects give nearly degenerate contributions to the power spectra. We use a generalization of the univariate bias approach, making the hypothesis that the number density of halos forming at a given position is a function of the local matter density contrast and of its local higher-order statistics. Using this approach, we calculate the halo-matter bispectra and analyze their properties. We determine a connection between the sign of the halo bispectrum on large scales and the parameter gNL. We also construct a combination of halo and matter bispectra that is sensitive to fNL, with little contamination from gNL. We study both the case of single and multiple sources to the primordial gravitational potential, discussing how to extend the concept of stochastic halo bias to the case of bispectra. We use a specific halo mass-function to calculate numerically the bispectra in appropriate squeezed limits, confirming our theoretical findings.