TQBTHEQUANTBATEMAN
TQB/ learn/ foundations/ brownian motion ito calculusEN · DARK
Probability & measures · intermediate

Brownian motion & Itô calculus

Why continuous paths accumulate quadratic variation and change the chain rule.

BY THE END, YOU CAN

01State the defining increment properties of Brownian motion.

02Explain why quadratic variation converges to elapsed time.

03Construct the Itô integral for adapted elementary processes.

04Apply Itô's lemma and identify the second-order correction.

01
INTUITION

Observe the object before formalizing it.

Brownian motion is continuous but not differentiable. Its increments are order √dt, so squared increments are order dt and survive when ordinary second-order terms would disappear.

01

Independent Gaussian increments encode a time-homogeneous noise clock.

02

Quadratic variation—not visual roughness—is the mathematical signature that changes calculus.

03

Adapted integrands prevent future information entering a stochastic integral.

02
WHY MARKETS CARE

Connect the mathematical object to a pricing question.

Diffusion models, hedging arguments, simulation schemes and measure changes all depend on Brownian scaling and Itô's correction.

INSTRUMENTS

European and barrier options

stochastic-volatility models

short-rate models

dynamic hedges

QUOTE CONVENTION

Time is measured in years; Brownian increments have √year units, and diffusion coefficients restore the financial unit of the state variable.

03
MATHEMATICS

Construct the definition and its invariants.

Formula · Short derivation

Brownian increments

Wt−Ws∼N(0,t−s),0≤s<tW_t-W_s\sim\mathcal N(0,t-s),\qquad 0\le s<t

Increment variance equals elapsed time and disjoint increments are independent.

Formula · Full derivation

Quadratic variation

[W]T=lim⁡∥Πn∥→0∑i(Wti+1−Wti)2=T[W]_T=\lim_{\lVert\Pi_n\rVert\to0}\sum_i(W_{t_{i+1}}-W_{t_i})^2=T

The squared path increments converge to time even though the ordinary variation diverges.

Formula · Full derivation

Itô's lemma

df(t,Xt)=(ft+atfx+12bt2fxx)dt+btfxdWtdf(t,X_t)=\left(f_t+a_tf_x+\tfrac12b_t^2f_{xx}\right)dt+b_tf_xdW_t

The half-variance curvature term is the surviving second-order contribution.

Full derivation
Full derivation

Why stochastic calculus keeps a second-order term

Expand over a short interval and classify terms by Brownian order: ΔW is √Δt, so (ΔW)² is Δt.

  1. 01

    Write a second-order expansion

    For ΔX=aΔt+bΔW, retain the Taylor terms that can contribute at order Δt.

    Δf≈ftΔt+fxΔX+12fxx(ΔX)2\Delta f\approx f_t\Delta t+f_x\Delta X+\tfrac12f_{xx}(\Delta X)^2
  2. 02

    Apply Brownian scaling

    Terms in (Δt)² and ΔtΔW vanish faster than Δt; the squared Brownian increment does not.

    (ΔX)2=b2(ΔW)2+o(Δt)(\Delta X)^2=b^2(\Delta W)^2+o(\Delta t)
  3. 03

    Use quadratic variation

    Across refining partitions, the accumulated squared Brownian increments converge to elapsed time.

    (dWt)2=dt,dWtdt=(dt)2=0(dW_t)^2=dt,\qquad dW_tdt=(dt)^2=0
  4. 04

    Collect drift and diffusion

    Substitute the differential identities and separate finite-variation from martingale terms.

    df=(ft+afx+12b2fxx)dt+bfxdWdf=(f_t+af_x+\tfrac12b^2f_{xx})dt+bf_xdW

Itô's lemma is the ordinary chain rule plus the deterministic contribution created by quadratic variation.

Inputs
  • Wₜ: standard Brownian motion
  • Πₙ: partition of [0,T]
  • [W]ₜ: quadratic variation
  • Xₜ: Itô process with drift a and diffusion b
Assumptions and limits
  • Brownian paths exclude jumps and microstructure discontinuities.
  • Continuous-time identities become discretization choices in code.
  • An adapted integrand and square integrability are not optional technicalities.
05
MODEL / PRICING

Fit, compute, then challenge the assumptions.

METHOD

Use analytical scaling to classify terms, then simulate only to illustrate path behavior and convergence—not to define the theorem.

CALIBRATION

Standard Brownian motion has no calibration parameter; model diffusion coefficients and correlations are the calibrated objects.

06PYTHON IMPLEMENTATIONOpen the implementation and checks.
ARCHITECTURE
  • Seeded normal-increment generator
  • Path and quadratic-variation calculations kept separate
  • Convergence diagnostic over refining partitions
PYTHON 3 · NUMPY / SCIPY

Brownian path and quadratic variation

Show that path increments scale with √dt while their squared sum approaches T.

REUSABLE EXAMPLE
01from __future__ import annotations
02
03import numpy as np
04
05def brownian_path(steps: int, horizon: float, seed: int) -> tuple[np.ndarray, float]:
06 if steps <= 0 or horizon <= 0:
07 raise ValueError("positive steps and horizon required")
08 rng = np.random.default_rng(seed)
09 increments = np.sqrt(horizon/steps)*rng.standard_normal(steps)
10 path = np.concatenate(([0.0], np.cumsum(increments)))
11 return path, float(np.dot(increments, increments))
12
13path, qv = brownian_path(100_000, 1.0, 7)
14assert path.shape == (100_001,)
15assert abs(qv-1.0) < 0.02
16print(f"terminal={path[-1]:.6f} quadratic_variation={qv:.6f}")
EXPECTED OUTPUTterminal=<seeded value> quadratic_variation≈1.000000
SANITY CHECKS

✓ Path begins at zero

✓ Quadratic variation approaches the horizon

✓ Invalid grid inputs are rejected

07
INTERACTIVE LAB

Run the thought experiment.

NUMERICAL ERROR LAB

Brownian motion & Itô calculus

Use a fixed seed to separate discretization, sampling and truncation effects.

SYNTHETIC · EDUCATIONAL
Euler terminal error0.2987M=32
Milstein terminal error0.0083same shocks
Monte Carlo SE0.2108N=4,096
spot level by time (years)

Exact, Euler and Milstein spot paths using the same shocks.

  • exact
  • Euler
  • Milstein
time (years): 0.0Y. exact: 100.00. Euler: 100.00. Milstein: 100.00.

Use Left/Right or Up/Down arrows to inspect values; Home and End jump to the bounds.

ERROR BUDGET

Model → algorithm → evidence

Each numerical approximation has an independent control and a reference diagnostic.

01Target lawQ / GBM

Measure and SDE fixed

02Algorithm32 steps

Discretization or transform control

03DiagnosticSE=0.211

Visible convergence evidence

MODEL BOUNDARY

GBM and Gaussian transform examples are analytical references, not production calibration engines.

08
FRONT OFFICE

Carry the abstraction into valuation.

ON THE DESK
“The noise is not the model; the diffusion coefficient decides what one Brownian shock means in market units.”
VISIBLE INPUTS

time grid and calendar

diffusion and correlation

seed and random stream

discretization scheme

CALIBRATION

Calibrate model diffusion parameters to market instruments; validate the simulation scheme independently against moments or analytical prices.

RISK

discretization bias

correlation construction

barrier monitoring

DAILY WORKFLOW
  1. fix measure and SDE
  2. choose scheme and grid
  3. reuse shocks for comparisons
  4. verify moments and convergence
Production failure modes
  • unseeded regression tests
  • negative states from a naive scheme
  • mis-scaled annual time
09MACRO CONNECTIONOpen the transmission channel.
MACRO CONNECTION

From uncertainty clock to market dispersion

Brownian time is a modelling clock; volatility and correlation translate it into asset-specific uncertainty rather than economic causality.

01Information arrivaltransmits

sets a modelling clock

02Diffusiontransmits

scales shocks into returns

03Hedge gridoutput

leaves discretization residuals

10COMMON PITFALLSOpen the failure checklist.
01

Treating dW/dt as an ordinary derivative.

02

Dropping the Itô correction.

03

Using future information in the integrand.

04

Concluding model realism from a visually plausible path.

11SOURCES / FURTHER READINGOpen sources and continue the track.
research

Ch. 1 §§1.2–1.3, printed pp. 9–24 and Ch. 2 §2.1.2, printed pp. 29–34 (PDF pp. 28–43, 49–55)

The formulation was checked against the cited sections; prose, examples and code are original to the Academy.

Source
Mathematical Modeling and Computation in Finance
Author
Cornelis W. Oosterlee and Lech A. Grzelak
Ref
First edition (2020)
OPEN ORIGINAL SOURCE ↗