Skip to main content

Chapter 4 — SDE Definition and Existence

Guiding question

Given a drift function b(t,x)b(t,x) and a diffusion function σ(t,x)\sigma(t,x), does the stochastic differential equation dXt=b(t,Xt) dt+σ(t,Xt) dWtdX_t = b(t,X_t)\,dt + \sigma(t,X_t)\,dW_t have a solution? Is it unique? What conditions on bb and σ\sigma are needed?

An SDE is a prescription for how a process evolves in time: a deterministic drift plus a random diffusion. The existence question asks whether such a prescription is self-consistent — whether a process with the prescribed dynamics can actually be constructed. Uniqueness asks whether the dynamics determine the process completely.


4.1 The SDE as an integral equation​

The shorthand dXt=b(t,Xt) dt+σ(t,Xt) dWtdX_t = b(t,X_t)\,dt + \sigma(t,X_t)\,dW_t is notation for the integral equation:

Xt=X0+∫0tb(s,Xs) ds+∫0tσ(s,Xs) dWs,0≤t≤T.X_t = X_0 + \int_0^t b(s,X_s)\,ds + \int_0^t \sigma(s,X_s)\,dW_s, \qquad 0\le t\le T.

The first integral is a Lebesgue integral (pathwise, since bb is a drift), and the second is an Itô integral. Finding XX means finding a continuous adapted process that satisfies this identity for all tt.

The SDE is implicit

The unknown process XX appears inside the integrals on the right side. This makes the SDE an implicit equation. The existence theorem constructs XX as the fixed point of a contraction map in the L2L^2 sense.


4.2 Strong and weak solutions​

There are two notions of solution.

A strong solution is a process XX defined on the given probability space (Ω,F,P)(\Omega,\mathcal F,\mathbb P), adapted to the filtration generated by WW, that satisfies the integral equation pathwise.

A weak solution is a triple (Ω,F,P,W,X)(\Omega,\mathcal F,\mathbb P, W, X) — possibly a new probability space with a new Brownian motion — on which the integral equation holds.

Every strong solution is a weak solution. The converse is false in general.

Analogy with ODEs

For an ODE x˙=f(x)\dot x = f(x), a strong solution constructs the trajectory from the given vector field ff and starting point. A weak solution allows changing the underlying space (redefining what "random" means) as long as the marginal dynamics are preserved.

In finance, the distinction matters for hedging: a strong solution gives a trading strategy that replicates under the original measure, while a weak solution only gives pricing without a hedging prescription.


4.3 Lipschitz conditions​

The standard hypothesis for existence and uniqueness of strong solutions is:

Standard conditions (Lipschitz and linear growth)

Suppose b,σ:[0,T]×R→Rb,\sigma:[0,T]\times\mathbb R\to\mathbb R satisfy:

  1. Lipschitz in space: there exists K>0K>0 such that for all t∈[0,T]t\in[0,T] and x,y∈Rx,y\in\mathbb R,
∣b(t,x)−b(t,y)∣+∣σ(t,x)−σ(t,y)∣≤K∣x−y∣.|b(t,x)-b(t,y)| + |\sigma(t,x)-\sigma(t,y)| \le K|x-y|.
  1. Linear growth: there exists C>0C>0 such that
∣b(t,x)∣2+∣σ(t,x)∣2≤C(1+x2).|b(t,x)|^2 + |\sigma(t,x)|^2 \le C(1+x^2).

Under these conditions and any F0\mathcal F_0-measurable X0X_0 with E[X02]<∞\mathbb E[X_0^2]<\infty, the SDE has a unique strong solution satisfying E[sup⁡t≤TXt2]<∞\mathbb E[\sup_{t\le T}X_t^2]<\infty.


4.4 Picard iteration: constructing the solution​

The proof uses the same idea as the classical Picard–Lindelöf theorem for ODEs.

Define the sequence Xt(0)=X0X^{(0)}_t = X_0 (constant), and

Xt(n+1)=X0+∫0tb(s,Xs(n)) ds+∫0tσ(s,Xs(n)) dWs.X^{(n+1)}_t = X_0 + \int_0^t b(s,X^{(n)}_s)\,ds + \int_0^t \sigma(s,X^{(n)}_s)\,dW_s.

Each X(n+1)X^{(n+1)} is well-defined because X(n)X^{(n)} is adapted and square-integrable. We show the sequence converges.

Convergence of Picard iterates

Let Zt(n)=Xt(n+1)−Xt(n)Z^{(n)}_t = X^{(n+1)}_t - X^{(n)}_t and εn(t)=E[sup⁡s≤t∣Zs(n)∣2]\varepsilon_n(t) = \mathbb E[\sup_{s\le t}|Z^{(n)}_s|^2].

By the Lipschitz condition and the Burkholder–Davis–Gundy inequality (Section 2.9):

εn+1(t)≤C∫0tεn(s) ds.\varepsilon_{n+1}(t) \le C\int_0^t \varepsilon_n(s)\,ds.

Starting from ε0(t)=E[sup⁡s≤t∣Xs(1)−X0∣2]≤C0t\varepsilon_0(t) = \mathbb E[\sup_{s\le t}|X^{(1)}_s - X_0|^2] \le C_0 t (from linear growth), iterating:

εn(t)≤Cntnn!C0.\varepsilon_n(t) \le \frac{C^n t^n}{n!}C_0.

Therefore ∑nεn(T)1/2<∞\sum_n\varepsilon_n(T)^{1/2} < \infty, and the Borel–Cantelli argument gives almost sure and L2L^2 convergence.

Picard iteration converging to OU solution dX = −X dt + dW
X⁽⁰⁾ = x₀X⁽1⁾X⁽2⁾X⁽3⁾True solution0T=1
Picard iteration builds X⁽ⁿ⁺¹⁾ by plugging X⁽ⁿ⁾ into the drift integral. Starting from the constant x₀, each iterate adds one more term of the exponential expansion of e^−t. Under Lipschitz conditions the iterates converge in L² — this is the constructive proof of existence.

4.5 Uniqueness via Gronwall's inequality​

Pathwise uniqueness

Under the Lipschitz condition, if XX and YY both solve the SDE with the same initial condition and the same Brownian motion, then Xt=YtX_t = Y_t almost surely for all tt.

Proof

Let Zt=Xt−YtZ_t = X_t - Y_t. Then

Zt=∫0t[b(s,Xs)−b(s,Ys)] ds+∫0t[σ(s,Xs)−σ(s,Ys)] dWs.Z_t = \int_0^t [b(s,X_s)-b(s,Y_s)]\,ds + \int_0^t [\sigma(s,X_s)-\sigma(s,Y_s)]\,dW_s.

By Lipschitz and the Itô isometry/BDG inequality:

E ⁣[sup⁡s≤tZs2]≤C∫0tE[Zs2] ds≤C∫0tE ⁣[sup⁡r≤sZr2]ds.\mathbb E\!\left[\sup_{s\le t}Z_s^2\right] \le C\int_0^t \mathbb E[Z_s^2]\,ds \le C\int_0^t \mathbb E\!\left[\sup_{r\le s}Z_r^2\right]ds.

Let u(t)=E[sup⁡s≤tZs2]u(t) = \mathbb E[\sup_{s\le t}Z_s^2]. Then u(t)≤C∫0tu(s) dsu(t)\le C\int_0^t u(s)\,ds with u(0)=0u(0)=0.

Gronwall's inequality: If u(t)≤A+C∫0tu(s) dsu(t)\le A + C\int_0^t u(s)\,ds with A=0A=0, then u(t)≤0u(t)\le 0. Therefore u≡0u\equiv 0, so Zt=0Z_t = 0 a.s.

Gronwall as an exponential estimate

Gronwall's inequality says: if a non-negative function grows no faster than a linear fraction of its own past accumulation, then it must be zero (when starting from zero). It is the integral analogue of the comparison principle for ODEs: if u˙≤Cu\dot u \le Cu, then u(t)≤u(0)eCtu(t)\le u(0)e^{Ct}. When u(0)=0u(0)=0, u≡0u\equiv 0.


4.6 The Markov property of SDE solutions​

Under the standard conditions, the strong solution of

dXt=b(t,Xt) dt+σ(t,Xt) dWtdX_t = b(t,X_t)\,dt + \sigma(t,X_t)\,dW_t

inherits the Markov property from Brownian motion: the future evolution depends on the past only through the current value XtX_t.

More precisely, for s≤ts\le t and measurable ff:

E[f(Xt)∣Fs]=E[f(Xt)∣Xs].\mathbb E[f(X_t)\mid\mathcal F_s] = \mathbb E[f(X_t)\mid X_s].

This follows from pathwise uniqueness and the fact that the increments Wt−WsW_t - W_s are independent of Fs\mathcal F_s.


4.7 Beyond Lipschitz: non-standard cases​

The Lipschitz condition is sufficient but not necessary. Several important models fail it but still have unique solutions:

ModelIssueResolution
CIR: dX=κ(θ−X)dt+σXdWdX = \kappa(\theta-X)dt + \sigma\sqrt{X}dWσ(x)=σx\sigma(x) = \sigma\sqrt{x} is not Lipschitz at x=0x=0Yamada–Watanabe uniqueness theorem
GBM: dS=μS dt+σS dWdS = \mu S\,dt + \sigma S\,dWLipschitz fails globally but holds locallyLocal Lipschitz + linear growth suffice
Absolute value: $dX =X^\alpha dW$

4.8 A worked example: verification​

Verify that Xt=x0e(μ−σ2/2)t+σWtX_t = x_0 e^{(\mu-\sigma^2/2)t + \sigma W_t} solves dX=μX dt+σX dWdX = \mu X\,dt + \sigma X\,dW.

Apply Itô's formula to f(t,x)=x0e(μ−σ2/2)t+σxf(t,x) = x_0 e^{(\mu-\sigma^2/2)t+\sigma x}:

∂tf=(μ−σ22)f,∂xf=σf,∂xxf=σ2f.\partial_t f = \left(\mu-\tfrac{\sigma^2}{2}\right)f,\quad \partial_x f = \sigma f,\quad \partial_{xx}f = \sigma^2 f.

Therefore

dXt=[(μ−σ22)f+12σ2f]dt+σf dWt=μXt dt+σXt dWt.dX_t = \left[\left(\mu-\tfrac{\sigma^2}{2}\right)f + \tfrac12\sigma^2 f\right]dt + \sigma f\,dW_t = \mu X_t\,dt + \sigma X_t\,dW_t.

The SDE is satisfied. Combined with uniqueness, this is the unique strong solution.


4.9 Problems​

Problem 1: Lipschitz check

Check whether b(x)=x2b(x) = x^2 satisfies the Lipschitz condition on R\mathbb R. Does it satisfy linear growth? What does this mean for global existence of dX=X2 dt+dWdX = X^2\,dt + dW?

Problem 2: Gronwall's lemma

Prove Gronwall's lemma: if u:[0,T]→[0,∞)u:[0,T]\to[0,\infty) is continuous and u(t)≤A+C∫0tu(s) dsu(t)\le A + C\int_0^t u(s)\,ds, then u(t)≤A eCtu(t)\le A\,e^{Ct}. (Hint: define v(t)=e−Ct(A+C∫0tu(s) ds)v(t) = e^{-Ct}(A+C\int_0^t u(s)\,ds) and differentiate.)

Problem 3: Additive noise SDE

Find the strong solution to dXt=−αXt dt+σ dWtdX_t = -\alpha X_t\,dt + \sigma\,dW_t with X0=x0X_0 = x_0. (Hint: use the integrating factor eαte^{\alpha t} to reduce to an exact differential.)

Problem 4: Pathwise uniqueness for the OU process

Verify directly (without citing the general theorem) that any two solutions of dX=−αX dt+dWdX = -\alpha X\,dt + dW with the same WW and same X0X_0 are indistinguishable, by writing the difference equation and applying Gronwall.


4.10 Selected solutions​

Solution 1

∣b(x)−b(y)∣=∣x2−y2∣=∣x+y∣⋅∣x−y∣|b(x)-b(y)| = |x^2-y^2| = |x+y|\cdot|x-y|. For large ∣x∣|x| or ∣y∣|y|, the factor ∣x+y∣|x+y| is unbounded, so Lipschitz fails globally. Also ∣b(x)∣2=x4|b(x)|^2 = x^4 grows faster than C(1+x2)C(1+x^2), so linear growth fails. The SDE dX=X2 dt+dWdX = X^2\,dt + dW has solutions that may explode (reach ±∞\pm\infty) in finite time; global existence is not guaranteed.

Solution 2

Let v(t)=e−Ct[A+C∫0tu(s) ds]v(t) = e^{-Ct}[A + C\int_0^t u(s)\,ds]. Then v(0)=Av(0) = A and

v′(t)=−Ce−Ct[A+C∫0tu]+Ce−Ctu(t)=Ce−Ct[u(t)−(A+C∫0tu)]≤0,v'(t) = -Ce^{-Ct}\left[A+C\int_0^t u\right] + Ce^{-Ct}u(t) = Ce^{-Ct}[u(t)-(A+C\int_0^t u)] \le 0,

since u(t)≤A+C∫0tuu(t)\le A + C\int_0^t u. So vv is decreasing, hence v(t)≤v(0)=Av(t)\le v(0) = A, giving A+C∫0tu≤AeCtA + C\int_0^t u\le Ae^{Ct}. Since u(t)≤A+C∫0tu≤AeCtu(t)\le A + C\int_0^t u\le Ae^{Ct}, Gronwall is proved.

Solution 3

Let Yt=eαtXtY_t = e^{\alpha t}X_t. By Itô's product rule:

dYt=αeαtXt dt+eαtdXt=αeαtXt dt+eαt(−αXt dt+σ dWt)=σeαt dWt.dY_t = \alpha e^{\alpha t}X_t\,dt + e^{\alpha t}dX_t = \alpha e^{\alpha t}X_t\,dt + e^{\alpha t}(-\alpha X_t\,dt + \sigma\,dW_t) = \sigma e^{\alpha t}\,dW_t.

Integrating: Yt=Y0+σ∫0teαs dWs=x0+σ∫0teαs dWsY_t = Y_0 + \sigma\int_0^t e^{\alpha s}\,dW_s = x_0 + \sigma\int_0^t e^{\alpha s}\,dW_s. Therefore

Xt=e−αtx0+σ∫0te−α(t−s) dWs.X_t = e^{-\alpha t}x_0 + \sigma\int_0^t e^{-\alpha(t-s)}\,dW_s.

This is the Ornstein–Uhlenbeck process. It is Gaussian, mean-reverting to zero with rate α\alpha.


4.11 What this chapter gives you​

ConceptContent
SDE as integral equationXt=X0+∫b ds+∫σ dWsX_t = X_0 + \int b\,ds + \int\sigma\,dW_s
Strong solutionadapted to filtration of WW, pathwise
Lipschitz + linear growthsufficient conditions for existence and uniqueness
Picard iterationconstructive proof via L2L^2 contraction
Gronwall's inequalityuniqueness via integral comparison
Markov propertySDE solutions inherit Brownian Markov structure

The next chapter solves a collection of important SDEs explicitly using Itô's formula and integrating factors.