Case Study 1: Molecular Vibrations — The Quantum Harmonic Oscillator in Chemistry

The Physical Situation

Every molecule vibrates. Even at absolute zero temperature, the atoms in a molecule oscillate around their equilibrium positions with zero-point motion that can never be eliminated. Understanding these vibrations is not merely an academic exercise — it is the foundation of infrared spectroscopy, one of the most powerful tools in analytical chemistry, atmospheric science, and materials characterization.

Consider the simplest case: a diatomic molecule such as hydrogen chloride (HCl). The two atoms are bound by a shared electron cloud that acts as an effective spring. When the bond is stretched or compressed from its equilibrium length $r_0$, the potential energy rises. For small displacements, this potential is parabolic — a harmonic oscillator.

Setting Up the Model

The internuclear potential of a diatomic molecule is well described by the Morse potential:

$$V(r) = D_e \left(1 - e^{-a(r - r_0)}\right)^2$$

where $D_e$ is the well depth (dissociation energy from the equilibrium position), $a$ controls the width, and $r_0$ is the equilibrium bond length.

Expanding around $r = r_0$ to second order:

$$V(r) \approx D_e a^2 (r - r_0)^2 + O((r - r_0)^3)$$

This is a harmonic oscillator with effective spring constant $k = 2D_e a^2$. The vibrational frequency is:

$$\omega = \sqrt{\frac{k}{\mu}} = a\sqrt{\frac{2D_e}{\mu}}$$

where $\mu = m_1 m_2/(m_1 + m_2)$ is the reduced mass of the two-atom system.

The quantum energy levels (in the harmonic approximation) are:

$$E_v = \left(v + \frac{1}{2}\right)\hbar\omega, \qquad v = 0, 1, 2, \ldots$$

Experimental Data: HCl as a Case Study

For ${}^1$H${}^{35}$Cl:

The harmonic model predicts equally spaced energy levels with spacing $\hbar\omega_e \approx 0.358$ eV. Let us compare with experiment.

Observed Vibrational Transitions

The harmonic model does well for the fundamental transition but becomes progressively worse for higher transitions. The observed levels get closer together as $v$ increases — the molecule is anharmonic.

Beyond the Harmonic Approximation

The full Morse potential gives the exact energy levels:

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

where $x_e = \hbar\omega_e / (4D_e)$ is the anharmonicity constant. For HCl, $x_e \approx 0.0174$. The anharmonic correction:

• Reduces the spacing between levels as $v$ increases (levels converge toward $D_e$).
• Limits the number of bound states to a finite value: $v_{\max} \approx 1/(2x_e) - 1/2$.
• Allows overtone transitions ($\Delta v = 2, 3, \ldots$), which are forbidden in the harmonic model but weakly allowed in the anharmonic case.

For HCl with the Morse potential, the anharmonic predictions match experiment to better than 0.1% for all observed transitions.

Isotope Effects: The Mass Dependence of Vibrations

Since $\omega = \sqrt{k/\mu}$ and the electronic potential (hence $k$) is essentially the same for different isotopes of the same element, the vibrational frequency scales as $\omega \propto 1/\sqrt{\mu}$.

Example: HCl vs. DCl

Replacing hydrogen with deuterium roughly doubles the reduced mass: $\mu_{\text{DCl}} \approx 1.904\mu_{\text{HCl}}$ (not exactly double because Cl is much heavier). Therefore:

$$\frac{\omega_{\text{DCl}}}{\omega_{\text{HCl}}} = \sqrt{\frac{\mu_{\text{HCl}}}{\mu_{\text{DCl}}}} \approx \frac{1}{\sqrt{1.904}} \approx 0.724$$

The predicted DCl fundamental frequency is $0.724 \times 2886 = 2090$ cm$^{-1}$. The observed value is 2091 cm$^{-1}$ — a triumph of the model.

Zero-point energy difference: The ZPE of HCl is $\frac{1}{2}\hbar\omega_{\text{HCl}} \approx 0.179$ eV, while for DCl it is $\frac{1}{2}\hbar\omega_{\text{DCl}} \approx 0.130$ eV. The difference (0.049 eV) means:

• The D-Cl bond is slightly more stable than the H-Cl bond.
• The DCl molecule has a marginally shorter effective bond length.
• Reactions involving bond breaking are slower for deuterium (the kinetic isotope effect), which is crucial in biochemistry and pharmaceutical research.

The Greenhouse Effect Connection

Molecules absorb infrared radiation at their vibrational frequencies. This is the molecular basis of the greenhouse effect:

1. Solar radiation (mostly visible/UV) passes through the atmosphere.
2. The Earth's surface absorbs and re-emits as infrared (thermal) radiation.
3. Greenhouse gases (CO$_2$, H$_2$O, CH$_4$) absorb this IR radiation at their vibrational frequencies.
4. The absorbed energy is re-emitted in all directions, trapping heat.

The CO$_2$ molecule has a bending mode at 667 cm$^{-1}$ (15 $\mu$m) that falls exactly in the atmospheric IR window. This single vibrational transition is the primary mechanism of CO$_2$-driven climate change. Understanding it requires the quantum harmonic oscillator.

Notably, symmetric diatomic molecules (N$_2$, O$_2$) have no dipole moment and therefore do not absorb IR radiation — they are transparent to thermal radiation and are not greenhouse gases. The QHO selection rule ($\Delta v = \pm 1$ with a changing dipole moment) explains this directly.

Discussion Questions

1. Why does the harmonic approximation work better for low vibrational quantum numbers than for high ones? Relate your answer to the Taylor expansion argument of Section 4.1.

2. The kinetic isotope effect (slower reaction rates for heavier isotopes) is used in pharmaceutical design to create drugs with longer half-lives ("deuterated drugs"). Explain the physics behind this application in terms of zero-point energy and bond dissociation.

3. The QHO model predicts that CO$_2$ absorbs at specific discrete frequencies. But the actual absorption spectrum shows broad bands, not sharp lines. What additional physics produces the broadening? (Consider rotation, collisions, and temperature effects.)

4. If you could measure the vibrational frequency of a diatomic molecule to arbitrary precision, could you determine both the force constant $k$ and the reduced mass $\mu$ independently? What additional measurement would you need?

Computational Exploration

Using the code from code/example-01-oscillator.py:

1. Plot the first six energy levels and wavefunctions for HCl.
2. Overlay the Morse potential with the harmonic approximation and observe where they diverge.
3. Calculate and plot the probability density $|\psi_v(r)|^2$ for $v = 0, 5, 10$. At what point does the classically forbidden penetration become negligible?
4. Compute the isotope shift for all hydrogen halides (HF, HCl, HBr, HI) upon deuterium substitution and plot $\omega_{\text{DX}}/\omega_{\text{HX}}$ vs. reduced mass ratio.