Technopedia

The technologist’s encyclopedia

Brownian motion

Brownian motion (the Wiener process) is the canonical continuous-time, continuous-state stochastic process: a random path WtW_t that starts at zero, moves continuously, and has independent, stationary Gaussian increments Wt+s−Ws∼N(0,s)W_{t+s} - W_s \sim \mathcal{N}(0, s). It is the scaling limit of random walks (Donsker’s invariance principle) and the building block of stochastic calculus (Itô integration, stochastic differential equations).

Where it is relevant.

Builds on Random variables and expectation, Conditional probability, and Limit theorems.

Authoritative sources

Definition and axioms

A standard Brownian motion (Wiener process) W=(Wt)t≥0W = (W_t)_{t \ge 0} is a stochastic process satisfying:

  1. W0=0W_0 = 0 almost surely.
  2. Continuous paths: t↦Wtt \mapsto W_t is continuous almost surely.
  3. Independent increments: for 0≤t0<t1<⋯<tn0 \le t_0 < t_1 < \dots < t_n, the increments Wti−Wti−1W_{t_i} - W_{t_{i-1}} are mutually independent.
  4. Stationary Gaussian increments: Wt−Ws∼N(0,t−s)W_t - W_s \sim \mathcal{N}(0, t-s) for 0≤s<t0 \le s < t.

Consequences (cheatsheet). Wt∼N(0,t)W_t \sim \mathcal{N}(0, t), so E[Wt]=0\mathbb{E}[W_t] = 0 and Var(Wt)=t\mathrm{Var}(W_t) = t; the covariance is Cov(Ws,Wt)=min⁡(s,t)\mathrm{Cov}(W_s, W_t) = \min(s, t); and WW is simultaneously a martingale, a Gaussian process (kernel min⁡(s,t)\min(s,t)), and a Markov process.

Details

The four axioms each pin down one aspect of the process, and together they force the entire law of WW.

Step 1 — Start at zero and move continuously. Axioms 1 and 2 say the path begins at the origin and never jumps: t↦Wtt \mapsto W_t is a continuous curve. This is what separates Brownian motion from a jump process such as the Poisson process. Continuity alone, however, does not fix the distribution — a deterministic drift Wt=μtW_t = \mu t is also continuous.

Step 2 — Independent, stationary increments. Axioms 3 and 4 say the “steps” of the path are memoryless and time-homogeneous: the increment over any interval of length hh has the same law, N(0,h)\mathcal{N}(0, h), and increments over disjoint intervals are independent. This is the continuous-time analogue of the i.i.d. increments of a random walk. It makes WW a Markov process (the future depends only on the present) and a martingale (the best prediction of WtW_t given the past is WsW_s).

Step 3 — The Gaussian law is forced. Because the increments are independent and stationary, WtW_t is a sum of tt worth of independent Gaussian increments; the only law consistent with all four axioms is Wt∼N(0,t)W_t \sim \mathcal{N}(0, t). The covariance follows by conditioning on the overlap of the two intervals:

Cov(Ws,Wt)=Cov(Ws, Ws+(Wt−Ws))=Var(Ws)=s=min⁡(s,t),s≤t.\mathrm{Cov}(W_s, W_t) = \mathrm{Cov}\big(W_s,\, W_s + (W_t - W_s)\big) = \mathrm{Var}(W_s) = s = \min(s,t), \quad s \le t.

Example. Take s=2s = 2, t=5t = 5. Then W2∼N(0,2)W_2 \sim \mathcal{N}(0,2), W5∼N(0,5)W_5 \sim \mathcal{N}(0,5), and Cov(W2,W5)=2\mathrm{Cov}(W_2, W_5) = 2. The correlation is Corr(W2,W5)=2/2⋅5≈0.63\mathrm{Corr}(W_2, W_5) = 2/\sqrt{2 \cdot 5} \approx 0.63: the two values are positively correlated because the path from 00 to 22 is shared.

t W_t one sample path, t ∈ [0, 10]

Authoritative sources

Example problem

Compute E[W3W7]\mathbb{E}[W_3 W_7] and Var(W7−W3)\mathrm{Var}(W_7 - W_3).

Solution

E[W3W7]=Cov(W3,W7)=min⁡(3,7)=3\mathbb{E}[W_3 W_7] = \mathrm{Cov}(W_3, W_7) = \min(3,7) = 3 (both means are zero). And W7−W3∼N(0,7−3)=N(0,4)W_7 - W_3 \sim \mathcal{N}(0, 7-3) = \mathcal{N}(0,4), so Var(W7−W3)=4\mathrm{Var}(W_7 - W_3) = 4.

Construction as the scaling limit of random walks

Donsker’s invariance principle (functional central limit theorem). Let X1,X2,…X_1, X_2, \dots be i.i.d. with E[X1]=0\mathbb{E}[X_1] = 0 and Var(X1)=σ2<∞\mathrm{Var}(X_1) = \sigma^2 < \infty, and let Sn=X1+⋯+XnS_n = X_1 + \dots + X_n. The rescaled random-walk process

W(n)(t)=S⌊nt⌋σn,t≥0,W^{(n)}(t) = \frac{S_{\lfloor n t \rfloor}}{\sigma \sqrt{n}}, \quad t \ge 0,

converges in distribution (in C[0,∞)C[0,\infty)) to a standard Brownian motion WW.

Details

Step 1 — The CLT is a statement about one time. The central limit theorem says Sn/(σn)→dN(0,1)S_n/(\sigma\sqrt{n}) \xrightarrow{d} \mathcal{N}(0,1): at a single time t=1t = 1, the rescaled walk looks Gaussian. Donsker’s theorem upgrades this to a statement about the whole path: the entire rescaled trajectory converges, not just its value at one point.

Step 2 — Why the n\sqrt{n} scaling. The walk SnS_n has standard deviation σn\sigma\sqrt{n}, so dividing by σn\sigma\sqrt{n} keeps the vertical scale fixed. Horizontally, nn steps are compressed into the unit interval [0,1][0,1], so each step spans 1/n1/n in time. The vertical step size is σ/n\sigma/\sqrt{n} and the horizontal step size is 1/n1/n; their ratio is σ/n:1/n=σn:1\sigma/\sqrt{n} : 1/n = \sigma\sqrt{n} : 1, which grows like n\sqrt{n}. This is why the limit path is continuous yet wildly rough — the steps shrink in both directions, but the vertical steps shrink more slowly.

Step 3 — The limit inherits the axioms. The rescaled process has (approximately) independent, stationary increments and a Gaussian marginal at each time; Donsker’s theorem makes this rigorous and shows the limit is continuous, hence a Brownian motion.

Example. A simple symmetric random walk (Xi=±1X_i = \pm 1 with probability 1/21/2) has σ2=1\sigma^2 = 1. After n=104n = 10^4 steps the walk is typically at height ±100\pm 100; rescaled by 1/n=1/1001/\sqrt{n} = 1/100 and plotted over [0,1][0,1], it looks like a Brownian path.

t S_i / √i rescaled random walk → Brownian path (Donsker)

Authoritative sources

Example problem

Let XiX_i be i.i.d. with E[Xi]=0\mathbb{E}[X_i] = 0 and Var(Xi)=4\mathrm{Var}(X_i) = 4. What is the correct rescaling of S⌊nt⌋S_{\lfloor n t \rfloor} so that the limit is a standard Brownian motion?

Solution

σ=4=2\sigma = \sqrt{4} = 2, so divide by 2n2\sqrt{n}: W(n)(t)=S⌊nt⌋/(2n)W^{(n)}(t) = S_{\lfloor nt \rfloor}/(2\sqrt{n}). Dividing only by n\sqrt{n} would give a Brownian motion scaled by 22, i.e. 2Wt2W_t, whose variance at time tt is 4t4t rather than tt.

Path properties

Details

Step 1 — Continuity does not imply smoothness. A continuous function can still be nowhere differentiable (think of the Weierstrass function). Brownian paths are even rougher: the derivative dWt/dtdW_t/dt does not exist anywhere. Intuitively, zooming into any interval reveals the same jagged structure at every scale — the path is self-similar.

Step 2 — Self-similarity. For any c>0c > 0, the process c−1/2Wctc^{-1/2} W_{ct} is again a standard Brownian motion. This scaling symmetry is the continuous-time echo of the n\sqrt{n} scaling in the random-walk limit, and it is why the roughness is the same at every magnification.

Step 3 — Quadratic variation. For a smooth function ff, the sum ∑(fti+1−fti)2\sum (f_{t_{i+1}} - f_{t_i})^2 over a fine partition tends to 00 (the squared increments are of order (Δt)2(\Delta t)^2). For Brownian motion the increments are of order Δt\sqrt{\Delta t}, so the squared increments are of order Δt\Delta t, and the sum converges to the length of the interval:

[W]t=lim⁡∣Π∣→0∑i(Wti+1−Wti)2=ta.s.[W]_t = \lim_{|\Pi| \to 0} \sum_{i} (W_{t_{i+1}} - W_{t_i})^2 = t \quad \text{a.s.}

This is the single fact that makes Itô calculus different from ordinary calculus: the “infinitesimal” dWtdW_t satisfies (dWt)2=dt(dW_t)^2 = dt, not 00.

Example. Over [0,1][0,1], refine the partition from 1010 to 10001000 intervals. The sum of squared increments stays near 11 (the value of tt), not near 00 as it would for a smooth curve.

t = 1 10 20 50 100 200 500 1000 number of intervals in the partition of [0,1] Σ (ΔW)²

Authoritative sources

Example problem

For the smooth function f(t)=tf(t) = t on [0,1][0,1], what is lim⁡∣Π∣→0∑(fti+1−fti)2\lim_{|\Pi| \to 0} \sum (f_{t_{i+1}} - f_{t_i})^2? Contrast with [W]1[W]_1.

Solution

For f(t)=tf(t) = t, each increment is Δti\Delta t_i and ∑(Δti)2≤(max⁡Δti)∑Δti=max⁡Δti→0\sum (\Delta t_i)^2 \le (\max \Delta t_i) \sum \Delta t_i = \max \Delta t_i \to 0. So the quadratic variation of a smooth function is 00, whereas [W]1=1[W]_1 = 1. The nonzero quadratic variation is precisely what forces the extra 12f′′(Wt) dt\tfrac{1}{2} f''(W_t)\,dt term in Itô’s lemma.

Itô calculus and stochastic differential equations

Details

Step 1 — Why the ordinary chain rule fails. The chain rule df(Wt)=f′(Wt) dWtdf(W_t) = f'(W_t)\,dW_t is wrong for Brownian motion because WW has nonzero quadratic variation. The missing term is 12f′′(Wt) dt\tfrac{1}{2} f''(W_t)\,dt, which comes from the rule (dWt)2=dt(dW_t)^2 = dt.

Step 2 — Itô’s lemma. For f∈C2f \in C^2, a Taylor expansion f(Wt+Δt)−f(Wt)≈f′(Wt)ΔW+12f′′(Wt)(ΔW)2f(W_{t+\Delta t}) - f(W_t) \approx f'(W_t)\Delta W + \tfrac{1}{2} f''(W_t)(\Delta W)^2 plus the fact that (ΔW)2≈Δt(\Delta W)^2 \approx \Delta t (in the quadratic-variation sense) gives

df(Wt)=f′(Wt) dWt+12f′′(Wt) dt.df(W_t) = f'(W_t)\,dW_t + \tfrac{1}{2} f''(W_t)\,dt.

Example. Take f(x)=x2f(x) = x^2. Then d(Wt2)=2Wt dWt+dtd(W_t^2) = 2 W_t\,dW_t + dt. Integrating: Wt2=2∫0tWs dWs+tW_t^2 = 2\int_0^t W_s\,dW_s + t, so ∫0tWs dWs=12(Wt2−t)\int_0^t W_s\,dW_s = \tfrac{1}{2}(W_t^2 - t). The −t-t term is the Itô correction that has no classical analogue.

Step 3 — SDEs and geometric Brownian motion. An SDE dXt=μ dt+σ dWtdX_t = \mu\,dt + \sigma\,dW_t is solved by Xt=X0+μt+σWtX_t = X_0 + \mu t + \sigma W_t. The geometric SDE dSt=μSt dt+σSt dWtdS_t = \mu S_t\,dt + \sigma S_t\,dW_t keeps St>0S_t > 0 (multiplicative noise) and is the Black–Scholes price model. Applying Itô’s lemma to f(x)=ln⁡xf(x) = \ln x gives dln⁡St=(μ−12σ2) dt+σ dWtd\ln S_t = (\mu - \tfrac{1}{2}\sigma^2)\,dt + \sigma\,dW_t, so

St=S0exp⁡((μ−12σ2)t+σWt).S_t = S_0 \exp\big((\mu - \tfrac{1}{2}\sigma^2)t + \sigma W_t\big).

The −12σ2-\tfrac{1}{2}\sigma^2 in the exponent is again the Itô correction: it is why the expected price E[St]=S0eμt\mathbb{E}[S_t] = S_0 e^{\mu t} grows at rate μ\mu even though the log-drift is μ−12σ2\mu - \tfrac{1}{2}\sigma^2.

Authoritative sources

Example problem

Apply Itô’s lemma to f(x)=exf(x) = e^{x} to find d(eWt)d(e^{W_t}).

Solution

f′(x)=exf'(x) = e^x, f′′(x)=exf''(x) = e^x. So d(eWt)=eWt dWt+12eWt dt=eWt(dWt+12dt)d(e^{W_t}) = e^{W_t}\,dW_t + \tfrac{1}{2} e^{W_t}\,dt = e^{W_t}(dW_t + \tfrac{1}{2}dt). The 12dt\tfrac{1}{2}dt term is the Itô correction; without it the exponential would not be a martingale.

Quizz

Quizz

A stock price follows geometric Brownian motion dSt=μSt dt+σSt dWtdS_t = \mu S_t\,dt + \sigma S_t\,dW_t with S0=100S_0 = 100, μ=0.08\mu = 0.08, and σ=0.2\sigma = 0.2.

(a) Find E[S1]\mathbb{E}[S_1].

(b) Find the drift of ln⁡St\ln S_t.

Solution

(a) For geometric Brownian motion, E[St]=S0eμt\mathbb{E}[S_t] = S_0 e^{\mu t}, so E[S1]=100 e0.08≈108.33\mathbb{E}[S_1] = 100\,e^{0.08} \approx 108.33.

(b) Applying Itô’s lemma to f(x)=ln⁡xf(x) = \ln x gives dln⁡St=(μ−12σ2) dt+σ dWt=(0.08−0.02) dt+0.2 dWt=0.06 dt+0.2 dWtd\ln S_t = (\mu - \tfrac{1}{2}\sigma^2)\,dt + \sigma\,dW_t = (0.08 - 0.02)\,dt + 0.2\,dW_t = 0.06\,dt + 0.2\,dW_t. The drift of ln⁡St\ln S_t is 0.060.06 — smaller than μ\mu by the Itô correction 12σ2=0.02\tfrac{1}{2}\sigma^2 = 0.02.