Constraining the primordial power spectrum using a differentiable likelihood

Abstract

The simplest inflationary models predict the primordial power spectrum (PPS) of curvature perturbations to be nearly scale-invariant. However, various other models of inflation predict deviations from this behaviour, motivating a data-driven approach to reconstruct the PPS and constrain its shape. In this work, we present a novel method that employs a fully differentiable pipeline to reconstruct the PPS using Gaussian Processes and uses neural network emulators for fast and differentiable theoretical predictions. By leveraging gradient-based sampling techniques, such as Hamiltonian Monte Carlo, our approach efficiently samples the high-dimensional parameter space of cosmological parameters and the free-form PPS, enabling joint constraints on both. Applying this framework to Planck 2018 Cosmic Microwave Background (CMB) temperature anisotropy data we find our reconstructed PPS to be consistent with near scale-invariance on small scales, while exhibiting large uncertainties at large scales, driven mostly by cosmic variance. Our results show an overestimation of the PPS amplitude compared to $\Lambda$CDM predictions from the Planck 2018 analysis, which we attribute to our choice of a wider prior on the optical depth $\tau$ based on Planck 2015 measurements. Adopting a prior consistent with Planck 2018 measurements brings our results into full agreement with previous work. To ensure robustness of our results, we validate our differentiable pipeline against a non-differentiable framework, and also demonstrate that our results are insensitive to the choice of Gaussian process hyperparameters. These promising results and the flexibility of our pipeline make it ideally suited for application to additional data sets such as CMB polarisation as well as Large-Scale Structure probes, thus moving towards multi-probe primordial power spectrum reconstruction.

31 pages, 15 figures, 2 tables. Version accepted for publication in JCAP

References in corpus (63)

• emcee: The MCMC Hammer
• Efficient Computation of CMB anisotropies in closed FRW models
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Preliminary Maps and Basic Results
• Planck 2018 results. X. Constraints on inflation
• The No-U-Turn Sampler: Adaptively Setting Path Lengths in Hamiltonian Monte Carlo
• Planck 2018 results. I. Overview and the cosmological legacy of Planck
• The Simons Observatory: Science goals and forecasts
• Planck 2018 results. V. CMB power spectra and likelihoods
• Importance Nested Sampling and the MultiNest Algorithm
• Planck 2015 results. XI. CMB power spectra, likelihoods, and robustness of parameters
• Pyro: Deep Universal Probabilistic Programming
• Planck 2018 results. III. High Frequency Instrument data processing and frequency maps
• Inflationary perturbations from a potential with a step
• Gaussian Process Cosmography
• Reconstructing the primordial power spectrum
• Composable Effects for Flexible and Accelerated Probabilistic Programming in NumPyro
• The Coyote Universe Extended: Precision Emulation of the Matter Power Spectrum
• COSMOPOWER: emulating cosmological power spectra for accelerated Bayesian inference from next-generation surveys
• Primordial power spectrum from WMAP
• New Universal Local Feature in the Inflationary Perturbation Spectrum
• The Atacama Cosmology Telescope: a measurement of the primordial power spectrum
• Model-Independent Reconstruction of the Primordial Power Spectrum from WMAP Data
• Reconstructing the primordial power spectrum - a new algorithm
• Searching for Oscillations in the Primordial Power Spectrum: Constraints from Planck (Paper II)
• On Minimally-Parametric Primordial Power Spectrum Reconstruction and the Evidence for a Red Tilt
• Model selection applied to reconstruction of the Primordial Power Spectrum
• Implications of a Running Spectral Index for Slow Roll Inflation
• How to use GP: Effects of the mean function and hyperparameter selection on Gaussian Process regression
• Reconstruction of the Primordial Power Spectrum using Temperature and Polarisation Data from Multiple Experiments
• Is a step in the primordial spectral index favored by CMB data ?
• The Shape of the Primordial Power Spectrum: A Last Stand Before Planck
• Reconstructing the Primordial Spectrum from WMAP Data by the Cosmic Inversion Method
• {\sc CosmoNet}: fast cosmological parameter estimation in non-flat models using neural networks
• Smoothing spline primordial power spectrum reconstruction
• Searching for Oscillations in the Primordial Power Spectrum: Perturbative Approach (Paper I)
• CONNECT: A neural network based framework for emulating cosmological observables and cosmological parameter inference
• Bayesian inflationary reconstructions from Planck 2018 data
• Reconstructing the primordial power spectrum from the CMB
• Cosmic Inversion: Reconstructing primordial spectrum from CMB anisotropy
• Fast Cosmological Parameter Estimation from Microwave Background Temperature and Polarization Power Spectra
• Data compression in cosmology: A compressed likelihood for Planck data
• WMAP 7 Constraints on Oscillations in the Primordial Power Spectrum
• Accelerating Large-Scale-Structure data analyses by emulating Boltzmann solvers and Lagrangian Perturbation Theory
• Running Spectral Index as a Probe of Physics at High Scales
• Model Independent Primordial Power Spectrum from Maxima, Boomerang, and DASI Data
• Model-independent constraints on ${\Omega }_m$ and $H(z)$ from the link between geometry and growth
• CosmoPower-JAX: high-dimensional Bayesian inference with differentiable cosmological emulators
• Direct Wavelet Expansion of the Primordial Power Spectrum
• Capse.jl: efficient and auto-differentiable CMB power spectra emulation
• Euclid-era cosmology for everyone: Neural net assisted MCMC sampling for the joint 3x2 likelihood
• Parameter Inference for Weak Lensing using Gaussian Processes and MOPED
• candl: Cosmic Microwave Background Analysis with a Differentiable Likelihood
• Bayesian reconstruction of the inflaton's speed of sound using CMB data
• On the running of the spectral index to all orders: a new model dependent approach to constrain inflationary models
• The future of cosmological likelihood-based inference: accelerated high-dimensional parameter estimation and model comparison
• Red, Straight, no bends: primordial power spectrum reconstruction from CMB and large-scale structure
• The impact of the Universe's expansion rate on constraints on modified growth of structure
• PRISM: Recovery of the primordial spectrum from Planck data
• LimberJack.jl: auto-differentiable methods for angular power spectra analyses
• First step toward matter power spectrum reconstruction with Stage III weak gravitational lensing surveys
• Bayesian inference methodology for Primordial Power Spectrum reconstructions from Large Scale Structure
• KiDS-1000: Detection of deviations from a purely cold dark matter power spectrum with tomographic weak gravitational lensing
• A model-independent reconstruction of the matter power spectrum

Cited by in corpus (1)

• J-PAS: forecast on the primordial power spectrum reconstruction