Quijote-PNG: Quasi-maximum Likelihood Estimation of Primordial Non-Gaussianity in the Nonlinear Halo Density Field
We study primordial non-Gaussian signatures in the redshift-space halo field on nonlinear scales, using a quasi-maximum likelihood estimator based on optimally compressed power spectrum and modal bispectrum statistics. We train and validate the estimator on a suite of halo catalogs constructed from...
Main Authors: | , , , , , , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
IOP Publishing
2023-01-01
|
Series: | The Astrophysical Journal |
Subjects: | |
Online Access: | https://doi.org/10.3847/1538-4357/acc4bd |
_version_ | 1797697041956077568 |
---|---|
author | Gabriel Jung Dionysios Karagiannis Michele Liguori Marco Baldi William R. Coulton Drew Jamieson Licia Verde Francisco Villaescusa-Navarro Benjamin D. Wandelt |
author_facet | Gabriel Jung Dionysios Karagiannis Michele Liguori Marco Baldi William R. Coulton Drew Jamieson Licia Verde Francisco Villaescusa-Navarro Benjamin D. Wandelt |
author_sort | Gabriel Jung |
collection | DOAJ |
description | We study primordial non-Gaussian signatures in the redshift-space halo field on nonlinear scales, using a quasi-maximum likelihood estimator based on optimally compressed power spectrum and modal bispectrum statistics. We train and validate the estimator on a suite of halo catalogs constructed from the Quijote-png N -body simulations, which we release to accompany this paper. We verify its unbiasedness and near-optimality for the three main types of primordial non-Gaussianity (PNG): local, equilateral, and orthogonal. We compare the modal bispectrum expansion with a k -binning approach, showing that the former allows for faster convergence of numerical derivatives in the computation of the score function, thus leading to better final constraints. We find, in agreement with previous studies, that the local PNG signal in the halo field is dominated by the scale-dependent bias signature on large scales and saturates at k ∼ 0.2 h Mpc ^−1 , whereas the small-scale bispectrum is the main source of information for equilateral and orthogonal PNG. Combining the power spectrum and bispectrum on nonlinear scales plays an important role in breaking degeneracies between cosmological and PNG parameters; such degeneracies, however, remain strong for equilateral PNG. We forecast that PNG parameters can be constrained with ${\rm{\Delta }}{f}_{\mathrm{NL}}^{\mathrm{local}}=45$ , ${\rm{\Delta }}{f}_{\mathrm{NL}}^{\mathrm{equil}}=570$ , and ${\rm{\Delta }}{f}_{\mathrm{NL}}^{\mathrm{ortho}}=110$ on a cubic volume of $1{\left(\mathrm{Gpc}\,{h}^{-1}\right)}^{3}$ at z = 1, considering scales up to ${k}_{\max }=0.5\,h\,{\mathrm{Mpc}}^{-1}$ . |
first_indexed | 2024-03-12T03:34:37Z |
format | Article |
id | doaj.art-4649ebe62a9a446398ab15cc78b03d87 |
institution | Directory Open Access Journal |
issn | 1538-4357 |
language | English |
last_indexed | 2024-03-12T03:34:37Z |
publishDate | 2023-01-01 |
publisher | IOP Publishing |
record_format | Article |
series | The Astrophysical Journal |
spelling | doaj.art-4649ebe62a9a446398ab15cc78b03d872023-09-03T13:20:01ZengIOP PublishingThe Astrophysical Journal1538-43572023-01-01948213510.3847/1538-4357/acc4bdQuijote-PNG: Quasi-maximum Likelihood Estimation of Primordial Non-Gaussianity in the Nonlinear Halo Density FieldGabriel Jung0https://orcid.org/0000-0002-2133-6254Dionysios Karagiannis1https://orcid.org/0000-0002-4927-0816Michele Liguori2https://orcid.org/0000-0002-0383-2254Marco Baldi3https://orcid.org/0000-0003-4145-1943William R. Coulton4https://orcid.org/0000-0002-1297-3673Drew Jamieson5https://orcid.org/0000-0001-5044-7204Licia Verde6https://orcid.org/0000-0003-2601-8770Francisco Villaescusa-Navarro7https://orcid.org/0000-0002-4816-0455Benjamin D. Wandelt8https://orcid.org/0000-0002-5854-8269Dipartimento di Fisica e Astronomia “G. Galilei,” Università degli Studi di Padova , Via Marzolo 8, I-35131, Padova, Italy; INFN, Sezione di Padova, Via Marzolo 8, I-35131, Padova, ItalyDepartment of Physics and Astronomy, University of the Western Cape , Cape Town 7535, South AfricaDipartimento di Fisica e Astronomia “G. Galilei,” Università degli Studi di Padova , Via Marzolo 8, I-35131, Padova, Italy; INFN, Sezione di Padova, Via Marzolo 8, I-35131, Padova, ItalyDipartimento di Fisica e Astronomia, Alma Mater Studiorum—University of Bologna , Via Piero Gobetti 93/2, I-40129 Bologna BO, Italy; INAF—Osservatorio Astronomico di Bologna , Via Piero Gobetti 93/3, I-40129 Bologna BO, Italy; Istituto Nazionale di Fisica Nucleare , Sezione di Bologna, Viale Berti Pichat 6/2, I-40127 Bologna BO, ItalyCenter for Computational Astrophysics, Flatiron Institute , 162 5th Avenue, New York, NY 10010, USAMax-Planck-Institut für Astrophysik , Karl-Schwarzschild-Straße 1, D-85748 Garching, GermanyICREA, Pg. Lluís Companys 23, Barcelona, E-08010, Spain; ICC, University of Barcelona , IEEC-UB, Martí i Franquès, 1, E-08028 Barcelona, SpainCenter for Computational Astrophysics, Flatiron Institute , 162 5th Avenue, New York, NY 10010, USA; Department of Astrophysical Sciences, Princeton University , 4 Ivy Lane, Princeton, NJ 08544, USACenter for Computational Astrophysics, Flatiron Institute , 162 5th Avenue, New York, NY 10010, USA; Sorbonne Université , CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98 bis boulevard Arago, F-75014 Paris, FranceWe study primordial non-Gaussian signatures in the redshift-space halo field on nonlinear scales, using a quasi-maximum likelihood estimator based on optimally compressed power spectrum and modal bispectrum statistics. We train and validate the estimator on a suite of halo catalogs constructed from the Quijote-png N -body simulations, which we release to accompany this paper. We verify its unbiasedness and near-optimality for the three main types of primordial non-Gaussianity (PNG): local, equilateral, and orthogonal. We compare the modal bispectrum expansion with a k -binning approach, showing that the former allows for faster convergence of numerical derivatives in the computation of the score function, thus leading to better final constraints. We find, in agreement with previous studies, that the local PNG signal in the halo field is dominated by the scale-dependent bias signature on large scales and saturates at k ∼ 0.2 h Mpc ^−1 , whereas the small-scale bispectrum is the main source of information for equilateral and orthogonal PNG. Combining the power spectrum and bispectrum on nonlinear scales plays an important role in breaking degeneracies between cosmological and PNG parameters; such degeneracies, however, remain strong for equilateral PNG. We forecast that PNG parameters can be constrained with ${\rm{\Delta }}{f}_{\mathrm{NL}}^{\mathrm{local}}=45$ , ${\rm{\Delta }}{f}_{\mathrm{NL}}^{\mathrm{equil}}=570$ , and ${\rm{\Delta }}{f}_{\mathrm{NL}}^{\mathrm{ortho}}=110$ on a cubic volume of $1{\left(\mathrm{Gpc}\,{h}^{-1}\right)}^{3}$ at z = 1, considering scales up to ${k}_{\max }=0.5\,h\,{\mathrm{Mpc}}^{-1}$ .https://doi.org/10.3847/1538-4357/acc4bdNon-GaussianityCosmological parameters from large-scale structureFisher’s Information |
spellingShingle | Gabriel Jung Dionysios Karagiannis Michele Liguori Marco Baldi William R. Coulton Drew Jamieson Licia Verde Francisco Villaescusa-Navarro Benjamin D. Wandelt Quijote-PNG: Quasi-maximum Likelihood Estimation of Primordial Non-Gaussianity in the Nonlinear Halo Density Field The Astrophysical Journal Non-Gaussianity Cosmological parameters from large-scale structure Fisher’s Information |
title | Quijote-PNG: Quasi-maximum Likelihood Estimation of Primordial Non-Gaussianity in the Nonlinear Halo Density Field |
title_full | Quijote-PNG: Quasi-maximum Likelihood Estimation of Primordial Non-Gaussianity in the Nonlinear Halo Density Field |
title_fullStr | Quijote-PNG: Quasi-maximum Likelihood Estimation of Primordial Non-Gaussianity in the Nonlinear Halo Density Field |
title_full_unstemmed | Quijote-PNG: Quasi-maximum Likelihood Estimation of Primordial Non-Gaussianity in the Nonlinear Halo Density Field |
title_short | Quijote-PNG: Quasi-maximum Likelihood Estimation of Primordial Non-Gaussianity in the Nonlinear Halo Density Field |
title_sort | quijote png quasi maximum likelihood estimation of primordial non gaussianity in the nonlinear halo density field |
topic | Non-Gaussianity Cosmological parameters from large-scale structure Fisher’s Information |
url | https://doi.org/10.3847/1538-4357/acc4bd |
work_keys_str_mv | AT gabrieljung quijotepngquasimaximumlikelihoodestimationofprimordialnongaussianityinthenonlinearhalodensityfield AT dionysioskaragiannis quijotepngquasimaximumlikelihoodestimationofprimordialnongaussianityinthenonlinearhalodensityfield AT micheleliguori quijotepngquasimaximumlikelihoodestimationofprimordialnongaussianityinthenonlinearhalodensityfield AT marcobaldi quijotepngquasimaximumlikelihoodestimationofprimordialnongaussianityinthenonlinearhalodensityfield AT williamrcoulton quijotepngquasimaximumlikelihoodestimationofprimordialnongaussianityinthenonlinearhalodensityfield AT drewjamieson quijotepngquasimaximumlikelihoodestimationofprimordialnongaussianityinthenonlinearhalodensityfield AT liciaverde quijotepngquasimaximumlikelihoodestimationofprimordialnongaussianityinthenonlinearhalodensityfield AT franciscovillaescusanavarro quijotepngquasimaximumlikelihoodestimationofprimordialnongaussianityinthenonlinearhalodensityfield AT benjamindwandelt quijotepngquasimaximumlikelihoodestimationofprimordialnongaussianityinthenonlinearhalodensityfield |