HMCODE-2020: improved modelling of non-linear cosmological power spectra with baryonic feedback

We present an updated version of the hmcode augmented halo model that can be used to make accurate predictions of the non-linear matter power spectrum over a wide range of cosmologies. Major improvements include modelling of baryon-acoustic oscillation (BAO) damping in the power spectrum and an updated treatment of massive neutrinos. We fit our model to simulated power spectra and show that we can match the results with an root mean square (RMS) error of 2.5 per cent across a range of cosmologies, scales , and redshifts z < 2. The error rarely exceeds 5 per cent and never exceeds 16 per cent. The worst-case errors occur at z similar or equal to 2, or for cosmologies with unusual dark energy equations of state. This represents a significant improvement over previous versions of hmcode, and over other popular fitting functions, particularly for massive-neutrino cosmologies with high neutrino mass. We also present a simple halo model that can be used to model the impact of baryonic feedback on the power spectrum. This six-parameter physical model includes gas expulsion by active galactic nuclei (AGN) feedback and encapsulates star formation. By comparing this model to data from hydrodynamical simulations, we demonstrate that the power spectrum response to feedback is matched at the <1 per cent level for z < 1 and . We also present a single-parameter variant of this model, parametrized in terms of feedback strength, which is only slightly less accurate. We make code available for our non-linear and baryon models at https://github.com/alexander-mead/HMcode and it is also available within camb and soon within class.

Automated summary

The updated hmcode is able to predict the non-linear matter power spectrum accurately, showing significant improvements in matching simulation results with low error rates across various cosmologies and scales. A simple halo model is also introduced to model the impact of baryonic feedback on the power spectrum, providing a physical understanding of gas expulsion and star formation processes.