The Harmonic Oscillator as a Model for Vibration

Every infrared spectrum recorded in an organic chemistry laboratory is interpreted, at the first level of approximation, through this exact model: a carbonyl stretch near 1700 cm⁻¹, an O–H stretch near 3300 cm⁻¹, a C–H bend near 1450 cm⁻¹ — each assigned to a characteristic $k$ and $\mu$ combination via $\omega = \sqrt{k/\mu}$, precisely as computed in the worked examples above. Force-field molecular mechanics, used throughout computational chemistry and drug design to model protein folding and small-molecule conformations cheaply, represents essentially every bond and angle as a harmonic (or near-harmonic) spring, with force constants fitted to reproduce experimental vibrational frequencies or high-level quantum chemical calculations. In materials science, the same parabolic-well idea underlies the Einstein and Debye models of specific heat in solids, where atoms vibrating about lattice sites are treated as independent (or coupled) harmonic oscillators to predict how heat capacity depends on temperature. In astrochemistry, the vibrational partition function built from harmonic oscillator energy levels feeds directly into equilibrium constants and reaction rate predictions for interstellar molecules, where direct experiments are impossible and spectroscopic identification relies entirely on matching predicted vibrational frequencies to observed emission lines. Modern vibrational instrumentation — Fourier-transform infrared (FTIR) spectrometers, Raman microscopes used in pharmaceutical quality control, and even the vibrational sum-frequency generation spectroscopy used to probe interfacial water structure — all report data whose first-pass interpretation is a harmonic force constant and a reduced mass. Even as increasingly sophisticated anharmonic corrections (Morse potentials, cubic and quartic force fields, explicit oscillator wavefunctions) refine the picture for precision spectroscopy, the parabola from a Taylor expansion remains the conceptual entry point for essentially every vibrational calculation done in chemistry.

Ask the tutor: How do the exact quantized energy levels and wavefunctions emerge once this parabolic Hamiltonian is actually solved?

Picture a diatomic molecule, say carbon monoxide, with its two nuclei connected by an electronic bonding potential $V(r)$ that depends on the internuclear separation $r$. This potential has a characteristic shape: it plunges steeply at short range (nuclear and electronic repulsion), passes through a minimum at the equilibrium bond length $r_e$, and climbs gently toward a flat asymptote at large $r$ where the bond has fully dissociated. This is the real, anharmonic, Morse-like curve that governs the true physics of the bond. Now Taylor expand $V(r)$ around $r_e$. Write $x = r - r_e$ as the displacement from equilibrium. Then

Two simplifications happen immediately. First, $V(0)$ is just a constant offset in energy — it can be set to zero without changing any physics, since only energy differences matter. Second, because $r_e$ is defined as the minimum of the potential, the first derivative $V'(0)$ vanishes by definition of a minimum. What remains, to leading nontrivial order, is

where $k \equiv V''(0)$ is the curvature of the potential at the bottom of the well — the force constant. This is the entire justification for the harmonic oscillator model: it is not a guess about bond physics, it is the unavoidable leading behavior of any smooth potential minimum, truncated at the first nonvanishing, non-constant term. The cubic and higher terms that were dropped are the anharmonic corrections. They are small when $x$ is small — that is, for a bond near its equilibrium length, executing modest-amplitude vibrations — and they become significant only as the bond stretches far enough to "feel" the true asymmetry of the potential (soft on the dissociation side, hard on the compression side). For the low vibrational states of most bonds at room temperature, the parabolic approximation is excellent. Once the potential is parabolic, the physical picture is exactly the mass-on-a-spring problem familiar from mechanics, but now applied to two nuclei connected by a bond rather than a block on a table. The two-body problem of nuclei with masses $m_1$ and $m_2$ moving under a mutual potential $V(x)$ separates cleanly into center-of-mass motion (translation of the whole molecule, uninteresting for vibration) and relative motion, described by a single fictitious particle of reduced mass

moving in the potential $V(x)$. This reduction is exact for any two-body central potential, not an extra approximation layered on top of the harmonic assumption. The vibrational problem is therefore, in full generality, "one particle of mass $\mu$ in the potential $\tfrac{1}{2}kx^2$." Classically, Newton's second law for this particle reads $\mu\ddot x = -kx$, whose solution is sinusoidal motion at angular frequency $\omega = \sqrt{k/\mu}$, independent of the amplitude — the hallmark of a harmonic (linear-restoring-force) system. Quantum mechanically, the same $k$ and $\mu$ enter the Hamiltonian operator, but the allowed energies turn out to be quantized rather than continuous, and the particle can never have exactly zero energy, even at absolute zero. Building that Hamiltonian rigorously, and understanding precisely how classical and quantum pictures relate, is the heart of the sections that follow.

The classical starting point is Hooke's law, the linear restoring force that defines a harmonic spring:

where:

• $F(x)$ is the restoring force on the particle at displacement $x$ from equilibrium,
• $k$ is the force constant (units of energy/length², e.g. N/m or mdyn/Å),
• $x = r - r_e$ is displacement from the equilibrium bond length.

Reading it: force is always directed back toward equilibrium and grows linearly with displacement — stretch the bond twice as far, and it pulls back twice as hard. This is exactly the leading behavior extracted from the Taylor expansion above. The classical equation of motion and its frequency follow from Newton's second law:

where:

• $\mu$ is the reduced mass of the two-nucleus system,
• $\omega$ is the classical angular frequency of oscillation (rad/s), related to the spectroscopic frequency by $\omega = 2\pi\nu$.

Reading it: stiffer bonds (larger $k$) vibrate faster; heavier reduced masses vibrate slower. This single relation is the reason C–H stretches (light $\mu$) appear at much higher wavenumber than C–Cl stretches (heavy $\mu$, comparable $k$). The quantum mechanical Hamiltonian is built by the standard correspondence: take the classical energy (kinetic plus potential), and replace momentum by its operator.

where:

• $\hat H$ is the vibrational Hamiltonian operator acting on wavefunctions $\psi(x)$,
• $-\dfrac{\hbar^2}{2\mu}\dfrac{d^2}{dx^2}$ is the kinetic energy operator, obtained from $\hat p^2/2\mu$ with $\hat p = -i\hbar\, d/dx$,
• $\tfrac{1}{2}kx^2$ is the parabolic potential energy operator (multiplicative, since $V$ depends only on position).

Reading it: this is the time-independent Schrödinger equation $\hat H \psi = E\psi$ for a single quantum particle confined not by hard walls (as in the particle in a box) but by a smoothly rising quadratic wall on both sides. The absence of infinite walls, and the smooth quadratic rise instead, is what will generate an infinite ladder of evenly spaced levels rather than the box's $n^2$-spaced levels — a qualitative difference worth anticipating even before solving the equation explicitly. It is often convenient to rewrite the Hamiltonian using $\omega$ rather than $k$, since $k = \mu\omega^2$:

where the substitution makes explicit that only one classical parameter, $\omega$, together with $\mu$ and $\hbar$, is needed to fully specify the quantum problem — a simplification that becomes essential when the exact energy levels and Hermite-polynomial wavefunctions are constructed.

Example 1 — Foundational: force constant and frequency of HCl. The H–Cl bond has an experimentally measured vibrational wavenumber of $\tilde\nu = 2886\ \text{cm}^{-1}$. Take atomic masses $m_H = 1.008$ u and $m_{Cl} = 34.97$ u (chlorine-35, the dominant isotope). Find the reduced mass and the force constant. First, the reduced mass:

Converting to kilograms using $1\ \text{u} = 1.6605\times10^{-27}\ \text{kg}$:

Convert the wavenumber to angular frequency. Since $\tilde\nu = 1/\lambda$ and $E = h c\tilde\nu = \hbar\omega$,

Now solve $\omega = \sqrt{k/\mu}$ for $k$:

This value, roughly $4.8\ \text{mdyn/Å}$ in the older spectroscopic unit, is a genuinely stiff bond — comparable in stiffness to many single bonds — and matches tabulated literature force constants for HCl closely. As a limiting-case sanity check, note that if the bond were far weaker, say $k \to 0$, then $\omega \to 0$ and the vibrational frequency would collapse to zero — an unphysically "floppy" bond, which is not what HCl's steep, well-defined potential minimum implies, consistent with the large $k$ actually obtained.

Example 2 — Applied: isotope shift from H to D in HCl. Predict the vibrational wavenumber of DCl (deuterium chloride) given the HCl force constant found above, assuming $k$ is unchanged by isotopic substitution (a very good approximation, since $k$ depends on the electronic structure, not on nuclear mass). The reduced mass of DCl, with $m_D = 2.014$ u:

The new angular frequency, using the same $k = 480.3\ \text{N/m}$:

Converting back to wavenumber:

This predicted red-shift, from $2886\ \text{cm}^{-1}$ down to about $2069\ \text{cm}^{-1}$, is a hallmark test of the harmonic model, and it agrees closely with the experimental DCl fundamental (observed near $2145\ \text{cm}^{-1}$, the small residual gap being attributable to anharmonicity, which affects H and D slightly differently). The predicted ratio,

compares to the naive limit of $\sqrt{1/2} = 0.7071$ expected if chlorine were infinitely heavy (so $\mu \approx m_H$ or $m_D$ exactly) — the small deviation from $0.7071$ reflects chlorine's finite mass entering the reduced mass, a useful limiting case to keep in mind: as $m_{Cl}\to\infty$, $\mu \to m_H$ (or $m_D$) and the ratio would approach exactly $\sqrt{1/2}$.

Example 3 — Exam-level: force constant from a Morse-type expansion. A diatomic's potential is modeled (for illustrative purposes) by a Morse function,

with dissociation energy $D_e = 7.50\times10^{-19}\ \text{J}$ and range parameter $\alpha = 1.90\ \text{Å}^{-1}$. Derive the effective harmonic force constant this potential implies near $x=0$, and compare the harmonic and true cubic-corrected potential at $x = 0.15$ Å. First expand $e^{-\alpha x}$ for small $x$:

Squaring and keeping terms through $x^3$:

So

Comparing to $\tfrac12 kx^2$ identifies the harmonic force constant as

Plugging in numbers (convert $\alpha$ to SI: $1.90\ \text{Å}^{-1} = 1.90\times10^{10}\ \text{m}^{-1}$):

This is a stiffer bond than HCl in Example 1 — consistent with a moderately deep, narrow Morse well. Now compare the harmonic prediction to the true Morse potential (with its leading cubic correction) at $x = 0.15\ \text{Å} = 1.50\times10^{-11}\ \text{m}$. The harmonic term:

The cubic correction term:

Computing $\alpha^3 = (1.90\times10^{10})^3 = 6.859\times10^{30}\ \text{m}^{-3}$ and $x^3 = (1.50\times10^{-11})^3 = 3.375\times10^{-33}\ \text{m}^3$:

The cubic term is about $-28.5\%$ of the quadratic term at this displacement — a substantial anharmonic correction, showing that even at a fairly modest stretch of 0.15 Å, the harmonic approximation is already beginning to noticeably overestimate the true restoring energy (since the true potential softens on the stretching side). This is precisely the expected trend: the harmonic model is trustworthy near $x=0$ (where, as $x\to0$, the cubic-to-quadratic ratio $\to 0$ exactly), but its reliability degrades as $|x|$ grows, always in the direction of the true potential being softer (lower energy) than the parabola predicts on the outward side, and stiffer on the inward side.

The mathematics of the harmonic oscillator predates quantum chemistry by centuries — Robert Hooke described the linear restoring force in 1676, and the classical mechanics of springs and pendulums was fully worked out by the eighteenth century. What is distinctly a twentieth-century achievement is recognizing that quantizing this exact same potential, inside the new Schrödinger equation, gives a model that maps almost perfectly onto real molecular vibrations. That synthesis belongs to the founders of quantum mechanics applying their new wave equation to molecular problems in the mid-1920s, but it is worth naming the figure whose name is permanently attached to the restoring-force law at the heart of the model.