|q⟩ Bad Qubits

intermediate · Physics · The Harmonic Oscillator: Ladder Operators

The Ground State

The ground state 0|0\rangle is the bottom rung of the ladder — the state that the lowering operator sends to zero. That single condition, a0=0a|0\rangle = 0, determines the entire ground-state wavefunction without solving the Schrödinger equation.

The defining condition

The ladder terminates at the bottom: there is no 1|{-1}\rangle, so

a0=0.a|0\rangle = 0.

This is the algebraic fingerprint of the ground state. Everything about 0|0\rangle follows from it. Its energy, for instance, is immediate:

H^0=ω(aa+12)0=ω(a0+120)=12ω0.\hat{H}|0\rangle = \hbar\omega\left(a^\dagger a + \tfrac{1}{2}\right)|0\rangle = \hbar\omega\left(a^\dagger \cdot 0 + \tfrac{1}{2}|0\rangle\right) = \tfrac{1}{2}\hbar\omega\,|0\rangle.

So E0=12ωE_0 = \tfrac{1}{2}\hbar\omega, the zero-point energy, with no integration required.

Recovering the Gaussian wavefunction

To find ψ0(x)=x0\psi_0(x) = \langle x|0\rangle, write the condition a0=0a|0\rangle = 0 in position space. In natural units a=(x^+ip^)/2a = (\hat{x} + i\hat{p})/\sqrt{2} with p^=id/dx\hat{p} = -i\,d/dx, so

a0=0  12(x+ddx)ψ0(x)=0  dψ0dx=xψ0.a|0\rangle = 0 \ \Longrightarrow\ \frac{1}{\sqrt{2}}\left(x + \frac{d}{dx}\right)\psi_0(x) = 0 \ \Longrightarrow\ \frac{d\psi_0}{dx} = -x\,\psi_0.

This is a first-order differential equation — far simpler than the original second-order Schrödinger equation — and it integrates at once to a Gaussian:

ψ0(x)=(1π)1/4ex2/2,\psi_0(x) = \left(\frac{1}{\pi}\right)^{1/4} e^{-x^2/2},

where the prefactor is fixed by normalization ψ02dx=1\int |\psi_0|^2\,dx = 1. The algebraic method turns finding the ground state into solving one first-order equation.

Spread of the ground state

Because ψ0\psi_0 is a Gaussian centred at the origin, x^=p^=0\langle\hat{x}\rangle = \langle\hat{p}\rangle = 0. Its variances follow from the ladder operators. Using x^=(a+a)/2\hat{x} = (a + a^\dagger)/\sqrt{2} and a0=0a|0\rangle = 0,

x^2=120(a+a)20=120aa0=12.\langle\hat{x}^2\rangle = \tfrac{1}{2}\langle 0|(a + a^\dagger)^2|0\rangle = \tfrac{1}{2}\langle 0|a a^\dagger|0\rangle = \tfrac{1}{2}.

An identical calculation gives p^2=12\langle\hat{p}^2\rangle = \tfrac{1}{2}. Hence the uncertainty product is

ΔxΔp=1212=12(restoring ),\Delta x\,\Delta p = \sqrt{\tfrac{1}{2}}\cdot\sqrt{\tfrac{1}{2}} = \tfrac{1}{2}\hbar \quad(\text{restoring }\hbar),

which saturates the Heisenberg bound ΔxΔp/2\Delta x\,\Delta p \ge \hbar/2. The ground state is a minimum-uncertainty state — the tightest a quantum particle can be squeezed into a parabolic well.

Try it

This is a numerical exercise — your code should return a number. Using x^=(a+a)/2\hat{x} = (a + a^\dagger)/\sqrt{2} and the truncated number-basis matrices, compute the ground-state spread 0x^20\langle 0|\hat{x}^2|0\rangle (which equals 12\tfrac{1}{2}).

Run your code to see the quantum state.

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