{"id":"adc811a3-7f49-4d58-84bc-df9ebe5d8644","arxiv_id":"1908.02273","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Random nonlinear elliptic equations homogenize at the optimal rate: order ε in d≥3 and order ε|log ε|^{1/2} in d=2, matching the sharp linear-elliptic rates.","lead":"This paper proves that solutions to nonlinear elliptic equations with tiny random, quickly decorrelating coefficients differ from their averaged effective solutions by only one microscopic length scale in three or more dimensions. It is the first sharp error bound for this class, and it also bounds how well small numerical samples estimate the effective material law.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified; the optimal-rate proof is internally consistent under (R) and the spectral-gap assumptions, with the conditionality already correctly tied to the deferred proofs in [21].","rationale":"The reader's CONDITIONAL verdict is appropriate. I found no mathematical contradiction, no fitted-parameter circularity, and the optimality argument via the linear subclass is sound. The stochastic estimates are all routed through the spectral gap inequality as the paper states; the dependence of the optimal rate on (R) is explicit. The only issue serious enough to keep the paper from ACCEPT is the deferred proof of localization and the two-scale expansion, which are load-bearing. Agreement is partial because the reader's weakest assumption was the spectral gap, whereas I weight the deferred proof components more heavily; both are real but the latter is the one that would block verification of the main theorems today.","tokens_in":85492,"tokens_out":32325,"duration_ms":385774,"concrete_test":"Independently prove Lemma 18 and the residual representation (14)/Proposition 36 from assumptions (A1)-(A3) and (P1)-(P2) without citing [21]; if either step fails or requires additional hypotheses, Theorems 2 and 4 must be amended or their assumptions strengthened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"I could not identify a new load-bearing flaw. The central claims follow from the stated assumptions (P1)-(P2), (R), and the corrector/multiscale-expansion lemmas, and the rate-matching argument from the linear subclass is legitimate. The weakest point remains the one the reader flagged: Lemma 18 (existence of localized correctors) and the formal two-scale expansion derivation are only sketched or deferred to the forthcoming paper [21], and the main theorems depend on them. This does not indicate an internal inconsistency, but it keeps the proof conditional until those components are available in published form.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a quantitative stochastic homogenization theory for uniformly elliptic monotone equations and systems with stationary random coefficients that decorrelate on scale ε. Under a spectral-gap inequality and an additional small-scale C^{1,α} regularity condition, it proves homogenization errors of order ε for d≥3 and of order ε|log ε|^{1/2} for d=2 (with a massive term), with random constants having bounded stretched exponential moments. Without the small-scale regularity condition, half rates are claimed, including on bounded domains. The paper also proves RVE approximation error rates of order (L/ε)^{-d/2} for the effective operator, with refined higher-order systematic error estimates. The proof strategy is based on localized correctors and linearized correctors, spectral-gap fluctuation estimates, a two-scale expansion with piecewise affine macroscopic slope, and flux-corrector residual estimates.","tokens_in":85401,"tokens_out":13073,"duration_ms":152073,"significance":"If correct, these results constitute the first optimal-order homogenization rates for a nonlinear stochastic homogenization problem, and the rate-matching to the linear elliptic subfamily makes the optimality argument transparent. The paper has substantial technical strengths: explicit stretched-exponential stochastic moment bounds, a modular proof structure that separates corrector estimates from the two-scale expansion and RVE analysis, and careful treatment of low dimensions and massive terms. The main caveat is that several load-bearing proof ingredients are deferred to a companion paper [21], which makes the manuscript, as submitted, not fully self-contained.","major_comments":[{"comment":"Lemma 18 is load-bearing: it provides existence, uniqueness, and continuity of the localized correctors and flux correctors on which every subsequent corrector estimate, including Propositions 19 and 20 and hence Theorems 2, 4, 7, and 14, depends. The proof is not contained in the manuscript; the text states 'For a detailed proof see [21]'. This makes the central theorems conditional on an unpublished companion paper. A journal version should either include this proof in full or state the main results as conditional on [21] with the exact dependence made explicit.","section":"3.1, Lemma 18"},{"comment":"The passage from the localized correctors with finite T to the actual correctors is handled by Lemma 33, whose proof is only sketched and points to 'the beginning of the proof of Lemma 40'; the argument there is described as a quantitative proof for φ that 'extends to σ', but it is not written out. Since Corollary 21 and the final error estimates rely on this T→∞ limit, the manuscript should provide a complete convergence argument for both the corrector and the flux corrector, or give a precise published reference.","section":"3.1-3.3, two-scale expansion"},{"comment":"The estimate for the vector potential θ_T is obtained by omitting the massive regularization of the equation Δθ_T = ∇φ_T; the authors write 'we omit this additional technicality'. This estimate feeds into Proposition 19 and therefore into the d≥3 rates, so the omitted argument is not purely cosmetic. The proof should be included or the status of the omitted regularization should be made precise.","section":"5.1, proof of Lemma 24, Part c"}],"minor_comments":[{"comment":"The displayed moment bound 'EL[exp(C1/C/C)]≤2' appears garbled; it should presumably read E_L[exp((C/C)^{1/C})]≤2 or an equivalent form. Please correct the notation throughout the theorem.","section":"Theorem 14(a)"},{"comment":"In the displayed estimate of Lemma 50, the sum on the right-hand side is written with ∇φ_T but the lemma concerns an arbitrary function v; it should be ∇v. As written, the statement is internally inconsistent, although the proof makes the intended meaning clear.","section":"Lemma 50"},{"comment":"There are several places where u is used instead of u_ε in expressions such as 'The error u−uhom' in Section 1.1. Please standardize the notation.","section":"1.1 and 2.1"},{"comment":"The symbol C is used both for deterministic constants and for random constants; in formulas such as E[exp(C^{ν}/C)]≤2 this creates ambiguity. Using a distinct notation such as C(ω) for the random constant would improve readability.","section":"Throughout"},{"comment":"The companion paper [21] is cited as forthcoming. If it is already available, the reference should be updated with preprint or publication data; otherwise the dependence of the main theorems on unpublished work should be addressed as in the major comments.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is mathematically serious and the central results appear likely correct, but it is not self-contained: Lemma 18, the T→∞ limit for flux correctors, and a regularization argument for the vector potential are deferred or omitted, and the main theorems depend on them. My major-revision recommendation is driven by completeness rather than by any detected internal contradiction. The editor may wish to verify that the companion paper [21] is publicly available or will be submitted with the revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Fischer and Neukamm claim the first optimal-order convergence rates for stochastic homogenization of nonlinear monotone elliptic equations and systems: ε in d≥3, ε√|log ε| in d=2, and a (L/ε)^{-d/2} bound for representative volume approximations. I read enough of the manuscript to believe the claim is real. The rates are genuinely new compared with the ε^δ rates of Armstrong–Smart, Armstrong–Mourrat, and Armstrong–Ferguson–Kuusi, and matching the linear-elliptic rates is a valid optimality argument because linear problems are a subclass.\n\nWhat the paper does well: the proof is built from explicit objects — localized correctors, flux correctors, linearized correctors, a piecewise-affine two-scale expansion, spectral-gap sensitivity estimates, and moment bounds on minimal radii. I found no circular step and no fitted parameters; the rates follow from the stated assumptions (P1)–(P2) and (R). The paper is also transparent about limitations: without (R) it only gets half rates, and it explains why the spectral gap approach does not cover strongly correlated or discrete-valued fields.\n\nThe soft spots are exactly what the reader flagged. Lemma 18, which establishes existence of localized correctors, and the formal derivation of the two-scale expansion residual are both deferred to the authors' forthcoming paper [21]. The main theorems (2, 4, 7, 14) depend on those two pieces, so the proof is conditional until they are available. That is a real gap, not a cosmetic one. The RVE result also requires an L-periodic approximation PL whose existence is verified only in cases; the paper states that the scope is restricted accordingly. These are specific, checkable items, so a referee can pin them down. A minor point: the stochastic integrability is only stretched exponential rather than Gaussian, but the paper acknowledges this.\n\nWho this is for: anyone working on quantitative homogenization or random media. It deserves peer review. I would make final acceptance contingent on the companion paper's proofs being available or replaced by in-paper arguments, but the present work is substantial enough to warrant serious referee time.","headline":"A strong and largely explicit proof of the first optimal-order rates in nonlinear stochastic homogenization, but the main theorems lean on two key lemmas deferred to a forthcoming paper, so the proof is conditional until those appear.","tokens_in":86136,"tokens_out":2558,"would_cite":true,"duration_ms":30300,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B27","35J60","60H25","35R60"],"pacs":[],"model":"deepseek-v4-flash","headline":"Optimal homogenization rates for nonlinear random elliptic PDEs are O(ε) in d≥3, matching linear theory.","keywords":["stochastic homogenization","monotone nonlinear elliptic PDE","optimal convergence rate","spectral gap inequality","corrector estimates","representative volume","random coefficient field","two-scale expansion"],"falsifier":"Construct a stationary random field ω_ε for which the spectral gap inequality (P2) holds but the corrector φ_ξ fails to satisfy the bound E[|∇φ_ξ|^2]^{1/2} ≲ |ξ| ε in d=3; then the paper's central homogenization-rate theorem would be false. Concretely, one could numerically compute the homogenization error for a scalar monotone equation with A(ω, ξ) = (1 + η(ω_ε(x))) ξ for a smooth Gaussian field with covariance decaying like (1 + |x−y|/ε)^{-(d+κ)}, κ>0, and check whether the error is indeed O(ε) in $L^{{2d/(d−2)}}$. A slower rate, such as ε|log ε|^{1/2}, would contradict the theorem.","tokens_in":85103,"feed_emoji":"📐","tokens_out":2700,"duration_ms":31305,"temperature":0.7,"pith_summary":"This paper establishes that solutions of random nonlinear uniformly elliptic equations with monotone nonlinearity converge to the homogenized solution at the optimal rate: order ε in d≥3, and order ε|log ε|^{1/2} in d=2 with a massive term. These are the same rates known for linear elliptic equations, which are a special case, so the rates cannot be improved in general. Previous nonlinear stochastic homogenization results only gave a small algebraic rate ε^δ. The paper also proves optimal-order error estimates for approximating the homogenized operator by representative volumes, with rate (L/ε)^{-d/2}. If correct, this is the first optimal-order quantitative homogenization result for any nonlinear stochastic homogenization problem.","feed_headline":"Optimal error rates for nonlinear random materials","feed_subtitle":"Homogenization error scales as ε in d≥3, matching the linear case; RVE approximations converge at (L/ε)^{-d/2}.","key_machinery":"The central objects are the homogenization corrector φ_ξ (a field correcting the macroscopic gradient ξ to the microscopic gradient ξ + ∇φ_ξ) and the flux corrector σ_ξ (a d−1-form solving ∇·σ_ξ = A(ω_ε, ξ + ∇φ_ξ) − A_hom(ξ)). Alongside these stand their localized versions φ_ξ^T, σ_ξ^T with a massive term 1/T that provides exponential localization, and the correctors for the linearized PDE around ξ + ∇φ_ξ, denoted φ_{ξ,Ξ}^T and σ_{ξ,Ξ}^T. The key quantization mechanism is a spectral gap inequality on the random field ω_ε: for every random variable F(ω_ε), the variance is bounded by ε^d times the $L^{2}$ norm of a spatial average of the sensitivity ∂F/∂ω_ε. Feeding corrector sensitivities into this inequality yields optimal fluctuation bounds, and a minimal-radius construction with a hole-filling estimate turns these into L^p corrector bounds. The optimal rates emerge from combining corrector bounds with a piecewise-affine two-scale expansion of u_hom.","core_discovery":"The central claim is that for a random monotone operator A(ω_ε(x), ξ) with stationary law (P1) and spectral-gap decorrelation on scale ε (P2), and under a small-scale $C^{{1,α}}$ regularity condition (R), the solution u_ε of the random PDE and the solution u_hom of the homogenized PDE satisfy ||u_ε − u_hom||_{$L^{{2d/(d−2)}}$(R^d)} ≤ C Ĉ(∇u_hom) ε for d≥3, and ||u_ε − u_hom||_{$L^{2}$(R^d)} ≤ C Ĉ(∇u_hom) ε |log ε|^{1/2} for d=2, with C a random constant having bounded stretched exponential moments. The same framework yields optimal RVE error bounds. The argument works for scalar equations, two-dimensional systems, and systems with Uhlenbeck structure; without (R) only half the rate is obtained. The proof identifies the correctors of the nonlinear problem with those of the linearized problem and exploits the spectral gap inequality to control their stochastic fluctuations.","pith_inferences":["If the spectral gap assumption is relaxed to slower-than-integrable correlations, the paper's strategy predicts that the ε and ε|log ε|^{1/2} rates fail; the rates likely degrade to those governed by the tail behavior of the corrector, analogous to the linear case.","The linearized-corrector identification suggests that derivative bounds for the effective operator A_hom(ξ) with respect to ξ might be obtainable by the same machinery, opening a route toward quantitative homogenization of nonlinear elasticity in the small-deformation regime.","A natural testable extension is to scalar monotone equations with p-growth (p≠2): the spectral-gap and two-scale-expansion framework would predict some rate ε^α(p), but the present paper does not cover it, and whether α(p) remains 1 for all p is an open question.","The RVE error estimates rely on a periodization of the field; for random fields where a matching periodization is hard to construct, the same fluctuation estimates suggest a coupling error that could be controlled by directly comparing local laws, possibly extending the result beyond periodic RVEs."],"forward_implications":["Numerical simulations of nonlinear random materials can be justified at the same resolution as linear ones: the RVE approximation of the effective material law converges at the central-limit rate (L/ε)^{-d/2} in d≤4.","The homogenization error in d≥3 is linear in the microstructural scale, so doubling the scale resolution halves the error—a sharp quantitative target for nonlinear stochastic homogenization.","In d=2 a logarithmic correction is unavoidable, matching the critical dimensional behavior of linear elliptic PDEs.","For systems with Uhlenbeck structure, or for scalar equations, the same optimal rates hold without assuming extra structure beyond the small-scale regularity condition.","The two-scale expansion with correctors also yields an H^1 error estimate for ∇u_ε − ∇û_ε of the same optimal order, making the microscopic gradient approximation quantitative."],"supporting_citations":[{"why":"Provide the spectral-gap-based optimal-rate framework for linear stochastic homogenization that the present proof adapts.","marker":"[32, 33]"},{"why":"Supply the localized corrector and flux-corrector estimates for linear elliptic equations that the nonlinear corrector estimates extend.","marker":"[31, 30]"},{"why":"Establishes the previous algebraic rate ε^δ for monotone nonlinear stochastic homogenization and introduces the linearized corrector concept that is central here.","marker":"[2]"},{"why":"Give the optimal-rate linear results (with Gaussian-type stochastic moments) whose rates the paper matches and uses as the optimality benchmark.","marker":"[4, 34]"},{"why":"Provide the qualitative homogenization theory for monotone systems that the paper makes quantitative.","marker":"[17, 18]"},{"why":"Supplies the first optimal periodic-RVE result in a linear discrete setting, which the RVE part of the paper generalizes to nonlinear monotone operators.","marker":"[29]"}],"fun_headline_variants":["Optimal ε homogenization for nonlinear random PDEs","Matching linear-case rates in nonlinear stochastic homogenization","Sharp convergence rates for random monotone elliptic operators","From ε^δ to ε: optimal rates in random nonlinear homogenization","Linear-rate convergence for nonlinear random materials"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the spectral gap inequality (P2): every random variable built from the coefficient field must have variance controlled by the ε-scaled integral of its local sensitivity, which encodes that correlations decay fast enough on scales larger than ε.","fun_headline_variants_meta":{"raw":{"variants":["Optimal ε homogenization for nonlinear random PDEs","Matching linear-case rates in nonlinear stochastic homogenization","Sharp convergence rates for random monotone elliptic operators","From ε^δ to ε: optimal rates in random nonlinear homogenization","Linear-rate convergence for nonlinear random materials"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000284,"raw_usage":{"total_tokens":1682,"prompt_tokens":958,"completion_tokens":724,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":574,"completion_tokens_details":{"reasoning_tokens":648}},"tokens_in":574,"tokens_out":724,"duration_ms":54980,"temperature":1.0,"reasoning_tokens":648,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:49:31.878888+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Construct a stationary random field ω_ε for which the spectral gap inequality (P2) holds but the corrector φ_ξ fails to satisfy the bound E[|∇φ_ξ|^2]^{1/2} ≲ |ξ| ε in d=3; then the paper's central homogenization-rate theorem would be false. Concretely, one could numerically compute the homogenization error for a scalar monotone equation with A(ω, ξ) = (1 + η(ω_ε(x))) ξ for a smooth Gaussian field with covariance decaying like (1 + |x−y|/ε)^{-(d+κ)}, κ>0, and check whether the error is indeed O(ε) in $L^{{2d/(d−2)}}$. A slower rate, such as ε|log ε|^{1/2}, would contradict the theorem.","supporting_citations":[{"cited_title":"Homogenization, linearization and large-scale regularity for nonlinear elliptic equations","cited_arxiv_id":"1805.00467","evidence_quote":"Establishes the previous algebraic rate ε^δ for monotone nonlinear stochastic homogenization and introduces the linearized corrector concept that is central here."},{"cited_title":"Gloria, S","cited_arxiv_id":null,"evidence_quote":"Supplies the first optimal periodic-RVE result in a linear discrete setting, which the RVE part of the paper generalizes to nonlinear monotone operators."}],"review_version":1}