{"id":"e842734c-8a25-4f16-9fc6-44643ffc8edc","arxiv_id":"2504.18450","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Temporal quadratic and higher-order variations of the nonlinear fractional stochastic heat equation converge to explicit integrals of the diffusion coefficient, yielding consistent estimators for the drift and the fractional Laplacian index.","lead":"A mathematics paper proves that time-based variation sums of the solution to a nonlinear stochastic heat equation converge to explicit limits, and uses them to estimate two unknown parameters. A smart generalist might read it because it extends statistical inference for random partial differential equations from linear additive noise to multiplicative nonlinear noise.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Proof gap: 'far' modified increments are not independent of u(t_i(δ),x); the factorization in l^(2,2)_{2,N} is unjustified.","rationale":"The reader's verification was careful and the acceptance is understandable: the limit statements are plausible, the rates match known α=2 results, and the even-integer restriction is disclosed. However, the stress-test found an internal gap in the proof of both main theorems that the reader did not flag. In the estimate of L_{2,N} (Theorem 1) and A_{2,N} (Theorem 2), the 'far' pairs are defined by i−j ≥ [N^{1−β}+1]+2, so t_j+δ ≤ t_i(δ). The paper claims this makes both modified increments independent of σ²(u(t_i(δ),x)) and σ²(u(t_j(δ),x)). But ~Δ_1u0(t_j,δ) is built from the original white noise on [t_j(δ),t_j+δ], and u(t_i(δ),x) depends on that same noise because t_j+δ ≤ t_i(δ). The covariance is nonzero already in the linear case, so the factorization used to apply Lemma 2 is not justified. Without an estimate for the additional cross-covariance term, the stated L1 rates do not follow from the written proof. This is a repairable gap rather than a demonstrated false theorem: the same bound might be recoverable by a direct Green-kernel estimate of the shared-noise covariance. I therefore recommend CONDITIONAL acceptance, with the condition being the addition of a rigorous treatment of this term. If the missing estimate is supplied and the rates survive, the paper's central claims and statistical applications stand.","tokens_in":25035,"tokens_out":23464,"duration_ms":212578,"concrete_test":"Check the independence claim directly in the linear case σ≡1. Fix i,j with i−j of order N^{1−β} so that t_j+δ ≤ t_i(δ), and compute Cov(~Δ_1u0(t_j,δ), u0(t_i(δ),x)) = ∫_{t_j(δ)}^{t_j+δ}∫_R G(t_i(δ)−a,x−y)(G(t_j+δ−a,x−y)−G(t_j−a,x−y)) dy da. Since the integrand is nonnegative and the time interval is nonempty, this covariance is positive, disproving the asserted independence. Then re-derive the contribution of the 'far' pairs to E|L_{2,N}|^2 without the factorization, bounding E[σ²(u(t_i(δ)))σ²(u(t_j(δ))) ~Δ_i ~Δ_j] conditionally on the shared noise on [t_j(δ),t_j+δ]. If the resulting bound is O(N^{−1}) or O(N^{−(α+1)/(2α)}), the theorem survives; if it is O(1), the stated rate (35) is not established.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"In the estimate of L_{2,N} in Theorem 1 (and A_{2,N} in Theorem 2), the authors split the double sum at i−j ≥ [N^{1−β}+1]+2, so that t_j+δ ≤ t_i(δ), and assert that this 'implies that ~Δ_1u0(t_i,δ) and ~Δ_1u0(t_j,δ) are independent of σ²(u(t_i(δ),x)) and σ²(u(t_j(δ),x))'. This implication is false. By construction, ~Δ_1u0(t_j,δ) uses the original noise W on [t_j(δ), t_j+δ]. Since t_j+δ ≤ t_i(δ), the random variable u(t_i(δ),x) also depends on W on exactly that interval. Already in the linear case σ≡1, Cov(~Δ_1u0(t_j,δ), u0(t_i(δ),x)) = ∫_{t_j(δ)}^{t_j+δ}∫_R G(t_i(δ)−a,x−y)(G(t_j+δ−a,x−y)−G(t_j−a,x−y)) dy da, which is not identically zero. Hence the factorization of l^(2,2)_{2,N} into E[σ²(u(t_i))σ²(u(t_j))]·E[~Δ_i~Δ_j] is unjustified, and the subsequent application of Lemma 2 bound (15) only controls the second factor. The cross-covariance term created by the shared noise is missing from the proof; without a bound on it, the L1 rate (35) does not follow from the written argument. The same gap affects a^(2,2)_{2,N} in Theorem 2 and rate (52). The limit theorems may still be true, but this step requires an additional estimate.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the temporal power variations of the solution to the nonlinear stochastic heat equation with a fractional Laplacian, ∂u/∂t = -(-Δ)^{α/2} u + σ(u) Ẇ, for α∈(1,2], observed at a fixed spatial point over an equidistant time partition. The main theoretical results are Theorem 1, which states that the renormalized quadratic variation N^{-1/α} Σ (u(t_{i+1},x)-u(t_i,x))² converges in L¹ to C_{0,α}² ∫_0^1 σ²(u(s,x)) ds with explicit rate (35), and Theorem 2, which states that, when p=2α/(α-1) is an even integer, the p-th power variation converges to B_{0,α} ∫_0^1 σ^p(u(s,x)) ds with rate (52). The proofs compare the increments of the nonlinear solution with those of the linear equation via a modified increment ~Δ₁u₀ that is independent of the past after a time shift, and use the representation of the linear solution as a perturbed fractional Brownian motion. The last section applies these limits to construct a consistent estimator for the anomality parameter α and two consistent estimators for the drift θ. The paper is clearly written and the rates are internally consistent, reproducing the known N^{-3/20} rate at α=2. However, the central decoupling step in the proof of Theorem 1, and its analogue in Theorem 2, contains a false independence assertion; the proofs of the limit theorems are therefore incomplete as written.","tokens_in":25443,"tokens_out":11505,"duration_ms":105015,"significance":"If the gap identified below is repaired, the paper would be a valuable contribution: it gives the first estimator for the fractional-Laplacian index α from discrete temporal observations at one spatial point, and it extends exact-variation results from the linear or standard-Laplacian setting to the nonlinear fractional setting with explicit convergence rates. The constants C_{0,α} and B_{0,α} are taken from the known linear theory and are not fitted, and the estimators are genuine functions of the data, so the statistical claims are not circular. The appendix estimates, in particular Proposition 4 and Lemma 4, appear sound and document the kernel estimates carefully. The main obstacle is a single but load-bearing independence argument in the proofs of Theorems 1 and 2; if that step cannot be justified, the stated L¹ rates and the consistency conclusions do not follow from the written proof.","major_comments":[{"comment":"The decoupling of l^(2,2)_{2,N} is not justified. The text asserts that when i-j ≥ [N^{-β}+1]+2 (so that t_j+δ ≤ t_i(δ)), the modified increments ~Δ₁u₀(t_i,δ) and ~Δ₁u₀(t_j,δ) are independent of σ²(u(t_i(δ),x)) and σ²(u(t_j(δ),x)). This is false. By construction, ~Δ₁u₀(t_j,δ) uses the original noise W on [t_j(δ), t_j+δ], and because t_j+δ ≤ t_i(δ), the random variable u(t_i(δ),x) also depends on W on exactly that interval. Already for the linear field, Cov(~Δ₁u₀(t_j,δ), u₀(t_i(δ),x)) = ∫_{t_j(δ)}^{t_j+δ}∫_R G_α(t_i(δ)-a,x-y)(G_α(t_j+δ-a,x-y)-G_α(t_j-a,x-y)) dy da, which is not identically zero; the same holds for the nonlinear solution. Consequently, E[σ²(u(t_i(δ),x))σ²(u(t_j(δ),x)) Z_i Z_j] does not factor as E[σ²(u(t_i(δ),x))σ²(u(t_j(δ),x))] E[Z_i Z_j], and the bound (38) does not follow from Lemma 2. Since this term feeds directly into E|L_{2,N}|² and therefore into the rate (35), Theorem 1 is not proved as written. The identical defect appears in the treatment of a^(2,2)_{2,N} in the proof of Theorem 2 and affects (52). The limit statements may well be true, but a supplementary estimate controlling the cross-covariance between ~Δ₁u₀(t_j,δ) and the σ(u(t_i(δ),x))-dependent factors is required.","section":"Section 3.2, proof of Theorem 1, estimate of l^(2,2)_{2,N}"}],"minor_comments":[{"comment":"The bound E|~Δ₁u₀(t_i,δ)|² ≤ C N^{-(α-1)/(2α)} has the exponent off by a factor of 2 for a second moment; the next line uses C N^{-(α-1)/α}, which is the correct order for the second moment of a temporal increment of u₀. Please correct the displayed bound.","section":"Section 3.2, Eq. (36)"},{"comment":"The condition for t_j+δ ≤ t_i(δ) is i-j ≥ N^{1-β}+1, not i-j ≥ [N^{-β}+1]+2; as printed, the 'far' set includes pairs that do not have the stated time separation. Since the subsequent estimate passes to the full double sum, this typo alone is not fatal, but it should be corrected.","section":"Sections 3.2 and 3.3, split threshold"},{"comment":"In the displayed formula for a^(2,2)_{2,N}, both factors in the product are written with index j: [ (~Δ₁u₀(t_j,δ))^p - B_{0,α}/N ][ (~Δ₁u₀(t_j,δ))^p - B_{0,α}/N ]. One factor should involve ~Δ₁u₀(t_i,δ).","section":"Section 3.3, proof of Theorem 2, a^(2,2)_{2,N}"},{"comment":"There are several typographical errors: 'spce-time' in the first paragraph of Section 3.3, 'we need need' before Lemma 3, 'Altough' in the proof of Proposition 1, and a missing parenthesis in reference [14]. These do not affect the mathematics.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read this paper carefully. The genuinely new result is an estimator for the fractional Laplacian index alpha in the nonlinear multiplicative-noise equation, built from a renormalized temporal quadratic variation. That estimation problem is new as far as I know, and the estimator is natural: a log-ratio of normalized variation to a Riemann sum of sigma^2. Theorems 1 and 2 also extend temporal variation limits from the linear to the nonlinear setting, and the rates check out at alpha=2 against the known Pospisil-Tribe result. The comparison-to-linear-solution strategy is the right high-level idea, and Proposition 4 gives a clean L2 bound on the replacement error.\n\nBut there is a load-bearing gap. In the proof of Theorem 1, the 'far' part l(2,2) is handled by asserting that ~Delta_i and ~Delta_j are independent of the sigma^2(u(t_i(delta))) and sigma^2(u(t_j(delta))) factors when t_j+delta <= t_i(delta). That is false. The modified increment ~Delta_j is defined using the original noise W on [t_j(delta), t_j+delta], and u(t_i(delta)) depends on W on exactly that interval. Using an independent copy ~W on the past does not remove that shared W. So the factorization E[sigma^2(u(t_i))sigma^2(u(t_j))] * E[Delta_i Delta_j] is unjustified. The same flaw appears in the a(2,2) part of Theorem 2. Without an extra estimate bounding the cross-covariance between sigma^2(u(t_i(delta))) and the modified increment of an earlier block, the stated L1 rates (35) and (52) do not follow from the written argument. This is not a typo; it is the step that converts the comparison bound into the limit theorem. It might be fixable—the cross-covariance should be small for large i-j—but the paper does not do it.\n\nOther weaknesses are minor by comparison. The even-integer condition on 2alpha/(alpha-1) is restrictive but the authors disclose it. The heavy reliance on earlier work from the same group is noticeable but the citations are relevant and not padding. The writing is clear and the appendix contains the necessary kernel estimates.\n\nThe paper is for researchers working on high-frequency statistics for SPDEs, especially anyone trying to estimate the fractional Laplacian index from discrete temporal observations. I would send it to a serious referee, because the problem is interesting and the gap may be closable, but I would not accept it as is; the referee should focus on the independence claim and request the missing estimate.","headline":"New and useful estimator for alpha in the nonlinear fractional heat equation, but a load-bearing independence claim in the proof of the variation limits is not justified as written; the paper deserves a serious referee if that gap can be closed.","tokens_in":25916,"tokens_out":5021,"would_cite":false,"duration_ms":48356,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60G15","60H05","60G18"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that temporal power variations of the nonlinear fractional heat equation converge to explicit functionals of the solution, and uses these limits to build consistent estimators for the drift $\\theta$ and the…","keywords":["stochastic heat equation","fractional Laplacian","fractional Brownian motion","power variation","quadratic variation","parameter estimation","nonlinear multiplicative noise","consistency"],"falsifier":"Simulate the linear solution $u_0$ for a fixed $\\alpha\\in(1,2)$ and directly measure $E|Y_{i+1}-Y_i|^2$ for the perturbation term in (28); if the decay is slower than $i^{-(\\alpha+1)/\\alpha}$, the rates in Theorems 1–2 and the consistency proofs collapse. Alternatively, for $\\alpha=3/2$, check numerically that $N^{-1/\\alpha}\\sum_{i=0}^{N-1}(u_\\theta(t_{i+1},x)-u_\\theta(t_i,x))^2$ approaches $C_{0,\\alpha}^2\\theta^{-1/\\alpha}\\int_0^1\\sigma^2(u_\\theta(s,x))\\,ds$; a systematic mismatch would refute the main limit theorem.","tokens_in":24830,"feed_emoji":"🧮","tokens_out":6954,"duration_ms":64563,"temperature":0.7,"pith_summary":"The paper establishes limit theorems for the temporal power variations of the solution to the nonlinear stochastic heat equation with a fractional Laplacian of order $\\alpha\\in(1,2]$, driven by spacetime white noise. Observing the solution at one fixed spatial point on an equidistant grid, it proves that a suitably renormalized quadratic variation converges in $L^1$ to an explicit constant times $\\theta^{-1/\\alpha}\\int_0^1\\sigma^2(u_\\theta(s,x))\\,ds$, and that, when the exponent $2\\alpha/(\\alpha-1)$ is an even integer, the corresponding power variation converges to a similar explicit integral. From these limits, the authors construct consistent estimators for the drift $\\theta$ and for the fractional-Laplacian order $\\alpha$. The point is that unknown parameters of a nonlinear SPDE can be recovered from the roughness of a single temporal path observed at discrete times.","feed_headline":"Fractional heat equation variations converge and reveal its parameters","feed_subtitle":"From discrete time samples at one spatial point, the drift and the fractional Laplacian order can be estimated consistently.","key_machinery":"The engine is the representation of the linear solution's time process as a perturbed fractional Brownian motion (equation (28)): $u_0(t,x)=C_{0,\\alpha}U_t+Y_t$, where $U$ is a fractional Brownian motion with Hurst index $H=(\\alpha-1)/(2\\alpha)$ and $Y$ is a self-similar Gaussian perturbation whose increments satisfy $E|Y_{i+1}-Y_i|^2\\le C i^{-(\\alpha+1)/\\alpha}$ (bound (29)). Around this, Proposition 4 compares true nonlinear increments with modified increments of the linear solution at shifted times $t_i(\\delta)=t_i-\\delta^{\\beta}$, $\\beta=2\\alpha/(2\\alpha+1)$, with error $O(\\delta^{4(\\alpha-1)/(2\\alpha+1)})$. The variation theorems then split the error into a nonlinear comparison term, a Gaussian fluctuation term controlled by the perturbed-fBm bounds in Lemmas 2 and 3, and a Riemann-sum term controlled by the solution's Hölder regularity.","core_discovery":"For the solution $u_\\theta$ of $\\partial_t u_\\theta = -\\theta(-\\Delta)^{\\alpha/2}u_\\theta + \\sigma(u_\\theta)\\dot W$ with zero initial condition, observed at times $t_i=i/N$ and a fixed spatial point, the renormalized temporal quadratic variation $N^{-1/\\alpha}\\sum_{i=0}^{N-1}(u_\\theta(t_{i+1},x)-u_\\theta(t_i,x))^2$ converges in $L^1$ to $C_{0,\\alpha}^2\\theta^{-1/\\alpha}\\int_0^1\\sigma^2(u_\\theta(s,x))\\,ds$. When $2\\alpha/(\\alpha-1)$ is an even integer, the power variation $\\sum_{i=0}^{N-1}|u_\\theta(t_{i+1},x)-u_\\theta(t_i,x)|^{2\\alpha/(\\alpha-1)}$ converges to $B_{0,\\alpha}\\theta^{-1/(\\alpha-1)}\\int_0^1\\sigma^{2\\alpha/(\\alpha-1)}(u_\\theta(s,x))\\,ds$. These two limits are the basis for consistent in-probability estimators of $\\alpha$ and $\\theta$: $\\hat\\alpha_N = \\log N/\\log A_N$, $\\hat\\theta_{N,1}$ from the quadratic variation, and $\\hat\\theta_{N,2}$ from the higher-order variation.","pith_inferences":["The even-integer condition on $2\\alpha/(\\alpha-1)$ appears to be technical; a central-limit or noninteger-power extension would make $\\hat\\theta_{N,2}$ available for all $\\alpha\\in(1,2)$, not only the discrete set where the exponent is even.","Because $\\hat\\alpha_N$ is independent of $\\theta$, a natural two-step inference procedure—estimate $\\alpha$ first, then $\\theta$—would decouple the parameters; the asymptotic joint distribution of such a procedure is not addressed in the paper.","If $\\sigma$ vanishes on a portion of the observed path, the denominators in (62)–(64) could degenerate; handling such cases would require a modified normalization that the paper does not discuss.","The same perturbed-fractional-Brownian strategy could be applied to spatial power variations at a fixed time, potentially yielding estimators that identify $\\alpha$ without relying on temporal sampling rates."],"forward_implications":["The estimator $\\hat\\alpha_N=\\log N/\\log A_N$ recovers the fractional-Laplacian order $\\alpha$ in probability without requiring knowledge of the drift $\\theta$.","The drift estimators $\\hat\\theta_{N,1}$ and $\\hat\\theta_{N,2}$ are consistent in probability for $\\theta$ when $\\alpha$ is known.","When $\\alpha=2$, the higher-order variation theorem yields a nontrivial quartic variation with error rate $N^{-3/20}$, matching the known standard-Laplacian case.","The limits extend to any time interval $[A_1,A_2]$, with the quadratic limit acquiring a factor $(A_2-A_1)^{2H-1}$.","The quadratic and higher-order variations converge at the same $L^1$ rate, so either observable can serve as the basis for parameter inference."],"supporting_citations":[{"why":"Supplies the increment bound (29) for the Gaussian perturbation in the linear solution's temporal decomposition, which sets the exponents used in Lemmas 2–3 and Theorems 1–2.","marker":"[4]"},{"why":"Cited for the representation (28) of the linear solution as a perturbed fractional Brownian motion and for the linear power-variation limit (42).","marker":"[29]"},{"why":"Provides the fBm quadratic-variation bounds in Lemma 1 and Proposition 1 and the spatial analogue whose arguments are adapted to the temporal setting.","marker":"[16]"},{"why":"Supplies the $\\alpha=2$ version of the nonlinear increment comparison and the quartic-variation rate cited in Remark 1.3.","marker":"[26]"},{"why":"Cited for the linear higher-order power-variation limit (42) and for the drift estimator in the linear case.","marker":"[23]"},{"why":"Establishes existence, moment bounds, and Hölder continuity (7)–(8) for the mild solution, which are used throughout the proofs.","marker":"[3]"},{"why":"Shows why integer powers enter the fBm power-variation analysis and motivates the even-integer restriction (44).","marker":"[24]"}],"fun_headline_variants":["Variation limits estimate fractional heat drift and Laplacian order","Convergent power variations reveal stochastic heat parameters","Quadratic and power variation limits identify heat equation parameters","From one point, variation limits recover fractional heat order and drift"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"All stated convergence rates inherit their exponents from the perturbation bound $E|Y_{i+1}-Y_i|^2\\le C i^{-(\\alpha+1)/\\alpha}$ in (29), and the higher-order variation theorem additionally assumes that $2\\alpha/(\\alpha-1)$ is an even integer.","fun_headline_variants_meta":{"raw":{"variants":["Variation limits estimate fractional heat drift and Laplacian order","Convergent power variations reveal stochastic heat parameters","Quadratic and power variation limits identify heat equation parameters","From one point, variation limits recover fractional heat order and drift"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000187,"raw_usage":{"total_tokens":1327,"prompt_tokens":941,"completion_tokens":386,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":320}},"tokens_in":557,"tokens_out":386,"duration_ms":4234,"temperature":1.0,"reasoning_tokens":320,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:16:35.592116+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate the linear solution $u_0$ for a fixed $\\alpha\\in(1,2)$ and directly measure $E|Y_{i+1}-Y_i|^2$ for the perturbation term in (28); if the decay is slower than $i^{-(\\alpha+1)/\\alpha}$, the rates in Theorems 1–2 and the consistency proofs collapse. Alternatively, for $\\alpha=3/2$, check numerically that $N^{-1/\\alpha}\\sum_{i=0}^{N-1}(u_\\theta(t_{i+1},x)-u_\\theta(t_i,x))^2$ approaches $C_{0,\\alpha}^2\\theta^{-1/\\alpha}\\int_0^1\\sigma^2(u_\\theta(s,x))\\,ds$; a systematic mismatch would refute the main limit theorem.","supporting_citations":[{"cited_title":"Bajja, K","cited_arxiv_id":null,"evidence_quote":"Supplies the increment bound (29) for the Gaussian perturbation in the linear solution's temporal decomposition, which sets the exponents used in Lemmas 2–3 and Theorems 1–2."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Cited for the representation (28) of the linear solution as a perturbed fractional Brownian motion and for the linear power-variation limit (42)."},{"cited_title":"Gamain and C","cited_arxiv_id":null,"evidence_quote":"Provides the fBm quadratic-variation bounds in Lemma 1 and Proposition 1 and the spatial analogue whose arguments are adapted to the temporal setting."},{"cited_title":"Pospisil and R","cited_arxiv_id":null,"evidence_quote":"Supplies the $\\alpha=2$ version of the nonlinear increment comparison and the quartic-variation rate cited in Remark 1.3."},{"cited_title":"Mahdi Khalil and C","cited_arxiv_id":null,"evidence_quote":"Cited for the linear higher-order power-variation limit (42) and for the drift estimator in the linear case."},{"cited_title":"Assaad, D","cited_arxiv_id":null,"evidence_quote":"Establishes existence, moment bounds, and Hölder continuity (7)–(8) for the mild solution, which are used throughout the proofs."},{"cited_title":"Nourdin, D","cited_arxiv_id":null,"evidence_quote":"Shows why integer powers enter the fBm power-variation analysis and motivates the even-integer restriction (44)."}],"review_version":1}