---
title: Modified Tamed Scheme Techniques for SDEs
url: https://www.emergentmind.com/topics/modified-tamed-scheme
type: topic
---

# Modified Tamed Scheme Techniques for SDEs

Searching arXiv for recent papers on modified/tamed schemes to ground the article in current literature.
arxiv_search(query="modified tamed scheme stochastic differential equations", max_results=10)
to=arxiv_search  荣富?
{"query":"modified tamed scheme stochastic differential equations", "max_results": 10}
Modified tamed scheme denotes a class of numerical discretizations in which the basic taming device used for nonglobally Lipschitz stochastic dynamics is altered so that explicit or semi-implicit methods remain stable under superlinear growth, low regularity, singular interactions, algebraic constraints, or long-time sampling requirements. In the cited literature, the modification may take the form of adaptive time stepping, a tamed Milstein correction, a smooth cut-off that activates taming only in the far field, a truncation near a singularity set, or a problem-specific implicit-explicit split. These constructions are used for strong approximation of SDEs, SDAEs, McKean–Vlasov equations, and related Langevin samplers, with convergence orders ranging from the Euler-type rate \(1/2\) to Milstein-type order \(1\), and with additional results on infinite-horizon error bounds, geometric ergodicity, and relative-entropy control [2411.01849] [2507.09475] [2412.19377].

## 1. Origins in taming and explicit stabilization

The conceptual starting point is the failure of classical explicit discretizations for SDEs with superlinearly growing coefficients. For such equations, the explicit Euler method can fail to converge strongly, and standard explicit Milstein analysis typically requires substantially more regularity than is available in many nonglobally Lipschitz models. Taming was introduced to preserve explicit implementability while suppressing increments that would otherwise become too large. In the tamed Euler setting, this restores strong order \(1/2\) in nonstandard problems such as stochastic differential equations with piecewise continuous arguments, while Milstein-type taming can recover higher-order strong accuracy in commutative-noise settings [1502.00058] [1102.0662].

Within this development, the adjective “modified” does not designate a single canonical algorithm. Rather, it refers to alterations of the standard taming mechanism that are tailored to a specific obstruction: low differentiability of coefficients, singular or constrained geometry, superlinear drift combined with long-time ergodicity, or loss of weak order under uniform denominator-based taming. A general abstract version appears in the framework of modified Euler recursions
\[
X_{(n+1)\delta}=T^{(\delta)}(X_{n\delta})+b^{(\delta)}\!\bigl(T^{(\delta)}(X_{n\delta})\bigr)\delta+\sqrt{\delta}\,\xi_{n+1},
\]
where \(T^{(\delta)}\) is contractive and \(b^{(\delta)}\to b\) as \(\delta\to 0\). This framework explicitly encompasses the tamed Euler scheme and the truncated Euler scheme [2412.19377].

## 2. Principal forms of modification

The literature uses several distinct mechanisms to modify a tamed scheme. They differ not only in implementation but in the analytical role played by the modification.

| Variant | Defining modification | Representative paper |
|---|---|---|
| Adaptive tamed Milstein | Adaptive step size and tamed Milstein derivative term \(q_\Delta\) | [2411.01849] |
| Far-field cut-off taming | \(b^h=b/(1+\psi(\gamma h^\alpha |b|))\), with no taming for small \(|b|\) | [2507.09475] |
| Semi-implicit taming | Linear drift implicit, nonlinear drift explicitly tamed | [2509.09032] |
| Singular-set cut-off | Drift turned off when \(\rho(x)\) is too small | [2408.08118] |
| Contractive modified Euler | Contractive state map \(T^{(\delta)}\) and modified drift \(b^{(\delta)}\) | [2412.19377] |

These variants share the same strategic objective: preserve explicitness or near-explicitness while controlling excursions that invalidate classical convergence theory. Their differences are structural. Adaptive Milstein schemes modify both local step length and Milstein correction. Semi-implicit schemes use taming only on the nonlinear component, leaving linear stiff structure implicit. Singular-drift schemes replace denominator taming by a geometric cut-off tied to distance from a singular set. Far-field schemes use a smooth cut-off so that the discretization coincides exactly with the untamed method on the region where the drift is moderate.

This diversity has an important methodological consequence. A modified tamed scheme is best understood not as a single formula but as a design principle: alter the taming only where instability or lack of regularity actually enters the proof.

## 3. Adaptive Milstein under low regularity

A particularly sharp realization appears in the tamed-adaptive Milstein scheme for the one-dimensional SDE
\[
X_t=x_0+\int_0^t \mu(X_s)\,ds+\int_0^t \sigma(X_s)\,dW_s,\qquad x_0\in\mathbb R,\ t\ge 0.
\]
The novelty is that the coefficients are assumed to satisfy only
\[
\mu,\sigma\in C_l^{1+\alpha},
\]
meaning that the first derivatives are locally Hölder continuous of order \(\alpha\in(0,1]\) with polynomial growth, rather than globally Lipschitz. The scheme is defined by
\[
t_0=0,\qquad \widehat X_0=x_0,\qquad t_{k+1}=t_k+h(\widehat X_{t_k})\Delta,
\]
and
\[
\widehat X_{t_{k+1}}
=
\widehat X_{t_k}
+\mu(\widehat X_{t_k})(t_{k+1}-t_k)
+\sigma(\widehat X_{t_k})(W_{t_{k+1}}-W_{t_k})
+\frac12 q_\Delta(\widehat X_{t_k})\Big((W_{t_{k+1}}-W_{t_k})^2-(t_{k+1}-t_k)\Big),
\]
with tamed Milstein coefficient
\[
q_{\Delta}(x)=\frac{\sigma'(x)\sigma(x)}{1+\Delta^{1/2}|\sigma'(x)\sigma(x)|}.
\]
The adaptive step-size function is
\[
h(x)=\frac{h_0}{(1+|\mu(x)|^2+|\mu'(x)|+|\sigma(x)|^4+|\sigma'(x)|^4+|q_\Delta(x)|+|x|^{l_0})^2}.
\]
Relative to classical Milstein, the method therefore changes both the step size and the derivative correction. Relative to standard tamed Milstein, it is explicitly tamed-adaptive rather than merely coefficient-tamed [2411.01849].

The assumptions are organized as \(A0\)–\(A2\). Besides \(C_l^{1+\alpha}\)-regularity, they impose a coercivity-type moment condition
\[
x\mu(x)+\frac{p_0-1}{2}\sigma^2(x)\le \gamma x^2+\eta,
\]
and a one-sided stability condition
\[
(x-y)(\mu(x)-\mu(y))+\frac12|\sigma(x)-\sigma(y)|^2\le \lambda |x-y|^2.
\]
The local low-regularity structure enters through the remainder estimate
\[
|f(x)-f(y)-f'(y)(x-y)|\le C(1+|x|^l+|y|^l)|x-y|^{1+\alpha}.
\]

Under these assumptions, if
\[
l_0 \ge \max\Big\{2,\frac{4l}{3(1+\alpha)}\Big\},\qquad p_0\ge 4(l+\alpha+1),
\]
then for every \(T>0\),
\[
\sup_{0\le t\le T}\mathbb E\big[|\widehat X_t-X_t|^2\big]\le C_T\Delta^{1+\alpha}.
\]
Hence the \(L_2\)-strong convergence rate is
\[
\frac{1+\alpha}{2}.
\]
If \(\lambda<0\) and \(\gamma<0\), the same estimate holds uniformly on the infinite horizon,
\[
\sup_{t\ge 0}\mathbb E\big[|\widehat X_t-X_t|^2\big]\le C\Delta^{1+\alpha},
\]
with \(C\) independent of \(T\). The same rate is therefore valid on both finite intervals and \((0,\infty)\) [2411.01849].

The proof combines moment bounds, adaptive-step control showing \(t_k\uparrow\infty\) almost surely, increment estimates, a Yamada–Watanabe smoothing argument, and a Gronwall-type estimate for the error process. A central estimate is
\[
|q_\Delta(x)-\sigma'(x)\sigma(x)|\le C\Delta^{1/2}|\sigma'(x)\sigma(x)|^2,
\]
which stabilizes the Milstein correction under only local Hölder continuity of derivatives. The resulting rate interpolates directly with the smoothness index \(\alpha\): smoother derivatives improve the strong order, but global Lipschitz differentiability is no longer required.

## 4. Far-field taming and preservation of Euler order

A different modification is designed to remove a specific weakness of standard denominator taming. In the usual tamed Euler construction,
\[
b^h_{\mathrm{tame}}(t,x)=\frac{b(t,x)}{1+h^\alpha |b(t,x)|},\qquad \alpha\in(0,1/2],
\]
the drift perturbation is present everywhere, including regions where \(|b|\) is moderate. The modified tamed scheme replaces this by
\[
b^h(t,x)=\frac{b(t,x)}{1+\psi(\gamma h^\alpha |b(t,x)|)},
\]
where \(\psi:\mathbb R_+\to\mathbb R_+\) is smooth and satisfies
\[
\psi(r)=
\begin{cases}
0, & r\le 1,\\
r, & r\ge 2,
\end{cases}
\]
with smooth increasing interpolation for \(r\in(1,2)\). Consequently, if \(\gamma h^\alpha |b(t,x)|\le 1\), then
\[
b^h(t,x)=b(t,x),
\]
so no taming occurs in the bulk region. In the far field one recovers a standard tamed-type denominator, and the bound
\[
|b^h(t,x)|\le \min\!\left(2\gamma^{-1}h^{-\alpha},\,|b(t,x)|\right)
\]
still provides stability [2507.09475].

For the Euler discretization,
\[
\hat X_{k+1}=\hat X_k+h\,b^h(T_k,\hat X_k)+\sigma(T_k,\hat X_k)\Delta W_k,
\]
this selective taming changes the error mechanism. The paper proves the high-order estimate
\[
\mathbb{E}\left|b^\xi(\kappa(t),\hat{X}_h(\kappa(t))) -b^{\xi,h}(\kappa(t),\hat{X}_h(\kappa(t)))\right|^p \le Ch^q
\]
for any \(p,q\ge 1\). The analytical reason is that taming is activated only on the event \(|b^\xi|\gtrsim h^{-\alpha}\), whose probability is made arbitrarily small by moment bounds. This removes the built-in low-order perturbation that standard taming introduces everywhere [2507.09475].

The resulting convergence theory recovers the underlying Euler orders. Under the strong assumptions,
\[
\sup_{0\le t\le T}\mathbb{E}|X(t)-\hat{X}_h(t)|^2 \le Ch,
\]
so the modified tamed Euler scheme has strong order \(1/2\). Under the weak assumptions, for any \(f\in C_p^\infty(\mathbb{R}^d)\),
\[
\sup_{0\le t\le T}\left|\mathbb{E}[f(X(t))]-\mathbb{E}[f(\hat{X}_h(t))]\right|\le Ch,
\]
so the weak order is \(1\). The paper explicitly contrasts this with standard tamed schemes, which are described as introducing an \(O(h^\alpha)\)-type perturbation that can lower weak order [2507.09475].

The same construction is extended to random-batch approximations and to stochastic gradient Langevin dynamics. For the Langevin application, the paper proves the uniform-in-time relative-entropy bound
\[
\sup_{t\ge 0}\mathcal{H}(\bar{\rho}_t^{(h)}\mid \rho_t)\le Ch^{2-\epsilon},
\]
for any \(\epsilon\in(0,1)\) with \(\alpha=\epsilon/2\). This yields a uniform-in-time near-sharp error estimate under relative entropy for super-linear sampling dynamics [2507.09475].

## 5. Structured, constrained, and singular variants

Problem-specific modifications become especially pronounced once the underlying stochastic model is no longer a standard SDE. For index-1 SDAEs with a singular matrix \(A(t)\), the modified construction is semi-implicit rather than fully explicit. The update
\[
A(t_n)X_{n+1}^N-A(t_n)X_n^N = h\,B(t_n)X_{n+1}^N + \frac{h\,f(X_n^N)}{1+h\|f(X_n^N)\|} + g(X_n^N)\Delta W_n^N
\]
treats the linear drift part implicitly, the nonlinear drift part explicitly but tamed, and the diffusion explicitly. A dual tamed scheme in reduced \((u,\hat v)\)-variables is then used for analysis, and the continuous-time interpolation satisfies
\[
\left(\mathbb{E}\sup_{t\in[0,T]}\|X_t-\bar X_t^N\|^p\right)^{1/p} \le C h^{1/2},
\qquad p\in[1,\infty).
\]
The modification is therefore dictated by singular structure and algebraic constraints rather than by low regularity alone [2509.09032].

For additive-noise SDEs with drift that is not locally integrable near a singularity set \(S\), denominator taming is replaced by a geometric cut-off. Writing \(\rho(x)=\operatorname{dist}(x,S)\), the modified drift is
\[
b_n(x):=b(x)\,\mathbf 1_{\{x\in D\}\cap\{\rho(x)\ge \delta n^{-w}\}},
\qquad
w\in\Bigl[0,\frac{1}{3l}\Bigr].
\]
Thus the drift is turned off whenever the scheme approaches the singular set too closely. The resulting explicit Euler-type approximation converges strongly in \(L^p\) with rate \(1/2\),
\[
E\sup_{t\in[0,T]}|X_t^n-X_t|^p\le c\,n^{-p/2},
\]
and the classical Euler scheme is shown not even to possess finite first moments in this regime [2408.08118].

Further extensions show that the same design principle survives in more structured stochastic systems. For McKean–Vlasov equations with common noise, explicit tamed Euler and tamed Milstein schemes are built for the interacting particle system, with denominator taming on both drift and diffusion and with measure-derivative Milstein corrections; the strong orders are \(1/2\) and \(1\), respectively. For regime-switching SDEs, an explicit tamed Milstein-type scheme with an additional switching correction achieves strong \(L^2\) order \(1\). For neutral stochastic differential delay equations, tamed EM and tamed \(\theta\)-EM schemes provide strong convergence and, in some settings, exponential mean-square stability and almost sure exponential stability [2006.00463] [1909.07886] [1603.06747] [1612.02800] [1909.03155].

## 6. Ergodicity, sampling, and recurrent misconceptions

Modified taming is not limited to finite-time pathwise approximation. In the study of invariant measures and long-time dynamics, modified Euler recursions have been shown to be geometrically ergodic under a mixed probability distance and under weighted total variation distance. For the tamed Euler specialization,
\[
X_{(n+1)\delta}=X_{n\delta}+b^{(\delta)}(X_{n\delta})\delta+\sqrt{\delta}\,\xi_{n+1},
\]
one also obtains strict \(L^1\)-Wasserstein contractivity,
\[
W_1\bigl(\mathcal L_x(X_{n\delta}),\mathcal L_y(X_{n\delta})\bigr) \le C e^{-\lambda_3 n\delta}|x-y|,
\]
and the invariant measures satisfy
\[
W_1(\pi,\pi^\delta)\le C\delta^\theta.
\]
These results use refined basic coupling for the general modified framework and coupling by reflection for the tamed Euler scheme [2412.19377].

In Langevin sampling, the modified Tamed Unadjusted Langevin Algorithm replaces the drift by
\[
h_\lambda(\theta)=\frac{h(\theta)}{(1+\lambda|\theta|^{2r})^{1/2}},
\]
and the chain
\[
\theta_{n+1}^\lambda=\theta_n^\lambda-\lambda h_\lambda(\theta_n^\lambda)+\sqrt{2\lambda\beta^{-1}}\,\xi_{n+1}
\]
admits non-asymptotic convergence bounds in Wasserstein distance under non-convexity and super-linear growth. The discretization bias is \(\mathcal O(\lambda)\) in \(W_1\) and \(\mathcal O(\lambda^{1/2})\) in \(W_2\), with the analysis relying on dissipativity, Lyapunov drift estimates, and contractivity of a weighted semimetric rather than global strong convexity [2207.02600].

Two recurrent misconceptions are corrected by this body of work. First, taming does not mean only the classical denominator \(1+h^\alpha|b|\). In the cited literature it may mean a tamed Milstein derivative correction, a smooth far-field cut-off, an implicit-explicit split, a contractive state map, or a singular-set truncation. Second, taming does not inherently force a loss of the underlying discretization order. The adaptive Milstein construction achieves rate \((1+\alpha)/2\) under only \(C_l^{1+\alpha}\) regularity, the selective far-field modification retains weak order \(1\) for Euler, and several Milstein-type tamed methods recover order \(1\) in commutative, mean-field, or switching settings [2411.01849] [2507.09475] [1102.0662] [2006.00463] [1909.07886].

Taken together, these results define the modified tamed scheme as a flexible analytical and algorithmic template rather than a single method. Its essential feature is selective stabilization: the discretization is altered exactly where superlinearity, irregularity, constraint structure, or long-time sampling analysis requires intervention, and left unchanged where the original explicit method already has the correct local behavior.

Source: https://www.emergentmind.com/topics/modified-tamed-scheme