Chapter 3 • Theory & Derivations

Unit 3: Exact Quantum Models II: The Harmonic Oscillator & Rigid Rotor

Rigorous quantum mechanical treatment of molecular vibrations and rotations: the quantum simple harmonic oscillator, power series solutions, Hermite polynomials, Gaussian ground states, ladder operator algebra (creation and annihilation operators), vibrational zero-point energy, quantum tunneling at classical turning points, the 3D rigid rotor, spherical harmonics, angular momentum operators, and space quantization.

§3.1 The Quantum Harmonic Oscillator: Power Series & Dimensionless Formulation

The harmonic oscillator is the quintessential model for molecular vibrations, describing parabolic potential wells about equilibrium geometry.

The Schrödinger Equation for a Harmonic Oscillator

A particle of mass \(m\) moving in a one-dimensional parabolic potential \(V(x) = \frac{1}{2} k x^2\) (where \(k\) is the force constant) obeys:

-\frac{\hbar^2}{2m} \frac{d^2\psi(x)}{dx^2} + \frac{1}{2} k x^2 \psi(x) = E \psi(x)

Introducing the classical harmonic oscillator frequency \(\omega = \sqrt{\frac{k}{m}}\):

\frac{d^2\psi(x)}{dx^2} + \frac{2m}{\hbar^2} \left( E - \frac{1}{2} m \omega^2 x^2 \right) \psi(x) = 0

Dimensionless Coordinate Transformation

Define the dimensionless coordinate \(y\) and dimensionless energy eigenvalue \(\epsilon\):

\alpha = \frac{m\omega}{\hbar}, \quad y = \sqrt{\alpha} x = \left(\frac{m\omega}{\hbar}\right)^{1/2} x

\epsilon = \frac{2E}{\hbar\omega}

Transforming derivatives using the chain rule (\(\frac{d}{dx} = \sqrt{\alpha} \frac{d}{dy}\), \(\frac{d^2}{dx^2} = \alpha \frac{d^2}{dy^2}\)):

\alpha \frac{d^2\psi}{dy^2} + \left( \alpha \epsilon - \alpha y^2 \right) \psi = 0

Dividing throughout by \(\alpha\):

\frac{d^2\psi(y)}{dy^2} + (\epsilon - y^2) \psi(y) = 0

Asymptotic Analysis at Large Distances

As \(y \rightarrow \pm\infty\), \(y^2 \gg \epsilon\):

\frac{d^2\psi}{dy^2} - y^2 \psi \approx 0

The asymptotic solution that remains square-integrable (\(\psi \rightarrow 0\) as \(|y| \rightarrow \infty\)) is:

\psi_{\text{asympt}}(y) \sim e^{-y^2 / 2}

We therefore introduce the ansatz:

\psi(y) = H(y) e^{-y^2 / 2}

Substituting this ansatz into the dimensionless differential equation:

\frac{d\psi}{dy} = \left( H' - y H \right) e^{-y^2 / 2}

\frac{d^2\psi}{dy^2} = \left( H'' - 2y H' + (y^2 - 1) H \right) e^{-y^2 / 2}

Substituting into \(\psi'' + (\epsilon - y^2)\psi = 0\) yields Hermite's differential equation:

H''(y) - 2y H'(y) + (\epsilon - 1) H(y) = 0

Truncation of the Power Series & Energy Quantization

Seeking a power series solution \(H(y) = \sum_{j=0}^\infty c_j y^j\), substituting into Hermite's equation gives the recurrence relation:

c_{j+2} = \frac{2j + 1 - \epsilon}{(j + 1)(j + 2)} c_j

For the wavefunction to remain square-integrable as \(y \rightarrow \infty\), the infinite series must terminate into a finite polynomial of degree \(v\). Requiring the numerator to vanish at \(j = v\):

2v + 1 - \epsilon = 0 \implies \epsilon = 2v + 1

Substituting \(\epsilon = \frac{2E}{\hbar\omega}\):

\frac{2E}{\hbar\omega} = 2v + 1 \implies E_v = \left( v + \frac{1}{2} \right) \hbar\omega = \left( v + \frac{1}{2} \right) h \nu_0, \quad v = 0, 1, 2, 3, \dots

The energy levels are equally spaced by \(\hbar\omega\), with a ground-state zero-point energy of \(E_0 = \frac{1}{2}\hbar\omega\).

§3.2 Hermite Polynomials, Wavefunction Parity & Classical Turning Points

The solutions \(H_v(y)\) of Hermite's differential equation are the classical orthogonal Hermite polynomials.

Rodrigues' Formula and Generating Function

The Hermite polynomials \(H_v(y)\) can be generated systematically via Rodrigues' formula:

H_v(y) = (-1)^v e^{y^2} \frac{d^v}{dy^v} \left( e^{-y^2} \right)

Alternatively, from the generating function \(G(y, s) = e^{2ys - s^2} = \sum_{v=0}^\infty \frac{H_v(y)}{v!} s^v\).

The first six Hermite polynomials are:

\(H_0(y) = 1\)
\(H_1(y) = 2y\)
\(H_2(y) = 4y^2 - 2\)
\(H_3(y) = 8y^3 - 12y\)
\(H_4(y) = 16y^4 - 48y^2 + 12\)
\(H_5(y) = 32y^5 - 160y^3 + 120y\)

Normalized Vibrational Wavefunctions

The normalized harmonic oscillator wavefunctions are:

\psi_v(x) = N_v H_v(\sqrt{\alpha} x) \exp\left( -\frac{\alpha x^2}{2} \right)

where the normalization constant is derived from \(\int_{-\infty}^\infty H_v^2(y) e^{-y^2} dy = \sqrt{\pi} 2^v v!\):

N_v = \left( \frac{\alpha}{\pi} \right)^{1/4} \frac{1}{\sqrt{2^v v!}} = \left( \frac{m\omega}{\pi\hbar} \right)^{1/4} \frac{1}{\sqrt{2^v v!}}

Symmetry and Nodal Properties

Wavefunction Parity: Because \(V(-x) = V(x)\), the states possess definite parity:

\psi_v(-x) = (-1)^v \psi_v(x)

Even \(v\) states are even (symmetric); odd \(v\) states are odd (antisymmetric).

Nodes: State \(\psi_v(x)\) possesses exactly \(v\) nodes (real roots of \(H_v\)).

Classical Turning Points & Quantum Penetration

In classical mechanics, an oscillator with energy \(E_v\) cannot exceed the amplitude where potential energy equals total energy:

E_v = \frac{1}{2} k x_c^2 \implies x_c = \pm \sqrt{\frac{2E_v}{k}} = \pm \sqrt{\frac{(2v + 1)\hbar\omega}{m\omega^2}} = \pm \sqrt{\frac{2v + 1}{\alpha}}

In quantum mechanics:

For \(|x| < x_c\), the kinetic energy is positive (\(E > V(x)\)), and the wavefunction is oscillatory.
For \(|x| > x_c\), the kinetic energy is negative (\(E < V(x)\)). The Gaussian factor \(\exp(-\alpha x^2 / 2)\) decays exponentially into the classically forbidden barrier.
For the ground state (\(v = 0\)), integrating \(|\psi_0(x)|^2\) beyond \(x_c = 1/\sqrt{\alpha}\) reveals that there is a 15.7% probability of finding the quantum particle outside the classical turning points!

§3.3 Ladder Operator Formalism: Raising & Lowering Algebra

Paul Dirac introduced an elegant algebraic method to solve the harmonic oscillator without differential equations using non-Hermitian creation and annihilation operators.

Definition of Ladder Operators

Define the dimensionless annihilation (lowering) operator \(\hat{a}\) and creation (raising) operator \(\hat{a}^\dagger\):

\hat{a} = \sqrt{\frac{m\omega}{2\hbar}} \left( \hat{x} + \frac{i}{m\omega} \hat{p} \right) = \frac{1}{\sqrt{2}} \left( \hat{y} + i \hat{p}_y \right)

\hat{a}^\dagger = \sqrt{\frac{m\omega}{2\hbar}} \left( \hat{x} - \frac{i}{m\omega} \hat{p} \right) = \frac{1}{\sqrt{2}} \left( \hat{y} - i \hat{p}_y \right)

Inverting for position and momentum operators:

\hat{x} = \sqrt{\frac{\hbar}{2m\omega}} (\hat{a} + \hat{a}^\dagger)

\hat{p} = -i \sqrt{\frac{m\hbar\omega}{2}} (\hat{a} - \hat{a}^\dagger)

The Fundamental Commutation Relation

Using \([\hat{x}, \hat{p}] = i\hbar\):

[\hat{a}, \hat{a}^\dagger] = \frac{m\omega}{2\hbar} \left[ \hat{x} + \frac{i}{m\omega}\hat{p}, \hat{x} - \frac{i}{m\omega}\hat{p} \right] = \frac{m\omega}{2\hbar} \left( -\frac{i}{m\omega}[\hat{x}, \hat{p}] + \frac{i}{m\omega}[\hat{p}, \hat{x}] \right)

[\hat{a}, \hat{a}^\dagger] = \frac{m\omega}{2\hbar} \left( -\frac{i}{m\omega}(i\hbar) - \frac{i}{m\omega}(i\hbar) \right) = \frac{m\omega}{2\hbar} \left( \frac{2\hbar}{m\omega} \right) = 1

Thus:

[\hat{a}, \hat{a}^\dagger] = \hat{I}

The Number Operator and Hamiltonian

The product \(\hat{a}^\dagger \hat{a}\) is the Hermitian number operator \(\hat{N}\):

\hat{N} = \hat{a}^\dagger \hat{a}

Expressing the Hamiltonian in terms of ladder operators:

\hat{a}^\dagger \hat{a} = \frac{m\omega}{2\hbar} \left( \hat{x}^2 + \frac{\hat{p}^2}{m^2\omega^2} + \frac{i}{m\omega}(\hat{x}\hat{p} - \hat{p}\hat{x}) \right) = \frac{1}{\hbar\omega} \left( \frac{\hat{p}^2}{2m} + \frac{1}{2}m\omega^2\hat{x}^2 \right) - \frac{1}{2}

Therefore:

\hat{H} = \hbar\omega \left( \hat{a}^\dagger \hat{a} + \frac{1}{2} \right) = \hbar\omega \left( \hat{N} + \frac{1}{2} \right)

Action on Number Eigenstates \(|v\rangle\)

Let \(|v\rangle\) be an eigenstate of \(\hat{N}\) with eigenvalue \(v\): \(\hat{N}|v\rangle = v|v\rangle\). From the commutator \([\hat{N}, \hat{a}] = [\hat{a}^\dagger \hat{a}, \hat{a}] = [\hat{a}^\dagger, \hat{a}]\hat{a} = -\hat{a}\):

\hat{N} (\hat{a}|v\rangle) = (\hat{a}\hat{N} - \hat{a}) |v\rangle = (v - 1) (\hat{a}|v\rangle)

Similarly, \([\hat{N}, \hat{a}^\dagger] = +\hat{a}^\dagger \implies \hat{N} (\hat{a}^\dagger|v\rangle) = (v + 1) (\hat{a}^\dagger|v\rangle)\). Evaluating normalization:

\hat{a}|v\rangle = \sqrt{v} |v - 1\rangle

\hat{a}^\dagger|v\rangle = \sqrt{v + 1} |v + 1\rangle

Because the norm \(\langle v | \hat{a}^\dagger \hat{a} | v \rangle = v \ge 0\), the spectrum must terminate at \(v = 0\) via:

\hat{a}|0\rangle = 0

This algebraic condition directly yields the ground-state differential equation \(\left(x + \frac{\hbar}{m\omega}\frac{d}{dx}\right)\psi_0 = 0 \implies \psi_0(x) \propto e^{-\alpha x^2 / 2}\). Any state \(|v\rangle\) can be constructed by repeated application of the creation operator:

|v\rangle = \frac{(\hat{a}^\dagger)^v}{\sqrt{v!}} |0\rangle

§3.4 Three-Dimensional Rigid Rotor & Spherical Harmonics

The rigid rotor describes the rotational motion of diatomic and polyatomic molecules with fixed internuclear bond lengths.

Hamiltonian in Spherical Polar Coordinates

Consider two masses \(m_1\) and \(m_2\) separated by a fixed bond distance \(R_0\). Reducing to a one-body problem of reduced mass \(\mu = \frac{m_1 m_2}{m_1 + m_2}\) with moment of inertia \(I = \mu R_0^2\):

\hat{H}_{\text{rot}} = \frac{\hat{\mathbf{L}}^2}{2I}

where \(\hat{\mathbf{L}}^2\) is the square of the orbital angular momentum operator:

\hat{\mathbf{L}}^2 = -\hbar^2 \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 TISE is:

\hat{\mathbf{L}}^2 Y(\theta, \phi) = 2 I E Y(\theta, \phi) = \hbar^2 l(l + 1) Y(\theta, \phi)

with \(E_l = \frac{\hbar^2 l(l + 1)}{2I}\).

Separation of Variables & Associated Legendre Polynomials

Using separation \(Y(\theta, \phi) = \Theta(\theta) \Phi(\phi)\):

1. Azimuthal Equation:

\frac{d^2\Phi}{d\phi^2} = -m^2 \Phi \implies \Phi_m(\phi) = \frac{1}{\sqrt{2\pi}} e^{i m \phi}, \quad m \in \{0, \pm 1, \pm 2, \dots\}

2. Polar Equation:

Substituting \(\Phi\) yields the Associated Legendre differential equation for \(\Theta(\theta)\). Solutions that remain finite at the poles (\(\theta = 0, \pi\)) require:

l \in \{0, 1, 2, 3, \dots\}, \quad m \in \{-l, -l+1, \dots, +l\}

giving the Associated Legendre functions \(P_l^{|m|}(\cos\theta)\).

Orthonormal Spherical Harmonics

The joint angular eigenfunctions are the spherical harmonics \(Y_l^m(\theta, \phi)\):

Y_l^m(\theta, \phi) = (-1)^m \sqrt{\frac{2l + 1}{4\pi} \frac{(l - m)!}{(l + m)!}} P_l^m(\cos\theta) e^{i m \phi}

They satisfy simultaneous eigenvalue equations:

\hat{\mathbf{L}}^2 Y_l^m(\theta, \phi) = l(l + 1) \hbar^2 Y_l^m(\theta, \phi)

\hat{L}_z Y_l^m(\theta, \phi) = m \hbar Y_l^m(\theta, \phi)

and orthonormal integration:

\int_0^{2\pi} d\phi \int_0^\pi d\theta \sin\theta Y_l^{m*}(\theta, \phi) Y_{l'}^{m'}(\theta, \phi) = \delta_{ll'} \delta_{mm'}

Each rotational energy level \(E_l\) is \((2l + 1)\)-fold degenerate corresponding to the allowed orientations of \(m \in \{-l, \dots, +l\}\) in space (space quantization).

§3.5 Angular Momentum Commutation Algebra & Space Quantization

Angular momentum in quantum mechanics is defined fundamentally by its commutation relations rather than by classical cross products.

Canonical Commutators

The Cartesian components of orbital angular momentum \(\hat{\mathbf{L}} = \hat{\mathbf{r}} \times \hat{\mathbf{p}}\) satisfy:

[\hat{L}_x, \hat{L}_y] = i\hbar \hat{L}_z

[\hat{L}_y, \hat{L}_z] = i\hbar \hat{L}_x

[\hat{L}_z, \hat{L}_x] = i\hbar \hat{L}_y

Compactly written using the Levi-Civita permutation symbol \(\epsilon_{ijk}\):

[\hat{L}_i, \hat{L}_j] = i\hbar \sum_k \epsilon_{ijk} \hat{L}_k

Because no two Cartesian components commute, they cannot be measured simultaneously. By the Robertson-Schrödinger uncertainty relation:

\Delta L_x \cdot \Delta L_y \ge \frac{\hbar}{2} |\langle \hat{L}_z \rangle|

Commutation with the Total Angular Momentum Operator

The scalar operator \(\hat{\mathbf{L}}^2 = \hat{L}_x^2 + \hat{L}_y^2 + \hat{L}_z^2\) commutes with all three components:

[\hat{\mathbf{L}}^2, \hat{L}_x] = [\hat{\mathbf{L}}^2, \hat{L}_y] = [\hat{\mathbf{L}}^2, \hat{L}_z] = 0

Proof for \(\hat{L}_z\):

[\hat{\mathbf{L}}^2, \hat{L}_z] = [\hat{L}_x^2, \hat{L}_z] + [\hat{L}_y^2, \hat{L}_z] + [\hat{L}_z^2, \hat{L}_z]

= \hat{L}_x [\hat{L}_x, \hat{L}_z] + [\hat{L}_x, \hat{L}_z]\hat{L}_x + \hat{L}_y [\hat{L}_y, \hat{L}_z] + [\hat{L}_y, \hat{L}_z]\hat{L}_y

= \hat{L}_x (-i\hbar \hat{L}_y) + (-i\hbar \hat{L}_y)\hat{L}_x + \hat{L}_y (i\hbar \hat{L}_x) + (i\hbar \hat{L}_x)\hat{L}_y = 0

Therefore, \(\hat{\mathbf{L}}^2\) and exactly one Cartesian projection (conventionally \(\hat{L}_z\)) form a Complete Set of Commuting Observables (CSCO).

Angular Momentum Ladder Operators \(\hat{L}_\pm\)

Define the ladder operators:

\hat{L}_+ = \hat{L}_x + i \hat{L}_y, \quad \hat{L}_- = \hat{L}_x - i \hat{L}_y

They satisfy:

[\hat{L}_z, \hat{L}_\pm] = \pm \hbar \hat{L}_\pm

[\hat{L}_+, \hat{L}_-] = 2\hbar \hat{L}_z

\hat{\mathbf{L}}^2 = \hat{L}_- \hat{L}_+ + \hat{L}_z^2 + \hbar \hat{L}_z = \hat{L}_+ \hat{L}_- + \hat{L}_z^2 - \hbar \hat{L}_z

Acting on mutual eigenstates \(|l, m\rangle\):

\hat{L}_\pm |l, m\rangle = \hbar \sqrt{l(l + 1) - m(m \pm 1)} |l, m \pm 1\rangle

Vector Model & Precessional Cones of Space Quantization

In state \(|l, m\rangle\), the magnitude of the angular momentum vector is \(|\mathbf{L}| = \sqrt{l(l + 1)} \hbar\), while its projection along the \(z\)-axis is \(L_z = m\hbar\). Because \(|m| \le l < \sqrt{l(l+1)}\), the vector \(\mathbf{L}\) can never align purely along the \(z\)-axis! The angle \(\theta\) between \(\mathbf{L}\) and the quantization axis is:

\cos\theta = \frac{L_z}{|\mathbf{L}|} = \frac{m}{\sqrt{l(l + 1)}}

The perpendicular components \(L_x\) and \(L_y\) remain completely indeterminate, causing \(\mathbf{L}\) to precess on a cone of half-angle \(\theta\) around the \(z\)-axis.

§3.6 Matrix Representations of Angular Momentum Operators

Using the orthonormal basis states \(\{|l, m\rangle\}\), angular momentum operators can be cast into explicit finite-dimensional matrix representations.

General Matrix Elements in the \(|l, m\rangle\) Basis

For a fixed angular momentum quantum number \(l\), the Hilbert subspace has dimension \(2l + 1\).

1. Matrix elements of \(\hat{L}_z\):

\langle l, m' | \hat{L}_z | l, m \rangle = m \hbar \delta_{m, m'}

\(\hat{L}_z\) is represented by a diagonal matrix.

2. Matrix elements of \(\hat{L}_+\) and \(\hat{L}_-\):

\langle l, m' | \hat{L}_+ | l, m \rangle = \hbar \sqrt{l(l + 1) - m(m + 1)} \delta_{m', m+1}

\langle l, m' | \hat{L}_- | l, m \rangle = \hbar \sqrt{l(l + 1) - m(m - 1)} \delta_{m', m-1}

3. Cartesian Components \(\hat{L}_x\) and \(\hat{L}_y\):

\hat{L}_x = \frac{1}{2}(\hat{L}_+ + \hat{L}_-), \quad \hat{L}_y = \frac{1}{2i}(\hat{L}_+ - \hat{L}_-)

Explicit Matrices for \(l = 1\) (Three-Dimensional Representation)

In the basis ordered as \(\{|1, +1\rangle, |1, 0\rangle, |1, -1\rangle\}\):

\mathbf{L}_z = \hbar \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & -1 \end{pmatrix}

Evaluating ladder operators:

\hat{L}_+|1, 0\rangle = \hbar\sqrt{1(2) - 0(1)}|1, 1\rangle = \sqrt{2}\hbar |1, 1\rangle

\hat{L}_+|1, -1\rangle = \hbar\sqrt{1(2) - (-1)(0)}|1, 0\rangle = \sqrt{2}\hbar |1, 0\rangle

\mathbf{L}_+ = \sqrt{2}\hbar \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix}, \quad \mathbf{L}_- = \mathbf{L}_+^\dagger = \sqrt{2}\hbar \begin{pmatrix} 0 & 0 & 0 \\ 1 & 0 & 0 \\ 0 & 1 & 0 \end{pmatrix}

Consequently:

\mathbf{L}_x = \frac{\hbar}{\sqrt{2}} \begin{pmatrix} 0 & 1 & 0 \\ 1 & 0 & 1 \\ 0 & 1 & 0 \end{pmatrix}, \quad \mathbf{L}_y = \frac{\hbar}{\sqrt{2}} \begin{pmatrix} 0 & -i & 0 \\ i & 0 & -i \\ 0 & i & 0 \end{pmatrix}

Evaluating \(\mathbf{L}_x^2 + \mathbf{L}_y^2 + \mathbf{L}_z^2\):

\mathbf{L}^2 = 2\hbar^2 \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} = l(l+1)\hbar^2 \mathbf{I}

which rigorously confirms the operator algebra.

§3.7 Rotational Transitions & Microwave Absorption Selection Rules

The quantum rigid rotor model directly governs pure rotational transitions observed in microwave spectroscopy.

Electric Dipole Transition Operator

The interaction between rotating molecules and electromagnetic radiation is mediated by the electric dipole transition moment operator \(\hat{\boldsymbol{\mu}} = \mu_0 \hat{\mathbf{u}}\), where \(\mu_0\) is the permanent molecular dipole moment:

\mathbf{M}_{l'm', lm} = \langle l', m' | \hat{\boldsymbol{\mu}} | l, m \rangle = \mu_0 \int Y_{l'}^{m'*}(\theta, \phi) \mathbf{u} Y_l^m(\theta, \phi) d\Omega

The Gross Selection Rule

If a molecule possesses zero permanent electric dipole moment (\(\mu_0 = 0\), as in homonuclear diatomics \(\text{N}_2, \text{O}_2, \text{H}_2\) or spherical rotors \(\text{CH}_4, \text{SF}_6\)), the transition dipole matrix element vanishes identically for all states:

\mu_0 = 0 \implies \mathbf{M}_{l'm', lm} = 0

Gross Selection Rule: A molecule must possess a non-zero permanent electric dipole moment (\(\mu_0 \neq 0\)) to exhibit an electric-dipole allowed pure rotational microwave spectrum.

The Specific Selection Rules

For linear rotors with \(\mu_0 \neq 0\), the dipole component along the space-fixed \(z\)-axis is \(\hat{\mu}_z = \mu_0 \cos\theta = \mu_0 \sqrt{\frac{4\pi}{3}} Y_1^0(\theta, \phi)\). The transition matrix element is:

\langle l', m' | \cos\theta | l, m \rangle \propto \int_0^{2\pi} e^{i (m - m')\phi} d\phi \int_0^\pi P_{l'}^{m'}(\cos\theta) \cos\theta P_l^m(\cos\theta) \sin\theta d\theta

Using the recurrence relation \((2l + 1) x P_l^m(x) = (l - m + 1) P_{l+1}^m(x) + (l + m) P_{l-1}^m(x)\):

The integral is strictly non-zero if and only if:

\Delta l = l' - l = \pm 1

\Delta m = m' - m = 0, \pm 1

For absorption transitions (\(\Delta l = +1\), denoted \(J \rightarrow J + 1\)):

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

where \(B = \frac{h}{8\pi^2 c I}\) is the rotational constant in \(\text{cm}^{-1}\). This yields equally spaced spectral absorption lines separated by \(2B\).

§3.8 Centrifugal Distortion, Rovibrational Coupling & P/R Branches

In real diatomic molecules, rotational and vibrational motions are coupled because stretching the chemical bond increases the moment of inertia, altering rotational energy levels.

The Vibrating Rotor Energy Term Values

In high-resolution infrared spectroscopy, the energy term value \(T(v, J) = E(v, J) / hc\) in wavenumbers (\(\text{cm}^{-1}\)) is expanded as:

\[T(v, J) = G(v) + F_v(J)\]

where the vibrational term includes anharmonicity:

G(v) = \tilde{\omega}_e \left( v + \frac{1}{2} \right) - \tilde{\omega}_e x_e \left( v + \frac{1}{2} \right)^2 + \dots

and the effective rotational term for vibrational level \(v\) includes centrifugal distortion:

F_v(J) = B_v J(J + 1) - D_v [J(J + 1)]^2 + \dots

The effective rotational constant \(B_v\) decreases in higher vibrational states due to increased average bond length:

B_v = B_e - \alpha_e \left( v + \frac{1}{2} \right)

where \(\alpha_e > 0\) is the vibration-rotation coupling constant, and \(D_v \approx \frac{4 B_e^3}{\tilde{\omega}_e^2}\) is the centrifugal distortion constant.

Selection Rules and Rovibrational Branches

For electric dipole absorption in a gas-phase diatomic molecule, the simultaneous selection rules are:

\Delta v = +1 \quad (\text{Fundamental band } v = 0 \rightarrow 1), \quad \Delta J = \pm 1

1. R-Branch (\(\Delta J = J' - J'' = +1\)): Transition from \(J''\) to \(J' = J'' + 1\):

\tilde{\nu}_R(J'') = T(1, J'' + 1) - T(0, J'') = \tilde{\nu}_0 + (2 B_1) + (3 B_1 - B_0) J'' + (B_1 - B_0) J''^2

Lines appear at frequencies higher than the band origin \(\tilde{\nu}_0\).

2. P-Branch (\(\Delta J = J' - J'' = -1\)): Transition from \(J''\) to \(J' = J'' - 1\):

\tilde{\nu}_P(J'') = T(1, J'' - 1) - T(0, J'') = \tilde{\nu}_0 - (B_1 + B_0) J'' + (B_1 - B_0) J''^2

Lines appear at frequencies lower than the band origin \(\tilde{\nu}_0\).

3. Band Origin and Missing Q-Branch:

For \(^1\Sigma\) diatomics (zero orbital angular momentum along the internuclear axis), \(\Delta J = 0\) is strictly forbidden. The central line at \(\tilde{\nu}_0\) is completely absent, producing the characteristic gap between the P and R branches.

Band Head Formation

Because \(B_1 < B_0\), the quadratic term \((B_1 - B_0) J''^2\) is negative. In the R-branch, line spacing decreases as \(J''\) increases until the spacing reverses sign, producing a maximum frequency pileup termed the band head.

## Advanced Mathematical Supplement: Lie Algebra SU(2) & Wigner D-Matrices

Lie Algebraic Formulation of Angular Momentum & SU(2) Representations

Orbital and spin angular momentum operators are the infinitesimal generators of spatial rotations, forming the Lie algebra \(\mathfrak{su}(2) \cong \mathfrak{so}(3)\).

The Lie Algebra Commutation Relations

The components satisfy the Lie bracket relations:

[\hat{J}_i, \hat{J}_j] = i\hbar \sum_{k} \epsilon_{i j k} \hat{J}_k

Defining ladder operators \(\hat{J}_\pm = \hat{J}_x \pm i \hat{J}_y\):

[\hat{J}_z, \hat{J}_\pm] = \pm \hbar \hat{J}_\pm, \quad [\hat{J}_+, \hat{J}_-] = 2\hbar \hat{J}_z

The Casimir operator of the algebra is \(\hat{J}^2 = \hat{J}_x^2 + \hat{J}_y^2 + \hat{J}_z^2\), which commutes with all generators: \([\hat{J}^2, \hat{J}_k] = 0\).

Finite Rotations and Wigner D-Matrices

A finite rotation of a quantum state by Euler angles \((\alpha, \beta, \gamma)\) in the \(z\)-\(y'\)-\(z''\) convention is represented by the unitary rotation operator:

\hat{\mathcal{D}}(\alpha, \beta, \gamma) = \exp\left( -\frac{i \alpha \hat{J}_z}{\hbar} \right) \exp\left( -\frac{i \beta \hat{J}_y}{\hbar} \right) \exp\left( -\frac{i \gamma \hat{J}_z}{\hbar} \right)

In the orthonormal angular momentum basis \(|j, m\rangle\), the matrix elements are the Wigner D-matrices:

D_{m' m}^j(\alpha, \beta, \gamma) = \langle j, m' | \hat{\mathcal{D}}(\alpha, \beta, \gamma) | j, m \rangle = e^{-i m' \alpha} d_{m' m}^j(\beta) e^{-i m \gamma}

where Wigner's small \(d\)-matrix is given by the Jacobi polynomial formula:

d_{m' m}^j(\beta) = \sqrt{\frac{(j+m')!(j-m')!}{(j+m)!(j-m)!}} \left(\sin\frac{\beta}{2}\right)^{m'-m} \left(\cos\frac{\beta}{2}\right)^{m'+m} P_{j-m'}^{(m'-m, m'+m)}(\cos\beta)

The Wigner-Eckart Theorem

For any irreducible spherical tensor operator \(\hat{T}_q^{(k)}\) of rank \(k\), all spatial matrix elements factor into a geometric Clebsch-Gordan coefficient and a physical reduced matrix element:

\langle j', m' | \hat{T}_q^{(k)} | j, m \rangle = (-1)^{j' - m'} \begin{pmatrix} j' & k & j \\ -m' & q & m \end{pmatrix} \langle j' || \hat{\mathbf{T}}^{(k)} || j \rangle

The Wigner-Eckart theorem proves that all spectroscopic selection rules (\(\Delta J, \Delta M\)) are governed strictly by the rotational symmetry of the transition operator, independent of the internal radial Hamiltonian.

## Research Monograph: The Geometric Berry Phase & Born-Oppenheimer Conical Intersections

Adiabatic Evolution & The Geometric Phase

When a quantum system with Hamiltonian \(\hat{H}(\mathbf{R})\) is transported slowly (adiabatically) along a closed contour \(\mathcal{C}\) in parameter space \(\mathbf{R}(t)\), Michael Berry discovered in 1984 that upon returning to the starting point \(\mathbf{R}(T) = \mathbf{R}(0)\), the state acquires an invariant geometric phase in addition to the standard dynamical phase:

|\psi(T)\rangle = \exp\left( -\frac{i}{\hbar} \int_0^T E_n(\mathbf{R}(t)) dt \right) \exp\left( i \gamma_n[\mathcal{C}] \right) |\psi(0)\rangle

The Berry Phase \(\gamma_n[\mathcal{C}]\) depends strictly on the geometry of the path in parameter space:

\gamma_n[\mathcal{C}] = \oint_{\mathcal{C}} \mathbf{A}_n(\mathbf{R}) \cdot d\mathbf{R} = \iint_{\mathcal{S}} \boldsymbol{\Omega}_n(\mathbf{R}) \cdot d\mathbf{S}

where the Berry Connection is \(\mathbf{A}_n(\mathbf{R}) = i \langle n(\mathbf{R}) | \nabla_{\mathbf{R}} | n(\mathbf{R}) \rangle\) and the Berry Curvature is:

\boldsymbol{\Omega}_n(\mathbf{R}) = \nabla_{\mathbf{R}} \times \mathbf{A}_n(\mathbf{R}) = i \sum_{m \ne n} \frac{\langle n | \nabla \hat{H} | m \rangle \times \langle m | \nabla \hat{H} | n \rangle}{(E_m(\mathbf{R}) - E_n(\mathbf{R}))^2}

Conical Intersections as Monopoles in Nuclear Space

In molecular chemistry, two adiabatic electronic potential energy surfaces of the same spatial and spin symmetry can intersect at a point in a two-dimensional branching plane \((X_1, X_2)\), forming a Conical Intersection (CoIn). Near the intersection point, the electronic Hamiltonian is:

\mathbf{H}_{\text{elec}}(\mathbf{R}) = \begin{pmatrix} c_1 X_1 & c_2 X_2 \\ c_2 X_2 & -c_1 X_1 \end{pmatrix}

The adiabatic energy surfaces form a double cone:

E_\pm(\mathbf{R}) = \pm \sqrt{c_1^2 X_1^2 + c_2^2 X_2^2}

The conical intersection acts as an effective Dirac magnetic monopole in nuclear coordinate space! When the nuclei traverse a closed loop encircling the conical intersection, the electronic wavefunction changes sign:

\gamma = \oint \mathbf{A} \cdot d\mathbf{R} = \pm \pi \implies e^{i \gamma} = -1

To keep the total electron-nuclear wavefunction single-valued, the nuclear vibrational wavefunction must also change sign, acquiring a half-integer angular momentum quantum number: the Molecular Aharonov-Bohm Effect.

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.

Foundational Example 3.1: Problem 1: Carbon Monoxide Vibrational Zero-Point Energy & Classical Turning Point

The fundamental vibrational wavenumber of carbon monoxide (\(^{12}\text{C}^{16}\text{O}\)) is \(\tilde{\nu}_0 = 2170\text{ cm}^{-1}\).

Calculate the effective force constant \(k\) of the \(\text{C}\equiv\text{O}\) triple bond (in \(\text{N/m}\)) using reduced mass \(\mu = \frac{m_{\text{C}} m_{\text{O}}}{m_{\text{C}} + m_{\text{O}}}\) with atomic masses \(m(^{12}\text{C}) = 12.000\text{ u}\) and \(m(^{16}\text{O}) = 15.995\text{ u}\).
Calculate the zero-point vibrational energy \(E_0\) in Joules and \(\text{kJ/mol}\).
Calculate the classical turning point amplitude \(x_c\) for the ground state (\(v = 0\)).

Comprehensive Multi-Step Solution:

Step 1: Reduced Mass and Force Constant Calculation

Calculate reduced mass \(\mu\):

\mu = \frac{12.000 \times 15.995}{12.000 + 15.995}\text{ u} = \frac{191.94}{27.995}\text{ u} \approx 6.85622\text{ u}

Convert to kilograms (\(1\text{ u} = 1.66054 \times 10^{-27}\text{ kg}\)):

\mu = 6.85622 \times 1.66054 \times 10^{-27}\text{ kg} \approx 1.1385 \times 10^{-26}\text{ kg}

The vibrational frequency is:

\nu_0 = c \tilde{\nu}_0 = (2.99792 \times 10^{10}\text{ cm/s}) \times 2170\text{ cm}^{-1} \approx 6.5055 \times 10^{13}\text{ s}^{-1}

The angular frequency is \(\omega = 2\pi \nu_0 = 2\pi (6.5055 \times 10^{13}) \approx 4.0875 \times 10^{14}\text{ rad/s}\). The force constant is:

k = \mu \omega^2 = (1.1385 \times 10^{-26}\text{ kg}) \times (4.0875 \times 10^{14}\text{ s}^{-1})^2 \approx 1902.1\text{ N/m}

This large force constant (\(\sim 1902\text{ N/m}\)) reflects the strong carbon-oxygen triple bond.

Step 2: Zero-Point Energy Calculation

E_0 = \frac{1}{2} h \nu_0 = \frac{1}{2} (6.62607 \times 10^{-34}\text{ J}\cdot\text{s})(6.5055 \times 10^{13}\text{ s}^{-1}) \approx 2.1553 \times 10^{-20}\text{ J}

Per mole:

E_{0, \text{molar}} = E_0 N_A = (2.1553 \times 10^{-20}\text{ J}) \times (6.02214 \times 10^{23}\text{ mol}^{-1}) \approx 12979\text{ J/mol} \approx 12.98\text{ kJ/mol}

Step 3: Classical Turning Point Amplitude \(x_c\)

At the turning point for \(v = 0\):

\frac{1}{2} k x_c^2 = E_0 \implies x_c = \sqrt{\frac{2 E_0}{k}} = \sqrt{\frac{\hbar\omega}{k}} = \sqrt{\frac{\hbar}{\mu\omega}}

x_c = \sqrt{\frac{2 \times 2.1553 \times 10^{-20}\text{ J}}{1902.1\text{ N/m}}} = \sqrt{2.2662 \times 10^{-23}\text{ m}^2} \approx 4.76 \times 10^{-12}\text{ m} = 0.0476\text{ Å}

The zero-point oscillation stretches or compresses the equilibrium bond length (\(R_e = 1.128\text{ Å}\)) by only \(\sim 4.2\%\).

Foundational Example 3.2: Problem 2: Matrix Elements of Position and Momentum using Ladder Operators

Using the ladder operator definitions:

\hat{x} = \sqrt{\frac{\hbar}{2m\omega}} (\hat{a} + \hat{a}^\dagger), \quad \hat{p} = -i\sqrt{\frac{m\hbar\omega}{2}} (\hat{a} - \hat{a}^\dagger)

Calculate the matrix elements \(\langle v' | \hat{x} | v \rangle\) and \(\langle v' | \hat{p} | v \rangle\) between arbitrary vibrational states.
Derive the selection rule for electric dipole vibrational transitions.
Calculate \(\langle v | \hat{x}^2 | v \rangle\) and prove that the average potential energy \(\langle V \rangle\) equals the average kinetic energy \(\langle T \rangle = \frac{1}{2} E_v\) (Virial Theorem).

Comprehensive Multi-Step Solution:

Step 1: Matrix Elements of \(\hat{x}\) and \(\hat{p}\)

Recall that:

\hat{a}|v\rangle = \sqrt{v}|v - 1\rangle, \quad \hat{a}^\dagger|v\rangle = \sqrt{v + 1}|v + 1\rangle

For position \(\hat{x}\):

\hat{x}|v\rangle = \sqrt{\frac{\hbar}{2m\omega}} \left( \sqrt{v}|v - 1\rangle + \sqrt{v + 1}|v + 1\rangle \right)

Taking the inner product with \(\langle v'|\):

\langle v' | \hat{x} | v \rangle = \sqrt{\frac{\hbar}{2m\omega}} \left( \sqrt{v} \delta_{v', v-1} + \sqrt{v+1} \delta_{v', v+1} \right)

Similarly for momentum \(\hat{p}\):

\langle v' | \hat{p} | v \rangle = -i\sqrt{\frac{m\hbar\omega}{2}} \left( \sqrt{v} \delta_{v', v-1} - \sqrt{v+1} \delta_{v', v+1} \right)

Step 2: Vibrational Dipole Selection Rule

In the harmonic approximation, the molecular dipole moment is expanded to first order in displacement \(x\):

\mu(x) = \mu_0 + \left(\frac{d\mu}{dx}\right)_0 x

The transition matrix element is:

\langle v' | \mu(x) | v \rangle = \left(\frac{d\mu}{dx}\right)_0 \langle v' | \hat{x} | v \rangle

From Step 1, \(\langle v' | \hat{x} | v \rangle\) is non-zero only when:

v' = v \pm 1 \iff \Delta v = \pm 1

Transitions with \(\Delta v = \pm 2, \pm 3, \dots\) (overtones) are strictly forbidden in a pure harmonic oscillator.

Step 3: Evaluation of \(\langle x^2 \rangle\) and the Virial Theorem

Expand \(\hat{x}^2\):

\hat{x}^2 = \frac{\hbar}{2m\omega} (\hat{a} + \hat{a}^\dagger)^2 = \frac{\hbar}{2m\omega} (\hat{a}^2 + \hat{a}\hat{a}^\dagger + \hat{a}^\dagger\hat{a} + (\hat{a}^\dagger)^2)

Because \(\hat{a}^2\) and \((\hat{a}^\dagger)^2\) change the quantum number by \(\pm 2\), their expectation values in state \(|v\rangle\) vanish:

\langle v | \hat{x}^2 | v \rangle = \frac{\hbar}{2m\omega} \langle v | (\hat{a}\hat{a}^\dagger + \hat{a}^\dagger\hat{a}) | v \rangle

Using \(\hat{a}\hat{a}^\dagger = \hat{a}^\dagger\hat{a} + 1\):

\hat{a}\hat{a}^\dagger + \hat{a}^\dagger\hat{a} = 2\hat{a}^\dagger\hat{a} + 1 = 2\hat{N} + 1

Therefore:

\langle x^2 \rangle_v = \frac{\hbar}{2m\omega} (2v + 1) = \left(v + \frac{1}{2}\right) \frac{\hbar}{m\omega}

Now calculate the average potential energy:

\langle V \rangle = \frac{1}{2} k \langle x^2 \rangle = \frac{1}{2} (m\omega^2) \left(v + \frac{1}{2}\right) \frac{\hbar}{m\omega} = \frac{1}{2} \left(v + \frac{1}{2}\right) \hbar\omega = \frac{1}{2} E_v

Because \(\langle T \rangle + \langle V \rangle = E_v\):

\langle T \rangle = E_v - \langle V \rangle = \frac{1}{2} E_v = \langle V \rangle

This proves the Virial Theorem (\(\langle T \rangle = \langle V \rangle\)) for a harmonic potential.

Intermediate Example 3.3: Problem 3: Rotational Constant and Bond Length of Hydrogen Chloride from Microwave Transitions

The pure rotational spectrum of gaseous \(^1\text{H}^{35}\text{Cl}\) exhibits a series of equidistant absorption lines in the far-infrared region with an average adjacent spacing \(\Delta\tilde{\nu} = 20.68\text{ cm}^{-1}\). Atomic masses: \(m(^1\text{H}) = 1.007825\text{ u}\), \(m(^{35}\text{Cl}) = 34.96885\text{ u}\).

Calculate the rotational constant \(B\) of \(^1\text{H}^{35}\text{Cl}\) in \(\text{cm}^{-1}\) and Joules.
Calculate the moment of inertia \(I\) (in \(\text{kg}\cdot\text{m}^2\)).
Calculate the equilibrium internuclear bond distance \(R_0\) (in Angstroms).
Predict the wavenumber of the \(J = 3 \rightarrow J = 4\) rotational transition.

Comprehensive Multi-Step Solution:

Step 1: Rotational Constant Calculation

For a rigid diatomic rotor, adjacent rotational lines are separated by:

\Delta\tilde{\nu} = 2B \implies B = \frac{\Delta\tilde{\nu}}{2} = \frac{20.68\text{ cm}^{-1}}{2} = 10.34\text{ cm}^{-1}

In energy units:

B(\text{Joules}) = h c B(\text{cm}^{-1}) = (6.62607 \times 10^{-34}\text{ J}\cdot\text{s})(2.99792 \times 10^{10}\text{ cm/s})(10.34\text{ cm}^{-1}) \approx 2.054 \times 10^{-22}\text{ J}

Step 2: Moment of Inertia

From the definition of the rotational constant:

B = \frac{h}{8\pi^2 c I} \implies I = \frac{h}{8\pi^2 c B}

I = \frac{6.62607 \times 10^{-34}\text{ J}\cdot\text{s}}{8 \pi^2 (2.99792 \times 10^{10}\text{ cm/s})(10.34\text{ cm}^{-1})} = \frac{6.62607 \times 10^{-34}}{2.4475 \times 10^{-11}} \approx 2.7073 \times 10^{-47}\text{ kg}\cdot\text{m}^2

Step 3: Internuclear Bond Length \(R_0\)

Calculate the reduced mass \(\mu\):

\mu = \frac{1.007825 \times 34.96885}{1.007825 + 34.96885}\text{ u} = \frac{35.2424}{35.9767}\text{ u} \approx 0.97959\text{ u}

In kilograms:

\mu = 0.97959 \times 1.66054 \times 10^{-27}\text{ kg} \approx 1.6266 \times 10^{-27}\text{ kg}

Since \(I = \mu R_0^2\):

R_0 = \sqrt{\frac{I}{\mu}} = \sqrt{\frac{2.7073 \times 10^{-47}\text{ kg}\cdot\text{m}^2}{1.6266 \times 10^{-27}\text{ kg}}} = \sqrt{1.6644 \times 10^{-20}\text{ m}^2} \approx 1.290 \times 10^{-10}\text{ m} = 1.290\text{ Å}

This precisely reproduces the experimental gas-phase bond length of \(\text{HCl}\) (\(1.275 - 1.29\text{ Å}\)).

Step 4: Wavenumber of \(J = 3 \rightarrow 4\) Transition

The transition wavenumber for \(J \rightarrow J + 1\) is:

\tilde{\nu}(3 \rightarrow 4) = 2B (3 + 1) = 8B = 8 \times 10.34\text{ cm}^{-1} = 82.72\text{ cm}^{-1}

Intermediate Example 3.4: Problem 4: Most Populated Rotational Energy Level at Thermal Equilibrium

At thermal equilibrium at temperature \(T = 300\text{ K}\), molecules are distributed among rotational states according to Boltzmann statistics with degeneracy \(g_J = 2J + 1\):

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

Derive the general formula for the rotational quantum number \(J_{\text{max}}\) corresponding to the most populated rotational level by treating \(J\) as a continuous variable.
For carbon monoxide (\(B = 1.931\text{ cm}^{-1}\)), calculate \(J_{\text{max}}\) at \(T = 300\text{ K}\).
Explain why the most intense spectral line in microwave absorption corresponds to transition from \(J_{\text{max}}\).

Comprehensive Multi-Step Solution:

Step 1: Derivation of \(J_{\text{max}}\)

Let the population function be:

f(J) = (2J + 1) \exp(-\beta J(J + 1))

where \(\beta = \frac{B h c}{k_B T}\). To find the maximum, set the first derivative to zero:

\frac{df}{dJ} = 2 \exp(-\beta J(J + 1)) + (2J + 1) [-\beta(2J + 1)] \exp(-\beta J(J + 1)) = 0

Dividing by the non-zero exponential factor:

2 - \beta (2J + 1)^2 = 0 \implies (2J + 1)^2 = \frac{2}{\beta} = \frac{2 k_B T}{B h c}

Taking the positive square root:

2J + 1 = \sqrt{\frac{2 k_B T}{B h c}} \implies J_{\text{max}} = \sqrt{\frac{k_B T}{2 B h c}} - \frac{1}{2}

Step 2: Numerical Calculation for CO at \(300\text{ K}\)

First calculate the dimensionless ratio \(\frac{k_B T}{B h c}\):

k_B T = (1.38065 \times 10^{-23}\text{ J/K})(300\text{ K}) = 4.14195 \times 10^{-21}\text{ J}

B h c = (1.931\text{ cm}^{-1})(6.62607 \times 10^{-34}\text{ J}\cdot\text{s})(2.99792 \times 10^{10}\text{ cm/s}) = 3.8358 \times 10^{-23}\text{ J}

\frac{k_B T}{B h c} = \frac{4.14195 \times 10^{-21}}{3.8358 \times 10^{-23}} \approx 107.98

Now evaluate \(J_{\text{max}}\):

J_{\text{max}} = \sqrt{\frac{107.98}{2}} - 0.5 = \sqrt{53.99} - 0.5 \approx 7.348 - 0.5 = 6.85

Rounding to the nearest integer:

J_{\text{max}} = 7

At \(300\text{ K}\), the rotational state \(J = 7\) is the most heavily populated state in gas-phase \(\text{CO}\).

Step 3: Connection to Spectral Intensity

Microwave absorption transition intensity is directly proportional to the population difference between states:

I(J \rightarrow J + 1) \propto N_J - N_{J+1} \approx N_J \frac{\Delta E}{k_B T}

Because \(N_J\) peaks at \(J = J_{\text{max}}\) while the transition energy \(2B(J+1)\) varies gently, the absorption spectrum envelope mirrors the rotational population distribution, displaying maximum absorption intensity around \(J = 7 \rightarrow 8\).

Advanced Example 3.5: Problem 5: Non-Commuting Angular Momentum and Uncertainty Products

Consider an angular momentum eigenstate \(|l, m\rangle\) where \(\hat{\mathbf{L}}^2|l, m\rangle = l(l+1)\hbar^2|l, m\rangle\) and \(\hat{L}_z|l, m\rangle = m\hbar|l, m\rangle\).

Prove that \(\langle \hat{L}_x \rangle = 0\) and \(\langle \hat{L}_y \rangle = 0\).
Calculate the expectation values \(\langle \hat{L}_x^2 \rangle\) and \(\langle \hat{L}_y^2 \rangle\) in terms of \(l, m, \hbar\).
Evaluate the uncertainty product \(\Delta L_x \cdot \Delta L_y\) and show that it satisfies the Robertson-Schrödinger uncertainty inequality:

\Delta L_x \cdot \Delta L_y \ge \frac{\hbar}{2} |\langle \hat{L}_z \rangle|

Identify under what condition \(\Delta L_x \cdot \Delta L_y\) achieves its minimum value.

Comprehensive Multi-Step Solution:

Step 1: Expectation Values of \(\hat{L}_x\) and \(\hat{L}_y\)

Using ladder operators:

\hat{L}_x = \frac{1}{2}(\hat{L}_+ + \hat{L}_-), \quad \hat{L}_y = \frac{1}{2i}(\hat{L}_+ - \hat{L}_-)

Because \(\hat{L}_\pm |l, m\rangle \propto |l, m \pm 1\rangle\) and \(\langle l, m | l, m \pm 1 \rangle = 0\) by orthogonality:

\langle \hat{L}_x \rangle = \frac{1}{2} (\langle l, m | \hat{L}_+ | l, m \rangle + \langle l, m | \hat{L}_- | l, m \rangle) = 0

\langle \hat{L}_y \rangle = \frac{1}{2i} (\langle l, m | \hat{L}_+ | l, m \rangle - \langle l, m | \hat{L}_- | l, m \rangle) = 0

Step 2: Evaluation of \(\langle \hat{L}_x^2 \rangle\) and \(\langle \hat{L}_y^2 \rangle\)

By rotational symmetry about the \(z\)-axis:

\langle \hat{L}_x^2 \rangle = \langle \hat{L}_y^2 \rangle

Since \(\hat{\mathbf{L}}^2 = \hat{L}_x^2 + \hat{L}_y^2 + \hat{L}_z^2\):

\langle \hat{L}_x^2 \rangle + \langle \hat{L}_y^2 \rangle = \langle \hat{\mathbf{L}}^2 \rangle - \langle \hat{L}_z^2 \rangle = l(l + 1)\hbar^2 - m^2 \hbar^2

Equating the two equal transverse components:

2 \langle \hat{L}_x^2 \rangle = \hbar^2 [l(l + 1) - m^2] \implies \langle \hat{L}_x^2 \rangle = \langle \hat{L}_y^2 \rangle = \frac{\hbar^2}{2} [l(l + 1) - m^2]

Step 3: Uncertainty Product Evaluation

Since \(\langle \hat{L}_x \rangle = \langle \hat{L}_y \rangle = 0\):

\Delta L_x = \sqrt{\langle \hat{L}_x^2 \rangle} = \hbar \sqrt{\frac{l(l + 1) - m^2}{2}}

\Delta L_y = \sqrt{\langle \hat{L}_y^2 \rangle} = \hbar \sqrt{\frac{l(l + 1) - m^2}{2}}

The uncertainty product is:

\Delta L_x \cdot \Delta L_y = \frac{\hbar^2}{2} [l(l + 1) - m^2]

Testing against the uncertainty inequality:

\Delta L_x \cdot \Delta L_y \ge \frac{\hbar}{2} |\langle \hat{L}_z \rangle| = \frac{\hbar^2}{2} |m|

Subtracting the two expressions:

\frac{\hbar^2}{2} [l(l + 1) - m^2 - |m|] = \frac{\hbar^2}{2} [l(l + 1) - |m|(|m| + 1)]

Since \(|m| \le l\), \(l(l + 1) - |m|(|m| + 1) \ge 0\), proving the inequality holds universally.

Step 4: Minimum Uncertainty Condition

The minimum occurs when \(|m| = l\) (maximum alignment along \(z\)):

\[l(l + 1) - l^2 = l\]

\Delta L_x \cdot \Delta L_y = \frac{\hbar^2}{2} l = \frac{\hbar}{2} |\langle \hat{L}_z \rangle|

States with \(m = \pm l\) are minimum uncertainty states for angular momentum.

Advanced Example 3.6: Problem 6: Anharmonic Morse Potential & Vibrational Dissociation Limit

A real diatomic molecule is modeled by the Morse potential:

V(x) = D_e [1 - e^{-a x}]^2

where \(D_e\) is the equilibrium dissociation energy and \(x = R - R_e\). The quantized energy levels are given by:

E_v = \hbar\omega_e \left(v + \frac{1}{2}\right) - \hbar\omega_e x_e \left(v + \frac{1}{2}\right)^2

where the anharmonicity constant is \(\omega_e x_e = \frac{\hbar\omega_e^2}{4 D_e}\).

For \(\text{HCl}\), \(\tilde{\omega}_e = 2990.9\text{ cm}^{-1}\) and \(\tilde{\omega}_e x_e = 52.8\text{ cm}^{-1}\). Calculate the depth \(D_e\) of the Morse potential well (in \(\text{eV}\)).
Calculate the fundamental transition wavenumber \(\tilde{\nu}_1 = E_1 - E_0\) and the first overtone \(\tilde{\nu}_2 = E_2 - E_0\).
Determine the maximum vibrational quantum number \(v_{\text{max}}\) before dissociation occurs.

Comprehensive Multi-Step Solution:

Step 1: Equilibrium Dissociation Energy \(D_e\)

From Morse theory:

\tilde{\omega}_e x_e = \frac{\tilde{\omega}_e^2}{4 D_e} \implies D_e(\text{cm}^{-1}) = \frac{\tilde{\omega}_e^2}{4 \tilde{\omega}_e x_e}

Substitute \(\tilde{\omega}_e = 2990.9\text{ cm}^{-1}\) and \(\tilde{\omega}_e x_e = 52.8\text{ cm}^{-1}\):

D_e = \frac{(2990.9)^2}{4 \times 52.8} = \frac{8945482.8}{211.2} \approx 42355.5\text{ cm}^{-1}

Convert to electron-volts (\(1\text{ eV} = 8065.54\text{ cm}^{-1}\)):

D_e = \frac{42355.5\text{ cm}^{-1}}{8065.54\text{ cm}^{-1}/\text{eV}} \approx 5.2514\text{ eV} \approx 5.25\text{ eV}

In \(\text{kJ/mol}\):

D_e = 5.2514 \times 96.485\text{ kJ/mol} \approx 506.7\text{ kJ/mol}

Step 2: Fundamental and First Overtone Transitions

Using \(E_v = \tilde{\omega}_e (v + 1/2) - \tilde{\omega}_e x_e (v + 1/2)^2\):

1. Fundamental Transition (\(v = 0 \rightarrow 1\)):

\Delta E(0 \rightarrow 1) = \tilde{\omega}_e - 2 \tilde{\omega}_e x_e

\tilde{\nu}_1 = 2990.9 - 2(52.8) = 2990.9 - 105.6 = 2885.3\text{ cm}^{-1}

2. First Overtone Transition (\(v = 0 \rightarrow 2\)):

\Delta E(0 \rightarrow 2) = 2 \tilde{\omega}_e - 6 \tilde{\omega}_e x_e

\tilde{\nu}_2 = 2(2990.9) - 6(52.8) = 5981.8 - 316.8 = 5665.0\text{ cm}^{-1}

Notice that \(\tilde{\nu}_2 < 2 \tilde{\nu}_1\) (\(5665.0\text{ cm}^{-1} < 5770.6\text{ cm}^{-1}\)) due to anharmonic level crowding.

Step 3: Maximum Vibrational Quantum Number \(v_{\text{max}}\)

Dissociation occurs when the spacing between adjacent levels vanishes:

\Delta E_v = E_{v+1} - E_v = \tilde{\omega}_e - 2 \tilde{\omega}_e x_e (v + 1) = 0

v_{\text{max}} + 1 = \frac{\tilde{\omega}_e}{2 \tilde{\omega}_e x_e} = \frac{2990.9}{2 \times 52.8} = \frac{2990.9}{105.6} \approx 28.32

v_{\text{max}} = 28.32 - 1 = 27.32 \implies v_{\text{max}} = 27

The \(\text{HCl}\) potential well supports 28 bound vibrational states (\(v = 0, 1, \dots, 27\)) before reaching the continuous dissociation continuum.

Advanced Example 3.7: Problem 7: Centrifugal Distortion & Non-Rigid Rotor Correction

A real diatomic molecule stretches as it rotates due to centrifugal force. The effective term values including the quartic centrifugal distortion constant \(D_J\) are:

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

where \(D_J \approx \frac{4 B^3}{\tilde{\omega}_e^2}\).

For \(^{12}\text{C}^{16}\text{O}\), \(B = 1.93128\text{ cm}^{-1}\) and \(\tilde{\omega}_e = 2169.8\text{ cm}^{-1}\). Calculate \(D_J\) in \(\text{cm}^{-1}\).
Calculate the transition wavenumber for \(J = 9 \rightarrow 10\) with and without the centrifugal distortion correction.
Determine the percentage error incurred if the rigid rotor approximation is used for this high-\(J\) transition.

Comprehensive Multi-Step Solution:

Step 1: Centrifugal Distortion Constant \(D_J\)

Using the Kratzer relation:

D_J = \frac{4 B^3}{\tilde{\omega}_e^2}

Substitute \(B = 1.93128\text{ cm}^{-1}\) and \(\tilde{\omega}_e = 2169.8\text{ cm}^{-1}\):

B^3 = (1.93128)^3 \approx 7.2033\text{ cm}^{-3}

\tilde{\omega}_e^2 = (2169.8)^2 \approx 4708032\text{ cm}^{-2}

D_J = \frac{4 \times 7.2033}{4708032} = \frac{28.8132}{4708032} \approx 6.120 \times 10^{-6}\text{ cm}^{-1}

Step 2: Transition Wavenumber for \(J = 9 \rightarrow 10\)

For non-rigid rotor:

\tilde{\nu}(J \rightarrow J + 1) = F(J + 1) - F(J) = 2 B (J + 1) - 4 D_J (J + 1)^3

For \(J = 9 \implies J + 1 = 10\):

1. Rigid Rotor Approximation (\(D_J = 0\)):

\tilde{\nu}_{\text{rigid}} = 2 B (10) = 20 B = 20 \times 1.93128\text{ cm}^{-1} = 38.6256\text{ cm}^{-1}

2. Centrifugal Correction:

\Delta\tilde{\nu}_{\text{centrif}} = -4 D_J (10)^3 = -4000 D_J = -4000 \times (6.120 \times 10^{-6}\text{ cm}^{-1}) = -0.02448\text{ cm}^{-1}

3. Corrected Transition Wavenumber:

\tilde{\nu}_{\text{actual}} = 38.6256 - 0.02448 = 38.6011\text{ cm}^{-1}

Step 3: Error Analysis

The absolute difference is \(0.0245\text{ cm}^{-1}\). The relative percentage error is:

\% \text{ Error} = \frac{0.02448\text{ cm}^{-1}}{38.6011\text{ cm}^{-1}} \times 100\% \approx 0.0634\%

In high-resolution microwave and millimeter-wave Fourier-transform spectroscopy (measurement precision \(\sim 10^{-4}\text{ cm}^{-1}\)), an error of \(0.0245\text{ cm}^{-1}\) is massive (over 200 times the instrumental resolution), demonstrating that centrifugal distortion corrections are essential for high-accuracy molecular structure determination.

Advanced Example 3.8: Rovibrational P-Branch and R-Branch Transitions for HCl

For the fundamental infrared band (\(v = 0 \rightarrow 1\)) of \(^1\text{H}^{35}\text{Cl}\): The spectroscopic constants are:

Band origin: \(\tilde{\nu}_0 = 2885.98\text{ cm}^{-1}\)
Ground state rotational constant: \(B_0 = 10.440\text{ cm}^{-1}\)
Excited state rotational constant: \(B_1 = 10.137\text{ cm}^{-1}\)

Write the explicit formulas for the transition wavenumbers of the R-branch \(\tilde{\nu}_R(J'')\) and P-branch \(\tilde{\nu}_P(J'')\) ignoring centrifugal distortion.
Calculate the wavenumbers for the first three lines of the R-branch (\(R(0), R(1), R(2)\)) and P-branch (\(P(1), P(2), P(3)\)).
Calculate the spacing between adjacent lines in the R and P branches and explain why lines converge in the R-branch while spreading in the P-branch.

Comprehensive Multi-Step Solution:

Step 1: Formulas for Rovibrational Transitions

The term values are \(T(v, J) = G(v) + B_v J(J + 1)\). The transition wavenumbers from ground state \(J''\) are:

1. R-branch (\(J' = J'' + 1\)):

\tilde{\nu}_R(J'') = \tilde{\nu}_0 + B_1(J'' + 1)(J'' + 2) - B_0 J''(J'' + 1)

Expanding:

\tilde{\nu}_R(J'') = \tilde{\nu}_0 + 2 B_1 + (3 B_1 - B_0) J'' + (B_1 - B_0) J''^2

2. P-branch (\(J' = J'' - 1\)):

\tilde{\nu}_P(J'') = \tilde{\nu}_0 + B_1(J'' - 1) J'' - B_0 J''(J'' + 1)

Expanding:

\tilde{\nu}_P(J'') = \tilde{\nu}_0 - (B_1 + B_0) J'' + (B_1 - B_0) J''^2

Step 2: Numerical Calculation of Transition Wavenumbers

Given \(\tilde{\nu}_0 = 2885.98\text{ cm}^{-1}\), \(B_0 = 10.440\text{ cm}^{-1}\), \(B_1 = 10.137\text{ cm}^{-1}\):

\(B_1 - B_0 = 10.137 - 10.440 = -0.303\text{ cm}^{-1}\)
\(3 B_1 - B_0 = 3(10.137) - 10.440 = 30.411 - 10.440 = 19.971\text{ cm}^{-1}\)
\(2 B_1 = 20.274\text{ cm}^{-1}\)
\(B_1 + B_0 = 10.137 + 10.440 = 20.577\text{ cm}^{-1}\)

1. R-Branch Lines:

\(R(0)\) (\(J'' = 0\)):

\tilde{\nu}_R(0) = 2885.98 + 20.274 = 2906.25\text{ cm}^{-1}

\(R(1)\) (\(J'' = 1\)):

\tilde{\nu}_R(1) = 2885.98 + 20.274 + 19.971(1) - 0.303(1)^2 = 2906.254 + 19.668 = 2925.92\text{ cm}^{-1}

\(R(2)\) (\(J'' = 2\)):

\tilde{\nu}_R(2) = 2885.98 + 20.274 + 19.971(2) - 0.303(4) = 2906.254 + 39.942 - 1.212 = 2944.98\text{ cm}^{-1}

2. P-Branch Lines:

\(P(1)\) (\(J'' = 1\)):

\tilde{\nu}_P(1) = 2885.98 - 20.577(1) - 0.303(1)^2 = 2885.98 - 20.880 = 2865.10\text{ cm}^{-1}

\(P(2)\) (\(J'' = 2\)):

\tilde{\nu}_P(2) = 2885.98 - 20.577(2) - 0.303(4) = 2885.98 - 41.154 - 1.212 = 2843.61\text{ cm}^{-1}

\(P(3)\) (\(J'' = 3\)):

\tilde{\nu}_P(3) = 2885.98 - 20.577(3) - 0.303(9) = 2885.98 - 61.731 - 2.727 = 2821.52\text{ cm}^{-1}

Step 3: Spacing and Asymmetry Analysis

Spacings in R-branch:
\(R(1) - R(0) = 2925.92 - 2906.25 = 19.67\text{ cm}^{-1}\)
\(R(2) - R(1) = 2944.98 - 2925.92 = 19.06\text{ cm}^{-1}\)

The spacing decreases by \(\approx 0.61\text{ cm}^{-1} \approx 2(B_0 - B_1)\).

Spacings in P-branch:
\(P(1) - P(2) = 2865.10 - 2843.61 = 21.49\text{ cm}^{-1}\)
\(P(2) - P(3) = 2843.61 - 2821.52 = 22.09\text{ cm}^{-1}\)

The spacing increases by \(\approx 0.60\text{ cm}^{-1} \approx 2(B_0 - B_1)\). Physical Reason: In higher vibrational states (\(v = 1\)), anharmonic stretching increases the average bond length \(\langle R \rangle\), which increases the moment of inertia \(I\) and decreases the rotational constant (\(B_1 < B_0\)). Because \((B_1 - B_0) < 0\), the quadratic correction shifts all transitions toward lower wavenumbers, compressing lines in the R-branch while spreading them apart in the P-branch.

Advanced Example 3.9: Spherical Harmonics and Rigid Rotor Transition Dipole Moment

Consider a rigid rotor molecule with permanent electric dipole moment \(\mu_0\) oriented along its internuclear axis \(\mathbf{r}/r\).

Write the electric dipole moment component \(\hat{\mu}_z = \mu_0 \cos\theta\) in terms of the Spherical Harmonic \(Y_1^0(\theta, \phi)\).
Evaluate the transition dipole matrix element \(\langle J', M' | \hat{\mu}_z | J, M \rangle\) between rotational states \(|J, M\rangle\) and \(|J', M'\rangle\).
Prove mathematically that rotational transitions require \(\Delta J = \pm 1\) and \(\Delta M = 0\) for \(z\)-polarized radiation.

Comprehensive Multi-Step Solution:

Step 1: Dipole Moment in Spherical Harmonics

The \(z\)-component of the dipole moment is:

\hat{\mu}_z = \mu_0 \cos\theta

Recall that the Spherical Harmonic for \(l = 1, m = 0\) is:

Y_1^0(\theta, \phi) = \sqrt{\frac{3}{4\pi}} \cos\theta \implies \cos\theta = \sqrt{\frac{4\pi}{3}} Y_1^0(\theta, \phi)

Therefore:

\hat{\mu}_z = \mu_0 \sqrt{\frac{4\pi}{3}} Y_1^0(\theta, \phi)

Step 2: Evaluation of the Transition Dipole Matrix Element

The transition dipole integral is:

\langle J', M' | \hat{\mu}_z | J, M \rangle = \mu_0 \sqrt{\frac{4\pi}{3}} \int_0^{2\pi} d\phi \int_0^\pi \sin\theta \, d\theta \, [Y_{J'}^{M'}(\theta, \phi)]^* Y_1^0(\theta, \phi) Y_J^M(\theta, \phi)

The integral over three Spherical Harmonics evaluates via Clebsch-Gordan coefficients / Wigner \(3j\)-symbols:

\int Y_{J'}^{M'*} Y_{J_2}^{M_2} Y_J^M d\Omega = (-1)^{M'} \sqrt{\frac{(2J'+1)(2J_2+1)(2J+1)}{4\pi}} \begin{pmatrix} J' & J_2 & J \\ 0 & 0 & 0 \end{pmatrix} \begin{pmatrix} J' & J_2 & J \\ -M' & M_2 & M \end{pmatrix}

Here \(J_2 = 1\) and \(M_2 = 0\):

1. Selection rule for \(M\):

The second \(3j\)-symbol requires \(-M' + 0 + M = 0 \implies M' = M \implies \Delta M = 0\).

2. Selection rule for \(J\):

By angular momentum addition triangle rules for coupling \(J\) and \(1\):

|J - 1| \le J' \le J + 1 \implies J' \in \{J - 1, J, J + 1\}

Furthermore, the parity condition embodied in \(\begin{pmatrix} J' & 1 & J \\ 0 & 0 & 0 \end{pmatrix}\) requires that \(J' + 1 + J\) must be an even integer.

If \(J' = J\): \(J + 1 + J = 2J + 1\) is odd, so the \(3j\)-symbol is strictly zero!
Therefore, \(J' = J\) is forbidden: \(\Delta J \ne 0\).
Only \(J' = J + 1\) and \(J' = J - 1\) survive.

Step 3: Selection Rule Proof

The explicit analytical result for the transition \(J \rightarrow J + 1\) with \(\Delta M = 0\) is:

|\langle J + 1, M | \hat{\mu}_z | J, M \rangle|^2 = \mu_0^2 \frac{(J + 1)^2 - M^2}{(2J + 1)(2J + 3)}

When \(J' = J \pm 1\), the transition dipole moment is non-zero (allowed electric dipole transition).
When \(J' = J\) or \(|J' - J| > 1\), the integral vanishes identically by parity and angular momentum orthogonality.

Thus, pure rotational absorption in linear molecules requires:

\Delta J = \pm 1, \quad \Delta M = 0 \quad (\text{for } z\text{-polarized radiation})

and \(\Delta M = \pm 1\) for circularly polarized radiation. This fundamental selection rule gives rise to the classic microwave spectra with lines evenly spaced by \(2B\).