Rotational Energy Levels - Physical Chemistry 2

SECTION 1

Historical Context & Motivation

The study of how molecules rotate and absorb electromagnetic radiation has deep roots in the development of quantum mechanics itself. In the late nineteenth century, physicists recognized that the heat capacities of diatomic gases could not be explained by classical mechanics alone—rotational degrees of freedom appeared to "freeze out" at low temperatures, a phenomenon that demanded a quantum mechanical explanation. The quantization of rotational energy levels provided one of the earliest and most elegant confirmations of quantum theory applied to molecular systems, bridging the gap between abstract wave mechanics and directly observable spectral lines in the microwave region of the electromagnetic spectrum.

1913 • Bohr's Quantization Postulate

Niels Bohr introduced quantized angular momentum for the hydrogen atom, establishing the principle that angular momentum—and by extension rotational motion—is restricted to discrete values.

1926 • Schrödinger's Wave Equation

Erwin Schrödinger formulated the wave equation, from which the rigid rotor solution was derived, providing exact energy eigenvalues for a rotating diatomic molecule in terms of the quantum number J.

1934 • Cleeton & Williams: First Microwave Absorption

C. E. Cleeton and N. H. Williams observed the first microwave absorption spectrum of ammonia (NH₃), directly confirming the quantization of molecular rotational energy levels.

1950s • Microwave Spectroscopy Matures

Post-war advances in microwave technology enabled precise determination of bond lengths, bond angles, and molecular geometries through high-resolution rotational spectroscopy, transforming structural chemistry.

The central question that rotational spectroscopy addresses is deceptively simple: how does a molecule rotate, and what does its rotational spectrum reveal about its structure? Classical mechanics predicts a continuous distribution of rotational energies, yet experimental microwave spectra show a series of discrete, evenly spaced absorption lines. Resolving this discrepancy requires treating the molecule as a quantum mechanical rigid rotor, a model that yields quantized energy levels and selection rules governing allowed transitions.

SECTION 2

Core Principles & Definitions

Before diving into the mathematics, it is essential to establish several foundational concepts. The treatment of molecular rotation in spectroscopy begins with the idealization of a molecule as a rigid rotor—a system in which the bond length remains fixed during rotation. This model, while approximate, captures the essential physics and serves as the starting point from which more realistic models (centrifugal distortion, non-rigid rotors) are developed. The key parameters are the moment of inertia I, the rotational constant B, the rotational quantum number J, and the selection rules that determine which transitions are spectroscopically active.

Rigid Rotor Model

A diatomic molecule is treated as two point masses separated by a fixed distance r₀. All potential energy is ignored; only kinetic energy of rotation is considered. The Schrödinger equation for this system yields quantized energy levels characterized by the quantum number J.

Moment of Inertia (I)

Defined as I = μr₀², where μ is the reduced mass of the diatomic pair and r₀ is the equilibrium bond length. The moment of inertia governs how "easily" a molecule rotates—larger I means greater resistance to angular acceleration and more closely spaced energy levels.

Rotational Constant (B)

The rotational constant B = ħ²/(2I) (in energy units) or B̃ = h/(8π²Ic) (in wavenumber units, cm⁻¹) encapsulates molecular geometry into a single spectroscopic parameter. Measuring B from a spectrum directly yields the bond length.

Selection Rules

For a pure rotational transition to occur, the molecule must possess a permanent dipole moment (μ ≠ 0), and the quantum number must change by exactly one unit: ΔJ = ±1. Homonuclear diatomics (H₂, N₂, O₂) are therefore rotationally inactive in absorption.

Degeneracy (2J + 1)

Each rotational level J is (2J + 1)-fold degenerate, corresponding to the 2J + 1 allowed values of the magnetic quantum number

M_J

(ranging from −J to +J). This degeneracy affects the Boltzmann population distribution and thus the intensities of spectral lines.

✦

KEY TAKEAWAY

Think of the rigid rotor like a figure skater spinning with arms held at a fixed distance. The skater cannot spin at any arbitrary speed—only at certain quantized angular velocities. A lighter skater with shorter arms (smaller I) spins faster for the same angular momentum, just as a lighter molecule with a shorter bond has a larger rotational constant B and more widely spaced energy levels. The selection rule ΔJ = ±1 is analogous to a constraint that the skater can only speed up or slow down by one "notch" at a time.

SECTION 3

Visual Explanation: Energy Level Diagram

Energy level diagram for a rigid rotor. Horizontal lines represent the allowed energy levels labeled by quantum number J. The energy of each level is E_J = BJ(J + 1), and the degeneracy g = 2J + 1 is shown on the right. Dashed vertical arrows indicate allowed transitions (ΔJ = +1), with transition energies increasing as 2B, 4B, 6B, 8B—producing equally spaced lines in the microwave spectrum.

The diagram above captures the most distinctive feature of rotational energy levels: the spacing between adjacent levels is not constant but increases linearly with J. Specifically, the gap between level J and level J + 1 is 2B(J + 1). However, when we apply the selection rule ΔJ = +1 and look at absorption transitions from level J to J + 1, the transition frequencies ν̃ = 2B̃(J + 1) are equally spaced by 2B̃ in wavenumber units. This elegant pattern—a "picket fence" of evenly spaced lines—is the hallmark of a pure rotational spectrum and allows immediate extraction of the rotational constant from experimental data.

SECTION 4

Mathematical Framework

The quantum mechanical treatment of molecular rotation begins with solving the Schrödinger equation for the rigid rotor Hamiltonian. For a diatomic molecule with reduced mass μ and fixed bond length r₀, the Hamiltonian contains only the kinetic energy of rotation, expressed in terms of the angular momentum operator. The resulting eigenvalue equation yields exact, closed-form solutions—the spherical harmonics Y_J^M(θ, φ)—with corresponding energy eigenvalues that depend on a single quantum number J.

REDUCED MASS

μ = m₁m₂ / (m₁ + m₂)

where m₁ and m₂ are the atomic masses of the two atoms in the diatomic molecule. The reduced mass captures the effective inertia of the two-body system in its center-of-mass frame.

MOMENT OF INERTIA

I = μr₀²

where r₀ is the equilibrium bond length. The moment of inertia I determines the rotational inertia of the molecule and directly connects molecular geometry to spectroscopic observables.

ENERGY EIGENVALUES

E_J = (ħ² / 2I) × J(J + 1) = BJ(J + 1), J = 0, 1, 2, …

Here B = ħ²/(2I) is the rotational constant in energy units (joules). In spectroscopy, it is conventional to express B in wavenumber units: B̃ = h/(8π²Ic), with units of cm⁻¹. The quantum number J takes non-negative integer values, and each level has a degeneracy of (2J + 1) due to the magnetic quantum number M_J.

TRANSITION FREQUENCIES (ΔJ = +1)

ν̃(J → J + 1) = 2B̃(J + 1), J = 0, 1, 2, …

This expression predicts that absorption lines in the pure rotational spectrum appear at ν̃ = 2B̃, 4B̃, 6B̃, … —an evenly spaced series separated by 2B̃. From any two adjacent lines, B̃ can be extracted, and thence the bond length r₀ can be calculated.

ℹ

Derivation Note

The transition frequency is obtained by computing ΔE = E_J+1 − E_J = B[(J+1)(J+2) − J(J+1)] = B[J² + 3J + 2 − J² − J] = 2B(J + 1). Converting to wavenumbers via ν̃ = ΔE/(hc) gives ν̃ = 2B̃(J + 1). This derivation is fundamental to analyzing rotational spectra.

SECTION 5

Rotational Spectrum & Intensity Distribution

While the positions of rotational spectral lines are determined by the energy level spacing, their intensities are governed by the population of each rotational level, which depends on both the Boltzmann distribution and the degeneracy factor (2J + 1). The population of level J relative to J = 0 is given by N_J/N₀ = (2J + 1) × exp[−BJ(J + 1)/(k_BT)]. At low J, the degeneracy factor dominates, causing populations to increase with J. At high J, the exponential Boltzmann factor overwhelms the degeneracy, and populations decline. The result is a population maximum at some intermediate J_max ≈ (k_BT / 2B)^1/2 − 1/2, and the rotational spectrum exhibits a characteristic intensity envelope that rises, reaches a maximum, and then falls off.

Simulated pure rotational absorption spectrum of a diatomic molecule. Lines are equally spaced by 2B̃ along the wavenumber axis. The dashed curve represents the intensity envelope, which peaks near J_max due to the interplay between increasing degeneracy (2J + 1) and decreasing Boltzmann population at higher J values.

The spectrum above illustrates a critical experimental observation: although energy level spacings grow with J, the transition frequencies ν̃ = 2B̃(J + 1) form an arithmetic progression with a common difference of 2B̃. The intensity pattern is not symmetric—it rises steeply at low J due to the rapid increase in degeneracy, reaches a maximum that depends on temperature and the rotational constant, and then decays exponentially. At room temperature, the most populated level for a typical small molecule such as HCl lies around J = 3 to J = 4, meaning the strongest absorption line appears near 8B̃ to 10B̃.

Electromagnetic Spectrum: Rotational Spectroscopy Region

Low energy High energy

Rotational transitions (25%)

SECTION 6

Worked Example: Bond Length of ¹²C¹⁶O from Microwave Spectrum

The pure rotational spectrum of carbon monoxide (¹²C¹⁶O) shows equally spaced absorption lines separated by 3.8626 cm⁻¹. From this single experimental measurement, we can determine the rotational constant, the moment of inertia, and ultimately the C–O bond length.

Determining the Bond Length of CO from Its Rotational Spectrum

1
Step 1 — Extract the Rotational Constant

The spacing between adjacent lines in a rigid rotor spectrum is 2B̃. Given that the line spacing is 3.8626 cm⁻¹, we have 2B̃ = 3.8626 cm⁻¹, so B̃ = 3.8626 / 2 = 1.9313 cm⁻¹.

B̃ = 1.9313 cm⁻¹
2
Step 2 — Calculate the Moment of Inertia

Using B̃ = h / (8π²Ic), we rearrange to find I = h / (8π²cB̃). Substituting h = 6.626 × 10⁻³⁴ J·s, c = 2.998 × 10¹⁰ cm/s, and B̃ = 1.9313 cm⁻¹: I = (6.626 × 10⁻³⁴) / (8 × π² × 2.998 × 10¹⁰ × 1.9313) I = (6.626 × 10⁻³⁴) / (4.5730 × 10¹¹) I = 1.449 × 10⁻⁴⁶ kg·m²

I = 1.449 × 10⁻⁴⁶ kg·m²
3
Step 3 — Compute the Reduced Mass

For ¹²C (m₁ = 12.000 u) and ¹⁶O (m₂ = 15.995 u), the reduced mass is: μ = (12.000 × 15.995) / (12.000 + 15.995) = 191.940 / 27.995 = 6.856 u Converting to kg: μ = 6.856 × 1.6605 × 10⁻²⁷ kg = 1.139 × 10⁻²⁶ kg

μ = 1.139 × 10⁻²⁶ kg
4
Step 4 — Determine the Bond Length

From I = μr₀², we solve for r₀: r₀ = $(I / μ)^{1/2}$ = $(1.449 × 10⁻⁴⁶ / 1.139 × 10⁻²⁶)^{1/2}$ r₀ = $(1.272 × 10⁻²⁰)^{1/2}$ = 1.128 × 10⁻¹⁰ m = 112.8 pm

r₀(C–O) = 112.8 pm
5
Step 5 — Verify Against Literature

The accepted bond length of CO is 112.8 pm, in excellent agreement with our calculated value. This demonstrates the remarkable precision of microwave spectroscopy for determining molecular geometry—bond lengths can be measured to within ±0.1 pm using this technique.

Literature value: 112.8 pm ✓

SECTION 7

Strengths & Limitations of the Rigid Rotor Model

The rigid rotor model provides a remarkably useful first approximation to the rotational behavior of molecules, but it is essential to understand where it succeeds and where it falls short. Real molecules are not perfectly rigid—bonds stretch under centrifugal force during rotation, and interactions between vibrational and rotational motion introduce additional corrections. Recognizing these limitations prepares us for the more sophisticated models encountered in advanced spectroscopy courses.

Feature	Strength	Limitation
Equal line spacing	Predicts evenly spaced lines (2B̃ apart), matching low-J experimental data very well.	At high J, centrifugal distortion causes lines to become slightly closer together than 2B̃ predicts.
Bond length determination	Yields accurate equilibrium bond lengths from a single spectroscopic constant B̃.	The measured B̃ is an effective constant averaged over vibrational motion; the true equilibrium $B_e$ requires vibrational corrections.
Applicability	Works for any linear molecule or diatomic, and extends to symmetric/asymmetric tops with additional quantum numbers.	Homonuclear diatomics (H₂, N₂) have no dipole moment and are invisible to pure rotational spectroscopy.
Mathematical simplicity	Closed-form energy expression $E_J$ = BJ(J + 1) makes calculations straightforward and transparent.	Ignores centrifugal distortion (D), vibration-rotation coupling (α_e), and anharmonicity—all of which are measurable in high-resolution spectra.

Comparison of rigid rotor strengths and limitations

✦

KEY TAKEAWAY

The rigid rotor model is to rotational spectroscopy what the ideal gas law is to thermodynamics: a foundational approximation that captures the essential physics and provides quantitative predictions accurate to within a few percent for most molecules. Just as real gases deviate from PV = nRT at high pressures, real molecules deviate from the rigid rotor at high J values due to centrifugal distortion. The correction term −DJ²(J + 1)² is added to produce the non-rigid rotor model, analogous to the van der Waals correction to the ideal gas equation.

SECTION 8

Connection to Advanced Theory: Non-Rigid Rotor & Polyatomic Molecules

The rigid rotor model naturally leads to several important extensions in molecular spectroscopy. Two immediate generalizations are the non-rigid rotor (which accounts for centrifugal distortion) and the treatment of polyatomic molecules (which requires classification into spherical, symmetric, and asymmetric tops based on their principal moments of inertia). Understanding these connections is critical for interpreting real spectroscopic data and for advanced courses in molecular structure determination.

Property	Rigid Rotor	Non-Rigid Rotor
Energy expression	$E_J$ = BJ(J + 1)	$E_J$ = BJ(J + 1) − DJ²(J + 1)²
Line spacing	Constant: 2B̃ between all adjacent lines	Decreases slightly at high J: 2B̃ − 4D̃(J + 1)²
Bond length assumption	Fixed at equilibrium value r₀	Effective bond length increases with J due to centrifugal stretching
Constants fitted	One: B̃	Two: B̃ and D̃ (centrifugal distortion constant)
Typical D̃/B̃ ratio	Not applicable (D = 0)	~10⁻⁴ to 10⁻⁶, confirming D̃ is a small correction

Rigid rotor vs. non-rigid rotor comparison

For polyatomic molecules, the rotational energy level structure becomes considerably richer. Linear polyatomics (e.g., CO₂, HCN) behave like diatomics with a single moment of inertia and quantum number J. Symmetric tops (e.g., CH₃Cl, NH₃) have two distinct moments of inertia and require a second quantum number K to describe the projection of angular momentum onto the molecular symmetry axis. Asymmetric tops (e.g., H₂O) have three distinct moments of inertia and generally lack closed-form energy expressions, requiring numerical diagonalization. The classification scheme based on principal moments—I_A ≤ I_B ≤ I_C—underpins the entire field of rotational spectroscopy of complex molecules.

💡

Looking Ahead

In subsequent topics, you will encounter vibration-rotation coupling, where each vibrational transition is accompanied by rotational fine structure (P, Q, and R branches), and Raman rotational spectroscopy, which circumvents the dipole moment requirement and enables the study of homonuclear diatomics. The rigid rotor energy level framework developed here provides the foundation for all of these extensions.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain why homonuclear diatomic molecules such as N₂ and O₂ do not exhibit pure rotational absorption spectra, even though they clearly rotate. What physical property must a molecule possess for a rotational transition to be spectroscopically active?

Show answer

A pure rotational transition requires the molecule to possess a permanent electric dipole moment (μ ≠ 0). When a polar molecule rotates, its oscillating dipole couples to the electric field of electromagnetic radiation, enabling absorption. Homonuclear diatomics have zero permanent dipole moment due to perfect symmetry—the center of charge coincides with the center of mass—so there is no oscillating dipole to drive the transition. This is a direct consequence of the selection rule for electric dipole transitions. These molecules can, however, be studied via Raman spectroscopy, which depends on the anisotropy of the molecular polarizability rather than the permanent dipole moment.

Question 2: BASIC CALCULATION

The rotational constant of HCl is B̃ = 10.593 cm⁻¹. Calculate the wavenumber of the J = 3 → J = 4 absorption line in the pure rotational spectrum.

Show answer

Using ν̃ = 2B̃(J + 1) for the transition from J to J + 1, with J = 3: ν̃ = 2 × 10.593 × (3 + 1) = 2 × 10.593 × 4 = 84.744 cm⁻¹. This line falls in the far-infrared/microwave region. The equal spacing between successive lines is 2B̃ = 21.186 cm⁻¹, meaning adjacent transitions differ by this constant amount—a signature of the rigid rotor approximation. The J = 0→1 line would appear at 21.186 cm⁻¹, J = 1→2 at 42.372 cm⁻¹, and so on.

Question 3: INTERMEDIATE

The first two absorption lines in the pure rotational spectrum of ¹H¹⁹F appear at 41.11 cm⁻¹ and 82.22 cm⁻¹. (a) Determine the rotational constant B̃. (b) Calculate the bond length of HF. Use m(¹H) = 1.0078 u, m(¹⁹F) = 18.998 u.

Show answer

(a) The first line corresponds to J = 0 → 1 at ν̃ = 2B̃, so B̃ = 41.11 / 2 = 20.555 cm⁻¹. As a check, the second line at 82.22 cm⁻¹ = 4B̃ = 4 × 20.555 = 82.22 cm⁻¹ ✓. (b) From B̃ = h/(8π²Ic), I = h/(8π²cB̃) = (6.626 × 10⁻³⁴)/(8π² × 2.998 × 10¹⁰ × 20.555) = 1.363 × 10⁻⁴⁷ kg·m². The reduced mass is μ = (1.0078 × 18.998)/(1.0078 + 18.998) = 0.9571 u = 1.590 × 10⁻²⁷ kg. Then r₀ = √(I/μ) = √(1.363 × 10⁻⁴⁷ / 1.590 × 10⁻²⁷) = √(8.574 × 10⁻²¹) = 9.26 × 10⁻¹¹ m = 92.6 pm. This agrees well with the accepted HF bond length of 91.7 pm; the small discrepancy arises from the rigid rotor approximation.

Question 4: APPLIED

In radio astronomy, the J = 1 → 0 emission line of ¹²C¹⁶O at 115.271 GHz is one of the most important tracers of molecular gas in the interstellar medium. (a) Convert this frequency to wavenumber (cm⁻¹). (b) Calculate B̃ for CO from this measurement. (c) Predict the frequency of the J = 2 → 1 transition.

Show answer

(a) ν̃ = ν/c = (115.271 × 10⁹ Hz) / (2.998 × 10¹⁰ cm/s) = 3.845 cm⁻¹. (b) The J = 1 → 0 transition has ν̃ = 2B̃(0 + 1) = 2B̃, so B̃ = 3.845 / 2 = 1.9225 cm⁻¹. This is very close to the laboratory value of 1.9313 cm⁻¹; the slight difference reflects that the astronomical measurement samples a specific isotopologue and vibrational state. (c) The J = 2 → 1 transition occurs at ν̃ = 2B̃ × 2 = 4B̃ = 4 × 1.9225 = 7.690 cm⁻¹, corresponding to ν = ν̃ × c = 7.690 × 2.998 × 10¹⁰ = 230.54 GHz. This is indeed the well-known CO(2–1) line at 230.538 GHz observed in millimeter-wave astronomy.

Question 5: CRITICAL THINKING

High-resolution spectroscopy of HCl reveals that the spacing between successive rotational absorption lines is not perfectly constant but decreases very slightly at higher J. The J = 0→1 line appears at 20.878 cm⁻¹ and the J = 9→10 line at 198.504 cm⁻¹. (a) What is the predicted position of the J = 9→10 line using the rigid rotor model with B̃ determined from the J = 0→1 line? (b) What is the discrepancy, and what physical phenomenon accounts for it? (c) Using the non-rigid rotor formula ν̃ = 2B̃(J + 1) − 4D̃(J + 1)³, estimate the centrifugal distortion constant D̃.

Show answer

(a) From the J = 0→1 line: 2B̃ = 20.878 cm⁻¹, so B̃ = 10.439 cm⁻¹. The rigid rotor prediction for J = 9→10 is ν̃ = 2B̃(10) = 20 × 10.439 = 208.78 cm⁻¹. (b) The observed value is 198.504 cm⁻¹, a discrepancy of 208.78 − 198.504 = 10.28 cm⁻¹. This arises from centrifugal distortion: as J increases, the molecule rotates faster, the centrifugal force stretches the bond, increasing I and effectively decreasing B. This lowers the energy levels relative to rigid rotor predictions. (c) Using ν̃ = 2B̃(J+1) − 4D̃(J+1)³ for J = 9: 198.504 = 2(10.439)(10) − 4D̃(10)³ → 198.504 = 208.78 − 4000D̃ → 4000D̃ = 10.276 → D̃ ≈ 2.57 × 10⁻³ cm⁻¹. Note: this is a simplified estimate; in practice, B̃ and D̃ are fitted simultaneously to all observed lines using least-squares regression, giving more precise values (literature D̃ for HCl ≈ 5.3 × 10⁻⁴ cm⁻¹). The discrepancy suggests our single-line B̃ value is itself slightly affected by the approximation.

SUMMARY

Summary

The rigid rotor model treats a diatomic molecule as two point masses separated by a fixed bond length, yielding quantized rotational energy levels given by E_J = BJ(J + 1), where the rotational constant B = ħ²/(2I) encodes the moment of inertia I = μr₀². Each level carries a degeneracy of (2J + 1), and the selection rule ΔJ = ±1 (requiring a permanent dipole moment) produces equally spaced absorption lines separated by 2B̃ in the microwave region of the electromagnetic spectrum.

From a single measured line spacing, one can extract B̃, compute I, and determine the bond length r₀ = (I/μ)^1/2 with picometer precision. The intensity pattern of the spectrum reflects the interplay between degeneracy and the Boltzmann distribution, producing a characteristic rise-and-fall envelope. At high J, centrifugal distortion causes deviations from equal spacing, leading to the non-rigid rotor correction E_J = BJ(J + 1) − DJ²(J + 1)². This framework forms the foundation for all molecular rotational spectroscopy, from simple diatomics to complex polyatomic symmetric and asymmetric top molecules.