Quantum Harmonic Oscillator

Why the harmonic oscillator matters

Any smooth potential, expanded around a stable minimum, is locally parabolic: V(x) ≈ V₀ + ½mω²(x − x₀)². The harmonic oscillator is the universal model for any system near equilibrium — the leading non-trivial term in the Taylor series. Solve it once, and you have a template for chemical bonds, crystal vibrations, photon modes, fields in spacetime, and trapped atoms. No other model in physics earns its keep across so many domains.

• Molecular vibrations. Diatomic molecules vibrate at IR frequencies. H–Cl stretches at 2990 cm⁻¹ (ℏω ≈ 371 meV); zero-point energy ≈ 186 meV. The vibrational spectrum (visible in IR spectroscopy as absorption peaks) is the QHO ladder in action.
• Phonons in solids. Lattice vibrations are collective harmonic motions. Each normal mode of a crystal is a QHO; its quanta are phonons. Heat capacity, thermal conductivity, sound velocity, and electron-phonon coupling all flow from this picture.
• Photons and free fields. Quantize the electromagnetic field mode-by-mode and each mode becomes a QHO. Photons are its quanta. Same for gluons, gravitons, and every free quantum field: free QFT is just an infinite collection of harmonic oscillators.
• Trapped ions and cold atoms. Ion traps confine atoms in approximately harmonic potentials. Sideband cooling, motional coherent states, and quantum logic gates on motional qubits all rely on QHO physics.
• Vacuum energy and Casimir effect. Each EM mode contributes ½ℏω to the vacuum. Mode-density differences between gaps and walls produce the measurable Casimir attraction. The unresolved cosmological constant problem is the same zero-point sum, taken seriously across all field modes.

The Hamiltonian and its spectrum

H = p²/(2m) + ½mω²x²

Two ways to solve: directly (substitute, get the Hermite differential equation, find polynomial solutions), or algebraically (introduce ladder operators and bypass differential equations entirely). The algebraic route, due to Dirac, is the cleanest.

Define a = (mωx + ip) / √(2mℏω) and its adjoint a†. Then H = ℏω(a†a + ½), and from [x, p] = iℏ we get [a, a†] = 1. Let |n⟩ be an eigenstate of n̂ = a†a with eigenvalue n. The commutator implies a†|n⟩ ∝ |n+1⟩ and a|n⟩ ∝ |n−1⟩. Normalization gives a†|n⟩ = √(n+1)|n+1⟩, a|n⟩ = √n|n−1⟩. The chain must terminate (energies must be bounded below), forcing a|0⟩ = 0 — the ground state — and n ∈ {0, 1, 2, ...}.

Energy spectrum:

E_n = ℏω(n + ½), n = 0, 1, 2, 3, ...

E₀ = ½ℏω    (zero-point energy)
E₁ = (3/2)ℏω
E₂ = (5/2)ℏω
...
ΔE = E_{n+1} − E_n = ℏω    (constant spacing)

Wave functions

Solve a|0⟩ = 0 in position space — it's a first-order ODE giving a Gaussian. Apply (a†)^n to climb the ladder. Result: ψ_n(x) = (mω/πℏ)^{1/4} · (1/√(2^n n!)) · H_n(ξ) · exp(−ξ²/2), where ξ = x/x₀ and x₀ = √(ℏ/mω) is the natural length scale; H_n are Hermite polynomials.

n	E_n	ψ_n(x) (unnormalized form)	Nodes
0	½ℏω	exp(−ξ²/2)	0
1	(3/2)ℏω	2ξ · exp(−ξ²/2)	1
2	(5/2)ℏω	(4ξ² − 2) · exp(−ξ²/2)	2
3	(7/2)ℏω	(8ξ³ − 12ξ) · exp(−ξ²/2)	3
n	(n+½)ℏω	H_n(ξ) · exp(−ξ²/2)	n

Each ψ_n has exactly n nodes (zeros), and the wave function leaks into the classically forbidden region |x| > x_classical. Probability oscillates inside the well, then decays Gaussianly outside.

Zero-point energy — atom-vibration estimate

For a diatomic molecule with ℏω ≈ 200 meV (typical mid-IR vibration), the zero-point energy is E₀ = ½ × 200 = 100 meV per mode. Three vibrational modes in a non-linear triatomic give ~0.3 eV of irreducible vibrational energy at zero kelvin. Thermal energy at 300 K (kT ≈ 25 meV) is much smaller, so most molecules sit in the vibrational ground state at room temperature — but with the n=0 zero-point energy locked in, not in the classical "frozen at x=0" sense.

This is measurable. Isotope substitution (D for H) reduces ω by 1/√2 because reduced mass doubles. Zero-point energy drops by the same factor. Reaction rates, bond lengths, and equilibrium constants shift accordingly — measurable kinetic isotope effects.

Common mistakes

• Saying the ground state has zero energy. E₀ = ½ℏω, not zero. The lowest classical energy is at x=0, p=0 — quantum mechanics forbids this exact configuration.
• Confusing the QHO with a real spring at any amplitude. Real springs and bonds are only parabolic near equilibrium. Beyond a few x₀, anharmonicity kicks in (overtone shifts, dissociation, level crowding). QHO is the leading approximation, not the full story.
• Forgetting that "ω" is angular frequency. The level spacing is ℏω = hν with ν = ω/(2π). Mid-IR ν ~ 3×10¹³ Hz → ω ~ 2×10¹⁴ rad/s → ℏω ~ 0.13 eV. Get the factor of 2π right.
• Treating a and a† as Hermitian. They are not — they are conjugates of each other. Only a†a (the number operator) and a + a† ∝ x or i(a − a†) ∝ p are Hermitian.
• Picturing only Hermite functions. Eigenstates are stationary — they don't slosh in the well. The classical-like sloshing motion comes from coherent superpositions of eigenstates (coherent states), which oscillate with the classical period 2π/ω.
• Ignoring degeneracy in higher dimensions. A 2D isotropic QHO has E = ℏω(n_x + n_y + 1) with n_x + n_y = N defining an N+1-fold degenerate level. 3D has ½(N+1)(N+2) degeneracy. Important in shell models of nuclei and quark confinement.