|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, a∣0⟩=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

a∣0⟩=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⟩=ℏω(a†a+12)∣0⟩=ℏω(a†⋅0+12∣0⟩)=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)=⟨x∣0⟩\psi_0(x) = \langle x|0\rangle, write the condition a∣0⟩=0a|0\rangle = 0 in position space. In natural units a=(x^+ip^)/2a = (\hat{x} + i\hat{p})/\sqrt{2} with p^=−i d/dx\hat{p} = -i\,d/dx, so

a∣0⟩=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/4e−x2/2,\psi_0(x) = \left(\frac{1}{\pi}\right)^{1/4} e^{-x^2/2},

where the prefactor is fixed by normalization ∫∣ψ0∣2 dx=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 a∣0⟩=0a|0\rangle = 0,

⟨x^2⟩=12⟨0∣(a+a†)2∣0⟩=12⟨0∣aa†∣0⟩=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=12⋅12=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 ⟨0∣x^2∣0⟩\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.