Segment convergence for super-linear stochastic functional differential equations by the truncated Euler-Maruyama method
Pith reviewed 2026-05-09 21:10 UTC · model grok-4.3
The pith
The truncated Euler-Maruyama scheme delivers strong convergence of numerical segments for SFDEs with super-linear coefficients.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For SFDEs with super-linear drift and diffusion coefficients, the explicit truncated Euler-Maruyama scheme produces a numerical segment process that remains uniformly bounded in moments over any finite time interval, satisfies an L2-error estimate between its continuous and step-process versions, and converges strongly to the true segment process at a positive order.
What carries the argument
The truncated Euler-Maruyama scheme, which applies a truncation function to the drift and diffusion to prevent explosion while remaining fully explicit.
If this is right
- The numerical segment inherits the ability to approximate invariant measures of the underlying SFDE.
- Ergodicity of the numerical segment follows from the established strong convergence.
- The method supplies a practical tool for pricing path-dependent financial options that depend on the full path segment.
- A supporting numerical example illustrates that the predicted convergence rates are observed in simulation.
Where Pith is reading between the lines
- Segment convergence may allow direct transfer of stability properties from the continuous SFDE to its numerical approximation over longer horizons than the finite-interval proofs cover.
- The truncation technique could be adapted to other explicit schemes for SFDEs, such as Milstein-type methods, to obtain segment results without losing computational simplicity.
- In practice, one could test the method on equations with polynomial growth of varying degrees to measure how the observed order depends on the growth exponent.
Load-bearing premise
The drift and diffusion coefficients obey super-linear growth bounds that still allow a suitable truncation level to produce uniform moment bounds for the numerical segment on finite intervals.
What would settle it
Take any concrete super-linear SFDE, fix a truncation level, and compute the L2 distance between the true segment and the numerical segment at a sequence of step sizes going to zero; if this distance fails to decay at the claimed positive order, the convergence result is false.
Figures
read the original abstract
Most existing literature focuses on pointwise convergence (i.e., convergence at a fixed time point) of numerical solutions for Stochastic functional differential equations (SFDEs). In contrast, this paper investigates the strong segment convergence (i.e., the strong order of convergence of the numerical segment process). For SFDEs with super-linear drift and diffusion coefficients, we employ the explicit truncated Euler-Maruyama (EM) scheme. First, we establish the uniform moment boundedness of the truncated EM solution over a finite time interval. Second, we derive the $L^2$-error estimate between the continuous numerical segment and the step numerical segment. Finally, we prove the strong convergence order of the numerical segment generated by the truncated EM. The results can be used to analyze invariant measures and ergodicity of numerical segment, and have important applications in practical problems such as path-dependent financial options. We also provide a numerical example to support the theoretical results.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript establishes strong segment convergence for the truncated Euler-Maruyama method applied to stochastic functional differential equations (SFDEs) with super-linear drift and diffusion coefficients. It proceeds in three steps: proving uniform moment boundedness of the truncated numerical solution over finite time intervals, deriving the L²-error estimate between the continuous numerical segment and the piecewise-constant step segment, and obtaining the strong convergence order of the numerical segment process. The results are illustrated with a numerical example and noted for applications to invariant measures, ergodicity, and path-dependent problems such as financial options.
Significance. If the central claims hold, the work meaningfully extends the literature on numerical methods for SFDEs from pointwise convergence to segment convergence, which is required for path-dependent functionals. The truncated EM approach is a standard technique for super-linear growth, and the explicit moment bounds plus segment error analysis provide a reusable foundation for studying numerical ergodicity. The three-step architecture is clearly delineated and supports reproducibility when the truncation threshold is fixed as described.
major comments (2)
- [§3] §3 (moment bounds): the uniform moment estimate for the truncated process requires the truncation threshold to be chosen independently of the step-size h; the manuscript should state explicitly how the threshold function is calibrated to the super-linear growth constants so that the bound remains uniform in h (otherwise the subsequent L² segment error may not pass to the limit).
- [§4] §4 (L² segment error): the passage from the continuous numerical segment to the step-function segment invokes a Lipschitz-type estimate on the functional term; the dependence of the constant on the delay length or the modulus of continuity of the functional must be tracked to confirm it does not reduce the final convergence order below the expected rate.
minor comments (2)
- [Abstract] The abstract states that a strong convergence order is proved but does not record the precise order (e.g., ½ or 1); adding this information would improve readability.
- [Numerical example] In the numerical example, the error tables or plots should include a log-log scale or explicit order computation to allow direct visual verification of the theoretical rate.
Simulated Author's Rebuttal
We thank the referee for the thorough review and valuable suggestions. We believe the comments can be addressed with minor revisions to enhance clarity, as detailed in our point-by-point responses below.
read point-by-point responses
-
Referee: [§3] §3 (moment bounds): the uniform moment estimate for the truncated process requires the truncation threshold to be chosen independently of the step-size h; the manuscript should state explicitly how the threshold function is calibrated to the super-linear growth constants so that the bound remains uniform in h (otherwise the subsequent L² segment error may not pass to the limit).
Authors: We appreciate this comment. Upon re-examination, our truncation threshold is indeed chosen independently of h, calibrated solely based on the super-linear growth constants appearing in the drift and diffusion coefficients and the length of the time interval. Specifically, the threshold R is selected larger than a constant determined by these parameters to ensure the moment bounds are uniform with respect to h. We will revise the manuscript in §3 to include an explicit statement of this calibration procedure, including the dependence on the growth constants, to make the uniformity clear and facilitate the passage to the limit in the error estimates. revision: yes
-
Referee: [§4] §4 (L² segment error): the passage from the continuous numerical segment to the step-function segment invokes a Lipschitz-type estimate on the functional term; the dependence of the constant on the delay length or the modulus of continuity of the functional must be tracked to confirm it does not reduce the final convergence order below the expected rate.
Authors: We agree that tracking the dependence is important for rigor. The functional term satisfies a global Lipschitz condition with a constant that depends on the fixed delay length but is independent of the step-size h. Since the delay is fixed and the modulus of continuity is controlled by the Lipschitz property, the constant in the L²-error estimate between the continuous and step segments does not depend on h and thus does not affect the convergence order. We will add a detailed remark in §4 explicitly stating this independence from h and the dependence on the delay length, confirming that the expected strong convergence order is preserved. revision: yes
Circularity Check
No significant circularity; standard three-step stochastic convergence proof
full rationale
The derivation follows the conventional architecture for truncated Euler-Maruyama on super-linear SFDEs: (1) uniform moment bounds on the truncated process over finite intervals under explicit growth hypotheses, (2) L2 error between continuous and piecewise-constant numerical segments, and (3) passage to the limit yielding segment convergence order. These steps are supported by standard stochastic analysis techniques and concrete truncation functions calibrated to the coefficients, without any reduction to self-defined quantities, fitted inputs renamed as predictions, or load-bearing self-citations. The abstract and described structure confirm the chain is self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- truncation threshold
axioms (2)
- domain assumption Existence and uniqueness of solutions to the SFDE under super-linear growth
- standard math Standard Ito calculus and stochastic integral properties
Reference graph
Works this paper leans on
-
[1]
F. Wu, S. Hu, Y . Liu, Positive solution and its asymptotic behaviour of stochastic functional Kolmogorov-type system, Journal of Mathematical Analysis and Applications 364 (2010) 104–118. 22 10 −4 10 −3 10 −2 10 −3 10 −2 10 −1 log (∆ ) log ( √ E∥XT − YT ∥2) Truncated EM, slope=0.5234 Reference Line, slope=0.5 Fig. 1. Rate of strong convergence for Example 5.1
work page 2010
-
[2]
Q. Luo, X. Mao, Y . Shen, Generalised theory on asymptotic stability and boundedness of stochastic functional di fferential equations, Automatica 47 (2011) 2075–2081
work page 2011
-
[3]
J. Bao, G. Yin, C. Y uan, Ergodicity for functional stocha stic di fferential equations and applications, Nonlinear Analysis: Theory, Methods and App lications 98 (2014) 66–82
work page 2014
-
[4]
F. Wu, F. Xi, C. Zhu, On a class of McKean-Vlasov stochasti c functional di fferential equations with applications, Journal of Di fferential Equations 371 (2023) 31–49
work page 2023
-
[5]
B. Shi, Y . Wang, F. Wu, Ergodicity of regime-switching fu nctional diffusions with infinite delay and application to a numerical algorithm for stochast ic optimization, SIAM Journal on Control and Optimization 50 (2022) 2658–2683
work page 2022
-
[6]
Mao, Numerical solutions of stochastic functional di fferential equations, LMS J
X. Mao, Numerical solutions of stochastic functional di fferential equations, LMS J. Com- put. Math 6 (2003) 1461–1570
work page 2003
-
[7]
J. Zhao, Y . Yi, Y . Xu, Strong convergence and stability of the split-step theta method for highly nonlinear neutral stochastic delay integro di fferential equation, Applied Numerical Mathematics 157 (2020) 385–404
work page 2020
-
[8]
X. Li, X. Mao, G. Song, An explicit approximation for supe r-linear stochastic functional differential equations, Stochastic Processes and their Applic ations 169 (2024) 104275
work page 2024
-
[9]
B. Shi, Y . Wang, X. Mao, F. Wu, Approximation of invariant measures of a class of back- ward Euler-Maruyama scheme for stochastic functional di fferential equations, Journal of Differential Equations 389 (2024) 415–456
work page 2024
- [10]
-
[11]
S. Zhou, Invariant measure of approximating method for stochastic functional di fferential equation, Communications on Applied Mathematics and Compu tation (2025) 1–21. 23
work page 2025
-
[12]
E. Buckwar, The Theta-Maruyama scheme for stochastic f unctional differential equations with distributed memory term, Monte Carlo Methods Appl. 10 ( 2004) 235–244
work page 2004
- [13]
-
[14]
S. Zhou, H. Jin, Strong convergence of implicit numeric al methods for nonlinear stochastic functional differential equations, Journal of Computational and Applied M athematics 324 (2017) 241–257
work page 2017
- [15]
-
[16]
J. Zhou, H. Xue, X. Mao, Truncated Euler-Maruyama metho d for hybrid stochastic func- tional differential equations with infinite time delay, Journal of Comp utational and Applied Mathematics 472 (2026) 116773
work page 2026
-
[17]
X. Mao, The truncated Euler-Maruyama method for stocha stic differential equations, Jour- nal of Computational and Applied Mathematics 290 (2015) 370 –384
work page 2015
-
[18]
F. Sancier, S. Mohammed, An option pricing model with me mory, Communications on Stochastic Analysis 11 (2017) 505–532
work page 2017
-
[19]
M. Hutzenthaler, A. Jentzen, Numerical approximation s of stochastic differential equations with non-globally Lipschitz continuous coe fficients, Memoirs of the American Mathemat- ical Society 236 (2015) 1–104
work page 2015
-
[20]
F. Wu, G. Yin, H. Mei, Stochastic functional di fferential equations with infinite delay: Existence and uniqueness of solutions, solution maps, Mark ov properties, and ergodicity, Journal of Differential Equations 262 (2017) 1226–1252
work page 2017
-
[21]
X. Li, X. Mao, G. Yin, Explicit numerical approximation s for stochastic di fferential equa- tions in finite and infinite horizons: truncation methods, co nvergence in pth moment and stability, IMA Journal of Numerical Analysis 39 (2019) 847– 892
work page 2019
-
[22]
D. J. Higham, X. Mao, A. M. Stuart, Strong convergence of Euler-type methods for non- linear stochastic di fferential equations, SIAM Journal on Numerical Analysis 40 ( 2003) 1041–1063. 24
work page 2003
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.