pith. sign in

arxiv: 2408.00953 · v2 · submitted 2024-08-01 · 🧮 math.NA · cs.NA

Approximation of the invariant measure for stochastic Allen-Cahn equation via an explicit fully discrete scheme

Yibo Wang , Wanrong Cao This is my paper

Pith reviewed 2026-05-23 22:14 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords stochastic Allen-Cahninvariant measureweak error analysisMalliavin calculusfully discrete schemespectral Galerkinexponential Euler
0
0 comments X

The pith

An explicit fully discrete scheme approximates the invariant measure of the stochastic Allen-Cahn equation.

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

The authors introduce an explicit fully discrete numerical scheme for the stochastic Allen-Cahn equation using spectral Galerkin in space and tamed accelerated exponential Euler in time. They establish time-independent boundedness of moments for the numerical solutions. This boundedness allows weak error analysis over infinite time intervals using Malliavin calculus. The result offers a method to numerically approximate the invariant measure of the equation.

Core claim

The paper claims that the proposed scheme yields numerical solutions whose moments are bounded independently of time, which in turn enables the application of Malliavin calculus to obtain weak error estimates on infinite time intervals and thereby approximate the invariant measure.

What carries the argument

Tamed accelerated exponential Euler scheme after spectral Galerkin discretization, whose moment bounds support infinite-time weak convergence analysis via Malliavin calculus.

If this is right

  • The weak error remains controlled as time tends to infinity.
  • The scheme can simulate the long-term behavior of the stochastic Allen-Cahn equation.
  • The approach provides convergence to the invariant measure under the discretization parameters.

Where Pith is reading between the lines

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

  • The same moment bound technique might apply to other stochastic evolution equations with similar dissipative properties.
  • Alternative error analysis methods could be explored if Malliavin calculus is not preferred.
  • Practical implementations could verify the moment boundedness numerically for specific parameters.

Load-bearing premise

The moments of the numerical solutions remain bounded independently of the time horizon.

What would settle it

Numerical experiments demonstrating that higher moments of the discrete solutions grow unboundedly with time would invalidate the infinite-time analysis.

read the original abstract

In this paper we propose an explicit fully discrete scheme to numerically solve the stochastic Allen-Cahn equation. The spatial discretization is done by a spectral Galerkin method, followed by the temporal discretization by a tamed accelerated exponential Euler scheme. Based on the time-independent boundedness of moments of numerical solutions, we present the weak error analysis in an infinite time interval by using Malliavin calculus. This provides a way to numerically approximate the invariant measure for the stochastic Allen-Cahn equation.

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

2 major / 1 minor

Summary. The manuscript proposes an explicit fully discrete scheme for the stochastic Allen-Cahn equation consisting of spectral Galerkin spatial discretization followed by a tamed accelerated exponential Euler time discretization. It establishes time-independent boundedness of moments for the numerical solutions and applies Malliavin calculus to obtain weak error estimates over the infinite time interval, thereby providing a numerical approximation to the invariant measure.

Significance. If the central claims hold, the work supplies a rigorous explicit scheme and infinite-horizon weak analysis via Malliavin calculus for approximating invariant measures of a nonlinear SPDE, which is a technically demanding task. The combination of taming for explicitness and Malliavin tools for long-time weak convergence is a substantive contribution to computational stochastic dynamics.

major comments (2)
  1. [moment bounds section] The section establishing time-independent moment bounds (presumably §3 or §4): these bounds are invoked to justify the Malliavin analysis on [0,∞), but must be shown to be uniform in the spectral cutoff N and time step Δt; non-uniformity would prevent passage to the limit in the integration-by-parts formula and undermine the invariant-measure approximation.
  2. [weak error analysis section] The weak error analysis on the infinite interval (presumably §5): the taming term must be controlled so that its effect on the invariant measure vanishes uniformly as N→∞ and Δt→0; otherwise the limiting object approximated by the scheme may differ from the true invariant measure of the continuous equation.
minor comments (1)
  1. [§2] Notation for the taming parameter and the accelerated exponential Euler scheme should be introduced with a clear reference to the continuous equation (1.1) to avoid ambiguity in the error analysis.

Simulated Author's Rebuttal

2 responses · 0 unresolved

Dear Editor, We thank the referee for the thorough review and valuable comments on our manuscript. We address each major comment point by point below, providing clarifications on the uniformity of estimates and the control of the taming term while remaining faithful to the content of the paper.

read point-by-point responses

  1. Referee: [moment bounds section] The section establishing time-independent moment bounds (presumably §3 or §4): these bounds are invoked to justify the Malliavin analysis on [0,∞), but must be shown to be uniform in the spectral cutoff N and time step Δt; non-uniformity would prevent passage to the limit in the integration-by-parts formula and undermine the invariant-measure approximation.

    Authors: We appreciate the referee's emphasis on uniformity. In Section 3, the time-independent moment bounds for the fully discrete solutions are derived using the one-sided Lipschitz property of the Allen-Cahn nonlinearity combined with the taming mechanism and the spectral projection. The resulting constants depend only on the equation coefficients, the noise intensity, and the domain, and are independent of both the Galerkin dimension N and the time step Δt. This independence follows directly from the a priori estimates that close without invoking any discretization-specific constants that grow with refinement. We will add an explicit remark after the main theorem in Section 3 to state this uniformity clearly, facilitating the subsequent passage to the limit in the Malliavin integration-by-parts formula. revision: partial

  2. Referee: [weak error analysis section] The weak error analysis on the infinite interval (presumably §5): the taming term must be controlled so that its effect on the invariant measure vanishes uniformly as N→∞ and Δt→0; otherwise the limiting object approximated by the scheme may differ from the true invariant measure of the continuous equation.

    Authors: We agree that uniform control of the taming term is essential. In Section 5 the weak error analysis proceeds by comparing the numerical invariant measure to the continuous one via Malliavin calculus on [0,∞). The taming is activated only on a set whose probability is controlled uniformly by the moment bounds of Section 3; consequently, the contribution of the taming correction to the weak error tends to zero as N→∞ and Δt→0, uniformly in time. This is obtained by splitting the expectation into the region where taming is inactive (where the scheme coincides with the untamed exponential Euler) and the complementary region (whose measure vanishes). We will strengthen the exposition in the revised Section 5 by adding a dedicated lemma that quantifies the uniform vanishing of the taming effect with respect to the discretization parameters. revision: yes

Circularity Check

0 steps flagged

No circularity; analysis rests on external Malliavin calculus and separate moment bounds

full rationale

The paper proposes an explicit fully discrete scheme (spectral Galerkin + tamed accelerated exponential Euler) and invokes time-independent boundedness of moments of the numerical solutions as the foundation for weak error analysis on [0,∞) via Malliavin calculus to approximate the invariant measure. No equations, fitted parameters, or self-citations are shown that reduce the claimed convergence or invariant-measure approximation to a definition or input by construction. The moment bounds are treated as an independent prerequisite rather than derived from the target result, and Malliavin calculus is an external tool. This is the normal case of a self-contained derivation without circular steps.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Review performed from abstract only; no explicit free parameters, axioms, or invented entities are stated in the provided text.

pith-pipeline@v0.9.0 · 5600 in / 1208 out tokens · 17855 ms · 2026-05-23T22:14:50.737796+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages

  1. [1]

    Br´ ehier, Approximation of the invariant distribution for a class of ergodic SPDEs using an explicit tamed exponential Euler scheme, ESAIM: M2AN, 56 (2022 ) 151-175

    C.E. Br´ ehier, Approximation of the invariant distribution for a class of ergodic SPDEs using an explicit tamed exponential Euler scheme, ESAIM: M2AN, 56 (2022 ) 151-175

    work page 2022

  2. [2]

    Br´ ehier, Approximation of the invariant measure with an Eule r scheme for stochastic PDEs driven by space-time white noise, Potential Anal., 40 (2014) 1– 40

    C.E. Br´ ehier, Approximation of the invariant measure with an Eule r scheme for stochastic PDEs driven by space-time white noise, Potential Anal., 40 (2014) 1– 40

    work page 2014

  3. [3]

    Br´ ehier, J.B

    C.E. Br´ ehier, J.B. Cui, J.L. Hong, Strong convergence rates of semi-discrete splitting ap- proximations for stochastic Allen–Cahn equation, IMA J. Numer. An al., 39 (2019) 2096- 2134

    work page 2019

  4. [4]

    Br´ ehier, M

    C.E. Br´ ehier, M. Kopec, Approximation of the invariant law of SPD Es: error analysis using a Poisson equation for a full-discretization scheme, IMA J. Num er. Anal., 37 (2017) 1375-1410

    work page 2017

  5. [5]

    Cai, S.Q

    M. Cai, S.Q. Gan, X.J. Wang, Weak convergence rates for an explic it full-discretization of stochastic Allen–Cahn equation with additive noise, J. Sci. Comput., 8 6 (2021) 34

    work page 2021

  6. [6]

    Cerrai, Second order Pde’s in finite and infinite dimension: a prob abilistic approach, Lecture Notes in Mathematics, Springer, 2001

    S. Cerrai, Second order Pde’s in finite and infinite dimension: a prob abilistic approach, Lecture Notes in Mathematics, Springer, 2001. 24

    work page 2001

  7. [7]

    Dang, J.L

    C.C.Chen, T.H. Dang, J.L. Hong, T. Zhou, CLT for approximating er godic limit of SPDEs via a full discretization, Stoch. Process. Appl., 157 (2023) 1-41

    work page 2023

  8. [8]

    Chen, J.L

    C.C. Chen, J.L. Hong, X. Wang, Approximation of invariant measur e for damped stochastic nonlinear Schr¨ odinger equation via an ergodic numerical scheme, Potential Anal., 46 (2017) 323-367

    work page 2017

  9. [9]

    Chen, S.Q

    Z.H. Chen, S.Q. Gan, X.J. Wang, A full-discrete exponential Euler a pproximation of the in- variant measure for parabolic stochastic partial differential equa tions, Appl. Numer. Math., 157 (2020) 135-158

    work page 2020

  10. [10]

    Cui, J.L

    J.B. Cui, J.L. Hong, L.Y. Sun, Weak convergence and invariant me asure of a full discretiza- tion for parabolic SPDEs with non-globally Lipschitz coefficients, Stoc h. Process. Appl., 134 (2021) 55-93

    work page 2021

  11. [11]

    Da Prato, An introduction to infinite-dimensional analysis, Sp ringer, 2006

    G. Da Prato, An introduction to infinite-dimensional analysis, Sp ringer, 2006

    work page 2006

  12. [12]

    J.L. Hong, X. Wang, Invariant measures for stochastic nonline ar Schr¨ odinger equations, numerical approximations and symplectic structures, Springer, 2 019

    work page

  13. [13]

    J.L. Hong, X. Wang, L.Y. Zhang, Numerical analysis on ergodic limit of approximations for stochastic NLS equation via multi-symplectic scheme, SIAM J. Nu mer. Anal., 55 (2017) 305-327

    work page 2017

  14. [14]

    Kruse, Optimal error estimates of Galerkin finite element met hods for stochastic partial differential equations with multiplicative noise, IMA J

    R. Kruse, Optimal error estimates of Galerkin finite element met hods for stochastic partial differential equations with multiplicative noise, IMA J. Numer. Anal., 34 (2014) 217-251

    work page 2014

  15. [15]

    Kruse, Strong and weak approximation of semilinear stochas tic evolution equations, Lecture Notes in Mathematics, Springer, 2014

    R. Kruse, Strong and weak approximation of semilinear stochas tic evolution equations, Lecture Notes in Mathematics, Springer, 2014

    work page 2014

  16. [16]

    Kruse, S

    R. Kruse, S. Larsson, Optimal regularity for semilinear stocha stic partial differential equa- tions with multiplicative noise, Electron. J. Probab., 17 (2012) 1-19

    work page 2012

  17. [17]

    Lin, R.S

    Q. Lin, R.S. Qi, Optimal weak order and approximation of the invar iant measure with a fully-discrete Euler scheme for semilinear stochastic parabolic equa tions with additive noise, Mathematics, 12 (2024) 112

    work page 2024

  18. [18]

    Liu, Z.H

    Z.H. Liu, Z.H. Qiao, Strong approximation of monotone stochast ic partial differential equa- tions driven by white noise, IMA J. Numer. Anal., 40 (2020) 1074-109 3

    work page 2020

  19. [19]

    Nualart, The Malliavin calculus and related topics, Probability an d its Applications, Springer, 2006

    D. Nualart, The Malliavin calculus and related topics, Probability an d its Applications, Springer, 2006

    work page 2006

  20. [20]

    Rothe, Global solutions of reaction-diffusion systems, Lect ure Notes in Mathematics, Springer, 1984

    F. Rothe, Global solutions of reaction-diffusion systems, Lect ure Notes in Mathematics, Springer, 1984

    work page 1984

  21. [21]

    Wang, An efficient explicit full-discrete scheme for strong ap proximation of stochastic Allen–Cahn equation, Stoch

    X.J. Wang, An efficient explicit full-discrete scheme for strong ap proximation of stochastic Allen–Cahn equation, Stoch. Process. Appl., 130 (2020) 6271-629 9

    work page 2020

  22. [22]

    Wang, R.S

    X.J. Wang, R.S. Qi, A note on an accelerated exponential Euler me thod for parabolic SPDEs with additive noise, Appl. Math. Lett., 46 (2015) 31-37. 25

    work page 2015