Boltzmann Distribution - Physical Chemistry 2

SECTION 1

Historical Context & Motivation

By the mid-nineteenth century, thermodynamics had matured into a powerful macroscopic framework—yet it offered no microscopic explanation for why gases behave the way they do. The pressure, temperature, and entropy of a system were empirically well-characterized, but the link between these bulk observables and the motion of individual molecules remained elusive. Ludwig Boltzmann set out to bridge that gap by asking a deceptively simple question: if a gas contains an enormous number of molecules sharing a fixed total energy, how is that energy distributed among them at equilibrium? His answer—grounded in probability theory and combinatorics—gave birth to statistical mechanics and produced one of the most consequential equations in all of physical science.

1860 • Maxwell's Speed Distribution

James Clerk Maxwell derives the distribution of molecular speeds in an ideal gas, providing the first statistical description of a thermodynamic system and setting the stage for Boltzmann's more general treatment.

1868 • Boltzmann's Generalization

Ludwig Boltzmann extends Maxwell's work to arbitrary energy degrees of freedom using the concept of the ergodic hypothesis, deriving the exponential energy distribution that now bears his name.

1877 • S = k ln W

Boltzmann publishes his celebrated entropy formula, connecting the macroscopic entropy S to the logarithm of the number of microstates W. This equation, later inscribed on his tombstone, provides the statistical foundation for the second law of thermodynamics.

1902 • Gibbs' Ensemble Theory

Josiah Willard Gibbs formalizes ensemble theory in his monograph, placing the Boltzmann distribution within the broader canonical ensemble framework and extending the formalism to interacting systems.

1920s • Quantum Corrections

The advent of quantum mechanics introduces Fermi-Dirac and Bose-Einstein statistics for indistinguishable particles, refining the classical Boltzmann distribution into a high-temperature limiting case of more general quantum statistics.

The central question the Boltzmann distribution answers is this: given a large collection of particles in thermal equilibrium at temperature T, what fraction of those particles will be found in any particular energy state? Understanding the answer is prerequisite to virtually every calculation in statistical thermodynamics—from predicting heat capacities and equilibrium constants to interpreting spectroscopic intensities and reaction rates.

SECTION 2

Core Principles & Definitions

The Boltzmann distribution rests on several foundational ideas drawn from probability theory, classical mechanics, and quantum mechanics. Before diving into the mathematics, it is essential to establish a clear conceptual vocabulary. The following principles form the intellectual scaffolding upon which the entire distribution is built, and each one addresses a specific aspect of how energy and probability interrelate in a system of many particles.

Microstates & Macrostates

A microstate is a complete specification of every particle's quantum numbers (or classical coordinates and momenta). A macrostate is defined by bulk properties (E, V, N). Many microstates map to the same macrostate.

Equal A Priori Probabilities

The fundamental postulate of statistical mechanics: every accessible microstate of an isolated system is equally probable at equilibrium. This postulate replaces dynamic equations of motion with a probabilistic framework.

The Partition Function (q or Z)

The partition function is the normalization constant that ensures all occupation probabilities sum to one. It encodes the complete thermodynamic information of the system and serves as a generating function for all ensemble averages.

Boltzmann Factor

The Boltzmann factor exp(−εᵢ / k_BT) gives the relative probability of occupying state i with energy εᵢ. Higher-energy states are exponentially less populated.

Temperature as a Distribution Parameter

Temperature T controls the width of the energy distribution. As T → 0, all particles collapse to the ground state; as T → ∞, all states become equally populated. Temperature is the Lagrange multiplier associated with the energy constraint.

✦

KEY TAKEAWAY

Think of the Boltzmann distribution as an energy "budget allocation" for molecules. Imagine a university distributing a fixed pool of scholarship money among students: most students receive modest awards (low energy), a smaller number receive generous scholarships (moderate energy), and only a very few receive the top prize (high energy). The temperature is analogous to the total budget—when it is large, the distribution flattens and more students can receive larger awards. The partition function is simply the accounting ledger that ensures every dollar is tracked.

SECTION 3

Visual Explanation

The Boltzmann distribution is best understood visually by examining how the population of energy levels changes with temperature. The following diagram shows a set of equally spaced energy levels and the relative number of particles occupying each level at two different temperatures. At low temperature, the population is strongly skewed toward the ground state. At high temperature, the distribution broadens and higher-energy states become appreciably populated.

Comparison of Boltzmann populations across six equally spaced energy levels at low temperature T₁ (left, blue) and high temperature T₂ = 3T₁ (right, red). Bar lengths represent the population ratio n_i/n₀ = exp(−iε/k_BT). Notice how the high-temperature distribution is substantially flatter, with appreciable population even at ε₅.

Several features of the diagram deserve careful attention. First, the ground state (ε₀) is always the most populated level regardless of temperature—this is a direct consequence of the exponential decay of the Boltzmann factor. Second, the ratio between successive populations is constant for equally spaced levels: n_i+1/n_i = exp(−Δε/k_BT), a quantity that approaches unity as T increases and approaches zero as T decreases. Third, the characteristic temperature Θ = ε/k_B provides a natural scale: when T ≪ Θ, only the ground state matters; when T ≫ Θ, many states are thermally accessible.

SECTION 4

Mathematical Framework

The Boltzmann distribution can be derived rigorously by maximizing the number of microstates W subject to constraints on total particle number N and total energy E, using the method of Lagrange multipliers. Equivalently, it emerges naturally from the canonical ensemble by considering a small subsystem in thermal contact with a large heat bath. Both approaches yield the same result, confirming the internal consistency of statistical mechanics. We present the key equations and their physical interpretation below.

BOLTZMANN DISTRIBUTION LAW

pᵢ = (gᵢ × exp(−εᵢ / k_B T)) / q

p_i = probability of occupying state i; g_i = degeneracy of level i; ε_i = energy of state i; k_B = Boltzmann constant (1.381 × 10⁻²³ J K⁻¹); T = absolute temperature (K); q = molecular partition function.

MOLECULAR PARTITION FUNCTION

q = Σᵢ gᵢ × exp(−εᵢ / k_B T)

The sum runs over all energy levels i. For non-degenerate states (g_i = 1), the partition function reduces to q = Σ exp(−ε_i / k_BT). The partition function counts the effective number of thermally accessible states.

POPULATION RATIO

nⱼ / nᵢ = (gⱼ / gᵢ) × exp(−(εⱼ − εᵢ) / k_B T)

This form is particularly useful in spectroscopy. It gives the ratio of populations in two levels directly, without needing the partition function. When j is the upper level and i is the lower, the ratio is always ≤ g_j/g_i at finite temperature.

AVERAGE ENERGY

⟨E⟩ = −(∂ ln q / ∂β)_V where β = 1 / k_B T

This expression demonstrates the power of the partition function as a generating function. By taking appropriate derivatives of ln q with respect to β (or T), one can extract the average energy, heat capacity, entropy, and other thermodynamic quantities without enumerating individual microstate probabilities.

ℹ

Derivation Sketch

Start from the number of ways to distribute N distinguishable particles among energy levels: W = N! / Πᵢ nᵢ!. Maximize ln W using Stirling's approximation (ln n! ≈ n ln n − n) subject to constraints Σ nᵢ = N and Σ nᵢεᵢ = E. Introduce Lagrange multipliers α and β. Setting ∂(ln W − αΣnᵢ − βΣnᵢεᵢ)/∂nᵢ = 0 yields nᵢ = exp(−α) × exp(−βεᵢ). Identifying β = 1/k_BT from thermodynamic consistency and using the normalization constraint to fix α recovers the Boltzmann distribution.

SECTION 5

Partition Functions by Degree of Freedom

For molecular systems, the total partition function factorizes into contributions from different degrees of freedom, provided those motions are approximately separable. This decomposition is one of the most powerful practical consequences of the Boltzmann framework, because it allows us to treat translation, rotation, vibration, and electronic excitation independently. The following diagram and table summarize the key partition functions encountered in physical chemistry.

The total molecular partition function factorizes into translational, rotational, vibrational, and electronic contributions (top). Each component has a characteristic temperature scale (Θ_rot, Θ_vib) that determines when that mode is thermally excited. The lower box shows how macroscopic thermodynamic quantities are extracted from q.

Mode	Characteristic Temperature	Typical Range	High-T Contribution to $C_V$
Translation	Θ_trans ≈ 10⁻¹⁵ K	Always classical	³⁄₂ k_B per molecule
Rotation (linear)	Θ_rot ≈ 2–85 K	Classical above ~100 K	k_B per molecule
Vibration	Θ_vib ≈ 300–6000 K	Often frozen at 300 K	k_B per mode per molecule
Electronic	Θ_elec ≈ 10⁴–10⁵ K	Usually only ground state	≈ 0 (except for low-lying states)

Characteristic temperatures and classical contributions for molecular degrees of freedom

The hierarchy of characteristic temperatures explains why the heat capacity of a diatomic gas increases stepwise with temperature. At room temperature, translational and rotational modes are fully excited (contributing ⁵⁄₂ k_B per molecule to C_V), vibrational modes are largely frozen out, and electronic contributions are negligible. This is the equipartition theorem in action—or rather, the Boltzmann distribution revealing exactly when and why equipartition fails.

SECTION 6

Worked Example

Consider the vibrational mode of carbon monoxide (CO), which has a fundamental vibrational frequency ν̃ = 2170 cm⁻¹. We wish to calculate the fraction of CO molecules in the first excited vibrational state (v = 1) relative to the ground state (v = 0) at T = 1000 K, and the vibrational partition function at that temperature.

Vibrational Population and Partition Function of CO at 1000 K

1
Step 1 — Convert Vibrational Frequency to Energy

The vibrational energy spacing is Δε = hcν̃, where h = 6.626 × 10⁻³⁴ J·s and c = 2.998 × 10¹⁰ cm/s. Substituting: Δε = (6.626 × 10⁻³⁴)(2.998 × 10¹⁰)(2170) = 4.312 × 10⁻²⁰ J. Alternatively, compute the characteristic vibrational temperature: Θ_vib = hcν̃ / k_B = 4.312 × 10⁻²⁰ / 1.381 × 10⁻²³ = 3123 K.

Θ_vib = 3123 K
2
Step 2 — Compute the Boltzmann Factor

The population ratio n₁/n₀ for a harmonic oscillator (all g_i = 1) is exp(−Θ_vib/T) = exp(−3123/1000) = exp(−3.123) = 0.04405.

n₁/n₀ = 0.0441
3
Step 3 — Evaluate the Vibrational Partition Function

For the quantum harmonic oscillator, the vibrational partition function (measured from the zero-point level) is the geometric series q_vib = 1 / (1 − exp(−Θ_vib/T)). Substituting: q_vib = 1 / (1 − 0.04405) = 1 / 0.9560 = 1.046.

q_vib = 1.046
4
Step 4 — Determine Absolute Population Fraction

The fraction of molecules in v = 0 is p₀ = 1/q_vib = 1/1.046 = 0.9560 (95.6%). The fraction in v = 1 is p₁ = 0.04405/1.046 = 0.04211 (4.2%). Only about 4% of CO molecules occupy the first excited vibrational state at 1000 K.

p₀ = 95.6%, p₁ = 4.2%
5
Step 5 — Interpret the Result

A partition function of q_vib ≈ 1.05 tells us that effectively only one vibrational state is significantly populated. This makes physical sense: at 1000 K, the thermal energy k_BT = 1.381 × 10⁻²⁰ J is roughly three times smaller than the vibrational quantum Δε = 4.3 × 10⁻²⁰ J, so the vibration is only partially 'thawed.' At T = 6000 K, the partition function would rise to approximately 2.9, indicating nearly three accessible states.

SECTION 7

Strengths, Limitations, and Scope

The Boltzmann distribution is extraordinarily powerful, but it is not universal. Understanding its domain of validity is essential to applying it correctly and recognizing when more sophisticated statistical frameworks are needed. The following table contrasts the strengths and limitations of the classical Boltzmann approach.

Strengths	Limitations
Exact for distinguishable particles (e.g., localized spins, nuclear states in NMR) and as the high-temperature limit of quantum statistics.	Fails for indistinguishable particles at low T or high density where quantum exchange effects become important (need Fermi-Dirac or Bose-Einstein statistics).
Yields closed-form partition functions for many model systems (harmonic oscillator, rigid rotor, particle in a box).	Assumes thermal equilibrium—inapplicable to transient, driven, or far-from-equilibrium systems without modification.
Connects microscopic quantum states directly to measurable thermodynamic quantities (U, S, $C_V$ , K_eq) through derivatives of ln q.	Requires knowledge of the energy level spectrum εᵢ; for complex many-body systems, this may not be analytically tractable.
Provides physical intuition: temperature controls the 'softness' of energy distribution; the partition function counts effective states.	Neglects inter-particle interactions in the ideal-gas formulation; corrections (virial expansion, configuration integrals) are needed for real gases and condensed phases.

Strengths and limitations of the Boltzmann distribution framework

✦

KEY TAKEAWAY

The Boltzmann distribution occupies a central position in statistical mechanics analogous to Newton's second law in classical mechanics: it is the default starting point for virtually every problem, and its limitations signal the need for more advanced frameworks (quantum statistics, nonequilibrium theory, or molecular simulation). In practice, for molecular gases above roughly 50 K, the Boltzmann distribution is an excellent approximation for translational, rotational, and often vibrational degrees of freedom.

SECTION 8

Connection to Quantum Statistics

The Boltzmann distribution assumes that particles are distinguishable and that no single quantum state has a macroscopic occupation number. When these assumptions break down—as they do for photons in a cavity, electrons in a metal, or atoms in a Bose-Einstein condensate—one must turn to the full quantum statistics developed by Fermi, Dirac, Bose, and Einstein. The table below highlights how the Boltzmann distribution relates to its quantum generalizations.

Feature	Boltzmann	Fermi-Dirac	Bose-Einstein
Particle type	Distinguishable (or dilute limit)	Fermions (half-integer spin)	Bosons (integer spin)
Mean occupation number ⟨n⟩	exp(−(ε − μ)/k_BT)	1 / [exp((ε − μ)/k_BT) + 1]	1 / [exp((ε − μ)/k_BT) − 1]
Max occupation per state	Unlimited (but assumed ≪ 1)	0 or 1 (Pauli exclusion)	Unlimited (bunching allowed)
Classical limit	Always valid	When exp((ε−μ)/k_BT) ≫ 1	When exp((ε−μ)/k_BT) ≫ 1
Example system	Molecular gas at moderate T	Electrons in metals, white dwarfs	Photons (Planck radiation), ⁴He superfluid

Comparison of classical (Boltzmann) and quantum (Fermi-Dirac, Bose-Einstein) distribution functions

The essential insight is that the Boltzmann distribution emerges as the dilute limit of both Fermi-Dirac and Bose-Einstein statistics. When the number of available quantum states vastly exceeds the number of particles—as is typically the case for molecules in a gas at ordinary temperatures—the probability of any two particles competing for the same state is negligible, and the ±1 in the denominator of the quantum distributions becomes irrelevant. This condition is quantified by requiring that the thermal de Broglie wavelength Λ be much smaller than the mean inter-particle spacing (V/N)¹ᐟ³, i.e., nΛ³ ≪ 1, where n = N/V is the number density.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain why the partition function q can be interpreted as the "effective number of thermally accessible states." Under what temperature conditions does q approach 1, and under what conditions does it grow large? What physical insight does each limit provide?

Show answer

The partition function q = Σ gᵢ exp(−εᵢ/k_BT) sums the Boltzmann weights of all states. When T → 0, every exponential except the ground-state term vanishes, so q → g₀ (typically 1), meaning only the ground state is accessible. When T → ∞, all exponentials approach 1, so q → Σ gᵢ, the total number of states, indicating every state is equally accessible. At intermediate temperatures, q counts how many states have significant Boltzmann weight. A q value of, say, 3.5 means the system's energy is effectively spread over about 3–4 states.

Question 2: BASIC CALCULATION

A two-level system has a ground state at ε₀ = 0 and an excited state at ε₁ = 4.0 × 10⁻²¹ J, both non-degenerate. Calculate the probability of finding the system in the excited state at T = 300 K. (Use k_B = 1.381 × 10⁻²³ J/K.)

Show answer

First compute the exponent: ε₁/(k_BT) = 4.0 × 10⁻²¹ / (1.381 × 10⁻²³ × 300) = 4.0 × 10⁻²¹ / 4.143 × 10⁻²¹ = 0.9655. The Boltzmann factor is exp(−0.9655) = 0.3808. The partition function is q = 1 + 0.3808 = 1.3808. Therefore, p₁ = 0.3808 / 1.3808 = 0.2758, or about 27.6%. This result indicates that at 300 K, the thermal energy k_BT is comparable to the energy gap, so the excited state is appreciably populated.

Question 3: INTERMEDIATE

The rotational constant of HCl is B̃ = 10.59 cm⁻¹. The rotational energy levels are ε_J = hcB̃ J(J+1) with degeneracy g_J = 2J+1. Determine which rotational level J has the maximum population at T = 300 K. (Θ_rot = hcB̃/k_B = 15.24 K.)

Show answer

The population of level J is proportional to f(J) = (2J+1) exp(−Θ_rot J(J+1)/T). To find the maximum, treat J as a continuous variable and differentiate: df/dJ = 2 exp(−Θ_rot J(J+1)/T) − (2J+1)² (Θ_rot/T) exp(−Θ_rot J(J+1)/T) = 0. This gives J_max = √(T/(2Θ_rot)) − ½ = √(300/(2 × 15.24)) − 0.5 = √(9.843) − 0.5 = 3.137 − 0.5 = 2.64. Since J must be an integer, the most populated level is J = 3. This illustrates how the interplay between increasing degeneracy (2J+1) and decreasing Boltzmann factor creates a population maximum at J > 0.

Question 4: APPLIED

In a flame at T = 2500 K, sodium atoms can be in the ground ²S_1/2 state (g₀ = 2, ε₀ = 0) or the first excited ²P_1/2,3/2 states (combined g₁ = 6, ε₁ corresponds to the 589 nm D-line). Calculate the fraction of sodium atoms in the excited state and discuss the implications for atomic emission spectroscopy.

Show answer

The energy gap is ε₁ = hc/λ = (6.626 × 10⁻³⁴)(2.998 × 10⁸)/(589 × 10⁻⁹) = 3.373 × 10⁻¹⁹ J. The dimensionless ratio ε₁/(k_BT) = 3.373 × 10⁻¹⁹ / (1.381 × 10⁻²³ × 2500) = 9.77. The Boltzmann factor is exp(−9.77) = 5.67 × 10⁻⁵. Including degeneracies: g₁ exp(−ε₁/k_BT) = 6 × 5.67 × 10⁻⁵ = 3.40 × 10⁻⁴. The partition function q ≈ 2 + 3.40 × 10⁻⁴ ≈ 2.00034. The excited-state fraction is p₁ = 3.40 × 10⁻⁴ / 2.00034 ≈ 1.7 × 10⁻⁴, or about 0.017%. Although this seems tiny, it is sufficient to produce the intense yellow emission of sodium in a flame. This is because even a small absolute population in a radiatively efficient state produces a detectable signal, illustrating how the Boltzmann distribution underpins quantitative flame photometry.

Question 5: CRITICAL THINKING

Consider a system of N independent quantum harmonic oscillators (frequency ν) in thermal equilibrium. Derive an expression for the average energy per oscillator ⟨ε⟩ as a function of T using the vibrational partition function. Show that your result reduces to k_BT in the high-temperature limit and to ½hν (zero-point energy contribution only) in the low-temperature limit. Discuss the physical significance of these limits in the context of the historical ultraviolet catastrophe.

Show answer

Starting from q_vib = exp(−Θ_v/(2T)) / [1 − exp(−Θ_v/T)] (including zero-point energy), where Θ_v = hν/k_B, use ⟨ε⟩ = −∂(ln q)/∂β to obtain ⟨ε⟩ = ½hν + hν/[exp(hν/k_BT) − 1]. In the high-T limit (k_BT ≫ hν), expand the exponential to first order: exp(x) ≈ 1 + x, giving ⟨ε⟩ ≈ ½hν + k_BT − ½hν = k_BT, recovering the classical equipartition result. In the low-T limit (k_BT ≪ hν), the exponential dominates and the second term vanishes, leaving ⟨ε⟩ → ½hν. The high-T limit corresponds to the Rayleigh-Jeans regime where classical equipartition assigns k_BT to every mode, leading to the ultraviolet catastrophe (divergent total energy). The quantized Boltzmann distribution automatically 'freezes out' high-frequency modes at finite T, resolving the catastrophe and providing one of the earliest confirmations that energy is quantized.

SUMMARY

Summary

The Boltzmann distribution is the cornerstone of statistical thermodynamics, providing the probability p_i = g_iexp(−ε_i/k_BT)/q that a molecule occupies energy state i at thermal equilibrium. The Boltzmann factor exp(−ε/k_BT) encodes the exponential suppression of high-energy states, while the partition function q = Σ g_iexp(−ε_i/k_BT) normalizes probabilities and serves as a thermodynamic generating function. Temperature controls the width of the distribution: low T concentrates population in the ground state, while high T spreads it across many levels.

For molecular systems, the partition function factorizes into translational, rotational, vibrational, and electronic contributions, each governed by a characteristic temperature that determines whether that mode is classically excited or quantum-mechanically frozen. The distribution is exact for distinguishable particles and serves as the high-temperature classical limit of both Fermi-Dirac and Bose-Einstein quantum statistics. Mastery of the Boltzmann distribution provides the essential bridge from microscopic quantum mechanics to macroscopic thermodynamics.