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.

arXiv comment22 pages, 12 figures

References in corpus (56)

• The cosmological simulation code GADGET-2
• Efficient Computation of CMB anisotropies in closed FRW models
• The Baryon Oscillation Spectroscopic Survey of SDSS-III
• Cosmology and Fundamental Physics with the Euclid Satellite
• Power Spectrum Analysis of Three-Dimensional Redshift Surveys
• 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
• Constraints on local primordial non-Gaussianity from large scale structure
• The effect of primordial non-Gaussianity on halo bias
• The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: RSD measurement from the power spectrum and bispectrum of the DR12 BOSS galaxies
• The Quijote simulations
• Massive Lossless Data Compression and Multiple Parameter Estimation from Galaxy Spectra
• Parameter inference with estimated covariance matrices
• Accurate Estimators of Correlation Functions in Fourier Space
• The power spectrum and bispectrum of SDSS DR11 BOSS galaxies I: bias and gravity
• Primordial non-gaussianity, statistics of collapsed objects, and the Integrated Sachs-Wolfe effect
• Large-scale Bias and Efficient Generation of Initial Conditions for Non-Local Primordial Non-Gaussianity
• General CMB and Primordial Bispectrum Estimation I: Mode Expansion, Map-Making and Measures of f_NL
• Constraining Primordial non-Gaussianity with Bispectrum and Power Spectum from Upcoming Optical and Radio Surveys
• Constraining $M_{\nu }$ with the Bispectrum I: Breaking Parameter Degeneracies
• Detection of Baryon Acoustic Oscillation Features in the Large-Scale 3-Point Correlation Function of SDSS BOSS DR12 CMASS Galaxies
• Detectability of the effect of Inflationary non-Gaussianity on halo bias
• 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
• The large-scale 3-point correlation function of the SDSS BOSS DR12 CMASS galaxies
• Generalized massive optimal data compression
• The CMB Bispectrum
• Precision analysis of the redshift-space galaxy bispectrum
• N-body simulations with generic non-Gaussian initial conditions I: Power Spectrum and halo mass function
• The power spectrum and bispectrum of SDSS DR11 BOSS galaxies II: cosmological interpretation
• 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
• Primordial non-Gaussianity without tails -- how to measure fNL with the bulk of the density PDF
• Rapid Separable Analysis of Higher Order Correlators in Large Scale Structure
• 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
• Enhancing BOSS bispectrum cosmological constraints with maximal compression
• The Persistence of Large Scale Structures I: Primordial non-Gaussianity
• Matter bispectrum of large-scale structure: Three-dimensional comparison between theoretical models and numerical simulations
• Matter trispectrum: theoretical modelling and comparison to N-body simulations
• Fast Estimation of Gravitational and Primordial Bispectra in Large Scale Structures
• The Covariance of Squeezed Bispectrum Configurations
• Hunting high and low: Disentangling primordial and late-time non-Gaussianity with cosmic densities in spheres
• Universal Non-Gaussian Initial Conditions for N-body Simulations
• The Power of Locality: Primordial Non-Gaussianity at the Map Level
• CARPool Covariance: Fast, unbiased covariance estimation for large-scale structure observables
• Constraining Primordial Non-Gaussianity with Post-reconstructed Galaxy Bispectrum in Redshift Space
• 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
• Matter bispectrum of large-scale structure with Gaussian and non-Gaussian initial conditions: Halo models, perturbation theory, and a three-shape model
• 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