13.1: Variational Principle

Suppose that we wish to solve the time-independent Schrödinger equation \[H\,\psi = E\,\psi,\] where \(H\) is a known (presumably complicated) time-independent Hamiltonian. Let \(\psi\) be a properly normalized trial solution to the previous equation. The variational principle states, quite simply, that the ground-state energy, \(E_0\), is always less than or equal to the expectation value of \(H\) calculated with the trial wavefunction: that is,

\[E_0 \leq \langle\psi|H|\psi\rangle.\] Thus, by varying \(\psi\) until the expectation value of \(H\) is minimized, we can obtain approximations to the wavefunction and the energy of the ground-state.

Let us prove the variational principle. Suppose that the \(\psi_n\) and the \(E_n\) are the true eigenstates and eigenvalues of \(H\): that is,

\[\label{e14.3} H\,\psi_n = E_n\,\psi_n.\]

Furthermore, let

\[\label{e14.4} E_0 < E_1 < E_2 < \cdots,\] so that \(\psi_0\) is the ground-state, \(\psi_1\) the first excited state, et cetera. The \(\psi_n\) are assumed to be orthonormal: that is,

\[\label{e14.5} \langle \psi_n|\psi_m\rangle = \delta_{nm}.\] If our trial wavefunction \(\psi\) is properly normalized then we can write \[\psi = \sum_n c_n\,\psi_n,\] where

\[\label{e14.7} \sum_n |c_n|^{\,2} = 1.\] Now, the expectation value of \(H\), calculated with \(\psi\), takes the form

\[\begin{aligned} \langle\psi|H|\psi\rangle & = \left.\left\langle \sum_n c_n\,\psi_n\right| H\left|\sum_m\,c_m\,\psi_m\right\rangle\right. = \sum_{n,m} c_n^{\,\ast}\,c_m\,\langle \psi_n|H|\psi_m\rangle\nonumber\\[0.5ex] &= \sum_n\,c_n^{\,\ast}\,c_m\,E_m\,\langle \psi_n|\psi_m\rangle= \sum_n E_n\,|c_n|^{\,2},\end{aligned}\]

where use has been made of Equations \ref{e14.3} and \ref{e14.5}. So, we can write

\[\langle \psi|H|\psi\rangle = |c_0|^{\,2}\,E_0 + \sum_{n>0} |c_n|^{\,2}\,E_n.\]

However, Equation \ref{e14.7} can be rearranged to give

\[|c_0|^{\,2} = 1-\sum_{n>0}|c_n|^{\,2}.\]

Combining the previous two equations, we obtain

\[\langle \psi|H|\psi\rangle = E_0 + \sum_{n>0} |c_n|^{\,2}\,(E_n-E_0).\]

The second term on the right-hand side of the previous expression is positive definite, because \(E_n-E_0>0\) for all \(n>0\) (Equation \ref{e14.4}). Hence, we obtain the desired result

\[\langle \psi|H|\psi\rangle \geq E_0.\]

Excited States

Suppose that we have found a good approximation, \(\tilde{\psi}_0\), to the ground-state wavefunction. If \(\psi\) is a normalized trial wavefunction that is orthogonal to \(\tilde{\psi}_0\) (i.e., \(\langle \psi|\tilde{\psi}_0\rangle=0\)) then, by repeating the previous analysis, we can easily demonstrate that

\[\langle \psi |H|\psi\rangle \geq E_1.\]

Thus, by varying \(\psi\) until the expectation value of \(H\) is minimized, we can obtain approximations to the wavefunction and the energy of the first excited state. Obviously, we can continue this process until we have approximations to all of the stationary eigenstates. Note, however, that the errors are clearly cumulative in this method, so that any approximations to highly excited states are unlikely to be very accurate. For this reason, the variational method is generally only used to calculate the ground-state and first few excited states of complicated quantum systems.

Contributors and Attributions

• Richard Fitzpatrick (Professor of Physics, The University of Texas at Austin)

  \( \newcommand {\ltapp} {\stackrel {_{\normalsize<}}{_{\normalsize \sim}}}\) \(\newcommand {\gtapp} {\stackrel {_{\normalsize>}}{_{\normalsize \sim}}}\) \(\newcommand {\btau}{\mbox{\boldmath$\tau$}}\) \(\newcommand {\bmu}{\mbox{\boldmath$\mu$}}\) \(\newcommand {\bsigma}{\mbox{\boldmath$\sigma$}}\) \(\newcommand {\bOmega}{\mbox{\boldmath$\Omega$}}\) \(\newcommand {\bomega}{\mbox{\boldmath$\omega$}}\) \(\newcommand {\bepsilon}{\mbox{\boldmath$\epsilon$}}\)