Brownian Motion and Quadratic Variation

Medium·18 min read·Interactive lab
Stochastic CalculusBrownian MotionQuadratic Variation

Setup

Let (Ω,F,P)(\Omega, \mathcal{F}, \mathbb{P}) be a complete probability space equipped with a filtration F=(Ft)t0\mathbb{F} = (\mathcal{F}_t)_{t \geq 0} satisfying the usual conditions: right-continuity (Ft=Ft+\mathcal{F}_t = \mathcal{F}_{t+}) and completeness (null sets of P\mathbb{P} are in F0\mathcal{F}_0).

All processes are defined on this space. Time horizon is [0,T][0, T] for a fixed T>0T > 0. Expectations are under P\mathbb{P} unless stated otherwise.


Definition

A standard Brownian motion (or Wiener process) is a stochastic process (Wt)t0(W_t)_{t \geq 0} satisfying:

  1. Initial condition: W0=0W_0 = 0 almost surely.
  2. Independent increments: for 0s<t0 \leq s < t, the increment WtWsW_t - W_s is independent of Fs\mathcal{F}_s.
  3. Gaussian increments: WtWsN(0,ts)W_t - W_s \sim \mathcal{N}(0,\, t - s) for all 0st0 \leq s \leq t.
  4. Continuous paths: tWt(ω)t \mapsto W_t(\omega) is continuous P\mathbb{P}-almost surely.

These axioms are not redundant. Axiom 4 is a topological regularity condition, independent of the distributional structure in axioms 2–3. Without it, the process would be defined only up to modification, and pathwise integration would be meaningless.

The covariance structure follows immediately: for sts \leq t, Cov(Ws,Wt)=E[WsWt]=s.\mathrm{Cov}(W_s, W_t) = \mathbb{E}[W_s W_t] = s.


Sample Path Properties

Hölder Continuity

By the Kolmogorov–Chentsov continuity theorem, a process with E[WtWsp]Cts1+α\mathbb{E}[|W_t - W_s|^p] \leq C |t - s|^{1 + \alpha} for some p,α,C>0p, \alpha, C > 0 admits a version with Hölder-continuous paths of any exponent γ<α/p\gamma < \alpha/p.

For Brownian motion, E[WtWs2k]=(2k1)!!tsk\mathbb{E}[|W_t - W_s|^{2k}] = (2k-1)!! \cdot |t-s|^k. Taking k=2k = 2: E[WtWs4]=3(ts)2\mathbb{E}[|W_t - W_s|^4] = 3(t-s)^2, so we may take p=4p = 4, α=1\alpha = 1, giving Hölder exponent up to γ<1/4\gamma < 1/4. Sharper analysis yields Hölder exponent γ<1/2\gamma < 1/2.

Brownian motion paths are not Hölder-1/21/2: the modulus of continuity is limh0suptshWtWs2hlog(1/h)=1a.s.\lim_{h \to 0} \sup_{|t - s| \leq h} \frac{|W_t - W_s|}{\sqrt{2h \log(1/h)}} = 1 \quad \text{a.s.} (the Lévy modulus). The hlog(1/h)\sqrt{h \log(1/h)} factor, not just h\sqrt{h}, captures the exact regularity.

Nowhere Differentiability

With probability one, Brownian motion is nowhere differentiable. To see why: if the path were differentiable at some point tt, then (Wt+hWt)/hWt(W_{t+h} - W_t)/h \to W'_t as h0h \to 0. But this quotient has standard deviation h1/2h^{-1/2}, which diverges. A formal proof uses Paley–Wiener–Zygmund or a direct Borel–Cantelli argument.

This is the mathematical source of the informal statement "Brownian motion fluctuates on every scale."

Infinite Total Variation

Define the total variation of WW on [0,T][0, T] over a partition π={0=t0<t1<<tn=T}\pi = \{0 = t_0 < t_1 < \cdots < t_n = T\}: TV(W,π)=i=1nWtiWti1.TV(W, \pi) = \sum_{i=1}^{n} |W_{t_i} - W_{t_{i-1}}|.

The total variation of Brownian motion is infinite almost surely: supπTV(W,π)=+a.s.\sup_{\pi} TV(W, \pi) = +\infty \quad \text{a.s.}

This has a fundamental consequence: Lebesgue–Stieltjes integration of the form fdW\int f \, dW cannot be defined pathwise. The standard integration-by-parts formula requires finite variation. For Brownian motion, a new theory — the Itô integral — is required.


Quadratic Variation

The quadratic variation of a process XX on [0,t][0, t] is defined as the limit in probability: [X]t=limπ0i=1n(XtiXti1)2,[X]_t = \lim_{|\pi| \to 0} \sum_{i=1}^{n} (X_{t_i} - X_{t_{i-1}})^2, where the limit is taken as the mesh π=maxi(titi1)0|\pi| = \max_i(t_i - t_{i-1}) \to 0 over arbitrary partitions.

Theorem: [W]t=t[W]_t = t

For standard Brownian motion, [W]t=t[W]_t = t almost surely for all t0t \geq 0.

Proof. Fix an equal-mesh partition with nn intervals and mesh h=t/nh = t/n. Let Δi=WtiWti1\Delta_i = W_{t_i} - W_{t_{i-1}} and define: Qn=i=1nΔi2.Q_n = \sum_{i=1}^{n} \Delta_i^2.

Since ΔiN(0,h)\Delta_i \sim \mathcal{N}(0, h) and all Δi\Delta_i are independent: E[Δi2]=h,Var(Δi2)=E[Δi4]h2=3h2h2=2h2.\mathbb{E}[\Delta_i^2] = h, \qquad \mathrm{Var}(\Delta_i^2) = \mathbb{E}[\Delta_i^4] - h^2 = 3h^2 - h^2 = 2h^2.

Therefore: E[Qn]=nh=t,Var(Qn)=n2h2=2t2/nn0.\mathbb{E}[Q_n] = nh = t, \qquad \mathrm{Var}(Q_n) = n \cdot 2h^2 = 2t^2/n \xrightarrow[n \to \infty]{} 0.

By Chebyshev's inequality, QntQ_n \to t in L2L^2, and hence in probability. Convergence a.s. follows by a subsequence argument extended to general (non-equal) partitions via an L2L^2 monotonicity argument. \square

Differential Notation

The result is commonly written as: (dWt)2=dt.(dW_t)^2 = dt.

This topic requires Premium

Only today's featured topic is free. Unlock the full Today's Focus archive with Premium.

Read the theory? Now explore it interactively.

Open the Lab