Read · Brownian motion · Explanation preview

§3 · The diffusion of small suspended spheres

Read this section in the whole paper →

Einstein's letters change 1 of this paper's 18 formulas. 17 stay in today's letters, with a note saying so.

§3 · The diffusion of small suspended spheres

A derivation, to an approximation

Why viscosity changes the spread

What fixes D for a small spherical tracer in a liquid?

The 9 printed paragraphs this passage explains
  • §3, paragraph 1: Suspended particles are spread irregularly through a liquid, and each is acted on by a force K that depends on position but not on time; for simplicity the force points along the X axis everywhere.
  • §3, paragraph 2: In equilibrium the number of particles per unit volume varies along x so that the free energy does not change under any small virtual displacement; working out the energy and entropy changes gives equation (1), in which the force K is held in balance by osmotic pressure. The passage that explains §3 states this balance directly; it does not follow Einstein's route through the variation of the free energy.
  • §3, paragraph 3: Einstein uses equation (1) to find the diffusion coefficient, reading the equilibrium as two opposed processes: the particles drifting under the force K, and diffusion driven by their irregular thermal motion.
  • §3, paragraph 4: For spheres of radius P in a liquid of viscosity k, the force K gives each particle the speed K divided by 6πkP, so that many particles times that speed cross unit area in unit time.
  • §3, footnote 1: The footnote refers to Kirchhoff's lectures on mechanics, lecture 26, section 4, for the drag on a sphere moving slowly through a liquid.
  • §3, paragraph 5: With D the diffusion coefficient, diffusion carries particles across unit area at a rate D times the fall of their number per unit volume along x; in dynamic equilibrium the drift and the diffusion cancel, which is equation (2).
  • §3, paragraph 6: Equations (1) and (2) together give D as RT/N divided by 6πkP: apart from universal constants and the absolute temperature, the diffusion coefficient depends only on the liquid's viscosity and the size of the particles.
  • §5, paragraph 1: Combining the diffusion coefficient found in §3 for small spheres of radius P with the displacement law of §4 and eliminating D gives λx in terms of the time, the temperature, the viscosity k and the radius P.
  • §5, paragraph 2: Taking N as 6 · 10²³ from the kinetic theory of gases, water at 17 °C as the liquid, and particles 0.001 mm across, λx for one second comes to 8 · 10⁻⁵ cm, or 0.8 micron, and about 6 microns in one minute.

The displacement law tells us what a given DD predicts. A separate model connects DD to a tracer’s physical surroundings. Let bb be mobility, so a small force FF produces mean drift bFbF. Stokes drag for a sphere gives b=1/(6πηa)b = 1/(6\pi\eta a).

Let cc be number density. At isothermal balance the force density cFcF balances the osmotic-pressure gradient. With ideal osmotic pressure cRT/NcRT/N, the drift flux cbFcbF becomes b (RT/N)b\,(RT/N) times the density gradient. Equating it with the opposite diffusive flux gives D=bRT/ND = bRT/N.

D=R T6 π η a NAD = \frac{R\,T}{6\,\pi\,\eta\,a\,N_A}
D=kB T6 π η aD = \frac{k_B\,T}{6\,\pi\,\eta\,a}

Shown in modern letters; Einstein's are in §3 of the German source face.

Stokes–Einstein diffusivity is the gas constant times temperature, divided by six pi times viscosity times radius times Avogadro's number; equivalently, Boltzmann's constant times temperature, divided by six pi times viscosity times radius.

The equality kB = R/N relates the two forms. For a prediction using modern constants it is convenient. For an inference of N, using a value of kB derived from that same N would defeat the point.

Show every step here: Why viscosity changes the spread
  1. Name the quantities. TT is the absolute temperature, in kelvin. η\eta is the liquid’s dynamic viscosity, its resistance to flow, in pascal seconds. aa is the radius of the spherical tracer. RR is the molar gas constant, and NN is the number of molecules in a mole. Einstein prints the viscosity as kk and the radius as PP; his kk is the viscosity, not Boltzmann’s constant.
  2. DD is the diffusivity of the earlier passages: it says how fast the mean square grows, ⟨x2⟩=2Dt\langle x^2\rangle = 2Dt. This passage asks what fixes DD. The route has two halves: how a tracer answers a steady push, and how the crowding of many tracers pushes back.
  3. First half, the push. Suppose a small steady force FF acts on each tracer. In a viscous liquid the tracer soon moves at a steady mean drift speed in proportion to the force. Call the ratio of speed to force the mobility bb, so the drift speed is bFbF.
  4. Stokes’s law gives the drag on a sphere moving slowly at speed vv through a liquid as 6πηav6\pi\eta av. At steady drift the drag balances the push, F=6πηavF = 6\pi\eta av, so v=F/(6πηa)v = F/(6\pi\eta a) and the mobility is b=1/(6πηa)b = 1/(6\pi\eta a). Stokes’s law is brought in from the theory of fluids; the random-walk argument does not supply it.
  1. Second half, the crowding. Let cc be the number of tracers per unit volume, which Einstein writes ν\nu. Few, small tracers in a liquid exert an osmotic pressure, written here Π, of the same form as an ideal gas’s pressure: Π = cRT/N. Dividing RR by NN turns the gas constant per mole into a constant per molecule.
  2. Let the density vary along xx. A gradient is the rate at which a quantity changes with position. With TT held fixed, only cc varies in cRT/N, so the constant factor comes outside: the gradient of Π is RT/NRT/N times the gradient of cc.
  3. In a steady balance, the force on the tracers in a thin slab is held by the difference in osmotic pressure across it. Per unit volume that reads: cFcF equals the gradient of Π, which is RT/NRT/N times the gradient of cc. This osmotic law, like Stokes’s, is brought in from outside; conservation of tracers alone does not give it.
  1. Now count flows. A flux is the number of tracers crossing unit area in unit time. Tracers drifting at speed bFbF with density cc carry the drift flux c bFc\,bF. Multiply the balance of the last step by bb: the drift flux is b (RT/N)b\,(RT/N) times the gradient of cc.
  2. Fick’s law gives the diffusive flux: tracers spread from crowded to sparse regions, and the flux is −D-D times the gradient of cc, with the minus sign because the flow runs down the gradient.
  3. In the balance nothing flows overall, so the two fluxes add to zero: b (RT/N)b\,(RT/N) times the gradient, minus DD times the gradient, is zero. Divide by the gradient, which is not zero where the density varies: D=bRT/ND = bRT/N.
  4. Substitute the Stokes mobility b=1/(6πηa)b = 1/(6\pi\eta a): D=(RT/N)×1/(6πηa)D = (RT/N) \times 1/(6\pi\eta a), which is RTRT divided by 6πηaN6\pi\eta aN. Einstein found this in §3 and quotes it at the start of §5 as D=RTN16πkPD = \frac{RT}{N}\frac{1}{6\pi kP}, with his kk for η\eta and his PP for aa.
  5. Boltzmann’s constant kBk_B equals R/NR/N, the gas constant per molecule. Replacing R/NR/N by kBk_B gives the second form. Both are one relation:
D=R T6 π η a NAD = \frac{R\,T}{6\,\pi\,\eta\,a\,N_A}
D=kB T6 π η aD = \frac{k_B\,T}{6\,\pi\,\eta\,a}

Shown in modern letters; Einstein's are in §3 of the German source face.

Stokes–Einstein diffusivity is the gas constant times temperature, divided by six pi times viscosity times radius times Avogadro's number; equivalently, Boltzmann's constant times temperature, divided by six pi times viscosity times radius.

  1. Which form to use depends on the question. To predict a spread from modern constants, the kBk_B form is convenient. To infer NN from a spread, it defeats the point: a value of kBk_B derived from a known NN would put the answer into the data. The RR form keeps NN separate, because RR can be measured on gases without counting molecules.
  1. Einstein’s own numbers, from §5: water at 17 °C, so TT = 290.15 K; the viscosity printed as kk = 1,35⋅10−21{,}35\cdot10^{-2}, a bare number that in the centimetre-gram-second units of the time means poise, and is 1.35×10−31.35\times10^{-3} Pa·s; a particle diameter printed as 0,001 mm, so aa = 0.0005 mm = 5×10−75\times10^{-7} m; and NN = 6×10236\times10^{23} from the kinetic theory of gases. The paper does not print RR; this edition supplies RR = 8.31 J/(mol·K).
  2. The denominator: 6π × (1.35×10−31.35\times10^{-3} Pa·s) × (5×10−75\times10^{-7} m) = 1.272×10−81.272\times10^{-8} kg/s. The numerator: RT/NRT/N = 8.31 × 290.15 / (6×10236\times10^{23}) = 4.019×10−214.019\times10^{-21} J. Their quotient is DD = 3.16×10−133.16\times10^{-13} m2/s\mathrm{m}^2/\mathrm{s}, since a joule divided by a kilogram per second is a square metre per second.
  3. With λx=2Dt\lambda_x = \sqrt{2Dt} from §4: at tt = 1 s the mean square is 2Dt2Dt = 6.32×10−136.32\times10^{-13} m2\mathrm{m}^2, whose square root is 7.95×10−77.95\times10^{-7} m, about 0.8 μm. Einstein prints “8⋅10−58\cdot10^{-5} cm = 0,8 Mikron”. At tt = 60 s the mean square is 60 times larger, so λx\lambda_x is 60\sqrt{60}, about 7.75, times larger: 6.16 μm, which Einstein gives as about 6 micrometres in one minute.
  4. Now double the viscosity, holding TT, aa and the constants fixed. η\eta sits in the denominator, so DD halves, to 1.58×10−131.58\times10^{-13} m2/s\mathrm{m}^2/\mathrm{s}, and the mean square at a fixed time halves with it. The RMS is a square root, so it changes by a factor 1/21/\sqrt{2}, about 0.707: at 1 s it falls from 0.795 μm to 0.562 μm. Doubling the viscosity divides the RMS displacement by 2\sqrt{2}, not by two.
Explore the equations in this step

Read each operation, check its units and assumptions, or open the mathematical step behind it. These are modern teaching equations, not the equations as printed.

Explore the equation · Modern model notationEinstein's letters

The same diffusivity, written with the gas constant

D=R T6 π η a NAD = \frac{R\,T}{6\,\pi\,\eta\,a\,N_A}
D=R T6 π k P ND = \frac{R\,T}{6\,\pi\,k\,P\,N}

is times over six pi times , and .

is times over six pi times , and .

Select a term or operation. In the formula, Down enters an operation, Up returns to its parent, and Left/Right move between siblings. Escape clears selection. Tab leaves the formula.

Each term and operation, in words

No live values

This equation states a relation. No instrument on this page computes its terms.

Read the equation aloud in words

D equals R T over six pi eta a N sub A.

Einstein's form: the gas constant per molecule, R over N, stands where the modern form writes Boltzmann's constant. A larger N would make each molecule's kick smaller and the spreading slower.

Model assumptions and every term’s meaning
  • A dilute suspension of spheres, each moving slowly through a Newtonian liquid (Stokes drag).
  • Avogadro's number is the quantity the paper sets out to measure.
  • These are modern teaching equations in SI notation, not a transcription of the printed paper. The paper writes k for the viscosity.
Six pi eta a N sub A
Stokes's drag coefficient for one sphere, times Avogadro's number. Read the prerequisite
R times T
Thermal energy per mole. Read the prerequisite
Thermal energy over drag
Kicks drive spreading; drag resists it. Read the prerequisite
What it asserts
The spreading rate is fixed by temperature, viscosity, size and the number of molecules in a mole. Read the prerequisite
Avogadro's number
Molecules per mole: what the measurement is meant to find. Read the prerequisite
Diffusion coefficient
How fast the mean square displacement grows: half its rate of growth. Not a speed. Read the prerequisite
Gas constant
Energy per mole per kelvin. Read the prerequisite
Particle radius
The radius of the suspended sphere, not its diameter. Read the prerequisite
Absolute temperature
Temperature in kelvin. Read the prerequisite
Viscosity
The liquid's resistance to shear. The paper writes it k, which is not Boltzmann's constant. Read the prerequisite

Every term’s units were checked when this page was built. That checks the units, not the model. This equation is written for this edition in modern notation, not transcribed from the paper.

Explore the equation · Modern model notation

Resistance to motion controls spreading

D=kB T6 π η aD = \frac{k_B\,T}{6\,\pi\,\eta\,a}

Shown in modern letters; Einstein's are in §3 of the German source face.

is thermal energy divided by the viscous drag coefficient.

Select a term or operation. In the formula, Down enters an operation, Up returns to its parent, and Left/Right move between siblings. Escape clears selection. Tab leaves the formula.

Each term and operation, in words

No live values

Symbolic here. Open the linked laboratory for a worked example and live values.

Read the equation aloud in words

The diffusion coefficient equals the Boltzmann constant times the absolute temperature, divided by six times pi times the dynamic viscosity times the particle radius.

Thermal energy competes with viscous drag. Doubling viscosity halves the diffusion coefficient while reducing RMS displacement only by the square root of two. Radius means radius, not diameter.

Model assumptions and every term’s meaning
  • Dilute, approximately spherical tracers with no-slip Stokes drag in a homogeneous Newtonian liquid.
  • Wall corrections, interactions, inertia, and observation noise are not included.
  • The modern SI 2019 constant set is used; this is not a historical inversion exercise.
An ideal model
This equation assumes the dilute spherical tracer and Stokes-drag model; dimensional consistency alone does not establish those assumptions. Read the prerequisite
Diffusion coefficient
A squared-distance-per-time coefficient. This value comes from the accepted model calculation. Read the prerequisite
Why divide by drag?
Increasing the denominator reduces diffusivity at fixed thermal energy. This is a model dependence, not a rule that every larger quantity must reduce motion. Read the prerequisite
Thermal energy scale
Boltzmann constant times absolute temperature supplies energy per particle. Read the prerequisite
A known modern constant
The SI 2019 constant set supplies this value. It must not be treated as independent historical evidence when trying to infer molecular number. Read the prerequisite
Absolute temperature
Use kelvin. The laboratory holds viscosity as an independently supplied parameter. Read the prerequisite
The drag coefficient
Six pi times viscosity times radius is Stokes drag per speed for an isolated sphere in the assumed low-Reynolds regime. Read the prerequisite
Dynamic viscosity
Viscosity is in the denominator. Doubling it at fixed temperature and radius halves D, not the RMS displacement. Read the prerequisite
Particle radius
The input is a radius, not a diameter. Confusing them changes the predicted coefficient by a factor of two. Read the prerequisite

Every term’s units were checked when this page was built. That checks the units, not the model. This equation is written for this edition in modern notation, not transcribed from the paper.

Modern qualifications

Einstein reaches this relation in §§2–3 through osmotic pressure and Stokes’s law, writing kk for the viscosity and PP for the radius: D=RTN16πkPD = \frac{RT}{N}\frac{1}{6\pi kP}. His kk is not Boltzmann’s constant. William Sutherland published the same relation independently in 1905, and Marian Smoluchowski reached the displacement law by another route in 1906. The laboratories use modern SI constants and say so.

Assumptions and limits: Why viscosity changes the spread

Assumed here

  • Few, small, spherical particles in a uniform liquid, with Stokes’s drag and ideal osmotic pressure.
  • R, the gas constant, is known separately; N is the number of molecules in a mole.

What this does not establish

  • Slip at the particle’s surface, its inertia, interactions between particles, unusual liquids, and gases are all left out.
  • Stokes’s law and the osmotic-pressure law are brought in from outside; conservation alone does not give them.

Earlier step: What the spreading curve predicts

Source context: Annalen der Physik (4), 17, 549–560 (1905), §§1–5. Bibliographic pointer; this preview is not a source transcription or translation.

Explanation

Einstein's displacement argument: the typical distance grows with the square root of time, not with time itself. Every explanation stays on the page when reading-only is on.

Static worked case

For radius 0.5 μm, viscosity 1.35×10⁻³ Pa·s, and T = 290.15 K, the RMS displacement is about 0.8 μm in one second. This static worked case stays in the markup; loading the live ensemble does not replace it.

Try the displacement argument yourself

Once open, the tracer ensemble stays open while you change the detail or read a lesson. It runs a new trial only when you apply its settings; opening an explanation never starts one.

Open the tracer ensemble in this reading

Open the tracer ensemble on its own page, where its worked example is part of the page.

Notes and laboratory
Page 554 of Annalen der Physik, volume 17, where §3 begins.
§3 begins on page 554 of Annalen der Physik, volume 17. Read it in Einstein’s German

Symbols in §3

  • Diffusion coefficient
  • Gas constant
  • Absolute temperature
  • Dynamic viscosity
  • Particle radius
  • Avogadro's number
  • Boltzmann constant

Keep the tracer ensemble in view

Read in another form

Download the full explanation as Markdown · The same explanation as data (JSON) · Browse the foundation library

The question we were answering:

Read this offline

Each chapter is one HTML file that opens in any browser without a connection. It holds the explanation at every level of detail, the lessons it links to, its source references and the numbers worked out when the site was built. It has no notes of yours and no running simulations. It is an explanation, not an edition of the paper.

Links to online sources still need a connection, and the interactive plots are not included. Every chapter you can keep