REVIEW 3 major objections 4 minor 1 cited by
Discrete reaction-diffusion system with stochastic dynamical boundary conditions: convergence results
T0 review · 3 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read This paper proves that a space-discrete approximation of a random reaction-diffusion system with stochastic boundary condition converges to the unique continuum solution.
desk verdict A real new convergence theorem for a semi-discrete scheme with stochastic dynamical boundary, but Proposition 43 has a fixable circular step and the fully discrete theorem overclaims. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The argument rests on three load-bearing pieces. First, the splitting $s=u+v$ (Proposition 21): $u$ is the solution of the discrete heat equation on the half-lattice with zero initial data and boundary value $\psi_t$, which absorbs the irregularity of the boundary path, while $v$ solves the nonlinear nonlocal system with deterministic initial and zero boundary data, coupled to $u$. Second, the regularization estimates for the discrete heat semigroup $e^{t\Delta_h}$ in Besov spaces (Lemmas 30-32), which transfer time regularity of $\psi$ into spatial regularity of $u$: the solution $u$ is bounded in $L^\infty_t W^{1,2}_x$ uniformly in $h$ (Proposition 34). Third, the linearized system (53)-(55) with a generic function $f$ in place of $s$, for which Propositions 38-42 give uniform $L^2$-type bounds; the nonlinear estimate (Proposition 43) is then obtained by formally taking $f=s$. Finally, the piecewise-constant extension operator $E_h$ (Definition 47) maps discrete Besov spaces continuously into continuum Besov spaces (Theorem 48), so the compactness argument of Theorem 56 yields convergence in $L^p([0,T],B^k_{p,p}(\mathbb{R}_+))$.
What would settle it
A concrete way to test the central claim is to compute $\|s_h\|_{L^2([0,T],H^1(\Omega_h^+))}$ for decreasing $h$ on the fully discrete scheme with a Pearson-process boundary (8) satisfying the paper's assumptions; if it is unbounded as $h\to0$, then Proposition 43 fails and the compactness argument of Theorem 56 has no uniform bound to work with. Conversely, a rigorous counterexample with initial data and a boundary path satisfying (7) for which this norm diverges would disprove the convergence theorem.
Extended reading notes
Core claim
The central claim is Theorem 56: let $(s_h,c_h)$ be the unique solution of the space-discrete system (72) with initial data converging in $H^1$ to $(s_0,c_0)$ and a boundary process $\psi$ satisfying the Hölder condition (7). Then the piecewise-constant extensions $(E_h s_h, E_h c_h)$ converge in $L^p([0,T], B^k_{p,p}(\mathbb{R}_+))$ for every $p\in[1,2]$ and every $k<1/p$ to the unique weak solution $(s,c)$ of the continuum system (67)-(68). The proof uses the splitting $s=u+v$ from Proposition 21, where $u$ solves the space-discrete heat equation with the random boundary condition and $v$ solves a nonlinear, nonlocal space-discrete system with deterministic boundary data; uniform-in-$h$ a priori estimates for both components are obtained through discrete heat-kernel regularization in Besov spaces and a linearization procedure, and compact embeddings of time-space Besov spaces on the lattice transfer the estimates to the continuum. The paper further proves that a fully discrete FTCS scheme converges to the space-discrete system with rate $O(k)$ for fixed $h$, and therefore to the continuum solution as $h\to0$ under the stability conditions (73)-(74).
Load-bearing premise
The proof of the uniform-in-$h$ a priori bound for the nonlinear system closes a circular step: Proposition 43 applies Proposition 42 with $f=s$, but condition (52) in Proposition 42 requires a norm bound on $f$ that is exactly the bound Proposition 43 is meant to produce, and the bootstrap from local existence (Theorem 46) to a global uniform estimate is not written out.
Editorial extensions
If this is right
- The semi-discrete scheme (25)-(26) is a convergent numerical method: its piecewise-constant extensions converge in $L^p([0,T], B^k_{p,p}(\mathbb{R}_+))$ to the unique continuum weak solution, giving a rigorous basis for simulations of the sulphation model.
- The fully discrete FTCS scheme, with time step $k(h)\to0$ chosen under the stability conditions (73)-(74), converges to the continuum solution in $L^\infty([0,T],H^1(\mathbb{R}_+))$ (Theorem 63).
- The convergence holds for every boundary process with the stated Hölder regularity, including the Pearson process (8) and fractional Brownian motion with Hurst index $>1/4$, so the result is not tied to one noise model.
- Time-space Besov spaces $B^k_{p,p}$ are the correct convergence spaces at this low boundary regularity: classical $L^2(0,T;H^1)$ convergence would generally fail because the solution's spatial regularity is tied to the Hölder exponent of the boundary path.
- The a priori estimates are uniform in $h$ for $0<h\le1$, which is what allows the extension operator and Besov compactness to transfer the discrete estimates to the continuum limit.
Reading between the lines
- The circular step in Proposition 43 (taking $f=s$ against condition (52)) is a gap in the written proof: a complete argument would need an intermediate bootstrap showing the local-in-time solution of Theorem 46 exists on $[0,T]$ with a norm bound uniform in $h$ before applying the linearized estimate. This may be fixable, but as written the uniform bound that Theorem 56 relies on is not fully esta
- The same splitting plus discrete-heat-kernel regularization in Besov spaces could plausibly handle other boundary-driven singular SPDEs on lattices, such as equations with interior multiplicative noise or with boundary noise of lower Hölder regularity; the paper's Lemma 32 already shows how far the heat kernel can compensate.
- Remark 62 notes the $O(k)$ rate for fixed $h$ is not uniform in $h$; tracking how the Lipschitz constants of the discrete nonlinearity grow as $h\to0$ could yield an explicit rate in $h$, which the paper does not provide.
- A cleaner convergence proof might avoid the $f=s$ bootstrap entirely by proving the uniform bound directly from a maximum principle or energy estimate for the nonlinear system (28), which would also strengthen Theorem 46 from local to global well-posedness in one step.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper studies a finite-difference semidiscretization of the random reaction-diffusion-ODE system (1)-(6) on the half-line, with a Holder-continuous stochastic Dirichlet boundary condition psi of regularity beta in (1/4,1/2). The authors split s=u+v into a random discrete heat equation with boundary psi and a nonlinear nonlocal discrete system for v with zero boundary, derive Besov-type a priori estimates for u and v, and then prove that piecewise-constant extensions of the discrete solutions converge in L^p([0,T],B^k_{p,p}(R+)) to the unique continuum mild/weak solution. A fully discrete scheme is also claimed to converge. The main result is Theorem 56; the proof relies on a uniform a priori bound (Proposition 43), compact embeddings of discrete Besov spaces, and a weak formulation of the discrete system.
Significance. If valid, Theorem 56 would be a meaningful first convergence theorem for this class of strongly nonlinear PDE-ODE systems with stochastic dynamical boundary conditions, and the use of discrete Besov spaces is well matched to the low boundary regularity. The paper's explicit discrete heat-kernel representation, the Feynman-Kac formula in Appendix A, and the reliance on the heat-kernel regularization from [8] are concrete strengths. However, the central uniform-in-h estimate is not closed as written, and the final strong-convergence arguments for both the semidiscrete and fully discrete schemes have gaps; these issues must be repaired before the theorem can be accepted.
major comments (3)
- [§6, Proposition 43] The proof of the uniform bound (65) is circular as written. Proposition 42 is stated for any f satisfying (52), and the third condition in (52) is exactly ||f||^2_{C([0,T],L^2(Omega_h^+))}+||D_h^+ f||^2_{L^2 L^2}<=K. Taking f=s, this is the left-hand side of (65), the estimate Proposition 43 is supposed to prove. Moreover, Theorem 46 invokes Proposition 43 to extend the local solution to time T_fin, so there is no independent global existence result to supply (52). A standard bootstrap repair would consist in first obtaining the structural conditions in (52) for the local nonlinear solution (positivity, s(.,0)=psi, ||s||_{L^infty}<=eta), then applying (63) with f=s and absorbing the term (1/2)||D_h^+ s||^2 to obtain L<=2mu; the paper does not provide these steps.
- [§7.2, proof of Theorem 56] The convergence proof only establishes weak convergence of E_h s_h in L^2([0,T],H^1) (and weak convergence of its difference quotients), whereas the theorem asserts strong convergence in L^p([0,T],B^k_{p,p}(R+)) for p in [1,2] and k<1/p. No uniform bound on the time derivative partial_t s_h (for instance in L^2_t H^{-1} or in a time-Besov space) is given, so the compactness needed to upgrade the weak convergence to the claimed strong convergence is absent. Strong convergence of c_h is proved separately, but it does not imply the missing time compactness for s_h; this gap is load-bearing because the final conclusion of Theorem 56 is exactly the strong convergence in L^p_t B^k_{p,p}.
- [§7.3.2, Theorem 63] As stated, Theorem 63 cannot hold: the fully discrete interpolant is piecewise constant in space, so it is not an H^1(R+) function, and Theorem 56 only provides convergence in L^p([0,T],B^k_{p,p}(R+)), not in L^infty([0,T],H^1(R+)). In the proof the first term ||s-s_h||_{L^infty H^1} is therefore not available. The theorem and its proof need to be reformulated in a norm compatible with the interpolant (for example L^p_t B^k_{p,p}), or the authors must introduce an H^1-conforming interpolation together with the required estimates.
minor comments (4)
- [Proposition 34] The condition 'r<2beta+p' appears to be a typo; the proof bounds u in W^{p,r} whenever r<alpha+k, with alpha<1/p-1 and k constrained by (49). The printed range should be corrected and the notation W^{p,r}_x made explicit (integrability p, differentiability r).
- [Theorem 56 proof] In the proof, 'weakly in L^2([0,T] x Omega_h^+)' should read 'weakly in L^2([0,T] x R+)' after extension by E_h; as printed the limit space is the discrete lattice.
- [Theorem 63] Theorem 63 contains two occurrences of 'lim_{h->+infty}' that should be 'lim_{h->0}', and the statement opens with 'Le us suppose' instead of 'Let us suppose'.
- [Lemma 30] Lemma 30 is imported from [8] without proof; since the heat-kernel regularization estimate is used repeatedly (Lemmas 31-32, Proposition 34), stating the precise Lemma 2.10 assumption and its proof or a self-contained derivation would improve verifiability.
Circularity Check
Proposition 43's uniform bound is derived by assuming the very norm bound it is meant to prove, making the convergence theorem's key compactness input circular as written.
-
self definitional
[Section 6, Proposition 43, proof of inequality (65), via Proposition 42, conditions (52) and inequality (64)]
"By Proposition 42 one can easily derive that the same estimate is true for the nonlinear system, taking f=s and applying inequality (64)."
Proposition 42 is stated only for f satisfying (52), whose first condition is ||f||^2_{C([0,T],L^2)} + ||D_h^+ f||^2_{L^2 L^2} ≤ K. Taking f=s, this hypothesis is exactly the left-hand side of the desired bound (65), i.e. sup_t ||s||^2_{L^2} + ∫ ||D_h^+ s||^2_{L^2} dt ≤ K. Thus the proof of Proposition 43 assumes (65) in order to apply Proposition 42 and conclude (65). Inequality (64) is not an independent estimate: it is obtained from (63) only after choosing K ≥ 2μ and after imposing the K-bound in (52). No separate bootstrap or continuation argument is written; Theorem 46 invokes Proposition 43 to extend local existence to T_fin, and Theorem 56 invokes Proposition 43 to obtain the uniform-in-h compactness bound.
full rationale
The central convergence proof is a derivation rather than a fit, so most of the paper is not circular. The well-posedness of the continuum system is imported from [35], an overlapping-author paper, but it is an independent existence/uniqueness theorem with its own hypotheses, and it is not derived from the discrete scheme; this self-citation is not itself load-bearing circularity. The heat-kernel Besov estimates cited from [8] are also overlapping-author results but are stated as independent lemmas used to regularize the boundary data. The genuinely circular step is Proposition 43: its proof sets f=s in Proposition 42, but Proposition 42's hypothesis (52) includes the norm bound that Proposition 43 is supposed to establish. Consequently the uniform estimate (65) used by Theorem 46 and Theorem 56 is not obtained from the paper's written arguments. Because this gap is load-bearing but localized and plausibly repairable by a bootstrap from the local existence theorem, the overall circularity score is 6 rather than higher.
Assumptions & free parameters
assumptions (6)
- domain assumption Existence, uniqueness, and regularity of the continuum mild solution to (9)-(10) with boundary (8), as stated in Theorem 3 of [35].
- domain assumption The boundary process psi is Hölder continuous with exponent beta in (1/4,1/2), bounded in [0,eta], with psi(0)=0, condition (7).
- domain assumption Porosity bound 0 < phi_min <= phi <= phi_max < 1, with B > 0 requiring 1/eta > B, Assumption 37 and condition (5).
- standard math Heat kernel regularization estimate on the positive lattice, Lemma 30, citing [8, Lemma 2.10].
- standard math Well-posedness of the martingale problem for the discrete generator with time-dependent coefficients, Lemma 68, citing [21, Theorem 7.3].
- standard math Equivalence of the weak formulation (67)-(68) with the mild formulation of Definition 2, via [7] and the mollifier argument in Lemma 55.
Cite this review
Pith. "Pith review of Discrete reaction-diffusion system with stochastic dynamical boundary conditions: convergence results." pith.science (2026). https://pith.science/paper/IEDUBBIA
@misc{pith2026250709278,
author = {Pith},
title = {Pith review of: Discrete reaction-diffusion system with stochastic dynamical boundary conditions: convergence results},
year = {2026},
howpublished = {\url{https://pith.science/paper/IEDUBBIA}},
note = {Machine review of arXiv:2507.09278}
}
read the original abstract
A space discrete approximation to a highly nonlinear reaction-diffusion system endowed with a stochastic dynamical boundary condition is analyzed and the convergence of the discrete scheme to the solution to the corresponding continuum random system is established. A splitting strategy allows us to decompose the random system into a space-discrete heat equation with a stochastic boundary condition, and a nonlinear and nonlocal space-discrete differential system coupled with the first one and with deterministic initial and boundary conditions. The convergence result is obtained by first establishing some a priori estimates for both space-discrete splitted variables and then exploiting compact embedding theorems for time-space Besov spaces on the positive lattice. The convergence of a fully discrete approximation of the random system is also discussed.
Forward citations
Cited by 1 Pith paper
-
Killed path-dependent McKean-Vlasov SDEs for a probabilistic representation of non-conservative McKean PDEs
A killed McKean-Vlasov SDE is shown to be well posed, and the density of its sub-probability law is a weak solution of the associated non-conservative path-dependent McKean PDE.
Reference graph
Works this paper leans on
-
[35]
M. Maurelli, D. Morale, and S. Ugolini. Well-posedness of a reaction-diffusion model with stochastic dynamical boundary conditions.Stochastic Processes and their Applications, 2025
work page 2025
- [4]
-
[8]
Elliptic stochastic quantization of Sinh-Gordon QFT
N. Barashkov and F. C. De Vecchi. Elliptic stochastic quantization of sinh-gordon qft.arXiv preprint arXiv:2108.12664, 2021
work page Pith review arXiv 2021
- [1]
-
[2]
F. Arceci, L.M. Giordano, M. Maurelli, D. Morale, and S. Ugolini. Randomness in a nonlinear model of sulphation phenomena. In G. Bretti, C. Cavaterra, M. Solci, and M. Spagnuolo, editors,Mathematical Modeling in Cultural Heritage, pages 31–50, Singapore, 2023. Springer Nature Singapore
work page 2023
- [3]
-
[5]
D. Aregba-Driollet, F. Diele, and R. Natalini. A mathematical model for the sulphur dioxide aggression to calcium carbonate stones: Numerical approximation and asymptotic analysis. SIAM Journal of Applied Mathematics, 64:1636–1667, 01 2004
work page 2004
-
[6]
Hajer Bahouri, Jean-Yves Chemin, and Raphael Danchin.Fourier Analysis and Nonlin- ear Partial Differential Equations, volume 343 ofGrundlehren der mathematischen Wis- senschaften. Springer, Berlin, 2011
work page 2011
Show all 41 references
-
[7]
Strongly continuous semigroups, weak solutions, and the variation of constants formula.Proceedings of the American Mathematical Society, 63:370–373, 04 1977
John Ball. Strongly continuous semigroups, weak solutions, and the variation of constants formula.Proceedings of the American Mathematical Society, 63:370–373, 04 1977
1977
-
[9]
Barbu, S
V. Barbu, S. Bonaccorsi, and L. Tubaro. Stochastic parabolic equations with nonlinear dy- namical boundary conditions.Journal of Mathematical Analysis and Applications, 427, 07 2015
2015
-
[10]
Bonaccorsi and G
S. Bonaccorsi and G. Ziglio. A semigroup approach to dynamic boundary value problems. Semigroup Forum, 69:159–183, 08 2004. 32
2004
-
[11]
Bonetti, C
E. Bonetti, C. Cavaterra, F. Freddi, M. Grasselli, and R. Natalini. A nonlinear model for marble sulphation including surface rugosity and mechanical damage.Nonlinear Analysis: Real World Applications, 73:103886, 2023
2023
-
[12]
Bonetti, C
E. Bonetti, C. Cavaterra, F. Freddi, M. Grasselli, and Natalini R. A nonlinear model for marble sulphation including surface rugosity: theoretical and numerical results.Communications on Pure & Applied Analysis, 18(2):977–998, 2019
2019
-
[13]
Brze´ zniak, B
Z. Brze´ zniak, B. Goldys, S. Peszat, and F. Russo. Second order pdes with dirichlet white noise boundary condition.Journal of Evolution Equations, 15, 05 2013
2013
-
[14]
Chalupecky and A
V. Chalupecky and A. Muntean. Semi-discrete finite difference multiscale scheme for a concrete corrosion model: a priori estimates and convergence.Japan Journal of Industral and Applied Mathematics, 29(2):289–316, 2012
2012
-
[15]
L Chen and S. Gan. Strong convergence and stationary distribution of an explicit scheme for the wright-fisher model.Journal of Computational and Applied Mathematics, 424, 12 2022
2022
-
[16]
First order strong convergence of an explicit scheme for the stochastic sis epidemic model.Journal of Computational and Applied Mathematics, 392, 02 2021
Lin Chen, Siqing Gan, and Xiaojie Wang. First order strong convergence of an explicit scheme for the stochastic sis epidemic model.Journal of Computational and Applied Mathematics, 392, 02 2021
2021
-
[17]
Chueshov and B
I. Chueshov and B. Schmalfuss. Parabolic stochastic partial differential equations with dy- namical boundary conditions.Differential and Integral Equations, 17, 07 2004
2004
-
[18]
Da Prato and J
G. Da Prato and J. Zabczyk. Evolution equations with white-noise boundary conditions. Stochastics and Stochastics Reports, 42:167–182, 1993
1993
-
[19]
Debussche, M
A. Debussche, M. Fuhrman, and G. Tessitore. Optimal control of a stochastic heat equation with boundary-noise and boundary-control.ESAIM: COCV, 13, 02 2007
2007
-
[20]
Discretisation of regularity structures.Ann
Dirk Erhard and Martin Hairer. Discretisation of regularity structures.Ann. Inst. Henri Poincar´ e Probab. Stat., 55(4):2209–2248, 2019
2019
-
[21]
Ethier and Thomas G
Stewart N. Ethier and Thomas G. Kurtz.Markov processes: Characterization and conver- gence. Wiley Series in Probability and Mathematical Statistics: Probability and Mathematical Statistics. John Wiley & Sons, Inc., New York, 1986
1986
-
[22]
F. R. Guarguaglini and R. Natalini. Global existence of solutions to a nonlinear model of sulphation phenomena in calcium carbonate stones.Nonlinear Analysis: Real World Applica- tions, 6(3):477–494, 2005
2005
-
[23]
F. R. Guarguaglini and R. Natalini. Fast reaction limit and large time behavior of solutions to a nonlinear model of sulphation phenomena.Communications in Partial Differential Equations, 32, 02 2007
2007
-
[24]
Gubinelli and M
M. Gubinelli and M. Hofmanov´ a. A PDE construction of the Euclideanϕ 4 3 quantum field theory.Comm. Math. Phys., 384(1):1–75, 2021
2021
-
[25]
Hairer and K
M. Hairer and K. Matetski. Discretisations of rough stochastic PDEs.Ann. Probab., 46(3):1651–1709, 2018
2018
-
[26]
The Φ 4 3 measure has sub-Gaussian tails.J
Martin Hairer and Rhys Steele. The Φ 4 3 measure has sub-Gaussian tails.J. Stat. Phys., 186(3):Paper No. 38, 25, 2022
2022
-
[27]
Hutzenthaler, A
M. Hutzenthaler, A. Jentzen, and P. Kloeden. Strong convergence of an explicit numerical method for sdes with nonglobally lipschitz continuous coefficients.The Annals of Applied Probability, 22, 10 2010
2010
-
[28]
J¨ averg˚ ard, D
J. J¨ averg˚ ard, D. Morale, G. Muntean, A. Rui, and S. Ugolini. A hybrid model of sulphation reactions: stochastic particles in a random continuum environment, 02 2025
2025
-
[29]
Jordan.Calculus of Finite Differences
C. Jordan.Calculus of Finite Differences. Chelsea Publishing Company, New York, NY, 1960. 33
1960
-
[30]
Jovanovi´ c and E
B. Jovanovi´ c and E. S¨ uli.Analysis of Finite Difference Schemes: For Linear Partial Differ- ential Equations with Generalized Solutions, volume 46. Springer, 01 2014
2014
-
[31]
P. E. Kloeden and E. Platen.Numerical Solution of Stochastic Differential Equations, vol- ume 23 ofStochastic Modelling and Applied Probability. Springer-Verlag, 1992
1992
-
[32]
Ladyˇ zenskaja, V.A
O.A. Ladyˇ zenskaja, V.A. Solonnikov, and N.N. Ural’ceva.Linear and Quasi-linear Equa- tions of Parabolic Type, volume 23 ofTranslations of Mathematical Monographs. American Mathematical Society, Providence, RI, 1968
1968
-
[33]
J. L. Lions.Quelques m´ ethodes de r´ esolution des probl` emes aux limites non lin´ eaires. Gauthier- Villars, Paris, 1969
1969
-
[34]
Martin and N
J. Martin and N. Perkowski. Paracontrolled distributions on bravais lattices and weak univer- sality of the 2d parabolic anderson model.Ann. Inst. H. Poincar´ e Probabilit´ es et Statistiques, 55(4):2058–2110, 2019
2019
-
[36]
Morale, G
D. Morale, G. Rui, and S. Ugolini. A stochastic interacting particle model for the marble sulphation process. InMathematical Modeling in Cultural Heritage, Singapore, 2025. Springer Nature Singapore
2025
-
[37]
Morale, L
D. Morale, L. Tarquini, and S. Ugolini. A probabilistic interpretation of a non-conservative and path-dependent nonlinear reaction-diffusion system for the marble sulphation in cultural heritage, 07 2024
2024
-
[38]
Springer, 01 2007
Alfio Quarteroni, Riccardo Sacco, and Fausto Saleri.Numerical Mathematics, volume 37. Springer, 01 2007
2007
-
[39]
L. C. G. Rogers and D. Williams.Diffusions, Markov Processes and Martingales, volume 2. Cambridge University Press, 2000
2000
-
[40]
Modern Birkh¨ auser Classics
Hans Triebel.Theory of Function Spaces. Modern Birkh¨ auser Classics. Birkh¨ auser Basel, 1 edition, 1983
1983
-
[41]
Turra.Infinitesimal analysis of singular stochastic partial differential equations
M. Turra.Infinitesimal analysis of singular stochastic partial differential equations. PhD thesis, University of Bonn, 2022. 34
2022
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.