We study the impact of different bias and redshift-space models on the halo power spectrum, quantifying their effect by comparing the fit to a subset of realizations taken from the WizCOLA suite. These provide simulated power spectrum measurements between $k_{\rm min}$ = 0.03 h/Mpc and $k_{\rm max}$ = 0.29 h/Mpc, constructed using the comoving Lagrangian acceleration method. For the bias prescription we include (i) simple linear bias; (ii) the McDonald & Roy model and (iii) its coevolution variant introduced by Saito et al.; and (iv) a very general model including all terms up to one-loop and corrections from advection. For the redshift-space modelling we include the Kaiser formula with exponential damping and the power spectrum provided by (i) tree-level perturbation theory and (ii) the Halofit prescription; (iii) one-loop perturbation theory, also with exponential damping; and (iv) an effective field theory description, also at one-loop, with damping represented by the EFT subtractions. We quantify the improvement from each layer of modelling by measuring the typical improvement in chi-square when fitting to a member of the simulation suite. We attempt to detect overfitting by testing for compatibility between the best-fit power spectrum per realization and the best-fit over the entire WizCOLA suite. For both bias and the redshift-space map we find that increasingly permissive models yield improvements in chi-square but with diminishing returns. The most permissive models show modest evidence for overfitting. Accounting for model complexity using the Bayesian Information Criterion, we argue that standard perturbation theory up to one-loop, or a related model such as that of Taruya, Nishimichi & Saito, coupled to the coevolution bias model, is likely to provide a good compromise for near-future galaxy surveys operating with comparable $k_{\rm max}$.

The use of Eulerian 'standard perturbation theory' to describe mass assembly in the early universe has traditionally been limited to modes with k $\lesssim$ 0.1 h/Mpc at z = 0. At larger k the SPT power spectrum deviates from measurements made using N-body simulations. Recently, there has been progress in extending the reach of perturbation theory to larger k using ideas borrowed from effective field theory. We revisit the computation of the redshift-space matter power spectrum within this framework, including for the first time the full one-loop time dependence. We use a resummation scheme proposed by Vlah et al. to account for damping of baryonic acoustic oscillations due to large-scale random motions and show that this has a significant effect on the multipole power spectra. We renormalize by comparison to a suite of custom N-body simulations matching the MultiDark MDR1 cosmology. At z = 0 and for scales k $\lesssim$ 0.4 h/Mpc we find that the EFT furnishes a description of the real-space power spectrum up to $\sim$2%, for the $\ell$ = 0 mode up to $\sim$5%, and for the $\ell$ = 2, 4 modes up to $\sim$25%. We argue that, in the MDR1 cosmology, positivity of the $\ell$ = 0 mode gives a firm upper limit of k $\approx$ 0.75 h/Mpc for the validity of the one-loop EFT prediction in redshift space using only the lowest-order counterterm. We show that replacing the one-loop growth factors by their Einstein-de Sitter counterparts is a good approximation for the $\ell$ = 0 mode, but can induce deviations as large as 2% for the $\ell$ = 2, 4 modes. An accompanying software bundle, distributed under open source licenses, includes Mathematica notebooks describing the calculation, together with parallel pipelines capable of computing both the necessary one-loop SPT integrals and the effective field theory counterterms.