---
title: Classical Statistics and Equipartition
module: Microstates, Phase Space, and Statistical Entropy
moduleNumber: 2
lessonNumber: 1
order: 201
summary: >
  A liter of gas holds on the order of a trillion trillion molecules, far too
  many to track by their equations of motion. Classical statistical mechanics
  replaces the trajectories with a single probability law, the Boltzmann
  distribution, and reads the measurable properties of matter off it: the
  Maxwell speed distribution, the average energy per degree of freedom, and the
  heat capacities of gases and solids — together with the low-temperature
  failures that forced the quantum revision.
topics: ["Microstates, Phase Space, and Statistical Entropy"]
sources:
  - book: Tipler & Llewellyn
    ref: "Ch. 8 Statistical Physics; §8-1 Classical Statistics: A Review, Temperature and Entropy, A Derivation of the Equipartition Theorem"
  - book: Tipler & Mosca
    ref: "Ch. 17 Temperature and Kinetic Theory; §17-6 The Kinetic Theory of Gases"
draft: false
---

A macroscopic system is assembled from a number of atoms of order $10^{22}$.
Classically, its future is fixed by the equations of motion for every
constituent given the state of motion at one instant, but solving $10^{22}$
coupled equations is not a program anyone can run. The measurable properties of
bulk matter can nonetheless be predicted, because conservation of energy and
momentum applied to a large ensemble fixes the _probable_ behavior of the
system, and the probable behavior is what a measurement reports.[^tl-intro] This
is **statistical mechanics**: the observable properties of a system are averages
over a probability distribution of microscopic states, and the distribution,
not the trajectories, is the object of study.

The central task is to find how a fixed amount of energy distributes itself
among the particles of a system in thermal equilibrium. The particles exchange
energy through collisions, so any one particle's energy fluctuates above and
below the mean over time. Classical statistical mechanics asserts that the
values the energy takes on follow one specific probability law, the Boltzmann
distribution, and that from it the properties of the whole system follow.

## The Boltzmann distribution

Classically the system is a large ensemble of identical but **distinguishable**
particles: all alike, yet in principle trackable through a collision, like
billiard balls with numbers painted on them. Boltzmann derived the law giving
the probable number of particles occupying each available energy state in such
an ensemble at equilibrium.

> **Definition (Boltzmann distribution).** The probability that a state of
> energy $E$ is occupied at absolute temperature $T$ is
> $$
> f_B(E) = A\,e^{-E/kT},
> $$
> where $A$ is a normalization constant set by the particular system and $k$ is
> the Boltzmann constant. The exponential $e^{-E/kT}$ is the **Boltzmann
> factor**.[^tl-boltz]

The Boltzmann constant is the gas constant per molecule,

$$
k = 1.381 \times 10^{-23}\ \mathrm{J/K} = 8.617 \times 10^{-5}\ \mathrm{eV/K},
$$

so that at room temperature $kT \approx 0.026\ \mathrm{eV}$, a number worth
memorizing: it is the energy scale that separates "easily excited" from
"frozen out." A state costing much more than $kT$ above the ground state is
almost never occupied; a state within $kT$ is populated freely.

A single state is rarely the object of interest. The number of _distinct_
states at energy $E$ — the **statistical weight** or **degeneracy** $g(E)$ —
multiplies the occupation probability, so the number of particles with energy
$E$ is

$$
n(E) = g(E)\,f_B(E) = A\,g(E)\,e^{-E/kT}.
$$

When the spectrum is dense, $E$ is treated as continuous and $g(E)$ becomes the
**density of states**, defined so that $g(E)\,\d E$ counts the states with energy
between $E$ and $E + \d E$. The exponential decay of $f_B$ with energy is the
single most important feature: it is drawn below as the falling population of a
ladder of levels.

$$
% caption: In equilibrium the number of particles on each level falls by the
% same Boltzmann factor with every step up in energy; the population of a level
% is set by its energy measured in units of kT.
\begin{tikzpicture}[scale=1.0, font=\small]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,-0.3) -- (0,4.4) node[left, font=\footnotesize] {energy};
\foreach \y/\len in {0/3.0, 0.9/1.82, 1.8/1.10, 2.7/0.67, 3.6/0.40} {
  \draw[black, thick] (-0.05,\y) -- (0.28,\y);
  \fill[acc!14, draw=acc] (0.4,\y-0.12) rectangle (0.4+\len,\y+0.12);
}
\node[anchor=west, font=\footnotesize] at (3.55,0) {ground level, most populated};
\node[anchor=west, font=\footnotesize] at (1.35,3.6) {sparsely populated};
\node[anchor=west, font=\footnotesize, text=acc] at (0.9,2.15) {bar length = population};
\end{tikzpicture}
$$

Two examples show the distribution at work.

**The law of atmospheres.** For an ideal gas in a uniform gravitational field,
a molecule of mass $m$ at height $z$ has energy $E = p^2/2m + mgz$. Integrating
the Boltzmann factor over the momentum components leaves the height dependence
alone,

$$
f_B(z)\,\d z = A\,e^{-mgz/kT}\,\d z,
$$

so the density falls off exponentially with altitude. For air ($M = 28.6$) at
$T = 300\ \mathrm{K}$, the ratio of densities at $1\ \mathrm{km}$ and at ground
is $e^{-mg(1000)/kT} = e^{-0.113} = 0.893$, giving
$\rho(1000) = 1.154\ \mathrm{kg/m^3}$ from a ground value of
$1.292\ \mathrm{kg/m^3}$.[^tl-atmos]

**Populations of atomic levels.** The ratio of populations of two levels is the
ratio of Boltzmann factors weighted by degeneracies,

$$
\frac{n_2}{n_1} = \frac{g_2}{g_1}\,e^{-(E_2 - E_1)/kT}.
$$

For the first excited state of [hydrogen](/atomic-physics/early-models-and-old-quantum-theory/bohr-model-hydrogen),
$E_2 - E_1 = 10.2\ \mathrm{eV}$ with $g_2/g_1 = 4$. At room temperature
$kT \approx 0.026\ \mathrm{eV}$ and

$$
\frac{n_2}{n_1} = 4\,e^{-10.2/0.026} = 4\,e^{-392} \approx 0,
$$

which is why a jar of hydrogen at room temperature does not glow with the
Balmer series. At the Sun's surface, $T = 5800\ \mathrm{K}$ and
$kT \approx 0.50\ \mathrm{eV}$, giving $n_2/n_1 \approx 4\,e^{-20.4} \approx
10^{-8}$: about $10^{15}$ atoms per mole sit in the first excited state at any
instant, enough to produce the observed absorption lines.

## Temperature and entropy

Temperature acquires a precise microscopic meaning through the number of
microstates. A **macrostate** (fixed total energy, volume, particle number) is
realized by many **microstates** — detailed assignments of energy to individual
particles. The number of microstates $W$ consistent with a macrostate is its
multiplicity, and Boltzmann's relation defines the entropy as its logarithm.

> **Definition (Boltzmann entropy).** The entropy of a macrostate with
> multiplicity $W$ is
> $$
> S = k \ln W.
> $$
> Equilibrium is the macrostate of largest $W$, hence of largest $S$: the system
> is found, overwhelmingly, in the arrangement that the most microstates
> realize.[^tl-entropy]

Absolute temperature is then the rate at which entropy responds to added energy
at fixed volume,

$$
\frac{1}{T} = \frac{\partial S}{\partial U}.
$$

A low-entropy macrostate concentrates the energy in a few particles; the
high-entropy macrostate spreads it out. There are vastly more ways to spread
energy than to concentrate it, so an isolated system drifts toward the spread
arrangement, and this drift is the second law.

$$
% caption: The same total energy distributed among particles: concentrating it
% on one particle is realized by few microstates, spreading it among many is
% realized by overwhelmingly more, so equilibrium is the spread arrangement.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
% left panel: concentrated
\draw[black] (-0.2,-0.2) rectangle (3.0,2.4);
\node[font=\footnotesize, anchor=south] at (1.4,2.45) {concentrated};
\foreach \x/\y in {0.4/0.3, 1.0/0.3, 1.6/0.3, 2.2/0.3, 0.7/0.9, 1.3/0.9, 1.9/0.9} \fill[black] (\x,\y) circle (2pt);
\fill[acc] (1.3,1.7) circle (2pt);
\draw[acc, thick] (1.3,1.7) circle (5pt);
\node[text=acc, anchor=west, font=\scriptsize] at (1.55,1.7) {all energy here};
\node[anchor=north, font=\scriptsize] at (1.4,-0.25) {small $W$, low $S$};
% right panel: spread
\begin{scope}[xshift=4.4cm]
\draw[black] (-0.2,-0.2) rectangle (3.0,2.4);
\node[font=\footnotesize, anchor=south] at (1.4,2.45) {spread};
\foreach \x/\y in {0.4/0.4, 1.0/0.4, 1.6/0.4, 2.2/0.4, 0.7/1.1, 1.3/1.1, 1.9/1.1, 0.5/1.8, 1.1/1.8, 1.7/1.8, 2.3/1.8} {\fill[black] (\x,\y) circle (2pt); \draw[black] (\x,\y) circle (3.4pt);}
\node[anchor=north, font=\scriptsize] at (1.4,-0.25) {large $W$, high $S$};
\end{scope}
\draw[->, acc, very thick] (3.25,1.1) -- (4.15,1.1);
\end{tikzpicture}
$$

## The Maxwell speed distribution

The Boltzmann distribution applied to the kinetic energy of gas molecules
reproduces Maxwell's velocity and speed distributions, derived a few years
before Boltzmann's general law. Assuming the three velocity components are
independent and the distribution depends only on speed, the normalized velocity
distribution is a product of three Gaussians,

$$
F(v_x, v_y, v_z) = \left(\frac{m}{2\pi kT}\right)^{3/2}
  e^{-m(v_x^2 + v_y^2 + v_z^2)/2kT}.
$$

Each component is symmetric about zero, so $\langle v_x \rangle = 0$: on average
a gas goes nowhere. The distribution of _speeds_ $v = |\mathbf{v}|$ is obtained
by multiplying the density in velocity space by the volume $4\pi v^2\,\d v$ of the
spherical shell of radius $v$,

> **Theorem (Maxwell speed distribution).** The number of molecules with speed
> between $v$ and $v + \d v$ in a gas of $N$ molecules at temperature $T$ is
> $$
> n(v)\,\d v = 4\pi N\left(\frac{m}{2\pi kT}\right)^{3/2} v^2\,e^{-mv^2/2kT}\,\d v.
> $$

The $v^2$ factor pulls $n(v)$ to zero at the origin and the exponential pulls
it to zero at large $v$, so the distribution peaks at an intermediate speed and
is asymmetric, with a long high-speed tail.

$$
% caption: The Maxwell speed distribution at two temperatures. Heating broadens
% and flattens the curve and shifts its peak to higher speed, while the area
% (the total number of molecules) is unchanged.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,0) -- (7.0,0) node[right] {speed};
\draw[->, black] (0,0) -- (0,3.0) node[above] {n(v)};
\draw[acc, thick, smooth] plot coordinates
  {(0,0)(0.5,0.55)(1,1.5)(1.5,2.25)(2,2.5)(2.5,2.3)(3,1.85)(3.5,1.35)(4,0.92)(4.5,0.6)(5,0.38)(5.5,0.23)(6,0.14)(6.5,0.08)};
\draw[black, thick, dashed, smooth] plot coordinates
  {(0,0)(0.5,0.28)(1,0.75)(1.5,1.2)(2,1.5)(2.5,1.68)(3,1.72)(3.5,1.66)(4,1.5)(4.5,1.3)(5,1.08)(5.5,0.86)(6,0.66)(6.5,0.49)};
\node[text=acc, anchor=west] at (2.0,2.55) {cooler};
\node[text=black, anchor=west] at (4.6,1.5) {hotter};
\end{tikzpicture}
$$

Three characteristic speeds summarize the curve. The **most probable speed**
$v_m$ is at the peak, the **average speed** $\langle v \rangle$ is the mean, and
the **rms speed** $v_\text{rms}$ is the square root of the mean square, tied to
the average kinetic energy. They are computed from the distribution by
integration and ordered $v_m < \langle v \rangle < v_\text{rms}$.

| Speed | Value | Definition |
| --- | --- | --- |
| Most probable $v_m$ | $\sqrt{2kT/m}$ | maximum of $n(v)$ |
| Average $\langle v \rangle$ | $\sqrt{8kT/\pi m}$ | $\tfrac{1}{N}\int v\,n(v)\,\d v$ |
| Rms $v_\text{rms}$ | $\sqrt{3kT/m}$ | $\sqrt{\langle v^2 \rangle}$ |

All three scale as $\sqrt{T/m}$: heavier molecules move slower at the same
temperature, and any gas speeds up as $\sqrt{T}$. For nitrogen
($m = 4.68 \times 10^{-26}\ \mathrm{kg}$) at $300\ \mathrm{K}$,

$$
\langle v \rangle = \sqrt{\frac{8kT}{\pi m}}
= \left[\frac{8(1.38\times 10^{-23})(300)}{\pi(4.68\times 10^{-26})}\right]^{1/2}
= 475\ \mathrm{m/s} \approx 1700\ \mathrm{km/h},
$$

about eight percent below $v_\text{rms}$.

## The distribution of kinetic energy

Changing variables from speed to kinetic energy $E = \tfrac{1}{2}mv^2$ in the
Maxwell distribution gives the distribution of molecular kinetic energy,

$$
n(E)\,\d E = \frac{2N}{\sqrt{\pi}\,(kT)^{3/2}}\,E^{1/2}\,e^{-E/kT}\,\d E,
$$

where the $E^{1/2}$ is the density of states carried over from the $v^2$ shell
factor. Multiplying by $E$ and integrating gives the average translational
kinetic energy,

$$
\langle E \rangle = \frac{1}{N}\int_0^\infty E\,n(E)\,\d E = \frac{3}{2}kT,
$$

independent of the molecular mass. This is the microscopic content of
temperature: absolute temperature is a direct measure of the mean translational
kinetic energy per molecule.

$$
% caption: Maxwell distribution of molecular kinetic energy. The rising density
% of states and the falling Boltzmann factor combine to a peaked curve whose
% mean sits at three-halves kT.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,0) -- (6.4,0) node[right] {E};
\draw[->, black] (0,0) -- (0,2.9) node[above] {n(E)};
\draw[acc, thick, smooth] plot coordinates
  {(0,0)(0.35,1.35)(0.7,1.7)(1.0,1.75)(1.5,1.55)(2,1.25)(2.5,0.95)(3,0.7)(3.5,0.5)(4,0.35)(4.5,0.24)(5,0.16)(5.5,0.1)(6,0.06)};
\draw[black, dashed] (1.5,0) -- (1.5,1.55);
\node[anchor=south, font=\scriptsize] at (1.5,1.58) {mean at $\tfrac{3}{2}kT$};
\node[anchor=north, font=\scriptsize] at (1.5,-0.05) {$\tfrac{3}{2}kT$};
\end{tikzpicture}
$$

That $\langle E \rangle$ does not depend on mass has a consequence for planetary
atmospheres. A gas escapes a planet over $\sim 10^8$ years if its average speed
reaches about one-sixth of the escape speed. At $T = 300\ \mathrm{K}$ every
molecule has $\langle E \rangle = \tfrac{3}{2}kT = 6.21 \times 10^{-21}\
\mathrm{J}$; for hydrogen ($m = 3.34 \times 10^{-27}\ \mathrm{kg}$) this gives
$v \approx 1.93\ \mathrm{km/s}$, above one-sixth of Earth's
$11.2\ \mathrm{km/s}$ escape speed. Hydrogen therefore leaks away, and its
absence from the atmosphere puts a lower bound on the age of Earth of $10^8$
years.[^tl-escape]

## Equipartition of energy

The result $\langle E \rangle = \tfrac{3}{2}kT$ splits into $\tfrac{1}{2}kT$ for
each of the three translational velocity components. This is one instance of a
general theorem.

> **Theorem (Equipartition).** In thermal equilibrium each degree of freedom
> contributes $\tfrac{1}{2}kT$ to the average energy per molecule, where a
> **degree of freedom** is any coordinate or velocity component that appears
> squared in the total energy.[^tl-equip]

The theorem follows because every quadratic term in the energy is a variable
integrated against the same Gaussian Boltzmann factor, and each such integral
returns $\tfrac{1}{2}kT$. A monatomic gas molecule has three translational
degrees of freedom and $\langle E \rangle = \tfrac{3}{2}kT$. A one-dimensional
harmonic oscillator has two — the coordinate $x$ (potential energy
$\tfrac{1}{2}\kappa x^2$) and the velocity $v_x$ (kinetic energy
$\tfrac{1}{2}mv_x^2$) — and averages $kT$.

## Heat capacities of gases

The molar heat capacity at constant volume, $C_V = (\partial U/\partial T)_V$
with $U$ the internal energy of a mole, is the equipartition theorem's most
direct experimental test. A rigid diatomic molecule modeled as a dumbbell can
translate along three axes and rotate about the two axes perpendicular to the
bond, for five degrees of freedom:

$$
E = \tfrac{1}{2}mv_x^2 + \tfrac{1}{2}mv_y^2 + \tfrac{1}{2}mv_z^2
  + \tfrac{1}{2}I_x\omega_x^2 + \tfrac{1}{2}I_y\omega_y^2.
$$

Rotation about the bond axis is neglected. The average energy is
$\tfrac{5}{2}kT$ per molecule, so $U = \tfrac{5}{2}RT$ per mole and
$C_V = \tfrac{5}{2}R$. That $C_V \approx \tfrac{5}{2}R$ for nitrogen and oxygen
led Clausius to infer that these gases are diatomic rotors.

$$
% caption: A rigid diatomic molecule has three translational degrees of freedom
% and two rotational ones; rotation about the bond axis carries negligible
% moment of inertia and does not count.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
% two atoms and bond
\fill[acc!18, draw=acc, thick] (-1.4,0) circle (10pt);
\fill[acc!18, draw=acc, thick] (1.4,0) circle (10pt);
\draw[black, thick] (-1.05,0) -- (1.05,0);
\node[anchor=north, font=\scriptsize] at (0,-0.5) {bond axis};
% translation arrow (upper left, clear of rotation)
\draw[->, black, thick] (-1.7,1.2) -- (-0.5,1.2);
\node[black!70, font=\scriptsize, anchor=south] at (-1.1,1.25) {translate};
% rotation arc (upper right, separated)
\draw[->, acc, thick] (1.5,0.95) arc (-25:250:0.5 and 0.32);
\node[text=acc, font=\scriptsize, anchor=south] at (1.25,1.35) {rotate};
% count
\node[font=\scriptsize, anchor=west] at (1.95,0) {5 DOF};
\end{tikzpicture}
$$

A nonrigid molecule vibrates along the bond, adding a kinetic and a potential
term, two more degrees of freedom, predicting $C_V = \tfrac{7}{2}R$. Yet the
measured values for most diatomic gases match $\tfrac{5}{2}R$ with no
vibrational contribution, and equipartition offers no reason why some degrees of
freedom should be inactive.

| Gas | $C_V\ (\mathrm{cal/mol\!\cdot\! deg})$ | $C_V/R$ | Active DOF |
| --- | --- | --- | --- |
| Ar, He (monatomic) | 2.98 | 1.50 | 3 (translation) |
| $N_2$, $O_2$, CO, NO | $\approx 4.9$–5.0 | $\approx 2.5$ | 5 (translation + rotation) |
| $\mathrm{Cl}_2$ | 5.93 | 2.98 | between 5 and 7 |
| $\mathrm{CO}_2$, $\mathrm{SO}_2$ | $> 6.7$ | $> 3.4$ | polyatomic |

The failure sharpens with temperature. The heat capacity of $\mathrm{H}_2$
depends on $T$, contradicting equipartition, which predicts a constant. Below
about $60\ \mathrm{K}$ hydrogen behaves as a monatomic gas with
$C_V = \tfrac{3}{2}R$; between $250$ and $700\ \mathrm{K}$ it is a rigid rotor
with $\tfrac{5}{2}R$; only near dissociation does it approach $\tfrac{7}{2}R$.
Degrees of freedom switch on one at a time as $T$ rises, in a staircase that
classical mechanics cannot explain.

$$
% caption: Molar heat capacity of hydrogen against temperature. Rotational and
% vibrational degrees of freedom activate in stages, giving plateaus at
% three-halves, five-halves, and seven-halves R rather than a constant value.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,0) -- (7.0,0) node[right] {T (log scale)};
\draw[->, black] (0,0) -- (0,4.1) node[above] {$\frac{C_V}{R}$};
\foreach \y/\lab in {1.5/{$\frac{3}{2}$}, 2.5/{$\frac{5}{2}$}, 3.5/{$\frac{7}{2}$}} {
  \draw[black, dashed] (0,\y) -- (6.6,\y);
  \node[anchor=east, font=\scriptsize] at (-0.05,\y) {\lab};
}
\draw[acc, thick, smooth] plot coordinates
  {(0.2,1.5)(1.4,1.5)(1.9,1.6)(2.3,1.95)(2.8,2.35)(3.3,2.5)(4.4,2.5)(4.9,2.62)(5.4,3.0)(5.9,3.35)(6.4,3.5)};
\node[text=acc, font=\scriptsize, anchor=south] at (0.9,1.5) {translation};
\node[text=acc, font=\scriptsize, anchor=south] at (3.6,2.5) {rotation};
\node[text=acc, font=\scriptsize, anchor=south] at (5.9,3.5) {vibration};
\end{tikzpicture}
$$

## Heat capacities of solids

For solids, equipartition gives the Dulong-Petit law. Modeling each atom as a
three-dimensional oscillator bound by springs, the vibrational energy has six
quadratic terms — three kinetic, three potential,

$$
E = \tfrac{1}{2}mv_x^2 + \tfrac{1}{2}mv_y^2 + \tfrac{1}{2}mv_z^2
  + \tfrac{1}{2}\kappa_1 x^2 + \tfrac{1}{2}\kappa_2 y^2 + \tfrac{1}{2}\kappa_3 z^2,
$$

so $\langle E \rangle = 6 \times \tfrac{1}{2}kT = 3kT$ and $U = 3RT$ per mole.

> **Definition (Dulong-Petit law).** The molar heat capacity of a solid is
> $$
> C_V = 3R \approx 6\ \mathrm{cal/K\!\cdot\! mol},
> $$
> independent of the material and of temperature.[^tl-dulong]

At high temperature every solid obeys it. But below a material-dependent
critical temperature $C_V$ drops toward zero as $T \to 0$, the critical
temperature being lower for soft solids such as lead and higher for hard ones
such as diamond. Equipartition, temperature-independent by construction, has no
account of this fall.

A second failure is quantitative. The classical free-electron picture of a
metal treats roughly one conduction electron per atom as a gas, which by
equipartition should add $\tfrac{3}{2}R$ to the heat capacity. Metals show no
such extra contribution: their heat capacities scarcely exceed those of
insulators. Both failures — the low-temperature collapse of $C_V$ and the
missing electron contribution — trace to the same source. Classical mechanics
is the wrong mechanics for atoms. The energy of an oscillator or a rotor is
quantized, degrees of freedom whose quantum spacing exceeds $kT$ are frozen out,
and the electrons obey an exclusion principle that keeps almost all of them from
sharing thermal energy. The repair requires the [quantum
distributions](/statistical-mechanics/quantum-statistics/quantum-statistics-bose-einstein-and-fermi-dirac), and the
[fermion gas](/statistical-mechanics/bose-systems/bose-einstein-condensation-and-the-fermion-gas) explains the
metals. The search for an understanding of specific heats was, historically, one
of the roads into the quantum theory itself.

[^tl-intro]: **Tipler & Llewellyn**, _Modern Physics_, Ch. 8 introduction — the statistical approach to systems of order $10^{22}$ particles, predicting bulk properties from probability rather than from individual trajectories.
[^tl-boltz]: **Tipler & Llewellyn**, _Modern Physics_, §8-1 (Boltzmann Distribution), Eqs. 8-1 and 8-2 — the distribution $f_B(E) = A e^{-E/kT}$, the Boltzmann factor and constant, and the statistical weight $g(E)$.
[^tl-atmos]: **Tipler & Llewellyn**, _Modern Physics_, §8-1, Example 8-1 (The Law of Atmospheres) and Example 8-2 (H atoms in the first excited state), Eq. 8-3.
[^tl-entropy]: **Tipler & Llewellyn**, _Modern Physics_, §8-1 "Temperature and Entropy" (Eqs. 8-4a–d) — the statistical definition of temperature and entropy through the multiplicity of microstates, $S = k \ln W$.
[^tl-escape]: **Tipler & Llewellyn**, _Modern Physics_, §8-1, Example 8-4 (Escape of $\mathrm{H}_2$ from Earth's Atmosphere), Eqs. 8-13 and 8-14.
[^tl-equip]: **Tipler & Llewellyn**, _Modern Physics_, §8-1 (Heat Capacities of Gases and Solids) and "A Derivation of the Equipartition Theorem" — each squared coordinate or velocity term contributes $\tfrac{1}{2}kT$.
[^tl-dulong]: **Tipler & Llewellyn**, _Modern Physics_, §8-1 ($C_V$ for Solids) — the Dulong-Petit result $C_V = 3R$ from six quadratic terms per atom, and its low-temperature failure.