Brownian Motion
In 1827 the botanist Robert Brown pointed his microscope at pollen grains suspended in water
and watched them jiggle — ceaselessly, erratically, with no current to blame and no sign of
life driving them. In 1900, in Paris, a doctoral student named Louis Bachelier put the
same jiggle to a very different use: he modelled the prices on the Paris Bourse with
it and derived the first mathematical theory of option pricing — five years before Einstein
wrote the physics paper that finally explained the pollen. Two worlds, one motion.
That motion, made rigorous, is Brownian motion — also called the
Wiener process and written (W_t)_{t \ge 0}. It is
randomness distilled to four axioms: the canonical continuous-time stochastic
process, the continuous cousin of the random walk, and the single object on which essentially
all of quantitative finance stands. The noise that drives stock prices in the Black–Scholes
world, the diffusions behind interest-rate models, the shocks in every Monte Carlo risk
engine — all are built from W_t. Learn this one process deeply and
the rest of the course is variations on it.
The four defining properties
Brownian motion is pinned down completely by four properties — nothing else is assumed, and
everything else we will ever prove about it follows from these:
-
Starts at zero: W_0 = 0.
-
Independent increments: over disjoint time intervals the increments
W_{t_2} - W_{t_1},\, W_{t_4} - W_{t_3}, \dots are mutually
independent — what the process does next does not depend on its past. This is the
mathematical form of "the market has no memory".
-
Gaussian stationary increments: for s < t,
W_t - W_s \sim N(0,\, t - s),
a normal
with mean 0 and variance equal to the elapsed time — the law
depends only on the gap t - s, not on where it sits. A week of
noise in March is statistically identical to a week of noise in November.
-
Continuous paths: almost every sample path
t \mapsto W_t(\omega) is continuous (no jumps).
A stochastic process (W_t)_{t \ge 0} is a
standard Brownian motion (Wiener process) if:
- W_0 = 0;
- for any 0 \le t_1 < t_2 \le t_3 < t_4, the increments
W_{t_2} - W_{t_1} and W_{t_4} - W_{t_3}
are independent (and likewise for any finite collection over disjoint intervals);
- for all 0 \le s < t,
W_t - W_s \sim N(0,\, t - s);
- the map t \mapsto W_t is continuous with probability one.
These four axioms determine the process uniquely (in law). Such a process exists — a
genuine theorem, not an assumption; see the construction vignette below.
Taking s = 0 gives the headline fact
W_t \sim N(0, t), so
\mathbb{E}[W_t] = 0 and
\operatorname{Var}(W_t) = t.
The spread grows like \sqrt{t}: uncertainty fans out with time —
but slower than time itself, a point we will make precise below.
One path of the noise
The axioms are austere; the object they describe is anything but. Below is a single sample
path of W_t over t \in [0,1], built as
a running sum of many tiny independent normal increments
\Delta W \sim N(0, \Delta t). Study its texture: it is continuous
— you can trace it without lifting your pen — yet violently jagged at every scale,
and it shows no preferred direction. It looks, not coincidentally, exactly like an intraday
price chart with the drift stripped out. Refresh to draw a fresh path — a different
\omega — and notice that while every path is different, the
statistics (the typical spread, the roughness) are always the same.
Worked examples: computing with the axioms
The four properties are not just a definition to memorise — they are a calculator. Three
computations every quant does in their head:
Example 1 — P(W_2 > 0). By the Gaussian-increments
property with s = 0, W_2 \sim N(0, 2): a
normal centred at zero. A centred normal is symmetric about 0, and
being a continuous distribution it puts zero probability on the single point
\{0\}, so the two halves split the probability evenly:
P(W_2 > 0) = \tfrac{1}{2}.
Nothing about the time 2 mattered: P(W_t > 0) = \tfrac{1}{2}
for every t > 0. Pure noise is equally likely to be above
or below where it started, at any horizon.
Example 2 — the distribution of W_3 - W_1. This is
an increment over the interval [1, 3], of length
3 - 1 = 2, so directly from the axiom
W_3 - W_1 \sim N(0, 2).
Note it has the same law as W_2 — that is stationarity of
increments in action: only the length of the window matters, not where it sits. Its standard
deviation is \sqrt{2} \approx 1.41.
Example 3 — P(W_4 > 2). Here we standardise.
W_4 \sim N(0, 4) has standard deviation
\sqrt{4} = 2, so Z = W_4 / 2 \sim N(0, 1)
and
P(W_4 > 2) = P(Z > 1) \approx 0.159.
A one-standard-deviation move — about a one-in-six event. Every "chance the price ends above
the strike" calculation in this course is a dressed-up version of this three-line computation.
Mean and variance, derived from the axioms
Let us not take those numbers on faith — every one of them falls straight out of the four
properties, with no skipped steps. Start with a single time t > 0 and
write W_t as an increment off the start:
W_t = W_t - W_0,
which is allowed because W_0 = 0 by the first property. The
increment W_t - W_0 runs over the interval
[0, t] of length t - 0 = t, so by the
Gaussian-increments property
W_t - W_0 \sim N(0,\, t - 0) = N(0, t).
Reading the mean and variance straight off this normal law:
\mathbb{E}[W_t] = 0,
\operatorname{Var}(W_t) = \operatorname{Var}(W_t - W_0) = t.
So the variance at time t is exactly t —
nothing more was needed than "starts at zero" plus "increments are
N(0, \text{gap})".
The covariance: \operatorname{Cov}(W_s, W_t) = \min(s, t)
Values at two different times are linked through the history they share. Fix two times and,
without loss of generality, label them so that s < t. The idea is to
split the later value into "what had happened by s" plus "what
happened afterwards", because those two parts are independent. Write
W_t = W_s + (W_t - W_s),
a plain algebraic identity. Now expand the covariance using its bilinearity (covariance is
linear in each argument):
\operatorname{Cov}(W_s, W_t) = \operatorname{Cov}\big(W_s,\, W_s + (W_t - W_s)\big),
= \operatorname{Cov}(W_s, W_s) + \operatorname{Cov}(W_s,\, W_t - W_s).
Take the two terms in turn. The first is just a variance,
\operatorname{Cov}(W_s, W_s) = \operatorname{Var}(W_s) = s,
by the result we just derived. For the second term, notice that
W_s = W_s - W_0 is the increment over [0, s]
and W_t - W_s is the increment over the disjoint interval
[s, t]. By the independent-increments property these
two increments are independent, and the covariance of independent random variables is zero:
\operatorname{Cov}(W_s,\, W_t - W_s) = 0.
Adding the two pieces back together,
\operatorname{Cov}(W_s, W_t) = s + 0 = s.
We assumed s < t, so s = \min(s, t); had we
labelled them the other way the same argument would give
t = \min(s, t). Either way the answer is the smaller of the two
times:
\operatorname{Cov}(W_s, W_t) = \min(s, t).
For a standard Brownian motion (W_t)_{t \ge 0} and any times
s, t \ge 0,
\mathbb{E}[W_t] = 0, \qquad \operatorname{Var}(W_t) = t, \qquad \operatorname{Cov}(W_s, W_t) = \min(s, t).
In particular W_t \sim N(0, t). These three facts follow from the
defining properties alone — W_0 = 0, independent increments, and
W_t - W_s \sim N(0, t - s) — and they characterise the
covariance structure of the process completely. (For a Gaussian process the mean and
covariance fix the law, so they fix Brownian motion itself.)
The four properties are easy to state, but it is not obvious that any process
satisfies all of them at once — continuous paths and independent Gaussian increments pull in
different directions. That a Brownian motion actually exists is a theorem, and there
are two famous ways to see it.
Direct construction (Lévy / Wiener). One can build W
explicitly on [0, 1] by the
Lévy–Ciesielski method: fill in the value at the midpoint with the correct
conditional Gaussian (a "Brownian bridge" interpolation), then the quarter-points, then the
eighths, and so on. The series converges uniformly with probability one, so the limit has
continuous paths by construction, and one checks it has the right Gaussian increments. This is
essentially Wiener's original existence proof, which is why the process bears his name.
Scaling limit of a random walk (Donsker). Take a simple symmetric random walk
S_n = \xi_1 + \cdots + \xi_n with
\pm 1 steps, speed it up n-fold in time and
shrink it by \sqrt{n} in space:
W^{(n)}_t = \frac{1}{\sqrt{n}}\, S_{\lfloor n t \rfloor}.
Donsker's invariance principle says these rescaled walks converge in
distribution (as processes on path space) to Brownian motion. The
1/\sqrt{n} is exactly the central-limit scaling, which is why the
limiting increments are Gaussian — Brownian motion is the continuous-time central limit
theorem made into a process, the precise sense in which it is "the continuous cousin of the
random walk".
The \sqrt{t} law: risk grows slower than time
\operatorname{Var}(W_t) = t looks innocent, but its square root is
one of the most consequential facts in finance. The typical size of the noise — its
standard deviation — is \sqrt{t}, so
quadrupling the horizon only doubles the spread:
\operatorname{sd}(W_{4t}) = \sqrt{4t} = 2\sqrt{t} = 2\,\operatorname{sd}(W_t).
To make the uncertainty 10\times bigger you must wait
100\times longer. Watch the law in the figure below: a bundle of
independent Brownian paths, fanning out inside the curves \pm\sqrt{t}
(one standard deviation — about 68\% of paths inside at any fixed
time) and \pm 2\sqrt{t} (about 95\%). The
envelope is a square-root, not a straight line — steep at first, then flattening.
This is exactly why traders quote volatility scaled by the square root of time.
An annualised volatility of 16\% corresponds to a daily volatility of
roughly 16\% / \sqrt{252} \approx 1\% (there are about
252 trading days in a year), and a "10-day value-at-risk" is the
1-day figure scaled by \sqrt{10} \approx 3.16, not by
10. Whenever you see a \sqrt{T} in a
pricing or risk formula, this axiom is where it came from.
Independent increments make variances add: split
[0, t] into n equal pieces and
W_t is the sum of n independent
increments, each of variance t/n, so
\operatorname{Var}(W_t) = n \cdot (t/n) = t. The
standard deviation is therefore \sqrt{t} — typical
displacement scales with the square root of time, the signature of diffusion. Doubling the
elapsed time only multiplies the typical wander by \sqrt{2}, not
by 2. If the spread grew like t instead,
the increments would have to conspire — always pushing the same way — which independence
forbids: half the shocks cancel, and cancellation is precisely what the square root measures.
Brownian motion is a martingale
A process is a martingale when its best forecast of any future value, given
everything seen so far, is simply its current value — the mathematical definition of a fair
game. Brownian motion is the prototype, and the proof is two lines of the same
split-the-increment trick we used for the covariance. Fix s < t and
let \mathcal{F}_s stand for the history of the path up to time
s. Split:
\mathbb{E}[W_t \mid \mathcal{F}_s] = \mathbb{E}[W_s + (W_t - W_s) \mid \mathcal{F}_s] = W_s + \mathbb{E}[W_t - W_s \mid \mathcal{F}_s].
The first term is W_s exactly because W_s
is already known at time s — conditioning on the history leaves it
alone. For the second term, the increment W_t - W_s is
independent of everything in \mathcal{F}_s (that is
the independent-increments axiom), so conditioning changes nothing and the conditional
expectation collapses to the plain one:
\mathbb{E}[W_t - W_s \mid \mathcal{F}_s] = \mathbb{E}[W_t - W_s] = 0,
since the increment is N(0, t-s), a centred normal. Therefore
\mathbb{E}[W_t \mid \mathcal{F}_s] = W_s.
Knowing the whole past of the noise gives you no edge in predicting where it goes
next — the expected future position is exactly where you stand. This fair-game property is the
seed of the deepest idea in the course: modern pricing theory says that (suitably discounted,
under the right probability measure) asset prices themselves must be
martingales,
or free money would exist. When you meet that principle, remember it started here.
Self-similarity: zoom in, same statistics
Take a Brownian motion, speed up its clock by a factor c > 0, and
shrink the result by \sqrt{c}:
B_t := \frac{1}{\sqrt{c}}\, W_{ct}.
Run the new process through the four axioms. It starts at
B_0 = W_0/\sqrt{c} = 0. Its increments over disjoint intervals come
from increments of W over disjoint intervals, so they are
independent. Each increment is Gaussian with
\operatorname{Var}(B_t - B_s) = \frac{1}{c}\operatorname{Var}(W_{ct} - W_{cs}) = \frac{ct - cs}{c} = t - s,
and its paths are continuous because W's are. All four axioms hold:
B is again a standard Brownian motion. This is
self-similarity: Brownian motion has no natural time scale. Zoom into any
stretch of the sample path above — milliseconds or decades — and after rescaling the axes by
the square-root law you see a curve with identical statistics. The path is a
statistical fractal, equally jagged at every magnification, which is why the figure looks like
a stock chart whether the axis reads minutes or years. (Real markets break this perfect
scaling — volatility clusters, jumps — and the Brownian benchmark is exactly what makes that
deviation measurable.)
Here is the trap: the paths are continuous, so it is tempting to treat
W_t like any calculus function and write
dW/dt. You cannot. With probability one, a Brownian
path is differentiable at no point whatsoever — continuous everywhere, differentiable
nowhere. The axioms show you why: the difference quotient over a small window
h is
\frac{W_{t+h} - W_t}{h} \sim N\!\left(0, \frac{h}{h^2}\right) = N\!\left(0, \frac{1}{h}\right),
whose variance blows up as h \to 0. The slope does not
settle down to a limit — it oscillates ever more wildly the closer you look. "Velocity of the
noise" simply does not exist.
Stranger still: over [0, t] the path's total variation
(the distance your pen travels tracing it) is infinite, yet its
quadratic variation — the sum of squared increments — is finite and
completely deterministic:
\sum_{i} \big(W_{t_{i+1}} - W_{t_i}\big)^2 \;\longrightarrow\; t \qquad \text{as the mesh} \to 0.
A random, nowhere-differentiable curve whose squared wiggles always add up to exactly the
elapsed time. This is the strangest fact on the page and the engine of everything later: it is
the rigorous meaning of the mnemonic (dW)^2 = dt, it is why ordinary
calculus fails on W, and it is why Itô had to invent a
new calculus — with a new chain rule carrying an extra second-order term — to
differentiate anything driven by Brownian motion. That invention, stochastic calculus, is
where this course goes next; the quadratic variation you just met is its beating heart.
A note on its history
The name honours the botanist Robert Brown,
who in 1827 watched pollen grains jitter erratically in water. The mathematics that turned that
jitter into a rigorous object was supplied by
Norbert Wiener in the 1920s — which is why
the process bears his name too. Remarkably, it had already been used to model the stock market
in 1900 by Louis Bachelier, five years
before Einstein's celebrated physical account — the first time chance was given a continuous,
time-evolving shape in finance.
Bachelier defended Théorie de la spéculation at the Sorbonne in 1900, with Henri
Poincaré on his jury. In it he modelled prices with Brownian motion, connected the process
to the heat equation, and priced options — a working mathematical theory of finance, from
scratch. The grade? Mention honorable — respectable, but short of the
très honorable an academic career required. The stock exchange was not considered a
fit subject for mathematics, and the thesis sank almost without trace.
The physicists caught up from the other side: in 1905
Einstein showed that Brown's jiggling
pollen was being bombarded by invisible molecules, derived the same
\sqrt{t} diffusion law, and turned the microscope into a tool
that let Jean Perrin count atoms — settling that atoms are real. Wiener made the
process rigorous in 1923. Bachelier stayed buried until the 1950s, when the statistician
Jimmie Savage mailed postcards to economists asking, in effect, "has anyone read this?"
Paul Samuelson had not — then did, was electrified, and put Bachelier's process
(exponentiated, to keep prices positive) at the centre of the theory that became
Black–Scholes. The man graded "honorable" had founded the field, half a century early.