Back

The Variational Method

Most quantum systems cannot be solved exactly, but there is a remarkably simple way to bound the ground-state energy from above: guess a wavefunction. The variational theorem guarantees that for any normalized trial state, the expectation $\langle H\rangle$ is greater than or equal to the true ground-state energy $E_0$, with equality only if your guess happens to be the exact ground state. So you build a trial wavefunction with an adjustable knob, compute the energy, and slide the knob to push $\langle H\rangle$ as low as it will go; the minimum is your best estimate, and you know it is an honest upper bound. This playground does it for the hydrogen atom using a Gaussian trial $e^{-a r^2}$, the same shape that underlies the Gaussian basis sets of every quantum-chemistry code. The catch is visible in the scene: the true ground state $e^{-r}$ has a sharp cusp at the nucleus, a kink that a smooth Gaussian simply cannot reproduce. Because of that mismatch the best Gaussian lands at $\langle H\rangle = -0.424$ Hartree, a respectable but clearly imperfect approximation to the exact $-0.5$. The lower panel plots $\langle H\rangle$ against the trial width $a$: the curve dips to a minimum at $a^\ast = 8/9\pi$ and, no matter how you tune the parameter, never crosses the forbidden line at $E_0$. That floor is the whole content of the variational principle.

Figure 1. The variational method. Left: the trial Gaussian (blue) against the exact 1s state (orange), which has a cusp at the nucleus the Gaussian cannot match. Right: the energy against the exact floor E0 = -0.5 Ha, split into kinetic and potential parts. Bottom: the variational curve (a), minimized at a* = 8/9 pi and never crossing E0. Method: the energy functional (a) = (3/2)a - 2 sqrt(2a/pi). Source: Griffiths, Introduction to Quantum Mechanics, 3rd ed., Ch. 7.
0.283

WHAT TO TRY

  • Sweep the width $a$: the energy $\langle H\rangle$ slides along its curve but never drops below $E_0$.
  • Jump to the best $a^\ast$: the minimum sits at $-0.424$ Ha, the closest a single Gaussian can get.
  • Look at $r=0$: the exact state has a cusp the Gaussian flattens, which is exactly why the bound is not tight.
  • Watch the kinetic and potential split: at the optimum the virial relation $2\langle T\rangle = -\langle V\rangle$ holds.