Skip to content

Oscillator as a Universal Local Model

The harmonic oscillator appears everywhere because smooth systems near stable equilibrium are locally quadratic. The exact oscillator is one solvable Hamiltonian, but the oscillator approximation is much broader: it is the first local model of small quantum fluctuations around a stable configuration.

This is the reason oscillator language appears in molecular vibrations, trapped ions, lattice phonons, electromagnetic modes, small field fluctuations, and many normal-mode decompositions.

Let x0x_0 be a stable minimum of a smooth potential. Then

V′(x0)=0,V′′(x0)>0.V'(x_0)=0, \qquad V''(x_0)\gt 0.

For a small displacement

η=x−x0,\eta=x-x_0,

Taylor expansion gives

V(x0+η)=V(x0)+12V′′(x0)η2+16V′′′(x0)η3+⋯ .V(x_0+\eta) = V(x_0) +\frac12V''(x_0)\eta^2 +\frac16V'''(x_0)\eta^3 +\cdots.

For sufficiently small η\eta, the quadratic term is the leading nonconstant term. Define

K=V′′(x0),ω=Km.K=V''(x_0), \qquad \omega=\sqrt{\frac{K}{m}}.

Then the local Hamiltonian is

H^≈p^22m+V(x0)+12mω2(x^−x0)2.\hat H \approx \frac{\hat p^2}{2m} +V(x_0) +\frac12m\omega^2(\hat x-x_0)^2.

Thus the low-lying energies are approximately

En≈V(x0)+ℏω(n+12),n=0,1,2,…E_n \approx V(x_0) +\hbar\omega \left( n+\frac12 \right), \qquad n=0,1,2,\ldots

The constant V(x0)V(x_0) sets the energy zero for that local well. The spacing is controlled by the curvature V′′(x0)V''(x_0) and the mass.

The oscillator approximation comes from two durable facts:

  • smooth functions have Taylor expansions;
  • a stable equilibrium has no linear term and has positive quadratic curvature.

The linear term vanishes because the force is zero at equilibrium. The quadratic term is the first restoring force. Quantum mechanics then turns this local quadratic Hamiltonian into a ladder of states with spacing ℏω\hbar\omega and a finite ground-state width.

This universality is local. It does not claim that the whole potential is parabolic, only that sufficiently small fluctuations around a stable point are governed first by the parabolic part.

The oscillator approximation is strongest for low-lying states whose wavefunctions live mostly near the minimum. The natural ground-state width is

Δx0=ℓ2,ℓ=ℏmω.\Delta x_0=\frac{\ell}{\sqrt2}, \qquad \ell=\sqrt{\frac{\hbar}{m\omega}}.

Anharmonic terms are small when the typical displacement sampled by the state is small compared with the length scale over which the potential curvature changes.

For example, comparing the cubic correction with the quadratic term at displacement η\eta gives the rough dimensionless ratio

ϵ3(η)∼∣V′′′(x0)∣ ∣η∣3V′′(x0).\epsilon_3(\eta) \sim \frac{\lvert V'''(x_0)\rvert\,\lvert\eta\rvert} {3V''(x_0)}.

If ϵ3(ℓ)≪1\epsilon_3(\ell)\ll1, the ground state is well described by the harmonic approximation. Higher excited states sample larger displacements, so the approximation usually worsens as nn increases.

For coordinates collected into a vector q\mathbf q, expand around an equilibrium q0\mathbf q_0:

q=q0+η.\mathbf q=\mathbf q_0+\boldsymbol\eta.

The potential has the local form

V(q0+η)=V(q0)+12ηTKη+⋯ ,V(\mathbf q_0+\boldsymbol\eta) = V(\mathbf q_0) +\frac12\boldsymbol\eta^T K\boldsymbol\eta +\cdots,

where KK is the Hessian matrix,

Kij=∂2V∂qi∂qj∣q=q0.K_{ij} = \left. \frac{\partial^2 V}{\partial q_i\partial q_j} \right\rvert_{\mathbf q=\mathbf q_0}.

With a mass matrix MM, the small-oscillation problem reduces to normal modes. In mass-weighted coordinates, the relevant symmetric matrix is

D=M−1/2KM−1/2.D=M^{-1/2}KM^{-1/2}.

Its eigenvalues are the squared normal-mode frequencies:

Deα=ωα2eα.D\mathbf e_\alpha=\omega_\alpha^2\mathbf e_\alpha.

After diagonalization, the quadratic Hamiltonian becomes a sum of independent oscillators:

H≈V(q0)+∑α(Pα22+12ωα2Qα2).H \approx V(\mathbf q_0) + \sum_\alpha \left( \frac{P_\alpha^2}{2} +\frac12\omega_\alpha^2Q_\alpha^2 \right).

Quantization gives

E{nα}≈V(q0)+∑αℏωα(nα+12).E_{\{n_\alpha\}} \approx V(\mathbf q_0) + \sum_\alpha \hbar\omega_\alpha \left( n_\alpha+\frac12 \right).

This is the general normal-mode pattern. The two-oscillator version is worked out explicitly in Coupled Oscillators: First Encounter.

SystemLocal oscillator coordinateWhat sets the frequencyMain caution
Diatomic moleculebond-length displacementcurvature of the molecular potential at equilibriumanharmonicity matters near dissociation
Polyatomic moleculenormal-mode coordinateHessian of the potential energy surfacerotations, translations, and mode coupling need separation
Crystal latticecollective displacement modeforce constants and mass densityphonons are normal-mode quanta, not atoms
Trapped iondisplacement in the traptrap curvature and ion massmicromotion and coupling to internal states can matter
Electromagnetic modefield-mode amplitudecavity or free-space mode frequencyfield quantization and gauge conventions enter
Small field fluctuationfluctuation around a backgroundquadratic part of the field Hamiltonian or actioninteractions and renormalization go beyond the free modes

These examples are not the same physical system. They share the local quadratic structure.

The oscillator becomes especially powerful when a system has many small-amplitude degrees of freedom. The quadratic approximation is first diagonalized classically into normal modes. Then each normal mode is quantized as an oscillator.

For a finite collection of modes this gives a finite sum. For a field in a box, the mode labels become wavevectors. In the continuum limit the sum becomes an integral, and the oscillator algebra becomes the starting point for field quantization.

The bridge is developed in Harmonic Oscillator to Fields. The important restriction is that free fields are quadratic. Interactions are the higher-order terms, just as anharmonicity is the higher-order correction to a molecular vibration.

The oscillator approximation is a local statement, not a global solution.

It does not determine:

  • tunneling between separated minima;
  • dissociation of a molecule at large bond length;
  • hard-wall boundary conditions far from the minimum;
  • highly excited states that explore anharmonic regions;
  • spectra near unstable equilibria;
  • zero modes associated with translations, rotations, or continuous symmetries;
  • interactions between normal-mode quanta.

If the quadratic term is absent or has zero curvature, the leading local model changes. If the curvature is negative, the point is unstable rather than an oscillator minimum. If the potential is nonsmooth, Taylor expansion may not be the right tool.

  • Saying every bound potential is globally a harmonic oscillator.
  • Using the harmonic approximation for states that sample far outside the quadratic region.
  • Ignoring the constant offset V(x0)V(x_0) when comparing different minima.
  • Treating normal-mode quanta as particles attached to original coordinates.
  • Forgetting zero modes when a symmetry leaves the energy unchanged.
  • Calling an unstable maximum an oscillator just because the potential has a quadratic term.
  • Presenting field modes as independent oscillators without stating that this is the free, quadratic approximation.
  • L. D. Landau and E. M. Lifshitz, Mechanics, 3rd ed., Butterworth-Heinemann, 1976.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018.
  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
  • M. Born and K. Huang, Dynamical Theory of Crystal Lattices, Oxford University Press, 1954.
  • N. W. Ashcroft and N. D. Mermin, Solid State Physics, Saunders, 1976.
  1. Let
V(x)=a(x2−b2)2,a>0,b>0.V(x)=a(x^2-b^2)^2, \qquad a\gt 0, \qquad b\gt 0.

Find the harmonic frequency near the minimum at x=bx=b.

Solution

Compute derivatives:

V′(x)=4ax(x2−b2),V'(x)=4ax(x^2-b^2),

and

V′′(x)=4a(3x2−b2).V''(x)=4a(3x^2-b^2).

At x=bx=b,

V′′(b)=8ab2.V''(b)=8ab^2.

Thus the local oscillator frequency is

ω=8ab2m.\omega=\sqrt{\frac{8ab^2}{m}}.
  1. A potential near x0x_0 is
V(x0+η)=V(x0)+12Kη2+⋯ .V(x_0+\eta) = V(x_0)+\frac12K\eta^2+\cdots.

Classify the local behavior for K>0K\gt 0, K=0K=0, and K<0K\lt 0.

Solution

If K>0K\gt 0, the point is a stable quadratic minimum and the local model is a harmonic oscillator. If K=0K=0, the quadratic term does not supply a restoring force; higher terms or symmetry directions must be examined. If K<0K\lt 0, the point is locally unstable, like an inverted oscillator rather than a bound harmonic oscillator.

  1. A two-dimensional system has mass mm in both coordinates and local potential
V=V0+12K1η12+12K2η22.V = V_0 +\frac12K_1\eta_1^2 +\frac12K_2\eta_2^2.

What are the approximate quantum energies?

Solution

The two coordinates are already normal modes. Their frequencies are

ω1=K1m,ω2=K2m.\omega_1=\sqrt{\frac{K_1}{m}}, \qquad \omega_2=\sqrt{\frac{K_2}{m}}.

The energies are

En1,n2≈V0+ℏω1(n1+12)+ℏω2(n2+12).E_{n_1,n_2} \approx V_0 +\hbar\omega_1 \left( n_1+\frac12 \right) +\hbar\omega_2 \left( n_2+\frac12 \right).
  1. Why does an exact double-well potential require more than one local oscillator approximation?
Solution

Each minimum has its own local quadratic approximation, which describes low-lying states localized near that minimum. The full double-well problem also includes tunneling through the barrier and global parity combinations of localized states. Those effects depend on the barrier between minima and are not determined by either local oscillator alone.