UC Berkeley

We present the Aemulus ν simulations: a suite of 150 (1.05 h-1 Gpc)3 N-body simulations with a mass resolution of 3.51 × 1010 Ωcb /0.3 h-1 M ⊙ in a wνCDM cosmological parameter space. The simulations have been explicitly designed to span a broad range in σ 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 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-105 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 ≤ 1% accuracy for k ≤ 1 hMpc-1 and 0 ≤ z ≤ 3, and ≤ 2% accuracy for k ≤ 4 hMpc-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.