Harmonic Oscillator Model - Physical Chemistry 2

SECTION 1

Historical Context & Motivation

The concept of oscillatory motion has been central to physics since the earliest studies of pendula and vibrating strings, but it was the crisis of blackbody radiation at the turn of the twentieth century that forced a quantum re-examination of the harmonic oscillator. Classical mechanics treated an oscillating particle as capable of possessing any continuous energy, yet this assumption led to the infamous ultraviolet catastrophe—the prediction that a blackbody would radiate infinite energy at short wavelengths. Max Planck resolved this paradox in 1900 by postulating that the energy of each oscillator in the cavity wall is restricted to discrete multiples of a fundamental quantum, hν. This single hypothesis launched quantum theory and placed the harmonic oscillator at the very heart of the new physics.

1900 • Planck's Quantum Hypothesis

Max Planck proposes that the energy of oscillators in a blackbody cavity is quantized in units of hν, resolving the ultraviolet catastrophe and introducing the quantum of action, h.

1926 • Schrödinger's Wave Equation

Erwin Schrödinger publishes his wave equation and solves it exactly for the harmonic oscillator potential, revealing equally spaced energy levels and Hermite-polynomial wave functions.

1927 • Dirac's Operator Method

Paul Dirac introduces the algebraic ladder-operator (creation and annihilation operator) approach, providing an elegant route to the harmonic oscillator spectrum without solving differential equations.

1929 • Morse Potential & Anharmonicity

Philip Morse proposes an anharmonic potential that accounts for bond dissociation, establishing the harmonic oscillator as the zeroth-order approximation to real molecular vibrations.

1960s • Coherent States & Quantum Optics

Roy Glauber and others develop coherent states of the harmonic oscillator to describe laser radiation, linking the quantum oscillator to modern photonics and quantum information science.

The central question that the quantum harmonic oscillator answers is deceptively simple: what are the allowed energies and spatial probability distributions of a particle confined by a restoring force proportional to its displacement? The answer reveals quantized energy levels, zero-point energy, and wave functions that extend into classically forbidden regions—features that have no counterpart in Newtonian mechanics and that directly manifest in infrared spectroscopy, lattice dynamics, and quantum field theory.

SECTION 2

Core Principles & Definitions

Before diving into the mathematics, it is essential to internalize the foundational principles that distinguish the quantum harmonic oscillator from its classical counterpart. The classical oscillator—a mass on a spring, for instance—can vibrate with any amplitude and therefore any energy. Quantum mechanics imposes a discretization on these energies and introduces phenomena such as zero-point energy and tunneling into classically forbidden regions. The following grid distills the four pillars upon which the entire model rests.

Parabolic Potential

The potential energy is V(x) = ½kx², where k is the force constant and x the displacement from equilibrium. This quadratic form is the leading-order Taylor expansion of any potential near a stable minimum.

Quantized Energy Levels

The allowed energies are E_v = (v + ½)ℏω, where v = 0, 1, 2, … is the vibrational quantum number. The levels are equally spaced by ℏω.

Zero-Point Energy

Even in the ground state (v = 0), the oscillator possesses a residual energy of ½ℏω. This non-vanishing energy is a direct consequence of the Heisenberg uncertainty principle: the particle cannot be simultaneously at rest and at the equilibrium position.

Hermite-Gaussian Wave Functions

The stationary-state wave functions are products of Hermite polynomials and a Gaussian envelope, ψ_v(x) = N_v H_v(αx) exp(−α²x²/2). They exhibit v nodes and non-negligible probability density beyond the classical turning points.

✦

KEY TAKEAWAY

Think of the quantum harmonic oscillator as a ladder whose rungs are evenly spaced by an amount ℏω. The lowest rung is not at the floor (zero energy) but half a rung above it—this is the zero-point energy. A molecule vibrating in the ground state is like a plucked guitar string that can never stop vibrating entirely: the uncertainty principle guarantees a minimum amplitude of motion, and that irreducible vibration carries real, measurable energy. Every rung you climb corresponds to absorbing one quantum of vibrational energy, which is precisely the photon energy seen in an infrared absorption spectrum.

SECTION 3

Visual Explanation — Energy Levels & Wave Functions

The parabolic potential V(x) = ½kx² is shown in purple. Horizontal lines mark the equally spaced energy levels E_v = (v + ½)ℏω for v = 0 through 4. Superimposed on the v = 0, 1, and 2 levels are schematic wave functions ψ_v(x), illustrating the increasing number of nodes and the spatial extension beyond the classical turning points.

The diagram above captures the defining features of the quantum harmonic oscillator in a single picture. The parabolic potential rises symmetrically about the equilibrium position x = 0, and the horizontal lines intersecting the parabola represent the equally spaced energy eigenvalues. Notice that the lowest level (v = 0) sits at ½ℏω, not at zero—this is the zero-point energy. The wave functions sketched on the first three levels reveal several important trends: each successive state gains one additional node (zero crossing), and the probability density increasingly concentrates near the classical turning points—the positions where the kinetic energy would be zero classically—even as the Gaussian tails extend well beyond them.

SECTION 4

Mathematical Framework

The quantum harmonic oscillator begins with writing the time-independent Schrödinger equation for a particle of mass m in a potential V(x) = ½kx². In what follows, we introduce the angular frequency ω = √(k/m) and the characteristic length α = (mω/ℏ)^1/2, which sets the spatial scale of the wave functions.

TIME-INDEPENDENT SCHRÖDINGER EQUATION

−(ℏ²/2m) d²ψ/dx² + ½kx²ψ = Eψ

ℏ = reduced Planck constant, m = particle mass, k = force constant, ψ = wave function, E = energy eigenvalue.

Introducing the dimensionless variable ξ = αx and defining ε = 2E/(ℏω) transforms the equation into the standard Hermite differential equation. Requiring that solutions remain normalizable (i.e., square-integrable) restricts ε to odd positive integers, yielding the celebrated energy eigenvalue formula.

ENERGY EIGENVALUES

E_v = (v + ½)ℏω, v = 0, 1, 2, …

v = vibrational quantum number, ω = √(k/m) = angular frequency. The uniform spacing ΔE = ℏω is the selection-rule transition energy in IR spectroscopy.

NORMALIZED WAVE FUNCTIONS

ψ_v(x) = (α/π^{1/2})^{1/2} (1/2^v v!)^{1/2} H_v(αx) exp(−α²x²/2)

H_v = Hermite polynomial of degree v, α = (mω/ℏ)^1/2. The Gaussian factor ensures rapid decay at large |x|, while the Hermite polynomial introduces v nodes.

LADDER OPERATORS (ALGEBRAIC METHOD)

â = √(mω/2ℏ)(x̂ + ip̂/mω), â† = √(mω/2ℏ)(x̂ − ip̂/mω)

â (annihilation) lowers the quantum number by one: â|v⟩ = √v |v−1⟩. â† (creation) raises it: â†|v⟩ = √(v+1) |v+1⟩. The Hamiltonian becomes Ĥ = ℏω(â†â + ½), making the spectrum derivable without solving differential equations.

💡

Why the Algebraic Approach Matters

The ladder-operator formalism is not merely an alternative derivation technique. It generalizes directly to quantum field theory, where photons are created and destroyed by identical operators acting on an oscillator vacuum state. Mastering this algebra here lays the groundwork for second quantization, which you will encounter in advanced quantum mechanics and statistical mechanics courses.

SECTION 5

Probability Densities & the Correspondence Principle

One of the most instructive aspects of the quantum harmonic oscillator is how its probability density |ψ_v(x)|² evolves as the quantum number v increases. In the ground state, the probability density is a simple Gaussian centered at the equilibrium position, implying that the particle is most likely to be found at x = 0. This is the opposite of the classical prediction, where the oscillating mass spends the least time at the center because it is moving fastest there. As v grows large, the quantum probability density develops many oscillations and its envelope increasingly resembles the classical time-averaged distribution, which peaks near the turning points. This convergence is a beautiful illustration of the correspondence principle—the requirement that quantum mechanics reproduce classical mechanics in the appropriate limit.

Top left: the ground-state (v = 0) probability density is a Gaussian peaked at the equilibrium position, with red dots marking the classical turning points x_tp. Top right: at v = 10, rapid oscillations appear and the envelope (dashed gold) approaches the classical U-shaped distribution. Bottom: direct comparison of the classical probability (dashed red, highest near turning points) and the quantum ground-state Gaussian (green).

The bottom panel in the figure highlights a central paradox of quantum mechanics: in the ground state, the particle is most likely to be found at the equilibrium position where the classical oscillator spends the least time. The classical distribution (dashed red) is U-shaped because the mass decelerates near the turning points and lingers there. As v → ∞, the locally averaged quantum density converges to this classical curve—an elegant quantitative demonstration of the correspondence principle first articulated by Niels Bohr.

Quantum Number v	Number of Nodes	Parity	Dominant Feature of \|ψ\|²
0	0	Even	Single Gaussian peak at x = 0
1	1	Odd	Node at x = 0; two symmetric lobes
2	2	Even	Three lobes; central lobe slightly smaller
3	3	Odd	Four lobes; amplitude shifts toward turning points
v (large)	v	(−1)ᵛ	Envelope approaches classical U-shaped distribution

Properties of the first few harmonic oscillator eigenstates

SECTION 6

Worked Example — IR Absorption of HCl

The harmonic oscillator model finds its most direct experimental application in infrared absorption spectroscopy. Diatomic molecules such as HCl absorb infrared photons whose energy matches the gap ℏω between adjacent vibrational levels. The following example demonstrates how to predict the fundamental absorption frequency from tabulated force-constant data.

Predicting the IR Fundamental Band of ¹H³⁵Cl

1
Step 1 — Identify Given Values

The force constant of the H–Cl bond is k = 516 N·m⁻¹. The atomic masses are m_H = 1.008 u and m_Cl = 34.97 u. We must use the reduced mass μ = m_Hm_Cl/(m_H + m_Cl).
2
Step 2 — Compute the Reduced Mass

μ = (1.008 × 34.97)/(1.008 + 34.97) u = 35.25/35.98 u = 0.9797 u. Converting to SI: μ = 0.9797 × 1.6605 × 10⁻²⁷ kg = 1.627 × 10⁻²⁷ kg.

μ = 1.627 × 10⁻²⁷ kg
3
Step 3 — Calculate the Angular Frequency ω

ω = √(k/μ) = √(516 / 1.627 × 10⁻²⁷) = √(3.172 × 10²⁹) = 5.632 × 10¹⁴ rad·s⁻¹.

ω = 5.632 × 10¹⁴ rad·s⁻¹
4
Step 4 — Convert to Wavenumber

The spectroscopic wavenumber is ν̃ = ω/(2πc), where c = 2.998 × 10¹⁰ cm·s⁻¹. ν̃ = 5.632 × 10¹⁴ / (2π × 2.998 × 10¹⁰) = 5.632 × 10¹⁴ / 1.884 × 10¹¹ = 2990 cm⁻¹.

ν̃ ≈ 2990 cm⁻¹
5
Step 5 — Compare with Experiment

The experimentally observed fundamental band of HCl is at approximately 2886 cm⁻¹. Our harmonic-oscillator prediction of ≈ 2990 cm⁻¹ overestimates the frequency by about 3.6%, which reflects the neglect of anharmonicity. The Morse oscillator model corrects for this by including higher-order terms in the potential, which slightly lower the transition frequency and also permit overtone transitions (Δv > 1).

Harmonic model: 2990 cm⁻¹ vs. Experiment: 2886 cm⁻¹ (≈ 3.6% error)

SECTION 7

Strengths & Limitations of the Harmonic Approximation

Like every model in physical chemistry, the harmonic oscillator is a simplification of reality. Its power lies in its analytic solvability and broad applicability, but its limitations become apparent when high vibrational excitation or bond dissociation are involved. The table below provides a systematic comparison.

Feature	Strength	Limitation
Energy Spacing	Equally spaced levels provide a clean, analytic formula for transition energies.	Real molecules show convergent spacing as v increases, deviating from equal gaps.
Selection Rules	Predicts the fundamental Δv = ±1 transition observed in IR spectra.	Cannot account for overtones (Δv = ±2, ±3, …) or hot bands without anharmonic corrections.
Dissociation	N/A — model does not include dissociation.	The parabolic potential rises to infinity; real bonds break at finite energy.
Zero-Point Energy	Correctly predicts non-zero ground-state energy, verified by isotopic substitution experiments.	The precise value is approximate; anharmonicity shifts the true ZPE.
Mathematical Tractability	Exact, closed-form solutions; ladder-operator algebra generalizes to many-body theory.	Perturbation or variational methods needed once anharmonic, rotation–vibration coupling, or multi-mode terms are added.

Strengths and limitations of the quantum harmonic oscillator model

✦

KEY TAKEAWAY

The harmonic oscillator is the hydrogen atom of vibrational physics: it is the exactly solvable benchmark against which all more realistic models are judged. Just as you perturb the hydrogen atom to study multi-electron atoms, you perturb the harmonic oscillator (adding cubic and quartic terms) to study anharmonic molecular vibrations. Its algebraic structure—ladder operators, number states, and equally spaced eigenvalues—recurs throughout quantum optics, condensed-matter physics, and quantum field theory, making it arguably the single most transferable model in all of quantum mechanics.

SECTION 8

Connection to Anharmonic & Advanced Models

The harmonic oscillator is rarely the final word in a physical chemistry analysis, but it is almost always the first. To model real diatomic vibrations more accurately, the Morse potential replaces the parabola with V(x) = D_e(1 − e^−βx)², where D_e is the well depth and β controls the width. This potential asymptotically approaches a finite dissociation limit and produces energy levels that converge rather than remaining equally spaced.

Property	Harmonic Oscillator	Morse Oscillator
Potential	V = ½kx² (parabolic)	V = Dₑ(1 − e⁻ᵝˣ)² (asymmetric well)
Energy Levels	Eᵥ = (v + ½)ℏω — equally spaced	Eᵥ = (v + ½)ℏω − (v + ½)²ℏωxₑ — spacing decreases
Dissociation	None; V → ∞ as x → ±∞	Finite number of bound states; V → Dₑ as x → ∞
Selection Rules	Δv = ±1 only	Δv = ±1, ±2, ±3, … (overtones allowed)
Analytic Solution	Exact (Hermite polynomials)	Exact (associated Laguerre functions)

Harmonic vs. Morse oscillator comparison

Beyond the Morse model, polyatomic molecules require the treatment of 3N − 6 (or 3N − 5 for linear molecules) normal modes, each of which is approximated as an independent harmonic oscillator in the zeroth-order picture. Coupling between these modes—Fermi resonance, Coriolis interaction, and Darling–Dennison resonance—introduces further anharmonic corrections. In solid-state physics, the lattice vibrations of a crystal are decomposed into phonons, each described by a harmonic oscillator Hamiltonian. The Debye and Einstein models of heat capacity are direct applications of the quantized oscillator energy levels to macroscopic thermodynamics. Finally, in quantum electrodynamics each mode of the electromagnetic field is a harmonic oscillator whose quanta are photons—making the harmonic oscillator the conceptual backbone of modern quantum field theory.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain, using the Heisenberg uncertainty principle, why the ground-state energy of the quantum harmonic oscillator cannot be zero. What would a zero ground-state energy imply about the position and momentum of the particle simultaneously?

Show answer

If E₀ were zero, the particle would need to be at rest (zero kinetic energy) at the equilibrium position (zero potential energy) simultaneously. This would require Δx = 0 and Δp = 0, violating the Heisenberg uncertainty principle ΔxΔp ≥ ℏ/2. The compromise between localizing the particle (lowering potential energy) and allowing it to spread (lowering kinetic energy) produces a minimum energy E₀ = ½ℏω, which is the zero-point energy. This residual energy is a purely quantum effect with no classical analogue and has measurable consequences, such as the Lamb shift and the Casimir effect.

Question 2: BASIC CALCULATION

A diatomic molecule modeled as a harmonic oscillator has a force constant k = 1860 N·m⁻¹ and a reduced mass μ = 1.139 × 10⁻²⁶ kg. Calculate the fundamental vibrational frequency ν̃ in cm⁻¹ and the zero-point energy in joules and electron-volts.

Show answer

First, ω = √(k/μ) = √(1860/1.139 × 10⁻²⁶) = √(1.633 × 10²⁹) = 4.041 × 10¹⁴ rad·s⁻¹. The wavenumber is ν̃ = ω/(2πc) = 4.041 × 10¹⁴/(2π × 2.998 × 10¹⁰) = 2143 cm⁻¹. This value is characteristic of carbon monoxide (CO). The zero-point energy is E₀ = ½ℏω = ½ × (1.055 × 10⁻³⁴)(4.041 × 10¹⁴) = 2.132 × 10⁻²⁰ J = 0.133 eV. The factor-of-two relationship—ν̃ gives the spacing, and ½ν̃ gives the zero-point contribution—is a hallmark of the equally spaced harmonic spectrum.

Question 3: INTERMEDIATE

Using the ladder-operator formalism, show that ⟨v|x̂²|v⟩ = (ℏ/2mω)(2v + 1). Then use this result to compute the mean square displacement of a CO molecule (ω = 4.041 × 10¹⁴ rad·s⁻¹, μ = 1.139 × 10⁻²⁶ kg) in its ground state and compare it with a typical bond length of 1.128 Å.

Show answer

Express x̂ = √(ℏ/2mω)(â + â†). Then x̂² = (ℏ/2mω)(â + â†)² = (ℏ/2mω)(â² + ââ† + â†â + â†²). Taking the diagonal matrix element, â²|v⟩ and â†²|v⟩ change v by ±2, so their expectation values vanish. The surviving terms give ⟨v|ââ†|v⟩ = v+1 and ⟨v|â†â|v⟩ = v, so ⟨x²⟩ = (ℏ/2mω)(2v+1). For v = 0: ⟨x²⟩ = ℏ/(2mω) = (1.055 × 10⁻³⁴)/(2 × 1.139 × 10⁻²⁶ × 4.041 × 10¹⁴) = 1.145 × 10⁻²³ m². The root-mean-square displacement is √(⟨x²⟩) = 3.39 × 10⁻¹² m = 0.034 Å, which is about 3% of the bond length. This small fraction confirms that the harmonic approximation (small-displacement expansion) is well justified for the ground vibrational state.

Question 4: APPLIED

Deuterium chloride (DCl) has the same force constant as HCl (k = 516 N·m⁻¹) but a different reduced mass because m_D = 2.014 u. Predict the ratio ν̃(HCl)/ν̃(DCl) and use it to estimate the fundamental absorption wavenumber of DCl. Discuss how this isotope shift is used experimentally to assign vibrational bands.

Show answer

Since ν̃ ∝ √(1/μ) and k is unchanged upon isotopic substitution, the ratio is ν̃(HCl)/ν̃(DCl) = √(μ_DCl/μ_HCl). μ_HCl = (1.008 × 34.97)/35.98 = 0.9797 u. μ_DCl = (2.014 × 34.97)/36.98 = 1.9044 u. The ratio = √(1.9044/0.9797) = √(1.944) = 1.394. Therefore ν̃(DCl) = 2990/1.394 ≈ 2145 cm⁻¹ (harmonic prediction). Experimentally, ν̃(DCl) ≈ 2091 cm⁻¹. The isotope shift is a powerful diagnostic: if a band shifts by a factor close to √2 upon H→D substitution, it is confirmed as an X–H stretching mode. This technique is widely used in assigning complex IR spectra of organic and inorganic molecules.

Question 5: CRITICAL THINKING

The harmonic oscillator selection rule Δv = ±1 arises from the transition dipole moment integral ⟨v'|μ̂|v⟩, where μ̂ is expanded to first order in the displacement x. Derive this selection rule using the ladder-operator expression for x̂ and discuss under what physical circumstances the rule breaks down, giving rise to overtone absorptions.

Show answer

The dipole moment is expanded as μ(x) ≈ μ₀ + (dμ/dx)₀ x + ½(d²μ/dx²)₀ x² + …. The zeroth-order term μ₀ gives no transitions (orthogonality). The first-order term is proportional to ⟨v'|x̂|v⟩. Since x̂ = √(ℏ/2mω)(â + â†), we have ⟨v'|x̂|v⟩ = √(ℏ/2mω)[√v δ_{v',v−1} + √(v+1) δ_{v',v+1}], which is non-zero only for Δv = ±1. This is the harmonic selection rule, and it also requires (dμ/dx)₀ ≠ 0, which is why homonuclear diatomics are IR-inactive. Overtones arise from two sources: (1) mechanical anharmonicity, where the wave functions are no longer pure Hermite-Gaussian and orthogonality with respect to x is relaxed; and (2) electrical anharmonicity, where the higher-order terms ½(d²μ/dx²)₀ x² etc. contribute matrix elements with Δv = ±2, ±3, …. In practice, overtone intensities decrease rapidly (typically by a factor of 10–20 per Δv increment), consistent with the weakness of these correction terms relative to the linear dipole approximation.

SUMMARY

Summary

The quantum harmonic oscillator models a particle in a parabolic potential V(x) = ½kx² and yields equally spaced energy levels given by E_v = (v + ½)ℏω. The ground state carries a zero-point energy of ½ℏω, mandated by the Heisenberg uncertainty principle. The stationary-state wave functions are products of Hermite polynomials and a Gaussian envelope, exhibiting v nodes and non-negligible probability density beyond the classical turning points.

The ladder-operator algebra (â and â†) provides an elegant alternative to solving the Schrödinger equation directly and generalizes to quantum field theory. In molecular spectroscopy, the harmonic model predicts the fundamental IR absorption frequency and the Δv = ±1 selection rule. Its limitations—inability to describe dissociation and overtones—are addressed by the Morse potential and higher-order anharmonic corrections. The correspondence principle is beautifully illustrated as the quantum probability density at large v converges to the classical distribution, reinforcing the oscillator's role as the foundational exactly solvable problem of quantum mechanics.