Tutorials › Stochastic Calculus › Auditing Ito Derivations with Explicit Moment Checks

Proof Strategies in Stochastic Calculus · Tutorial 574 of 600

Auditing Ito Derivations with Explicit Moment Checks

Learn to verify two common Ito-derived solutions by comparing their explicit formulas with independently calculated moments.

Intermediate 9 min read

What You'll Learn

  • Derive the explicit solutions of geometric Brownian motion and the Ornstein–Uhlenbeck process.
  • Compute their unconditional means and variances from Brownian motion’s normal law and the Ito isometry.
  • Check the geometric Brownian motion second moment using Ito’s formula and justify the stochastic-integral expectation.
  • Check the Ornstein–Uhlenbeck variance by applying Ito’s formula to the square of the process.
  • Detect the error caused by omitting Ito’s second-derivative correction.

What a Moment Check Can Detect

An explicit solution to a stochastic differential equation (SDE) can look plausible and still have a missing drift term, a wrong diffusion coefficient, or an incorrect Ito correction. A useful audit is to compute the solution’s mean and variance by a route independent of the proposed formula, then compare those moments with the SDE. This tutorial carries out that check for geometric Brownian motion and the Ornstein–Uhlenbeck process.

All calculations are under \(\mathbb{P}\) on a filtered probability space \((\Omega,\mathcal{F},(\mathcal{F}_t)_{t\ge0},\mathbb{P})\) carrying a standard Brownian motion \(B_t\). An Ito integral \(\int_0^t H_s\,dB_s\) below has a progressively measurable integrand \(H\) satisfying \(\mathbb{E}\int_0^t H_s^2\,ds<\infty\); this class is denoted by \(\mathcal{L}^2\). In that case the integral is a square-integrable martingale, and the Ito isometry gives \(\mathbb{E}[(\int_0^t H_s\,dB_s)^2]=\mathbb{E}\int_0^t H_s^2\,ds\). No risk-neutral measure \(\mathbb{Q}\) is used here.

We will also use the normal moment-generating formula. If \(B_t\sim N(0,t)\), then completing the square in its density gives \(\mathbb{E}[e^{aB_t}]=e^{a^2t/2}\) for every real \(a\). This provides a direct way to calculate moments of an exponential solution without taking expectations of an SDE.

Geometric Brownian Motion: The Independent Moment Calculation

Consider the SDE with initial value \(s>0\), drift \(\mu\in\mathbb{R}\), and diffusion coefficient \(\sigma\in\mathbb{R}\):

$$ dS_t=\mu S_t\,dt+\sigma S_t\,dB_t,\qquad S_0=s. $$

The drift coefficient is \(\mu S_t\), while the diffusion coefficient is \(\sigma S_t\). To obtain its candidate explicit solution, apply Ito’s formula to \(f(t,x)=\log x\) for \(x>0\). Its derivatives are \(f_t=0\), \(f_x=1/x\), and \(f_{xx}=-1/x^2\). The second-order term is therefore \(\tfrac12 f_{xx}(\sigma x)^2dt=-\tfrac12\sigma^2dt\), so \(d\log S_t=(\mu-\tfrac12\sigma^2)dt+\sigma\,dB_t\). Integrating and exponentiating gives the candidate below. Equivalently, applying Ito’s formula directly verifies it:

$$ S_t=s\exp\!\left((\mu-\tfrac12\sigma^2)t+\sigma B_t\right). $$

For that direct verification, set \(g(t,b)=s\exp((\mu-\tfrac12\sigma^2)t+\sigma b)\). Then \(g_t=(\mu-\tfrac12\sigma^2)g\), \(g_b=\sigma g\), and \(g_{bb}=\sigma^2g\). Ito’s formula gives \(dg=(g_t+\tfrac12g_{bb})dt+g_b\,dB_t=\mu g\,dt+\sigma g\,dB_t\), as required, and \(g(0,0)=s\). The coefficients are globally Lipschitz and of linear growth, so the standard existence-and-uniqueness theorem for SDEs ensures this solution agrees with the unique strong solution.

Theorem: For \(t\ge0\), the geometric Brownian motion solution has \(\mathbb{E}[S_t]=se^{\mu t}\), \(\mathbb{E}[S_t^2]=s^2e^{(2\mu+\sigma^2)t}\), and \(\operatorname{Var}(S_t)=s^2e^{2\mu t}(e^{\sigma^2t}-1)\).

Proof. The displayed explicit solution gives \(S_t^k=s^k\exp(k(\mu-\tfrac12\sigma^2)t+k\sigma B_t)\) for \(k=1,2\). Applying the normal moment-generating formula with \(a=k\sigma\) yields \(\mathbb{E}[S_t^k]=s^k\exp(k(\mu-\tfrac12\sigma^2)t+\tfrac12k^2\sigma^2t)\). Substituting \(k=1\) gives \(se^{\mu t}\), and substituting \(k=2\) gives \(s^2e^{(2\mu+\sigma^2)t}\). Subtracting the square of the mean from the second moment gives \(s^2e^{2\mu t}(e^{\sigma^2t}-1)\). These moments are finite, including the fourth moment, which will justify the Ito-integral check below. \(\square\)

Worked Example: Auditing Geometric Brownian Motion Moments

Take \(s=2\), \(\mu=\tfrac14\), \(\sigma=\tfrac12\), and \(t=2\). Compute the mean, second moment, and variance, and check the variance against the formula.

The exponent in the explicit solution is \((\tfrac14-\tfrac12\cdot\tfrac14)2+\tfrac12B_2=(\tfrac18)2+\tfrac12B_2=\tfrac14+\tfrac12B_2\). The mean from the theorem is \(2e^{\mu t}=2e^{(1/4)2}=2e^{1/2}\). The second moment is \(2^2e^{(2\mu+\sigma^2)t}=4e^{(2(1/4)+(1/2)^2)2}=4e^{(1/2+1/4)2}=4e^{3/2}\).

The variance obtained by subtraction is \(4e^{3/2}-(2e^{1/2})^2=4e^{3/2}-4e=4e(e^{1/2}-1)\). The direct variance formula gives \(s^2e^{2\mu t}(e^{\sigma^2t}-1)=4e^{2(1/4)2}(e^{(1/2)^2 2}-1)=4e(e^{1/2}-1)\). The two calculations agree. As a special case, setting \(\sigma=0\) in the formula gives variance \(0\), consistent with the deterministic solution \(S_t=se^{\mu t}\).

A Second Audit Through Ito’s Formula

There is another check for the geometric Brownian motion second moment. Apply Ito’s formula to \(f(x)=x^2\). Its derivatives are \(f'(x)=2x\) and \(f''(x)=2\); the correction term is \(\tfrac12 f''(S_t)(\sigma S_t)^2dt=\sigma^2S_t^2dt\). Thus

$$ d(S_t^2)=(2\mu+\sigma^2)S_t^2\,dt+2\sigma S_t^2\,dB_t. $$

Taking expectations is justified here, rather than automatic: the explicit formula has finite fourth moment, so \(\mathbb{E}\int_0^t(2\sigma S_s^2)^2ds=4\sigma^2\int_0^t\mathbb{E}[S_s^4]ds<\infty\). The stochastic integral is consequently a true square-integrable martingale with mean zero. Taking expectations in the integrated identity gives \(m_2(t)=s^2+(2\mu+\sigma^2)\int_0^t m_2(u)\,du\), where \(m_2(t)=\mathbb{E}[S_t^2]\). The explicit moment already derived satisfies this equation: its derivative is \((2\mu+\sigma^2)s^2e^{(2\mu+\sigma^2)t}\) and its initial value is \(s^2\). This agreement audits both the Ito correction and the second-moment exponent.

The equation \(dB_t\,dB_t=dt\) accounts for the positive \(\sigma^2S_t^2dt\) term. The other multiplication rules are \(dB_t\,dt=0\) and \(dt\,dt=0\); they do not remove the Brownian quadratic-variation term. The equality \(dB_t\,dB_t=dt\) is an Ito-calculus rule, not an ordinary chain rule for differentiable paths.

Ornstein–Uhlenbeck: Mean and Variance from the Integral

Consider an Ornstein–Uhlenbeck process with initial value \(x\), mean-reversion rate \(\lambda>0\), and constant diffusion coefficient \(\eta\in\mathbb{R}\):

$$ dX_t=-\lambda X_t\,dt+\eta\,dB_t,\qquad X_0=x. $$

Here the drift is \(-\lambda X_t\) and the diffusion coefficient is \(\eta\). Multiplying by the deterministic integrating factor \(e^{\lambda t}\), using the Ito product rule, and integrating gives the candidate solution

$$ X_t=xe^{-\lambda t}+\eta\int_0^t e^{-\lambda(t-s)}\,dB_s. $$

For fixed \(t\), the integrand is deterministic and \(\int_0^t e^{-2\lambda(t-s)}ds=(1-e^{-2\lambda t})/(2\lambda)<\infty\), so it belongs to \(\mathcal{L}^2\). The integral has mean zero and, by the Ito isometry, variance \(\int_0^t e^{-2\lambda(t-s)}ds\). Therefore

Theorem: For \(\lambda>0\), the Ornstein–Uhlenbeck solution has \(\mathbb{E}[X_t]=xe^{-\lambda t}\) and \(\operatorname{Var}(X_t)=\dfrac{\eta^2}{2\lambda}(1-e^{-2\lambda t})\).

Proof. The stochastic integral in the explicit solution is a square-integrable martingale at its upper time limit, so its expectation is zero. Its covariance with the deterministic constant \(xe^{-\lambda t}\) is zero. Applying the Ito isometry to its variance gives \(\eta^2\int_0^t e^{-2\lambda(t-s)}ds\). Evaluating the integral yields \(\eta^2(1-e^{-2\lambda t})/(2\lambda)\), proving both formulas. \(\square\)

Worked Example: An Exact Ornstein–Uhlenbeck Check

Let \(x=3\), \(\lambda=1\), \(\eta=2\), and \(t=\ln 2\). Find the mean and variance, and check them using Ito’s formula for \(X_t^2\).

Since \(e^{-\lambda t}=e^{-\ln2}=\tfrac12\), the mean is \(3\cdot\tfrac12=\tfrac32\). Also \(e^{-2\lambda t}=e^{-2\ln2}=\tfrac14\). The variance from the integral formula is \(\dfrac{2^2}{2\cdot1}(1-\tfrac14)=2\cdot\tfrac34=\tfrac32\). Hence the second moment is \((\tfrac32)^2+\tfrac32=\tfrac94+\tfrac64=\tfrac{15}{4}\).

For the independent check, apply Ito’s formula to \(f(x)=x^2\), whose derivatives are \(f'(x)=2x\) and \(f''(x)=2\). The correction term is \(\tfrac12(2)\eta^2dt=4dt\), giving \(d(X_t^2)=(-2X_t^2+4)dt+4X_t\,dB_t\). The last integral has mean zero: the explicit second moment is finite on the bounded time interval \([0,\ln2]\), so \(\mathbb{E}\int_0^{\ln2}(4X_s)^2ds<\infty\). Thus \(m_2'(t)=-2m_2(t)+4\), \(m_2(0)=9\). Its solution is \(m_2(t)=9e^{-2t}+2(1-e^{-2t})\); at \(t=\ln2\), this is \(9/4+2(3/4)=15/4\), matching the independently computed second moment.

Common Mistake: Dropping the Ito Correction

Worked Example: A Missing Term Changes the Mean

Set \(s=1\), \(\mu=0\), \(\sigma=1\), and \(t=1\) in the geometric Brownian motion SDE. A mistaken ordinary-chain-rule calculation might write \(d\log S_t=dB_t\), giving the incorrect candidate \(\widetilde S_t=e^{B_t}\).

For the actual solution, the derivatives of \(g(t,b)=e^{-t/2+b}\) are \(g_t=-g/2\), \(g_b=g\), and \(g_{bb}=g\). The Ito correction is \(\tfrac12g_{bb}dt=\tfrac12gdt\), so \(dg=(-g/2+g/2)dt+g\,dB_t=g\,dB_t\), as the SDE with zero drift requires. Its mean at \(t=1\) is \(e^{0}=1\), and its variance is \(e^0(e^1-1)=e-1\).

The incorrect candidate has mean \(\mathbb{E}[e^{B_1}]=e^{1/2}\), not \(1\), and variance \(\mathbb{E}[e^{2B_1}]-(\mathbb{E}[e^{B_1}])^2=e^2-e\). It also fails the SDE: applying Ito’s formula to \(e^{B_t}\), with derivatives \(f'=f\) and \(f''=f\), gives \(d(e^{B_t})=\tfrac12e^{B_t}dt+e^{B_t}dB_t\), which contains a drift absent from the proposed equation. The moment discrepancy exposes the same error as direct substitution.

How to Use the Audit

For an Ito-derived explicit formula, separate two questions. First, does Ito’s formula verify that the candidate satisfies the SDE, including the second-derivative correction? Second, do independent calculations of its mean and variance agree with the drift and diffusion in that SDE? The checks are complementary: substitution detects a pathwise algebra error, while moments can reveal a wrong exponential factor or a stochastic integral whose expectation has been assumed to vanish without adequate integrability.

Do not confuse a local martingale with a true martingale. A stochastic integral with an \(\mathcal{L}^2\) integrand is a true square-integrable martingale, so its expectation and conditional expectation increments are zero. Without that integrability (or another sufficient condition), an Ito integral may only be a local martingale, and taking expectations may be invalid. In the checks above, the needed square-integrability was verified explicitly.

Check Your Understanding

For each question, state the integrability or correction-term justification needed for your calculation.

  1. For geometric Brownian motion with \(s=1\), \(\mu=0\), \(\sigma=2\), and \(t=1\), compute the mean and variance from the moment formulas.
  2. Apply Ito’s formula to \(X_t^2\) when \(dX_t=-\lambda X_tdt+\eta dB_t\). Identify the drift correction and the stochastic-integral integrand.
  3. For the Ornstein–Uhlenbeck process, take the limit of the variance formula as \(\lambda\) decreases to zero, with \(\eta\) and \(t\) fixed.
  4. Why does finite \(\mathbb{E}\int_0^t H_s^2ds\) justify taking the expectation of \(\int_0^tH_s\,dB_s\)?
  5. For \(dS_t=\sigma S_t\,dB_t\), explain why the candidate \(S_t=se^{\sigma B_t}\) does not satisfy the equation when \(\sigma\ne0\).