Unit 2: Rotational & Microwave Spectroscopy of Diatomic & Polyatomic Rotors

A rigid diatomic rotator consists of two point masses \(m_1\) and \(m_2\) held at fixed bond length \(r_0\). In center-of-mass coordinates, this two-body problem reduces to the motion of a single fictitious particle with reduced mass \(\mu\) at distance \(r_0\) from the origin:

\[ \mu = \frac{m_1 m_2}{m_1 + m_2}, \quad I = \mu r_0^2 \]

where \(I\) is the moment of inertia. Classically, the kinetic energy of rotation with angular velocity \(\omega\) is \(E_{\text{rot}} = \frac{1}{2} I \omega^2 = \frac{J^2}{2I}\), where \(J = I\omega\) is the magnitude of the angular momentum vector.

Schrödinger Equation & Spherical Harmonics

Quantum mechanically, the Hamiltonian for a free rigid rotor is:

\[ \hat{H} = \frac{\hat{J}^2}{2I} = -\frac{\hbar^2}{2I} \left[ \frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial}{\partial\theta}\right) + \frac{1}{\sin^2\theta}\frac{\partial^2}{\partial\phi^2} \right] \]

The eigenfunctions of \(\hat{J}^2\) are the spherical harmonics \(Y_{J, M_J}(\theta, \phi)\), with eigenvalues:

\[ \hat{J}^2 Y_{J, M_J} = \hbar^2 J(J+1) Y_{J, M_J}, \quad \hat{J}_z Y_{J, M_J} = \hbar M_J Y_{J, M_J} \]

where \(J = 0, 1, 2, \dots\) is the rotational quantum number, and \(M_J = -J, -J+1, \dots, +J\) is the magnetic quantum number. The energy eigenvalues are:

\[ E_J = \frac{\hbar^2}{2I} J(J+1) \]

Each rotational energy level \(E_J\) is \((2J+1)\)-fold degenerate in the absence of external fields.

In spectroscopy, rotational energy levels are universally expressed in wavenumber units (\(\text{cm}^{-1}\)) as rotational term values \(F(J) = E_J / (hc)\):

\[ F(J) = B J(J+1) \]

where \(B\) is the rotational constant:

\[ B = \frac{h}{8\pi^2 c I} = \frac{\hbar}{4\pi c \mu r_0^2} \]

Gross & Specific Selection Rules

The interaction between the rotating molecule and microwave radiation requires the dipole matrix element to be non-zero:

1. Gross Selection Rule: The molecule must possess a permanent electric dipole moment: \(\mu_0 \neq 0\). Heteronuclear diatomic molecules (e.g., \(\text{CO}, \text{HCl}, \text{NO}\)) exhibit pure rotational microwave spectra, whereas homonuclear diatomics (e.g., \(\text{N}_2, \text{O}_2, \text{H}_2\)) have zero permanent dipole moment and are completely microwave-inactive.
2. Specific Selection Rule: For transitions within an electric dipole permitted manifold: \[ \Delta J = \pm 1, \quad \Delta M_J = 0, \pm 1 \]

Rotational Transition Frequencies

For an absorption transition from level \(J\) to \(J+1\):

\[ \tilde{\nu}_{J \to J+1} = F(J+1) - F(J) = B(J+1)(J+2) - BJ(J+1) = 2B(J+1) \]

This yields a sequence of equally spaced lines in the microwave spectrum:

\[ \tilde{\nu}_{0\to 1} = 2B, \quad \tilde{\nu}_{1\to 2} = 4B, \quad \tilde{\nu}_{2\to 3} = 6B, \quad \tilde{\nu}_{3\to 4} = 8B, \dots \]

The separation between any two adjacent rotational absorption lines is strictly constant and equal to \(2B\):

\[ \Delta \tilde{\nu} = \tilde{\nu}_{J+1\to J+2} - \tilde{\nu}_{J\to J+1} = 2B \]

Measuring \(\Delta \tilde{\nu}\) allows immediate calculation of \(B\), \(I\), and the equilibrium bond length \(r_0\) with extraordinary sub-picometer precision.

### Advanced Mathematical Formulation: Watson's Reduced Asymmetric Rotor Hamiltonians For a general asymmetric top molecule (\(I_A \neq I_B \neq I_C\)), the centrifugal distortion cannot be described by simple scalar constants \(D_J\). Instead, centrifugal distortion couples angular momentum operators of higher even powers. #### 1. The Quartic Centrifugal Distortion Hamiltonian The unreduced effective rotational Hamiltonian up to fourth order in angular momentum is: \[ \hat{H}_{\text{rot}} = X \hat{J}_x^2 + Y \hat{J}_y^2 + Z \hat{J}_z^2 + \frac{1}{4} \sum_{\alpha,\beta,\gamma,\delta} \tau_{\alpha\beta\gamma\delta} \hat{J}_\alpha \hat{J}_\beta \hat{J}_\gamma \hat{J}_\delta \] where \(\alpha, \beta, \gamma, \delta \in \{x, y, z\}\). Because the components of total angular momentum do not commute (\([\hat{J}_x, \hat{J}_y] = -i\hbar \hat{J}_z\)), there exist linear dependencies among the 81 \(\tau\) coefficients. Symmetrization leaves nine independent quartic centrifugal distortion coefficients. #### 2. Watson's Unitary Contact Transformation To eliminate indeterminacy and produce a hermitian Hamiltonian with minimal parameters that can be unambiguously fitted to observed microwave spectra, J. K. G. Watson applied an operator unitary transformation: \[ \tilde{H} = e^{i \hat{S}} \hat{H} e^{-i \hat{S}} = \hat{H} + i [\hat{S}, \hat{H}] - \frac{1}{2} [\hat{S}, [\hat{S}, \hat{H}]] + \dots \] where \(\hat{S}\) is an odd hermitian operator of third degree in \(\hat{J}_\alpha\). Depending on the molecular symmetry and prolate/oblate character, two standardized reductions are employed: **Watson's \(A\)-Reduction (Asymmetric Reduction):** \[ \hat{H}^{(A)} = \frac{1}{2}(B + C) \hat{\mathbf{J}}^2 + \left[A - \frac{1}{2}(B + C)\right] \hat{J}_z^2 + \frac{1}{2}(B - C) (\hat{J}_x^2 - \hat{J}_y^2) \] \[ - \Delta_J \hat{\mathbf{J}}^4 - \Delta_{JK} \hat{\mathbf{J}}^2 \hat{J}_z^2 - \Delta_K \hat{J}_z^4 - 2 \delta_J \hat{\mathbf{J}}^2 (\hat{J}_x^2 - \hat{J}_y^2) - \delta_K \left[ \hat{J}_z^2 (\hat{J}_x^2 - \hat{J}_y^2) + (\hat{J}_x^2 - \hat{J}_y^2) \hat{J}_z^2 \right] \] where the 5 Watson quartic constants \(\Delta_J, \Delta_{JK}, \Delta_K, \delta_J, \delta_K\) represent: - \(\Delta_J\): Isotropic centrifugal stretching. - \(\Delta_{JK}\): Cross-coupling between total rotation and axial rotation. - \(\Delta_K\): Centrifugal expansion along the principal symmetry axis. - \(\delta_J, \delta_K\): Asymmetry-splitting centrifugal modulations. **Watson's \(S\)-Reduction (Symmetric Reduction):** Preferred for near-spherical or accidental symmetric tops where \(\delta_K\) in the \(A\)-reduction becomes ill-conditioned: \[ \hat{H}^{(S)} = \frac{1}{2}(B + C) \hat{\mathbf{J}}^2 + \left[A - \frac{1}{2}(B + C)\right] \hat{J}_z^2 + \frac{1}{2}(B - C) (\hat{J}_x^2 - \hat{J}_y^2) \] \[ - D_J \hat{\mathbf{J}}^4 - D_{JK} \hat{\mathbf{J}}^2 \hat{J}_z^2 - D_K \hat{J}_z^4 + d_1 \hat{\mathbf{J}}^2 (\hat{J}_+^2 + \hat{J}_-^2) + d_2 (\hat{J}_+^4 + \hat{J}_-^4) \] #### 3. Matrix Representation in the Symmetric Top Basis \(|J, K, M\rangle\) The matrix elements are evaluated analytically: \[ \langle J, K | \hat{\mathbf{J}}^2 | J, K \rangle = J(J+1) \] \[ \langle J, K | \hat{J}_z^2 | J, K \rangle = K^2 \] \[ \langle J, K \pm 2 | (\hat{J}_x^2 - \hat{J}_y^2) | J, K \rangle = \frac{1}{2} \sqrt{ [J(J+1) - K(K \pm 1)] [J(J+1) - (K \pm 1)(K \pm 2)] } \] Diagonalization of this tridiagonal band matrix yields the exact rovibrational rotational energy eigenvalues with sub-kHz precision.

Because isotopic substitution changes the nuclear mass without altering the electronic potential energy surface, the equilibrium internuclear distance \(r_0\) remains identical between isotopologues to a very high approximation.

Consider two isotopologues with reduced masses \(\mu\) and \(\mu'\) (where \(\mu' > \mu\)):

\[ I = \mu r_0^2, \quad I' = \mu' r_0^2 \implies \frac{B}{B'} = \frac{I'}{I} = \frac{\mu'}{\mu} \]

Since \(\mu' > \mu\), \(B' < B\). The heavier isotopologue exhibits smaller rotational line spacing. For carbon monoxide:

• \(^{12}\text{C}^{16}\text{O}\): \(\mu = \frac{12.0000 \times 15.9949}{12.0000 + 15.9949} = 6.8562\text{ u} \implies B \approx 1.92118\text{ cm}^{-1}\)
• \(^{13}\text{C}^{16}\text{O}\): \(\mu' = \frac{13.0034 \times 15.9949}{13.0034 + 15.9949} = 7.1699\text{ u} \implies B' \approx 1.83797\text{ cm}^{-1}\)

Kraitchman's Equations for Polyatomics

For polyatomic molecules with multiple unknown bond lengths and angles, isotopic substitution at atom \(i\) shifts the center of mass. Kraitchman derived exact expressions determining the coordinate \(|z_i|\) of the substituted atom relative to the principal axes of the parent molecule:

\[ |z_i| = \sqrt{\frac{\Delta I_x}{\mu_s}}, \quad \mu_s = \frac{M \Delta m}{M + \Delta m} \]

where \(M\) is the total molecular mass and \(\Delta m\) is the mass increase. By systematically substituting each atomic site, the complete 3D structure is elucidated without assuming bond angles.

In real molecules, the chemical bond is not a rigid rod but an elastic spring with force constant \(k\). As the rotational quantum number \(J\) increases, the centrifugal force \(\mu \omega^2 r\) stretches the bond, increasing the moment of inertia \(I\) and slightly lowering the rotational energy levels below the rigid rotor expectation.

Treating centrifugal distortion with perturbation theory gives the modified rotational term formula:

\[ F(J) = B J(J+1) - D J^2(J+1)^2 \]

where \(D\) is the centrifugal distortion constant. For a harmonic bond with vibrational frequency \(\bar{\omega}\) (\(\text{cm}^{-1}\)), Kratzer's relation connects \(D\) directly to \(B\) and \(\bar{\omega}\):

\[ D = \frac{4 B^3}{\bar{\omega}^2} \]

Transition Frequencies for Non-Rigid Rotors

The absorption transition frequency becomes:

\[ \tilde{\nu}_{J \to J+1} = F(J+1) - F(J) = 2B(J+1) - 4D(J+1)^3 \]

The separation between consecutive lines is no longer strictly constant:

\[ \Delta \tilde{\nu} = \tilde{\nu}_{J+1\to J+2} - \tilde{\nu}_{J\to J+1} = 2B - 4D[ (J+2)^3 - (J+1)^3 ] = 2B - 6D(J+1)^2 - 6D(J+1) - 4D \]

Because \(D \ll B\) (typically \(D / B \sim 10^{-4} - 10^{-5}\)), centrifugal distortion is small at low \(J\) but becomes increasingly prominent at high rotational quantum numbers.

The intensity of a rotational absorption line \(J \to J+1\) is governed by the thermal population \(N_J\) of the initial state \(J\). According to the Maxwell-Boltzmann distribution:

\[ N_J \propto g_J \exp\left(-\frac{E_J}{k_B T}\right) = (2J+1) \exp\left(-\frac{h c B J(J+1)}{k_B T}\right) \]

Two opposing mathematical factors govern \(N_J\) as \(J\) increases:

1. The degeneracy factor \(g_J = (2J+1)\) increases linearly with \(J\), promoting population in higher states.
2. The Boltzmann exponential factor \(\exp\left(-\frac{hcBJ(J+1)}{k_B T}\right)\) decreases exponentially with \(J(J+1)\).

Finding the Most Populated State (\(J_{\max}\))

To find the value of \(J\) corresponding to maximum population, we treat \(J\) as a continuous variable and set \(\frac{d N_J}{d J} = 0\):

\[ \frac{d}{dJ}\left[ (2J+1) e^{-\frac{hcB J(J+1)}{k_B T}} \right] = 2 e^{-\frac{hcB J(J+1)}{k_B T}} - (2J+1) \left(\frac{hcB (2J+1)}{k_B T}\right) e^{-\frac{hcB J(J+1)}{k_B T}} = 0 \]

Factoring out the non-zero exponential term:

\[ 2 - \frac{hcB}{k_B T} (2J+1)^2 = 0 \implies (2J+1)^2 = \frac{2 k_B T}{hcB} \] \[ 2J+1 = \sqrt{\frac{2 k_B T}{hcB}} \implies J_{\max} = \sqrt{\frac{k_B T}{2 hcB}} - \frac{1}{2} \]

This explains why microwave spectra display a characteristic intensity envelope: spectral line intensities initially rise with \(J\), reach a smooth maximum at \(J \approx J_{\max}\), and then decay asymptotically to zero at large \(J\).

Polyatomic molecules possess three principal moments of inertia \(I_A \le I_B \le I_C\) along mutually orthogonal principal axes of inertia \(a, b, c\). Rotational constants are defined as:

\[ A = \frac{h}{8\pi^2 c I_A}, \quad B = \frac{h}{8\pi^2 c I_B}, \quad C = \frac{h}{8\pi^2 c I_C} \quad (A \ge B \ge C) \]

Taxonomy of Molecular Rotors

Symmetric Top Quantum Term Formula

For a prolate symmetric top, the projection of total angular momentum \(\vec{J}\) onto the principal molecular symmetry axis is quantized with quantum number \(K = -J, \dots, +J\):

\[ F(J, K) = B J(J+1) + (A - B) K^2 \]

The electric dipole selection rules for a prolate top with dipole along the symmetry axis are \(\Delta J = \pm 1, \Delta K = 0\). Consequently, the transition frequencies are \(\tilde{\nu} = 2B(J+1)\), identical to a linear rotor, meaning \(K\)-splittings are absent in rigid symmetric tops unless centrifugal distortion or the Stark effect is introduced.

When an external static electric field \(\vec{\mathcal{E}}\) is applied to a rotating polar molecule, the interaction between the field and the permanent dipole moment \(\vec{\mu}_0\) lifts the \((2J+1)\)-fold \(M_J\) degeneracy. This field-induced splitting is the Stark effect.

Linear Rotors: Second-Order Stark Shift

For a linear rotor in state \(J\), the first-order Stark energy perturbation vanishes by parity: \(\langle J, M_J | \hat{\vec{\mu}} \cdot \vec{\mathcal{E}} | J, M_J \rangle = 0\). The energy shifts arise from second-order perturbation theory:

\[ \Delta E_{\text{Stark}}^{(2)}(J, M_J) = \frac{\mu_0^2 \mathcal{E}^2}{2hcB} \left[ \frac{J(J+1) - 3 M_J^2}{J(J+1)(2J-1)(2J+3)} \right] \quad (J > 0) \]

For \(J = 0\), the second-order shift is \(\Delta E^{(2)}(0, 0) = -\frac{\mu_0^2 \mathcal{E}^2}{6 hc B}\). Notice that the shift depends on \(M_J^2\), splitting each \(J\) level into \(J+1\) distinct components (\(|M_J| = 0, 1, \dots, J\)).

Symmetric Tops: First-Order Stark Shift

In symmetric tops with \(K \neq 0\), the molecule possesses an average projection of dipole moment along the space-fixed \(Z\)-axis, producing a dramatic first-order Stark effect linear in electric field:

\[ \Delta E_{\text{Stark}}^{(1)}(J, K, M_J) = -\frac{\mu_0 \mathcal{E} K M_J}{J(J+1)} \]

Because the shift is proportional to \(\mu_0 \mathcal{E}\), measuring the frequency separation between Stark components as a function of calibrated field \(\mathcal{E}\) provides the most precise experimental method for determining permanent molecular electric dipole moments in the gas phase.

### Advanced Research Monograph: Chirped-Pulse Fourier Transform Microwave (CP-FTMW) Astrochemistry Historically, microwave spectroscopy was plagued by slow acquisition rates because cavity Fourier-transform spectrometers (Balle-Flygare design) could only sample bandwidths of \(\sim 1\text{ MHz}\) per cavity resonance tuning. #### 1. Broadband Chirped-Pulse Revolution Developed by B. H. Pate and coworkers, Chirped-Pulse Fourier Transform Microwave (CP-FTMW) spectroscopy utilizes high-speed arbitrary waveform generators (AWGs) to sweep microwave frequencies across an 11 GHz bandwidth (e.g., 7 to 18 GHz) in a single microsecond pulse (\(\Delta t = 1.0\text{ }\mu\text{s}\)): \[ E_{\text{chirp}}(t) = E_0 \cos\left( \omega_0 t + \frac{1}{2} \alpha t^2 \right) \] where \(\alpha = \frac{2\pi \Delta \nu}{\Delta t}\) is the linear chirp rate. #### 2. Macroscopic Polarization and Free Induction Decay (FID) - The broadband pulse simultaneously polarizes all electric-dipole allowed rotational transitions across thousands of molecular species in a supersonic expansion jet. - As the molecular ensemble macroscopically dephases, it radiates an 11-GHz broad Free Induction Decay (FID) lasting tens of microseconds. - Real-time digitizers sampling at \(40\text{ GS/s}\) capture the complete FID, which yields the entire high-resolution microwave spectrum upon digital FFT with sub-kHz precision. #### 3. Interstellar Astrochemistry and Prebiotic Molecule Discovery Because pure rotational transitions are the unique fingerprint of gas-phase polar molecules in cold interstellar space (\(T \sim 10 - 50\text{ K}\)): - CP-FTMW laboratory data directly calibrate ALMA (Atacama Large Millimeter/submillimeter Array) radio telescope observations. - This has led to the definitive identification of complex organic molecules (COMs) in the Taurus Molecular Cloud (TMC-1) and Sagittarius B2, including chiral propylene oxide, benzonitrile (the first aromatic ring detected by radio astronomy), and cyano-polyynes.

Microwave spectroscopy encompasses the frequency range from \(3\text{ GHz}\) to \(300\text{ GHz}\) (\(0.1 - 10\text{ cm}^{-1}\)). Because coaxial cables suffer unacceptable attenuation at these frequencies, microwave radiation is guided through hollow rectangular metal waveguides and resonant cavity structures.

Microwave Radiation Sources & Stark Cells

1. Solid-State Generators: Reflex klystrons and backward wave oscillators (BWOs) have been largely superseded by solid-state Gunn diodes and phase-locked microwave frequency synthesizers, providing narrow spectral purity (\(\Delta\nu / \nu < 10^{-8}\)).
2. The Stark Absorption Cell: A long rectangular waveguide (typically \(1 - 3\text{ meters}\)) containing a central isolated metal septum plate. A high-voltage square wave (\(0 - 2\text{ kV}\) at \(100\text{ kHz}\)) is applied to modulate the Stark splitting, enabling phase-sensitive lock-in detection that eliminates low-frequency source drift.
3. Detectors: Low-noise point-contact silicon-tungsten crystal diodes, Schottky barrier diodes, or cryogenic bolometers.

Dielectric Relaxation & Microwave Heating Physics

In condensed phases (liquids and solids), microwave radiation interacts with molecular dipoles not through discrete rotational transitions (which are quenched by collisions), but via collective dielectric relaxation. The response of the medium is characterized by the complex relative permittivity:

\[ \varepsilon^*(\omega) = \varepsilon'(\omega) - i \varepsilon''(\omega) \]

where \(\varepsilon'\) is the real dielectric constant (characterizing capacitive energy storage) and \(\varepsilon''\) is the imaginary dielectric loss factor (characterizing irreversible thermal dissipation).

The Debye Relaxation Formulation

Peter Debye modeled the rotational relaxation of molecular dipoles in a viscous solvent with rotational correlation time \(\tau_D\):

\[ \varepsilon'(\omega) = \varepsilon_\infty + \frac{\varepsilon_s - \varepsilon_\infty}{1 + \omega^2 \tau_D^2}, \quad \varepsilon''(\omega) = \frac{(\varepsilon_s - \varepsilon_\infty) \omega \tau_D}{1 + \omega^2 \tau_D^2} \]

where \(\varepsilon_s\) is the static (low-frequency) permittivity and \(\varepsilon_\infty\) is the high-frequency optical permittivity (\(\approx n^2\)). Dielectric loss \(\varepsilon''\) reaches its absolute maximum when \(\omega \tau_D = 1\), i.e., when the electric field frequency matches the reciprocal molecular reorientation time:

\[ f_{\max} = \frac{1}{2\pi \tau_D} \]

For liquid water at \(25^\circ\text{C}\), \(\tau_D \approx 8.3\text{ ps}\), yielding \(f_{\max} \approx 19\text{ GHz}\). Domestic microwave ovens operate at \(2.45\text{ GHz}\) (\(\lambda = 12.24\text{ cm}\)), where \(\varepsilon''\) is moderate, ensuring uniform penetration depth rather than superficial surface heating.