Primordial power spectrum and cosmology from black-box galaxy surveys

Abstract

We propose a new, likelihood-free approach to inferring the primordial matter power spectrum and cosmological parameters from arbitrarily complex forward models of galaxy surveys where all relevant statistics can be determined from numerical simulations, i.e. black-boxes. Our approach, which we call simulator expansion for likelihood-free inference (SELFI), builds upon approximate Bayesian computation using a novel effective likelihood, and upon the linearisation of black-box models around an expansion point. Consequently, we obtain simple "filter equations" for an effective posterior of the primordial power spectrum, and a straightforward scheme for cosmological parameter inference. We demonstrate that the workload is computationally tractable, fixed a priori, and perfectly parallel. As a proof of concept, we apply our framework to a realistic synthetic galaxy survey, with a data model accounting for physical structure formation and incomplete and noisy galaxy observations. In doing so, we show that the use of non-linear numerical models allows the galaxy power spectrum to be safely fitted up to at least $k_{\max }=0.5$ $h$/Mpc, outperforming state-of-the-art backward-modelling techniques by a factor of $\sim 5$ in the number of modes used. The result is an unbiased inference of the primordial matter power spectrum across the entire range of scales considered, including a high-fidelity reconstruction of baryon acoustic oscillations. It translates into an unbiased and robust inference of cosmological parameters. Our results pave the path towards easy applications of likelihood-free simulation-based inference in cosmology. We have made our code pySELFI and our data products publicly available at http://pyselfi.florent-leclercq.eu.

arXiv comment22 pages, 12 figures. Matches MNRAS published version. The main results are equations (25), (26), and figure 10. The pySELFI code is publicly available on the first author's website (currently http://pyselfi.florent-leclercq.eu) and on GitHub (https://github.com/florent-leclercq/pyselfi/)

References in corpus (50)

• emcee: The MCMC Hammer
• Planck 2015 results. XIII. Cosmological parameters
• Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC
• The cosmological simulation code GADGET-2
• Baryonic Features in the Matter Transfer Function
• The Clustering of the SDSS DR7 Main Galaxy Sample I: A 4 per cent Distance Measure at z=0.15
• Planck 2018 results. I. Overview and the cosmological legacy of Planck
• The Galaxy Luminosity Function and Luminosity Density at Redshift z=0.1
• Euclid Definition Study Report
• Why your model parameter confidences might be too optimistic -- unbiased estimation of the inverse covariance matrix
• Measuring the Baryon Acoustic Oscillation scale using the SDSS and 2dFGRS
• A 2% Distance to z=0.35 by Reconstructing Baryon Acoustic Oscillations - I : Methods and Application to the Sloan Digital Sky Survey
• The 2dF Galaxy Redshift Survey: The bias of galaxies and the density of the Universe
• Improving Cosmological Distance Measurements by Reconstruction of the Baryon Acoustic Peak
• On the Robustness of the Acoustic Scale in the Low-Redshift Clustering of Matter
• Solving Large Scale Structure in Ten Easy Steps with COLA
• Bayesian physical reconstruction of initial conditions from large scale structure surveys
• The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Analysis of potential systematics
• Baryonic Signatures in Large-Scale Structure
• Parameter inference with estimated covariance matrices
• Global, Exact Cosmic Microwave Background Data Analysis Using Gibbs Sampling
• Cosmological $N$-body simulations with suppressed variance
• Massive optimal data compression and density estimation for scalable, likelihood-free inference in cosmology
• ELUCID - Exploring the Local Universe with reConstructed Initial Density field I: Hamiltonian Markov Chain Monte Carlo Method with Particle Mesh Dynamics
• Physical Bayesian modelling of the non-linear matter distribution: new insights into the Nearby Universe
• Past and present cosmic structure in the SDSS DR7 main sample
• Bayesian power-spectrum inference for Large Scale Structure data
• Towards optimal extraction of cosmological information from nonlinear data
• Reconstruction within the Zeldovich approximation
• Dark Energy Survey Year 1 Results: Curved-Sky Weak Lensing Mass Map
• Approximate Bayesian Computation for Forward Modeling in Cosmology
• Reconstructing the Initial Density Field of the Local Universe: Method and Test with Mock Catalogs
• Unmasking the Masked Universe: the 2M++ catalogue through Bayesian eyes
• Reconstruction in Fourier space
• cosmoabc: Likelihood-free inference via Population Monte Carlo Approximate Bayesian Computation
• Bayesian optimisation for likelihood-free cosmological inference
• Reconstructing cosmological initial conditions from galaxy peculiar velocities. I. Reverse Zeldovich Approximation
• astroABC: An Approximate Bayesian Computation Sequential Monte Carlo sampler for cosmological parameter estimation
• Statistical properties of paired fixed fields
• Bayesian analysis of the dynamic cosmic web in the SDSS galaxy survey
• Methods for Bayesian power spectrum inference with galaxy surveys
• Probing deviations from General Relativity with the Euclid spectroscopic survey
• Matrix-free Large Scale Bayesian inference in cosmology
• The VIMOS Public Extragalactic Redshift Survey (VIPERS): Reconstruction of the redshift-space galaxy density field
• A Bayesian method for combining theoretical and simulated covariance matrices for large-scale structure surveys
• Bayesian cosmic density field inference from redshift space dark matter maps
• Bayesian power-spectrum inference with foreground and target contamination treatment
• Towards optimal estimation of the galaxy power spectrum
• Parameter inference and model comparison using theoretical predictions from noisy simulations
• Introductory Overview of Eulerian and Lagrangian Perturbation Theories

Cited by in corpus (17)

• AI-assisted super-resolution cosmological simulations
• Fast and Credible Likelihood-Free Cosmology with Truncated Marginal Neural Ratio Estimation
• ENCORE: An $\mathcal{O}(N_{\mathrm{g}}^2)$ Estimator for Galaxy $N$-Point Correlation Functions
• Information Content of Higher-Order Galaxy Correlation Functions
• Kernel-Based Emulator for the 3D Matter Power Spectrum from CLASS
• Baryon acoustic oscillations reconstruction using convolutional neural networks
• On the accuracy and precision of correlation functions and field-level inference in cosmology
• Systematic-free inference of the cosmic matter density field from SDSS3-BOSS data
• Parameter Inference for Weak Lensing using Gaussian Processes and MOPED
• The PAU Survey: Measurement of Narrow-band galaxy properties with Approximate Bayesian Computation
• Accurate X-ray Timing in the Presence of Systematic Biases With Simulation-Based Inference
• Reconstructing the Universe with Variational self-Boosted Sampling
• Perfectly parallel cosmological simulations using spatial comoving Lagrangian acceleration
• Exploring Early and Late Cosmology with Next Generation Surveys
• COmoving Computer Acceleration (COCA): $N$-body simulations in an emulated frame of reference
• Simulation-based inference of Bayesian hierarchical models while checking for model misspecification
• Diagnosing Systematic Effects Using the Inferred Initial Power Spectrum