Emulating the Lyman-Alpha forest 1D power spectrum from cosmological simulations: New models and constraints from the eBOSS measurement

TL;DR

This paper introduces high-resolution cosmological simulations and a Gaussian process emulator for the Lyman-alpha forest 1D power spectrum, enabling precise cosmological parameter constraints from eBOSS data.

Contribution

The paper presents the Lyssa simulation suite, a new emulator framework, and demonstrates its application to derive cosmological constraints from Lyman-alpha forest data.

Findings

Emulator achieves sub-percent accuracy in power spectrum predictions.

Constraints on matter power spectrum parameters from eBOSS data.

Tension with Planck results when using certain priors on mean transmission.

Abstract

We present the Lyssa suite of high-resolution cosmological simulations of the Lyman-$\alpha$ forest designed for cosmological analyses. These 18 simulations have been run using the Nyx code with $4096^3$ hydrodynamical cells in a 120 Mpc ($\sim$ 81Mpc/h) comoving box and individually provide sub-percent level convergence of the Lyman-$\alpha$ forest 1d flux power spectrum. We build a Gaussian process emulator for the Lyssa simulations in the lym1d likelihood framework to interpolate the power spectrum at arbitrary parameter values. We validate this emulator based on leave-one-out tests and based on the parameter constraints for simulations outside of the training set. We also perform comparisons with a previous emulator, showing a percent level accuracy and a good recovery of the expected cosmological parameters. Using this emulator we derive constraints on the linear matter power spectrum amplitude and slope parameters $A_{\mathrm{Ly}\alpha}$ and $n_{\mathrm{Ly}\alpha}$. While the best-fit Planck $\Lambda$CDM model has $A_{\mathrm{Ly}\alpha}=8.79$ and $n_{\mathrm{Ly}\alpha}=-2.363$, from DR14 eBOSS data we find that $A_{\mathrm{Ly}\alpha}<7.6$ (95\% CI) and $n_{\mathrm{Ly}\alpha}=-2.369 \pm 0.008$. The low value of $A_{\mathrm{Ly}\alpha}$, in tension with Planck, is driven by the correlation of this parameter with the mean transmission of the Lyman-$\alpha$ forest. This tension disappears when imposing a well-motivated external prior on this mean transmission, in which case we find $A_{\mathrm{Ly}\alpha}=9.8\pm1.1$ in accordance with Planck.