Euclid: Field-level inference of primordial non-Gaussianity and cosmic initial conditions

Abstract: A primary target of the \Euclid space mission is to constrain early-universe physics by searching for deviations from a primordial Gaussian random field. A significant detection of primordial non-Gaussianity would rule out the simplest models of cosmic inflation and transform our understanding of the origin of the Universe. This paper forecasts how well field-level inference of galaxy redshift surveys can constrain the amplitude of local primordial non-Gaussianity ($f_{NL}$), within a Bayesian hierarchical framework, in the upcoming \Euclid data. We design and simulate mock data sets and perform Markov chain Monte Carlo analyses using a full-field forward modelling approach. By including the formation history of the cosmic matter field in the analysis, the method takes into account all available probes of primordial non-Gaussianity, and goes beyond statistical summary estimators of $f_{NL}$. Probes include, for example, two-point and higher-order statistics, peculiar velocity fields, and scale-dependent galaxy biases. Furthermore, the method simultaneously handles systematic survey effects, such as selection effects, survey geometries, and galaxy biases. The forecast shows that the method can reach precision levels of up to $\sigma (f_{NL}) = 2.3$ (68.3\% CI, and at the grid resolution $\Delta L = 62.5\,h^{-1}$Mpc) with \Euclid data. We also provide data products, including realistic $N$-body simulations with nonzero values of $f_{NL}$ and maps of adiabatic curvature fluctuations. The results underscore the feasibility and advantages of field-level inference to constrain $f_{NL}$ in galaxy redshift surveys. Our approach consistently captures all the information available in the large-scale structure to constrain $f_{NL}$, and resolves the degeneracy between early-universe physics and late-time gravitational effects, while mitigating the impact of systematic and observational effects.