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

TL;DR

This paper demonstrates a comprehensive Bayesian field-level inference method for galaxy surveys, capable of constraining primordial non-Gaussianity with high precision, accounting for systematics and leveraging all available large-scale structure information.

Contribution

It introduces a novel Bayesian hierarchical framework that integrates multiple probes and systematics for field-level inference of primordial non-Gaussianity in Euclid data.

Findings

Achieves forecasted precision of σ(f_NL) = 2.3 with Euclid data.

Incorporates formation history and all large-scale structure probes.

Handles survey systematics and observational effects effectively.

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.