11.02.01 · stat-mech-physics / kinetic-theory

Maxwell-Boltzmann distribution from kinetic theory

draft3 tiersLean: nonepending prereqs

Anchor (Master): Landau & Lifshitz, *Statistical Physics*, Part 1, §4–5; Huang, *Statistical Mechanics*, 2e, Ch. 3

Intuition [Beginner]

In a gas at temperature $T$ , every molecule is in constant motion — but not at the same speed. Some drift slowly, some zip around. The Maxwell-Boltzmann distribution tells you what fraction of molecules have each speed. It looks like a lopsided bell curve: peaked near a most probable speed, with a long tail of fast movers.

Temperature controls the shape. Raise $T$ and the peak shifts right toward faster speeds and the curve flattens. Cool the gas and the peak shifts left toward slower speeds and sharpens.

Three characteristic speeds mark the distribution. The most probable speed is $v_{p} = 2 k T / m$ , the average speed is $v_{avg} = 8 k T / (π m)$ , and the root-mean-square speed is $v_{rms} = 3 k T / m$ , where $k$ is Boltzmann's constant and $m$ is the mass of one molecule.

For nitrogen ( $N_{2}$ , $m = 28 u$ ) at room temperature ( $300 K$ ): $v_{rms} \approx 517 m/s$ — faster than the speed of sound ( $343 m/s$ ). Sound propagates through air because the molecules are already moving faster than the sound wave; individual collisions carry the disturbance forward.

Visual [Beginner]

The Maxwell-Boltzmann speed distribution for three temperatures. At low $T$ the peak is tall, narrow, and leftward. As $T$ increases the peak shifts right, broadens, and drops. The three characteristic speeds ( $v_{p}$ , $v_{avg}$ , $v_{rms}$ ) appear as vertical lines, with $v_{p}$ leftmost and $v_{rms}$ rightmost.

The area under each curve between two speeds equals the fraction of molecules with speeds in that range. Because the total area must be 1 (every molecule has some speed), the curve at higher $T$ is both wider and shorter.

Worked example [Beginner]

Compute $v_{rms}$ for helium ( $^{4} He$ , $m = 4 u$ ) at $T = 300 K$ .

Step 1. Convert mass to kilograms. One atomic mass unit is $1 u = 1.66 \times 1 0^{- 27} kg$ , so $m = 4 \times 1.66 \times 1 0^{- 27} kg = 6.64 \times 1 0^{- 27} kg$ .

Step 2. Evaluate $3 k T / m$ with $k = 1.38 \times 1 0^{- 23} J/K$ :

\frac{3 k T}{m} = \frac{3 \times 1.38 \times 1 0 ^{- 23} \times 300}{6.64 \times 1 0 ^{- 27}} = \frac{1.242 \times 1 0 ^{- 20}}{6.64 \times 1 0 ^{- 27}} \approx 1.87 \times 1 0^{6} m^{2} / s^{2} .

Step 3. Take the square root:

v_{rms} = 1.87 \times 1 0^{6} \approx 1367 m/s .

Helium molecules at room temperature move at about $1367 m/s$ — more than twice as fast as nitrogen molecules ( $517 m/s$ ). The lighter the molecule, the faster it moves at the same temperature. This explains why helium balloons deflate quickly: the lightweight He atoms move so fast they slip through the microscopic pores in the rubber, escaping faster than heavier air molecules can diffuse in.

Check your understanding [Beginner]

Formal definition [Intermediate+]

The velocity distribution

Consider an ideal gas of $N$ identical particles of mass $m$ in thermal equilibrium at temperature $T$ . The Maxwell-Boltzmann velocity distribution is the probability density for the velocity components $(v_{x}, v_{y}, v_{z})$ of a single particle:

f (v_{x}, v_{y}, v_{z}) = (\frac{m}{2 π k T})^{3/2} exp (- \frac{m ( v _{x}^{2} + v _{y}^{2} + v _{z}^{2} )}{2 k T}) .

This factorises as $f (v_{x}) f (v_{y}) f (v_{z})$ where each component has the Gaussian form

f (v_{i}) = (\frac{m}{2 π k T})^{1/2} exp (- \frac{m v _{i}^{2}}{2 k T}) .

Derivation from symmetry. Assume the equilibrium velocity distribution of a gas in a container (i) depends on $v_{x}, v_{y}, v_{z}$ only through $v_{x}^{2}, v_{y}^{2}, v_{z}^{2}$ (isotropy), (ii) factorises as $g (v_{x}) g (v_{y}) g (v_{z})$ (independence of components), and (iii) is normalised. Write $g (v_{x}^{2}) = h (v_{x}^{2})$ . Independence of components and rotational invariance force $h$ to satisfy $h (a) h (b) h (c) = h (a + b + c)$ for all $a, b, c \geq 0$ , which (together with smoothness and normalisation) pins down $h (t) = C exp (- α t)$ for constants $C, α > 0$ . The constants follow from normalisation and matching the mean kinetic energy to $\frac{3}{2} k T$ .

The speed distribution

The speed distribution $f (v)$ (probability density for the magnitude $v = ∣ v ∣$ ) is obtained by integrating the velocity distribution over angular directions in velocity space. In spherical coordinates $(v, θ, φ)$ the solid-angle measure gives a geometric factor $4 π v^{2}$ (the surface area of a sphere of radius $v$ ), producing

f (v) = 4 π v^{2} (\frac{m}{2 π k T})^{3/2} exp (- \frac{m v ^{2}}{2 k T}), v \geq 0.

The competition between the geometric prefactor (growing as $v^{2}$ ) and the exponential (decaying as $e^{- m v^{2} /2 k T}$ ) produces the peak at $v_{p} = 2 k T / m$ . The moments of this distribution give the three characteristic speeds:

v_{p} = \frac{2 k T}{m}, ⟨ v ⟩ = \frac{8 k T}{π m}, v_{rms} = \frac{3 k T}{m} .

Equipartition theorem

The equipartition theorem states: in thermal equilibrium at temperature $T$ , each quadratic degree of freedom in the Hamiltonian contributes $\frac{1}{2} k T$ to the average energy.

For a monatomic ideal gas, the internal energy is purely translational kinetic energy with three quadratic terms ( $\frac{1}{2} m v_{x}^{2}$ , $\frac{1}{2} m v_{y}^{2}$ , $\frac{1}{2} m v_{z}^{2}$ ), so $U = \frac{3}{2} N k T$ . The molar specific heat at constant volume is $C_{V} = \frac{3}{2} R$ where $R = N_{A} k$ is the gas constant and $N_{A}$ is Avogadro's number.

For a diatomic gas (e.g. $N_{2}$ , $O_{2}$ ) at moderate temperatures, there are five quadratic degrees of freedom: three translational and two rotational (the axis along the bond has negligible moment of inertia). This gives $C_{V} = \frac{5}{2} R$ , in good agreement with experiment near room temperature. At higher temperatures vibrational modes unfreeze, adding two more quadratic terms (kinetic + potential) for $C_{V} = \frac{7}{2} R$ .

Mean free path and collision rate

The mean free path $λ$ — the average distance a molecule travels between collisions — for a gas of number density $n = N / V$ with molecular diameter $d$ is

λ = \frac{1}{2 π d ^{2} n} .

The collision rate (average number of collisions per molecule per unit time) is $z = ⟨ v ⟩ / λ$ . For air at standard conditions ( $n \approx 2.5 \times 1 0^{25} m^{- 3}$ , $d \approx 3.7 \times 1 0^{- 10} m$ ): $λ \approx 66 nm$ and $z \approx 7 \times 1 0^{9} s^{- 1}$ — each molecule suffers billions of collisions per second.

Counterexamples to common slips

The Maxwell-Boltzmann velocity distribution is a Gaussian (symmetric bell curve) in each component. The speed distribution is not Gaussian — the $4 π v^{2}$ prefactor skews it. Confusing the two is the most common error in first courses.
The equipartition theorem applies only to quadratic degrees of freedom. A particle in a non-harmonic potential $V (x) \propto x^{4}$ does not contribute $\frac{1}{2} k T$ from that coordinate.
The kinetic-theory derivation of the ideal gas law assumes non-interacting point particles. Real gases with intermolecular forces deviate; the van der Waals equation corrects for finite molecular volume and attractive forces.

Key theorem with proof [Intermediate+]

Theorem (Ideal gas law from kinetic theory). A gas of $N$ non-interacting point particles of mass $m$ in a container of volume $V$ at temperature $T$ exerts pressure $P$ satisfying $P V = N k T$ .

Proof. Consider a gas confined to a cubical box of side $L$ (volume $V = L^{3}$ ). A molecule with $x$ -component of velocity $v_{x}$ strikes the wall at $x = L$ with momentum $m v_{x}$ and rebounds with momentum $- m v_{x}$ , imparting impulse $2 m v_{x}$ to the wall per collision. The time between successive collisions with the same wall is $Δ t = 2 L / v_{x}$ (round trip to the opposite wall and back). The average force on this wall from one molecule's $x$ -motion is

F_{x} = \frac{2 m v _{x}}{Δ t} = \frac{m v _{x}^{2}}{L} .

Summing over all $N$ molecules and dividing by the wall area $L^{2}$ gives the pressure:

P = \frac{1}{L ^{2}} i = 1 \sum N \frac{m v _{x, i}^{2}}{L} = \frac{m}{L ^{3}} i = 1 \sum N v_{x, i}^{2} = \frac{N m ⟨ v _{x}^{2} ⟩}{V} .

By isotropy, $⟨ v_{x}^{2} ⟩ = ⟨ v_{y}^{2} ⟩ = ⟨ v_{z}^{2} ⟩ = \frac{1}{3} ⟨ v^{2} ⟩$ . Equipartition gives $\frac{1}{2} m ⟨ v^{2} ⟩ = \frac{3}{2} k T$ , so $⟨ v^{2} ⟩ = 3 k T / m$ . Substituting:

P = \frac{N m}{V} \cdot \frac{k T}{m} = \frac{N k T}{V} .

Rearranging: $P V = N k T$ . $□$

The result is independent of the box shape — a derivation using the stress tensor confirms this. It is also independent of the particle mass: heavier molecules move slower but carry more momentum per collision, and the two effects cancel exactly.

Corollary. Writing $N = n N_{A}$ where $n$ is the number of moles, $P V = n R T$ with $R = N_{A} k \approx 8.314 J/(mol K)$ . This is the familiar ideal gas law from thermodynamics 11.01.01 pending, now derived from microscopic mechanics rather than postulated.

Exercises [Intermediate+]

Exercise 7 (medium, symbolic).

A container holds a mixture of two ideal gases: $N_{1}$ molecules of mass $m_{1}$ and $N_{2}$ molecules of mass $m_{2}$ , at the same temperature $T$ . Derive the total pressure and show that the partial pressures satisfy Dalton's law: $P = P_{1} + P_{2}$ .

Hint

Apply the kinetic-theory pressure formula to each species separately. Each has its own $⟨ v^{2} ⟩$ but the same temperature.

Answer

For species $i$ : $P_{i} V = N_{i} m_{i} ⟨ v_{i}^{2} ⟩ /3$ . Equipartition gives $\frac{1}{2} m_{i} ⟨ v_{i}^{2} ⟩ = \frac{3}{2} k T$ regardless of mass, so $P_{i} V = N_{i} k T$ .

Total pressure: $(P_{1} + P_{2}) V = N_{1} k T + N_{2} k T = (N_{1} + N_{2}) k T$ . This is Dalton's law of partial pressures: the total pressure is the sum of the partial pressures, each given by $P_{i} = N_{i} k T / V$ , independent of the mass or identity of the other species.

Exercise 8 (hard, symbolic).

The fraction of molecules with speed exceeding a threshold $v_{0}$ is $F (v_{0}) = \int_{v_{0}}^{\infty} f (v) d v$ . Show that the fraction of molecules with kinetic energy exceeding $ϵ_{0} = \frac{1}{2} m v_{0}^{2}$ is

F (ϵ_{0}) = \frac{2}{π} u_{0} e^{- u_{0}^{2}} + erfc (u_{0}), u_{0} = \frac{ϵ _{0}}{k T},

where erfc is the complementary error function. Evaluate numerically for $ϵ_{0} = 5 k T$ (relevant to chemical kinetics).

Hint

Substitute $u^{2} = m v^{2} / (2 k T)$ and integrate by parts. Use $erfc (x) = \frac{2}{π} \int_{x}^{\infty} e^{- t^{2}} d t$ .

Answer

With $u = v m / (2 k T)$ , the fraction becomes $F = \frac{4}{π} \int_{u_{0}}^{\infty} u^{2} e^{- u^{2}} d u$ . Integrate by parts:

\int_{u_{0}}^{\infty} u^{2} e^{- u^{2}} d u = [- \frac{u}{2} e^{- u^{2}}]_{u_{0}}^{\infty} + \frac{1}{2} \int_{u_{0}}^{\infty} e^{- u^{2}} d u = \frac{u _{0}}{2} e^{- u_{0}^{2}} + \frac{π}{4} erfc (u_{0}) .

So $F = \frac{4}{π} (\frac{u _{0}}{2} e^{- u_{0}^{2}} + \frac{π}{4} erfc (u_{0})) = \frac{2 u _{0}}{π} e^{- u_{0}^{2}} + erfc (u_{0})$ , as claimed.

For $ϵ_{0} = 5 k T$ : $u_{0} = 5 \approx 2.236$ . Then $\frac{2 \times 2.236}{π} e^{- 5} \approx 0.0170$ and $erfc (2.236) \approx 0.0011$ , giving $F \approx 0.018$ . About 1.8% of molecules exceed $5 k T$ — this is the high-energy tail that drives activated chemical reactions.

Exercise 9 (hard, symbolic).

Derive the Maxwell-Boltzmann velocity distribution from scratch using only: (i) the velocity components $v_{x}, v_{y}, v_{z}$ are statistically independent; (ii) the distribution is isotropic — it depends only on $v^{2} = v_{x}^{2} + v_{y}^{2} + v_{z}^{2}$ ; (iii) the distribution is normalised. Show that these assumptions force a Gaussian.

Hint

Let $F (v) = g (v_{x}^{2}) g (v_{y}^{2}) g (v_{z}^{2})$ and demand $F$ depends on $v^{2}$ only. Take logarithms and use a functional equation.

Answer

Write $F = g (v_{x}^{2}) g (v_{y}^{2}) g (v_{z}^{2})$ . Isotropy requires $F$ to be a function of $v^{2}$ only. Taking logarithms: $ln F = ln g (v_{x}^{2}) + ln g (v_{y}^{2}) + ln g (v_{z}^{2}) = h (v^{2})$ for some function $h$ .

Set $v_{y} = v_{z} = 0$ : $h (v_{x}^{2}) = ln g (v_{x}^{2}) + 2 ln g (0)$ . So $ln g (a) = h (a) - 2 ln g (0)$ . Then $h (a + b + c) = h (a) + h (b) + h (c) - 2 c_{0}$ where $c_{0} = 2 ln g (0)$ . Define $H (a) = h (a) - c_{0}$ : $H (a + b + c) = H (a) + H (b) + H (c)$ — the Cauchy functional equation. Assuming smoothness, $H (a) = - α a$ for a constant $α$ . So $ln g (a) = - α a + const$ , giving $g (v_{i}^{2}) = A e^{- α v_{i}^{2}}$ .

Normalisation determines $A = (α / π)^{1/2}$ , and matching $⟨ \frac{1}{2} m v^{2} ⟩ = \frac{3}{2} k T$ gives $α = m / (2 k T)$ . This is the Maxwell 1860 derivation.

Exercise 10 (hard, symbolic).

A gas of particles of mass $m$ is in thermal equilibrium at temperature $T$ in a gravitational potential $V (z) = m g z$ near Earth's surface. Write down the full phase-space distribution $ρ (r, v)$ (the Boltzmann distribution). Then derive the barometric formula $n (z) = n (0) e^{- m g z / k T}$ for the number density as a function of height.

Hint

The distribution factorises: $ρ (r, v) \propto e^{- H (r, v) / k T}$ where $H = \frac{1}{2} m v^{2} + m g z$ . Integrate over velocities to get the spatial density.

Answer

The Boltzmann distribution in phase space is $ρ (r, v) \propto exp (- \frac{1}{2} m v^{2} / k T - m g z / k T)$ . This factorises as $ρ_{vel} (v) \cdot ρ_{spatial} (r)$ with $ρ_{spatial} \propto e^{- m g z / k T}$ .

For $V (z) = m g z$ : $n (z) = n (0) e^{- m g z / k T}$ . The velocity distribution at every height is the same Maxwell-Boltzmann distribution; only the density changes with altitude. At height $h = k T / (m g)$ (the scale height), the density drops by a factor of $e$ . For air at $300 K$ : $h \approx 8.6 km$ .

This is the barometric formula — the spatial part of the Maxwell-Boltzmann distribution in a gravitational field. It connects kinetic theory to the Boltzmann factor that underpins all of statistical mechanics 11.03.01 pending.

The Boltzmann equation and transport theory [Master]

The BBGKY hierarchy

The exact microscopic state of an $N$ -particle gas is described by the full $N$ -particle phase-space density $ρ^{(N)} (r_{1}, p_{1}, \dots, r_{N}, p_{N}, t)$ . The BBGKY hierarchy (Bogoliubov–Born–Green–Kirkwood–Yvon) is the tower of coupled equations obtained by integrating Liouville's equation over progressively more particles.

The one-particle distribution function $f_{1} (r, p, t)$ describes the probability of finding a particle near $(r, p)$ regardless of the other particles. Its equation of motion involves the two-particle correlation $f_{2}$ ; $f_{2}$ 's equation involves $f_{3}$ ; and so on. The hierarchy does not close on its own — each level requires knowledge of the next.

The Boltzmann equation

Boltzmann's insight (1872) was to close the BBGKY hierarchy at the first level by assuming molecular chaos (the Stosszahlansatz): before each binary collision, the two particles are uncorrelated, so $f_{2} (r_{1}, p_{1}, r_{2}, p_{2}) \approx f_{1} (r_{1}, p_{1}) f_{1} (r_{2}, p_{2})$ . This discards information about correlations and is the microscopic origin of the arrow of time in kinetic theory.

The resulting Boltzmann equation for $f (r, v, t)$ is

\frac{\partial f}{\partial t} + v \cdot \nabla_{r} f + \frac{F}{m} \cdot \nabla_{v} f = \frac{\partial f}{\partial t}_{coll} .

The left side is the streaming term (free transport plus external force). The right side is the collision integral, which for hard-sphere interactions reads

\frac{\partial f}{\partial t}_{coll} = \int d^{3} v_{2} \int d Ω ∣ v_{1} - v_{2} ∣ σ (Ω) [f (v_{1}^{'}) f (v_{2}^{'}) - f (v_{1}) f (v_{2})],

where the primed velocities are post-collision and $σ (Ω)$ is the differential cross section. The gain-loss structure makes the collision integral vanish when $f$ is the Maxwellian — the Maxwell-Boltzmann distribution is the unique stationary solution.

The H-theorem

Boltzmann defined the H-functional

H [f] (t) = \int d^{3} r \int d^{3} v f (r, v, t) ln f (r, v, t),

and proved $d H / d t \leq 0$ , with equality if and only if $f$ is the Maxwellian. The quantity $S = - k H$ is (up to additive constants) the Boltzmann entropy. The H-theorem is the first microscopic proof of the second law of thermodynamics: entropy increases monotonically toward the equilibrium (Maxwellian) state.

Loschmidt's and Zermelo's paradoxes

Loschmidt's reversibility paradox (1876). The microscopic equations of motion are time-reversible, yet the H-theorem claims irreversible behaviour. Resolution: the Stosszahlansatz is an irreversible assumption about initial conditions (uncorrelated incoming particles), not a consequence of reversible dynamics. Time-reversing a post-collision configuration violates molecular chaos, and the reversed state has measure-zero probability in a macroscopic gas.

Zermelo's recurrence paradox (1896). By the Poincare recurrence theorem, any bounded mechanical system returns arbitrarily close to its initial state given sufficient time. How can $H$ decrease monotonically? Resolution: Poincare recurrence times for macroscopic systems far exceed the age of the universe. The H-theorem holds on all physically relevant timescales.

Transport coefficients and the Chapman-Enskog expansion

Near equilibrium, the Boltzmann equation is solved perturbatively by writing $f = f_{MB} (1 + ϵ ϕ_{1} + \dots)$ where $f_{MB}$ is the local Maxwellian and $ϵ$ is proportional to the Knudsen number (mean free path divided by macroscopic length scale).

At first order, the Chapman-Enskog expansion yields Newton's law of viscosity ( $τ = η \nabla u$ ), Fourier's law of heat conduction ( $q = - κ \nabla T$ ), and Fick's law of diffusion ( $J = - D \nabla n$ ). The transport coefficients $η$ (viscosity), $κ$ (thermal conductivity), and $D$ (diffusion coefficient) are predicted as functions of molecular parameters, in remarkable agreement with experiment for dilute gases.

A key prediction: for a hard-sphere gas, viscosity $η \propto T$ — it increases with temperature, unlike liquids. This counterintuitive result arises because faster molecules transport momentum more effectively across shear layers, and it is confirmed by measurement.

Connections [Master]

The Maxwell-Boltzmann distribution is the hub connecting kinetic theory to the rest of physics and chemistry:

Statistical mechanics (Unit 11.03.01). The MB distribution is the equilibrium distribution of the microcanonical and canonical ensembles in the classical limit. It is the bridge from the microscopic postulates of statistical mechanics to the thermodynamic properties of ideal gases.
Chemical kinetics (Unit 14.08.01). The Arrhenius equation $k = A e^{- E_{a} / (R T)}$ derives from the fraction of the MB tail above the activation energy $E_{a}$ . The temperature dependence of reaction rates is a direct consequence of the MB speed distribution.
Biological energy metabolism (Unit 17.04.01). Enzyme-substrate collision frequencies in cells depend on molecular speeds governed by the MB distribution. Temperature regulation in homeotherms maintains optimal enzyme kinetics partly by keeping molecular speeds in a narrow range.
Quantum statistics (Unit 11.05.01). The MB distribution is the classical limit of both the Bose-Einstein and Fermi-Dirac distributions. The failure of classical MB statistics at low temperatures and high densities is what necessitated quantum statistics.
Fluid dynamics. The Navier-Stokes equations are the macroscopic limit of the Boltzmann equation. The Chapman-Enskog expansion makes this connection explicit.
Plasma physics. The MB distribution is the starting point for describing Maxwellian plasmas; deviations from Maxwellian (non-thermal tails, beam distributions) drive plasma instabilities.
Atmospheric science. The barometric formula (Exercise 10) and atmospheric escape (the Jeans escape mechanism, where the high-velocity tail of the MB distribution exceeds escape velocity) are direct applications.

Historical and philosophical context [Master]

Maxwell (1860). James Clerk Maxwell derived the velocity distribution in Illustrations of the Dynamical Theory of Gases (Philosophical Magazine, 1860). His derivation followed the symmetry argument reproduced in Exercise 9 — independence of components, isotropy, and the functional equation forcing a Gaussian. Maxwell was initially uncertain about the independence assumption; later derivations by Boltzmann put it on firmer ground.

Boltzmann (1872). Ludwig Boltzmann's H-theorem was the first demonstration that irreversible macroscopic behaviour (the second law) emerges from reversible microscopic dynamics — provided one makes the Stosszahlansatz. This was a revolutionary claim: time-asymmetry from time-symmetric laws.

Loschmidt (1876). Johann Josef Loschmidt's reversibility objection forced Boltzmann (and later the foundations-of-physics community) to confront the role of initial conditions and probability in the second law. The modern resolution — that entropy increase is statistical, not absolute, and that the arrow of time originates in the special low-entropy initial state of the universe — traces back to this debate.

Zermelo (1896). Ernst Zermelo's recurrence objection (drawing on Poincare's recurrence theorem) pushed Boltzmann toward a more explicitly probabilistic formulation, now understood as the statistical mechanics of the canonical ensemble.

The modern synthesis. The Ehrenfests (1911) clarified the distinction between the Boltzmann equation (valid for dilute hard-sphere gases under molecular chaos) and the coarser-grained Gibbs approach. The BBGKY hierarchy (1940s–1950s) placed the Boltzmann equation within a systematic framework. The Lanford theorem (1975) rigorously proved the validity of the Boltzmann equation for hard spheres in the Boltzmann-Grad limit (for short times), establishing the precise conditions under which the Stosszahlansatz is justified.

Bibliography [Master]

Maxwell, J. C. (1860). "Illustrations of the Dynamical Theory of Gases." Philosophical Magazine, 19, 19–32. The original derivation of the velocity distribution.
Boltzmann, L. (1872). "Weitere Studien über das Wärmegleichgewicht unter Gasmolekülen." Sitzungsberichte der Kaiserlichen Akademie der Wissenschaften, 66, 275–370. The H-theorem.
Loschmidt, J. (1876). "Über den Zustand des Wärmegleichgewichtes eines Systems von Körpern." Sitzungsberichte der Kaiserlichen Akademie der Wissenschaften, 73, 128–142. The reversibility objection.
Zermelo, E. (1896). "Über einen Satz der Dynamik und die mechanische Wärmetheorie." Annalen der Physik, 57, 485–494. The recurrence objection.
Chapman, S. & Cowling, T. G. (1970). The Mathematical Theory of Non-uniform Gases, 3rd ed. Cambridge University Press. The definitive treatment of the Chapman-Enskog expansion.
Cercignani, C., Illner, R. & Pulvirenti, M. (1994). The Mathematical Theory of Dilute Gases. Springer. Rigorous results including the Lanford theorem.
Huang, K. (1987). Statistical Mechanics, 2nd ed. Wiley. Ch. 3: careful treatment of the Boltzmann equation and transport theory.
Lifshitz, E. M. & Pitaevskii, L. P. (1981). Physical Kinetics. Course of Theoretical Physics Vol. 10, Pergamon. Comprehensive kinetic theory beyond the ideal gas.
Lanford, O. E. (1975). "Time Evolution of Large Classical Systems." In Dynamical Systems, Theory and Applications, Lecture Notes in Physics 38, Springer. Rigorous derivation of the Boltzmann equation from Hamiltonian dynamics.

Prerequisites

11.01.01 pending
09.01.03 pending
09.01.01 pending
02.12.01 pending

Used in

11.03.01

Tier anchors

beginner: Susskind, *The Theoretical Minimum: Statistical Mechanics* (2013), Lecture 4
intermediate: Schroeder, *Thermal Physics* (2000), Ch. 6
master: Landau & Lifshitz, *Statistical Physics*, Part 1, §4–5; Huang, *Statistical Mechanics*, 2e, Ch. 3

References

tong
raw/pdfs/statphys/statphys.pdf · §2 Kinetic theory — Maxwell-Boltzmann distribution, ideal gas, equipartition
TODO_REF
Schroeder — Thermal Physics (Addison-Wesley, 2000) · Ch. 6 Ideal Gas
TODO_REF
Landau & Lifshitz — Statistical Physics, Part 1, 3e (Pergamon, 1980) · §4 The Maxwellian distribution; §5 The Boltzmann distribution

Reviewer

Tyler (pending external stat-mech reviewer per PHYSICS_PLAN §6)

Estimated time

beginner: 15m
intermediate: 35m
master: 50m