Quantification of high dimensional non-Gaussianities and its implication to Fisher analysis in cosmology

Abstract

It is well known that the power spectrum is not able to fully characterize the statistical properties of non-Gaussian density fields. Recently, many different statistics have been proposed to extract information from non-Gaussian cosmological fields that perform better than the power spectrum. The Fisher matrix formalism is commonly used to quantify the accuracy with which a given statistic can constrain the value of the cosmological parameters. However, these calculations typically rely on the assumption that the likelihood of the considered statistic follows a multivariate Gaussian distribution. In this work we follow Sellentin & Heavens (2017) and use two different statistical tests to identify non-Gaussianities in different statistics such as the power spectrum, bispectrum, marked power spectrum, and wavelet scatering transform (WST). We remove the non-Gaussian components of the different statistics and perform Fisher matrix calculations with the \textit{Gaussianized} statistics using Quijote simulations. We show that constraints on the parameters can change by a factor of $\sim 2$ in some cases. We show with simple examples how statistics that do not follow a multivariate Gaussian distribution can achieve artificially tight bounds on the cosmological parameters when using the Fisher matrix formalism. We think that the non-Gaussian tests used in this work represent a powerful tool to quantify the robustness of Fisher matrix calculations and their underlying assumptions. We release the code used to compute the power spectra, bispectra, and WST that can be run on both CPUs and GPUs.

arXiv comment24 pages, 6 figures

References in corpus (54)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Karhunen-Loeve eigenvalue problems in cosmology: how should we tackle large data sets?
• Why your model parameter confidences might be too optimistic -- unbiased estimation of the inverse covariance matrix
• The Quijote simulations
• The Bispectrum: From Theory to Observations
• Rejuvenating the matter power spectrum: restoring information with a logarithmic density mapping
• A new approach to observational cosmology using the scattering transform
• Constraining $M_{\nu }$ with the Bispectrum I: Breaking Parameter Degeneracies
• Nuisance hardened data compression for fast likelihood-free inference
• Towards an Optimal Estimation of Cosmological Parameters with the Wavelet Scattering Transform
• Fisher for complements: Extracting cosmology and neutrino mass from the counts-in-cells PDF
• Constraining $M_{\nu }$ with the Bispectrum II: The Total Information Content of the Galaxy Bispectrum
• Automatic physical inference with information maximising neural networks
• Nearest Neighbor distributions: new statistical measures for cosmological clustering
• Using the Marked Power Spectrum to Detect the Signature of Neutrinos in Large-Scale Structure
• On the assumption of Gaussianity for cosmological two-point statistics and parameter dependent covariance matrices
• New Interpretable Statistics for Large Scale Structure Analysis and Generation
• A marked correlation function for constraining modified gravity models
• Detecting Neutrino Mass by Combining Matter Clustering, Halos, and Voids
• Cosmic Shear Cosmology Beyond 2-Point Statistics: A Combined Peak Count and Correlation Function Analysis of DES-Y1
• Constraining neutrino mass with tomographic weak lensing one-point probability distribution function and power spectrum
• Constraining neutrino mass with weak lensing Minkowski Functionals
• Constraining Neutrino Mass with the Tomographic Weak Lensing Bispectrum
• Constraining neutrino mass with tomographic weak lensing peak counts
• The RWST, a comprehensive statistical description of the non-Gaussian structures in the ISM
• Weak lensing scattering transform: dark energy and neutrino mass sensitivity
• Primordial non-Gaussianity without tails -- how to measure fNL with the bulk of the density PDF
• Joint analysis of anisotropic power spectrum, bispectrum and trispectrum: application to N-body simulations
• On the insufficiency of arbitrarily precise covariance matrices: non-Gaussian weak lensing likelihoods
• A fast estimator for the bispectrum and beyond - A practical method for measuring non-Gaussianity in 21-cm maps
• Clipping the Cosmos: The Bias and Bispectrum of Large Scale Structure
• Matter trispectrum: theoretical modelling and comparison to N-body simulations
• Cosmological cross-correlations and nearest neighbor distributions
• Lossless, Scalable Implicit Likelihood Inference for Cosmological Fields
• What does the Marked Power Spectrum Measure? Insights from Perturbation Theory
• Classification of Magnetohydrodynamic Simulations using Wavelet Scattering Transforms
• Statistical description of dust polarized emission from the diffuse interstellar medium -- A RWST approach
• Information Content of Higher-Order Galaxy Correlation Functions
• Likelihood Non-Gaussianity in Large-Scale Structure Analyses
• Minkowski Functionals of Convergence Maps and the Lensing Figure of Merit
• Cosmology and neutrino mass with the Minimum Spanning Tree
• Exploring KSZ velocity reconstruction with $N$-body simulations and the halo model
• What Can We Learn by Combining the Skew Spectrum and the Power Spectrum?
• Clipping the Cosmos II: Cosmological information from non-linear scales
• The void halo mass function: a promising probe of neutrino mass
• Weighing neutrinos with the halo environment
• Probing massive neutrinos with the Minkowski functionals of large-scale structure
• Towards cosmological constraints from the compressed modal bispectrum: a robust comparison of real-space bispectrum estimators
• The unequal-time matter power spectrum: impact on weak lensing observables
• Flow-Based Likelihoods for Non-Gaussian Inference
• The Information Content of Projected Galaxy Fields
• Cosmology with the kinetic Sunyaev-Zeldovich effect: Independent of the optical depth and ${\sigma }_8$
• Constraining the Neutrino Mass with the Drifting Coefficient of the Field Cluster Mass Function

Cited by in corpus (9)

• Precise Cosmological Constraints from BOSS Galaxy Clustering with a Simulation-Based Emulator of the Wavelet Scattering Transform
• Galaxy Clustering Analysis with SimBIG and the Wavelet Scattering Transform
• EFTofLSS meets simulation-based inference: ${\sigma }_8$ from biased tracers
• Probing massive neutrinos with the Minkowski functionals of the galaxy distribution
• Wavelet Based Statistics for Enhanced 21cm EoR Parameter Constraints
• Towards unveiling the large-scale nature of gravity with the wavelet scattering transform
• Hitting the mark: Optimising Marked Power Spectra for Cosmology
• Neural Network Reconstruction of Non-Gaussian Initial Conditions from Dark Matter Halos
• Comparing the morphology of molecular clouds without supervision