Abstract

This paper develops a renormalized perturbation theory framework for nonlinear structure formation in a broad class of modified gravity models that exhibit Vainshtein screening, with a focus on a viable subclass of Horndeski theories.We extend earlier perturbative methods, originally applied to the Dvali–Gabadadze– Porrati model, to construct a self-consistent treatment that captures both the linear modifications to gravity at large scales and the nonlinear screening effects at small scales. In the framework, the response of the gravitational potential to matter density fluctuations is characterized by renormalized propagators, leading to the definition of a nonlinear (or renormalized) effective gravitational constant. This paper details several numerical strategies to compute this renormalized gravitational constant. Numerical examples illustrate how the effective gravitational constant evolves with scale and redshift. These results are key to accurately predicting cosmological observables such as the matter power spectrum and bispectrum in modified gravity scenarios. © (2025), (American Physical Society). All rights reserved.