Rotational Partition Function

Drop a molecule of carbon monoxide into a warm gas cell and it does not just fly around. It also tumbles end over end, sampling a huge number of rotational states almost continuously. The rotational partition function, $q_{\text{rot}}$ , counts how many such tumbling states are thermally accessible at temperature $T$ . This quantity matters because rotations store energy and entropy long before vibrations or electronic excitations switch on. Rotational contributions dominate the heat capacity of gases near room temperature and set the intensity pattern of microwave and rotational Raman spectra. This article builds $q_{\text{rot}}$ from the rigid-rotor energy levels. You will see the exact sum for a linear molecule. Then you will see why, at ordinary temperatures, that sum collapses into a simple closed-form integral. Finally you will extend the result to nonlinear, asymmetric-top molecules.

Key takeaways

$q_{\text{rot}}$ sums Boltzmann factors over rotational energy levels, each with its correct degeneracy $2J+1$ .
For a linear rotor, the high-temperature limit gives $q_{\text{rot}} = k_BT/(hcB)$ , remarkably simple and exact to good approximation above ~50 K for most molecules.
Nonlinear molecules need three rotational constants $A$ , $B$ , $C$ and give $q_{\text{rot}} \propto T^{3/2}(ABC)^{-1/2}$ .
The high- $T$ approximation replaces a discrete sum with an integral, valid whenever level spacing is much smaller than $k_BT$ .
Small, light molecules like H $_2$ at low $T$ break the approximation and need the exact sum plus symmetry corrections.
Key idea: because rotational energy levels are closely spaced compared with $k_BT$ at ordinary temperatures, the partition function sum becomes an integral, giving $q_{\text{rot}}$ in closed form from the rotational constant alone.

The core insight

Rotational energy levels are so finely spaced relative to $k_BT$ for most molecules at room temperature that the partition function's discrete sum over $J$ can be replaced by a continuous integral, yielding the compact classical-limit formulas $q_{\text{rot}} = k_BT/(hcB)$ for linear rotors and an analogous $T^{3/2}$ law for nonlinear ones.

The physical picture

A rigid diatomic or linear molecule rotates like a dumbbell with a fixed bond length. Quantum mechanics restricts its rotational energy to discrete levels labeled by the quantum number $J = 0, 1, 2, \dots$ . Each level is $(2J+1)$ -fold degenerate, since the rotor can point its angular momentum in $2J+1$ different orientations (values of $M_J$ ) for the same energy. The rotational partition function was already anticipated when the total molecular partition function was factorized into translational, rotational, vibrational, and electronic pieces. Here we build the rotational factor explicitly. Because typical rotational energy spacings, expressed as $hcB$ , are only a few wavenumbers, thousands of levels are populated even at 298 K. That is precisely the condition under which a sum can be approximated by an integral — the same logic used for translational partition functions, but now applied to a discrete ladder of degenerate levels rather than a near-continuum from the start. The key physical trade-off is between two counteracting effects as $J$ increases: the energy $hcBJ(J+1)$ grows, suppressing the Boltzmann factor, while the degeneracy $2J+1$ grows, enhancing the population. Their product peaks at an intermediate $J$ , giving the population distribution its characteristic asymmetric hump — the origin of the intensity envelope seen in rotational spectra.

Figure

Energy level ladder and degeneracy-weighted population for a linear rotor

The second useful picture is how $q_{\text{rot}}$ itself grows with temperature. Since more levels become thermally accessible as $T$ rises, $q_{\text{rot}}$ increases without bound, tracking $T$ linearly for a linear rotor and as $T^{3/2}$ for a nonlinear one.

Figure

Growth of q_rot with temperature for CO, B = 1.931 cm⁻¹

The core relations

q_{\text{rot}} = \sum_{J=0}^{\infty} (2J+1)\, e^{-hcBJ(J+1)/k_BT}

where: $B$ is the rotational constant in wavenumbers, $J$ the rotational quantum number, $k_B$ Boltzmann's constant.

q_{\text{rot}} \approx \frac{k_BT}{hcB} \quad \text{(linear molecule, high-}T\text{ limit)}

where: valid when $hcB \ll k_BT$ , i.e. essentially all molecules above a few tens of kelvin except H $_2$ and its isotopologues.

q_{\text{rot}} \approx \frac{\sqrt{\pi}}{\sigma}\left(\frac{k_BT}{hc}\right)^{3/2}\frac{1}{\sqrt{ABC}} \quad \text{(nonlinear molecule)}

where: $A \ge B \ge C$ are the three rotational constants (cm⁻¹) about the principal axes, $\sigma$ is the symmetry number (previewed here, treated in full next).

E_J = hcBJ(J+1), \qquad B = \frac{h}{8\pi^2 c I}

where: $I$ is the moment of inertia, $I = \mu r^2$ for a diatomic with reduced mass $\mu$ and bond length $r$ .

The core relations

Derivation AFrom the rigid-rotor sum to q_rot = kT/hcB

Goal. Show that the discrete sum over rotational levels reduces, at ordinary temperatures, to the closed-form expression $q_{\text{rot}} = k_BT/(hcB)$ for a linear rigid rotor.

Step 1. Write the exact partition function. The allowed energies of a rigid linear rotor are $E_J = hcBJ(J+1)$ with degeneracy $g_J = 2J+1$ . The partition function is therefore

q_{\text{rot}} = \sum_{J=0}^{\infty}(2J+1)\, e^{-hcBJ(J+1)/k_BT}.

Step 2. Identify the small parameter. Define $\beta = hcB/k_BT$ . For CO, $B = 1.931\ \text{cm}^{-1}$ , and at $T = 298\ \text{K}$ , $k_BT/hc = 207.2\ \text{cm}^{-1}$ , so

\beta = \frac{1.931}{207.2} = 9.32\times10^{-3} \ll 1.

Because successive terms in the sum change only slightly (the exponent changes by roughly $2\beta J$ per step), the sum behaves like a smoothly varying function of $J$ sampled at unit intervals.

Step 3. Replace the sum with an integral. When a function varies slowly over unit steps in its index, the Euler–Maclaurin approximation justifies

\sum_{J=0}^{\infty} f(J) \approx \int_0^{\infty} f(J)\, dJ,

with corrections of order $\beta$ that vanish as $\beta \to 0$ . Applying this to $q_{\text{rot}}$ ,

q_{\text{rot}} \approx \int_0^{\infty} (2J+1)\, e^{-\beta J(J+1)}\, dJ.

Step 4. Perform a substitution to make the integral elementary. Let $x = J(J+1)$ , so $dx = (2J+1)\,dJ$ . The integral becomes exactly

q_{\text{rot}} \approx \int_0^{\infty} e^{-\beta x}\, dx = \frac{1}{\beta}.

Step 5. Substitute back for β. Since $\beta = hcB/k_BT$ ,

q_{\text{rot}} \approx \frac{1}{\beta} = \frac{k_BT}{hcB}.

Step 6. Extend to a nonlinear rotor. A nonlinear molecule has three rotational constants $A$ , $B$ , $C$ (one per principal axis) and energy levels that, for an asymmetric top, have no simple closed form — but the same high- $T$ integral trick applied to the classical rotational partition function (integrating over all three Euler angles and their conjugate momenta) yields

q_{\text{rot}} = \frac{\sqrt{\pi}}{\sigma}\left(\frac{k_BT}{hc}\right)^{3/2}(ABC)^{-1/2}.

The symmetry number $\sigma$ corrects for rotational configurations that are physically indistinguishable — treated fully in the next topic.

Boxed result.

\boxed{q_{\text{rot,linear}} = \frac{k_BT}{hcB}, \qquad q_{\text{rot,nonlinear}} = \frac{\sqrt{\pi}}{\sigma}\left(\frac{k_BT}{hc}\right)^{3/2}\frac{1}{\sqrt{ABC}}}

Worked examples

Example 1 — CO at 298 K (foundational). Carbon monoxide has rotational constant $B = 1.931\ \text{cm}^{-1}$ (spectroscopic value). We want $q_{\text{rot}}$ at 298 K, ignoring symmetry corrections (CO is heteronuclear, so $\sigma = 1$ ). Step 1: compute $k_BT/hc$ , the thermal energy expressed as an equivalent wavenumber.

\frac{k_BT}{hc} = \frac{(1.381\times10^{-23}\ \text{J K}^{-1})(298\ \text{K})}{(6.626\times10^{-34}\ \text{J s})(2.998\times10^{10}\ \text{cm s}^{-1})}

Numerator: $1.381\times10^{-23} \times 298 = 4.116\times10^{-21}\ \text{J}$ . Denominator: $6.626\times10^{-34}\times 2.998\times10^{10} = 1.986\times10^{-23}\ \text{J cm}$ .

\frac{k_BT}{hc} = \frac{4.116\times10^{-21}}{1.986\times10^{-23}} = 207.3\ \text{cm}^{-1}.

Step 2: divide by $B$ .

q_{\text{rot}} = \frac{207.3}{1.931} = 107.3.

So about 107 rotational states are effectively occupied per molecule at room temperature — consistent with the dense rotational fine structure seen in the CO vibration-rotation (IR) spectrum, where dozens of resolvable $J$ lines appear on each side of the band center. Sanity check: since $\beta = 1/107.3 = 9.3\times10^{-3} \ll 1$ , the high- $T$ approximation is excellent, matching Step 2 of the derivation.

Example 2 — Water vapor at 373 K (applied, nonlinear rotor). Water is a nonlinear (bent) asymmetric top with rotational constants (from microwave spectroscopy):

A = 27.88\ \text{cm}^{-1}, \quad B = 14.52\ \text{cm}^{-1}, \quad C = 9.285\ \text{cm}^{-1}.

Water's symmetry number is $\sigma = 2$ (a $C_2$ axis swaps the two equivalent hydrogens). We want $q_{\text{rot}}$ at the normal boiling point, 373 K. Step 1: compute $k_BT/hc$ at 373 K.

\frac{k_BT}{hc} = \frac{(1.381\times10^{-23})(373)}{1.986\times10^{-23}} = \frac{5.152\times10^{-21}}{1.986\times10^{-23}} = 259.4\ \text{cm}^{-1}.

Step 2: raise to the 3/2 power.

(259.4)^{3/2} = 259.4 \times \sqrt{259.4} = 259.4 \times 16.11 = 4178.

Step 3: compute $ABC$ and its square root.

ABC = 27.88 \times 14.52 \times 9.285 = 27.88 \times 134.85 = 3759.

\sqrt{ABC} = \sqrt{3759} = 61.31.

Step 4: assemble, with $\sqrt{\pi} = 1.7725$ .

q_{\text{rot}} = \frac{1.7725}{2}\times\frac{4178}{61.31} = 0.8862 \times 68.15 = 60.4.

So water vapor at 373 K has roughly 60 thermally accessible rotational states per molecule — smaller than a typical diatomic despite three rotational degrees of freedom, because water's small moment of inertia (light atoms, short bonds) gives large rotational constants and hence widely spaced levels. Sanity check: since $A$ , $B$ , $C \sim 10$ – $28\ \text{cm}^{-1}$ are still much smaller than $k_BT/hc \approx 259\ \text{cm}^{-1}$ , the classical high- $T$ formula remains valid, though water is the borderline case among common nonlinear molecules — the approximation would fail badly at, say, 20 K.

Limiting cases & common traps

Situation	What happens	Consequence
$T \to 0$	$\beta \to \infty$ , high- $T$ integral invalid	Must sum discrete $J$ terms directly; only $J=0$ populated
Very light molecule (H $_2$ , HD)	$B$ large (H $_2$ : $60.85\ \text{cm}^{-1}$ ), $\beta$ not small even at 298 K	High- $T$ formula overestimates $q_{\text{rot}}$ ; exact sum or Euler–Maclaurin correction terms needed
Homonuclear diatomic (N $_2$ , O $_2$ , H $_2$ )	Rotating molecule by 180° regenerates an indistinguishable configuration	Must divide $q_{\text{rot}}$ by symmetry number $\sigma = 2$ (details in the next topic)
Nonlinear vs linear molecule	Three rotational degrees of freedom vs two	$q_{\text{rot}} \propto T^{3/2}$ instead of $T^{1}$
Using $B$ in Hz instead of cm⁻¹	Mixing unit systems	Always confirm whether $B$ is quoted in cm⁻¹ (needs a factor of $hc$ ) or in Hz/MHz (needs a factor of $h$ only)

The chemist behind it

The rigid-rotor model and its use in statistical mechanics owe much to the synthesis of spectroscopy and thermodynamics pioneered in the early twentieth century. Applying quantum mechanics to molecular rotation to predict thermodynamic properties — heat capacities, entropies — from spectroscopic constants alone became one of the great triumphs of the emerging field of statistical thermodynamics, tying laboratory microwave and infrared spectra directly to bulk calorimetric data.

Origins · 1 minRichard C. Tolman and how this was discovered

Richard C. Tolman

1881–1948 · American physical chemist who helped establish the statistical-mechanical foundations connecting quantum energy levels, like the rigid rotor, to measurable thermodynamic functions

Ask the tutor: Want to see how the symmetry number changes q_rot for molecules like benzene or ammonia?

Why does q_rot use (2J+1) as a degeneracy factor?

Each rotational level with quantum number $J$ corresponds to $2J+1$ distinct values of the projection quantum number $M_J$ , which labels the orientation of the angular momentum vector along a chosen axis. All $2J+1$ orientations share exactly the same energy in the absence of an external field, so each contributes equally to the Boltzmann sum. Omitting this degeneracy factor would undercount the true number of accessible microstates and give an incorrect partition function.

How is the rotational constant B related to bond length?

The rotational constant is defined as $B = h/(8\pi^2 cI)$ , where $I$ is the moment of inertia. For a diatomic, $I = \mu r^2$ , with $\mu$ the reduced mass and $r$ the bond length. Because $B$ is measured with high precision from microwave or infrared spectra, this relation is routinely inverted to determine bond lengths to four or five significant figures, one of the most accurate structural tools in chemistry.

Does the high-temperature approximation ever fail badly for common gases?

Yes, most notably for hydrogen and its isotopologues. H $_2$ has $B = 60.85\ \text{cm}^{-1}$ , so even at 298 K, $hcB/k_BT \approx 0.29$ , far from the small-parameter limit assumed in the derivation. Exact summation, or Euler–Maclaurin correction terms beyond the leading $k_BT/hcB$ , is needed for quantitative work on H $_2$ , D $_2$ , and HD near and below room temperature.

Why do homonuclear diatomics need a symmetry number?

Rotating a homonuclear diatomic like N $_2$ by 180° produces a configuration physically identical to the original, since the two nuclei are indistinguishable. Naively summing over all $J$ states double-counts these equivalent orientations, so the partition function must be divided by a symmetry number $\sigma = 2$ . This correction, and its generalization to polyatomic point groups, is developed as the next topic.

How does q_rot feed into rotational heat capacity?

Once $q_{\text{rot}}(T)$ is known, thermodynamic functions follow by differentiation: the internal energy is $U = k_BT^2(\partial \ln q_{\text{rot}}/\partial T)$ . For a linear rotor in the high- $T$ limit, $q_{\text{rot}} \propto T$ , giving $U_{\text{rot}} = RT$ per mole and a constant heat capacity contribution $C_{V,\text{rot}} = R$ , in agreement with equipartition and observed heat capacities of diatomic gases near room temperature.

Rotational Partition Function

The physical picture

The core relations

The core relations

Worked examples

Limiting cases & common traps

The chemist behind it

Before this

Where this leads