In this paper, we have studied a new fractional reaction-diffusion two-species system as an extension to the Rosenzweig-MacArthur reaction-diffusion di-trophic food chain system which models the spatial interactions between a prey and predator. To guarantee good working guidelines when numerically simulating the model, we first show that the system is locally asymptotically stable, as it provides good conditions and correct choice of ecological parameters to enhance a biologically meaningful result. We propose a fast and accurate method for numerical solutions of space fractional reaction-diffusion equations. The technique is based on Fourier spectral method in space and exponential integrator scheme in time. The complexity of fractional derivative index in fractional reaction diffusion model is numerically formulated and graphically displayed in one-, two- and three-dimensions. https://doi.org/10.28919/cmbn/3364