7.1: The Variational Method
In this section we introduce the powerful and versatile variational method and use it to improve the approximate solutions we found for the helium atom using the independent electron approximation.
Experimental (True) Energy of the Helium Atom
The helium atom has two electrons bound to a nucleus with charge \(Z = 2\). The successive removal of the two electrons can be considered stepwise:
\[\ce{He} \xrightarrow {\textit{I}_1} \ce{He}^+ + e^\xrightarrow {\textit{I}_2}\ce{He}^{++}+2e^\label{7.1.1}\]
The first ionization energy \(I_1\) is the minimum energy required to remove the first electron from helium gas and is experimentally determined:
\[\textit{I}_1=\textit{E}_{1\textit{s}}(\ce{He}) = 24.59\;eV \label{7.1.2a}\]
The second ionization energy, \(I_2\) can experimentally determined, but also can be calculated exactly from the hydrogen atom solutions since \(He^+\) is a hydrogenlike ion with \(Z=2\). Hence, we have
\[\textit{I}_2=\textit{E}_{ 1\textit{s}}(\ce{He}^+)=\dfrac{Z^2}{2n^2}=54.42\mbox{ eV}\label{7.1.2b}\]
The energy of the three separated particles on the right side of Equation \(\ref{7.1.1}\) is (by definition) zero. Therefore the groundstate energy of helium atom is given by \(E_0=(\textit{I}_1+\textit{I}_2)=79.02\mbox{ eV}\). We will attempt to reproduce this true value, as close as possible, by theoretical analysis.
The "Ignorance is Bliss" Approximation
The Hamiltonian for the Helium atom is:
\[\hat{H} = \dfrac{\hbar^2}{2m_e}\bigtriangledown_{el_{1}}^2 \dfrac{\hbar^2}{2m_e}\bigtriangledown_{el_{2}}^2  \dfrac {Ze^2}{4\pi\epsilon_0 r_1}  \dfrac {Ze^2}{4\pi\epsilon_0 r_2} + \dfrac {e^2}{4\pi \epsilon_0 r_{12}}\label{7.1.3}\]
If we simply ignore the \(\dfrac {e^2}{4\pi\epsilon_0 r_{12}}\) electronelectron repulsion term, then electron 1 is separable from electron 2, so that the total wavefunction is a product to two oneelectron wavefunctions:
\[\Psi_{total} = \psi_{el_{1}}\psi_{el_{2}} \label{7.1.4a}\]
or in braket notation
\[  \Psi_{total} \rangle = \hat{H}  \psi_{el_1} \rangle  \psi_{el_2} \rangle \label{7.1.4b}\]
With some operator algebra, something important arises  the one electron energies are additive:
\[\hat{H} \Psi_{total} = (\hat{H}_{el_1} + \hat{H}_{el_2}) \psi_{n\ {el_1}} \psi_{n\ {el_2}} = (E_{n_1} + E_{n_2}) \psi_{n\ {el_1}} \psi_{n\ {el_2}} \label{7.1.5a}\]
or in braket notation
\[\hat{H}  \Psi_{total} \rangle = \hat{H}  \psi_{el_1} \rangle  \psi_{el_2} \rangle = (E_{n_1} + E_{n_2})  \psi_{1} \rangle  \psi_{1} \rangle \label{ 7.1.5b}\]
The energy for a ground state Helium atom (both electrons in lowest state) is then
\[E_{He_{1s}} = \underset{\text{energy of single electron in helium}}{E_{n_1}} + \underset{\text{energy of single electron in helium}}{E_{n_2}} = R\left(\dfrac{Z^2}{1}\right) R \left(\dfrac{Z^2}{1}\right) = 8R \label{7.1.6}\]
where \(R\) is the Rydberg constant (\(13.6 \; eV\) that also maps the lowest energy of the hydrogen atom and \(Z=2\) for the helium nucleus. Experimentally, we find that the total energy for a ground state Helium atom \(E_{He_{1s}} = 5.8066\,R\). The big difference in results is due to the approximation we used to get to the value of \(8\;R\). The "ignorance is bliss" approximation overestimates the energy of the helium atom greatly; this is a poor approximation and we need to address electronelectron repulsion properly (or better at least).
Shielding and Penetration
One way to take electronelectron repulsion into account is to modify the form of the wavefunction. A logical modification is to change the nuclear charge, \(Z\), in the wavefunctions to an effective nuclear charge, from +2 to a smaller value, \(\zeta\) (called zeta). The rationale for making this modification is that one electron partially shields the nuclear charge from the other electron, as shown in Figure \(\PageIndex{1}\).
Figure \(\PageIndex{1}\): Electronelectron shielding leading to a reduced effective nuclear charge. The attractive force of the nucleus on electron 2, \(V(r_2)\), is partially countered by the repulsive force between electron 1 and electron 2, \(V(r_{12})\).
A region of negative charge density between one of the electrons and the +2 nucleus makes the potential energy between them more positive (decreases the attraction between them). We can effect this change mathematically by using \(\zeta < 2\) in the wavefunction expression. If the shielding were complete, then \(\zeta\) would equal 1. If there is no shielding, then \(\zeta = 2\). So a way to take into account the electronelectron interaction is by saying it produces a shielding effect. The shielding is not zero, and it is not complete, so the effective nuclear charge varies between one and two.
Exercise \(\PageIndex{1}\)
Compute the binding energy for He using Equation (916) and a shielded nucleus with \(\zeta\) = 1.5. How does this result, relative to the one obtained with \(Z = 2\), compare to the experimental binding energy in Table \(\PageIndex{1}\)? Does using \(\zeta\) = 1.5 over estimate or under estimate the shielding effect?
In general, a theory should be able to make predictions in advance of knowledge of the experimental result. Consequently, a principle and method for choosing the best value for \(\zeta\) or any other adjustable parameter that is to be optimized in a calculation is needed. The Variational Principle provides the required criterion and method and says that the best value for any variable parameter in an approximate wavefunction is the value that gives the lowest energy for the ground state; i.e., the value that minimizes the energy. The variational method is the procedure that is used to find the lowest energy and the best values for the variable parameters.
A Better Approximation: The Variational Method
Tthe variational method is one way of finding approximations to the lowest energy eigenstate or ground state, and some excited states. This allows calculating approximate wavefunctions and is the variational principle. The method consists in choosing a "trial wavefunction" depending on one or more parameters, and finding the values of these parameters for which the expectation value of the energy is the lowest possible. The wavefunction obtained by fixing the parameters to such values is then an approximation to the ground state wavefunction, and the expectation value of the energy in that state is an upper bound to the ground state energy.
The variational principle means that the expectation value for the binding energy obtained using an approximate wavefunction and the exact Hamiltonian operator will be higher than or equal to the true energy for the system. This idea is really powerful. When implemented, it permits us to find the best approximate wavefunction from a given wavefunction that contains one or more adjustable parameters, called a trial wavefunction. A mathematical statement of the variational principle is
\[ E_{trial} \ge E_{true} \label {7.1.7}\]
where
\[ E_{trial} = \dfrac{ \langle \psi _{trial} \hat {H}  \psi _{trial} \rangle}{\langle \psi _{trial}  \psi _{trial} \rangle} = \dfrac {\int \psi _{trial} ^* \hat {H} \psi _{trial} d \tau}{\int \psi _{trial} ^* \psi _{trial} d\tau } \label {7.1.8}\]
Equation \(\ref{7.1.7}\) is call the variational theorem and states that for a timeindependent Hamiltonian operator, any trial wave function will have an variational energy expectation value that is greater than or equal to the true ground state wave function corresponding to the given Hamiltonian. Because of this, the variational energy is an upper bound to the true ground state energy of a given molecule.The general approach of this method consists in choosing a "trial wavefunction" depending on one or more parameters, and finding the values of these parameters for which the expectation value of the energy is the lowest possible (Figure \(\PageIndex{2}\)).
Figure \(\PageIndex{2}\): Simplified algorithmic flowchart of the Variational Method approximation
The variational energy \(E_{trial}\) is only equal to the true energy \(E_{true}\) when the the corresponding trial wavefunction \(\Psi_{trial}\) is equal to the true wavefunction \(\Psi_{true}\).
Application to the Helium atom Ground State
Often the expectation values (numerator) and normalization integrals (denominator) in Equation \(\ref{7.1.8}\) can be evaluated analytically. For the case of the He atom, the trial wavefunction is the product wavefunction given by Equation \(\ref{7.1.9}\), the adjustable or variable parameter in the trial wavefunction is the effective nuclear charge \(\zeta\) (would be equal to \(Z\) if fully unshielded), and the Hamiltonian is the complete form. When the expectation value for the trial energy is calculated for helium, the result is a function that depends on the adjustable parameter, \(\zeta\).
\[  E_{trial} (\zeta) \rangle = \dfrac {\mu e^4}{4 \epsilon ^2_0 h} \left ( \zeta ^2  \dfrac {27}{8} \zeta \right ) \label {7.1.9}\]
This function is shown in Figure \(\PageIndex{1}\). According to the variation principle, the minimum value of the energy on this graph is the best approximation of the true energy of the system, and the associated value of \(\zeta\) is the best value for the adjustable parameter.
Figure \(\PageIndex{2}\): Graph of trial energies for helium atom as a function of the adjustable parameter \(\zeta\) in Equation \(\ref{7.1.9}\) that represents the effective nuclear charge felt by the electrons.
Using the mathematical function for the energy of a system, the minimum energy with respect to the adjustable parameter can be found by taking the derivative of the energy with respect to that parameter, setting the resulting expression equal to zero, and solving for the parameter, in this case \(\zeta\). This is a standard method in calculus for finding maxima and minima.
Exercise \(\PageIndex{2}\)
Find the value for \(\zeta\) that minimizes the helium binding energy and compare the binding energy to the experimental value. What is the percent error in the calculated value?
When this procedure is carried out for He, we find \(\zeta\) = 1.6875 and the approximate energy we calculate using this approximation method, E_{approx} = 77.483 eV. Table \(\PageIndex{1}\) and Figure \(\PageIndex{2}\) show that a substantial improvement in the accuracy of the computed binding energy is obtained by using shielding to account for the electronelectron interaction. Including the effect of electron shielding in the wavefunction reduces the error in the binding energy to about 2%. This idea is very simple, elegant, and significant.











The improvement we have seen in the total energy calculations using a variable parameter \(\zeta\) indicates that an important contribution of electronelectron interaction or repulsion to the total binding energy arises from the fact that each electron shields the nuclear charge from the other electron. It is reasonable to assume the electrons are independent; i.e., that they move independently, but the shielding must be taken into account in order to finetune the wavefunctions. The inclusion of optimizable parameters in the wavefunction allows us to develop a clear physical image of the consequences of our variation calculation. Calculating energies correctly is important, and it is also important to be able to visualize electron densities for multielectron systems. In the next two sections, we take a temporary break from our consideration of approximation methods in order to examine multielectron wavefunctions more closely.
Contributors
Adapted from "Quantum States of Atoms and Molecules" by David M. Hanson, Erica Harvey, Robert Sweeney, Theresa Julia Zielinski