Chapter 8 • Theory & Derivations

Unit 8: Molecular Partition Functions & The Thermodynamic Bridge

Rigorous quantum-to-classical factorization of the molecular partition function: translational density of states and thermal de Broglie wavelength, rotational partition functions of linear and non-linear polyatomic rotors, symmetry numbers and nuclear spin statistics (ortho/para hydrogen), vibrational partition functions and zero-point energy, electronic degeneracies, temperature-dependent heat capacity C_V(T) equipartition transitions, and the Sackur-Tetrode formula for absolute translational entropy.

§8.1 Factorization of the Molecular Partition Function

For an ideal gas molecule, the total Hamiltonian can be approximated to high accuracy as a sum of uncoupled independent degrees of freedom:

\hat{H}_{\text{mol}} \approx \hat{H}_{\text{trans}} + \hat{H}_{\text{rot}} + \hat{H}_{\text{vib}} + \hat{H}_{\text{elec}} + \hat{H}_{\text{nuc}}

Because the Hamiltonian operators commute, the total energy of a single molecule is additive:

\varepsilon = \varepsilon_{\text{trans}} + \varepsilon_{\text{rot}} + \varepsilon_{\text{vib}} + \varepsilon_{\text{elec}} + \varepsilon_{\text{nuc}}

Factorization of the Partition Function

The single-molecule partition function is:

q = \sum_i e^{-\beta \varepsilon_i} = \sum_{\text{trans}} \sum_{\text{rot}} \sum_{\text{vib}} \sum_{\text{elec}} \sum_{\text{nuc}} e^{-\beta (\varepsilon_{\text{trans}} + \varepsilon_{\text{rot}} + \varepsilon_{\text{vib}} + \varepsilon_{\text{elec}} + \varepsilon_{\text{nuc}})}

Because the exponential of a sum is the product of exponentials, the multiple sum factors cleanly:

q = q_{\text{trans}} \times q_{\text{rot}} \times q_{\text{vib}} \times q_{\text{elec}} \times q_{\text{nuc}}

Taking natural logarithms converts this into an additive sum of terms:

\ln q = \ln q_{\text{trans}} + \ln q_{\text{rot}} + \ln q_{\text{vib}} + \ln q_{\text{elec}} + \ln q_{\text{nuc}}

Consequently, every thermodynamic property decomposes into additive components:

U = U_{\text{trans}} + U_{\text{rot}} + U_{\text{vib}} + U_{\text{elec}}

C_V = C_{V, \text{trans}} + C_{V, \text{rot}} + C_{V, \text{vib}} + C_{V, \text{elec}}

S = S_{\text{trans}} + S_{\text{rot}} + S_{\text{vib}} + S_{\text{elec}}

This factorization forms the computational foundation for predicting all macroscopic gas-phase thermochemical properties directly from molecular spectroscopic data.

§8.2 Translational Partition Function & The Thermal de Broglie Wavelength

Consider a molecule of mass \(m\) moving freely inside a three-dimensional rectangular box of dimensions \(a \times b \times c = V\). The quantized translational energy levels are:

\varepsilon_{n_x, n_y, n_z} = \frac{h^2}{8m} \left( \frac{n_x^2}{a^2} + \frac{n_y^2}{b^2} + \frac{n_z^2}{c^2} \right)

Because typical molecular masses and macroscopic volumes produce extraordinarily tiny energy level spacings (\(\Delta \varepsilon \sim 10^{-38}\text{ J} \ll k_B T \sim 4 \times 10^{-21}\text{ J}\)), the quantum sum can be replaced by an integral:

q_{\text{trans}, x} = \sum_{n_x=1}^\infty e^{-\frac{h^2 n_x^2}{8 m a^2 k_B T}} \approx \int_0^\infty e^{-\frac{h^2 n^2}{8 m a^2 k_B T}} dn

Using the standard Gaussian integral \(\int_0^\infty e^{-\alpha x^2} dx = \frac{1}{2} \sqrt{\frac{\pi}{\alpha}}\):

q_{\text{trans}, x} = \frac{1}{2} \sqrt{\frac{\pi \times 8 m a^2 k_B T}{h^2}} = \frac{a}{h} \sqrt{2\pi m k_B T}

For 3D translation:

q_{\text{trans}} = q_x q_y q_z = \frac{a b c}{h^3} (2\pi m k_B T)^{3/2} = \frac{V}{h^3} (2\pi m k_B T)^{3/2}

The Thermal de Broglie Wavelength \(\Lambda\)

We define the thermal de Broglie wavelength:

\Lambda = \frac{h}{\sqrt{2\pi m k_B T}} = \frac{h}{p_{\text{thermal}}}

This represents the quantum spatial coherence length of a thermal particle at temperature \(T\). In terms of \(\Lambda\), the translational partition function assumes the remarkably simple form:

q_{\text{trans}} = \frac{V}{\Lambda^3}

Physical Meaning of \(\Lambda\):

\(V / \Lambda^3\) is the ratio of macroscopic container volume to the effective quantum volume \(\Lambda^3\) of a particle.
When \(\Lambda \ll d\) (where \(d = (V/N)^{1/3}\) is the average interparticle distance), quantum wavepackets do not overlap, and the gas behaves strictly as a classical Maxwell-Boltzmann system.
When \(\Lambda \sim d\), wavepacket overlap becomes significant, requiring quantum statistics (Bose-Einstein or Fermi-Dirac).

§8.3 Rotational Partition Function, Symmetry Numbers & Spin Statistics

Linear Rotors (Diatomic & Linear Polyatomic Molecules)

The rigid rotor quantum energy levels with moment of inertia \(I\) are:

\varepsilon_J = \frac{\hbar^2}{2 I} J(J+1) = h c B J(J+1), \quad g_J = 2J + 1 \quad (J = 0, 1, 2, \dots)

where the rotational constant is \(B = \frac{\hbar}{4\pi c I}\). Defining the characteristic rotational temperature \(\theta_{\text{rot}} = \frac{h c B}{k_B} = \frac{\hbar^2}{2 I k_B}\):

q_{\text{rot}} = \sum_{J=0}^\infty (2J + 1) \exp\left( -J(J+1) \frac{\theta_{\text{rot}}}{T} \right)

At temperatures \(T \gg \theta_{\text{rot}}\) (which is true for almost all molecules at room temperature except \(\text{H}_2\)), the sum can be approximated by an integral:

q_{\text{rot}} \approx \int_0^\infty (2J + 1) e^{-\frac{\theta_{\text{rot}}}{T} (J^2 + J)} dJ = \left[ -\frac{T}{\theta_{\text{rot}}} e^{-\frac{\theta_{\text{rot}}}{T} (J^2 + J)} \right]_0^\infty = \frac{T}{\theta_{\text{rot}}} = \frac{8\pi^2 I k_B T}{h^2}

The Rotational Symmetry Number \(\sigma\)

For homonuclear or symmetrical molecules, rotating the molecule by \(180^\circ\) interchanges identical nuclei. To avoid double-counting identical indistinguishable configurations in phase space, we divide by the rotational symmetry number \(\sigma\):

q_{\text{rot}} = \frac{T}{\sigma \theta_{\text{rot}}} = \frac{8\pi^2 I k_B T}{\sigma h^2}

Values of \(\sigma\):

Heteronuclear diatomics (\(\text{HCl}, \text{CO}\)): \(\sigma = 1\)
Homonuclear diatomics (\(\text{N}_2, \text{O}_2, \text{H}_2\)): \(\sigma = 2\)
Water (\(\text{H}_2\text{O}\), \(C_{2v}\)): \(\sigma = 2\)
Ammonia (\(\text{NH}_3\), \(C_{3v}\)): \(\sigma = 3\)
Methane (\(\text{CH}_4\), \(T_d\)): \(\sigma = 12\)
Benzene (\(\text{C}_6\text{H}_6\), \(D_{6h}\)): \(\sigma = 12\)

Non-Linear Polyatomic Rotors

For asymmetric and spherical tops with three principal moments of inertia \(I_A, I_B, I_C\):

q_{\text{rot}} = \frac{\sqrt{\pi}}{\sigma} \left( \frac{T^3}{\theta_A \theta_B \theta_C} \right)^{1/2} = \frac{\sqrt{\pi}}{\sigma} \left( \frac{8\pi^2 k_B T}{h^2} \right)^{3/2} (I_A I_B I_C)^{1/2}

Nuclear Spin Statistics (Ortho vs Para Hydrogen)

In molecular hydrogen \(\text{H}_2\), the protons are fermions (\(I_{\text{nuc}} = 1/2\)). The total wavefunction must be antisymmetric with respect to proton exchange.

Para-Hydrogen (Singlet spin, \(I_{\text{tot}} = 0\), antisymmetric spin): Must pair with symmetric rotational states (\(J = 0, 2, 4, \dots\)). Nuclear spin statistical weight: \(g_{\text{nuc}} = 1\).
Ortho-Hydrogen (Triplet spin, \(I_{\text{tot}} = 1\), symmetric spin): Must pair with antisymmetric rotational states (\(J = 1, 3, 5, \dots\)). Nuclear spin statistical weight: \(g_{\text{nuc}} = 3\).

At high temperature, the equilibrium ratio is ortho:para = 3:1. At low temperature (\(T \rightarrow 0\)), hydrogen relaxes exclusively into the \(J = 0\) para-state.

§8.4 Vibrational Partition Function & Zero-Point Energy

A molecule with \(N_{\text{at}}\) atoms possesses \(3 N_{\text{at}} - 5\) normal vibrational modes if linear, and \(3 N_{\text{at}} - 6\) normal modes if non-linear. In the harmonic oscillator approximation, each normal mode \(i\) with fundamental frequency \(\nu_i\) is independent:

\varepsilon_{v_i} = \left( v_i + \frac{1}{2} \right) h \nu_i \quad (v_i = 0, 1, 2, \dots)

Defining the characteristic vibrational temperature \(\theta_{\text{vib}, i} = \frac{h \nu_i}{k_B} = \frac{h c \tilde{\nu}_i}{k_B}\):

Zero-Point Energy Conventions

1. Measured relative to the bottom of the potential well:

q_{\text{vib}, i} = \sum_{v=0}^\infty e^{-\beta (v + 1/2) h\nu_i} = e^{-\theta_{\text{vib}, i} / 2T} \sum_{v=0}^\infty e^{-v \theta_{\text{vib}, i} / T} = \frac{e^{-\theta_{\text{vib}, i} / 2T}}{1 - e^{-\theta_{\text{vib}, i} / T}}

2. Measured relative to the \(v = 0\) vibrational ground state:

Setting \(\varepsilon_0 = 0\) by factoring out the zero-point energy:

q_{\text{vib}, i}' = \sum_{v=0}^\infty e^{-v \theta_{\text{vib}, i} / T} = \frac{1}{1 - e^{-\theta_{\text{vib}, i} / T}}

For a polyatomic molecule with all normal modes:

q_{\text{vib}} = \prod_{i=1}^{3N_{\text{at}} - 6} \frac{1}{1 - e^{-\theta_{\text{vib}, i} / T}}

High and Low Temperature Limits

Because vibrational frequencies are typically high (\(\tilde{\nu} \sim 500 - 3500\text{ cm}^{-1} \implies \theta_{\text{vib}} \sim 700 - 5000\text{ K}\)):

At room temperature (\(T \approx 300\text{ K} \ll \theta_{\text{vib}}\)): \(e^{-\theta_{\text{vib}}/T} \ll 1\), so \(q_{\text{vib}} \approx 1\). Molecules are frozen in their ground vibrational state (\(v = 0\)).
At high temperature (\(T \gg \theta_{\text{vib}}\)): Taylor expand \(1 - e^{-\theta_{\text{vib}}/T} \approx \frac{\theta_{\text{vib}}}{T}\):

q_{\text{vib}} \approx \frac{T}{\theta_{\text{vib}}} = \frac{k_B T}{h \nu}

recovering the classical harmonic oscillator partition function.

§8.5 Electronic and Nuclear Spin Partition Functions

Electronic Partition Function

The electronic partition function is:

q_{\text{elec}} = \sum_i g_{e, i} e^{-\beta \varepsilon_{e, i}} = g_{e, 0} + g_{e, 1} e^{-\varepsilon_{e, 1} / k_B T} + \dots

For the vast majority of stable molecules (e.g., \(\text{N}_2, \text{H}_2, \text{CO}_2, \text{CH}_4, \text{H}_2\text{O}\)), the electronic ground state is a closed-shell singlet (\(^1\Sigma\) or \(^1A_1\)) with degeneracy \(g_{e, 0} = 1\), and the first electronically excited state lies in the ultraviolet (\(\Delta \varepsilon_1 \sim 3 - 8\text{ eV} \implies \theta_{\text{elec}} \sim 35,000 - 100,000\text{ K}\)). Therefore, at all accessible chemical temperatures:

q_{\text{elec}} \approx g_{e, 0}

Important exceptions:

1. Triplet Ground States:

Molecular oxygen \(\text{O}_2\) has a \(^3\Sigma_g^-\) ground state with \(g_{e, 0} = 3\), so \(q_{\text{elec}} = 3\).

2. Open-Shell Radicals:

Nitric oxide (\(\text{NO}\)): Has a \(^2\Pi_{1/2}\) ground state (\(g_0 = 2\)) and a low-lying \(^2\Pi_{3/2}\) excited state (\(g_1 = 2\)) only \(\Delta\tilde{\nu} = 121.1\text{ cm}^{-1}\) (\(\theta_{\text{elec}} \approx 174\text{ K}\)) above the ground state:

q_{\text{elec}}(\text{NO}) = 2 + 2 e^{-174 / T}

This gives rise to an electronic contribution to the heat capacity at low temperatures.

Atomic Halogens (e.g., \(\text{F}, \text{Cl}\)): \(^2P_{3/2}\) ground state and \(^2P_{1/2}\) excited state separated by spin-orbit coupling.

Nuclear Spin Partition Function

For a molecule containing nuclei with spin quantum numbers \(I_1, I_2, \dots\), the nuclear spin degeneracy is:

g_{\text{nuc}} = \prod_k (2 I_k + 1)

Because nuclear spin energy splittings are infinitesimal (\(\sim 10^{-6}\text{ K}\) in zero magnetic field), all nuclear spin states are equally populated:

q_{\text{nuc}} = \prod_k (2 I_k + 1)

Nuclear spin partition functions cancel identically in all chemical equilibrium and reaction rate calculations.

§8.6 Statistical Heat Capacities & The Classical Equipartition Limits

The molar isochoric heat capacity is the sum of contributions from all active degrees of freedom:

C_{V, \text{molar}} = C_{V, \text{trans}} + C_{V, \text{rot}} + C_{V, \text{vib}} + C_{V, \text{elec}}

Temperature Evolution of Heat Capacity in Diatomic Gases

1. Translational Contribution:

Always fully classical above a fraction of a Kelvin:

U_{\text{trans}} = \frac{3}{2} R T \implies C_{V, \text{trans}} = \frac{3}{2} R

2. Rotational Contribution:

At \(T \ll \theta_{\text{rot}}\) (e.g., \(T < 50\text{ K}\) for \(\text{H}_2\)): Rotational motion is quantum-mechanically frozen (\(C_{V, \text{rot}} \rightarrow 0\)).
At \(T \gg \theta_{\text{rot}}\): Fully active classical 2D rotation:

U_{\text{rot}} = R T \implies C_{V, \text{rot}} = R

3. Vibrational Contribution:

From the harmonic oscillator:

C_{V, \text{vib}} = R \left( \frac{\theta_{\text{vib}}}{T} \right)^2 \frac{e^{\theta_{\text{vib}}/T}}{(e^{\theta_{\text{vib}}/T} - 1)^2}

At \(T \ll \theta_{\text{vib}}\): Frozen zero-point motion (\(C_{V, \text{vib}} \rightarrow 0\)).
At \(T \gg \theta_{\text{vib}}\): Equipartition limit (1 kinetic + 1 potential quadratic term):

C_{V, \text{vib}} \rightarrow R

Summary of Plateaus for a Diatomic Molecule

Low \(T\) (\(T < \theta_{\text{rot}}\)): \(C_V = \frac{3}{2} R\) (monatomic-like behavior)
Intermediate \(T\) (\(\theta_{\text{rot}} < T \ll \theta_{\text{vib}}\)): \(C_V = \frac{3}{2} R + R = \frac{5}{2} R\) (room temperature for \(\text{N}_2, \text{O}_2, \text{CO}\))
High \(T\) (\(T \gg \theta_{\text{vib}}\)): \(C_V = \frac{5}{2} R + R = \frac{7}{2} R\)

This stepwise activation of degrees of freedom provided historical proof of energy quantization, resolving the failure of classical equipartition.

§8.7 Sackur-Tetrode Equation for Absolute Translational Entropy

The Sackur-Tetrode equation is the crowning achievement of early quantum statistical mechanics. Derived independently in 1912 by Otto Sackur and Hugo Tetrode, it predicts the absolute entropy of an ideal monatomic gas directly from fundamental physical constants (\(h, k_B, m\)).

Mathematical Derivation

For an indistinguishable monatomic gas, the total canonical partition function is:

Q = \frac{q_{\text{trans}}^N}{N!} g_{e, 0}^N = \frac{1}{N!} \left( \frac{V}{\Lambda^3} \right)^N g_{e, 0}^N

where \(\Lambda = \frac{h}{\sqrt{2\pi m k_B T}}\). Applying Stirling's approximation \(\ln N! \approx N \ln N - N\):

\ln Q = N \ln\left(\frac{V}{\Lambda^3}\right) - (N \ln N - N) + N \ln g_{e, 0} = N \left[ \ln\left(\frac{V}{N \Lambda^3}\right) + 1 + \ln g_{e, 0} \right]

Recall the statistical thermodynamic formula for entropy:

S = k_B \ln Q + k_B T \left( \frac{\partial \ln Q}{\partial T} \right)_{V, N}

Evaluate the temperature derivative:

\ln Q = N \ln V - 3 N \ln \Lambda - N \ln N + N + N \ln g_{e, 0}

Since \(\Lambda \propto T^{-1/2} \implies \ln \Lambda = -\frac{1}{2} \ln T + \text{const}\):

-3 N \ln \Lambda = \frac{3}{2} N \ln T + \text{const} \implies \left( \frac{\partial \ln Q}{\partial T} \right)_V = \frac{3 N}{2 T}

Substituting into \(S\):

S = k_B N \left[ \ln\left(\frac{V}{N \Lambda^3}\right) + 1 + \ln g_{e, 0} \right] + k_B T \left( \frac{3 N}{2 T} \right) = N k_B \left[ \ln\left(\frac{V}{N \Lambda^3}\right) + \frac{5}{2} + \ln g_{e, 0} \right]

Molar Form of the Sackur-Tetrode Equation

For 1 mole (\(N = N_A, N_A k_B = R\)) with ideal gas volume \(V/N_A = k_B T / P\):

S_m^\circ = R \left[ \ln\left( \frac{k_B T}{P^\circ \Lambda^3} \right) + \frac{5}{2} + \ln g_{e, 0} \right]

Substituting \(\Lambda = \frac{h}{\sqrt{2\pi m k_B T}}\) and expanding:

S_m^\circ(T, P^\circ) = R \left[ \frac{3}{2} \ln M + \frac{5}{2} \ln T - \ln P^\circ - 1.1517 + \ln g_{e, 0} \right]

where \(M\) is molar mass in \(\text{g/mol}\), \(T\) in Kelvin, and \(P^\circ\) in bar. The Sackur-Tetrode equation yields agreement to within \(0.01\%\) with experimental Third Law calorimetric entropy measurements for noble gases (He, Ne, Ar, Kr, Xe).

§8.8 Hindered Internal Rotations & Torsional Partition Functions

In flexible molecules (such as ethane \(\text{CH}_3-\text{CH}_3\), butane, or polymers), internal rotation around single \(\sigma\)-bonds is neither completely free nor completely rigid. It is a hindered (torsional) rotation.

The Torsional Potential Energy Function

For ethane, rotation of one methyl group relative to the other about the \(\text{C}-\text{C}\) axis encounters a 3-fold periodic potential barrier:

V(\phi) = \frac{V_3}{2} (1 - \cos 3\phi)

where \(\phi \in [0, 2\pi]\) is the torsional dihedral angle, \(V_3 \approx 12.0\text{ kJ/mol}\) is the barrier height, and the 3 minima correspond to the staggered conformations.

Limiting Behaviors

1. High Temperature or Low Barrier (\(k_B T \gg V_3\), Free Rotor Limit):

The potential barrier \(V_3\) is negligible compared to thermal energy. The internal rotation behaves as a 1D free rotor with reduced moment of inertia \(I_{\text{red}}\):

q_{\text{free}} = \frac{1}{\sigma_{\text{int}}} \left( \frac{8\pi^2 I_{\text{red}} k_B T}{h^2} \right)^{1/2} = \frac{1}{\sigma_{\text{int}}} \left( \frac{\pi T}{\theta_{\text{int}}} \right)^{1/2}

where \(\sigma_{\text{int}} = 3\) for ethane, and \(I_{\text{red}} = \frac{I_{A} I_{B}}{I_A + I_B}\). The molar heat capacity contribution is \(C_{V, \text{free}} = \frac{1}{2} R\).

2. Low Temperature or High Barrier (\(k_B T \ll V_3\), Torsional Oscillator Limit):

Molecules are trapped near the bottom of the potential wells (\(\phi \approx 0\)). Taylor expanding \(\cos(3\phi) \approx 1 - \frac{9\phi^2}{2}\):

V(\phi) \approx \frac{V_3}{2} \left[ 1 - \left( 1 - \frac{9\phi^2}{2} \right) \right] = \frac{9 V_3}{4} \phi^2 = \frac{1}{2} k_\phi \phi^2

The motion reduces to a simple harmonic torsional vibration with frequency:

\nu_{\text{tors}} = \frac{3}{2\pi} \sqrt{\frac{V_3}{2 I_{\text{red}}}}

The partition function is the standard 1D harmonic oscillator formula:

q_{\text{tors}} = \frac{1}{1 - e^{-h \nu_{\text{tors}} / k_B T}}

The heat capacity approaches \(C_{V, \text{tors}} = R \left( \frac{\theta}{T} \right)^2 \frac{e^{\theta/T}}{(e^{\theta/T} - 1)^2}\), which reaches \(R\) at intermediate temperatures.

The Pitzer-Gwinn Approximation for Intermediate Barriers

At ambient temperatures (\(k_B T \approx V_3\)), neither extreme is valid. Kenneth Pitzer and William Gwinn developed numerical tables and analytical approximations:

q_{\text{hindered}} = q_{\text{free}} \times \left( \frac{q_{\text{harm, QM}}}{q_{\text{harm, class}}} \right) I_0\left( \frac{V_3}{2 k_B T} \right) e^{-V_3 / 2 k_B T}

where \(I_0\) is the modified Bessel function. Accurate calculation of hindered rotation partition functions is critical for predicting conformational equilibria in pharmaceutical drug discovery and polymer physics.

## Advanced Mathematical Supplement: Dunham Coefficients & High-T Asymptotics

Dunham Expansions & High-Temperature Partition Function Asymptotics

In high-precision thermochemistry, idealized harmonic oscillator and rigid rotor models are insufficient. J. L. Dunham expanded the rovibrational energy levels in a bivariate double power series in \((v + 1/2)\) and \(J(J + 1)\):

E(v, J) = \sum_{l=0}^\infty \sum_{j=0}^\infty Y_{l, j} \left( v + \frac{1}{2} \right)^l [J(J + 1)]^j

where \(Y_{l, j}\) are the Dunham Coefficients:

\(Y_{1, 0} \approx \tilde{\omega}_e\) (fundamental harmonic frequency)
\(Y_{2, 0} \approx -\tilde{\omega}_e x_e\) (anharmonicity constant)
\(Y_{0, 1} \approx B_e\) (equilibrium rotational constant)
\(Y_{1, 1} \approx -\alpha_e\) (vibration-rotation coupling)
\(Y_{0, 2} \approx -D_e\) (centrifugal distortion)

High-Temperature Asymptotics via the Euler-Maclaurin Summation

To evaluate the partition function without truncation error, we apply the Euler-Maclaurin formula:

\sum_{k=a}^b f(k) = \int_a^b f(x) dx + \frac{f(a) + f(b)}{2} + \sum_{k=1}^m \frac{B_{2k}}{(2k)!} [f^{(2k-1)}(b) - f^{(2k-1)}(a)]

where \(B_{2k}\) are the Bernoulli numbers (\(B_2 = 1/6, B_4 = -1/30\)). For the rotational partition function of a linear rotor:

q_{\text{rot}}(T) = \sum_{J=0}^\infty (2J + 1) e^{-J(J+1) \theta_{\text{rot}} / T} = \frac{T}{\sigma \theta_{\text{rot}}} \left[ 1 + \frac{1}{3}\left(\frac{\theta_{\text{rot}}}{T}\right) + \frac{1}{15}\left(\frac{\theta_{\text{rot}}}{T}\right)^2 + \frac{4}{315}\left(\frac{\theta_{\text{rot}}}{T}\right)^3 + \dots \right]

The leading term \(\frac{T}{\sigma \theta_{\text{rot}}}\) is the classical integral.
The term \(+\frac{1}{3}\) is the first quantum correction.
Higher-order terms ensure six-figure accuracy down to cryogenic temperatures without explicit summation.

## Research Monograph: Path Integral Molecular Dynamics & Nuclear Quantum Effects

Nuclear Quantum Effects in Chemistry

In standard Born-Oppenheimer molecular dynamics, electrons are treated quantum mechanically, but nuclei are approximated as classical point masses following Newton's equations. However, for light nuclei (protons, deuterons, lithium), Nuclear Quantum Effects (NQEs)—including zero-point energy (ZPE) and quantum tunneling—dramatically alter hydrogen bond networks, proton transfer kinetics, and aqueous solvation.

Feynman-Kac Ring Polymer Isomorphism

David Chandler and Peter Wolynes showed that the quantum canonical partition function of a single quantum particle is mathematically isomorphic to the classical partition function of a ring polymer (necklace) consisting of \(P\) classical beads joined by harmonic springs:

Q_{\text{quantum}} = \lim_{P \rightarrow \infty} \left( \frac{m P k_B T}{2\pi \hbar^2} \right)^{P/2} \int d\mathbf{r}_1 \dots d\mathbf{r}_P \exp\left( -\beta_{\text{eff}} U_{\text{ring}}(\mathbf{r}_1, \dots, \mathbf{r}_P) \right)

where the ring polymer potential is:

U_{\text{ring}} = \sum_{k=1}^P \left[ \frac{1}{2} m \omega_P^2 (\mathbf{r}_{k+1} - \mathbf{r}_k)^2 + \frac{1}{P} V(\mathbf{r}_k) \right]

with cyclic boundary \(\mathbf{r}_{P+1} = \mathbf{r}_1\), effective spring constant \(\omega_P = \frac{\sqrt{P}}{\beta \hbar}\), and effective temperature \(T_{\text{eff}} = P T\).

Path Integral Molecular Dynamics (PIMD)

In PIMD, the \(P\) beads are propagated simultaneously using standard classical molecular dynamics algorithms:

1. Spatial Delocalization: The radius of gyration of the ring polymer measures the physical quantum spatial uncertainty (de Broglie wavepacket spread) of the nucleus.

2. Tunneling: In barrier crossing, the polymer stretches across the potential barrier, representing instanton tunneling trajectories.

PIMD simulations reveal that in liquid water, zero-point motion destabilizes hydrogen bonds while proton tunneling strengthens them; these two quantum effects cancel almost exactly at \(300\text{ K}\), explaining why classical models of water appear fortuitously successful.

Worked Problems & Step-by-Step Quantum Derivations

Multi-step solved problems covering Planck distribution, photoelectric kinetics, Compton shift, de Broglie wavelengths, uncertainty relations, and Hermitian operator commutation algebra.

Advanced Example 8.1: Thermal de Broglie Wavelength and q_trans for Argon at 300 K

For argon gas (\(\text{Ar}\), atomic mass \(M = 39.948\text{ g/mol}\)) at \(T = 300.0\text{ K}\) and \(P = 1.00\text{ bar}\):

Calculate the single-atom mass \(m\) and the thermal de Broglie wavelength \(\Lambda\).
Calculate the average volume per atom \(V/N\) and compare the thermal de Broglie wavelength to the average interatomic spacing \(d = (V/N)^{1/3}\).
Compute the single-particle translational partition function \(q_{\text{trans}}\) in a volume of \(1.00\text{ dm}^3\).

Comprehensive Multi-Step Solution:

Step 1: Mass and Thermal de Broglie Wavelength

1. Single-atom mass \(m\):

m = \frac{M}{N_A} = \frac{0.039948\text{ kg/mol}}{6.02214 \times 10^{23}\text{ mol}^{-1}} \approx 6.6335 \times 10^{-26}\text{ kg}

2. Thermal de Broglie wavelength:

\Lambda = \frac{h}{\sqrt{2\pi m k_B T}}

Evaluate the denominator:

2\pi m k_B T = 2 \pi (6.6335 \times 10^{-26}\text{ kg}) (1.38065 \times 10^{-23}\text{ J/K}) (300.0\text{ K})

= 2 \pi \times 6.6335 \times 1.38065 \times 300 \times 10^{-49} \approx 1.7258 \times 10^{-45}\text{ kg}^2\cdot\text{m}^2/\text{s}^2

Taking square root:

\sqrt{2\pi m k_B T} \approx 4.1543 \times 10^{-23}\text{ kg}\cdot\text{m/s}

Thus:

\Lambda = \frac{6.62607 \times 10^{-34}\text{ J}\cdot\text{s}}{4.1543 \times 10^{-23}\text{ kg}\cdot\text{m/s}} \approx 1.595 \times 10^{-11}\text{ m} = 0.01595\text{ nm} = 0.1595\text{ \AA}

Step 2: Average Interatomic Distance Comparison

At \(P = 1.00\text{ bar} = 10^5\text{ Pa}\) and \(T = 300.0\text{ K}\):

\frac{V}{N} = \frac{k_B T}{P} = \frac{(1.38065 \times 10^{-23}\text{ J/K})(300\text{ K})}{10^5\text{ Pa}} = 4.142 \times 10^{-26}\text{ m}^3

The average interatomic spacing is:

d = \left( \frac{V}{N} \right)^{1/3} = (41.42 \times 10^{-27}\text{ m}^3)^{1/3} \approx 3.46 \times 10^{-9}\text{ m} = 3.46\text{ nm} = 34.6\text{ \AA}

Comparing scales:

\frac{\Lambda}{d} = \frac{0.01595\text{ nm}}{3.46\text{ nm}} \approx 4.6 \times 10^{-3} \ll 1

Because \(\Lambda \ll d\), quantum wavepacket overlap is completely negligible. Classical Maxwell-Boltzmann statistics is accurate to within 1 part in \(10^7\).

Step 3: Translational Partition Function in \(V = 1.00\text{ dm}^3 = 10^{-3}\text{ m}^3\)

q_{\text{trans}} = \frac{V}{\Lambda^3} = \frac{1.00 \times 10^{-3}\text{ m}^3}{(1.595 \times 10^{-11}\text{ m})^3} = \frac{1.00 \times 10^{-3}}{4.057 \times 10^{-33}} \approx 2.46 \times 10^{29}

A single argon atom has over \(10^{29}\) thermally accessible translational quantum states in a 1-liter flask at room temperature.

Advanced Example 8.2: Rotational Partition Function and Characteristic Temperature for HCl and CO

Given spectroscopic constants:

Hydrogen chloride (\(^{1}\text{H}^{35}\text{Cl}\)): \(B = 10.593\text{ cm}^{-1}\), \(\sigma = 1\)
Carbon monoxide (\(^{12}\text{C}^{16}\text{O}\)): \(B = 1.931\text{ cm}^{-1}\), \(\sigma = 1\)

Calculate the characteristic rotational temperature \(\theta_{\text{rot}}\) for \(\text{HCl}\) and \(\text{CO}\).
Compute the rotational partition function \(q_{\text{rot}}\) for both molecules at \(T = 300.0\text{ K}\).
Determine the rotational quantum number \(J_{\text{max}}\) that corresponds to the most populated rotational level at \(300\text{ K}\) for both gases.

Comprehensive Multi-Step Solution:

Step 1: Characteristic Rotational Temperature \(\theta_{\text{rot}}\)

The characteristic rotational temperature is:

\theta_{\text{rot}} = \frac{h c B}{k_B}

Using \(\frac{h c}{k_B} \approx 1.43878\text{ cm}\cdot\text{K}\):

1. For \(\text{HCl}\):

\theta_{\text{rot}}(\text{HCl}) = 1.43878 \times 10.593\text{ cm}^{-1} \approx 15.24\text{ K}

2. For \(\text{CO}\):

\theta_{\text{rot}}(\text{CO}) = 1.43878 \times 1.931\text{ cm}^{-1} \approx 2.778\text{ K}

Step 2: Rotational Partition Function at \(300.0\text{ K}\)

Since \(T \gg \theta_{\text{rot}}\) for both molecules, the high-temperature approximation \(q_{\text{rot}} = \frac{T}{\sigma \theta_{\text{rot}}}\) is valid:

1. For \(\text{HCl}\) (\(\sigma = 1\)):

q_{\text{rot}} = \frac{300.0}{1 \times 15.24} \approx 19.685

(Including the first-order quantum correction \(q_{\text{rot}} = \frac{T}{\theta} + \frac{1}{3} = 19.685 + 0.333 \approx 20.02\)).

2. For \(\text{CO}\) (\(\sigma = 1\)):

q_{\text{rot}} = \frac{300.0}{1 \times 2.778} \approx 107.99

Step 3: Most Populated Rotational Level \(J_{\text{max}}\)

The population of rotational level \(J\) is:

N_J \propto (2J + 1) \exp\left( -J(J+1) \frac{\theta_{\text{rot}}}{T} \right)

Treating \(J\) as continuous and maximizing \(N(J)\):

\frac{d N_J}{dJ} = \left[ 2 - (2J + 1)^2 \frac{\theta_{\text{rot}}}{T} \right] e^{-J(J+1)\theta_{\text{rot}}/T} = 0

(2J + 1)^2 \frac{\theta_{\text{rot}}}{T} = 2 \implies 2J + 1 = \sqrt{\frac{2 T}{\theta_{\text{rot}}}} \implies J_{\text{max}} = \sqrt{\frac{T}{2 \theta_{\text{rot}}}} - \frac{1}{2}

1. For \(\text{HCl}\):

J_{\text{max}} = \sqrt{\frac{300}{2 \times 15.24}} - 0.5 = \sqrt{9.8425} - 0.5 = 3.137 - 0.5 \approx 2.64 \implies J = 3

2. For \(\text{CO}\):

J_{\text{max}} = \sqrt{\frac{300}{2 \times 2.778}} - 0.5 = \sqrt{53.996} - 0.5 = 7.348 - 0.5 \approx 6.85 \implies J = 7

Advanced Example 8.3: Nuclear Spin Statistics of Ortho and Para Molecular Hydrogen

For molecular hydrogen (\(\text{H}_2\), \(\theta_{\text{rot}} = 85.3\text{ K}\)):

Write the explicit rotational partition function for para-\(\text{H}_2\) (\(q_{\text{para}}\)) and ortho-\(\text{H}_2\) (\(q_{\text{ortho}}\)).
Calculate the equilibrium ratio \(N_{\text{ortho}} / N_{\text{para}}\) at \(T = 50.0\text{ K}\) and at \(T = 300.0\text{ K}\).
Explain why liquid hydrogen (\(T = 20.3\text{ K}\)) must be catalytically converted from normal hydrogen to para-hydrogen before long-term storage.

Comprehensive Multi-Step Solution:

Step 1: Explicit Rotational Partition Functions

For protons with spin \(I = 1/2\):

Para-\(\text{H}_2\) (spin singlet, \(g_{\text{spin}} = 1\)): restricted to even \(J \in \{0, 2, 4, \dots\}\)

q_{\text{para}} = \sum_{J = 0, 2, 4, \dots} (2J + 1) e^{-J(J+1)\theta_{\text{rot}}/T}

Ortho-\(\text{H}_2\) (spin triplet, \(g_{\text{spin}} = 3\)): restricted to odd \(J \in \{1, 3, 5, \dots\}\)

q_{\text{ortho}} = 3 \sum_{J = 1, 3, 5, \dots} (2J + 1) e^{-J(J+1)\theta_{\text{rot}}/T}

Step 2: Equilibrium Ratio at \(50\text{ K}\) and \(300\text{ K}\)

The equilibrium ratio of populations is:

\frac{N_{\text{ortho}}}{N_{\text{para}}} = \frac{q_{\text{ortho}}}{q_{\text{para}}} = \frac{3 \sum_{J \text{ odd}} (2J + 1) e^{-J(J+1)\theta_{\text{rot}}/T}}{\sum_{J \text{ even}} (2J + 1) e^{-J(J+1)\theta_{\text{rot}}/T}}

Given \(\theta_{\text{rot}} = 85.3\text{ K}\):

1. At \(T = 50.0\text{ K}\):

\(\theta / T = 85.3 / 50.0 = 1.706\)
Even terms:
\(J = 0\): \(1 \times e^0 = 1.000\)
\(J = 2\): \(5 \times e^{-6 \times 1.706} = 5 \times e^{-10.236} \approx 5 \times (3.58 \times 10^{-5}) \approx 0.00018\)
\(q_{\text{para}} \approx 1.00018\)
Odd terms:
\(J = 1\): \(3 \times e^{-2 \times 1.706} = 3 \times e^{-3.412} \approx 3 \times 0.03297 \approx 0.0989\)
\(J = 3\): \(7 \times e^{-12 \times 1.706} = 7 \times e^{-20.47} \approx 0.0000\)
\(q_{\text{ortho}} \approx 3 \times 0.0989 = 0.2967\)

Ratio:

\frac{N_{\text{ortho}}}{N_{\text{para}}} = \frac{0.2967}{1.00018} \approx 0.297

At \(50\text{ K}\), hydrogen is roughly 77% para and 23% ortho.

2. At \(T = 300.0\text{ K}\):

At high temperatures (\(T \gg \theta_{\text{rot}}\)), both odd and even sums converge to the identical integral value \(\frac{T}{2 \theta_{\text{rot}}}\):

\lim_{T \rightarrow \infty} \frac{q_{\text{ortho}}}{q_{\text{para}}} = \frac{3 \times \frac{T}{2\theta_{\text{rot}}}}{1 \times \frac{T}{2\theta_{\text{rot}}}} = 3.000

Normal hydrogen at room temperature is exactly 75% ortho and 25% para (3:1 ratio).

Step 3: Cryogenic Storage and Catalysis

When room-temperature hydrogen (75% ortho) is liquefied at \(20.3\text{ K}\), the nuclear spin conversion \(J = 1 \rightarrow J = 0\) is spin-forbidden and proceeds very slowly in the absence of a magnetic catalyst. However, the ortho-to-para conversion is exothermic:

\Delta H = E(J=1) - E(J=0) = 2 k_B \theta_{\text{rot}} = 2 \times (1.38 \times 10^{-23}) \times 85.3 \approx 2.35 \times 10^{-21}\text{ J/molecule} \approx 1.42\text{ kJ/mol}

This heat of conversion (\(1.42\text{ kJ/mol}\)) exceeds the heat of vaporization of liquid hydrogen (\(\Delta H_{\text{vap}} \approx 0.90\text{ kJ/mol}\)). Without a paramagnet catalyst (such as ferric oxide \(\text{Fe}_2\text{O}_3\) or activated carbon), uncatalyzed ortho-to-para relaxation in storage tanks releases enough heat to spontaneously boil off up to 50% of the liquid hydrogen within days!

Advanced Example 8.4: Vibrational Partition Function and Excited State Populations for N2 and I2

Given vibrational wavenumbers:

Nitrogen (\(\text{N}_2\)): \(\tilde{\nu} = 2358.6\text{ cm}^{-1}\)
Iodine vapor (\(\text{I}_2\)): \(\tilde{\nu} = 214.5\text{ cm}^{-1}\)

Calculate the characteristic vibrational temperatures \(\theta_{\text{vib}}\) for \(\text{N}_2\) and \(\text{I}_2\).
Calculate the vibrational partition function \(q_{\text{vib}}\) (relative to \(v = 0\)) for both molecules at \(300.0\text{ K}\) and at \(1000.0\text{ K}\).
Determine the percentage fraction of molecules in vibrationally excited states (\(v \ge 1\)) for both gases at \(300\text{ K}\) and \(1000\text{ K}\).

Comprehensive Multi-Step Solution:

Step 1: Characteristic Vibrational Temperatures

Using \(\theta_{\text{vib}} = \frac{h c \tilde{\nu}}{k_B} = 1.43878 \times \tilde{\nu}\):

1. For \(\text{N}_2\):

\theta_{\text{vib}}(\text{N}_2) = 1.43878 \times 2358.6\text{ cm}^{-1} \approx 3393.5\text{ K}

2. For \(\text{I}_2\):

\theta_{\text{vib}}(\text{I}_2) = 1.43878 \times 214.5\text{ cm}^{-1} \approx 308.6\text{ K}

Step 2: Vibrational Partition Functions

The partition function relative to the ground state is:

q_{\text{vib}} = \frac{1}{1 - e^{-\theta_{\text{vib}} / T}}

1. At \(T = 300.0\text{ K}\):

For \(\text{N}_2\): \(\theta / T = 3393.5 / 300 = 11.312\)

e^{-11.312} \approx 1.22 \times 10^{-5} \implies q_{\text{vib}} = \frac{1}{1 - 1.22 \times 10^{-5}} \approx 1.000012

For \(\text{I}_2\): \(\theta / T = 308.6 / 300 = 1.0287\)

e^{-1.0287} \approx 0.35747 \implies q_{\text{vib}} = \frac{1}{1 - 0.35747} = \frac{1}{0.64253} \approx 1.5563

2. At \(T = 1000.0\text{ K}\):

For \(\text{N}_2\): \(\theta / T = 3393.5 / 1000 = 3.3935\)

e^{-3.3935} \approx 0.03359 \implies q_{\text{vib}} = \frac{1}{1 - 0.03359} \approx 1.0348

For \(\text{I}_2\): \(\theta / T = 308.6 / 1000 = 0.3086\)

e^{-0.3086} \approx 0.73447 \implies q_{\text{vib}} = \frac{1}{1 - 0.73447} = \frac{1}{0.26553} \approx 3.766

Step 3: Fraction in Excited States (\(v \ge 1\))

The ground-state fraction is \(P_0 = \frac{1}{q_{\text{vib}}} = 1 - e^{-\theta_{\text{vib}}/T}\). The excited-state fraction is:

P(v \ge 1) = 1 - P_0 = e^{-\theta_{\text{vib}} / T}

1. For \(\text{N}_2\):

At \(300\text{ K}\): \(P(v \ge 1) = 1.22 \times 10^{-5} \approx 0.0012\%\) (Virtually 100% in \(v = 0\))
At \(1000\text{ K}\): \(P(v \ge 1) \approx 0.0336 \approx 3.36\%\)

2. For \(\text{I}_2\):

At \(300\text{ K}\): \(P(v \ge 1) \approx 0.3575 \approx 35.75\%\)
At \(1000\text{ K}\): \(P(v \ge 1) \approx 0.7345 \approx 73.45\%\)

Because iodine has a weak bond and heavy atoms, its vibrational frequency is low, resulting in significant vibrational excitation even at room temperature.

Advanced Example 8.5: Absolute Standard Entropy of Argon at 298.15 K via Sackur-Tetrode

For argon gas (\(M = 39.948\text{ g/mol}\), \(g_{e, 0} = 1\)) at \(T = 298.15\text{ K}\) and standard pressure \(P^\circ = 1.000\text{ bar} = 10^5\text{ Pa}\):

Compute the thermal de Broglie wavelength \(\Lambda\) at \(298.15\text{ K}\).
Use the Sackur-Tetrode equation to calculate the absolute molar standard entropy \(S_m^\circ\).
Compare the calculated value with the experimental Third Law calorimetric entropy \(S_m^\circ = 154.84\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}\).

Comprehensive Multi-Step Solution:

Step 1: Thermal de Broglie Wavelength

At \(T = 298.15\text{ K}\) with atomic mass \(m = \frac{0.039948}{6.02214 \times 10^{23}} = 6.6335 \times 10^{-26}\text{ kg}\):

\Lambda = \frac{h}{\sqrt{2\pi m k_B T}}

2\pi m k_B T = 2\pi (6.6335 \times 10^{-26})(1.38065 \times 10^{-23})(298.15) \approx 1.7152 \times 10^{-45}\text{ kg}^2\cdot\text{m}^2/\text{s}^2

\sqrt{2\pi m k_B T} \approx 4.1415 \times 10^{-23}\text{ kg}\cdot\text{m/s}

\Lambda = \frac{6.62607 \times 10^{-34}}{4.1415 \times 10^{-23}} \approx 1.5999 \times 10^{-11}\text{ m}

Step 2: Sackur-Tetrode Entropy Calculation

The Sackur-Tetrode equation is:

S_m^\circ = R \left[ \ln\left( \frac{k_B T}{P^\circ \Lambda^3} \right) + \frac{5}{2} + \ln g_{e, 0} \right]

Evaluate the argument inside the logarithm:

k_B T = (1.38065 \times 10^{-23})(298.15) \approx 4.1164 \times 10^{-21}\text{ J}

\Lambda^3 = (1.5999 \times 10^{-11}\text{ m})^3 \approx 4.0952 \times 10^{-33}\text{ m}^3

P^\circ \Lambda^3 = (10^5\text{ Pa})(4.0952 \times 10^{-33}\text{ m}^3) = 4.0952 \times 10^{-28}\text{ J}

Taking the ratio:

\frac{k_B T}{P^\circ \Lambda^3} = \frac{4.1164 \times 10^{-21}}{4.0952 \times 10^{-28}} \approx 1.00518 \times 10^7

Taking the natural logarithm:

\ln(1.00518 \times 10^7) = \ln(1.00518) + 7 \ln(10) = 0.00516 + 7(2.302585) = 0.00516 + 16.11810 = 16.12326

Adding the constant term \(\frac{5}{2} = 2.50000\):

\text{Bracket sum} = 16.12326 + 2.50000 = 18.62326

Multiplying by the universal gas constant \(R = 8.31446\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}\):

S_m^\circ = 8.31446 \times 18.62326 \approx 154.842\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}

Step 3: Comparison with Experiment

Calculated \(S_m^\circ = 154.84\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}\)
Experimental \(S_m^\circ = 154.84\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}\)

The agreement is exact to four significant figures, confirming the validity of quantum statistical mechanics.

Advanced Example 8.6: Electronic Partition Function and Heat Capacity of Nitric Oxide

Nitric oxide (\(\text{NO}\)) has an open-shell electronic ground configuration with two low-lying states:

Ground level \(^2\Pi_{1/2}\): degeneracy \(g_0 = 2\), \(\varepsilon_0 = 0\)
Excited level \(^2\Pi_{3/2}\): degeneracy \(g_1 = 2\), \(\varepsilon_1 = 121.1\text{ cm}^{-1}\) (\(\theta_{\text{elec}} = 174.2\text{ K}\))

Higher electronic states lie far above in the UV.

Write the electronic partition function \(q_{\text{elec}}(T)\).
Derive the electronic contribution to internal energy \(U_{\text{elec}}(T)\) and molar heat capacity \(C_{V, \text{elec}}(T)\).
Calculate \(C_{V, \text{elec}}\) at \(T = 50\text{ K}\), \(T = 174\text{ K}\), and \(T = 1000\text{ K}\).

Comprehensive Multi-Step Solution:

Step 1: Electronic Partition Function

The states are \(\varepsilon_0 = 0\) and \(\varepsilon_1 = k_B \theta_{\text{elec}}\) with \(g_0 = g_1 = 2\).

q_{\text{elec}}(T) = g_0 + g_1 e^{-\theta_{\text{elec}} / T} = 2 \left( 1 + e^{-\theta_{\text{elec}} / T} \right)

Step 2: Internal Energy and Heat Capacity

1. Internal Energy:

U_{\text{elec}} = R T^2 \left( \frac{\partial \ln q_{\text{elec}}}{\partial T} \right) = R T^2 \frac{1}{1 + e^{-\theta/T}} \left( \frac{\theta}{T^2} e^{-\theta/T} \right) = R \theta \frac{e^{-\theta/T}}{1 + e^{-\theta/T}} = \frac{R \theta}{e^{\theta/T} + 1}

2. Molar Heat Capacity:

C_{V, \text{elec}} = \frac{d U_{\text{elec}}}{dT} = R \theta \left( -\frac{1}{(e^{\theta/T} + 1)^2} \right) e^{\theta/T} \left( -\frac{\theta}{T^2} \right) = R \left( \frac{\theta}{T} \right)^2 \frac{e^{\theta/T}}{(e^{\theta/T} + 1)^2}

This has the exact mathematical form of a two-level Schottky heat capacity anomaly!

Step 3: Numerical Evaluation at Specified Temperatures

Given \(R = 8.3145\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}\) and \(\theta = 174.2\text{ K}\):

1. At \(T = 50.0\text{ K}\):

\(x = \theta / T = 174.2 / 50.0 = 3.484\)
\(e^x = e^{3.484} \approx 32.59\)
\(\frac{e^x}{(e^x + 1)^2} = \frac{32.59}{(33.59)^2} = \frac{32.59}{1128.3} \approx 0.02888\)
\(C_{V, \text{elec}} = 8.3145 \times (3.484)^2 \times 0.02888 = 8.3145 \times 12.138 \times 0.02888 \approx 2.91\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}\)

2. At \(T = 174.2\text{ K}\) (\(T = \theta\)):

\(x = 1.000\)
\(e^1 \approx 2.7183\)
\(\frac{e}{(e + 1)^2} = \frac{2.7183}{(3.7183)^2} = \frac{2.7183}{13.826} \approx 0.1966\)
\(C_{V, \text{elec}} = 8.3145 \times (1.0)^2 \times 0.1966 \approx 1.63\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}\)

3. At \(T = 1000.0\text{ K}\):

\(x = 174.2 / 1000 = 0.1742\)
\(e^x \approx 1 + 0.1742 = 1.1742\)
\(\frac{e^x}{(e^x + 1)^2} \approx \frac{1}{(2)^2} = 0.25\)
\(C_{V, \text{elec}} \approx 8.3145 \times (0.1742)^2 \times 0.25 = 8.3145 \times 0.03035 \times 0.25 \approx 0.063\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}\)

As \(T \rightarrow \infty\), both levels are equally populated (\(P_0 = P_1 = 0.5\)), so no further thermal energy can be absorbed, causing \(C_{V, \text{elec}} \rightarrow 0\).

Advanced Example 8.7: Temperature-Dependent Molar Heat Capacity of Carbon Dioxide

Carbon dioxide (\(\text{CO}_2\)) is a linear triatomic molecule (\(3N - 5 = 4\) normal modes) with vibrational wavenumbers:

Symmetric stretch (\(\tilde{\nu}_1\)): \(1388\text{ cm}^{-1}\) (non-degenerate)
Bending (\(\tilde{\nu}_2\)): \(667\text{ cm}^{-1}\) (doubly degenerate, \(g_2 = 2\))
Asymmetric stretch (\(\tilde{\nu}_3\)): \(2349\text{ cm}^{-1}\) (non-degenerate)

Calculate the characteristic vibrational temperatures \(\theta_1, \theta_2, \theta_3\).
Calculate the translational, rotational, and vibrational molar heat capacities at \(T = 300.0\text{ K}\).
Compute the total \(C_{V, \text{molar}}\) and heat capacity ratio \(\gamma = C_P / C_V\) at \(300\text{ K}\) and compare to the high-temperature equipartition limit.

Comprehensive Multi-Step Solution:

Step 1: Characteristic Vibrational Temperatures

Using \(\theta = 1.43878 \times \tilde{\nu}\):

\(\theta_1 = 1.43878 \times 1388\text{ cm}^{-1} \approx 1997\text{ K}\)
\(\theta_2 = 1.43878 \times 667\text{ cm}^{-1} \approx 960\text{ K}\) (doubly degenerate)
\(\theta_3 = 1.43878 \times 2349\text{ cm}^{-1} \approx 3380\text{ K}\)

Step 2: Heat Capacity Contributions at \(300.0\text{ K}\)

1. Translation (3 degrees of freedom):

C_{V, \text{trans}} = \frac{3}{2} R \approx 1.5 \times 8.3145 \approx 12.472\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}

2. Rotation (Linear molecule, 2 degrees of freedom):

Since \(T = 300\text{ K} \gg \theta_{\text{rot}} \approx 0.56\text{ K}\), rotation is fully classical:

C_{V, \text{rot}} = R \approx 8.3145\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}

3. Vibration (Einstein function \(C_{\text{vib}}(x) = R x^2 \frac{e^x}{(e^x - 1)^2}\) with \(x = \theta / T\)):

Mode 1 (\(\theta_1 = 1997\text{ K}\)): \(x_1 = 1997 / 300 = 6.657\)

C_{V, 1} = R (6.657)^2 \frac{e^{6.657}}{(e^{6.657} - 1)^2} \approx R (44.31) e^{-6.657} \approx R (44.31)(0.00128) \approx 0.057 R \approx 0.47\text{ J/mol}\cdot\text{K}

Mode 2 (Bending, \(\theta_2 = 960\text{ K}\), degeneracy 2): \(x_2 = 960 / 300 = 3.200\)

\frac{e^{3.20}}{(e^{3.20} - 1)^2} = \frac{24.53}{(23.53)^2} = \frac{24.53}{553.8} \approx 0.0443

C_{V, 2} = 2 \times R (3.20)^2 (0.0443) = 2 \times R (10.24)(0.0443) \approx 0.907 R \approx 7.54\text{ J/mol}\cdot\text{K}

Mode 3 (\(\theta_3 = 3380\text{ K}\)): \(x_3 = 3380 / 300 = 11.27\)

e^{-11.27} \approx 1.27 \times 10^{-5} \implies C_{V, 3} \approx 0.0016 R \approx 0.01\text{ J/mol}\cdot\text{K}

Total vibrational heat capacity at \(300\text{ K}\):

C_{V, \text{vib}} = 0.47 + 7.54 + 0.01 \approx 8.02\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1} \approx 0.965 R

Step 3: Total Heat Capacity and Ratio \(\gamma\)

Total \(C_V\):

C_V = C_{V, \text{trans}} + C_{V, \text{rot}} + C_{V, \text{vib}} = 12.47 + 8.31 + 8.02 = 28.80\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1} \approx 3.465 R

Molar \(C_P = C_V + R = 28.80 + 8.31 = 37.11\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1} \approx 4.465 R\). Heat capacity ratio:

\gamma = \frac{C_P}{C_V} = \frac{37.11}{28.80} \approx 1.289

Comparison with high-temperature equipartition limit: At \(T \rightarrow \infty\), all 4 vibrational modes are fully active (\(C_{V, \text{vib}} \rightarrow 4 R\)):

C_{V, \text{classical}} = \frac{3}{2} R + R + 4 R = \frac{13}{2} R = 6.5 R \approx 54.04\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}

\gamma_{\text{classical}} = \frac{7.5 R}{6.5 R} = \frac{15}{13} \approx 1.154

At room temperature, \(\text{CO}_2\) is midway between rigid linear rotor behavior (\(\gamma = 7/5 = 1.40\)) and fully active vibrational equipartition (\(\gamma = 1.154\)) because only the low-frequency bending mode is partially active.

Advanced Example 8.8: Torsional Partition Function and Hindered Internal Rotation in Ethane

For ethane (\(\text{CH}_3-\text{CH}_3\)) at \(T = 300.0\text{ K}\): The internal rotation barrier is \(V_3 = 12.0\text{ kJ/mol}\). The reduced moment of inertia for methyl group rotation is \(I_{\text{red}} = 2.65 \times 10^{-47}\text{ kg}\cdot\text{m}^2\), and the symmetry number is \(\sigma_{\text{int}} = 3\).

Calculate the free-rotor partition function \(q_{\text{free}}\) at \(300\text{ K}\).
Calculate the torsional vibrational frequency \(\nu_{\text{tors}}\) and the characteristic vibrational temperature \(\theta_{\text{tors}}\) in the harmonic oscillator limit.
Compute the harmonic torsional partition function \(q_{\text{tors}}\) and compare it with \(q_{\text{free}}\).

Comprehensive Multi-Step Solution:

Step 1: Free-Rotor Partition Function

The 1D free-rotor formula is:

q_{\text{free}} = \frac{1}{\sigma_{\text{int}}} \left( \frac{8\pi^2 I_{\text{red}} k_B T}{h^2} \right)^{1/2}

Evaluate the argument:

8\pi^2 I_{\text{red}} k_B T = 8\pi^2 (2.65 \times 10^{-47}\text{ kg}\cdot\text{m}^2)(1.38065 \times 10^{-23}\text{ J/K})(300\text{ K})

= 8\pi^2 \times 2.65 \times 1.38065 \times 300 \times 10^{-70} \approx 8.665 \times 10^{-67}\text{ J}^2\cdot\text{s}^2

Divide by \(h^2 = (6.626 \times 10^{-34})^2 \approx 4.3905 \times 10^{-67}\):

\frac{8\pi^2 I_{\text{red}} k_B T}{h^2} = \frac{8.665 \times 10^{-67}}{4.3905 \times 10^{-67}} \approx 1.9736

Taking square root:

\sqrt{1.9736} \approx 1.4048

Dividing by \(\sigma_{\text{int}} = 3\):

q_{\text{free}} = \frac{1.4048}{3} \approx 0.468

Step 2: Torsional Vibrational Frequency in the Harmonic Limit

The barrier per molecule is:

V_3 = \frac{12000\text{ J/mol}}{6.02214 \times 10^{23}\text{ mol}^{-1}} \approx 1.9926 \times 10^{-20}\text{ J}

The effective torsional spring constant is \(k_\phi = \frac{9 V_3}{2}\):

k_\phi = 4.5 \times (1.9926 \times 10^{-20}\text{ J}) \approx 8.967 \times 10^{-20}\text{ J/rad}^2

The torsional frequency is:

\nu_{\text{tors}} = \frac{1}{2\pi} \sqrt{\frac{k_\phi}{I_{\text{red}}}} = \frac{1}{2\pi} \sqrt{\frac{8.967 \times 10^{-20}\text{ J}}{2.65 \times 10^{-47}\text{ kg}\cdot\text{m}^2}} = \frac{1}{2\pi} \sqrt{3.3837 \times 10^{27}} \approx \frac{5.817 \times 10^{13}}{2\pi} \approx 9.258 \times 10^{12}\text{ Hz}

In wavenumber:

\tilde{\nu}_{\text{tors}} = \frac{\nu_{\text{tors}}}{c} = \frac{9.258 \times 10^{12}\text{ s}^{-1}}{2.9979 \times 10^{10}\text{ cm/s}} \approx 308.8\text{ cm}^{-1}

Characteristic temperature:

\theta_{\text{tors}} = \frac{h \nu_{\text{tors}}}{k_B} = 1.43878 \times 308.8 \approx 444.3\text{ K}

Step 3: Harmonic Torsional Partition Function

At \(T = 300.0\text{ K}\):

\frac{\theta_{\text{tors}}}{T} = \frac{444.3}{300.0} = 1.481

e^{-\theta_{\text{tors}} / T} = e^{-1.481} \approx 0.2274

The harmonic oscillator partition function is:

q_{\text{tors}} = \frac{1}{1 - e^{-\theta_{\text{tors}} / T}} = \frac{1}{1 - 0.2274} = \frac{1}{0.7726} \approx 1.294

Because \(V_3 / k_B T = \frac{12000}{8.314 \times 300} \approx 4.81 > 1\), the potential barrier is substantial compared to thermal energy. Ethane at room temperature behaves predominantly as a torsional harmonic oscillator (\(q \approx 1.29\)) rather than a free internal rotor.

Advanced Example 8.9: Rotational Relaxation and Ortho-to-Para Conversion Rates

Consider a cryogenic sample of molecular hydrogen \(\text{H}_2\) cooled from \(300\text{ K}\) to \(20.0\text{ K}\).

State the fraction of ortho-\(\text{H}_2\) and para-\(\text{H}_2\) at \(300\text{ K}\) (normal hydrogen) and at \(20.0\text{ K}\) at true thermal equilibrium.
In the absence of a catalyst, ortho-to-para conversion follows second-order kinetics driven by nuclear dipole-dipole interactions between colliding molecules:

-\frac{d x_{\text{ortho}}}{dt} = k_{\text{conv}} x_{\text{ortho}}^2

with rate constant \(k_{\text{conv}} = 0.019\text{ hr}^{-1}\) in the liquid phase. Calculate the time required for the ortho fraction to drop from \(0.75\) to \(0.50\).

If the enthalpy of ortho-to-para conversion is \(\Delta H_{\text{conv}} = 1.42\text{ kJ/mol}\) and the heat of vaporization of liquid \(\text{H}_2\) is \(\Delta H_{\text{vap}} = 0.904\text{ kJ/mol}\), calculate the mass of liquid hydrogen boiled off per mole of unconverted ortho-hydrogen that undergoes spontaneous relaxation in an unvented storage dewar.

Comprehensive Multi-Step Solution:

Step 1: Equilibrium Ortho and Para Fractions

At \(T = 300\text{ K}\) (high temperature limit):

x_{\text{ortho}} = \frac{3}{3 + 1} = 0.75 = 75\%, \quad x_{\text{para}} = \frac{1}{3 + 1} = 0.25 = 25\%

At \(T = 20.0\text{ K}\) (\(\theta_{\text{rot}} = 85.3\text{ K} \implies \theta / T = 4.265\)):

q_{\text{para}} \approx 1.000, \quad q_{\text{ortho}} = 3(3) e^{-2(4.265)} = 9 e^{-8.53} \approx 9(0.000197) \approx 0.00177

Equilibrium fractions:

x_{\text{ortho}} \approx \frac{0.00177}{1.00177} \approx 0.18\%, \quad x_{\text{para}} \approx 99.82\%

True equilibrium liquid hydrogen is virtually 100% pure para-\(\text{H}_2\).

Step 2: Uncatalyzed Conversion Kinetics

The rate law is:

-\frac{d x}{dt} = k_{\text{conv}} x^2 \implies \int_{x_0}^{x_t} \frac{dx}{x^2} = -k_{\text{conv}} \int_0^t dt

\frac{1}{x_t} - \frac{1}{x_0} = k_{\text{conv}} t \implies t = \frac{1}{k_{\text{conv}}} \left( \frac{1}{x_t} - \frac{1}{x_0} \right)

Substitute \(x_0 = 0.75\), \(x_t = 0.50\), and \(k_{\text{conv}} = 0.019\text{ hr}^{-1}\):

\frac{1}{0.50} - \frac{1}{0.75} = 2.000 - 1.333 = 0.667

t = \frac{0.667}{0.019\text{ hr}^{-1}} \approx 35.1\text{ hours}

Without a catalyst, it takes over 35 hours for the ortho concentration to decrease from 75% to 50%.

Step 3: Liquid Boil-Off Calculation

Relaxation of 1 mole of ortho-\(\text{H}_2\) to para-\(\text{H}_2\) releases:

Q_{\text{released}} = \Delta H_{\text{conv}} = 1.42\text{ kJ/mol}

Each mole of liquid hydrogen vaporized absorbs:

\Delta H_{\text{vap}} = 0.904\text{ kJ/mol}

The moles of liquid hydrogen boiled off per mole of relaxing ortho-\(\text{H}_2\) is:

n_{\text{boiled}} = \frac{Q_{\text{released}}}{\Delta H_{\text{vap}}} = \frac{1.42\text{ kJ}}{0.904\text{ kJ/mol}} \approx 1.571\text{ moles of }\text{H}_2

The mass of liquid hydrogen boiled off is:

m_{\text{boiled}} = n_{\text{boiled}} \times M(\text{H}_2) = 1.571\text{ mol} \times 2.016\text{ g/mol} \approx 3.17\text{ grams}

Because \(\Delta H_{\text{conv}} > \Delta H_{\text{vap}}\), every single mole of ortho-\(\text{H}_2\) that relaxes inside a cryogenic tank boils off more than 1.5 moles of liquid hydrogen! This underscores why industrial rocket propellants (such as liquid hydrogen for Saturn V and SLS) mandate complete catalytic para-conversion during liquefaction.