Why Newton's gravity is practically reliable in the large-scale   cosmological simulations
  Until now, it has been common to use Newton's gravity to study the non-linear clustering properties of the large-scale structures. Without confirmation from Einstein's theory, however, it has been unclear whether we can rely on the analysis, for example, near the horizon scale. In this work we will provide a confirmation of using Newton's gravity in cosmology based on relativistic analysis of weakly non-linear situations to the third order in perturbations. We will show that, except for the gravitational wave contribution, the relativistic zero-pressure fluid equations perturbed to the second order in a flat Friedmann background coincide exactly with the Newtonian results. We will also present the pure relativistic correction terms appearing in the third order. The third-order correction terms show that these are the linear-order curvature perturbation strength higher than the second-order relativistic/Newtonian terms. Thus, the pure general relativistic corrections in the third order are independent of the horizon scale and are small in the large-scale due to the low-level temperature anisotropy of the cosmic microwave background radiation. Since we include the cosmological constant, our results are relevant to currently favoured cosmology. As we prove that the Newtonian hydrodynamic equations are valid in all cosmological scales to the second order, and that the third-order correction terms are small, our result has a practically important implication that one can now use the large-scale Newtonian numerical simulation more reliably as the simulation scale approaches and even goes beyond the horizon. 
J. Hwang and H. Noh