Accelerating cosmological inference with Gaussian processes and neural networks -- an application to LSST Y1 weak lensing and galaxy clustering

Abstract

Studying the impact of systematic effects, optimizing survey strategies, assessing tensions between different probes and exploring synergies of different data sets require a large number of simulated likelihood analyses, each of which cost thousands of CPU hours. In this paper, we present a method to accelerate cosmological inference using emulators based on Gaussian process regression and neural networks. We iteratively acquire training samples in regions of high posterior probability which enables accurate emulation of data vectors even in high dimensional parameter spaces. We showcase the performance of our emulator with a simulated 3x2 point analysis of LSST-Y1 with realistic theoretical and systematics modelling. We show that our emulator leads to high-fidelity posterior contours, with an order of magnitude speed-up. Most importantly, the trained emulator can be re-used for extremely fast impact and optimization studies. We demonstrate this feature by studying baryonic physics effects in LSST-Y1 3x2 point analyses where each one of our MCMC runs takes approximately 5 minutes. This technique enables future cosmological analyses to map out the science return as a function of analysis choices and survey strategy.

arXiv comment13 pages, 8 figures, To be submitted to MNRAS

References in corpus (51)

• emcee: The MCMC Hammer
• Efficient Computation of CMB anisotropies in closed FRW models
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Cosmological parameters from CMB and other data: a Monte-Carlo approach
• The Cosmic Linear Anisotropy Solving System (CLASS) II: Approximation schemes
• First results from the IllustrisTNG simulations: matter and galaxy clustering
• The Simons Observatory: Science goals and forecasts
• Revising the Halofit Model for the Nonlinear Matter Power Spectrum
• Dancing in the dark: galactic properties trace spin swings along the cosmic web
• Cobaya: Code for Bayesian Analysis of hierarchical physical models
• Extragalactic Science, Cosmology and Galactic Archaeology with the Subaru Prime Focus Spectrograph (PFS)
• Efficient sampling of fast and slow cosmological parameters
• Cosmology from cosmic shear power spectra with Subaru Hyper Suprime-Cam first-year data
• CMB power spectrum parameter degeneracies in the era of precision cosmology
• The BAHAMAS project: Calibrated hydrodynamical simulations for large-scale structure cosmology
• Dark energy constraints from cosmic shear power spectra: impact of intrinsic alignments on photometric redshift requirements
• Towards a realistic population of simulated galaxy groups and clusters
• The MassiveBlack-II Simulation: The Evolution of Halos and Galaxies to z~0
• CosmoLike - Cosmological Likelihood Analyses for Photometric Galaxy Surveys
• The Coyote Universe II: Cosmological Models and Precision Emulation of the Nonlinear Matter Power Spectrum
• Concordance and Discordance in Cosmology
• Euclid preparation: II. The EuclidEmulator -- A tool to compute the cosmology dependence of the nonlinear matter power spectrum
• COSMOPOWER: emulating cosmological power spectra for accelerated Bayesian inference from next-generation surveys
• Deep optical galaxy counts with the Keck Telescope
• Cosmology Constraints from the Weak Lensing Peak Counts and the Power Spectrum in CFHTLenS
• The Aemulus Project II: Emulating the Halo Mass Function
• Cosmology with the Wide-Field Infrared Survey Telescope -- Multi-Probe Strategies
• Modeling baryonic physics in future weak lensing surveys
• Accounting for baryonic effects in cosmic shear tomography: Determining a minimal set of nuisance parameters using PCA
• Full-shape cosmology analysis of SDSS-III BOSS galaxy power spectrum using emulator-based halo model: a $5\%$ determination of ${\sigma }_8$
• The cosmology dependence of galaxy clustering and lensing from a hybrid $N$-body-perturbation theory model
• 2D-FFTLog: Efficient computation of real space covariance matrices for galaxy clustering and weak lensing
• Fast cosmological parameter estimation using neural networks
• Constraining neutrino mass with tomographic weak lensing one-point probability distribution function and power spectrum
• Bayesian optimisation for likelihood-free cosmological inference
• SPECULATOR: Emulating stellar population synthesis for fast and accurate galaxy spectra and photometry
• Bayesian emulator optimisation for cosmology: application to the Lyman-alpha forest
• Constraining neutrino mass with weak lensing Minkowski Functionals
• {\sc CosmoNet}: fast cosmological parameter estimation in non-flat models using neural networks
• Cosmology with the Wide-Field Infrared Survey Telescope -- Synergies with the Rubin Observatory Legacy Survey of Space and Time
• Dark Energy Survey Year 1 Results: Constraining Baryonic Physics in the Universe
• An emulator for the Lyman-$\alpha$ forest in beyond-$\Lambda$CDM cosmologies
• Fast Cosmological Parameter Estimation from Microwave Background Temperature and Polarization Power Spectra
• Cosmological parameter estimation via iterative emulation of likelihoods
• Kernel-Based Emulator for the 3D Matter Power Spectrum from CLASS
• Multi-Fidelity Emulation for the Matter Power Spectrum using Gaussian Processes
• Euclid-era cosmology for everyone: Neural net assisted MCMC sampling for the joint 3x2 likelihood
• $\text{matryoshka}$: Halo Model Emulator for the Galaxy Power Spectrum
• Parameter Inference for Weak Lensing using Gaussian Processes and MOPED
• Combined full shape analysis of BOSS galaxies and eBOSS quasars using an iterative emulator
• Interpreting Internal Consistency of DES Measurements

Cited by in corpus (4)

• Cosmological Constraints from Galaxy Clusters and Groups in the eROSITA Final Equatorial Depth Survey
• CONNECT: A neural network based framework for emulating cosmological observables and cosmological parameter inference
• Fast and robust Bayesian Inference using Gaussian Processes with GPry
• Growth and Geometry Split in Light of the DES-Y3 Survey