High precision accelerator for our hybrid model of the redshift space power spectrum

Abstract

Upcoming Large Scale Structure surveys aim to achieve an unprecedented level of precision in measuring galaxy clustering. However, accurately modeling these statistics may require theoretical templates that go beyond second-order perturbation theory, especially for achieving precision at smaller scales. In our previous work, we introduced a hybrid model for the redshift space power spectrum of galaxies. This model combines second-order templates with N-body simulations to capture the influence of scale-independent parameters on the galaxy power spectrum. However, the impact of scale-dependent parameters was addressed by precomputing a set of input statistics derived from computationally expensive N-body simulations. As a result, exploring the scale-dependent parameter space was not feasible in this approach. To address this challenge, we present an accelerated methodology that utilizes Gaussian processes, a machine learning technique, to emulate these input statistics. Our emulators exhibit remarkable accuracy, achieving reliable results with just 13 N-body simulations for training. We reproduce all necessary input statistics for a set of test simulations with an error of approximately 0.1 per cent in the parameter space within $5\sigma$ of the Planck predictions, specifically for scales around $k > 0.1$ $h$Mpc$^{-1}$. Following the training of our emulators, we can predict all inputs for our hybrid model in approximately 0.2,seconds at a specified redshift. Given that performing 13 N-body simulations is a manageable task, our present methodology enables us to construct efficient and highly accurate models of the galaxy power spectra within a manageable time frame.Comment: 22 pages, 9 figure