Stochastic Processes — Brownian Motion and Martingales

⚡ Key Takeaways

Brownian motion is the canonical continuous-time stochastic process with independent, stationary, Gaussian increments and continuous paths; a martingale is a process whose conditional expected future value equals its present value. Key result: Brownian motion is the prototypical martingale, and the two notions are linked by Lévy's characterization (a continuous martingale with quadratic variation tt is Brownian motion) and by the Itô stochastic calculus. Applies to: mathematical finance, statistical physics, PDE theory, filtering, population genetics. Originated with Brown (1827), Einstein (1905), Bachelier (1900), formalized by Wiener (1923), Lévy, Doob, and Itô (1940s–50s).

💡 The Core Insight

Randomness accumulated continuously over time, when made to have independent increments of variance proportional to elapsed time, forces a Gaussian limit (via the CLT) and a path so rough it has no derivative but a perfectly predictable "energy" ⟨W⟩t=t\langle W \rangle_t = t. The martingale property—conditional fairness—is what allows this roughness to be tamed: stochastic integrals against martingales are again martingales, and quadratic variation replaces the chain rule, giving the Itô formula. Brownian motion is simultaneously a Markov process, a Gaussian process, and a continuous martingale; these three identities, when forced to coexist, pin it down uniquely.

Formal Statement

DefinitionFiltered Probability Space

A filtered probability space (Ω,F,{Ft}t≥0,P)(\Omega, \mathcal{F}, \{\mathcal{F}_t\}_{t\ge 0}, \mathbb{P}) consists of a probability space together with a filtration: an increasing family of σ\sigma-algebras Fs⊆Ft⊆F\mathcal{F}_s \subseteq \mathcal{F}_t \subseteq \mathcal{F} for s≤ts \le t. It satisfies the usual conditions if it is right-continuous (Ft=⋂u>tFu\mathcal{F}_t = \bigcap_{u>t}\mathcal{F}_u) and F0\mathcal{F}_0 contains all P\mathbb{P}-null sets. A process XX is adapted if XtX_t is Ft\mathcal{F}_t-measurable for all tt.

DefinitionBrownian Motion

A real-valued process {Wt}t≥0\{W_t\}_{t\ge 0} is a standard Brownian motion (Wiener process) if: 1. W0=0W_0 = 0 almost surely; 2. independent increments: for 0≤t0<t1<⋯<tn0\le t_0 < t_1 < \cdots < t_n, the increments Wti−Wti−1W_{t_{i}}-W_{t_{i-1}} are independent; 3. stationary Gaussian increments: Wt−Ws∼N(0,t−s)W_t - W_s \sim \mathcal{N}(0, t-s) for s<ts<t; 4. continuous paths: t↦Wtt \mapsto W_t is continuous almost surely.

1
P(Wt−Ws∈dx)=12π(t−s)exp⁡ ⁣(−x22(t−s))dx\quad \mathbb{P}(W_t - W_s \in dx) = \frac{1}{\sqrt{2\pi(t-s)}}\exp\!\left(-\frac{x^2}{2(t-s)}\right)dx

Derivation hint: this is the transition density; it solves the heat equation ∂tp=12∂xxp\partial_t p = \tfrac12 \partial_{xx} p with p(0,⋅)=δ0p(0,\cdot)=\delta_0, encoding the Gaussian increment law.

DefinitionMartingale, Submartingale, Supermartingale

An adapted, integrable (E∣Mt∣<∞\mathbb{E}|M_t|<\infty) process is a martingale if

2
E[Mt∣Fs]=Msfor all s≤t.\quad \mathbb{E}[M_t \mid \mathcal{F}_s] = M_s \quad \text{for all } s \le t.

Derivation hint: replace "==" with "≥\ge" for a submartingale and "≤\le" for a supermartingale; the equality is the tower property applied to a fair-game increment with zero conditional mean. The defining relation says the best forecast of MtM_t given information up to ss is the current value MsM_s.

TheoremBrownian Motion is a Martingale; Associated Martingales

Let WW be a standard Brownian motion with its natural filtration. Then each of the following is a martingale:

3
Wt,Wt2−t,exp⁡ ⁣(λWt−12λ2t) (λ∈R).\quad W_t,\qquad W_t^2 - t,\qquad \exp\!\left(\lambda W_t - \tfrac12 \lambda^2 t\right)\ (\lambda\in\mathbb{R}).

Derivation hint: for WtW_t, use independent increments and zero mean; for Wt2−tW_t^2-t, expand Wt2=(Wt−Ws)2+2Ws(Wt−Ws)+Ws2W_t^2=(W_t-W_s)^2+2W_s(W_t-W_s)+W_s^2 and take conditional expectation; for the exponential, use the Gaussian MGF E[eλ(Wt−Ws)∣Fs]=eλ2(t−s)/2\mathbb{E}[e^{\lambda(W_t-W_s)}\mid\mathcal F_s]=e^{\lambda^2(t-s)/2}.

TheoremLévy's Characterization

Let MM be a continuous local martingale with M0=0M_0=0 and quadratic variation ⟨M⟩t=t\langle M\rangle_t = t. Then MM is a standard Brownian motion.

4
⟨M⟩t=t and M cont. local mart. ⟹ MtWt\quad \langle M\rangle_t = t \ \text{and } M \text{ cont. local mart.} \ \Longrightarrow\ M_t \stackrel{d}{=} W_t

Derivation hint: apply Itô's formula to exp⁡(iλMt+12λ2⟨M⟩t)\exp(i\lambda M_t + \tfrac12\lambda^2\langle M\rangle_t); the stochastic integral term is a martingale, forcing the conditional characteristic function to be that of N(0,t−s)\mathcal N(0,t-s), the law of an increment.

Why It Holds / Derivation

ProofQuadratic Variation of Brownian Motion

We show that over [0,t][0,t] with partition Πn={0=t0<⋯<tn=t}\Pi_n = \{0=t_0<\cdots<t_n=t\} of mesh ∥Πn∥→0\|\Pi_n\|\to 0, the sums of squared increments converge to tt.

Step 1. Define Qn=∑i=1n(Wti−Wti−1)2Q_n = \sum_{i=1}^n (W_{t_i}-W_{t_{i-1}})^2. By stationary increments, E[(Wti−Wti−1)2]=ti−ti−1\mathbb{E}[(W_{t_i}-W_{t_{i-1}})^2] = t_i - t_{i-1}, so

5
E[Qn]=∑i=1n(ti−ti−1)=t.\quad \mathbb{E}[Q_n] = \sum_{i=1}^n (t_i - t_{i-1}) = t.

Derivation hint: variance of a centered Gaussian increment equals its time span; sum telescopes.

Step 2. Increments are independent, so the variance of QnQ_n adds:

6
Var⁡(Qn)=∑i=1nVar⁡((Wti−Wti−1)2)=2∑i=1n(ti−ti−1)2.\quad \operatorname{Var}(Q_n) = \sum_{i=1}^n \operatorname{Var}\big((W_{t_i}-W_{t_{i-1}})^2\big) = 2\sum_{i=1}^n (t_i-t_{i-1})^2.

Derivation hint: for Z∼N(0,σ2)Z\sim\mathcal N(0,\sigma^2), Var⁡(Z2)=2σ4\operatorname{Var}(Z^2)=2\sigma^4; here σ2=ti−ti−1\sigma^2=t_i-t_{i-1}.

Step 3. Bound the sum by the mesh:

7
Var⁡(Qn)=2∑i=1n(ti−ti−1)2≤2 ∥Πn∥∑i=1n(ti−ti−1)=2 ∥Πn∥ t→0.\quad \operatorname{Var}(Q_n) = 2\sum_{i=1}^n (t_i-t_{i-1})^2 \le 2\,\|\Pi_n\|\sum_{i=1}^n (t_i-t_{i-1}) = 2\,\|\Pi_n\|\, t \to 0.

Derivation hint: factor out the largest interval; the remaining sum telescopes to tt.

Step 4. Since E[Qn]=t\mathbb{E}[Q_n]=t and Var⁡(Qn)→0\operatorname{Var}(Q_n)\to 0, Chebyshev gives Qn→tQ_n \to t in L2L^2 (and in probability). For nested partitions or summable meshes, the convergence is almost sure. Hence ⟨W⟩t=t\langle W\rangle_t = t. ■\blacksquare

∎
ProofWt2−tW_t^2 - t is a Martingale

Step 1. For s<ts<t write Wt=Ws+(Wt−Ws)W_t = W_s + (W_t - W_s), so

8
Wt2=Ws2+2Ws(Wt−Ws)+(Wt−Ws)2.\quad W_t^2 = W_s^2 + 2W_s(W_t - W_s) + (W_t - W_s)^2.

Derivation hint: binomial expansion of (Ws+Δ)2(W_s + \Delta)^2.

Step 2. Take E[⋅∣Fs]\mathbb{E}[\cdot\mid\mathcal F_s]. Since WsW_s is Fs\mathcal F_s-measurable and the increment is independent of Fs\mathcal F_s with mean 00 and variance t−st-s:

9
E[Wt2∣Fs]=Ws2+2Ws⋅0+(t−s)=Ws2+(t−s).\quad \mathbb{E}[W_t^2\mid\mathcal F_s] = W_s^2 + 2W_s\cdot 0 + (t-s) = W_s^2 + (t-s).

Derivation hint: pull out Fs\mathcal F_s-measurable factors; use E[Δ∣Fs]=0\mathbb{E}[\Delta\mid\mathcal F_s]=0 and E[Δ2∣Fs]=t−s\mathbb{E}[\Delta^2\mid\mathcal F_s]=t-s.

Step 3. Subtract tt from both sides:

10
E[Wt2−t∣Fs]=Ws2+(t−s)−t=Ws2−s.\quad \mathbb{E}[W_t^2 - t \mid \mathcal F_s] = W_s^2 + (t-s) - t = W_s^2 - s.

Derivation hint: rearrange; the result is the value at ss, confirming the martingale property. ■\blacksquare

∎
TheoremItô's Formula

Let f∈C2(R)f\in C^2(\mathbb{R}) and WW a Brownian motion. Then a.s.,

11
f(Wt)=f(W0)+∫0tf′(Ws) dWs+12∫0tf′′(Ws) ds.\quad f(W_t) = f(W_0) + \int_0^t f'(W_s)\,dW_s + \frac12\int_0^t f''(W_s)\,ds.

Derivation hint: Taylor-expand f(Wti)−f(Wti−1)f(W_{t_{i}})-f(W_{t_{i-1}}) to second order; the first-order terms build the Itô integral, while ∑f′′(⋅)(Wti−Wti−1)2→∫f′′ d⟨W⟩=∫f′′ ds\sum f''(\cdot)(W_{t_i}-W_{t_{i-1}})^2 \to \int f''\,d\langle W\rangle = \int f''\,ds because ( dW )2=dt(\,dW\,)^2 = dt.

Figure
Brownian sample path with its variance envelope ±t\pm\sqrt{t}
t W +√t −√t W_t (rough, nowhere differentiable)
mathematician · Norbert Wiener | 1894–1964

Wiener gave the first rigorous mathematical construction of Brownian motion (1923), proving that a probability measure—the Wiener measure—exists on the space C[0,∞)C[0,\infty) of continuous functions realizing Einstein's physical predictions. He established the continuity and nowhere-differentiability of almost every path, and developed the abstract machinery (the Wiener integral) that anticipated Itô's calculus. Beyond probability, Wiener founded cybernetics, contributed to harmonic analysis (the Paley–Wiener theorems), and shaped wartime work on prediction and filtering. The Wiener process bears his name precisely because he transformed a physical phenomenon into a well-defined object of pure mathematics.

Building Intuition

🧠 Intuition

Picture a drunkard taking nn steps of size ±1/n\pm 1/\sqrt n per unit time. After time tt, position is a sum of ∼nt\sim nt tiny independent steps; the CLT makes this Gaussian with variance tt. Send n→∞n\to\infty: the discrete walk becomes Brownian motion. The variance growing linearly in time, not quadratically, is the signature of diffusion—random cancellation means displacement scales like t\sqrt t, not tt. The martingale property is "no information advantage": knowing the entire past gives no edge in predicting the next increment's direction, because E[Wt−Ws∣Fs]=0\mathbb{E}[W_t - W_s \mid \mathcal F_s] = 0. The roughness (( dW )2=dt(\,dW\,)^2 = dt) is what distinguishes stochastic from ordinary calculus: a Brownian path wiggles so much that second-order effects survive in the limit, producing the extra 12f′′\tfrac12 f'' term in Itô's formula. Think of it as a deterministic "drift" generated purely by the path's roughness.

AE
Albert Einstein
1879–1955

In his 1905 annus mirabilis paper on Brownian motion, Einstein derived the diffusion equation for suspended particles and predicted that the mean-squared displacement grows linearly in time: ⟨x2⟩=2Dt\langle x^2\rangle = 2Dt, with D=kBT/(6πηr)D = k_B T/(6\pi\eta r) relating the diffusion constant to temperature, viscosity, and particle radius. This gave a route to measuring Avogadro's number, providing decisive evidence for the atomic hypothesis (confirmed experimentally by Perrin, Nobel 1926). Einstein's physical reasoning—random molecular bombardment producing a Gaussian spreading—is exactly mirrored by the mathematical transition density [1][1].

Worked Examples

ExampleHitting Time and Optional Stopping

Given: Standard Brownian motion WW starting at 00. Let a>0a>0, b>0b>0, and τ=inf⁡{t:Wt=a or Wt=−b}\tau = \inf\{t: W_t = a \text{ or } W_t = -b\} be the first exit time from (−b,a)(-b, a). Find P(Wτ=a)\mathbb{P}(W_\tau = a).

Steps: 1. WtW_t is a martingale; τ\tau is a stopping time with E[τ]<∞\mathbb{E}[\tau]<\infty (since the gap is bounded). By the Optional Stopping Theorem, E[Wτ]=W0=0\mathbb{E}[W_\tau] = W_0 = 0. 2. WτW_\tau takes value aa with probability pp and −b-b with probability 1−p1-p. So

12
0=E[Wτ]=a p+(−b)(1−p)=p(a+b)−b.\quad 0 = \mathbb{E}[W_\tau] = a\,p + (-b)(1-p) = p(a+b) - b.

Derivation hint: equate the martingale's expected stopped value to its start, then solve the two-point distribution. 3. Solve: p=ba+bp = \dfrac{b}{a+b}.

Boxed answer: P(Wτ=a)=ba+b\displaystyle \boxed{\mathbb{P}(W_\tau = a) = \frac{b}{a+b}}. With a=2, b=3a=2,\ b=3: probability =3/5=0.6= 3/5 = 0.6 (dimensionless).

ExampleExpected Exit Time via Wt2−tW_t^2 - t

Given: Same setup; compute E[τ]\mathbb{E}[\tau], the expected time to exit (−b,a)(-b,a).

Steps: 1. The process Wt2−tW_t^2 - t is a martingale (Eq. [10][10]). Optional stopping gives E[Wτ2−τ]=0\mathbb{E}[W_\tau^2 - \tau] = 0, hence E[τ]=E[Wτ2]\mathbb{E}[\tau] = \mathbb{E}[W_\tau^2]. 2. Using p=b/(a+b)p = b/(a+b) from the previous example:

13
E[Wτ2]=a2⋅ba+b+b2⋅aa+b=ab(a+b)a+b=ab.\quad \mathbb{E}[W_\tau^2] = a^2\cdot\frac{b}{a+b} + b^2\cdot\frac{a}{a+b} = \frac{ab(a+b)}{a+b} = ab.

Derivation hint: weight the squared exit values by their probabilities; factor the common abab. 3. Therefore E[τ]=ab\mathbb{E}[\tau] = ab.

Boxed answer: E[τ]=ab\displaystyle \boxed{\mathbb{E}[\tau] = ab}. With a=2, b=3a=2,\ b=3: E[τ]=6\mathbb{E}[\tau] = 6 (time units).

Common Misconceptions

⚠️ Common mistake

Brownian motion is not differentiable. It is a frequent error to write "dWt/dtdW_t/dt" as if it were an ordinary velocity. Almost every path is nowhere differentiable, has infinite total variation on every interval, and "white noise" W˙\dot W exists only as a distribution (generalized process). Consequently the Itô integral ∫H dW\int H\,dW is not a Riemann–Stieltjes integral; it must be defined as an L2L^2 limit of left-endpoint Riemann sums, and the choice of evaluation point matters—the Stratonovich integral (midpoint) obeys the ordinary chain rule but is not a martingale.

⚠️ Common mistake

Optional stopping requires hypotheses. The identity E[Mτ]=E[M0]\mathbb{E}[M_\tau]=\mathbb{E}[M_0] can fail without integrability/uniform-integrability conditions or a bounded stopping time. A classic trap: τa=inf⁡{t:Wt=a}\tau_a = \inf\{t: W_t = a\} for a>0a>0 has P(τa<∞)=1\mathbb{P}(\tau_a<\infty)=1 yet E[τa]=∞\mathbb{E}[\tau_a]=\infty, and E[Wτa]=a≠0=W0\mathbb{E}[W_{\tau_a}] = a \neq 0 = W_0—optional stopping fails because {Wt∧τa}\{W_{t\wedge\tau_a}\} is not uniformly integrable.

Counterexample

A martingale need not converge. One might guess every martingale settles down. But WtW_t itself is a martingale that oscillates forever: by the law of the iterated logarithm, lim sup⁡t→∞Wt/2tlog⁡log⁡t=1\limsup_{t\to\infty} W_t/\sqrt{2t\log\log t} = 1 and lim inf⁡=−1\liminf = -1 a.s., so WtW_t does not converge—it visits +∞+\infty and −∞-\infty in the limsup sense. Martingale convergence requires an L1L^1-bound: sup⁡tE∣Mt∣<∞\sup_t \mathbb{E}|M_t|<\infty. Brownian motion has E∣Wt∣=2t/π→∞\mathbb{E}|W_t| = \sqrt{2t/\pi}\to\infty, so the hypothesis fails.

Connections

connection

Brownian motion and martingales sit at the crossroads of several theories. Via the Feynman–Kac formula, expectations of functionals of WW solve parabolic PDEs, linking Partial Differential Equations to probability: u(t,x)=E[ϕ(Wt)∣W0=x]u(t,x) = \mathbb{E}[\phi(W_t)\mid W_0=x] solves the heat equation. The Itô integral and Itô's formula generate the theory of Stochastic Differential Equations, whose solutions are diffusions. The Martingale Representation Theorem states every L2L^2 Brownian martingale is a stochastic integral ∫H dW\int H\,dW, underpinning Mathematical Finance and market completeness. Lévy's characterization [4][4] shows Brownian motion is the only continuous martingale with quadratic variation tt, and the Dambis–Dubins–Schwarz theorem reveals every continuous martingale is a time-changed Brownian motion. The increment structure makes WW a Markov process and a Gaussian process simultaneously, connecting to Markov Processes and Gaussian field theory.

Real-World Applications

🔬 Where this appears

Historical Development

⏳ History
- 1827 — Robert Brown observes erratic motion of pollen grains suspended in water under a microscope, describing the phenomenon now bearing his name.
- 1900 — Louis Bachelier, in Théorie de la spéculation, models stock prices with Brownian motion—five years before Einstein—introducing the heat-equation connection.
- 1905 — Albert Einstein derives the diffusion law ⟨x2⟩=2Dt\langle x^2\rangle = 2Dt, linking microscopic molecular motion to macroscopic diffusion and Avogadro's number.
- 1906 — Marian Smoluchowski independently develops a kinetic theory of Brownian motion.
- 1908 — Paul Langevin introduces the stochastic differential equation mx¨=−γx˙+ξ(t)m\ddot x = -\gamma\dot x + \xi(t), the first SDE in physics.
- 1923 — Norbert Wiener constructs Brownian motion rigorously as a measure on continuous-function space, proving path continuity and nowhere-differentiability.
- 1933–1939 — Paul Lévy develops fine path properties (law of iterated logarithm, local time, characterization theorem) and the modulus of continuity.
- 1940s — Joseph Doob systematizes martingale theory: the convergence theorem, optional stopping, and maximal inequalities, providing the abstract backbone.
- 1944–1951 — Kiyosi Itô defines the stochastic integral and proves Itô's formula, founding stochastic calculus and the theory of diffusion processes.

Subtleties & Limits

open-problem

While Brownian motion itself is fully understood, surrounding questions remain deep and active. The construction and uniqueness of solutions to SDEs with irregular (merely measurable or distributional) coefficients—regularization by noise—is an ongoing program (Davie, Krylov–Röckner, and the rough-paths/regularity-structures approaches of Lyons and Hairer). The multidimensional and infinite-dimensional analogues—stochastic PDEs driven by space-time white noise—require renormalization, as in Hairer's Fields-medal-winning theory of regularity structures resolving the KPZ and Φ34\Phi^4_3 equations. Fine geometric properties of the Brownian zero set, the exact Hausdorff measure of paths, and intersection local times in critical dimensions continue to be studied. In finance, rough volatility models (Gatheral–Jaisson–Rosenbaum) replace Brownian volatility with fractional Brownian motion of Hurst index H≈0.1H\approx 0.1, raising new questions about non-Markovian, non-semimartingale calculus where Itô's formula no longer directly applies.

See Also

Martingale Convergence Theorem · Optional Stopping Theorem · Markov Processes · Mathematical Finance · Fractional Brownian Motion · Local Time · Doob's Inequalities