Improved recipes for peculiar velocity power spectra using Evolution Mapping

Abstract

Abstract We present new fitting functions for the velocity divergence auto- and cross-power spectra, Pθθ(k) and Pδθ(k), calibrated on gravity-only N-body simulations. By applying the Evolution Mapping framework, we revise existing prescriptions to introduce a physically motivated parametrisation in terms of the clustering amplitude σ12, the RMS density fluctuation smoothed at 12 Mpc. This approach improves robustness and extends the range of applicability beyond that of previous models. Our fits are calibrated using a suite of multi-resolution simulations, with numerical convergence carefully quantified and sampling artefacts mitigated through a conservative patching strategy. This yields converged measurements up to k ≃ 0.56 Mpc−1 and percent-level accuracy for both Pθθ(k) and Pδθ(k) over a wide range of σ12. Validation against independent simulations spanning a broad range of cosmological models confirms an accuracy of 1–2 per cent on scales where the measurements are robust, systematically outperforming existing prescriptions. We further assess the impact of deviations from the exact evolution mapping relation induced by differing growth histories. For most cosmologies of practical interest, we find that neglecting these effects introduces only subdominant errors. We show that expressing fitting functions in h-dependent units leads to spurious, unphysical dependencies on the Hubble parameter, even for models with identical linear clustering. This provides strong empirical support for parametrising non-linear evolution in terms of σ12 rather than σ8. Our fitting functions provide a robust description of velocity power spectra, with direct applications to redshift-space distortion modelling in galaxy redshift surveys.