The quantum harmonic oscillator (QHO) models molecular bond stretching and bending vibrations. Its energy levels are equally spaced: E_v = ℏω(v + 1/2), where v = 0, 1, 2, … The zero-point energy ℏω/2 persists even at absolute zero, reflecting the uncertainty principle. Wavefunctions are Hermite polynomials multiplied by a Gaussian envelope, and they extend into classically forbidden regions (tunneling). The anharmonic Morse potential is a more realistic model for real bonds, accounting for bond dissociation at high vibrational quantum numbers.
Verify the equally spaced energy levels first, then focus on the physical meaning of zero-point energy. Compare the QHO to the Morse oscillator to see how anharmonicity leads to overtone transitions in IR spectra.
The classical harmonic oscillator — a mass on a spring — has a continuous range of energies depending on how far you stretch it. Its quantum counterpart has something fundamentally different: energy comes in discrete packages. The allowed vibrational energies are E_v = ℏω(v + 1/2), where v = 0, 1, 2, … is the vibrational quantum number and ω = √(k/μ) is the angular frequency set by the force constant k and reduced mass μ. The energy levels are equally spaced by ℏω, much like the equally spaced levels you saw in the particle-in-a-box, but now the potential is curved rather than flat.
The term v + 1/2 rather than v contains a crucial message: even the lowest state (v = 0) has energy ℏω/2, not zero. This zero-point energy is a direct consequence of the uncertainty principle. A particle confined to a potential well cannot be simultaneously at rest at the minimum — that would require ΔxΔp = 0. Instead it must have some residual kinetic energy, which is the zero-point energy. This is not a curiosity: zero-point energy is physically real. It keeps helium liquid at atmospheric pressure even at 0 K, it contributes to isotope effects in chemical reactions (heavier isotopes have lower ω and lower zero-point energy, changing reaction rates), and it means molecular bonds are always vibrating.
The wavefunctions of the QHO are Hermite polynomials multiplied by a Gaussian, ψ_v(x) = N_v · H_v(αx) · e^(−α²x²/2). For v = 0, this is a simple Gaussian centered at equilibrium — the particle is most likely found near the center. Higher v states show more nodes and larger spatial spread. Critically, the wavefunctions extend into the classically forbidden regions beyond the classical turning points. This quantum tunneling has real spectroscopic consequences: it is partly why proton-transfer reactions can proceed faster than a classical analysis predicts.
In connecting this to molecular spectroscopy, the force constant k corresponds to the curvature of the potential energy surface at the bond's equilibrium length. Stiff bonds (like C≡C) have large k and high ω, absorbing IR light at high wavenumber. Weak bonds (like H-bonds) absorb at low wavenumber. Because the reduced mass μ also enters ω, isotopic substitution (e.g., H → D) shifts vibrational frequencies predictably — this is isotopic labeling, a powerful tool in structural chemistry.
The harmonic model is an approximation, valid only for small displacements. For large v, the real potential — better described by the Morse potential V(r) = D_e(1 − e^(−a(r−r_e)))² — deviates significantly. The Morse levels converge and ultimately the bond dissociates. This anharmonicity relaxes the strict Δv = ±1 selection rule, allowing overtone transitions (Δv = ±2, ±3) to appear with weaker intensity in IR spectra. The quantum harmonic oscillator is thus not just a model — it is the first rung of a ladder that leads to quantitative prediction of every molecular vibration in an IR or Raman spectrum.