REVIEW 3 major objections 3 minor 1 cited by
$W_{\bf d}$-convergence rate of EM schemes for invariant measures of supercritical stable SDEs
T0 review · 3 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read This paper proves that the invariant measure of an Euler–Maruyama discretization of an SDE driven by multiplicative α-stable noise converges to the true invariant measure at an explicit polynomial rate in the step size, for α∈(1/2,2).
desk verdict First quantitative EM invariant-measure rates for supercritical stable SDEs; two repairable but real gaps in the Stein-equation and regularity arguments. 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 central object is the nonlocal Stein/Poisson equation $\mathcal{L}f=\mu(g)-g$ for the stable generator $\mathcal{L}f=b\cdot\nabla f+\mathcal{L}^{(\alpha)}_{\sigma(x)}f$. The argument shows that $f=\int_0^\infty P_t\bar g\,dt$ belongs to weighted Hölder spaces with weight $(\log(e+|x|^2))^{-1}$, using two-sided heat-kernel estimates for the SDE density with unbounded drift, commutator estimates, and a third-difference regularity lemma. This regularity is what converts the generator error $\mathbb{E}[\mathcal{L}f(\xi)]$ for a stationary EM step ξ into the explicit powers of η that appear in the rate.
What would settle it
Take the one-dimensional OU process (B.1) for α∈(1/2,1] with the metric $d(x,y)=|x-y|^\gamma\wedge 1$. Proposition B.1 gives the exact invariant distributions of both the SDE and its EM scheme, so $W_d(\mu,\mu_\eta)$ can be computed or simulated to high precision; if its scaling in η is slower than $\eta^{\gamma\wedge(2\alpha-1-\varepsilon)}$ for all small ε, then Theorem 1.2 is false. Alternatively, finding a drift satisfying (Hγ) for which the two-sided heat-kernel bound (1.7)–(1.8) fails would invalidate the regularity lemmas and the rate.
Extended reading notes
Core claim
Under the dissipativity assumption $\langle x,b(x)\rangle\le -\kappa_1|x|^2+\kappa_2$ and the Hölder regularity assumption (Hγ), the invariant measure μη of the Euler–Maruyama scheme satisfies $W_d(\mu,\mu_\eta)\le C(\eta^{\gamma\wedge(2\alpha-1-\varepsilon)})$ for α∈(1/2,1] and $W_d(\mu,\mu_\eta)\le C\eta^{\gamma/\alpha}$ for α>1, where $d(x,y)=|x-y|^\gamma\wedge 1$ and γ∈((1-α)_+,1]. The discovery is that the Stein/Poisson equation $\mathcal{L}f=\mu(g)-g$ has a solution $f=\int_0^\infty P_t\bar g\,dt$ lying in weighted Hölder spaces with logarithmic growth, $C^1_{\log}$ and $C^{\beta+\gamma}_{\log}$, for every β∈(0,α); this regularity is exactly strong enough to control the one-step generator error along the stationary EM chain and to yield the first explicit discretization rate for supercritical stable SDEs with multiplicative noise.
Load-bearing premise
The proof requires the two-sided heat-kernel estimate (1.7)–(1.8) to hold for the whole coefficient class (Hγ) with unbounded drift; if that estimate fails for some admissible drift and diffusion, the Stein-solver regularity and the stated rate would collapse.
Editorial extensions
If this is right
- For α∈(1/2,1], the invariant measure of the EM chain converges under $W_d$ at an explicitly polynomial rate, so sampling from heavy-tailed targets by Euler–Maruyama is justified with a rate.
- For α∈(1,2) with monotone drift, the stronger Wasserstein-1 rate $\eta^{1/\alpha}$ holds, matching the time-scaling of stable increments.
- The Stein/Poisson regularity estimates apply to the whole supercritical class with multiplicative noise and unbounded drift, so the same machinery can bound discretization errors for other observables beyond the Wasserstein distance.
- Choosing γ=1 when α>1 recovers a bounded-Lipschitz-type metric, while choosing γ close to $(1-\alpha)_+$ covers the finest Hölder metric available to the method.
Reading between the lines
- Editorial inference: the logarithmic-growth weight mechanism suggests the same $W_d$ rate should extend to the Pareto-type EM scheme mentioned in Remark 1.4, since only the small-time scaling of the increments is used.
- Editorial inference: the $\varepsilon$ loss in the α∈(1/2,1] rate looks like an artifact of taking a limit β↑α in the proof, not an intrinsic obstruction; the one-dimensional OU example in Appendix B offers a testbed for sharpness.
- Editorial inference: for α≤1/2 the short-time gradient estimate diverges, so a genuinely new estimate would be needed; if a rate were sought in an unweighted metric such as $W_1$, the logarithmic-growth regularity would not suffice.
- Editorial inference: the same Stein/Poisson regularity could be used to prove rates for other discretizations or for time-averaged functionals, because the central bound is about the semigroup, not about the specific EM recursion.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the Euler-Maruyama (EM) approximation of the invariant measure of supercritical stable SDEs with multiplicative noise. Under Hölder regularity of the coefficients and a dissipativity condition, the authors prove Wasserstein convergence rates for the invariant measures of the EM scheme: for α∈(1/2,2) and the bounded-Hölder-type distance d(x,y)=|x-y|^γ∧1, Theorem 1.2 gives a rate of η^{γ∧(2α-1-ε)} for α∈(1/2,1] and η^{γ/α} for α>1, up to logarithmic factors. The proof is based on Stein's method: the authors establish weighted Hölder regularity for solutions of nonlocal Stein/Poisson equations, using two-sided heat-kernel bounds for stable SDEs with unbounded drift, and then estimate the discrepancy between the EM semigroup and the continuous semigroup. An appendix gives matching first-order rates for Ornstein-Uhlenbeck processes, showing that the rates are at least sharp in that model.
Significance. If correct, the main result is the first explicit discretization rate for invariant measures of supercritical stable SDEs with multiplicative noise, going substantially beyond the subcritical case treated by Chen, Deng, Schilling and Xu. The paper also provides a useful Stein-equation regularity theory for nonlocal generators with unbounded coefficients and only Hölder-regular data. Strengths include the explicit nature of the estimates, the use of external heat-kernel and ergodicity theorems rather than fitted constants, and the sharp OU example in Appendix B that calibrates the rate. The two technical gaps identified below concern the rigorous derivation of the Stein equation and one norm-equivalence step; both appear repairable without changing the main claims, but they currently prevent the proof from being fully valid as written.
major comments (3)
- [§2.4, Proposition 2.15] The Hille–Yosida argument proving (2.45) is not valid as written. Since A is the generator of a C0-semigroup on C0, one has D(A)⊂C0; however, Theorem 2.14, specifically (2.40), shows that f(x)=∫_0^∞ P_t ḡ(x)dt has logarithmic growth. For non-centered g∈C0, the function ḡ=g-μ(g) does not vanish at infinity, so f is not in C0 and the limits lim_{ε↓0}U_ε ḡ=f and lim_{ε↓0}g_ε=-ḡ are not C0-limits. Thus the closure argument cannot be applied. Because Theorem 2.16 invokes (2.45) for the compactly supported truncations g_n, the pointwise identity Lf=-ḡ — the foundation of the Stein method in Theorems 1.1 and 1.2 — is not rigorously established as written. This gap is likely repairable by differentiating under the integral and using the regularity of Theorem 2.14 together with exponential ergodicity (2.20), but the argument must be supplied.
- [§2.3, Lemma 2.13 and equation (2.32)] In the small-time part of Lemma 2.13, Theorem 2.3(i) is applied with β0=(β+γ)∧2. The equivalent-characterization statement in Theorem 2.3 is stated for non-integer β0 only. When β+γ≥2, which can occur for α>1 with γ=1 and β close to α, one has β0=2, an integer, and the inequality sup_{|v|≤1}‖δ^3_v P_t g‖∞/|v|^2 does not by itself justify the C^2-norm bound (2.32). Since (2.32) is used in the large-time estimate (2.34) and in Theorem 2.14, this is a load-bearing gap. The fix is local: replace β0=(β+γ)∧2 by β0=(β+γ)∧(2-ε) with ε small so that β0>α and non-integer, and adjust the subsequent constants and τ accordingly.
- [§2.3, Lemma 2.12] In the proof of Lemma 2.12, the estimate E|f(X^{x±v-z}_t)| ≲ (1+|x±v-z|^m)‖f‖_{L∞_m} for m∈[0,α) is asserted as 'standard' without proof or reference. This estimate is used with f equal to [δ_v,L]P_{s+t0}g in Lemma 2.13, so it is needed in a form that supplies the weighted L∞ norm of the commutator term. A short proof or a precise citation should be provided.
minor comments (3)
- [§3, Proof of Theorem 1.1] The references '(3.5)' after the estimates for I3 and I4 point to an equation that is introduced only later in the proof of Theorem 1.2; the intended reference appears to be (3.1).
- [§2.4, Theorem 2.14] The conclusion ‖f‖_{C^{β+γ}_{log}} is stated for every β∈(0,α), but the proof via Lemma 2.13 and Theorem 2.3(i) requires β+γ to be non-integer. Either state the result for non-integer β+γ only, or add an interpolation argument for the exceptional integer values such as β+γ=2.
- [Throughout] The displayed text contains numerous formatting artifacts, such as '/BD' and misplaced inequality signs, which should be cleaned in the production version. The abstract also contains a typo ('Stei n').
Circularity Check
No significant circularity: the W_d-rate theorem is a genuine derivation from external heat-kernel and ergodicity estimates; self-citations in [29] and [48] are published theorems with independent proofs and do not smuggle in the target result.
full rationale
The central claim, Theorem 1.2, is not an input or a renamed fit. The proof obtains the rate by combining the Stein/Poisson representation (1.13)-(1.15), the regularity estimate for the Stein solution in Theorem 2.14, and direct EM local-error estimates in Section 3. No parameter is fitted to the target quantity, and the Wasserstein metric d(x,y)=|x-y|^gamma ∧ 1 is chosen with the same gamma as the coefficient Holder exponent, but this is a metric choice, not a reduction of the conclusion to the assumption. The cited results [29] and [48] are external published theorems with proofs under assumptions (H_gamma) and (1.9); those assumptions do not include the W_d-convergence rate, so the citations are independent support rather than circular self-citation. The main derivation is self-contained given those external inputs. The reviewer-identified issue in Proposition 2.15 and Theorem 2.16 is a proof gap, not circularity: f = ∫ P_t g-bar dt is shown by the paper itself to have only logarithmic growth, so it need not lie in C0 = D(A), and the Hille-Yosida resolvent passage in C0 is not justified as written. This is a genuine correctness concern, but it does not equate the theorem with its inputs or relabel a fitted value as a prediction; therefore it does not raise the circularity score.
Assumptions & free parameters
assumptions (4)
- domain assumption Two-sided heat kernel estimates (1.7)-(1.8) from Menozzi and Zhang [29] hold for the SDE (1.2) under (Hγ).
- domain assumption Exponential ergodicity with m-th moment bounds (2.20) from Zhang and Zhang [48] holds for the semigroup under dissipativity (2.19).
- domain assumption Gradient estimate (2.12) from Liang and Wang [27] holds in the monotone case (H1) and (2.11).
- standard math Interpolation and norm equivalence for weighted Hölder spaces in Theorem 2.3, taken from Hao, Zhang, Zhu and Zhu [21].
Cite this review
Pith. "Pith review of $W_{\bf d}$-convergence rate of EM schemes for invariant measures of supercritical stable SDEs." pith.science (2026). https://pith.science/paper/HIJ3RK6I
@misc{pith2026241109949,
author = {Pith},
title = {Pith review of: $W_\bf d$-convergence rate of EM schemes for invariant measures of supercritical stable SDEs},
year = {2026},
howpublished = {\url{https://pith.science/paper/HIJ3RK6I}},
note = {Machine review of arXiv:2411.09949}
}
abstract
By establishing the regularity estimates for nonlocal Stein/Poisson equations under $\gamma$-order H\"older and dissipative conditions on the coefficients, we derive the $W_{\bf d}$-convergence rate for the Euler-Maruyama schemes applied to the invariant measure of SDEs driven by multiplicative $\alpha$-stable noises with $\alpha \in (\frac{1}{2}, 2)$, where $W_{\bf d}$ denotes the Wasserstein metric with ${\bf d}(x,y)=|x-y|^\gamma\wedge 1$ and $\gamma \in ((1-\alpha)_+, 1]$.
Forward citations
Cited by 1 Pith paper
-
Discretization, Uniform-in-Time Estimations and Approximation of Invariant Measures for Nonlinear Stochastic Differential Equations with Non-Uniform Dissipativity
The numerical invariant measure of an explicit truncated Euler–Maruyama scheme converges to the true invariant measure at order h^{1/2} in the L1-Wasserstein distance, uniformly in time, for drifts that are dissipativ...
Reference graph
Works this paper leans on
-
[8]
P. Chen, C. S. Deng, R.L. Schilling and L. Xu (2023): Appro ximation of the invariant measure of stable SDEs by an Euler-Maruyama scheme. Stochastic Processes and their Applications. [J]. 163, pp. 136-167
work page 2023
-
[29]
S. Menozzi, and X. Zhang (2022): Heat kernel of supercri tical nonlocal operators with unbounded drifts. Journal de l’ ´Ecole polytechnique — Math´ ematiques. [J]. 9, pp. 537-579
work page 2022
-
[48]
X. Zhang and X. Zhang (2023): Ergodicity of supercritic al SDEs driven by α -stable processes and heavy-tailed sampling. Bernoulli. [J]. 29(3), pp. 1933-1958. (P. Chen) School of Mathematics, Nanjing University of Aeronautics and As tronautics, Nanjing, 211106, China Email address : chenpengmath@nuaa.edu.cn (L. Xu) 1. Department of Mathematics, F aculty o...
work page 2023
-
[1]
B. C. Arnold (2015): Pareto Distributions, 2nd ed. Chapman and Hall/CRC
work page 2015
-
[2]
B. Arras and C. Houdr´ e (2019): On Stein’s method for mult ivariate self-decomposable laws. Electronic Journal of Probability. [J]. 24(128), pp. 1-63
work page 2019
-
[3]
V. Bally and D. Talay (1996): The law of the Euler scheme fo r stochastic differential equations. Probability theory and related fields. [J]. 104(1), pp. 43-60
work page 1996
-
[4]
J. Bao, X. Mao and C. Yuan (2009): Rate of convergence for n umerical solutions to SFDEs with jumps. arXiv preprint arXiv:0906.3455
work page Pith review arXiv 2009
-
[5]
R. F. Bass (2009): Regularity results for stable-like op erators. Journal of Functional Analysis. [J]. 8, pp. 2693-2722
work page 2009
Show all 48 references
-
[6]
N. H. Bingham, C. M. Goldie and J. L. Teugels (1987): Regul ar variation. Cambridge University Press, Cambridge . Wd-CONVERGENCE RATE OF EM SCHEMES FOR INV ARIANT MEASURES 23
1987
-
[7]
Bubeck, R
S. Bubeck, R. Eldan and J. Lehec (2018): Sampling from a lo g-concave distribution with projected Langevin Monte Carlo. Discrete & Computational Geometry. [J]. 4, pp. 757-783
2018
-
[9]
P. Chen, Q. M. Shao and L. Xu (2023): A probability approxi mation framework: Markov process approach. The Annals of Applied Probability. [J]. 33(2), pp. 1619-1659
2023
-
[10]
Z. Q. Chen, R. Song and X. Zhang (2018): Stochastic flows f or L´ evy processes with H¨ older drifts. Revista matem´atica iberoamericana. [J]. 34(4), pp. 1755-1788
2018
-
[11]
Z. Q. Chen, X. Zhang and G. Zhao (2021): Supercritical SD Es driven by multiplicative stable-like L´ evy processes. Transactions of the American Mathematical Society. [J]. 374(11), pp. 7621-7655
2021
-
[12]
A. S. Dalalyan (2017): Theoretical guarantees for appr oximate sampling from smooth and log-concave densities. Journal of the Royal Statistical Society: Series B (Statist ical Methodology). [J]. 79(3), pp. 651-676
2017
-
[13]
A. S. Dalalyan and L. Riou-Durand (2020): On sampling fr om a log-concave density using kinetic Langevin diffusions. Bernoulli. [J]. 26(3), pp. 1956-1988
2020
-
[14]
Durrett (2010): Probability: Theory and Examples , 4th ed
R. Durrett (2010): Probability: Theory and Examples , 4th ed. Cambridge University Press
2010
-
[15]
Eberle, A
A. Eberle, A. Guillin and R. Zimmer (2019): Couplings an d quantitative contraction rates for Langevin dy- namics. The Annals of Probability. [J]. 47(4), pp. 1982-2010
2019
-
[16]
Embrechts, C
P. Embrechts, C. Kl¨ uppelberg and T. Mikosch (1997): Mo delling extremal events: For insurance and finance. Springer-Verlag, Berlin
1997
-
[17]
Fang and M
W. Fang and M. B. Giles (2016): Adaptive Euler-Maruyama method for SDEs with non-globally Lipschitz drift. In: International Conference on Monte Carlo and Quasi-Monte Ca rlo Methods in Scientific Computing . Springer, Cham, pp. 217-234
2016
-
[18]
X. Fang, Q. M. Shao and L. Xu (2019): Multivariate approx imations in W asserstein distance by Stein’s method and Bismut’s formula. Probability Theory and Related Fields. [J]. 174(3), pp. 945-979
2019
-
[19]
S. Foss, D. Korshunov and S. Zachary (2011): An introduc tion to heavy-tailed and subexponential distributions. Springer Series in Operations Research and Financial Engin eering. Springer, New York
2011
-
[20]
Gilbarg and N
D. Gilbarg and N. Trudinger (1983): Elliptic partial di fferential equations of second order. 2nd ed. Grundlehren der mathematischen Wissenschaften [Fundamental Principl es of Mathematical Sciences]. Springer-Verlag, Berlin
1983
-
[21]
Z. Hao, X. Zhang, R. Zhu and X. Zhu (2024): Singular kinet ic equations and applications. The Annals of Probability. [J]. 52(2), pp. 576-657
2024
-
[22]
Huang and Z
X. Huang and Z. W. Liao (2018): The Euler-Maruyama metho d for S(F)DEs with H¨ older drift and α -stable noise. Stochastic Analysis and Applications. [J]. 36(1), pp. 28-39
2018
-
[23]
S. Kim, Q. Song and F. Liang (2022): Stochastic gradient Langevin dynamics with adaptive drifts. Journal of Statistical Computation and Simulation. [J]. 92(2), pp. 318-336
2022
-
[24]
K¨ uhn (2021): Schauder estimates for Poisson equati ons associated with non-local Feller generators
F. K¨ uhn (2021): Schauder estimates for Poisson equati ons associated with non-local Feller generators. Journal of Theoretical Probability [J]. 34(3), pp. 1506–1578
2021
-
[25]
K¨ uhn (2022): Interior Schauder estimates for ellip tic equations associated with L´ evy operators
F. K¨ uhn (2022): Interior Schauder estimates for ellip tic equations associated with L´ evy operators. Potential Analysis. [J]. 56(3), pp. 459-481
2022
-
[26]
K¨ uhn and R
F. K¨ uhn and R. L. Schilling (2019): Strong convergence of the Euler-Maruyama approximation for a class of L´ evy-driven SDEs.Stochastic Processes and their Applications. [J]. 129(8), pp. 2654-2680
2019
-
[27]
Liang and J
M. Liang and J. W ang (2020): Gradient estimates and ergo dicity for SDEs driven by multiplicative L´ evy noises via coupling. Stochastic Processes and their Applications. [J]. 130(5), pp. 3053-3094
2020
-
[28]
S. L. Lohr (2010): Sampling: Design and Analysis. 2nd ed . Advanced Series, Boston
2010
-
[30]
S. P. Meyn and R. L. Tweedie (1992): Stability of Markovi an processes I: Criteria for discrete-time chains. Advances in Applied Probability. [J]. 24(3), pp. 542-574
1992
-
[31]
T. H. Nguyen, U. Simsekli and G. Richard (2019): Non-asy mptotic analysis of fractional Langevin Monte Carlo for non-convex optimization. International Conference on Machine Learning. [M]. pp. 4810-4819
2019
-
[32]
Pag` es and F
G. Pag` es and F. Panloup (2023): Unadjusted Langevin al gorithm with multiplicative noise: Total variation and W asserstein bounds. The Annals of Applied Probability. [J]. 33(1), pp. 726-779
2023
-
[33]
Pardoux and A
E. Pardoux and A. Yu. Veretennikov (2001): On the Poisso n equation and diffusion approximation. I The Annals of Probability. [J]. 29, pp. 1061-1085
2001
-
[34]
Pardoux and A
E. Pardoux and A. Yu. Veretennikov (2003): On the Poisso n equation and diffusion approximation. II The Annals of Probability. [J]. 31, pp. 1166-1192
2003
-
[35]
Pardoux and A
E. Pardoux and A. Yu. Veretennikov (2005): On the Poisso n equation and diffusion approximation. III The Annals of Probability. [J]. 33, pp. 1111-1133
2005
-
[36]
Pazy (1983): Semigroups of linear operators and appl ications to partial differential equations
A. Pazy (1983): Semigroups of linear operators and appl ications to partial differential equations. Springer- Verlag, New York. 44, pp. viii+279
1983
-
[37]
Protter and D
P. Protter and D. Talay (1997): The Euler scheme for L´ ev y driven stochastic differential equations. The Annals of Probability. [J]. 25(1), pp. 393-423. 24 P. CHEN, L. XU, X. ZHANG, AND X. ZHANG
1997
-
[38]
Revuz and M
D. Revuz and M. Yor (1999): Continuous martingales and B rownian motion, 3th ed. Springer-Verlag, Berlin
1999
-
[39]
G. O. Roberts and R. L. Tweedie (1996): Exponential conv ergence of Langevin distributions and their discrete approximations. Bernoulli. [J]. 2(4), pp. 341-363
1996
-
[40]
R¨ ockner and L
M. R¨ ockner and L. Xie (2021): Diffusion approximation f or fully coupled stochastic differential equations. The Annals of Probability. [J]. 49(3), pp. 1205-1236
2021
-
[41]
Ros-Oton and J
X. Ros-Oton and J. Serra (2016): Regularity theory for g eneral stable operators. Journal of Differential Equa- tions. [J]. 260, pp. 8675–8715
2016
-
[42]
Y. R. Rubinstein and D. P. Kroese (2016): Simulation and the Monte Carlo Method. Wiley Series in Probability and Statistics
2016
-
[43]
E. M. Stein (1970): Singular integrals and differentiab ility properties of functions. Princeton University Press, Princeton, NJ
1970
-
[44]
Triebel (2006): Theory of function spaces
H. Triebel (2006): Theory of function spaces. III Birkh¨ auser Verlag, Basel
2006
-
[45]
Villani (2009): Optimal transport, Old and new
C. Villani (2009): Optimal transport, Old and new. Springer-Verlag, Berlin
2009
-
[46]
F. Y. W ang, L. Xu and X. Zhang (2015): Gradient estimates for SDEs driven by multiplicative L´ evy noise. Journal of Functional Analysis. [J]. 269(10), pp. 3195-3219
2015
-
[47]
Zhang (2010): Stochastic Volterra equations in Bana ch spaces and stochastic partial differential equations
X. Zhang (2010): Stochastic Volterra equations in Bana ch spaces and stochastic partial differential equations. Journal of Functional Analysis. [J]. 258, pp. 1361-1425
2010
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.