Quijote-PNG: Quasi-maximum likelihood estimation of Primordial Non-Gaussianity in the non-linear dark matter density field

Abstract

Future Large Scale Structure surveys are expected to improve over current bounds on primordial non-Gaussianity (PNG), with a significant impact on our understanding of early Universe physics. The level of such improvements will however strongly depend on the extent to which late time non-linearities erase the PNG signal on small scales. In this work, we show how much primordial information remains in the bispectrum of the non-linear dark matter density field by implementing a new, simulation-based, methodology for joint estimation of PNG amplitudes ($f_{\mathrm{NL}}$) and standard $\Lambda$CDM parameters. The estimator is based on optimally compressed statistics, which, for a given input density field, combine power spectrum and modal bispectrum measurements, and numerically evaluate their covariance and their response to changes in cosmological parameters. We train and validate the estimator using a large suite of N-body simulations (QUIJOTE-PNG), including different types of PNG (local, equilateral, orthogonal). We explicitly test the estimator's unbiasedness, optimality and stability with respect to changes in the total number of input realizations. While the dark matter power spectrum itself contains negligible PNG information, as expected, including it as an ancillary statistic increases the PNG information content extracted from the bispectrum by a factor of order $2$. As a result, we prove the capability of our approach to optimally extract PNG information on non-linear scales beyond the perturbative regime, up to $k_{\max }=0.5h{\mathrm{Mpc}}^{-1}$, obtaining marginalized $1$-$\sigma$ bounds of $\Delta f_{\mathrm{NL}}^{\mathrm{local}}\sim 16$, $\Delta f_{\mathrm{NL}}^{\mathrm{equil}}\sim 77$ and $\Delta f_{\mathrm{NL}}^{\mathrm{ortho}}\sim 40$ on a cubic volume of $1(\mathrm{Gpc}/h{)}^3$ at $z=1$. At the same time, we discuss the significant information on cosmological parameters contained on these scales.

22 pages, 12 figures

References in corpus (25)

• Why your model parameter confidences might be too optimistic -- unbiased estimation of the inverse covariance matrix
• The imprints of primordial non-gaussianities on large-scale structure: scale dependent bias and abundance of virialized objects
• Transients from Initial Conditions in Cosmological Simulations
• Primordial non-gaussianity, statistics of collapsed objects, and the Integrated Sachs-Wolfe effect
• Constraints on Multi-Field Inflation from the BOSS Galaxy Survey
• Nuisance hardened data compression for fast likelihood-free inference
• Constraining $M_{\nu }$ with the Bispectrum II: The Total Information Content of the Galaxy Bispectrum
• Constraints on Single-Field Inflation from the BOSS Galaxy Survey
• Precision analysis of the redshift-space galaxy bispectrum
• Primordial non-Gaussianity in the large scale structure of the Universe
• Joint analysis of anisotropic power spectrum, bispectrum and trispectrum: application to N-body simulations
• Maximal compression of the redshift space galaxy power spectrum and bispectrum
• CARPool: fast, accurate computation of large-scale structure statistics by pairing costly and cheap cosmological simulations
• The Persistence of Large Scale Structures I: Primordial non-Gaussianity
• Matter trispectrum: theoretical modelling and comparison to N-body simulations
• The Covariance of Squeezed Bispectrum Configurations
• Hunting high and low: Disentangling primordial and late-time non-Gaussianity with cosmic densities in spheres
• The Power of Locality: Primordial Non-Gaussianity at the Map Level
• Constraining Primordial Non-Gaussianity with Post-reconstructed Galaxy Bispectrum in Redshift Space
• CARPool Covariance: Fast, unbiased covariance estimation for large-scale structure observables
• Fisher Forecasts for Primordial non-Gaussianity from Persistent Homology
• Modal compression of the redshift-space galaxy bispectrum
• Towards cosmological constraints from the compressed modal bispectrum: a robust comparison of real-space bispectrum estimators
• Constraining Primordial Non-Gaussianity with Moments of the Large Scale Density Field
• Bayesian Control Variates for optimal covariance estimation with pairs of simulations and surrogates

Cited by in corpus (13)

• Quijote-PNG: Simulations of primordial non-Gaussianity and the information content of the matter field power spectrum and bispectrum
• GLASS: Generator for Large Scale Structure
• Cosmological constraints from the power spectrum and bispectrum of 21cm intensity maps
• Light Fields during Inflation from BOSS and Future Galaxy Surveys
• Taming assembly bias for primordial non-Gaussianity
• Improving constraints on primordial non-Gaussianity using neural network based reconstruction
• Quijote-PNG: The Information Content of the Halo Mass Function
• Quijote-PNG: Quasi-maximum likelihood estimation of Primordial Non-Gaussianity in the non-linear halo density field
• Constraining $f_{NL}$ using the Large-Scale Modulation of Small-Scale Statistics
• Squeezing information from radio surveys to probe the primordial Universe
• Bending the web: exploring the impact of modified gravity on the density field and halo properties within the cosmic web
• Power Spectrum Emulators from Neural Networks and Tree-Based Methods
• Quijote PNG: The information content of the halo power spectrum and bispectrum