The Physics Travel Guide

Harmonic Oscillator

Prerequisites: Newton's Second Law, Differential Equations. Full path (3 pages)

L=m2(dxdt)2−k2x2L= \frac{ \color{blue}{m}}{2} \left( \color{olive}{\frac{dx}{dt}}\right)^2 - \frac{ \color{red}{k}}{2} \color{magenta}{x}^2

Intuitive

A harmonic oscillator simply consists of an object attached to a spring and is one of the simplest physical systems we can study.

To describe the system we usually track the position\color{magenta}{\text{position}} of the object at the end of the spring. The continuous up and down movement is called an oscillation, which is where the name "oscillator" comes from.

The movement depends crucially on the properties of the spring\color{red}{\text{properties of the spring}}. When the spring is harder, the object will oscillate faster. In addition, the motion of the object depends on its mass\color{blue}{\text{mass}}. A heavier object will oscillate more slowly.

Concrete

Classical oscillator. From the Lagrangian L=m2x˙2−k2x2L = \frac{m}{2}\dot x^2 - \frac{k}{2}x^2 the Euler-Lagrange equation gives mx¨=−kx,orx¨+ω2x=0,ω=k/m.m\ddot x = -kx , \qquad \text{or} \qquad \ddot x + \omega^2 x = 0, \quad \omega = \sqrt{k/m} . The general solution is x(t)=Acos⁡(ωt+φ)x(t) = A\cos(\omega t + \varphi), where the amplitude AA and phase φ\varphi are fixed by the initial position and velocity. The period T=2πm/kT = 2\pi\sqrt{m/k} does not depend on the amplitude: a stiffer spring makes the oscillation faster, a heavier mass makes it slower. The energy E=m2x˙2+k2x2=12kA2E = \frac{m}{2}\dot x^2 + \frac{k}{2}x^2 = \frac{1}{2}kA^2 is conserved and moves back and forth between kinetic and potential form.

Why it is everywhere. Any system near a stable equilibrium x0x_0 of a smooth potential V(x)V(x) is approximately a harmonic oscillator. Expanding, V(x)=V(x0)+12V′′(x0)(x−x0)2+O((x−x0)3),V(x) = V(x_0) + \frac{1}{2}V''(x_0)(x - x_0)^2 + \mathcal O\big((x-x_0)^3\big) , since V′(x0)=0V'(x_0) = 0. For small displacements the effective spring constant is k=V′′(x0)>0k = V''(x_0) > 0. The same argument applied to a system with many degrees of freedom gives a set of decoupled oscillators, the normal modes, after diagonalizing the matrix of second derivatives. Molecular vibrations, phonons in crystals and the modes of a free field are all treated this way.

Quantum oscillator. With k=mω2k = m\omega^2 the Hamiltonian is H=p22m+12mω2x2,[x,p]=iℏ.H = \frac{p^2}{2m} + \frac{1}{2}m\omega^2 x^2, \qquad [x, p] = i\hbar . Define the ladder operators a=mω2ℏ(x+ipmω),a†=mω2ℏ(x−ipmω).a = \sqrt{\frac{m\omega}{2\hbar}}\left(x + \frac{i p}{m\omega}\right), \qquad a^\dagger = \sqrt{\frac{m\omega}{2\hbar}}\left(x - \frac{i p}{m\omega}\right) . They satisfy [a,a†]=1[a, a^\dagger] = 1, and the Hamiltonian becomes H=ℏω(a†a+12)H = \hbar\omega\left(a^\dagger a + \frac12\right). From [H,a]=−ℏω a[H, a] = -\hbar\omega\, a and [H,a†]=ℏω a†[H, a^\dagger] = \hbar\omega\, a^\dagger, the operator a†a^\dagger raises the energy by ℏω\hbar\omega and aa lowers it. Since HH is bounded below, there is a state with a∣0⟩=0a|0\rangle = 0, and the spectrum is En=ℏω(n+12),∣n⟩=(a†)nn!∣0⟩,n=0,1,2,…E_n = \hbar\omega\left(n + \frac12\right), \qquad |n\rangle = \frac{(a^\dagger)^n}{\sqrt{n!}}|0\rangle, \qquad n = 0, 1, 2, \dots In position space a∣0⟩=0a|0\rangle = 0 is a first-order differential equation with solution ψ0(x)=(mωπℏ)1/4e−mωx2/2ℏ\psi_0(x) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4} e^{-m\omega x^2/2\hbar}. The ground state energy 12ℏω\frac12\hbar\omega is not zero: the particle cannot sit at rest at the bottom of the well, as the uncertainty principle requires. The excited states are ψn∝Hn ⁣(mω/ℏ x)ψ0\psi_n \propto H_n\!\left(\sqrt{m\omega/\hbar}\,x\right)\psi_0, with HnH_n the Hermite polynomials.

Abstract

Classical. The phase space is R2\mathbb R^2 with coordinates (q,p)(q, p) and symplectic form Ω=dp∧dq\Omega = dp \wedge dq. In the rescaled coordinates Q=mω qQ = \sqrt{m\omega}\, q, P=p/mωP = p/\sqrt{m\omega} (still canonical), the Hamiltonian is H=ω2(P2+Q2)H = \frac{\omega}{2}(P^2 + Q^2) and its Hamiltonian flow is rigid rotation with angular velocity ω\omega. Thus H/ωH/\omega is the moment map of a U(1)U(1) action on phase space. Action-angle variables are I=H/ωI = H/\omega and the rotation angle θ\theta, with Ω=dI∧dθ\Omega = dI \wedge d\theta; the orbits are circles enclosing area ∮p dq=2πI\oint p\, dq = 2\pi I. For nn degrees of freedom with equal frequencies, H=ω2∣Z∣2H = \frac{\omega}{2}|Z|^2 with Z=Q+iP∈CnZ = Q + iP \in \mathbb C^n, and the group of linear symplectic maps preserving HH is U(n)U(n). Quadratic Hamiltonians on R2n\mathbb R^{2n} form the Lie algebra sp(2n,R)\mathfrak{sp}(2n, \mathbb R) under the Poisson bracket, and the positive-definite isotropic oscillator generates the center of its maximal compact subalgebra u(n)\mathfrak u(n).

Quantum. Quantization represents the Heisenberg algebra [a,a†]=1[a, a^\dagger] = 1. By the Stone-von Neumann theorem, for finitely many degrees of freedom the irreducible representation with aa and a†a^\dagger adjoint to each other (integrated to the Weyl relations) is unique up to unitary equivalence. A convenient model is the Bargmann-Fock space of holomorphic functions f(z)f(z) with ⟨f,g⟩=1π∫Cf(z)‾ g(z) e−∣z∣2 d2z,\langle f, g\rangle = \frac{1}{\pi}\int_{\mathbb C} \overline{f(z)}\, g(z)\, e^{-|z|^2}\, d^2 z , where a†a^\dagger acts as multiplication by zz, aa as ∂z\partial_z, and zn/n!z^n/\sqrt{n!} is the orthonormal eigenbasis ∣n⟩|n\rangle of a†aa^\dagger a.

The quadratic operators K0=14(a†a+aa†),K+=12(a†)2,K−=12a2K_0 = \frac14\left(a^\dagger a + a a^\dagger\right), \qquad K_+ = \frac12 (a^\dagger)^2, \qquad K_- = \frac12 a^2 satisfy [K0,K±]=±K±[K_0, K_\pm] = \pm K_\pm and [K−,K+]=2K0[K_-, K_+] = 2K_0, a realization of sp(2,R)≅sl(2,R)\mathfrak{sp}(2, \mathbb R) \cong \mathfrak{sl}(2, \mathbb R). The Hamiltonian is H=2ℏωK0H = 2\hbar\omega K_0. This is the metaplectic (oscillator) representation: it exponentiates to the double cover Mp(2,R)Mp(2,\mathbb R) of Sp(2,R)Sp(2,\mathbb R), not to Sp(2,R)Sp(2,\mathbb R) itself, and it splits into two irreducible lowest-weight representations on even and odd nn, with lowest K0K_0 weights 14\frac14 and 34\frac34. The zero-point energy 12ℏω\frac12\hbar\omega is tied to this double cover: in geometric quantization it appears through the half-form (metaplectic) correction, and the Bohr-Sommerfeld condition ∮p dq=2πℏ(n+12)\oint p\,dq = 2\pi\hbar\left(n + \frac12\right) including the Maslov correction reproduces the exact spectrum.

A free quantum field is a collection of such oscillators, one per mode, and its Fock space is built from the same algebra. For infinitely many degrees of freedom the Stone-von Neumann uniqueness fails, which is the origin of inequivalent representations in quantum field theory.

The phase space of a harmonic oscillator

Why is it interesting?

“The career of a young theoretical physicist consists of treating the harmonic oscillator in ever-increasing levels of abstraction.” Sidney Coleman

In a sense, the harmonic oscillator is to physics what the set of natural numbers is to mathematics. It is a simple, almost “trivial” system, but one which conceals much subtlety and beauty and from which a great deal of what is of interest in the subject evolves. Gregory L. Naber

The harmonic oscillator is an important model system pervading many areas in classical physics; it is likewise ubiquitous in quantum mechanics. The nonrelativistic Scrodinger equation with a harmonic oscillator potential is readily ¨ solved with standard analytic methods, whether in one or three dimensions. However, we will take a different tack here, and address the one-dimension problem more as an excuse to introduce the notion of “creation” and “annihilation” operators, or “step-up” and “step-down” operators. This is an example of a type of operation which will repeat itself in many contexts, including the theory of angular momentum.

The Simple Harmonic Oscillator by F. Porter

FAQ

Why do we introduce the ladder operators to treat the harmonic oscillator in quantum mechanics?

Path to this page

  1. Calculus
  2. Differential Equations
  3. Newton's Second Law

Leads to: Pendulum.