1.3 Brownian motion as fluctuation

In 1827 the botanist Robert Brown watched, under a microscope, pollen grains suspended in water. The grains jittered erratically, as if alive — but the same jitter was visible in inanimate dust particles, in soot, in anything small enough. Brown could not explain it. Three-quarters of a century later, Einstein could: the grains were being kicked by the unseen molecules of the surrounding fluid, themselves in vigorous thermal motion of the kind we just visualised. Einstein’s 1905 paper on Brownian motion was the first concrete evidence that molecules are real. The probabilistic structure underlying this jitter — random walks, mean-square displacement growing as N\sqrt{N}, the diffusion equation as continuum limit — is developed in Foundations 11.3.

The setup

The simulation below sets the stage Brown actually saw. One large red particle of mass MM sits in a bath of about 200 small black particles, each of mass mMm \ll M. Every collision between the bath and the big particle is perfectly elastic and time-reversible — same physics as the kinetic-theory gas. But the mass asymmetry means each hit barely perturbs the big particle, and the motion you observe is the sum of many small random kicks.

one heavy particle (M = 20m) in a bath of 5000 light onesMean squared displacement vs. timeConverges to Einstein's relation ⟨|Δr|²⟩ = 4Dt01234 Mpx²0306090120Time (seconds)straight-line motion: r² = ⟨V²⟩ · t²average of 0 runsthis run

Each 120-second run of the heavy particle's (t) is noisy (light blue). Averaged over many runs, the noise washes out (black). If the particle moved in a straight line at the thermal speed ⟨V²⟩ = 2kT/M, the average would grow as — the green parabola. It does not. Because each collision randomises direction, the average grows linearly in t, with slope 4D. Linearity is the signature of randomness.

The plot on the right is the squared displacement Δr(t)2|\Delta \mathbf{r}(t)|^2 of the big particle from its starting position, as a function of elapsed time. Watch what it does over a few tens of seconds: it grows, on average, linearly in tt.

Einstein’s relation

Einstein’s 1905 analysis ties together three things the chapter has now met: the random thermal kicks of the bath, the drag a fluid exerts on a moving body, and the thermal energy kBTk_B T. The result, the Einstein relation, is built in three steps — define how fast the particle spreads, define how hard the fluid resists being pushed through, and show that the two are the same physics.

The diffusion coefficient

The squared-displacement plot grows linearly in time, and that linear growth is what defines the diffusion coefficient. Averaged over many runs (the angle brackets \langle\,\cdot\,\rangle denote that average), in dd dimensions

Δr(t)2  =  2dDt,\bigl\langle |\Delta \mathbf{r}(t)|^2 \bigr\rangle \;=\; 2 d\, D\, t,

so 11-D gives 2Dt2Dt, 22-D gives 4Dt4Dt (the case plotted), and 33-D gives 6Dt6Dt. The diffusion coefficient DD (units m2/s\text{m}^2/\text{s}) measures how fast the particle’s probability cloud spreads: the typical distance wandered grows as Δr2Dt\sqrt{\langle|\Delta\mathbf r|^2\rangle} \sim \sqrt{Dt} — the square-root-of-time signature of a random walk, slower than the linear-in-time progress of directed motion.

Why the mean-square displacement grows as 2dDt Derivation

The big particle performs a random walk: each molecular kick nudges it a small step in a random direction, uncorrelated with the previous one. For a sum of independent steps the variances add, so the mean-square displacement grows in proportion to the number of steps — and therefore in proportion to time. Defining the coefficient of that growth to be 2D2D in one dimension,

Δx2  =  2Dt.\langle \Delta x^2\rangle \;=\; 2 D t.

The full random-walk derivation of this linear law, and the diffusion equation it becomes in the continuum limit, are in Foundations 11.3. The motions along xx, yy, zz are independent and each contributes 2Dt2Dt, so in dd dimensions the squared displacements add:

Δr2  =  Δx2+Δy2+  =  2dDt.\bigl\langle |\Delta\mathbf r|^2\bigr\rangle \;=\; \langle\Delta x^2\rangle + \langle\Delta y^2\rangle + \cdots \;=\; 2 d\, D\, t.

The prefactor 2,4,62, 4, 6 is twice the number of dimensions. ✓

Drag: the same collisions, biased by motion

When the big particle sits still, the kicks are random and cancel on average; their leftover is the wander measured by DD. Push the particle through the fluid at velocity v\mathbf v, though, and the kicks stop cancelling: it runs into more molecules on its leading face than its trailing one, and the imbalance is a net drag force opposing the motion,

Fdrag  =  γv,\mathbf{F}_\text{drag} \;=\; -\gamma\,\mathbf{v},

where γ\gamma is the friction (drag) coefficient, units kg/s\text{kg/s}. (Its inverse 1/γ1/\gamma is the mobility, the drift speed a unit force produces.) For a sphere of radius aa in a fluid of dynamic viscosity η\eta, the slow-motion drag is Stokes’s law γ=6πηa\gamma = 6\pi\eta a.

Why the drag scales as γ = 6πηa Derivation

At low speed a small sphere’s drag can depend only on its radius aa, the fluid’s dynamic viscosity η\eta (units Pa⋅s=kg/(ms)\text{Pa·s} = \text{kg}/(\text{m}\,\text{s})), and its speed vv; the fluid’s inertia — its density — drops out when motion is slow enough. The drag must be linear in vv (doubling the speed doubles the rate at which momentum is handed to the fluid), so Fdrag=γvF_\text{drag} = \gamma v with γ\gamma assembled from η\eta and aa alone. Only one combination has the units of γ\gamma:

[ηa]  =  kgmsm  =  kgs  =  [γ],[\eta\,a] \;=\; \frac{\text{kg}}{\text{m}\,\text{s}}\cdot\text{m} \;=\; \frac{\text{kg}}{\text{s}} \;=\; [\gamma],

so γηa\gamma \propto \eta a. Solving the fluid equations at low Reynolds number fixes the constant at 6π6\pi, giving γ=6πηa\gamma = 6\pi\eta a; the solution is in Physics → Viscosity & diffusion. ✓

The relation

Diffusion (DD) and drag (γ\gamma) are the same molecular collisions seen two ways: random and self-cancelling, they make the particle fluctuate; biased by the particle’s motion, they dissipate it. Einstein’s result is that the two are locked together by the thermal energy kBTk_B T — with kBk_B Boltzmann’s constant, introduced in 1.2:

  D  =  kBTγ.  \boxed{\;D \;=\; \frac{k_B T}{\gamma}.\;}

This is the first instance of the fluctuation–dissipation principle: the size of a system’s spontaneous fluctuations (DD) is fixed by the friction that damps its driven motion (γ\gamma), with temperature setting the scale. A fluid that resists motion strongly lets its particles wander slowly, and vice versa.

The Einstein relation from a balance of drift and diffusion Derivation

Einstein’s argument is an equilibrium balance. Subject the Brownian particles to a weak external force — gravity on a colloid, say — derived from a potential U(x)U(x), so the force is F=dU/dxF = -\,dU/dx. Two particle currents result, and in equilibrium they must cancel everywhere.

Drift. The force drives each particle through the fluid, which resists with the drag γv-\gamma v. The particle quickly reaches the terminal velocity where force balances drag, F=γvdriftF = \gamma v_\text{drift}, so vdrift=F/γv_\text{drift} = F/\gamma. A population of number density n(x)n(x) drifting at this speed is a flux

Jdrift  =  nvdrift  =  nFγ.J_\text{drift} \;=\; n\,v_\text{drift} \;=\; \frac{nF}{\gamma}.

Diffusion. Random kicks carry particles from high to low concentration — Fick’s law — a flux

Jdiff  =  Ddndx.J_\text{diff} \;=\; -D\,\frac{dn}{dx}.

Equilibrium. With nothing flowing on average, the particles settle into the Boltzmann distribution in the potential (the equilibrium law, Physics → Kinetic theory),

n(x)  =  n0eU(x)/kBT,n(x) \;=\; n_0\, e^{-U(x)/k_B T},

and the total flux vanishes, Jdrift+Jdiff=0J_\text{drift} + J_\text{diff} = 0. Differentiating the Boltzmann distribution,

dndx  =  1kBTdUdxn  =  FkBTn,\frac{dn}{dx} \;=\; -\frac{1}{k_B T}\frac{dU}{dx}\,n \;=\; \frac{F}{k_B T}\,n,

using F=dU/dxF = -dU/dx. Put the two fluxes equal and opposite:

nFγ  =  Ddndx  =  DnFkBT.\frac{nF}{\gamma} \;=\; D\,\frac{dn}{dx} \;=\; D\,\frac{nF}{k_B T}.

The common factor nFnF cancels, leaving

D  =  kBTγ.  D \;=\; \frac{k_B T}{\gamma}. \;✓

The external force was only a scaffold to set up the balance; it has dropped out, so the relation is a property of the free particle in the bath.

For the spherical particle, substituting Stokes’s law γ=6πηa\gamma = 6\pi\eta a gives the form Einstein used to read molecular sizes off a microscope:

D  =  kBT6πηa.D \;=\; \frac{k_B T}{6\pi\eta a}.

Higher temperature raises DD (more energetic kicks); a larger or more viscously-embedded particle lowers it (more drag). In the simulation, slide TT up and watch the squared-displacement plot grow faster; slide the mass ratio up and watch it grow slower.

Why this matters for sound

Brownian motion is not sound. It is the equilibrium fluctuation of the medium itself — the residual jitter that remains in a fluid even when nothing macroscopic is happening. The bath kicks the big particle, but the kicks are incoherent: positive and negative, in random directions, uncorrelated from one moment to the next. They do not propagate; they do not carry information.

A sound, by contrast, is a coherent deviation from equilibrium. The kicks are organised: the molecules at one location all push in the same direction at the same time, then return, and then push the other direction, in sync. The disturbance has structure both in space and in time. It carries something — energy, information, the contents of a sentence — from where it was generated to wherever the air arrives.

We are about to spend a chapter learning to derive the equation that governs such coherent deviations. Before we do, the next lesson states clearly what a “sound” is, so we know what we are deriving an equation for.