Third order solutions of cosmological density perturbations in Horndeski's most general scalar-tensor theory with the Vainshtein mechanism

We study the third order solutions of the cosmological density perturbations in Horndeski's most general scalar-tensor theory under the condition that the Vainshtein mechanism is at work. In this work, we thoroughly investigate the independence property of the functions describing the nonlinear mode couplings, which is also useful for models within general relativity. Then, we find that the solutions of the density contrast and the velocity divergence up to third order are characterized by six parameters. Furthermore, the one-loop order power spectra obtained with third order solutions are described by four parameters. We exemplify the behavior of the one-loop order power spectra assuming the kinetic gravity braiding model, which demonstrates that the effect of modified gravity appears more significantly in the power spectrum of the velocity divergence than the density contrast.