A free particle in 1 dimension is described by

H^=p^22m\hat{H} = \frac{\hat{p}^2}{2m}

What are eigenvectors of H^\hat{H}? Momentum basis states -- since

p^∣k⟩=ℏk∣k⟩\hat{p}\ket{k}=\hbar k\ket{k}

then

H^∣k⟩=(ℏk)22m∣k⟩\hat{H}\ket{k} = \frac{(\hbar k)^2}{2m}\ket{k}

Time evolution operator in spectral form yields

U^(t)=e−iH^t/ℏ=∫−∞∞e−iℏk22mt ∣k⟩⟨k∣ dk\hat{U}(t) = e^{-i\hat{H}t/\hbar} = \int_{-\infty}^{\infty} e^{-i\frac{\hbar k^2}{2m}t}\,\ket{k}\bra{k}\,dk

Plane wave packet

ψ(x,0)=eik0x⟶ψ(x,t)=ei(k0x−ω0t),ω0=ℏk022m\psi(x,0) = e^{ik_0 x} \longrightarrow \psi(x,t) = e^{i(k_0 x - \omega_0 t)}, \quad \omega_0 = \frac{\hbar k_0^2}{2m}

Gaussian wave packet

Recall Gaussian Wave Packet. This is the superposition of plane waves. We chose the Gaussian as the wave function ϕ(x)\phi(x) because it saturates the Heisenberg Uncertainty Relation.

ψ(x,0)=N e−(x−x0)24σ2 eik0x\psi(x,0) = N\,e^{-\frac{(x-x_0)^2}{4\sigma^2}}\,e^{ik_0 x}

Perform a FT

ψ~(k,0)=N~ e−σ2(k−k0)2\tilde\psi(k,0) = \tilde N\,e^{-\sigma^2(k-k_0)^2}

Time-evolve each component by eiℏk2t/(2m)e^{i\hbar k^2 t/(2m)}

ψ~(k,t)=N~ e−iℏk22mt e−σ2(k−k0)2\tilde\psi(k,t) = \tilde N\,e^{-i\frac{\hbar k^2}{2m}t}\,e^{-\sigma^2(k-k_0)^2}

Perform an inverse FT

ψ(x,t)=N~ e−[x−xc(t)]24σc(t)2 ei(⋯ )\psi(x,t) = \tilde N\,e^{-\frac{[x - x_c(t)]^2}{4\sigma_c(t)^2}}\,e^{i(\cdots)}

Where i(...)i(...) is some crazy boi

i(⋯ )=ik0x−iω0t+iτ [x−xc(t)]24σ2(1+τ2)−12iarctan⁡(τ)i(\cdots) = i k_0 x - i\omega_0 t + i\frac{\tau\,[x-x_c(t)]^2}{4\sigma^2(1+\tau^2)} - \tfrac{1}{2}i\arctan(\tau)

and

xc(t)=x0+ℏk0mt=x0+v0t,v0=ℏk0mx_c(t) = x_0 + \frac{\hbar k_0}{m}t = x_0 + v_0 t, \qquad v_0 = \frac{\hbar k_0}{m}

and

σc2(t)=σ2+ℏ24m2σ2t2\sigma_c^2(t) = \sigma^2 + \frac{\hbar^2}{4m^2\sigma^2}t^2

Moving Gaussian

If we square the wave function then we can see its probability density

∣ψ(x,t)∣2=∣N~∣2 e−(x−xc(t))22σc(t)2|\psi(x,t)|^2 = |\tilde N|^2\,e^{-\frac{(x - x_c(t))^2}{2\sigma_c(t)^2}}

We can see that the Gaussian probability cloud drifts and broadens. The center xc(t)x_c(t) slides at constant velocity and the width σc(t)\sigma_c(t) grows over time. The shape always remains Gaussian.

This is similar to diffusive brownian motion of a classical particle. A particle buffeted by random collision usually has a spread like:

σclassical(t)∼Dt\sigma_\text{classical}(t) \sim \sqrt{Dt}

but a quantum particle spreads at

σc(t)≈ℏt2mσ\sigma_c(t) \approx \frac{\hbar t}{2m\sigma}

The quantum particle spreads linear in tt, but the brownian one spreads at t\sqrt{t}. This is faster than a classical diffusing cloud

Wave packet spreads faster if the initial position is more localized/confined -- that means σ\sigma is small.