|q⟩ Bad Qubits

intermediate · Physics · The Harmonic Oscillator: Ladder Operators

Creation and Annihilation Operators

The breakthrough of the operator method is to factor the Hamiltonian. A sum of two squares like p2+(mωx)2p^2 + (m\omega x)^2 does not factor over the reals, but it nearly does over the complex numbers — and the small correction it leaves behind is exactly the zero-point energy.

Natural units

To keep the algebra clean we set =m=ω=1\hbar = m = \omega = 1 for the rest of the module; physical factors can be restored at the end by dimensional analysis. In these units the Hamiltonian is

H^=12(p^2+x^2),[x^,p^]=i.\hat{H} = \frac{1}{2}\left(\hat{p}^2 + \hat{x}^2\right), \qquad [\hat{x}, \hat{p}] = i.

Defining the ladder operators

We introduce two operators built from x^\hat{x} and p^\hat{p}:

a=x^+ip^2,a=x^ip^2.a = \frac{\hat{x} + i\hat{p}}{\sqrt{2}}, \qquad a^\dagger = \frac{\hat{x} - i\hat{p}}{\sqrt{2}}.

The operator aa is the annihilation (or lowering) operator and aa^\dagger the creation (or raising) operator. Because x^\hat{x} and p^\hat{p} are Hermitian, aa^\dagger is genuinely the adjoint of aa — taking the dagger swaps the sign of the ip^i\hat{p} term, as shown. Note that aa is not Hermitian: it is not an observable but a tool for moving between states.

Inverting the definitions expresses position and momentum back in terms of the ladder operators,

x^=a+a2,p^=aai2=i(aa)2,\hat{x} = \frac{a + a^\dagger}{\sqrt{2}}, \qquad \hat{p} = \frac{a - a^\dagger}{i\sqrt{2}} = \frac{i(a^\dagger - a)}{\sqrt{2}},

a substitution we will use repeatedly.

Why this factors the Hamiltonian

Compute the product aaa^\dagger a directly, being careful that x^\hat{x} and p^\hat{p} do not commute:

aa=(x^ip^)(x^+ip^)2=x^2+p^2+i(x^p^p^x^)2=x^2+p^2+i[x^,p^]2.a^\dagger a = \frac{(\hat{x} - i\hat{p})(\hat{x} + i\hat{p})}{2} = \frac{\hat{x}^2 + \hat{p}^2 + i(\hat{x}\hat{p} - \hat{p}\hat{x})}{2} = \frac{\hat{x}^2 + \hat{p}^2 + i[\hat{x}, \hat{p}]}{2}.

Using [x^,p^]=i[\hat{x}, \hat{p}] = i, the cross term becomes ii=1i \cdot i = -1, so

aa=x^2+p^2212=H^12.a^\dagger a = \frac{\hat{x}^2 + \hat{p}^2}{2} - \frac{1}{2} = \hat{H} - \frac{1}{2}.

Rearranging gives the central identity of the whole module:

 H^=aa+12 (in units of ω).\boxed{\ \hat{H} = a^\dagger a + \frac{1}{2}\ } \quad(\text{in units of }\hbar\omega).

The leftover 12\tfrac{1}{2} is the zero-point energy, and it arose entirely from the non-commutativity of x^\hat{x} and p^\hat{p}. Had the operators commuted (the classical case), the factorization would have been exact and the ground-state energy zero.

Try it

This is a numerical exercise — your code should return a number. Using the number-basis action an=nn1a|n\rangle = \sqrt{n}\,|n-1\rangle and an=n+1n+1a^\dagger|n\rangle = \sqrt{n+1}\,|n+1\rangle, build small matrices for aa and aa^\dagger, form aa+12a^\dagger a + \tfrac{1}{2}, and confirm its (0,0)(0,0) entry is the ground-state energy E0=12E_0 = \tfrac{1}{2}.

Run your code to see the quantum state.

Sign in on the full site to ask questions and join the discussion.