←Return

Stochastic Calculus: Conceptual Theory, FAQs, and Quantitative Finance Interview Problems

In quantitative finance, theoretical physics, and academic screens, stochastic calculus is both the ultimate language of continuous-time modeling and the primary filter in technical interviews. Questions on stochastic calculus generally split into two distinct categories:

  1. →Conceptual Theory: Why classical calculus fails, how quadratic variation fundamentally changes integration, the difference between Itô and Stratonovich, and the machinery of measure changes.
  2. →Practical Problem-Solving & Interview Puzzles: Direct calculations involving Itô's Lemma, solving Stochastic Differential Equations (SDEs), and evaluating stochastic integrals and expectations under pressure.

This guide provides a rigorous yet intuitive reference covering both dimensions.


Part I: Conceptual Theory & Foundational FAQs

1. Why does standard Newtonian calculus fail for Brownian motion?

The Core Reason: Standard Riemann-Stieltjes integration ∫f(t)dg(t)\int f(t) dg(t) requires the integrator g(t)g(t) to have bounded (finite) total variation.

A standard Brownian motion (or Wiener process) WtW_t:

  • →Is almost surely continuous everywhere, but
  • →Is almost surely nowhere differentiable, and
  • →Has unbounded (infinite) total variation over any finite time interval [0,t][0, t]: lim⁡Δt→0∑i=1n∣Wti−Wti−1∣=∞a.s.\lim_{\Delta t \to 0} \sum_{i=1}^n |W_{t_i} - W_{t_{i-1}}| = \infty \quad \text{a.s.}

Because Brownian motion oscillates infinitely wildly at infinitesimal scales, the ordinary chain rule: dfdt=f′(x)dxdt\frac{df}{dt} = f'(x) \frac{dx}{dt} cannot be evaluated because dWtdt\frac{dW_t}{dt} does not exist in the classical sense.

2. What is Quadratic Variation, and why is (dWt)2=dt(dW_t)^2 = dt?

While the total variation of Brownian motion is infinite, its quadratic variation is finite and deterministic.

Over a partition Π={0=t0<t1<⋯<tn=t}\Pi = \{0 = t_0 < t_1 < \dots < t_n = t\} with mesh size ∥Π∥→0\|\Pi\| \to 0: [W,W]t=lim⁡∥Π∥→0∑i=1n(Wti−Wti−1)2=tin L2[W, W]_t = \lim_{\|\Pi\| \to 0} \sum_{i=1}^n (W_{t_i} - W_{t_{i-1}})^2 = t \quad \text{in } L^2

In differential notation, this yields the fundamental multiplication rule of stochastic calculus: (dWt)2=dt(dW_t)^2 = dt

Intuition: Brownian increments scale with the square root of time: ΔWt∼N(0,Δt)  ⟹  ΔWt=O(Δt)\Delta W_t \sim \mathcal{N}(0, \Delta t) \implies \Delta W_t = \mathcal{O}(\sqrt{\Delta t}) In classical calculus, terms of order (Δt)2(\Delta t)^2 vanish when taking limits. But for Brownian motion: (ΔWt)2∼O(Δt)(\Delta W_t)^2 \sim \mathcal{O}(\Delta t) Because (ΔWt)2(\Delta W_t)^2 is of the first order in time, it cannot be discarded in a Taylor expansion.

Multiplication Rules of Stochastic Calculus:
  dt   × dt   = 0
  dt   × dW_t = 0
  dW_t × dt   = 0
  dW_t × dW_t = dt

3. What is Itô's Lemma, and where does the extra term come from?

Itô's Lemma is the stochastic counterpart of the chain rule. If XtX_t follows the Itô drift-diffusion process: dXt=μ(t,Xt)dt+σ(t,Xt)dWtdX_t = \mu(t, X_t) dt + \sigma(t, X_t) dW_t

And f(t,Xt)f(t, X_t) is a twice continuously differentiable function C1,2C^{1,2}, then the differential df(t,Xt)df(t, X_t) is:

df(t,Xt)=(∂f∂t+μ∂f∂x+12σ2∂2f∂x2)dt+σ∂f∂xdWtdf(t, X_t) = \left( \frac{\partial f}{\partial t} + \mu \frac{\partial f}{\partial x} + \frac{1}{2} \sigma^2 \frac{\partial^2 f}{\partial x^2} \right) dt + \sigma \frac{\partial f}{\partial x} dW_t

Where does the extra 12σ2∂2f∂x2dt\frac{1}{2}\sigma^2 \frac{\partial^2 f}{\partial x^2} dt come from? It arises directly from the second-order term of the Taylor expansion: Δf=∂f∂tΔt+∂f∂xΔXt+12∂2f∂x2(ΔXt)2+…\Delta f = \frac{\partial f}{\partial t} \Delta t + \frac{\partial f}{\partial x} \Delta X_t + \frac{1}{2} \frac{\partial^2 f}{\partial x^2} (\Delta X_t)^2 + \dots Substituting ΔXt=μΔt+σΔWt\Delta X_t = \mu \Delta t + \sigma \Delta W_t: (ΔXt)2=σ2(ΔWt)2+O((Δt)3/2)=σ2Δt(\Delta X_t)^2 = \sigma^2 (\Delta W_t)^2 + \mathcal{O}((\Delta t)^{3/2}) = \sigma^2 \Delta t This is Itô's correction term (often called the convexity or Jensen's inequality correction). If ff is convex (∂2f∂x2>0\frac{\partial^2 f}{\partial x^2} > 0), randomness inherently drifts the expected value upward.


4. What is the difference between Itô and Stratonovich integration?

The distinction lies in the evaluation point chosen when defining the Riemann-type sum: ∫0tYs∘dWs=lim⁡∥Π∥→0∑i=0n−1Yτi(Wti+1−Wti)\int_0^t Y_s \circ dW_s = \lim_{\|\Pi\| \to 0} \sum_{i=0}^{n-1} Y_{\tau_i} (W_{t_{i+1}} - W_{t_i})

  1. →

    Itô Integral (τi=ti\tau_i = t_i, left endpoint):

    • →Definition: Evaluated at the beginning of each interval: ∑YtiΔWi\sum Y_{t_i} \Delta W_i.
    • →Key Property: It is non-anticipative (causal). The integrator does not look into the future.
    • →Martingale Property: If YtY_t is square-integrable, E[∫0tYsdWs]=0\mathbb{E}\left[ \int_0^t Y_s dW_s \right] = 0.
    • →Consequence: The classical chain rule breaks; requires Itô's Lemma.
    • →Domain: Used universally in quantitative finance and economics, where agents cannot see future prices.
  2. →

    Stratonovich Integral (τi=ti+ti+12\tau_i = \frac{t_i + t_{i+1}}{2}, midpoint):

    • →Definition: Evaluated at the midpoint: ∑Yti+ti+12ΔWi\sum Y_{\frac{t_i + t_{i+1}}{2}} \Delta W_i, denoted with ∘\circ.
    • →Key Property: Obeys standard Newtonian chain rule without second-order correction terms!
    • →Consequence: Not a martingale; E[∫0tYs∘dWs]≠0\mathbb{E}\left[ \int_0^t Y_s \circ dW_s \right] \ne 0.
    • →Domain: Used in physics and engineering (e.g., Langevin equations with colored physical noise tending to white noise).

5. What is Girsanov's Theorem, and why is it essential in pricing?

Girsanov's Theorem describes how the drift of a stochastic process changes when switching from one equivalent probability measure to another.

In quantitative finance:

  • →Under the physical (real-world) measure P\mathbb{P}, an asset drifts at its expected return μ\mu: dSt=μStdt+σStdWtPdS_t = \mu S_t dt + \sigma S_t dW_t^{\mathbb{P}}
  • →Under the risk-neutral measure Q\mathbb{Q}, discounted asset prices must be martingales. Therefore, the asset must drift at the risk-free rate rr: dSt=rStdt+σStdWtQdS_t = r S_t dt + \sigma S_t dW_t^{\mathbb{Q}}

Girsanov's theorem defines the new Brownian motion WtQW_t^{\mathbb{Q}} via the Radon-Nikodym derivative process dQdP\frac{d\mathbb{Q}}{d\mathbb{P}}: WtQ=WtP+∫0tμ−rσdsW_t^{\mathbb{Q}} = W_t^{\mathbb{P}} + \int_0^t \frac{\mu - r}{\sigma} ds

Crucial Insight: Changing measures alters the drift (μ→r\mu \to r), but leaves volatility σ\sigma and quadratic variation completely invariant.


6. What is the Feynman-Kac Formula?

The Feynman-Kac Theorem bridges the gap between stochastic differential equations and deterministic partial differential equations (PDEs).

If a function V(t,x)V(t, x) solves the parabolic terminal-value PDE: ∂V∂t+μ(t,x)∂V∂x+12σ2(t,x)∂2V∂x2−rV=0,V(T,x)=g(x)\frac{\partial V}{\partial t} + \mu(t, x) \frac{\partial V}{\partial x} + \frac{1}{2} \sigma^2(t, x) \frac{\partial^2 V}{\partial x^2} - r V = 0, \quad V(T, x) = g(x)

Then V(t,x)V(t, x) can be represented as the conditional expectation under the risk-neutral measure: V(t,x)=EQ[e−r(T−t)g(XT)∣Xt=x]V(t, x) = \mathbb{E}^{\mathbb{Q}} \left[ e^{-r(T - t)} g(X_T) \mid X_t = x \right] where dXt=μ(t,Xt)dt+σ(t,Xt)dWtdX_t = \mu(t, X_t) dt + \sigma(t, X_t) dW_t.

This theorem proves that pricing a derivative via a Black-Scholes PDE or via Monte Carlo expectation are mathematically identical.


Part II: Practical Interview Problems & Quant Puzzles

Problem 1: Solve Geometric Brownian Motion (GBM)

Problem: Find the analytical solution to the SDE: dSt=μStdt+σStdWt,S0>0dS_t = \mu S_t dt + \sigma S_t dW_t, \quad S_0 > 0

Solution:

  1. →Define the transformation f(St)=ln⁡Stf(S_t) = \ln S_t.
  2. →Compute partial derivatives: ∂f∂S=1S,∂2f∂S2=−1S2,∂f∂t=0\frac{\partial f}{\partial S} = \frac{1}{S}, \quad \frac{\partial^2 f}{\partial S^2} = -\frac{1}{S^2}, \quad \frac{\partial f}{\partial t} = 0
  3. →Apply Itô's Lemma: d(ln⁡St)=1StdSt+12(−1St2)(dSt)2d(\ln S_t) = \frac{1}{S_t} dS_t + \frac{1}{2} \left( -\frac{1}{S_t^2} \right) (dS_t)^2
  4. →Substitute dSt=μStdt+σStdWtdS_t = \mu S_t dt + \sigma S_t dW_t and (dSt)2=σ2St2dt(dS_t)^2 = \sigma^2 S_t^2 dt: d(ln⁡St)=1St(μStdt+σStdWt)−12St2(σ2St2dt)d(\ln S_t) = \frac{1}{S_t} (\mu S_t dt + \sigma S_t dW_t) - \frac{1}{2 S_t^2} (\sigma^2 S_t^2 dt) d(ln⁡St)=(μ−12σ2)dt+σdWtd(\ln S_t) = \left( \mu - \frac{1}{2}\sigma^2 \right) dt + \sigma dW_t
  5. →Integrate both sides from 00 to tt: ln⁡St−ln⁡S0=(μ−12σ2)t+σWt\ln S_t - \ln S_0 = \left( \mu - \frac{1}{2}\sigma^2 \right) t + \sigma W_t
  6. →Exponentiate: St=S0exp⁡((μ−12σ2)t+σWt)\mathbf{S_t = S_0 \exp\left( \left( \mu - \frac{1}{2}\sigma^2 \right) t + \sigma W_t \right)}

Problem 2: Evaluate the Stochastic Integral ∫0tWsdWs\int_0^t W_s dW_s

Problem: Compute ∫0tWsdWs\int_0^t W_s dW_s and find its expectation.

Solution: In classical calculus, ∫xdx=12x2\int x dx = \frac{1}{2}x^2. Let's test f(Wt)=Wt2f(W_t) = W_t^2 using Itô's Lemma:

  1. →f(x)=x2  ⟹  f′(x)=2x,f′′(x)=2f(x) = x^2 \implies f'(x) = 2x, \quad f''(x) = 2.
  2. →Applying Itô's Lemma to Wt2W_t^2: d(Wt2)=2WtdWt+12(2)dt=2WtdWt+dtd(W_t^2) = 2 W_t dW_t + \frac{1}{2}(2) dt = 2 W_t dW_t + dt
  3. →Rearranging: WtdWt=12d(Wt2)−12dtW_t dW_t = \frac{1}{2} d(W_t^2) - \frac{1}{2} dt
  4. →Integrating from 00 to tt (noting W0=0W_0 = 0): ∫0tWsdWs=12Wt2−12t\mathbf{\int_0^t W_s dW_s = \frac{1}{2} W_t^2 - \frac{1}{2} t}

Expectation Verification: E[∫0tWsdWs]=E[12Wt2−12t]=12E[Wt2]−12t\mathbb{E}\left[ \int_0^t W_s dW_s \right] = \mathbb{E}\left[ \frac{1}{2} W_t^2 - \frac{1}{2} t \right] = \frac{1}{2} \mathbb{E}[W_t^2] - \frac{1}{2}t Since E[Wt2]=Var(Wt)=t\mathbb{E}[W_t^2] = \text{Var}(W_t) = t: 12(t)−12(t)=0\frac{1}{2}(t) - \frac{1}{2}(t) = 0 This perfectly confirms the Itô isometry property that the expectation of an Itô integral with respect to Brownian motion is zero.


Problem 3: Find the Fourth Moment E[Wt4]\mathbb{E}[W_t^4]

Problem: Calculate E[Wt4]\mathbb{E}[W_t^4] using Itô's Lemma.

Solution:

  1. →Let f(x)=x4  ⟹  f′(x)=4x3,f′′(x)=12x2f(x) = x^4 \implies f'(x) = 4x^3, \quad f''(x) = 12x^2.
  2. →By Itô's Lemma: d(Wt4)=4Wt3dWt+6Wt2dtd(W_t^4) = 4 W_t^3 dW_t + 6 W_t^2 dt
  3. →Write in integral form: Wt4=4∫0tWs3dWs+6∫0tWs2dsW_t^4 = 4 \int_0^t W_s^3 dW_s + 6 \int_0^t W_s^2 ds
  4. →Take expectations of both sides. Since E[∫0tWs3dWs]=0\mathbb{E}\left[\int_0^t W_s^3 dW_s\right] = 0: E[Wt4]=6∫0tE[Ws2]ds\mathbb{E}[W_t^4] = 6 \int_0^t \mathbb{E}[W_s^2] ds
  5. →Knowing E[Ws2]=s\mathbb{E}[W_s^2] = s: E[Wt4]=6∫0ts ds=6[t22]=3t2\mathbb{E}[W_t^4] = 6 \int_0^t s \, ds = 6 \left[ \frac{t^2}{2} \right] = \mathbf{3t^2}

This matches the standard normal kurtosis formula: for Z∼N(0,1)Z \sim \mathcal{N}(0, 1), E[Z4]=3  ⟹  E[(tZ)4]=3t2\mathbb{E}[Z^4] = 3 \implies \mathbb{E}[(\sqrt{t}Z)^4] = 3t^2.


Problem 4: The Ornstein-Uhlenbeck (Mean-Reverting) Process

Problem: Solve the Vasicek / Ornstein-Uhlenbeck SDE: dXt=θ(μ−Xt)dt+σdWt,θ>0dX_t = \theta(\mu - X_t) dt + \sigma dW_t, \quad \theta > 0 Find the conditional expectation E[Xt∣X0]\mathbb{E}[X_t \mid X_0] and conditional variance Var(Xt∣X0)\text{Var}(X_t \mid X_0).

Solution:

  1. →Multiply by the integrating factor eθte^{\theta t}: d(eθtXt)=θeθtXtdt+eθtdXtd(e^{\theta t} X_t) = \theta e^{\theta t} X_t dt + e^{\theta t} dX_t
  2. →Substitute dXtdX_t: d(eθtXt)=θeθtXtdt+eθt(θ(μ−Xt)dt+σdWt)d(e^{\theta t} X_t) = \theta e^{\theta t} X_t dt + e^{\theta t} (\theta(\mu - X_t) dt + \sigma dW_t) d(eθtXt)=θμeθtdt+σeθtdWtd(e^{\theta t} X_t) = \theta \mu e^{\theta t} dt + \sigma e^{\theta t} dW_t
  3. →Integrate from 00 to tt: eθtXt−X0=θμ∫0teθsds+σ∫0teθsdWse^{\theta t} X_t - X_0 = \theta \mu \int_0^t e^{\theta s} ds + \sigma \int_0^t e^{\theta s} dW_s eθtXt−X0=θμ(eθt−1θ)+σ∫0teθsdWse^{\theta t} X_t - X_0 = \theta \mu \left( \frac{e^{\theta t} - 1}{\theta} \right) + \sigma \int_0^t e^{\theta s} dW_s
  4. →Multiply through by e−θte^{-\theta t}: Xt=X0e−θt+μ(1−e−θt)+σ∫0te−θ(t−s)dWs\mathbf{X_t = X_0 e^{-\theta t} + \mu (1 - e^{-\theta t}) + \sigma \int_0^t e^{-\theta(t - s)} dW_s}

Moments:

  • →Mean: E[Xt]=X0e−θt+μ(1−e−θt)→t→∞μ\mathbb{E}[X_t] = X_0 e^{-\theta t} + \mu (1 - e^{-\theta t}) \xrightarrow{t \to \infty} \mu.
  • →Variance (via Itô Isometry): Var(Xt)=E[(σ∫0te−θ(t−s)dWs)2]=σ2∫0te−2θ(t−s)ds=σ22θ(1−e−2θt)\text{Var}(X_t) = \mathbb{E}\left[ \left( \sigma \int_0^t e^{-\theta(t - s)} dW_s \right)^2 \right] = \sigma^2 \int_0^t e^{-2\theta(t - s)} ds = \mathbf{\frac{\sigma^2}{2\theta} (1 - e^{-2\theta t})}
  • →Steady-State Variance: As t→∞t \to \infty, Var(X∞)=σ22θ\text{Var}(X_\infty) = \frac{\sigma^2}{2\theta}.

Summary Cheat Sheet for Interviews

ConceptKey Equation / Property
Brownian IncrementE[dWt]=0,Var(dWt)=dt\mathbb{E}[dW_t] = 0, \quad \text{Var}(dW_t) = dt
Quadratic Variation(dWt)2=dt,dt⋅dWt=0(dW_t)^2 = dt, \quad dt \cdot dW_t = 0
Itô 1D Ruledf=ftdt+fxdXt+12fxx(dXt)2df = f_t dt + f_x dX_t + \frac{1}{2} f_{xx} (dX_t)^2
Itô IsometryE[(∫0tHsdWs)2]=∫0tE[Hs2]ds\mathbb{E}\left[ \left(\int_0^t H_s dW_s\right)^2 \right] = \int_0^t \mathbb{E}[H_s^2] ds
Martingale ConditionE[Mt∣Fs]=Ms(∀s≤t)\mathbb{E}[M_t \mid \mathcal{F}_s] = M_s \quad (\forall s \le t)
Girsanov PropertyPreserves volatility σ\sigma; shifts drift by λ=μ−rσ\lambda = \frac{\mu - r}{\sigma}

Institutional Proof

Dive deeper into Stochastic Calculus

See the complete formal proof, animated visual derivations, and the full architectural breakdown in the library.

Enter the Library →

The Journal

Subscribe for bi-weekly deep dives into abstract mathematics and statistical inference.