Pith. sign in

REVIEW 1 major objections 2 minor 37 references

Numerical Approximation of the stochastic Cahn--Hilliard equation with singular potential

T0 review · 1 major / 2 minor · reviewed 2026-06-27 · grok-4.3

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

Pith's one-line read A regularized finite element scheme converges to the pathwise unique strong solution of the singular stochastic Cahn-Hilliard equation.

desk verdict A convergent regularized finite-element scheme for the singular stochastic Cahn-Hilliard equation, with the uniform H1 bound as the key unverified step. read the letter →

arxiv 2606.07164 v1 pith:BXOVAFHF submitted 2026-06-05 math.NA cs.NA

Ľubomír Baňas , Stefan Metzger This is my paper · ORCID
classification math.NAcs.NA
keywords stochasticCahn-Hilliardequationsingularpotentialfiniteelementapproximationnumericalconvergenceregularizationmultiplicativenoisemonotonicitymethod
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

The paper constructs a regularized fully discrete finite element scheme for the stochastic Cahn-Hilliard equation that includes a singular double-obstacle potential and multiplicative conservative noise. Stability estimates hold uniformly in the discretization parameters, allowing convergence first to a regularized equation via monotonicity arguments and then to the original singular equation via a uniform H^1 bound. This matters for obtaining reliable approximations when the singular potential prevents direct numerical treatment of the limiting problem.

What carries the argument

The regularised fully discrete finite element approximation scheme, whose uniform H^1-estimate independent of regularization and discretization parameters enables passage from the monotonicity limit to the singular limit.

What would settle it

Numerical evidence that the H^1 norm of solutions to the regularised problem grows unbounded as the regularization parameter tends to zero would disprove the passage to the singular limit.

Watch

Extended reading notes

Core claim

Thanks to a uniform H^1-estimate for the regularised problem the regularised solution converges towards the pathwise unique probabilistically strong solution of the original singular stochastic Cahn--Hilliard equation.

Load-bearing premise

The regularised problem admits a uniform H^1-estimate independent of the regularization and discretization parameters.

Editorial extensions

If this is right

  • The scheme produces stable approximations that remain controlled as mesh size and time step vanish.
  • Convergence holds pathwise to the unique probabilistically strong solution.
  • Numerical experiments can directly compare the regularised scheme against its unregularised counterpart and display the influence of the conservative noise.

Reading between the lines

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

  • The uniform stability technique could be tested on related singular stochastic phase-field models with different noise structures.
  • If the H^1 bound persists under weaker noise assumptions, the method might extend to equations with non-conservative multiplicative noise.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

1 major / 2 minor

Summary. The manuscript develops a regularised fully discrete finite element scheme for the stochastic Cahn-Hilliard equation with singular double-obstacle potential and multiplicative conservative noise. It establishes stability estimates uniform in the discretization parameters, proves convergence of the scheme to a regularised (still singular-potential) equation via monotonicity arguments, and then invokes a uniform H¹ estimate on the regularised solutions to pass to the limit as the regularization parameter vanishes, obtaining convergence to the pathwise-unique probabilistically strong solution of the original equation. Numerical simulations comparing the regularised and unregularised approximations are included.

Significance. If the uniform H¹ estimate is established independently of the regularization parameter, the result supplies a rigorous convergence theory for numerical approximation of a singular stochastic PDE that is otherwise difficult to treat directly; the combination of monotonicity-based convergence for the regularised problem with a uniform energy bound is a standard but technically demanding route in this area and would be a useful addition to the literature on numerical SPDEs.

major comments (1)
  1. [uniform H¹-estimate (abstract and the section establishing the singular limit)] The passage from the regularised to the original singular equation rests on a uniform H¹ estimate for the regularised problem that is independent of the regularization parameter ε (invoked after the monotonicity step). In the multiplicative-noise setting this bound is typically derived via Itô calculus applied to a regularised energy; the proof must explicitly control any ε-dependent contributions arising from the noise term and from the approximation of the singular potential, otherwise the limit may fail to satisfy the double-obstacle constraint in the required sense.
minor comments (2)
  1. The precise form of the regularization of the singular potential (e.g., the specific mollification or penalty function) should be stated explicitly in the abstract and in the statement of the regularised problem.
  2. Notation for the fully discrete scheme (time-step, spatial mesh size, regularization parameter) should be introduced once and used consistently throughout the stability and convergence statements.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive feedback. We address the single major comment below.

read point-by-point responses

  1. Referee: [uniform H¹-estimate (abstract and the section establishing the singular limit)] The passage from the regularised to the original singular equation rests on a uniform H¹ estimate for the regularised problem that is independent of the regularization parameter ε (invoked after the monotonicity step). In the multiplicative-noise setting this bound is typically derived via Itô calculus applied to a regularised energy; the proof must explicitly control any ε-dependent contributions arising from the noise term and from the approximation of the singular potential, otherwise the limit may fail to satisfy the double-obstacle constraint in the required sense.

    Authors: We appreciate this observation. The uniform H¹ estimate is derived in Section 4 by applying Itô's formula to the regularised energy. The multiplicative conservative noise term produces no ε-dependent blow-up because the noise is divergence-free and the regularisation preserves the zero-mean constraint; the resulting stochastic integral is controlled via Burkholder-Davis-Gundy and the uniform L² bound already obtained from the monotonicity argument. The approximation of the singular potential contributes only non-positive terms that vanish as ε→0 by construction of the regularisation. Consequently the double-obstacle constraint is recovered in the limit. We will add an explicit remark after the statement of the uniform estimate clarifying these controls. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: convergence proof relies on independent uniform estimates and monotonicity arguments.

full rationale

The derivation proceeds by first obtaining stability uniform in discretization parameters, then monotonicity-based convergence to a regularised equation, followed by passage to the singular limit via an invoked uniform H^1 bound. The target is an independently defined pathwise-unique strong solution; no step reduces by construction to a fitted parameter, self-referential definition, or load-bearing self-citation chain. The uniform estimate is presented as a separate analytical ingredient rather than being smuggled in via ansatz or renaming.

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

The paper rests on standard background results from stochastic PDE theory and numerical analysis; no free parameters or invented entities are introduced in the abstract.

assumptions (2)
  • domain assumption The regularised problem admits a uniform H^1-estimate independent of regularization and discretization parameters.
    Invoked to pass from the regularised limit to the singular equation.
  • domain assumption The singular stochastic Cahn-Hilliard equation possesses a pathwise unique probabilistically strong solution.
    The numerical scheme is shown to converge to this solution.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Numerical Approximation of the stochastic Cahn--Hilliard equation with singular potential." pith.science (2026). https://pith.science/paper/BXOVAFHF

@misc{pith2026260607164,
  author       = {Pith},
  title        = {Pith review of: Numerical Approximation of the stochastic Cahn--Hilliard equation with singular potential},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/BXOVAFHF}},
  note         = {Machine review of arXiv:2606.07164}
}
abstract

We discuss the numerical approximation of the stochastic Cahn--Hilliard equation with a singular double-obstacle potential and multiplicative conservative noise. We propose a regularised fully discrete finite element approximation scheme for the problem and show that it satisfies stability estimates which are uniform with respect to the discretization parameters. We show convergence of the approximation for vanishing discretization parameters towards a regularised version of the singular stochastic Cahn--Hilliard equation by monotonicity arguments. Hence, thanks to a uniform $H^1$-estimate for the regularised problem we show that the regularised solution converges towards the pathwise unique probabilistically strong solution of the original singular stochastic Cahn--Hilliard equation. We conclude by presenting numerical simulations where we compare the regularised numerical approximation to its unregularised counterpart and illustrate the effect of the conservative noise.

Figures

Figures reproduced from arXiv: 2606.07164 by the authors.

Figure 1
Figure 1. Initial condition (left), deterministic numerical solution (ν = 0) (middle) and stochastic numerical solution (ν = 1.6 −4 ) (right) at time t = 1.2 × 10−2 . of the square, i.e., dim Ve h = (L + 1)2 + L 2 and ψek and the nodal basis functions of Ve h. In the experiments below we choose L = 16. 8.1. Spinodal decomposition. We take T = 0.012, ε = 1/(12π), ν = 10−5 . The computations were performed with the time step si… view at source ↗
Figure 2
Figure 2. Deterministic numerical solution (ν = 0) at time t = 1.2, 1.5, 2, 2.5 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p046_2.png] view at source ↗
Figure 3
Figure 3. One path of the stochastic numerical solution (ν = 1.6 × 10−4 ) at time t = 0.6, 0.8, 0.9, 2.5 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p046_3.png] view at source ↗
Figures from the paper (8 more)
Figure 4
Figure 4. Figure 4: One path of the stochastic numerical solution (ν = 1.6 × 10−4 , halved squares, from lower left to upper right corner) at time t = 0.6, 0.8, 0.9, 2.5 × 10−3 . corresponding adaptive meshes are displayed in [PITH_FULL_IMAGE:figures/full_fig_p046_4.png]
Figure 5
Figure 5. Figure 5: One path of the stochastic numerical solution (ν = 1.6 × 10−4 , halved squares, from lower right to upper left corner) at time t = 0.6, 0.8, 0.9, 2.5 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p047_5.png]
Figure 6
Figure 6. Figure 6: Deterministic numerical solution (ν = 0) at time t = 0, 2.5, 6.5, 8 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p047_6.png]
Figure 7
Figure 7. Figure 7: Stochastic numerical solution (ν = 1.6) at time t = 0, 2.5, 6.5, 8 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p047_7.png]
Figure 8
Figure 8. Figure 8: Finite element mesh for one path of the stochastic numerical solution (ν = 1.6) at time t = 0, 2.5, 6.5, 8 × 10−3 [PITH_FULL_IMAGE:figures/full_fig_p047_8.png]
Figure 9
Figure 9. Figure 9: Evolution of the expected value of the energy (left) and of the area of the middle circle (right). 0 0.002 0.004 0.006 0.008 0.01 0.012 0.006 0.0065 0.007 0.0075 0.008 0.0085 0.009 nu=1.6 nu=1.6, delta=1e-3 nu=0 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.0065 0.007 0.0075 …
Figure 10
Figure 10. Figure 10: Histogram of closing times of the inner circle for ν = 1.6 (left) and for ν = 0.16, 0.8, 1.6 (right) along with the area of the inner circle for ν = 0. In the deterministic setting the inner circle disappears at the time t = 0.00674, the the stochastic setting the occ…
Figure 11
Figure 11. Figure 11: Evolution of the expected value of the energy (left) and of the area of the inner circle (right) for ν = 0, 0.16, 0.8, 1.6. 3. , Finite element approximation of a model for phase separation of a multi-component alloy with non-smooth free energy and a concentration dep…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

37 extracted references

  1. [1]

    Antonopoulou, Ľ

    D. Antonopoulou, Ľ. Baňas, R. Nürnberg, and A. Prohl,Numerical approximation of the stochastic Cahn-Hilliard equation near the sharp interface limit, Numer. Math.147(2021), no. 3, 505–551. MR 4224928

    2021

  2. [2]

    J. W. Barrett and J. F. Blowey,Finite element approximation of a model for phase separation of a multi-component alloy with non-smooth free energy, Numer. Math.77(1997), 1–34. APPROXIMATION OF THE NON-SMOOTH STOCHASTIC CAHN–HILLIARD EQUATION 49 5 5.2 5.4 5.6 5.8 6 6.2 0.0065 0.007 0.0075 0.008 0.0085 nu=1.6 nu=0.8 nu=0.16 nu=0 0 0.002 0.004 0.006 0.008 0....

    1997

  3. [3]

    ,Finite element approximation of a model for phase separation of a multi-component alloy with non-smooth free energy and a concentration dependent mobility matrix, M3AS9(1999), 627–663

    1999

  4. [4]

    J.W.Barrett, R.Nürnberg, andV.Styles,Finite element approximation of a phase field model for void electromigration, SIAM J. Numer. Anal.42(2004), no. 2, 738–772. MR MR2084234 (2005m:65206)

    2004

  5. [5]

    Bauzet, E

    C. Bauzet, E. Bonetti, G. Bonfanti, F. Lebon, and G. Vallet,A global existence and uniqueness result for a stochastic Allen–Cahn equation with constraint, Mathematical Methods in the Applied Sciences 40(2017), no. 14, 5241–5261

    2017

  6. [6]

    Bauzet, C

    C. Bauzet, C. Sultan, G. Vallet, and A. Zimmermann,Theoretical analysis of a finite-volume scheme for a stochastic Allen-Cahn problem with constraint, Nonlinear Anal.259(2025), Paper No. 113812, 21

    2025

  7. [7]

    Baňas, Ľ

    R. Baňas, Ľ. Nürnberg,Numerical approximation of a non-smooth phase-field model for multicompo- nent incompressible flow, ESAIM Math. Model. Numer. Anal.51(2017), no. 3, 1089–1117

    2017

  8. [8]

    Baňas, Z

    Ľ. Baňas, Z. Brzeźniak, M. Neklyudov, and A. Prohl,A convergent finite-element-based discretization of the stochastic landau–lifshitz–gilbert equation, IMA J. Numer. Anal.34(2013), no. 2, 502–549

    2013

Show all 37 references
  1. [9]

    Baňas and J

    Ľ. Baňas and J. D. Mukam,Numerical approximation of the stochastic cahn-hilliard equation with space-time white noise near the sharp interface limit, IMA J. Numer. Anal. (2025), draf064

    2025

  2. [10]

    Baňas, A

    Ľ. Baňas, A. Novick-Cohen, and R. Nürnberg,The degenerate and non-degenerate deep quench ob- stacle problem: A numerical comparison, Netw. Heterog. Media8(2013), no. 1, 37–64

    2013

  3. [11]

    Baňas and R

    Ľ. Baňas and R. Nürnberg,Adaptive finite element methods for Cahn-Hilliard equations, J. Comput. Appl. Math.218(2008), no. 1, 2–11

    2008

  4. [12]

    Baňas and R

    Ľ. Baňas and R. Nürnberg,Finite element approximation of a three dimensional phase field model for void electromigration, J. Sci. Comput.37(2008), no. 2, 202–232. MR 2453219

    2008

  5. [13]

    Baňas and R

    Ľ. Baňas and R. Nürnberg,A posteriori estimates for the Cahn-Hilliard equation with obstacle free energy, M2AN Math. Model. Numer. Anal.43(2009), no. 5, 1003–1026

    2009

  6. [14]

    Baňas and C

    Ľ. Baňas and C. Vieth,Robust a posteriori estimates for the stochastic Cahn-Hilliard equation, Math. Comp.92(2023), no. 343, 2025–2063. MR 4593209

    2023

  7. [15]

    J. F. Blowey and C. M. Elliott,The Cahn–Hilliard gradient theory for phase separation with non- smooth free energy. Part I: Mathematical analysis, European J. Appl. Math.2(1991), 233–279

    1991

  8. [16]

    Part II: Numerical analysis, European J

    ,The Cahn–Hilliard gradient theory for phase separation with non-smooth free energy. Part II: Numerical analysis, European J. Appl. Math.3(1992), 147–179

    1992

  9. [17]

    Bréhier and L

    C.-E. Bréhier and L. Goudenège,Weak convergence rates of splitting schemes for the stochastic Allen- Cahn equation, BIT60(2020), no. 3, 543–582

    2020

  10. [18]

    Breit and A

    D. Breit and A. Prohl,Weak error analysis for the stochastic Allen-Cahn equation, Stoch. Partial Differ. Equ. Anal. Comput.12(2024), no. 4, 2181–2245

    2024

  11. [19]

    S. C. Brenner and L. R. Scott,The mathematical theory of finite element methods (second edition), Springer-Verlag, New York, 2002. 50 Ľ. BAŇAS, S. METZGER

    2002

  12. [20]

    J. W. Cahn and J. E. Hilliard,Free energy of a non-uniform system. I. Interfacial free energy, J. Chem. Phys.28(1958), 258–267

    1958

  13. [21]

    Cook,Brownian motion in spinodal decomposition, Acta Metallurgica18(1970), no

    H.E. Cook,Brownian motion in spinodal decomposition, Acta Metallurgica18(1970), no. 3, 297–306

    1970

  14. [22]

    3, 1473– 1494

    Arnaud Debussche and Ludovic Goudenege,Stochastic Cahn-Hilliard Equation with Double Singular Nonlinearities and Two Reflections, Siam Journal on Mathematical Analysis43(2011), no. 3, 1473– 1494

    2011

  15. [23]

    Du and X

    Q. Du and X. Feng,The phase field method for geometric moving interfaces and their numeri- cal approximations, Geometric partial differential equations. Part I, Handb. Numer. Anal., vol. 21, Elsevier/North-Holland, Amsterdam, 2020, pp. 425–508. MR 4378430

    2020

  16. [24]

    Feireisl and A

    E. Feireisl and A. Novotný,Singular Limits in Thermodynamics of Viscous Fluids, Springer-Verlag GmbH, 2017

    2017

  17. [25]

    Flandoli and D

    F. Flandoli and D. Gatarek,Martingale and stationary solutions for stochastic Navier-Stokes equa- tions, Probab. Theory Related Fields102(1995), no. 3, 367–391

    1995

  18. [26]

    Furihata, M

    D. Furihata, M. Kovács, Larsson S., and F. Lindgren,Strong convergence of a fully discrete finite element approximation of the stochastic Cahn-Hilliard equation, SIAM J. Numer. Anal.56(2018), no. 2, 708–731

    2018

  19. [27]

    Giorgini, M

    A. Giorgini, M. Grasselli, and A. Miranville,The Cahn–Hilliard–Oono equation with singular poten- tial, Mathematical Models and Methods in Applied Sciences27(2017), no. 13, 2485–2510

    2017

  20. [28]

    Girault and P.-A

    V. Girault and P.-A. Raviart,Finite element methods for Navier–Stokes, Springer-Verlag, Berlin, 1986

    1986

  21. [29]

    Gräser and R

    C. Gräser and R. Kornhuber,Nonsmooth Newton methods for set-valued saddle point problems, SIAM J. Numer. Anal.47(2009), no. 2, 1251–1273. MR 2485452

    2009

  22. [30]

    Liu and M

    W. Liu and M. Röckner,Stochastic partial differential equations: an introduction, Springer, Cham; Heidelberg; New York; Dordrecht; London, 2015

    2015

  23. [31]

    Metzger,A convergent finite element scheme for a fourth-order liquid crystal model, IMA J

    S. Metzger,A convergent finite element scheme for a fourth-order liquid crystal model, IMA J. Nu- mer. Anal.42(2022), no. 1, 440–486, published online 2020

    2022

  24. [32]

    (2025), 64 pages

    S.Metzger,A convergent augmented SAV scheme for stochastic Cahn–Hilliard equations with dynamic boundary conditions describing contact line tension, Interfaces Free Bound. (2025), 64 pages

    2025

  25. [33]

    ,A convergent stochastic scalar auxiliary variable method, IMA J. Numer. Anal.45(2025), no. 4, 2404–2469, published online 2024

    2025

  26. [34]

    ,Strong error estimates for a fully discrete SAV scheme for the stochastic Allen–Cahn equation with multiplicative noise, ESAIM: Math. Model. Numer. Anal.59(2025), no. 4, 2305–2327

    2025

  27. [35]

    Di Primio, M

    A. Di Primio, M. Grasselli, and L. Scarpa,Stochastic Cahn–Hilliard and conserved Allen–Cahn equa- tions with logarithmic potential and conservative noise, Nonlinearity37(2024), no. 12, 125005

    2024

  28. [36]

    Scarpa,On the stochastic Cahn-Hilliard equation with a singular double-well potential, Nonlinear Anal.171(2018), 102–133

    L. Scarpa,On the stochastic Cahn-Hilliard equation with a singular double-well potential, Nonlinear Anal.171(2018), 102–133. MR 3778274

    2018

  29. [37]

    Simon,Sobolev, Besov and Nikolskii fractional spaces: Imbeddings and comparisons for vector valued spaces on an interval, Ann

    J. Simon,Sobolev, Besov and Nikolskii fractional spaces: Imbeddings and comparisons for vector valued spaces on an interval, Ann. di Mat. Pura ed Appl. (IV)157(1990), no. 1, 117–148

    1990

Pith tools

Reviewed June 27, 2026 · model on record in the stance chip above.