Nonlinear evolution of the matter power spectrum in modified theories of gravity

We present a formalism to calculate the nonlinear matter power spectrum in modified gravity models that explain the late-time acceleration of the Universe without dark energy. Any successful modified gravity models should contain a mechanism to recover general relativity (GR) on small scales in order to avoid the stringent constrains on deviations from GR at solar system scales. Based on our formalism, the quasi-nonlinear power spectrum in the Dvali-Gabadadze-Porratti braneworld models and f(R) gravity models are derived by taking into account the mechanism to recover GR properly. We also extrapolate our predictions to fully nonlinear scales using the parametrized post-Friedmann framework. In the Dvali-Gabadadze-Porratti and f(R) gravity models, the predicted nonlinear power spectrum is shown to reproduce N-body results. We find that the mechanism to recover GR suppresses the difference between the modified gravity models and dark energy models with the same expansion history, but the difference remains large at the weakly nonlinear regime in these models. Our formalism is applicable to a wide variety of modified gravity models and it is ready to use once consistent models for modified gravity are developed.