Aemulus $\hat{I}½$: Precise Predictions for Matter and Biased Tracer Power Spectra in the Presence of Neutrinos

Abstract

We present the Aemulus $\hat{I}½$ simulations: a suite of 150 $(1.05h^{-1}\mathrm{Gpc}{)}^3$ $N$-body simulations with a mass resolution of $3.51\times 10^{10}\frac{\hat{I}{\text{c◯}}_{cb}}{0.3}{h}^{-1}{M}_{\odot }$ in a $w\hat{I}½$CDM cosmological parameter space. The simulations have been explicitly designed to span a broad range in $\ddot{I}{ƒ}_8$ to facilitate investigations of tension between large scale structure and cosmic microwave background cosmological probes. Neutrinos are treated as a second particle species to ensure accuracy to $0.5\mathrm{eV}$, the maximum neutrino mass that we have simulated. By employing Zel'dovich control variates, we increase the effective volume of our simulations by factors of $10-10^5$ depending on the statistic in question. As a first application of these simulations, we build new hybrid effective field theory and matter power spectrum surrogate models, demonstrating that they achieve $\le 1\%$ accuracy for $k\le 1h\mathrm{Mp}{\mathrm{c}}^{-1}$ and $0\le z\le 3$, and $\le 2\%$ accuracy for $k\le 4h\mathrm{Mp}{\mathrm{c}}^{-1}$ for the matter power spectrum. We publicly release the trained surrogate models, and estimates of the surrogate model errors in the hope that they will be broadly applicable to a range of cosmological analyses for many years to come.

37 pages, 15 figures, matching version accepted by JCAP