pith. sign in

arxiv: 2202.09055 · v1 · submitted 2022-02-18 · 🧮 math.NA · cs.NA

Convergence analysis of a finite difference method for stochastic Cahn--Hilliard equation

Jialin Hong , Diancong Jin , Derui Sheng This is my paper

Pith reviewed 2026-05-24 12:25 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords stochastic Cahn-Hilliard equationfinite difference methodstrong convergenceMalliavin derivativemultiplicative noiseexponential Euler methoddensity convergencenegative moments
0
0 comments X

The pith

The finite difference method for the stochastic Cahn-Hilliard equation achieves strong spatial convergence of order 1 and temporal convergence of order nearly 3/8.

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

The paper proves that a spatial finite difference discretization of the stochastic Cahn-Hilliard equation with Lipschitz nonlinearity and multiplicative noise converges strongly with order 1. It establishes this rate for both the semi-discrete solution and its Malliavin derivative by means of estimates on the discrete Green function. Negative-moment bounds on the exact solution then yield L1 convergence of the numerical densities to the true density. An exponential Euler time discretization of the semi-discrete system produces strong convergence of order nearly 3/8 once Hölder continuity of the spatial approximation is controlled.

Core claim

Based on fine estimates of the discrete Green function, both the spatial semi-discrete numerical solution and its Malliavin derivative have strong convergence order 1. By showing the negative moment estimates of the exact solution, the density of the spatial semi-discrete numerical solution converges in L1(R) to the exact one. An exponential Euler method applied to the spatial semi-discrete solution yields temporal strong convergence order nearly 3/8 after the optimal Hölder continuity of that solution is derived.

What carries the argument

Fine estimates of the discrete Green function that produce uniform bounds on the Malliavin derivative of the semi-discrete solution.

If this is right

  • The spatial semi-discrete solution converges strongly with order 1 to the exact solution.
  • The Malliavin derivative of the spatial semi-discrete solution also converges strongly with order 1.
  • The probability density of the spatial semi-discrete solution converges to the exact density in L1(R).
  • The fully discrete scheme obtained by exponential Euler time stepping converges strongly with order nearly 3/8.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same Green-function technique could be tested on other stochastic fourth-order equations that admit a similar mild formulation.
  • Density convergence in L1 supplies a route to error bounds on expectations of bounded continuous functionals of the solution.
  • Improved temporal schemes might raise the 3/8 rate if they exploit the Hölder regularity already established for the spatial approximation.

Load-bearing premise

The nonlinearity is Lipschitz continuous and the exact solution satisfies negative moment bounds.

What would settle it

A sequence of numerical experiments with successively halved spatial mesh sizes in which the measured strong error fails to decrease proportionally to the mesh size would falsify the order-1 spatial claim.

read the original abstract

This paper presents the convergence analysis of the spatial finite difference method (FDM) for the stochastic Cahn--Hilliard equation with Lipschitz nonlinearity and multiplicative noise. Based on fine estimates of the discrete Green function, we prove that both the spatial semi-discrete numerical solution and its Malliavin derivative have strong convergence order $1$. Further, by showing the negative moment estimates of the exact solution, we obtain that the density of the spatial semi-discrete numerical solution converges in $L^1(\mathbb R)$ to the exact one. Finally, we apply an exponential Euler method to discretize the spatial semi-discrete numerical solution in time and show that the temporal strong convergence order is nearly $\frac38$, where a difficulty we overcome is to derive the optimal H\"older continuity of the spatial semi-discrete numerical solution.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The paper analyzes convergence of a finite-difference spatial semi-discretization for the stochastic Cahn-Hilliard equation with Lipschitz nonlinearity and multiplicative noise. Using discrete Green-function estimates, it proves strong order-1 convergence of both the semi-discrete solution and its Malliavin derivative. Negative-moment bounds on the exact solution yield L¹(ℝ) convergence of the numerical densities. An exponential Euler time discretization is then applied to the semi-discrete solution; the resulting temporal strong convergence rate is shown to be nearly 3/8 after establishing optimal Hölder continuity of the semi-discrete solution.

Significance. If the stated rates hold under the given hypotheses, the work supplies rigorous, explicit error bounds for a standard spatial discretization of an important SPDE, together with density convergence and a temporal rate obtained from Hölder regularity. The explicit construction of the discrete Green-function bounds, the negative-moment estimates, and the Hölder continuity argument are all strengths that make the analysis self-contained and potentially useful for related phase-field models with multiplicative noise.

minor comments (2)
  1. [Abstract, §1] Abstract and §1: the temporal rate is stated only as 'nearly 3/8'. The main theorem should record the precise form (e.g., 3/8−ε for arbitrary ε>0) together with the dependence on the Hölder exponent derived in the preceding section.
  2. [§2 (assumptions)] The statement that the nonlinearity is Lipschitz is used globally; a brief remark on whether the constant enters the final constants in the error bounds would clarify the dependence on the data.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the positive assessment of its contributions. We are pleased that the referee recommends minor revision. As the report contains no specific major comments, we have prepared no point-by-point revisions below.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper conducts a standard mathematical convergence analysis for a finite-difference discretization of the stochastic Cahn-Hilliard equation. All load-bearing ingredients (Lipschitz nonlinearity, discrete Green-function estimates for the biharmonic operator, negative-moment bounds, and Hölder continuity of the semi-discrete solution) are explicitly constructed or invoked under stated hypotheses within dedicated sections of the manuscript. No derivation step reduces by construction to a fitted parameter, self-referential definition, or load-bearing self-citation chain; the claimed spatial order-1 and temporal nearly-3/8 rates follow directly from the derived a-priori estimates rather than from renaming or smuggling prior results. The argument is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper relies on standard mathematical estimates and domain assumptions for stochastic PDEs rather than introducing new free parameters or entities.

axioms (2)
  • domain assumption Nonlinearity is Lipschitz continuous
    Explicitly stated in the abstract as a standing assumption on the equation.
  • standard math Discrete Green function admits fine estimates sufficient for order-1 convergence
    Invoked as the basis for the spatial convergence proof.

pith-pipeline@v0.9.0 · 5674 in / 1332 out tokens · 26313 ms · 2026-05-24T12:25:18.589843+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

25 extracted references · 25 canonical work pages

  1. [1]

    A fully discrete approximation of the one- dimensional stochastic heat equation

    Rikard Anton, David Cohen, and Lluis Quer-Sardanyons. A fully discrete approximation of the one- dimensional stochastic heat equation. IMA J. Numer. Anal. , 40(1):247–284, 2020

    work page 2020

  2. [2]

    Malliavin calculus for w hite noise driven parabolic SPDEs

    Vlad Bally and Etienne Pardoux. Malliavin calculus for w hite noise driven parabolic SPDEs. Potential Anal., 9(1):27–64, 1998

    work page 1998

  3. [3]

    The law of the Euler scheme for stochastic differential equations

    Vlad Bally and Denis Talay. The law of the Euler scheme for stochastic differential equations. II. Conver- gence rate of the density. Monte Carlo Methods Appl. , 2(2):93–128, 1996

    work page 1996

  4. [4]

    Cahn-Hilliard stochastic equa tion: existence of the solution and of its density

    Caroline Cardon-Weber. Cahn-Hilliard stochastic equa tion: existence of the solution and of its density. Bernoulli, 7(5):777–816, 2001

    work page 2001

  5. [5]

    Co nforming finite element methods for the stochastic Cahn-Hilliard-Cook equation

    Shimin Chai, Yanzhao Cao, Yongkui Zou, and Wenju Zhao. Co nforming finite element methods for the stochastic Cahn-Hilliard-Cook equation. Appl. Numer. Math. , 124:44–56, 2018

    work page 2018

  6. [6]

    Co nvergence of Density Approximations for Stochastic Heat Equation

    Chuchu Chen, Jianbo Cui, Jialin Hong, and Derui Sheng. Co nvergence of Density Approximations for Stochastic Heat Equation. arXiv:2007.12960

    work page arXiv 2007

  7. [7]

    A fully discrete approximation of the one-dimensional stochastic wave equation

    David Cohen and Llu ´ ıs Quer-Sardanyons. A fully discrete approximation of the one-dimensional stochastic wave equation. IMA J. Numer. Anal. , 36(1):400–420, 2016

    work page 2016

  8. [8]

    Absolute continuity and nume rical approximation of stochastic Cahn-Hilliard equation with unbounded noise diffusion

    Jianbo Cui and Jialin Hong. Absolute continuity and nume rical approximation of stochastic Cahn-Hilliard equation with unbounded noise diffusion. J. Differential Equations , 269(11):10143–10180, 2020

    work page 2020

  9. [9]

    Convergence in Density of Splitting A VF Scheme for Stochastic Langevin Equation

    Jianbo Cui, Jialin Hong, and Derui Sheng. Convergence in Density of Splitting A VF Scheme for Stochastic Langevin Equation. arXiv:1906.03439

    work page arXiv 1906

  10. [10]

    Strong converg ence of full discretization for stochastic Cahn- Hilliard equation driven by additive noise

    Jianbo Cui, Jialin Hong, and Liying Sun. Strong converg ence of full discretization for stochastic Cahn- Hilliard equation driven by additive noise. SIAM J. Numer. Anal. , 59(6):2866–2899, 2021. CONVERGENCE ANALYSIS OF FDM 23

    work page 2021

  11. [11]

    A. M. Davie and J. G. Gaines. Convergence of numerical sc hemes for the solution of parabolic stochastic partial differential equations. Math. Comp. , 70(233):121–134, 2001

    work page 2001

  12. [12]

    Qiang Du and R. A. Nicolaides. Numerical analysis of a co ntinuum model of phase transition. SIAM J. Numer. Anal. , 28(5):1310–1322, 1991

    work page 1991

  13. [13]

    Elliott and Stig Larsson

    Charles M. Elliott and Stig Larsson. Error estimates wi th smooth and nonsmooth data for a finite element method for the Cahn-Hilliard equation. Math. Comp. , 58(198):603–630, S33–S36, 1992

    work page 1992

  14. [14]

    Strong convergence of a fully discrete finite element approximation of the stochastic Cah n-Hilliard equation

    Daisuke Furihata, Mih´ aly Kov´ acs, Stig Larsson, and F redrik Lindgren. Strong convergence of a fully discrete finite element approximation of the stochastic Cah n-Hilliard equation. SIAM J. Numer. Anal. , 56(2):708–731, 2018

    work page 2018

  15. [15]

    Lattice approximations for stochastic quasi-linear parabolic partial differential equations driven by space-time white noise

    Istv´ an Gy¨ ongy. Lattice approximations for stochastic quasi-linear parabolic partial differential equations driven by space-time white noise. I. Potential Anal. , 9(1):1–25, 1998

    work page 1998

  16. [16]

    Hu and S

    Y. Hu and S. Watanabe. Donsker’s delta functions and app roximation of heat kernels by the time dis- cretization methods. J. Math. Kyoto Univ. , 36(3):499–518, 1996

    work page 1996

  17. [17]

    Analysis of Stochastic Partial Differential Equations , volume 119 of CBMS Regional Conference Series in Mathematics

    Davar Khoshnevisan. Analysis of Stochastic Partial Differential Equations , volume 119 of CBMS Regional Conference Series in Mathematics . Published for the Conference Board of the Mathematical Sci ences, Washington, DC; by the American Mathematical Society, Prov idence, RI, 2014

    work page 2014

  18. [18]

    Kohatsu-Higa

    A. Kohatsu-Higa. High order Itˆ o-Taylor approximations to heat kernels. J. Math. Kyoto Univ. , 37(1):129– 150, 1997

    work page 1997

  19. [19]

    Finite element approximation of the Cahn-Hilliard-Cook equation

    Mih´ aly Kov´ acs, Stig Larsson, and Ali Mesforush. Finite element approximation of the Cahn-Hilliard-Cook equation. SIAM J. Numer. Anal. , 49(6):2407–2429, 2011

    work page 2011

  20. [20]

    R. Kruse. Strong and Weak Approximation of Semilinear Stochastic Evo lution Equations , volume 2093 of Lecture Notes in Mathematics . Springer, Cham, 2014

    work page 2093

  21. [21]

    Finite-element approx imation of the linearized Cahn-Hilliard-Cook equa- tion

    Stig Larsson and Ali Mesforush. Finite-element approx imation of the linearized Cahn-Hilliard-Cook equa- tion. IMA J. Numer. Anal. , 31(4):1315–1333, 2011

    work page 2011

  22. [22]

    Convergence in total v ariation on Wiener chaos

    Ivan Nourdin and Guillaume Poly. Convergence in total v ariation on Wiener chaos. Stochastic Process. Appl., 123(2):651–674, 2013

    work page 2013

  23. [23]

    The Malliavin Calculus and Related Topics

    David Nualart. The Malliavin Calculus and Related Topics . Probability and its Applications (New York). Springer-Verlag, Berlin, second edition, 2006

    work page 2006

  24. [24]

    Error estimates of semidi screte and fully discrete finite element methods for the Cahn-Hilliard-Cook equation

    Ruisheng Qi and Xiaojie Wang. Error estimates of semidi screte and fully discrete finite element methods for the Cahn-Hilliard-Cook equation. SIAM J. Numer. Anal. , 58(3):1613–1653, 2020

    work page 2020

  25. [25]

    A stochastic wave equation in dimension 3: smoothness of the law

    Llu ´ ıs Quer-Sardanyons and Marta Sanz-Sol´ e. A stochastic wave equation in dimension 3: smoothness of the law. Bernoulli, 10(1):165–186, 2004. LSEC, ICMSEC, Academy of Mathematics and Systems Science, C hinese Academy of Sci- ences, and School of Mathematical Sciences, University of C hinese Academy of Sciences, Beijing 100049, China Email address : hj...

    work page 2004