Ladder Operators for the Harmonic Oscillator

Vibrational spectroscopy still uses the $\Delta v=\pm1$ harmonic selection rule derived above as the zeroth-order guide for interpreting infrared and Raman spectra of every diatomic and polyatomic molecule, with anharmonic corrections computed as perturbative mixing between number states via cubic and quartic terms in $\hat x$. Molecular dynamics force fields and normal-mode analyses represent a polyatomic molecule's small-amplitude vibrations as a set of independent harmonic oscillators (normal modes), each quantized with its own ladder operators, and thermodynamic partition functions for vibrations are built directly from the evenly spaced $E_n=(n+\tfrac12)\hbar\omega$ ladder. In condensed-matter physics, lattice vibrations are quantized as phonons using precisely this creation/annihilation algebra applied mode-by-mode to the crystal's normal coordinates; phonon population statistics, thermal conductivity, and Raman/infrared selection rules in solids all trace back to $[\hat a_k,\hat a_k^\dagger]=1$ for each phonon mode $k$. Quantum optics and cavity quantum electrodynamics use the identical algebra for photon number states $|n\rangle$ of a cavity mode, with coherent states $|\alpha\rangle$ (eigenstates of $\hat a$) describing laser light and squeezed states built by non-trivial combinations of $\hat a,\hat a^\dagger$ used in gravitational-wave detectors such as LIGO to beat the standard quantum limit. In trapped-ion quantum computing, a single ion's motional state in a harmonic trap is manipulated and cooled to its ground state ($n=0$) using exactly the raising/lowering language developed here, and sideband cooling exploits the $\Delta n=\pm1$ coupling to laser fields to reach the vibrational ground state needed for high-fidelity quantum logic gates.

:::quiz[What is the value of $[\hat a,\hat a^\dagger]$ for the harmonic oscillator ladder operators as defined in this article?]

• $i\hbar$
• $\hbar\omega$
• $1$ ✓
• $0$

Ask the tutor: How do the ladder operators change when treating a two-dimensional isotropic harmonic oscillator, and how does that foreshadow angular momentum quantization?

Think of the harmonic oscillator Hamiltonian as an unfactored quadratic form, $\hat H \propto \hat x^2+\hat p^2$ in suitable units. If $\hat x$ and $\hat p$ were ordinary numbers this would factor as a difference of squares, $(\hat x+i\hat p)(\hat x-i\hat p)$, using $i^2=-1$. Because $\hat x$ and $\hat p$ do not commute, the factorization is not exact — but it is exact up to an extra additive constant that comes entirely from the commutator $[\hat x,\hat p]=i\hbar$. This leftover constant is the origin of the famous zero-point energy $\tfrac12\hbar\omega$. Physically, $\hat a^\dagger$ adds one quantum of vibrational energy $\hbar\omega$ to the system, and $\hat a$ removes one. Repeated application of $\hat a$ on any state must eventually terminate — you cannot remove infinitely many quanta from a system whose energy is bounded below — and this termination condition is what pins the ground state energy to exactly $\tfrac12\hbar\omega$ rather than zero. Everything above the ground state is built by climbing the ladder with $\hat a^\dagger$, picking up a $\sqrt{n+1}$ normalization factor at each rung, exactly mirroring how angular momentum raising and lowering operators build the $|j,m\rangle$ ladder used later for the Particle on a Sphere and Spherical Harmonics problem. The position-space wavefunctions $\psi_n(x)$ built from Hermite polynomials are not abandoned by this approach; they are recovered as the unique solution of $\hat a\,\psi_0(x)=0$, a first-order differential equation, followed by repeated differentiation-and-multiplication via $\hat a^\dagger$. The ladder operators are thus a bridge: they give the spectrum algebraically and, when needed, regenerate the explicit real-space functions with far less labor than solving the second-order Hermite equation directly.

The starting Hamiltonian for a particle of mass $m$ on a spring of force constant $k=m\omega^2$ is

where: $\hat x,\hat p$ are the position and momentum operators satisfying $[\hat x,\hat p]=i\hbar$, and $\omega=\sqrt{k/m}$ is the classical angular frequency. Reading: this is exactly the same Hamiltonian solved by Hermite functions; nothing physical has changed, only the algebra used to diagonalize it. Define the dimensionless-in-spirit combinations

where: $\hat a$ is the annihilation (lowering) operator, $\hat a^\dagger$ its Hermitian-conjugate creation (raising) operator, and the prefactor $\sqrt{m\omega/2\hbar}$ is chosen precisely so the resulting commutator is unity. Reading: $\hat a$ is a complex combination of position and momentum, not itself Hermitian, so it is not an observable — but $\hat a^\dagger\hat a$ will be. The fundamental commutation relation, derived below, is

where: the right side is the identity operator (times $\hbar^0$, i.e. a pure number in the chosen units). Reading: this single algebraic fact, inherited entirely from $[\hat x,\hat p]=i\hbar$, is the engine that produces the entire energy ladder. The Hamiltonian rewritten in terms of the number operator $\hat N \equiv \hat a^\dagger\hat a$ is

where: $\hat N$ is Hermitian ($\hat N^\dagger = (\hat a^\dagger\hat a)^\dagger = \hat a^\dagger\hat a = \hat N$) so its eigenvalues are real, and in fact turn out to be the non-negative integers $n=0,1,2,\dots$. Reading: once you know $\hat N|n\rangle = n|n\rangle$, the energies fall out immediately as $E_n=(n+\tfrac12)\hbar\omega$. The ladder action on normalized eigenstates $|n\rangle$ is

where: the square-root prefactors are exactly what keep $|n\rangle$ normalized to unity at every step; they are not free choices but forced by the algebra, as shown in the derivation below. Reading: this is the operational heart of the method — every matrix element, every selection rule, and every explicit wavefunction follows from these two lines. Position and momentum recovered from the ladder operators are

where: these are just the inverted definitions of $\hat a,\hat a^\dagger$. Reading: because $\hat a,\hat a^\dagger$ only connect states differing by one quantum, $\hat x$ and $\hat p$ can only have nonzero matrix elements between $|n\rangle$ and $|n\pm1\rangle$ — the algebraic origin of the $\Delta v=\pm1$ vibrational selection rule (in the harmonic approximation).

:::derivation[The commutator $[\hat a,\hat a^\dagger]=1$ from $[\hat x,\hat p]=i\hbar$] Goal. Show that the specific prefactor $\sqrt{m\omega/2\hbar}$ chosen in the definitions of $\hat a,\hat a^\dagger$ makes their commutator exactly the number $1$, using only $[\hat x,\hat p]=i\hbar$.

Step 1 — Write the product $\hat a\hat a^\dagger$. Substitute the definitions:

Expand the product term by term, being careful not to commute $\hat x$ and $\hat p$ past each other:

Step 2 — Group the cross terms. The two middle terms combine as

Insert $[\hat x,\hat p]=i\hbar$:

So

Step 3 — Compute $\hat a^\dagger\hat a$ the same way. Swapping the order of the factors flips the sign of the cross term:

Notice the quadratic terms in $\hat x^2,\hat p^2$ are identical in both expressions — only the sign of the constant flips.

Step 4 — Subtract. All the $\hat x^2$ and $\hat p^2$ pieces cancel exactly:

This confirms $[\hat a,\hat a^\dagger]=1$ as a pure number, independent of $m,\omega,\hbar$ — those constants were absorbed entirely into the choice of prefactor in Step 1. Equivalently, adding the two boxed expressions from Steps 2–3 gives $\hat a\hat a^\dagger+\hat a^\dagger\hat a = \frac{m\omega}{\hbar}\hat x^2+\frac{1}{m\hbar\omega}\hat p^2$, and comparing to $\hat H=\frac{\hat p^2}{2m}+\frac12 m\omega^2\hat x^2$ shows

exactly the factorization promised.

Example 1 — Foundational: building $|2\rangle$ from the vacuum and checking its energy. Starting from $|0\rangle$, defined by $\hat a|0\rangle=0$ and $\langle 0|0\rangle=1$, construct $|2\rangle$ and verify $\hat H|2\rangle = \tfrac52\hbar\omega|2\rangle$. Apply $\hat a^\dagger$ once:

Apply $\hat a^\dagger$ a second time, using the ladder rule with $n=1$:

So $(\hat a^\dagger)^2|0\rangle = \sqrt{2}\,|2\rangle$, i.e.

Now check the energy directly using $\hat H=\hbar\omega(\hat a^\dagger\hat a+\tfrac12)$. First lower: $\hat a|2\rangle=\sqrt{2}\,|1\rangle$. Then raise again: $\hat a^\dagger(\sqrt2|1\rangle) = \sqrt2\cdot\sqrt{2}\,|2\rangle = 2|2\rangle$. So $\hat a^\dagger\hat a|2\rangle = \hat N|2\rangle = 2|2\rangle$, confirming $n=2$ directly by explicit operator action rather than by assumption. Therefore

matching $E_n=(n+\tfrac12)\hbar\omega$ at $n=2$ exactly, with no reference to Hermite polynomials.

Example 2 — Applied: vibrational excitation energy of gas-phase CO and the algebraic origin of its IR selection rule. Carbon monoxide has a harmonic vibrational wavenumber $\tilde\nu = 2170\ \text{cm}^{-1}$ (a genuine, oft-tabulated value; note it is not a round number, by design). Using $E=hc\tilde\nu$, the fundamental quantum is

Compute the product stepwise:

The energy of the $v=0\to v=1$ transition is exactly this quantum, since

regardless of which rung you start from — the ladder spacing is uniform, a direct consequence of $\hat N$ having integer eigenvalues one apart. Why does only $\Delta v=\pm1$ appear (in the harmonic/double-harmonic approximation)? The transition dipole matrix element for absorption involves $\langle v'|\hat x|v\rangle$ (since the dipole moment is expanded linearly in the bond-length displacement for the harmonic approximation). Using $\hat x=\sqrt{\hbar/2m\omega}\,(\hat a+\hat a^\dagger)$,

Every term on the right is nonzero only when $v'=v\pm1$; there is no way to get $v'=v\pm2$ or higher from a single power of $\hat a+\hat a^\dagger$. This is the operator-algebra proof of the harmonic selection rule, sharper and shorter than doing the corresponding Hermite-polynomial integral by hand.

Example 3 — Competition-level: matrix elements of $\hat x^2$ and the anharmonic correction structure. Compute $\langle 3|\hat x^2|1\rangle$ exactly (a nonzero, non-diagonal, non-adjacent-looking matrix element that tests whether the ladder algebra is applied correctly), for a particle of mass $m$ and frequency $\omega$. Write $\hat x^2$ in terms of ladder operators:

Act term by term on $|1\rangle$. Term $\hat a^2|1\rangle$: first $\hat a|1\rangle=\sqrt1|0\rangle=|0\rangle$, then $\hat a|0\rangle=0$. So $\hat a^2|1\rangle=0$, contributing nothing to $\langle3|\hat a^2|1\rangle$. Term $\hat a\hat a^\dagger|1\rangle$: first $\hat a^\dagger|1\rangle=\sqrt2|2\rangle$, then $\hat a(\sqrt2|2\rangle)=\sqrt2\cdot\sqrt2|1\rangle=2|1\rangle$. This produces $|1\rangle$, orthogonal to $|3\rangle$, so it contributes $0$ to $\langle3|\cdots|1\rangle$. Term $\hat a^\dagger\hat a|1\rangle$: $\hat a|1\rangle=|0\rangle$, then $\hat a^\dagger|0\rangle=|1\rangle$. Again produces $|1\rangle$, contributing $0$ to the $\langle3|$ overlap. Term $(\hat a^\dagger)^2|1\rangle$: first $\hat a^\dagger|1\rangle=\sqrt2|2\rangle$, then $\hat a^\dagger(\sqrt2|2\rangle)=\sqrt2\cdot\sqrt3|3\rangle=\sqrt6|3\rangle$. This is the only term with nonzero overlap onto $\langle3|$, giving $\langle3|(\hat a^\dagger)^2|1\rangle=\sqrt6$. Collecting only the surviving term,

This nonzero $\Delta v=\pm2$ matrix element of $\hat x^2$ (as opposed to $\hat x$) is exactly why overtone and combination-band-inducing terms in real anharmonic potentials — which contain $\hat x^3,\hat x^4,\dots$ — connect states two, three, or four rungs apart, and why Raman scattering (which couples to $\hat x^2$-type polarizability terms) shows different selection rules from IR absorption.

Paul Adrien Maurice Dirac introduced the operator method for the harmonic oscillator in the earliest years of matrix and operator quantum mechanics, treating $\hat x$ and $\hat p$ as abstract non-commuting objects rather than differential operators acting on wavefunctions. His approach — factor the Hamiltonian, extract the algebra, read off the spectrum from positivity of the norm — became the template not only for the oscillator but for angular momentum, spin, and eventually the creation/annihilation operator formalism of quantum electrodynamics, where photons themselves are literally treated as quanta of a harmonic-oscillator-like field mode.