Vibrational Energy Levels

Historical Context & Motivation

The study of molecular vibrations sits at the intersection of classical mechanics and quantum theory, a convergence that transformed how chemists understand the internal dynamics of molecules. Long before spectroscopists could assign individual absorption bands to specific bond stretches and bends, physicists grappled with a deceptively simple question: why do molecules absorb light only at discrete frequencies rather than across a continuous range? The resolution of this question drew upon the same quantum principles that Planck and Einstein had invoked to explain blackbody radiation and the photoelectric effect, and it ultimately gave rise to one of the most powerful analytical tools in modern chemistry—infrared (IR) spectroscopy. Understanding the historical pathway from classical harmonic motion to the quantum harmonic oscillator illuminates why vibrational energy levels are quantized and why that quantization has such far-reaching experimental consequences.

The central question that this lesson addresses is: How does the quantization of vibrational energy arise from the potential energy surface of a chemical bond, and what observable spectral features result? Answering this question requires us to move from the classical picture of a vibrating spring to the quantum-mechanical harmonic and anharmonic oscillator models, and then to connect those models to the infrared absorption spectra that chemists measure in the laboratory.

Core Principles & Definitions

Before diving into the mathematics, it is essential to establish several foundational concepts. A diatomic molecule can be modeled as two masses connected by a spring-like restoring force; this is the harmonic oscillator approximation. When the displacement from the equilibrium bond length is small, the potential energy is parabolic and the restoring force obeys Hooke's law. The quantum-mechanical treatment of this system yields equally spaced energy levels characterized by a vibrational quantum number v. Real molecules, however, deviate from this ideal behavior because the true potential energy curve is anharmonic: the repulsive wall at short internuclear distances is steeper than the attractive tail at long distances, and the molecule can ultimately dissociate. The Morse potential captures this asymmetry and predicts energy levels that converge as v increases, eventually reaching the dissociation limit.

Visualizing Vibrational Energy Levels

The most instructive way to understand vibrational quantization is to compare the harmonic and anharmonic (Morse) potential energy curves side by side. In the diagram below, the parabolic harmonic potential is shown alongside the Morse potential, with horizontal lines indicating the allowed energy levels for each model. Note how the harmonic oscillator predicts uniformly spaced levels, whereas the Morse oscillator shows levels that become progressively closer together as the vibrational quantum number increases, converging toward the dissociation energy Dₑ.

Several features of this diagram deserve careful attention. First, the zero-point energy is visible as the lowest horizontal line sitting above the bottom of each potential well—confirming that the molecule can never reside at the very minimum of the curve. Second, the harmonic oscillator's equal spacing (ΔE = hν₀ for every transition) stands in contrast with the Morse oscillator, where the spacing between adjacent levels decreases roughly linearly with v. Third, the Morse potential possesses a finite number of bound states, whereas the harmonic oscillator supports infinitely many. These distinctions have direct experimental consequences: the harmonic model predicts a single absorption frequency (the fundamental), while the anharmonic model predicts weaker overtone bands at frequencies slightly less than 2ν₀, 3ν₀, and so forth, exactly as observed in high-resolution IR spectra.

Mathematical Framework

The quantum-mechanical treatment of vibrational motion begins with the time-independent Schrödinger equation for a one-dimensional harmonic oscillator. We express the potential energy of a diatomic molecule in terms of the displacement coordinate x = r − rₑ, where rₑ is the equilibrium bond length. The Hamiltonian operator takes the familiar form Ĥ = −(ħ²/2μ)(d²/dx²) + ½kx², and solution by the power-series or algebraic (ladder-operator) method yields the following energy eigenvalues.

In spectroscopy, energies are conventionally expressed in wavenumber units (cm⁻¹) by dividing by hc. The vibrational term value is therefore written as follows.

To account for anharmonicity, the Morse oscillator introduces correction terms. The energy expression expands as a power series in (v + ½), where the first anharmonicity constant ω̃ₑxₑ causes the level spacing to decrease with increasing v.

Vibrational Transitions & Spectral Features

The observable IR spectrum of a molecule is governed by two complementary factors: the energies of the allowed transitions and the intensities (transition dipole moments) of those transitions. The gross selection rule requires that the electric dipole moment of the molecule change during the vibration for it to be IR-active. Homonuclear diatomics such as N₂ and O₂ have zero dipole moment at all bond lengths and therefore show no IR absorption, whereas heteronuclear diatomics like HCl and CO are strongly IR-active. The specific selection rule for the harmonic oscillator restricts transitions to Δv = ±1. Anharmonicity weakens this restriction, permitting overtone transitions (Δv = ±2, ±3, …) and combination bands in polyatomics, though these appear at progressively lower intensities.

The right panel of the diagram schematically illustrates the appearance of a typical IR spectrum. The fundamental band near ω̃ₑ dominates the spectrum because the transition dipole matrix element for Δv = 1 is largest. The first overtone near 2ω̃ₑ appears at roughly 1–5% of the fundamental's intensity, and the second overtone near 3ω̃ₑ is weaker still. At elevated temperatures, hot bands emerge because a non-negligible fraction of molecules populate v = 1 (governed by the Boltzmann distribution). These bands appear at slightly lower wavenumbers than the fundamental because ΔG(3/2) < ΔG(1/2) due to anharmonicity, and their intensity grows with temperature.

Worked Example: HCl Fundamental & Overtone

Let us apply the anharmonic oscillator model to hydrogen chloride (¹H³⁵Cl), one of the most well-studied diatomic molecules. The spectroscopic constants for the ground electronic state are ω̃ₑ = 2990.95 cm⁻¹ and ω̃ₑxₑ = 52.82 cm⁻¹. We will calculate the positions of the fundamental and first overtone transitions, as well as the dissociation energy D₀.

Harmonic vs. Anharmonic Oscillator: Strengths & Limitations

Both the harmonic and anharmonic oscillator models are valuable, but their domains of applicability differ substantially. The harmonic oscillator is analytically tractable, provides a first approximation to vibrational frequencies, and serves as the basis for normal-mode analysis in polyatomic molecules. However, it fails to predict overtone bands, does not account for thermal expansion of bonds, and implies infinite dissociation energy. The Morse oscillator corrects these shortcomings at the cost of introducing an empirical anharmonicity parameter and losing the clean algebraic structure of the harmonic solution.

Connection to Advanced Theory

The one-dimensional treatment of vibrational energy levels presented so far is only the beginning. In more advanced courses and research contexts, several extensions become necessary. For polyatomic molecules, the concept of normal modes provides a framework in which each of the 3N − 6 (or 3N − 5) independent vibrations is treated as a separate quantum harmonic oscillator, and the total vibrational energy is the sum of individual mode energies. Coupling between modes introduces Fermi resonance (accidental degeneracy between a fundamental and an overtone or combination band of similar symmetry) and Coriolis coupling (interaction between vibration and rotation). Computationally, ab initio methods like coupled-cluster theory (CCSD(T)) with large basis sets can predict harmonic frequencies to within 1–2% of experiment, while variational or perturbative VPT2 corrections add anharmonic contributions.

A particularly important extension is the inclusion of vibration–rotation coupling. In a gas-phase spectrum, each vibrational band is not a single line but a structured envelope of closely spaced rotational transitions. The combined energy expression E(v, J) = G(v) + F(J), where F(J) = B̃ₑJ(J+1) − D̃ₑ[J(J+1)]² and B̃ₑ depends weakly on v through the vibration–rotation interaction constant αₑ, produces the familiar P, Q, and R branch structure of rovibrational bands. Mastery of this coupled problem is the natural next step after understanding vibrational energy levels in isolation.

Practice Problems

Summary