Improved recipes for peculiar velocity power spectra using Evolution Mapping

TL;DR

This paper introduces improved, physically motivated fitting functions for velocity divergence power spectra, calibrated on simulations, extending applicability and accuracy for redshift-space distortion analysis.

Contribution

The authors develop new calibration methods for velocity power spectra using Evolution Mapping, enhancing robustness and extending the model's validity beyond previous approaches.

Findings

Achieved percent-level accuracy up to k ≈ 0.56 Mpc^{-1}

Validated against diverse cosmological models with 1-2% precision

Demonstrated superiority over existing prescriptions in robustness and accuracy

Abstract

We present new fitting functions for the velocity divergence auto- and cross-power spectra, $P_{\theta\theta}(k)$ and $P_{\delta\theta}(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 $\sigma_{12}$, the RMS density fluctuation smoothed at $12\,\text{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\simeq0.56\,\mathrm{Mpc}^{-1}$ and percent-level accuracy for both $P_{\theta\theta}(k)$ and $P_{\delta\theta}(k)$ over a wide range of $\sigma_{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 $\sigma_{12}$ rather than $\sigma_{8}$. Our fitting functions provide a robust description of velocity power spectra, with direct applications to redshift-space distortion modelling in galaxy redshift surveys.