Pith. sign in

REVIEW 2 major objections 3 minor 23 references

Time-dependent averages of a critical long-range stochastic heat equation

T0 review · 2 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read For a critical long-range stochastic heat equation in d≥3, ball averages of the solution are Gaussian when t≪R², a non-Gaussian law when t=cR², and extinct when t≫R².

desk verdict A clean, genuinely new phase diagram for spatial averages of a critical SHE; the non-Gaussian limit is real, and the main proof concern raised by the stress test turns out to be a misreading. read the letter →

arxiv 2411.09058 v2 pith:YJX6CRSC submitted 2024-11-13 math.PR

Sefika Kuzgun , Ran Tao This is my paper · ORCID
classification math.PR MSC 60H1560F05
keywords stochasticheatequationRieszkernelcriticalSPDEspatialaveragescentrallimittheoremnon-GaussianscalinginvarianceWienerchaos
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper establishes a three-regime limit law for the spatial average of the solution to the critical stochastic heat equation $\partial_t u=\frac12\Delta u+\kappa u\dot F$ in $d\ge 3$, where $\dot F$ is white in time and has spatial covariance $\|x-y\|^{-2}$. The solution is a singular random measure, so the ball average $u_t(B_R)$ is the natural observable. The main theorem shows that for $t_R=o(R^2)$, the normalized average converges to a Gaussian with variance $\sigma^2=\kappa^2\int_{B_1}\int_{B_1}\|x-y\|^{-2}\,dx\,dy$; for $t_R=cR^2$, the rescaled average is, in law, a non-degenerate positive random variable $u_c(B_1)$; and for $t_R^{-1}=o(R^{-2})$, the average converges to zero in probability. The interest is that the middle regime is a non-Gaussian spatial-average limit, produced by all Wiener chaos orders contributing at exactly the diffusive time scale.

What carries the argument

The load-bearing identity is the scaling law $u_t(B_R)=\varepsilon^{-d/2}u_{\varepsilon t}(B_{\sqrt{\varepsilon}R})$ in law, which follows from the noise scaling $\dot F(t,x)=\varepsilon^2\dot F(\varepsilon^2 t,\varepsilon x)$ together with uniqueness in law of the martingale-problem solution. This identity converts every time-dependent statement into a fixed-time statement: part (1) reduces to a central limit theorem for $u_1(B_{\tilde R})$ as $\tilde R\to\infty$, part (2) is the exact identity $R^{-d}u_{cR^2}(B_R)=u_c(B_1)$, and part (3) reduces to the known extinction of $u_t(B_1)$ as $t\to\infty$. The quantitative work is done in Fourier space: the variance of the $n$-th Wiener chaos is expressed through the Riesz spectral measure $c_d\|\xi\|^{-d+2}d\xi$, the Fourier transform of $\mathbf{1}_{B_R}$ in terms of the Bessel function $J_{d/2}$, and an auxiliary function $\varphi(\eta^n)$ encoding the time integrals; Lemmas 3.7 and 3.8 show that the summed higher-chaos variances vanish under $\kappa<(d-2)/2$.

What would settle it

For two different radii $R_1$ and $R_2$, compare the law of $R^{-d}u_{cR^2}(B_R)$: part (2) predicts these laws coincide for every $R$, so any detectable dependence on $R$ would disprove the scaling identity and with it the phase diagram. A second check is to evaluate the summed higher-chaos variances in Proposition 3.3 at $\kappa=(d-2)/2$; if they fail to vanish, the Gaussian regime genuinely stops at the stated upper bound.

Watch

Extended reading notes

Core claim

The paper's central claim is Theorem 1.1: for every $0<\kappa<(d-2)/2$, the time-dependent spatial averages of the solution constructed in [MT04] have three different limits according to how $t_R$ compares with $R^2$. If $t_R=o(R^2)$, then $(t_R R^{d-1})^{-1/2}(u_{t_R}(B_R)-\omega_d R^d)$ converges in law to a Gaussian with the variance $\sigma^2$ given above. If $t_R=cR^2$ for $c>0$, then $R^{-d}u_{t_R}(B_R)=u_c(B_1)$ in law for every $R>0$, where $u_c(B_1)$ is positive almost surely and non-degenerate. If $t_R^{-1}=o(R^{-2})$, then $R^{-d}u_{t_R}(B_R)$ converges to $0$ in probability. The Gaussian regime is driven by the first Wiener chaos being dominant, while the non-Gaussian regime appears because all chaos orders contribute at the critical space-time scale.

Load-bearing premise

The entire phase diagram rests on the scaling law $u_t(B_R)=\varepsilon^{-d/2}u_{\varepsilon t}(B_{\sqrt{\varepsilon}R})$ in law; if the rescaled random-measure process were not itself a solution of the same equation, or if the solution's law were not unique, the three regimes would not follow.

Editorial extensions

If this is right

  • For all $0<\kappa<(d-2)/2$, the full space-time phase diagram of the ball average is: sub-diffusive time gives Gaussian fluctuations, diffusive time gives a $c$-dependent non-Gaussian law, and super-diffusive time gives extinction.
  • The Gaussian regime is equivalent, by the scaling law, to a central limit theorem for the fixed-ball random variable $u_t(B_1)$ as $t\to 0$.
  • The limit in the diffusive regime is the first non-Gaussian spatial-average limit for an SPDE in the literature surveyed by the paper.
  • The contrast with Riesz index $p\in(2,d)$ says that $p=2$ is the only index at which all Wiener chaos orders survive the diffusive scaling, so the large-scale fluctuations remain multiplicative rather than becoming Edwards-Wilkinson.
  • The boundary case $\kappa=(d-2)/2$ is left open for the central limit theorem, while parts (2) and (3) are expected to hold for the chaos-expansion solution.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the scaling identity is as general as it looks, the same phase diagram should be testable for the chaos-expansion solution at $\kappa=(d-2)/2$, with the Gaussian regime the part most likely to fail first.
  • Between the three regimes, for example $t_R=R^{2-\delta}$ with $0<\delta<2$, one would expect a family of intermediate limits interpolating between the Gaussian law and $u_c(B_1)$, though the paper does not address this.
  • A quantitative central limit theorem for part (1) would require controlling the decay rate of the summed higher-chaos variances without pointwise Malliavin derivatives; the paper's Fourier-space estimates are a natural starting point.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 3 minor

Summary. The paper studies the time-dependent spatial averages of a critical stochastic heat equation in d≥3 with noise white in time and colored in space with covariance kernel ∥x−y∥^{−2}. The solution is the measure-valued process constructed by Mueller and Tribe, and the main result is a three-regime phase diagram for the average over a ball of radius R at time t_R: for t_R=o(R^2) the normalized average converges to a centered Gaussian with variance σ²; for t_R=cR^2 the average equals u_c(B_1) in law, a positive non-Gaussian random variable; for t_R^{−1}=o(R^{−2}) the average tends to zero in probability. The proof uses the scaling invariance of the equation to reduce parts (2) and (3) to previous results of Mueller and Tribe, and to reduce part (1) to the fixed-time t=1 case. The fixed-time CLT is proved via the Wiener chaos decomposition, with the first-chaos variance computed explicitly and the higher-chaos variance sum shown to vanish.

Significance. If the proof is completed, this is a valuable contribution: it gives a complete phase diagram for spatial averages of a critical, measure-valued SPDE, and it appears to provide the first example of a non-Gaussian spatial-average limit for an SPDE, obtained through an exact scaling identity instead of a fitted limit. The paper is also honest about its scope, emphasizing the restriction κ<(d−2)/2 and the open boundary case. The proof strategy is natural and the first-chaos computation is clean; the higher-chaos estimates are new and use explicit constants and dominated convergence rather than Dalang-type integrability conditions. The main weakness is a concrete gap in Lemma 3.8, which is load-bearing for the central limit theorem.

major comments (2)
  1. [Section 3.2.3, Lemma 3.8, display between (3.15) and (3.16)] The displayed inequality asserting that the integral in (3.15) is bounded by the integral with product ∏_{j=2}^{n}(∥η_j−η_{j−1}∥^{−d+2}∥η_j∥^{−2})∥η_n∥^{−d+2} is false for m>0: the right-hand side contains extra factors ∥η_j∥^{−2} for j=2,…,m+1 that are not present on the left, and these factors can be arbitrarily small. For example, with m=1 and η2 large, the asserted pointwise inequality fails. Since this bound is the starting point for the convolution argument that yields (3.16) and the summability of κ^{2n}c_d^n C′(n,d), Proposition 3.3 — and hence Theorem 1.1(1) — is not established as written. The argument appears repairable by retaining the factors ∏_{j=2}^{m+1}∥η_j−η_{j−1}∥^{−d+2} in the bound and absorbing the factors ∥η_i∥^{−2} (i=m+2,…,n−1) one at a time in the successive applications of Lemma 2.1; the constants in (3.16) are consistent with such a repaired derivation, but the correction must be supplied.
  2. [Section 3.2.3, Lemma 3.8, definition of γ′₀ after (3.17)] The formula γ′₀ = 3/2 + 4m₀ − (d−3)/2 can exceed 2 (e.g., d=7 gives γ′₀=7/2), contradicting the stated condition 0≤γ′₀<2 and the requirement 0<2−γ′₀. The correct expression implied by (3.17) is γ′₀ = 2m₀ + 3 − d/2 (equivalently, γ′₀ = 3/2 + (4m₀ − (d−3))/2). As printed, the verification that the convolution kernels are integrable and that the resulting series is summable for all κ<(d−2)/2 is not justified. This is part of the same load-bearing gap in Lemma 3.8.
minor comments (3)
  1. [Section 3.1, Lemma 3.1] The proof of the scaling identity omits the verification that the rescaled measure-valued process {v^ε_t} satisfies the martingale problem (2.2)–(2.3) with the same coupling constant κ and initial data. The claim is indeed straightforward from (1.6), but since uniqueness in law is what makes the identity exact, a one-line computation (or the precise statement of [MT04, Lemma 2]) should be included.
  2. [Section 3.2.1, Lemma 3.5] The phrase 'change-of-variable η_j → η_j/R' is ambiguous and can lead the reader to question the factor R^{2d−2n} and the argument of φ in the final formula. The authors should state explicitly that the new variables are η_j^{new} = R η_j^{old} (so that ∫ dη_old = R^{−nd} dη_new) and then display the resulting expression.
  3. [Section 3.2.3, equation (3.14)] Equation (3.14) applies the bound (3.8) to the block of variables η_{m+2},…,η_{n−1} rather than to the initial block. This is valid because the simplex S_n is symmetric under permutations of the variables, but the paper should state this symmetry explicitly when deriving (3.14).

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the derivation uses external results from [MT04] and exact scaling identities; no fitted parameter is renamed as a prediction and no self-citation is load-bearing.

full rationale

The paper's results are not circular. Theorem 1.1 part (2) is an exact distributional identity obtained from the scaling relation Lemma 3.1, namely u_t(B_R) = ε^{-d/2} u_{εt}(B_{√ε R}) in law, which yields u_tR(B_R)/R^d = u_{cR^2/R^2}(B_1) = u_c(B_1) when t_R = cR^2. This is a consequence of the scaling invariance of the equation and the uniqueness-in-law theorem proved in the prior external paper [MT04], not of a fitted parameter or of the present authors' own results. Part (3) follows from the externally established extinction result [MT04, Theorem 2 ii] together with the same scaling identity; no quantity is defined in terms of the limit it is supposed to explain. Part (1) is proved by reducing the time-dependent average to a fixed-time average at t=1 via the scaling lemma, then computing the variance of the first chaos explicitly (Proposition 3.2) and showing the sum of all higher-order chaos variances vanishes (Proposition 3.3), with the limiting Gaussian variance σ² computed from the noise covariance by Fourier analysis. No parameter is fitted to any subset of the data and no prediction is statistically forced by an input. The paper's reliance on [MT04] is a reference to independent prior work with no author overlap with the present paper; the existence, uniqueness in law, and extinction results are stated assumptions from that external source, not self-citations. The skeptical observation about a possible index inconsistency around equation (3.16) in Lemma 3.8 concerns the correctness of a technical estimate in the written proof, not circularity, and therefore does not affect the circularity score under the review rules.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central claim rests on the Mueller-Tribe solution theory and on the scaling property of the solution, not on fitted parameters or invented entities. All additional ingredients are standard Fourier analysis and Riesz potential estimates.

assumptions (4)
  • domain assumption Existence, uniqueness in law, and measure-valuedness of the martingale-problem solution for (1.1) when 0<kappa<(d-2)/2, from [MT04].
    The paper builds on the Mueller-Tribe solution theory; the spatial average u_t(B_R) presupposes the solution is a nonnegative Radon measure-valued process, and uniqueness in law is used to identify the scaling limit.
  • domain assumption Scaling invariance of the noise and of the solution process (Lemma 3.1), including the transfer of scaling via uniqueness in law.
    The proof of parts (2) and (3) and the reduction of part (1) rely on u_t(B_R)=epsilon^{-d/2}u_{epsilon t}(B_{sqrt(epsilon) R}) in law; the verification that the rescaled process is a solution is omitted and referenced to [MT04, Lemma 2].
  • standard math Riesz kernel composition formula (Lemma 2.1) and Bessel function asymptotics (Lemma 2.2).
    Used in the variance computations and higher-chaos estimates; the proof of Lemma 2.1 is omitted as an elementary gamma-beta computation.
  • standard math Finite-dimensional distribution and Brownian bridge computations leading to the Feynman-Kac type formula (2.5) for E[beta_t(x-y)^n].
    Used to derive the Fourier-space variance formula in Lemma 3.5; standard Ito isometry plus Brownian bridge arguments are invoked.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Time-dependent averages of a critical long-range stochastic heat equation." pith.science (2026). https://pith.science/paper/YJX6CRSC

@misc{pith2026241109058,
  author       = {Pith},
  title        = {Pith review of: Time-dependent averages of a critical long-range stochastic heat equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YJX6CRSC}},
  note         = {Machine review of arXiv:2411.09058}
}
abstract

We study the time-dependent spatial averages of a critical stochastic partial differential equation, namely the stochastic heat equation in dimension $d\geq 3$ with noise white in time and colored in space with covariance kernel $\|\cdot\|^{-2}$. The solution to this SPDE is a singular measure and was constructed by Mueller and Tribe in [MT04]. We show that the time-dependent spatial averages of this SPDE over a ball of radius $R$ at time $t$ have different limits under different space-time scales. In particular, when $t\ll R^2$, the central limit theorem holds; when $t=R^2$, the spatial average is a non-Gaussian random variable; when $t\gg R^2$, the spatial average becomes extinct.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

23 extracted references · 18 canonical work pages

  1. [1]

    Averaging 2d stochastic wave equation

    Raul Bola \ n os Guerrero, David Nualart, and Guangqu Zheng. Averaging 2d stochastic wave equation. Electronic Journal of Probability , 26:1--32, 2021

    work page 2021

  2. [2]

    Central limit theorems for spatial averages of the stochastic heat equation via M alliavin- S tein's method

    Le Chen, Davar Khoshnevisan, David Nualart, and Fei Pu. Central limit theorems for spatial averages of the stochastic heat equation via M alliavin- S tein's method. Stoch. Partial Differ. Equ. Anal. Comput. , 11(1):122--176, 2023

    work page 2023

  3. [3]

    Universality in marginally relevant disordered systems

    Francesco Caravenna, Rongfeng Sun, and Nikos Zygouras. Universality in marginally relevant disordered systems . The Annals of Applied Probability , 27(5):3050--3112, 2017

    work page 2017

  4. [4]

    The critical 2d stochastic heat flow

    Francesco Caravenna, Rongfeng Sun, and Nikos Zygouras. The critical 2d stochastic heat flow. Inventiones mathematicae , 233:325--460, 2023

    work page 2023

  5. [5]

    Robert C. Dalang. Extending the martingale measure stochastic integral with applications to spatially homogeneous s.p.d.e.'s. Electron. J. Probab. , 4:no. 6, 29, 1999

    work page 1999

  6. [6]

    Flow equation approach to singular stochastic pdes

    Pawe Duch. Flow equation approach to singular stochastic pdes. arXiv preprint arXiv:2109.11380 , 2021

    arXiv 2021

  7. [7]

    A Central Limit Theorem for the stochastic wave equation with fractional noise

    Francisco Delgado-Vences, David Nualart, and Guangqu Zheng. A Central Limit Theorem for the stochastic wave equation with fractional noise . Annales de l'Institut Henri Poincaré, Probabilités et Statistiques , 56(4):3020 -- 3042, 2020

    work page 2020

  8. [8]

    Fluctuations of stochastic pdes with long-range correlations

    Luca Gerolla, Martin Hairer, and Xue-Mei Li. Fluctuations of stochastic pdes with long-range correlations. arXiv preprint arXiv:2303.09811 , 2023

    arXiv 2023

Show all 23 references
  1. [9]

    Paracontrolled distributions and singular PDE s

    Massimiliano Gubinelli, Peter Imkeller, and Nicolas Perkowski. Paracontrolled distributions and singular PDE s. Forum Math. Pi , 3:e6, 75, 2015

    2015

  2. [10]

    Fluctuations of a nonlinear stochastic heat equation in dimensions three and higher

    Yu Gu and Jiawei Li. Fluctuations of a nonlinear stochastic heat equation in dimensions three and higher. SIAM Journal on Mathematical Analysis , 52(6):5422--5440, 2020

    2020

  3. [11]

    urgen G\

    J\"urgen G\"artner and Adrian Schnitzler. Stable limit laws for the parabolic A nderson model between quenched and annealed behaviour. Ann. Inst. Henri Poincar\'e Probab. Stat. , 51(1):194--206, 2015

    2015

  4. [12]

    A theory of regularity structures

    Martin Hairer. A theory of regularity structures. Inventiones mathematicae , 198(2):269--504, 2014

    2014

  5. [13]

    A central limit theorem for the stochastic heat equation

    Jingyu Huang, David Nualart, and Lauri Viitasaari. A central limit theorem for the stochastic heat equation. Stochastic Processes and Their Applications , 130(12):7170--7184, 2020

    2020

  6. [14]

    Gaussian fluctuations for the stochastic heat equation with colored noise

    Jingyu Huang, David Nualart, Lauri Viitasaari, and Guangqu Zheng. Gaussian fluctuations for the stochastic heat equation with colored noise. Stochastics and Partial Differential Equations: Analysis and Computations , 8:402--421, 2020

    2020

  7. [15]

    Foundations of Stochastic Differential Equations in Infinite Dimensional Spaces , volume 47 of CBMS-NSF Regional Conference Series in Applied Mathematics

    Kiyosi Ito. Foundations of Stochastic Differential Equations in Infinite Dimensional Spaces , volume 47 of CBMS-NSF Regional Conference Series in Applied Mathematics . SIAM, 1984

    1984

  8. [16]

    Limit theorems for time-dependent averages of nonlinear stochastic heat equations

    Kunwoo Kim and Jaeyun Yi. Limit theorems for time-dependent averages of nonlinear stochastic heat equations. Bernoulli , 28(1):214--238, 2022

    2022

  9. [17]

    Nikolai N. Lebedev. Special Functions and Their Applications . Dover Publications, New York, revised english edition edition, 1972

    1972

  10. [18]

    A singular parabolic A nderson model

    Carl Mueller and Roger Tribe. A singular parabolic A nderson model. Electron. J. Probab. , 9:no. 5, 98--144, 2004

    2004

  11. [19]

    Spatial averages for the parabolic A nderson model driven by rough noise

    David Nualart, Xiaoming Song, and Guangqu Zheng. Spatial averages for the parabolic A nderson model driven by rough noise. ALEA Lat. Am. J. Probab. Math. Stat. , 18(1):907--943, 2021

    2021

  12. [20]

    The M alliavin calculus and related topics

    David Nualart. The M alliavin calculus and related topics . Probability and its Applications (New York). Springer-Verlag, New York, 1995

    1995

  13. [21]

    Quantitative central limit theorems for the parabolic A nderson model driven by colored noises

    David Nualart, Panqiu Xia, and Guangqu Zheng. Quantitative central limit theorems for the parabolic A nderson model driven by colored noises. Electron. J. Probab. , 27:Paper No. 120, 43, 2022

    2022

  14. [22]

    Averaging G aussian functionals

    David Nualart and Guangqu Zheng. Averaging G aussian functionals. Electron. J. Probab. , 25:Paper No. 48, 54, 2020

    2020

  15. [23]

    A priori bounds for quasi-linear spdes in the full sub-critical regime

    Felix Otto, Jonas Sauer, Scott Smith, and Hendrik Weber. A priori bounds for quasi-linear spdes in the full sub-critical regime. arXiv preprint arXiv:2103.11039 , 2021

    2021 arXiv

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.