We study the third order solutions of the cosmological density perturbations in the 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 the general relativity. Then, we find that the solutions of the density contrast and the velocity divergence up to the third order ones are characterized by 6 parameters. Furthermore, the 1-loop order power spectra obtained with the third order solutions are described by 4 parameters. We exemplify the behavior of the 1-loop order power spectra assuming the kinetic gravity braiding model, which demonstrates that the effect of the modified gravity appears more significantly in the power spectrum of the velocity divergence than the density contrast.

We study the bispectrum of the matter density perturbations induced by the large scale structure formation in the most general second-order scalar-tensor theory that may possess the Vainshtein mechanism as a screening mechanism. On the basis of the standard perturbation theory, we derive the bispectrum being expressed by a kernel of the second order of the density perturbations. We find that the kernel at the leading order is characterized by one parameter, which is determined by the solutions of the linear density perturbations, the Hubble parameter and the other function specifying nonlinear interactions. This does not allow for varied behavior in the bispectrum of the matter density perturbations in the most general second-order scalar-tensor theory equipped with the Vainshtein mechanism. We exemplify the typical behavior of the bispectrum in a kinetic gravity braiding model.