MCMC generation of cosmological fields far beyond Gaussianity
arXiv:2107.05639 · doi:10.21105/astro.2107.05639
Abstract
Structure formation in our Universe creates non-Gaussian random fields that will soon be observed over almost the entire sky by the Euclid satellite, the Vera-Rubin observatory, and the Square Kilometre Array. An unsolved problem is how to analyze best such non-Gaussian fields, e.g. to infer the physical laws that created them. This problem could be solved if a parametric non-Gaussian sampling distribution for such fields were known, as this distribution could serve as likelihood during inference. We therefore create a sampling distribution for non-Gaussian random fields. Our approach is capable of handling strong non-Gaussianity, while perturbative approaches such as the Edgeworth expansion cannot. To imitate cosmological structure formation, we enforce our fields to be (i) statistically isotropic, (ii) statistically homogeneous, and (iii) statistically independent at large distances. We generate such fields via a Monte Carlo Markov Chain technique and find that even strong non-Gaussianity is not necessarily visible to the human eye. We also find that sampled marginals for pixel pairs have an almost generic Gauss-like appearance, even if the joint distribution of all pixels is markedly non-Gaussian. This apparent Gaussianity is a consequence of the high dimensionality of random fields. We conclude that vast amounts of non-Gaussian information can be hidden in random fields that appear nearly Gaussian in simple tests, and that it would be short-sighted not to try and extract it.
arXiv commentAccepted by the Open Journal of Astrophysics
References in corpus (25)
- Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results
- Non-Gaussian features of primordial fluctuations in single field inflationary models
- First results from the IllustrisTNG simulations: matter and galaxy clustering
- Karhunen-Loeve eigenvalue problems in cosmology: how should we tackle large data sets?
- Euclid Definition Study Report
- LSST Science Book, Version 2.0
- A Conceptual Introduction to Hamiltonian Monte Carlo
- KiDS+VIKING-450: Cosmic shear tomography with optical+infrared data
- The Atacama Cosmology Telescope: Cosmological parameters from three seasons of data
- Transients from Initial Conditions: A Perturbative Analysis
- Reionization and the Cosmic Dawn with the Square Kilometre Array
- Measurement of Cosmic Microwave Background Polarization Power Spectra from Two Years of BICEP Data
- Planck 2018 results. IX. Constraints on primordial non-Gaussianity
- Quantifying the non-Gaussianity in the EoR 21-cm signal through bispectrum
- Breaking the spell of Gaussianity: forecasting with higher order Fisher matrices
- Fast optimal CMB power spectrum estimation with Hamiltonian sampling
- Cosmological parameters, shear maps and power spectra from CFHTLenS using Bayesian hierarchical inference
- The 21cm bispectrum as a probe of non-Gaussianities due to X-ray heating
- Dark Energy Survey Year 1 Results: Joint Analysis of Galaxy Clustering, Galaxy Lensing, and CMB Lensing Two-point Functions
- A fast, always positive definite and normalizable approximation of non-Gaussian likelihoods
- Generalised Fisher Matrices
- On the use of the Edgeworth expansion in cosmology I: how to foresee and evade its pitfalls
- Optimizing parameter constraints: a new tool for Fisher matrix forecasts
- Non-Gaussian forecasts of weak lensing with and without priors
- Markov Chain Truncation for Doubly-Intractable Inference