A Multi-Fidelity Emulator for the Matter Power Spectrum using Gaussian Processes

We present methods for emulating the matter power spectrum which effectively combine information from cosmological $N$-body simulations at different resolutions. An emulator allows estimation of simulation output by interpolating across the parameter space of a handful of simulations. We present the first implementation of multi-fidelity emulation in cosmology, where many low-resolution simulations are combined with a few high-resolution simulations to achieve an increased emulation accuracy. The power spectrums dependence on cosmology is learned from the low-resolution simulations, which are in turn calibrated using high-resolution simulations. We show that our multi-fidelity emulator can achieve percent-level accuracy on average with only $3$ high-fidelity simulations and outperforms a single-fidelity emulator that uses $11$ simulations. With a fixed number of high-fidelity training simulations, we show that our multi-fidelity emulator is $simeq 100$ times better than a single-fidelity emulator at $k leq 2 ,htextrm{Mpc}{^{-1}}$, and $simeq 20$ times better at $3 leq k < 6.4 ,htextrm{Mpc}{^{-1}}$. Multi-fidelity emulation is fast to train, using only a simple modification to standard Gaussian processes. Our proposed emulator shows a new way to predict non-linear scales by fusing simulations from different fidelities.