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=kk\hat{p}\ket{k}=\hbar k\ket{k}

then

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

Time evolution operator in spectral form yields

U^(t)=eiH^t/=eik22mtkkdk\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)=Ne(xx0)24σ2eik0x\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(kk0)2\tilde\psi(k,0) = \tilde N\,e^{-\sigma^2(k-k_0)^2}

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

ψ~(k,t)=N~eik22mteσ2(kk0)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[xxc(t)]24σc(t)2ei()\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()=ik0xiω0t+iτ[xxc(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~2e(xxc(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.