{"id":"8b2894f0-0a28-423b-8d9d-5b2c6544ada5","arxiv_id":"2505.03999","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A renormalized perturbation theory is used to derive a scale-dependent effective gravitational constant for Horndeski gravity that recovers general relativity at small scales through Vainshtein screening.","lead":"This paper develops a method to compute how the effective strength of gravity changes from large to small scales in modified gravity theories with Vainshtein screening, where ordinary gravity is recovered in dense regions. The result is a scale-dependent gravitational constant that could sharpen predictions for how galaxies and large-scale structure cluster in such theories.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The master equation rests on an unvalidated closure: truncating the response expansion at Γ^(1) before substituting into the nonlinear Poisson equation (33) makes Eq. (85) an ansatz, not a derived resummation.","rationale":"The paper is a serious extension of DGP renormalized perturbation theory to a viable Horndeski subclass, and it contains real independent checks: the IR p² scaling is recovered in both the SPT and renormalized treatments, and Method 4's positivity-preserving iteration is a concrete algorithmic proposal. The reader's CONDITIONAL verdict is appropriate. My stress-test pass identifies one concern that is more specific and more load-bearing than the quasi-static/F-screening caveats: the master equation (85) is derived by substituting a first-order-only response ansatz into a nonlinear field equation, without establishing that the omitted connected response functions are subdominant in the equation that determines Γ^(1). Because the exact response expansion (39) includes a nonzero tree-level Γ^(2) from Eq. (56), the step from Eq. (33) to Eq. (81) is a closure assumption, not a derivation. The subsequent averaging over δm does not repair this; it converts the closure into a bispectrum-weighted self-consistency condition. The paper's own numerical experience (Methods 1–3 failing, only the positivity-enforcing Method 4 converging) further underscores that the equation's solution structure is delicate. The proposed concrete test—comparing the small-g expansion of the master equation with the exact one-loop propagator—would settle whether the closure is consistent with the perturbative response. Since the concern is addressable and does not by itself refute the framework, the verdict remains CONDITIONAL rather than moving to REJECT.","tokens_in":23953,"tokens_out":15545,"duration_ms":166706,"concrete_test":"Take the small-g (small µF−1) limit of Eq. (85) using the tree-level bispectrum (116) and compare the resulting Γ^(1)(p) with the exact one-loop propagator Γ^(1)_1-loop of Eq. (55) at the same order in G_NL_eff. If the two do not agree to within the estimated truncation error, or if the k→0 coefficients in Section 11 do not satisfy Eq. (165) numerically, the closure Φ≈Γ^(1)δm is not a valid reduction of Eq. (33) and the master equation is not the resummation of the response expansion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central formula (77) with M fixed by (85) inherits the closure ansatz Φ(p)=Γ^(1)(p)δm(p), i.e., the response functional (39) is truncated at first order before substitution into the nonlinear Poisson equation (33). This is not a controlled truncation: Eq. (56) gives a nonzero tree-level Γ^(2), and Section 5's estimate only shows Γ^(2)/Γ^(1) is small in a spherical/one-point sense, not that the contribution of Γ^(2) to the equation determining Γ^(1) is negligible. More importantly, Eq. (81) is presented as a field equation holding for each realization, but it asserts that δm(p) equals a deterministic convolution of δm with itself; for a genuinely stochastic δm this cannot hold pointwise. The statistical content is inserted only later by multiplying by δm(k) and averaging (Eqs. 82–85). That step replaces the field equation by a bispectrum-weighted closure, which is a new physical assumption rather than a consequence of Eq. (33). The paper acknowledges that F is not screened (Section 5, after Eq. 63) and defers the self-consistent iteration of P and B (Section 8), but the more basic issue is that no check is given that the M from (85) reproduces the perturbative Γ^(1) at one loop, Eq. (55). The IR comparison in Section 11 only verifies the p² scaling, not the coefficient; Eq. (165) is left as an unverified consistency condition. Thus the interpolation between µ at large scales and 1/F at small scales is plausible but currently rests on an unvalidated closure.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a renormalized perturbation theory for the quasi-static Horndeski subclass (1), with the goal of defining a scale- and time-dependent nonlinear effective gravitational constant G_NL_eff. After deriving the nonlinear Poisson equation (33), the authors introduce the response functional expansion (39), truncate it to the linear response Γ^(1), and write Γ^(1) in the form (76) in terms of an unknown function M(t,p). The main result is Eq. (77), with M determined by the integral equation (85) involving the matter power spectrum and bispectrum. Sections 5-8 develop a spherical approximation and a formal iteration scheme for the power spectrum and bispectrum. Section 10 presents numerical solutions for a model with F=1, cK=1, cB=0.5, using four methods of which Method 4 is selected. Section 11 compares the IR limit to one-loop SPT and finds the same p² scaling. The paper concludes that the framework captures Vainshtein screening and can be iterated in future work.","tokens_in":24338,"tokens_out":8684,"duration_ms":87555,"significance":"If the master equation (85) is a valid resummation, the paper provides a useful nonperturbative tool for estimating screening effects on the gravitational coupling in Horndeski models, with potential applications to power spectrum and bispectrum predictions. The derivation of the quasi-static field equations (13)-(33), the explicit nonlinear Poisson equation, and the spherical approximation are clear and carefully presented. Strengths include the absence of parameter fitting to the target quantity, the nonlinear integral-equation structure of Eq. (85), and the explicit comparison with one-loop SPT in the IR. However, the central equation rests on a closure that is not derived from a controlled truncation, and the numerical implementation inherits inconsistencies from the input spectra. The usefulness of the framework therefore depends on the additional validation requested below.","major_comments":[{"comment":"The master equation is an unvalidated closure, not a derived resummation. The substitution Φ(p)=Γ^(1)(p)δm(p) truncates the exact response expansion (42) at first order, while Eq. (56) gives a nonzero tree-level Γ^(2); the estimates in Eqs. (69)-(70) are based on a spherical one-point approximation and do not control the contribution of Γ^(2) to the equation for Γ^(1). Furthermore, Eq. (81) is presented as a pointwise field equation, but it equates the linear functional (M(p)-1)δm(p) to a quadratic convolution of δm with itself, which cannot hold for each realization of a stochastic density field; the statistical closure is inserted only when Eq. (82) is obtained by multiplying by δm(k) and averaging. The derivation should either start from functional differentiation of Eq. (33) and a systematic truncation of the propagator hierarchy, or the ansatz should be stated as such and tested, e.g., against the one-loop Γ^(1) of Eq. (55) and against N-body simulations.","section":"Section 7, Eqs. (78)-(85)"},{"comment":"The comparison with one-loop SPT establishes only the p² scaling of the IR correction, not the coefficient. Equation (165) is a nontrivial consistency condition that is left unchecked, and no demonstration is given that the solution M(p) of Eq. (85) reproduces the perturbative one-loop propagator Γ^(1)_1-loop in Eq. (55). Without this check, the statement that both approaches agree in the IR limit is premature. Please verify Eq. (165) numerically or analytically for the adopted parameter choices, or at least show the size of the residual.","section":"Section 11, Eqs. (161)-(165)"},{"comment":"The numerical results for G_NL_eff are presented without convergence criteria, residual measures, or error estimates. Methods 1-3 fail or oscillate, and Method 4 is selected because it converges in two iterations, but no quantitative stopping rule is provided. In addition, the power spectrum and bispectrum used in Eq. (85) come from halofit and SPT/BiHalofit and are not corrected by the Horndeski dynamics; the paper acknowledges this (§8) but does not estimate the resulting error in the converged M(p). Please report residuals, provide error bars or robustness tests, and quantify the sensitivity to the input P and B.","section":"Section 10, Figs. 6-7"}],"minor_comments":[{"comment":"The paper correctly acknowledges after Eq. (63) that F(ϕ) is not screened; for general F this means the small-scale limit of G_NL_eff is G/F, not G. Since the numerical model sets F=1, the abstract's claim of restoration of GR should be qualified to that choice.","section":"Section 3, Eq. (63)"},{"comment":"The notation 'p2γ' in Eq. (135) is ambiguous; it should read p1² γ to match Eq. (85), where p1 is the integration variable.","section":"Section 10.4, Eq. (135)"},{"comment":"In Eq. (82), the bispectrum argument is written as B(p1,p2,k), but the definition (83) uses B(p1,p2,p3) with p3=-k; use consistent notation and check the Dirac deltas in the averaging step.","section":"Section 7, Eq. (82)"},{"comment":"The replacement δm → ∆(k) in Eq. (67) mixes a dimensionless field with a dimensionless amplitude of the power spectrum; clarify that this is a heuristic identification and not an exact substitution.","section":"Section 5, Eq. (67)"},{"comment":"The 'elephant problem' is attributed to reference [1] in the introduction, but that reference is a general review; either provide the original source or cite a more specific reference.","section":"Introduction, paragraph 2"}],"recommendation":"major_revision","confidential_remarks":"The central issue is the validation of the closure leading to Eq. (85); the framework is potentially useful, but the derivation should be turned into a controlled approximation or explicitly tested against the one-loop limit and simulations before publication. The paper is within the scope of the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague], quick take on arXiv:2505.03999. The headline: this is a serious extension of Scoccimarro's DGP renormalized perturbation theory to a viable Horndeski subclass, and it gives a new master equation for the scale-dependent effective gravitational constant. The interesting physics is that G_NL_eff interpolates between the linear value μ at large scales and G/F at small scales, encoding Vainshtein screening. The numerics converge with a positivity-preserving iteration (Section 10.4) that is itself a useful contribution. The IR check against one-loop SPT in Section 11 is honest and passes at the level of the p² scaling.\n\nBut there is a real soft spot, and it is the one the stress-test note flags. The master equation is not derived by a controlled resummation. The authors truncate the response functional at Γ^(1) before substituting into the nonlinear Poisson equation (33). Equation (81) is presented as a pointwise identity for δm, which cannot hold for a stochastic field; the statistical content only enters later by multiplying by δm(k) and averaging. That step replaces a field equation with a bispectrum-weighted closure. It is a plausible physical assumption, and Scoccimarro's DGP approach has a similar flavor, but this paper does not show that the M(p) from (85) reproduces the one-loop Γ^(1) (Eq. 55) or that the neglected Γ^(2) is small in the equation that determines Γ^(1). The Section 5 estimate is spherical/one-point and doesn't control the full functional impact. So the central result, Eq. (77), is an ansatz whose regime of validity is not established.\n\nOther limitations are more minor and the paper is upfront about them: F is not screened (which the authors flag as a possible quasi-static artifact), the numerical pipeline uses unrenormalized P and B and defers the self-consistent iteration, and there are no error bars or released code. The coefficient-level IR comparison is left as the unverified consistency condition (165).\n\nWho this is for: people working on modified gravity predictions for the power spectrum and bispectrum. The formalism is genuinely new and the authors are honest about where it stops. As it stands, the numbers are illustrative, not ready for survey forecasts. With a serious revision that tests the closure (e.g., comparing M against the perturbative Γ^(1) at one loop) and either justifies or removes the pointwise-to-averaged step, this could be a solid method paper.\n\nMy recommendation: send it to peer review. A serious referee should spend time on Sections 7 and 11. The flaws are substantial but addressable, and the core idea deserves to be in the literature.\n\nBest.","headline":"A serious new formalism for a scale-dependent effective gravitational constant in screened Horndeski gravity, but the master equation is an ansatz and the numerics are not yet a full renormalization.","tokens_in":24823,"tokens_out":3814,"would_cite":true,"duration_ms":35669,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C25","83D05","85A40"],"pacs":["04.50.Kd","98.80.-k"],"model":"deepseek-v4-flash","headline":"A single function now sets the effective gravitational strength at every cosmic scale.","keywords":["Horndeski gravity","Vainshtein screening","renormalized perturbation theory","effective gravitational constant","quasi-static approximation","cosmological structure formation","modified gravity","matter power spectrum"],"falsifier":"Evaluate the same nonlinear Poisson equation in a direct numerical simulation of the $c_K=1$, $c_B=0.5$ model and measure $G_{\\mathrm{eff}}$ through $-p^2\\langle\\Phi\\delta_m\\rangle/P_{\\delta\\delta}(p)$; if the measured scale dependence does not match Eq. (77) with the converged $M(p)$ from Eq. (85), or if iterations of the full system move $M(p)$ outside $[0,1]$ at any scale where the quasi-static approximation is expected to hold, the central claim would be contradicted.","tokens_in":23740,"feed_emoji":"🌌","tokens_out":9117,"duration_ms":83270,"temperature":0.7,"pith_summary":"This paper develops a way to compute how the strength of gravity changes with scale in a broad class of modified gravity theories that hide their deviations from Einstein gravity at small scales through the Vainshtein mechanism. Extending an earlier renormalized perturbation approach, it defines a nonlinear effective gravitational constant $G_{\\mathrm{eff}}^{\\mathrm{NL}}/G$ as the response of the gravitational potential to a single Fourier mode of the matter density. The central result is that this coupling is controlled by one scale- and time-dependent function $M(t,p)$, which interpolates between $M=1$ at large scales (modified gravity with strength $\\mu$) and $M=0$ at small scales (strength $G/F$), with $M$ obeying a closed nonlinear integral equation. A sympathetic reader would care because this scale-dependent coupling is what must be put into growth and lensing equations to predict the matter power spectrum and bispectrum in these theories without assuming the screening is perfect.","feed_headline":"One function sets gravity's effective strength at every cosmic scale","feed_subtitle":"The new formula fixes how strong gravity is at every scale—key for clustering and lensing predictions.","key_machinery":"The machinery is a closure ansatz for the response function: the full linear response of the potential to the nonlinear density field is written as $\\Gamma^{(1)}(t,p) = -(3 a^2 H^2 \\Omega_m/2p^2)\\bigl(1 + (\\mu F-1)M(t,p)\\bigr)$, reducing all unknown nonlinearity to a single function $M$. Substituting this ansatz into the quasi-static, Fourier-space nonlinear Poisson equation, which keeps the Galileon-type kernel $\\gamma(p_1,p_2)=1-(p_1\\cdot p_2)^2/(p_1^2 p_2^2)$, and taking the bispectrum moment yields the master equation. The paper also introduces a positivity-preserving quadratic iteration for solving this integral equation, replacing a naive fixed-point iteration that drifts negative.","core_discovery":"On the paper's own terms, the central discovery is a resummation formula: after loop corrections are folded into the linear response, the renormalized effective gravitational constant is $G_{\\mathrm{eff}}^{\\mathrm{NL}}(t,p)/G = (1/F)\\bigl(1 + (\\mu F - 1) M(t,p)\\bigr)$, where $\\mu$ is the linear growth coupling, $F$ is the nonminimal coupling function, and $M(t,p)$ solves the master equation $$P(p)(1-M(p)) = \\frac{g}{(2\\pi)^2}\\int $p_1^{2}$\\, \\gamma(p_1,p-p_1)\\, B(p_1,p-p_1)\\, M(p_1)M(|p-p_1|)\\, dp_1\\, d\\mu_\\$\\theta$.$$ In the infrared $M\\to 1$ and one recovers the linear coupling $\\mu$; in the ultraviolet $M\\to 0$ and the coupling becomes $G/F$, showing that the derivative (Galileon-type) interactions screen the extra force but the conformal factor $F$ is not screened in this quasi-static treatment. The paper further exhibits a numerical iteration that converges in two steps and shows that the result reproduces the $p^2$ infrared behaviour of one-loop standard perturbation theory, while remaining finite where the one-loop expression fails.","pith_inferences":["Editorial extension: if $F\\neq 1$ is not screened, the small-scale limit is $G/F$, not Newton's $G$; in models where $F$ deviates appreciably at late times, this residual coupling would appear as a local fifth force unless a separate mechanism acts on $F$.","Editorial extension: the positivity-preserving quadratic iteration used here is not tied to the specific field theory; the same strategy could stabilize comparable nonlinear integral equations in other resummation schemes.","Editorial extension: because the master equation needs $P$ and $B$ as inputs, a sharp test of the framework would be a direct N-body measurement of $\\langle\\Phi\\delta_m\\rangle/P_{\\delta\\delta}$ in the $c_K=1$, $c_B=0.5$ model and a comparison with the iterated formula."],"forward_implications":["The effective gravitational coupling becomes a continuous function of scale and redshift, so Vainshtein screening is encoded as $M(t,p)\\to 0$ rather than requiring a sharp cutoff or a separate treatment of each object.","The matter power spectrum can be corrected by solving a single linear equation for a renormalized growth factor $D_{\\mathrm{renorm}}(t,p)$ with $G_{\\mathrm{eff}}^{\\mathrm{NL}}$ in place of Newton's constant.","Bispectrum corrections follow from the same coupling; the squeezed-limit second-order kernel is fixed by the Ward identity to $D_1 = [D_{\\mathrm{renorm}}]^2$, with no new soft poles, so non-Gaussianity changes only through the modified growth and finite multipole coefficients.","In the infrared the renormalized coupling recovers the one-loop SPT $p^2$ scaling, while at small scales it remains well behaved where one-loop SPT breaks down.","The converged result depends only weakly on the assumed input bispectrum, which suggests the method is stable enough for an iterative power-spectrum and bispectrum renormalization."],"supporting_citations":[{"why":"Supplies the renormalized perturbation response-function formalism that this paper extends from a braneworld model to Horndeski gravity.","marker":"[15]"},{"why":"Defines the $\\alpha_K, \\alpha_B$ parametrization used to specify the viable Horndeski subclass and evolve the background.","marker":"[30]"},{"why":"Provides the quasi-static scalar-field and metric perturbation equations, including the nonlinear Poisson equation, from which the master equation is built.","marker":"[31]"},{"why":"Gives the standard perturbation-theory kernels used for the power spectrum, bispectrum corrections, and the SPT comparison.","marker":"[33]"},{"why":"Establishes the extended Galilean-invariance and soft-limit properties of the density kernels used to fix the squeezed bispectrum.","marker":"[34]"},{"why":"Supplies the Ward identity that fixes the soft piece of the second-order kernel.","marker":"[35]"},{"why":"Provides the nonlinear matter power spectrum used as input to the numerical solution of the master equation.","marker":"[40]"},{"why":"Provides the fitting-form bispectrum used as one of the two inputs for the numerical iteration.","marker":"[41]"},{"why":"Gives the second-order SPT solutions for density and potential used to build the bispectrum and SPT comparison.","marker":"[45]"},{"why":"Gives the third-order SPT solutions whose cross-spectra define the one-loop effective coupling.","marker":"[46]"}],"fun_headline_variants":["One formula sets gravity's effective strength at all scales","New equation reveals gravity's screened strength in Horndeski","Master formula sets gravity's strength scale by scale","Resummed gravity coupling: one equation from IR to UV","How one resummation sets gravity's effective strength"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation's load-bearing premise is that the quasi-static approximation holds down to the small scales where Vainshtein screening operates, and that the only nonlinearities that matter there are the Galileon-type derivative terms retained in the truncated Poisson equation, so the unscreened coupling $F(\\phi)$ can be treated as a background function.","fun_headline_variants_meta":{"raw":{"variants":["One formula sets gravity's effective strength at all scales","New equation reveals gravity's screened strength in Horndeski","Master formula sets gravity's strength scale by scale","Resummed gravity coupling: one equation from IR to UV","How one resummation sets gravity's effective strength"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000515,"raw_usage":{"total_tokens":2516,"prompt_tokens":978,"completion_tokens":1538,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":594,"completion_tokens_details":{"reasoning_tokens":1459}},"tokens_in":594,"tokens_out":1538,"duration_ms":9477,"temperature":1.0,"reasoning_tokens":1459,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:40:19.113960+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the same nonlinear Poisson equation in a direct numerical simulation of the $c_K=1$, $c_B=0.5$ model and measure $G_{\\mathrm{eff}}$ through $-p^2\\langle\\Phi\\delta_m\\rangle/P_{\\delta\\delta}(p)$; if the measured scale dependence does not match Eq. (77) with the converged $M(p)$ from Eq. (85), or if iterations of the full system move $M(p)$ outside $[0,1]$ at any scale where the quasi-static approximation is expected to hold, the central claim would be contradicted.","supporting_citations":[{"cited_title":"Large-Scale Structure in Brane-Induced Gravity I. Perturbation Theory","cited_arxiv_id":"0906.4545","evidence_quote":"Supplies the renormalized perturbation response-function formalism that this paper extends from a braneworld model to Horndeski gravity."},{"cited_title":"Ward identities and consistency relations for the large scale structure with multiple species","cited_arxiv_id":"1310.7915","evidence_quote":"Supplies the Ward identity that fixes the soft piece of the second-order kernel."},{"cited_title":"Bispectrum of cosmological density perturbations in the most general second-order scalar-tensor theory","cited_arxiv_id":"1311.0281","evidence_quote":"Gives the second-order SPT solutions for density and potential used to build the bispectrum and SPT comparison."}],"review_version":1}