Part VI β€” Spectroscopy & Molecular Structure

Rotational Spectroscopy

1. Introduction: Microwave Spectroscopy and Molecular Rotation

Rotational spectroscopy probes the quantized angular momentum of molecules by observing absorption and emission in the microwave and far-infrared regions of the electromagnetic spectrum, typically at frequencies from 1 GHz to 1 THz (wavenumbers ~0.03–30 cm$^{-1}$).

A molecule can absorb microwave radiation only if it possesses a permanent electric dipole moment. The oscillating electric field of the radiation couples to the rotating dipole, driving transitions between quantized rotational energy levels. This is why homonuclear diatomics such as H$_2$, N$_2$, and O$_2$ are microwave-inactive β€” they lack a permanent dipole.

From the precisely measured frequencies of rotational transitions, one can extract bond lengths, bond angles, and force constants with extraordinary precision. The technique is among the most accurate methods for determining molecular geometry, with bond lengths routinely measured to $\pm 0.001$ Γ…. Rotational spectroscopy also plays a pivotal role in radio astronomy, where microwave emission lines identify molecules in the interstellar medium.

Key Concepts at a Glance

  • Rigid rotor model β€” the simplest quantum description of a rotating diatomic
  • Selection rules β€” $\Delta J = \pm 1$ for a polar molecule
  • Equally spaced lines β€” spacing of $2B$ in wavenumber
  • Centrifugal distortion β€” corrections to the rigid rotor for bond stretching
  • Population effects β€” Boltzmann distribution governs line intensities

2. Derivation 1: Rigid Rotor Energy Levels

Consider a diatomic molecule with atoms of masses $m_1$ and $m_2$separated by a fixed bond length $r_0$. In the center-of-mass frame, the two-body rotation reduces to a single particle of reduced mass:

$$\mu = \frac{m_1 m_2}{m_1 + m_2}$$

moving on the surface of a sphere of radius $r_0$. The moment of inertia is:

$$I = \mu r_0^2$$

The SchrΓΆdinger Equation for the Rigid Rotor

Since the bond length is fixed, the potential energy is zero and the Hamiltonian contains only kinetic (rotational) energy. In spherical coordinates, the rigid-rotor Hamiltonian is:

$$\hat{H} = \frac{\hat{L}^2}{2I}$$

where $\hat{L}^2$ is the square of the orbital angular momentum operator. The eigenfunctions of $\hat{L}^2$ are the spherical harmonics $Y_J^{M_J}(\theta, \phi)$, satisfying:

$$\hat{L}^2 \, Y_J^{M_J}(\theta, \phi) = \hbar^2 J(J+1) \, Y_J^{M_J}(\theta, \phi)$$

where $J = 0, 1, 2, \ldots$ is the rotational quantum number and$M_J = -J, -J{+}1, \ldots, J{-}1, J$. Substituting into the eigenvalue equation$\hat{H}\psi = E\psi$:

$$\frac{\hat{L}^2}{2I}\,Y_J^{M_J} = \frac{\hbar^2 J(J+1)}{2I}\,Y_J^{M_J}$$

Thus the energy eigenvalues are:

$$\boxed{E_J = \frac{\hbar^2}{2I}\,J(J+1) = B\,J(J+1)}$$

where we define the rotational constant:

$$B = \frac{\hbar^2}{2I} \quad \text{(in energy units)}, \qquad \tilde{B} = \frac{h}{8\pi^2 I c} \quad \text{(in cm}^{-1}\text{)}$$

Degeneracy

The energy $E_J$ depends only on $J$, not on $M_J$. Since $M_J$ runs from $-J$ to $+J$, there are$2J + 1$ degenerate states for each $J$:

$$g_J = 2J + 1$$

This degeneracy reflects the fact that the energy of a freely rotating molecule is independent of the orientation of its angular momentum in space.

Energy Level Spacing

The energy gap between adjacent levels grows linearly with $J$:

$$E_{J+1} - E_J = B\bigl[(J{+}1)(J{+}2) - J(J{+}1)\bigr] = 2B(J+1)$$

The levels are not equally spaced β€” the gap increases as $J$increases β€” but we shall see that the transition frequencies are nearly equally spaced under the selection rule $\Delta J = +1$.

3. Derivation 2: Selection Rules and Transition Frequencies

Transition Dipole Matrix Element

A transition between rotational states $|J, M_J\rangle$ and$|J', M_J'\rangle$ is driven by the interaction of the radiation field with the molecular dipole moment $\hat{\boldsymbol{\mu}}$. The transition rate is proportional to the square of the transition dipole matrix element:

$$\boldsymbol{\mu}_{fi} = \langle J', M_J' | \hat{\boldsymbol{\mu}} | J, M_J \rangle$$

For a diatomic molecule with permanent dipole moment $\mu_0$ along the internuclear axis, the dipole operator in the space-fixed frame involves direction cosines. The $z$-component of the dipole is$\hat{\mu}_z = \mu_0 \cos\theta$, and similarly for $x$ and $y$. Using the relation $\cos\theta \propto Y_1^0(\theta,\phi)$, the integral becomes:

$$\mu_z^{fi} = \mu_0 \int Y_{J'}^{M_J'*}(\theta,\phi)\;\cos\theta\;Y_J^{M_J}(\theta,\phi)\;\sin\theta\,d\theta\,d\phi$$

This is a standard integral over three spherical harmonics (a Gaunt integral). The properties of the Clebsch–Gordan coefficients dictate that it is nonzeroonly when:

$$\boxed{\Delta J = J' - J = \pm 1, \qquad \Delta M_J = 0, \pm 1}$$

and the molecule must have $\mu_0 \neq 0$ (permanent dipole moment). This is the gross selection rule: a polar molecule is required.

Absorbed Frequencies

For an absorption transition $J \to J + 1$, the photon energy must match the energy gap:

$$h\nu = E_{J+1} - E_J = B(J{+}1)(J{+}2) - BJ(J{+}1) = 2B(J+1)$$

In wavenumber units ($\tilde{\nu} = \nu / c$):

$$\boxed{\tilde{\nu}_J = 2\tilde{B}(J+1), \quad J = 0, 1, 2, \ldots}$$

Uniform Spacing of Lines

The first few transitions produce lines at:

  • $J=0 \to 1$: $\tilde{\nu} = 2\tilde{B}$
  • $J=1 \to 2$: $\tilde{\nu} = 4\tilde{B}$
  • $J=2 \to 3$: $\tilde{\nu} = 6\tilde{B}$
  • $J=3 \to 4$: $\tilde{\nu} = 8\tilde{B}$

The spacing between consecutive absorption lines is constant:

$$\Delta\tilde{\nu} = \tilde{\nu}_{J+1} - \tilde{\nu}_J = 2\tilde{B}(J{+}2) - 2\tilde{B}(J{+}1) = 2\tilde{B}$$

This equal spacing of $2\tilde{B}$ is the hallmark of a pure rotational spectrum in the rigid-rotor approximation.

Why Homonuclear Diatomics Have No Rotational Spectrum

Homonuclear diatomic molecules (H$_2$, N$_2$, O$_2$, etc.) possess a center of symmetry. By symmetry, their permanent electric dipole moment vanishes:$\mu_0 = 0$. Since the transition dipole is proportional to $\mu_0$,every rotational matrix element is zero. These molecules are therefore microwave inactive and do not exhibit a pure rotational absorption spectrum.

Note that homonuclear diatomics can still have rotational Raman spectra, as the Raman selection rule depends on the polarizability rather than the dipole moment.

4. Derivation 3: Non-Rigid Rotor (Centrifugal Distortion)

Real molecules are not perfectly rigid. As the molecule rotates faster (higher $J$), the centrifugal force stretches the bond, increasing the moment of inertia and lowering the rotational energy compared to the rigid-rotor prediction.

Classical Picture: Centrifugal Stretching

For a diatomic rotating with angular momentum $L$, the centrifugal force on the reduced mass is $F_{\text{cent}} = \mu\omega^2 r = L^2/(\mu r^3)$. In equilibrium, this is balanced by the restoring force of the bond, modeled as a spring with force constant $k$:

$$k(r - r_0) = \frac{L^2}{\mu r^3}$$

For small displacements $\delta r = r - r_0 \ll r_0$, we approximate$r \approx r_0$ in the centrifugal term:

$$\delta r \approx \frac{L^2}{\mu k r_0^3}$$

Corrected Energy Levels

Treating the centrifugal effect as a perturbation and replacing $L^2$ by its quantum eigenvalue $\hbar^2 J(J+1)$, one obtains the first-order corrected energy:

$$\boxed{E_J = BJ(J+1) - DJ^2(J+1)^2}$$

where $D$ is the centrifugal distortion constant.

Derivation of D

The stretching lowers the rotational energy by an amount proportional to $[J(J+1)]^2$. The detailed calculation yields:

$$D = \frac{4B^3}{\omega^2}$$

where $\omega$ is the harmonic vibrational frequency of the bond (in the same units as $B$). Since vibrational frequencies are much larger than rotational constants ($\omega \gg B$), we have $D \ll B$, typically$D/B \sim 10^{-4}$ to $10^{-6}$.

Effect on Spectral Lines

The transition frequencies for $J \to J + 1$ now become:

$$\tilde{\nu}_J = 2\tilde{B}(J+1) - 4\tilde{D}(J+1)^3$$

The centrifugal correction causes each line to shift slightly to lower frequency. The spacing between consecutive lines decreases:

$$\Delta\tilde{\nu} = 2\tilde{B} - 4\tilde{D}\bigl[3(J+1)^2 + 3(J+1) + 1\bigr]$$

This gradual decrease in spacing is clearly observable in high-resolution spectra. By fitting the measured frequencies to the non-rigid rotor formula, one extracts both $\tilde{B}$and $\tilde{D}$, giving information about both the bond length and the force constant.

5. Derivation 4: Symmetric Top Molecules

A symmetric top molecule has two equal principal moments of inertia. For a prolate symmetric top (e.g., CH$_3$Cl),$I_A < I_B = I_C$, while for an oblate top (e.g., benzene),$I_A = I_B < I_C$. We focus on the prolate case, which includes most common examples such as NH$_3$ and CH$_3$F.

Energy Levels

The rotational Hamiltonian for a symmetric top is:

$$\hat{H} = \frac{\hat{J}^2}{2I_B} + \left(\frac{1}{2I_A} - \frac{1}{2I_B}\right)\hat{J}_a^2$$

where $\hat{J}^2$ is the total angular momentum squared and$\hat{J}_a$ is its component along the molecular symmetry axis. The eigenvalues of $\hat{J}^2$ are $\hbar^2 J(J+1)$ and those of$\hat{J}_a$ are $\hbar K$ with$K = 0, \pm 1, \pm 2, \ldots, \pm J$. Defining rotational constants$A = \hbar^2/(2I_A)$ and $B = \hbar^2/(2I_B)$:

$$\boxed{E_{J,K} = BJ(J+1) + (A - B)K^2}$$

For a prolate top, $A > B$, so $(A - B) > 0$ and higher$|K|$ increases the energy. For an oblate top, replace $A$ with$C$: $E_{J,K} = BJ(J+1) + (C - B)K^2$ where$C - B < 0$.

Selection Rules

For a symmetric top with a dipole moment along the symmetry axis:

$$\Delta J = \pm 1, \qquad \Delta K = 0$$

The $\Delta K = 0$ rule means that the component of angular momentum about the symmetry axis does not change. Since the energy depends on $K^2$, and$K$ does not change, the $K$-dependent term cancels in the transition frequency. The result is the same expression as for a linear rotor:

$$\tilde{\nu} = 2\tilde{B}(J+1)$$

Thus, the pure rotational spectrum of a symmetric top looks like that of a diatomic molecule β€” equally spaced lines separated by $2\tilde{B}$ β€” but each line may carry different degeneracy and intensity patterns.

Application: NH$_3$ and CH$_3$Cl

Ammonia (NH$_3$): A prolate symmetric top with $\tilde{B} \approx 9.94$ cm$^{-1}$ and$\tilde{A} \approx 6.30$ cm$^{-1}$. Its famous inversion doubling at 23.87 GHz (the transition used in the first maser) overlays the rotational structure.

Methyl chloride (CH$_3$Cl): With$\tilde{B} \approx 0.443$ cm$^{-1}$, its rotational lines lie in the far infrared. Centrifugal distortion terms $D_J$, $D_{JK}$, and $D_K$ are needed for precise fitting.

6. Derivation 5: Rotational Temperature and Population Distribution

At thermal equilibrium, the population of rotational level $J$ is given by the Boltzmann distribution:

$$\boxed{\frac{N_J}{N} = \frac{(2J+1)\,g_{\text{ns}}\,\exp\!\bigl[-BJ(J+1)/(k_BT)\bigr]}{q_{\text{rot}}}}$$

where $g_{\text{ns}}$ is the nuclear spin statistical weight (important for molecules with equivalent nuclei), and the rotational partition function is:

$$q_{\text{rot}} = \sum_{J=0}^{\infty}(2J+1)\,g_{\text{ns}}\,\exp\!\bigl[-BJ(J+1)/(k_BT)\bigr]$$

In the high-temperature limit ($k_BT \gg B$), this sum can be replaced by an integral, giving $q_{\text{rot}} \approx k_BT/(\sigma B)$ where$\sigma$ is the symmetry number.

Most Populated Level

The population $N_J$ first increases with $J$ (due to the degeneracy factor $2J+1$) then decreases (due to the Boltzmann exponential). To find the maximum, differentiate $f(J) = (2J+1)\exp[-BJ(J+1)/(k_BT)]$ with respect to $J$ and set it to zero. Treating $J$ as continuous:

$$\frac{df}{dJ} = 2\,e^{-BJ(J+1)/k_BT} + (2J+1)\left(-\frac{B(2J+1)}{k_BT}\right)e^{-BJ(J+1)/k_BT} = 0$$

Simplifying:

$$2 = \frac{B(2J+1)^2}{k_BT}$$

Solving for $J$:

$$\boxed{J_{\max} = \sqrt{\frac{k_BT}{2B}} - \frac{1}{2}}$$

This result has a simple physical interpretation: higher temperature populates higher angular momentum states. For HCl at 300 K with $\tilde{B} \approx 10.6$ cm$^{-1}$,$J_{\max} \approx 3$.

Intensity Envelope

The absorption intensity of the $J \to J+1$ transition is proportional to the population of the lower state $N_J$ and to the transition probability, which scales as $(J+1)$ from the squared matrix element. Thus:

$$I_{J \to J+1} \propto (J+1)(2J+1)\,\exp\!\bigl[-BJ(J+1)/(k_BT)\bigr]$$

The resulting intensity envelope rises steeply at low $J$, peaks near$J_{\max}$, and then decays exponentially. At higher temperatures, the envelope broadens and the peak shifts to higher $J$, distributing intensity over more lines.

7. Applications of Rotational Spectroscopy

Bond Length Determination

Measuring the line spacing $2\tilde{B}$ directly yields the rotational constant $\tilde{B}$. Since $\tilde{B} = h/(8\pi^2 I c)$ and$I = \mu r_0^2$, knowing the atomic masses gives the bond length:

$r_0 = \sqrt{\frac{h}{8\pi^2 \mu c \tilde{B}}}$

For CO, $\tilde{B} = 1.9225$ cm$^{-1}$ yields$r_0 = 1.1283$ Γ….

Isotope Effects

Isotopic substitution changes the reduced mass but not the bond length (to first approximation). Comparing $\tilde{B}$ for two isotopologues:

$\frac{\tilde{B}'}{\tilde{B}} = \frac{\mu}{\mu'}$

For H$^{35}$Cl vs H$^{37}$Cl, the slightly different reduced masses produce doublets in the spectrum that are clearly resolved.

Interstellar Molecule Detection

Radio telescopes detect rotational emission lines from molecules in interstellar clouds. Over 200 molecular species have been identified via their microwave signatures, including CO (the most commonly used tracer of molecular gas), HCN, H$_2$O, NH$_3$, and complex organics. The $J = 1 \to 0$ line of CO at 115.27 GHz is the workhorse of molecular astrophysics.

Atmospheric Monitoring

Microwave and millimeter-wave spectroscopy monitors atmospheric trace gases. Rotational transitions of O$_3$, H$_2$O, ClO, and other species are measured from ground and satellite platforms to track ozone depletion, water vapor profiles, and pollutant distributions with high sensitivity.

8. Historical Context

Cleeton and Williams (1934)

C.E. Cleeton and N.H. Williams at the University of Michigan performed the first microwave spectroscopy experiment, observing the absorption of centimeter-wavelength radiation by ammonia gas. They measured the NH$_3$ inversion line at approximately 24 GHz (1.25 cm wavelength), demonstrating that molecular transitions could be probed in the microwave region. Their work laid the foundation for all of microwave spectroscopy and molecular physics.

Charles H. Townes

Townes exploited the ammonia inversion transition to build the first maser (Microwave Amplification by Stimulated Emission of Radiation) in 1954. His deep understanding of rotational spectroscopy was essential to selecting a suitable molecular transition for population inversion. Townes shared the 1964 Nobel Prize in Physics for this work, which also paved the way for the laser.

Walter Gordy

Gordy, working at Duke University from the 1940s onward, systematically extended microwave spectroscopy to hundreds of molecules. His group pioneered precision measurements of rotational constants, centrifugal distortion parameters, and nuclear quadrupole coupling constants. His textbook Microwave Molecular Spectra (with R.L. Cook) remains a standard reference in the field.

The development of microwave spectroscopy was greatly accelerated by radar technology from World War II, which produced high-power microwave sources (klystrons, magnetrons) and sensitive detection systems that were subsequently adapted for laboratory spectroscopy.

9. Python Simulation: Microwave Spectrum of HCl

This simulation computes and plots the rotational absorption spectrum of HCl at various temperatures, showing how the Boltzmann population distribution shapes the intensity envelope. The four panels display: (1) population $N_J/N$ vs $J$, (2) broadened spectra at different temperatures, (3) a stick spectrum with labeled transitions at 300 K, and (4) the most populated level $J_{\max}$ as a function of temperature.

Rotational Spectrum of HCl: Population & Intensity Analysis

Python

Simulates microwave absorption spectrum of HCl with rigid rotor + centrifugal distortion at multiple temperatures

script.py180 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

10. Fortran Simulation: Rotational Constants and Line Positions

High-precision computation of rotational constants, centrifugal distortion parameters, spectral line positions, and population distributions for a series of diatomic molecules (HCl, HBr, CO, HF, HI). The program also demonstrates bond length determination from $\tilde{B}$ and the isotope effect for H$^{35}$Cl vs H$^{37}$Cl.

Diatomic Rotational Spectroscopy Calculator

Fortran

Computes rotational constants, spectral lines, populations, and isotope effects for diatomic molecules

rotational_spectroscopy.f90140 lines

Click Run to execute the Fortran code

Code will be compiled with gfortran and executed on the server