Tutorials › Stochastic Calculus › Explosion Tests for Superlinear Drift Diffusions

Explosion, Comparison Theorems and Boundary Classification · Tutorial 289 of 600

Explosion Tests for Superlinear Drift Diffusions

Compute the scale and speed endpoint tests for cubic drift diffusions, and see why an outward cubic drift can explode while an inward one cannot.

Advanced 12 min read

What You'll Learn

  • Define the scale density, scale function, and speed density for constant-diffusion scalar SDEs.
  • Apply Feller’s endpoint test to the outward cubic drift.
  • Bound the scale-tail integral that determines whether the outward endpoint is accessible.
  • Show why the inward cubic drift has infinite scale at both endpoints.
  • Verify nonexplosion for the inward drift with a stopped Lyapunov estimate.
  • Identify the sign error that reverses the scale test’s conclusion.

Two Cubic Drifts, Two Different Outcomes

A drift that grows cubically can dominate Brownian fluctuations at large distances, but its sign matters. The drift \(x^3\) points away from the origin on both sides, whereas \(-x^3\) points toward it. Feller’s explosion test makes this distinction quantitative through the scale and speed densities. For the outward drift, a tail integral at each endpoint is finite; for the inward drift, the scale function diverges at both endpoints.

Work on a filtered probability space \((\Omega,\mathcal{F},(\mathcal{F}_t)_{t\ge0},\mathbb{P})\) carrying a standard Brownian motion \(B_t\). For either equation, the diffusion coefficient is the constant \(1\), and the drift is locally Lipschitz. Thus, for each initial value \(X_0=x\), there is a pathwise unique strong solution up to its lifetime. Concretely, before explosion it satisfies \(X_t=x+\int_0^t b(X_s)\,ds+\int_0^t 1\,dB_s\). The integrand \(1\) is progressively measurable and belongs to \(\mathcal{L}^2\) on every finite interval, since \(\mathbb{E}\int_0^T1^2\,ds=T\). Local Lipschitz continuity does not ensure a global solution: define \(\tau_n=\inf\{t\ge0:|X_t|\ge n\}\) and \(\tau_\infty=\lim_{n\to\infty}\tau_n\), the lifetime at which explosion may occur.

Definition: For \(dX_t=b(X_t)\,dt+dB_t\), define the scale density \(p(x)=\exp\!\left(-\int_0^x 2b(u)\,du\right)\), the scale function \(s(x)=\int_0^x p(u)\,du\), and the speed density \(m(x)=2/p(x)\). If \(s(+\infty)<\infty\), the Feller integral for the right endpoint, with reference point \(0\), is \(I_+=\int_0^\infty (s(+\infty)-s(y))m(y)\,dy\). The corresponding left-endpoint integral is \(I_-=\int_{-\infty}^0(s(y)-s(-\infty))m(y)\,dy\), when \(s(-\infty)>-\infty\).

We use the standard Feller explosion test from the earlier material on finite-time explosion: an endpoint with finite scale endpoint and finite corresponding integral is accessible as an exit endpoint, so there is positive probability of reaching it in finite time and hence of explosion there. If the scale endpoint is infinite, it cannot be reached in finite time. The test is an endpoint criterion; it does not say that every sample path explodes, nor does it give a deterministic explosion time.

Outward Cubic Drift: \(dX_t=X_t^3\,dt+dB_t\)

For \(b(x)=x^3\), integration gives \(\int_0^x2u^3\,du=x^4/2\). Therefore the scale density is \(p(x)=e^{-x^4/2}\), and the speed density is \(m(x)=2e^{x^4/2}\). Since \(p\) is even and decays rapidly, both scale endpoints are finite. The key point is to check that the scale tail multiplied by the speed density is integrable.

Worked Example: The Right Endpoint for Outward Cubic Drift

Compute and bound \(I_+\) for \(b(x)=x^3\), showing each step that establishes its finiteness.

For \(y\ge1\), the difference of fourth powers satisfies \(z^4-y^4\ge4y^3(z-y)\) whenever \(z\ge y\), because the derivative of \(z^4\) is \(4z^3\ge4y^3\) on that range. Thus

$$ \begin{aligned} s(+\infty)-s(y) &=\int_y^\infty e^{-z^4/2}\,dz\\ &\le e^{-y^4/2}\int_y^\infty e^{-2y^3(z-y)}\,dz\\ &=e^{-y^4/2}\int_0^\infty e^{-2y^3v}\,dv =\frac{e^{-y^4/2}}{2y^3}. \end{aligned} $$

Multiplying by \(m(y)=2e^{y^4/2}\) gives \((s(+\infty)-s(y))m(y)\le1/y^3\). Consequently, the part of \(I_+\) from \(1\) to infinity is bounded by \(\int_1^\infty y^{-3}\,dy=1/2\). On \([0,1]\), the integrand is continuous and finite, so its integral is finite as well. Hence \(I_+<\infty\).

Check: the tail bound integrates to \(\int_1^\infty y^{-3}\,dy=[-1/(2y^2)]_1^\infty=1/2\), a finite value. Feller’s test therefore classifies \(+\infty\) as an accessible exit endpoint for this diffusion.

The coefficients and the scale density are symmetric under \(x\mapsto -x\). In particular, changing variables \(y=-v\) in the left-endpoint integral gives the same tail bound as on the right. Both \(I_-\) and \(I_+\) are finite. Feller’s test therefore permits finite-time explosion through either endpoint, with positive probability from a finite starting point. The symmetry of the coefficients does not mean that a given path explodes in both directions; the two endpoints describe alternative ways the lifetime can end.

Result: For \(dX_t=X_t^3\,dt+dB_t\), the scale function has finite limits at both endpoints and both Feller endpoint integrals are finite. Thus both endpoints are accessible exits, and finite-time explosion has positive probability.

Justification. The calculation above establishes finiteness of the right integral and, by symmetry, of the left integral. The scale density is positive and integrable on each tail by the same exponential estimate, so both scale limits are finite. The stated conclusion then follows from Feller’s explosion test. \(\square\)

Inward Cubic Drift: \(dX_t=-X_t^3\,dt+dB_t\)

For \(b(x)=-x^3\), the sign in the scale exponent reverses: \(p(x)=e^{x^4/2}\), and \(m(x)=2e^{-x^4/2}\). Because \(p(x)\ge1\), its integral over either unbounded half-line diverges. Thus \(s(+\infty)=+\infty\) and \(s(-\infty)=-\infty\). The scale test already rules out finite-time explosion at either endpoint. A separate stopped estimate confirms the conclusion without relying on endpoint classification alone.

Worked Example: Infinite Scale for Inward Cubic Drift

Verify the scale divergence at each endpoint and contrast it with the outward case.

On the right, \(e^{x^4/2}\ge1\) for \(x\ge0\), so for every \(R>0\),

$$ s(R)=\int_0^R e^{x^4/2}\,dx\ge\int_0^R1\,dx=R. $$

As \(R\to\infty\), this lower bound tends to infinity, so \(s(+\infty)=+\infty\). On the left, evenness of \(p\) gives \(s(-R)=-\int_0^R e^{x^4/2}\,dx\le-R\), so \(s(-\infty)=-\infty\). The endpoint test uses infinite scale here; no finite Feller exit integral needs to be calculated.

Check: replacing \(b(x)=-x^3\) by \(b(x)=x^3\) changes the exponent from \(+x^4/2\) to \(-x^4/2\). In the outward case the scale density is integrable on the tails; in the inward case it is bounded below by \(1\), so its tail integral diverges.

Proposition: The solution of \(dX_t=-X_t^3\,dt+dB_t\) is nonexplosive: for every finite \(T\), \(\mathbb{P}(\tau_\infty\le T)=0\).

Proof. Apply Itô’s formula, with \(f(x)=x^2\), \(f'(x)=2x\), and \(f''(x)=2\), to the solution stopped at \(\tau_n\). The drift is \(-x^3\) and the diffusion coefficient is \(1\), so the correction term is \(\tfrac12 f''(x)\,dt=dt\). Therefore

$$ d(X_{t\wedge\tau_n})^2 =\mathbf{1}_{\{t<\tau_n\}}(1-2X_t^4)\,dt +2\mathbf{1}_{\{t<\tau_n\}}X_t\,dB_t. $$

Before \(\tau_n\), \(|X_t|<n\), so the stopped stochastic integrand is bounded on \([0,T]\) and belongs to \(\mathcal{L}^2\). Its Itô integral is a true martingale: it is square-integrable by the Itô isometry, and its conditional expectation of each future increment given the current filtration is zero. Taking expectations, and using \(1-2X_t^4\le1\), yields \(\mathbb{E}[X_{T\wedge\tau_n}^2]\le x^2+T\). On \(\{\tau_n\le T\}\), continuity of the stopped path gives \(|X_{T\wedge\tau_n}|=n\). Hence

$$ n^2\mathbb{P}(\tau_n\le T) \le \mathbb{E}[X_{T\wedge\tau_n}^2] \le x^2+T, \qquad \mathbb{P}(\tau_n\le T)\le\frac{x^2+T}{n^2}. $$

The stopping times \(\tau_n\) increase with \(n\), and \(\{\tau_\infty\le T\}\subseteq\{\tau_n\le T\}\) for every \(n\). Thus \(\mathbb{P}(\tau_\infty\le T)\le(x^2+T)/n^2\) for every \(n\), and letting \(n\to\infty\) gives \(\mathbb{P}(\tau_\infty\le T)=0\). This holds for each finite \(T\), proving nonexplosion. \(\square\)

Why the Scale Sign Matters

Worked Example: A Positive Cubic Coefficient

Let \(a>0\) and consider \(dX_t=aX_t^3\,dt+dB_t\). Recompute the scale and tail bound to see how the outward-drift calculation depends on \(a\).

Here \(\int_0^x2au^3\,du=ax^4/2\), so \(p(x)=e^{-ax^4/2}\) and \(m(x)=2e^{ax^4/2}\). For \(z\ge y\ge1\), \(z^4-y^4\ge4y^3(z-y)\), giving

$$ \begin{aligned} s(+\infty)-s(y) &=\int_y^\infty e^{-az^4/2}\,dz\\ &\le e^{-ay^4/2}\int_0^\infty e^{-2ay^3v}\,dv =\frac{e^{-ay^4/2}}{2ay^3}. \end{aligned} $$

Multiplication by \(m(y)=2e^{ay^4/2}\) bounds the Feller integrand by \(1/(ay^3)\). Its integral from \(1\) to infinity is \(1/(2a)\), finite for every \(a>0\). The compact part is finite by continuity, and symmetry gives the same conclusion at the left endpoint. Therefore every positive cubic coefficient gives finite endpoint integrals.

Check: substituting \(a=1\) gives the earlier tail bound \(1/y^3\) and tail integral \(1/2\). The factor \(1/(2a)\) comes directly from \(\int_1^\infty (ay^3)^{-1}\,dy=1/(2a)\).

A common mistake is to omit the minus sign in the definition of the scale density. For \(b(x)=x^3\), the incorrect calculation \(p(x)=\exp(\int_0^x2u^3\,du)=e^{x^4/2}\) would imply infinite scale at both endpoints and would wrongly suggest nonexplosion. The correct calculation is \(p(x)=\exp(-\int_0^x2u^3\,du)=e^{-x^4/2}\). Reversing the drift to \(-x^3\) reverses this exponent. The scale calculation is therefore sensitive to the sign of the drift, not just to its growth rate.

These tests are useful because local existence and pathwise uniqueness do not settle whether a solution lasts for all time. Both cubic equations have locally Lipschitz coefficients, yet the outward drift admits explosion and the inward drift does not. The Feller calculation identifies the endpoint behavior; the Lyapunov estimate supplies a direct independent check for the nonexplosive case. Neither calculation asserts uniqueness in law as a substitute for pathwise uniqueness: here the local strong solution is pathwise unique up to its lifetime, and the stopped estimate rules out a finite lifetime for the inward equation.

Check Your Understanding

Use the scale and speed definitions above, and distinguish an endpoint’s scale behavior from the finiteness of its Feller integral.

  1. For \(b(x)=2x^3\), compute \(p(x)\) and \(m(x)\), then state a bound for the right-tail Feller integrand for \(y\ge1\).
  2. Why does \(p(x)=e^{x^4/2}\) imply that the scale function diverges at both endpoints?
  3. In the stopped Itô estimate for \(X_t^2\) with drift \(-X_t^3\), identify the second-derivative correction term and explain why the stopped stochastic integral has mean zero.
  4. Does finite Feller integral at \(+\infty\) imply that every path explodes through \(+\infty\)? Explain what the endpoint test concludes.
  5. Find the sign error in the calculation \(p(x)=\exp(\int_0^x2u^3\,du)\) for drift \(b(x)=x^3\), and give the corrected scale density.