We assess the performance of the Taruya, Nishimichi and Saito (TNS) model for the halo redshift space power spectrum, focusing on utilising mildly non-linear scales to constrain the growth rate of structure f . Using simulations with volume and number density typical of forthcoming Stage IV galaxy surveys, we determine ranges of validity for the model at redshifts z = 0.5 and z = 1. We proceed to perform a Bayesian MCMC analysis utilising the monopole, quadrupole, and hexadecapole spectra, followed by an exploratory Fisher matrix analysis. As previously noted in other forecasts as well as in real data analyses, we find that including the hexadecapole can significantly improve the constraints. However, a restricted range of scales is required for the hexadecapole in order for the growth parameter estimation to remain unbiased, limiting the improvement. We consistently quantify these effects by employing the multipole expansion formalism in both our Fisher and MCMC forecasts.A way to attack this challenge is by using perturbation theory based models (Bernardeau et al. 2002;Kaiser 1987;Scoccimarro 2004) that can be easily extended to include non-standard theories of gravity and dark energy (see Koyama 2016, 2017; Bose et al. 2018a,b, for example). These can be combined with phenomenological ingredients to model non-linear physics (Taruya et al. 2010;Senatore and Zaldarriaga 2014;de la Bella et al. 2017). Additionally, a model for the galaxy bias is required to relate the dark matter and galaxy distributions .An important factor, which is the main focus of this paper, is our ability to use these prescriptions to model non-linear structure formation in an unbiased manner, and how this affects the growth of structure parameter estimation with Stage IV surveys. dida.markovic@port.ac.uk benjamin.bose@unige.ch a.pourtsidou@qmul.ac.uk 1 www.euclid-ec.org 2 https://wfirst.gsfc.nasa.gov/ 3 www.desi.lbl.gov