Emulating the non-linear matter power spectrum in mixed axion dark matter models

In order to constrain ultra light dark matter models with current and near future weak lensing surveys, we need the predictions for the non-linear dark matter power spectrum. This is commonly extracted from numerical simulations or from semi-analytical methods. For ultra light dark matter models, such numerical simulations are often very expensive due to the need of having a very high force resolution, often limiting the simulations to very small simulation boxes that do not contain very large scales. In this work, we take a different approach by relying on fast approximate N-body simulations. In these simulations, axion physics is only included in the initial conditions, allowing us to run a large number of simulations with varying axion and cosmological parameters. From our simulation suite, we used machine learning tools to create an emulator for the ratio of the dark matter power spectrum in mixed axion models – models where dark matter is a combination of CDM and axion – to that of ΛCDM. The resulting emulator only needs to be combined with existing emulators for ΛCDM to be able to be used in parameter constraints. We compared the emulator to semi-analytical methods, but a more thorough test to full simulations to verify the true accuracy of this approach is not possible at the present time and is left for future work.