9.5: Singleelectron Wavefunctions and Basis Functions
Finding the most useful singleelectron wavefunctions to serve as building blocks for a multielectron wavefunction is one of the main challenges in finding approximate solutions to the multielectron Schrödinger Equation. The functions must be different for different atoms because the nuclear charge and number of electrons are different. The attraction of an electron for the nucleus depends on the nuclear charge, and the electronelectron interaction depends upon the number of electrons.
As we saw in our initial approximation methods, the most straightforward place to start in finding reasonable singleelectron wavefunctions for multielectron atoms is with the atomic orbitals produced in the quantum treatment of hydrogen, the socalled “hydrogenic” spinorbitals. These traditional atomic orbitals, with a few modifications, give quite reasonable calculated results and are still in wide use for conceptually understanding multielectron atoms. In this section and in Chapter 10 we will explore some of the many other singleelectron functions that also can be used as atomic orbitals.
Hydrogenic spinorbitals used as components of multielectron systems are identified in the same way as they are for the hydrogen atom. Each spinorbital consists of a spatial wavefunction, specified by the quantum numbers (n, \(l , m_l\)) and denoted ls, 2s, 2p, 3s, 3p, 3d, etc, multiplied by a spin function, specified by the quantum number \(m_s\) and denoted \(\alpha\) or \(\beta\). In our initial approximation methods, we ignored the spin components of the hydrogenic orbitals, but they must be considered in order to develop a complete description of multielectron systems. The subscript on the argument of the spatial function reveals which electron is being described (\(r_1\) is a vector that refers to the coordinates of electron 1, for example.) No argument is given for the spin function. An example of a spinorbital for electron 2 in a \(3p_z\) orbital:
\[  \varphi _{3p_z} \alpha (r_2) \rangle = \varphi _{3,1,0}(r_2) \alpha \label {9.5.1}\]
In the alternative shorthand notation for this spinorbital shown below, the coordinates for electron 2 in the spatial function are abbreviated simply by the number “2,” and the spatial function is represented by “\(3p_z\)” rather than "\(\varphi _{3,1,0}\)". The argument “2” given for the spin function refers to the unknown spin variable for electron 2. Many slight variations on these shorthand forms are in use in this and other texts, so flexibility and careful reading are important.
\[  \varphi _{3p_z}\alpha (2) \rangle = 3p_z (2) \alpha (2) \label {9.5.2}\]
In this chapter we will continue the trend of moving away from writing specific mathematical functions and toward a more symbolic, condensed representation. Your understanding of the material in this and future chapters requires that you keep in mind the form and properties of the specific functions denoted by the symbols used in each equation.
Exercise \(\PageIndex{1}\)
Write the full mathematical form of \(\varphi _{3pz\alpha}\) using as much explicit functional detail as possible.
The basic mathematical functions and thus the general shapes and angular momenta for hydrogenic orbitals are the same as those for hydrogen orbitals. The differences between atomic orbitals for the hydrogen atom and those used as components in the wavefunctions for multielectron systems lie in the radial parts of the wavefunctions and in the energies. Specifically, the differences arise from the replacement of the nuclear charge Z in the radial parts of the wavefunctions by an adjustable parameter \(\zeta\) that is allowed to vary in approximation calculations in order to model the interactions between the electrons. We discussed such a procedure for helium The Variational Method previously. The result is that electrons in orbitals with different values for the angular momentum quantum number, \(l\), have different energies. Figure \(\PageIndex{1}\) shows the results of a quantum mechanical calculation on argon in which the degeneracy of the 2s and 2p orbitals is found to be removed, as is the degeneracy of the 3s, 3p, and 3d orbitals.
Figure \(\PageIndex{1}\): Ordering of energy levels for Ar. Energy level differences are not to scale.
The energy of each electron now depends not only on its principle quantum number, \(n\), but also on its angular momentum quantum number, \(l\).
The presence of \(\zeta\) in the radial portions of the wavefunctions also means that the electron probability distributions associated with hydrogenic atomic orbitals in multielectron systems are different from the exact atomic orbitals for hydrogen. Figure \(\PageIndex{2}\) compares the radial distribution functions for an electron in a 1s orbital of hydrogen (the ground state), a 2s orbital in hydrogen (an excited configuration of hydrogen) and a 1s orbital in helium that is described by the best variational value of \(\zeta\). Our use of hydrogenlike orbitals in quantum mechanical calculations for multielectron atoms helps us to interpret our results for multielectron atoms in terms of the properties of a system we can solve exactly.
Figure \(\PageIndex{2}\): Radial distribution functions for 1s of hydrogen (red, \(\zeta\) = 1), 2s of hydrogen (blue, \(\zeta\) = 1) and 1s of helium (black, \(\zeta\) = 1.6875).
Exercise \(\PageIndex{2}\)
Analyze Figure \(\PageIndex{2}\) and write a paragraph about what you can discern about the relative sizes of ground state hydrogen, excited state hydrogen and ground state helium atoms.
While they provide useful stepping off points for understanding computational results, nothing requires us to use the hydrogenic functions as the building blocks for multielectrons wavefunctions. In practice, the radial part of the hydrogenic atomic orbital presents a computational difficulty because the radial function has nodes, positive and negative lobes, and steep variations that make accurate evaluation of integrals by a computer slow. Consequently other types of functions are generally used in building multielectron functions. These usually are related to the hydrogenic orbitals to aid in the analysis of molecular electronic structure. For example, Slatertype atomic orbitals (STO’s), designated below as \(S_{nlm} (r, \theta , \varphi )\), avoid the difficulties imposed by the hydrogenic functions. The STO’s, named after their creator, John Slater, were the first alternative functions that were used extensively in computations. STO’s do not have any radial nodes, but still contain a variational parameter \(\zeta\) (zeta), that corresponds to the effective nuclear charge in the hydrogenic orbitals. In Equation \(\ref{936}\) and elsewhere in this chapter, the distance, \(r\), is measured in units of the Bohr radius, \(a_0\).
\[ S_{nlm} (r, \theta , \varphi ) = \dfrac {(2 \zeta )^{n+1/2}}{[(2n)!]^{1/2}} r^{n1} e^{\zeta r } Y^m_l (\theta , \varphi ) \label {936}\]
Exercise \(\PageIndex{3}\)
 Write the radial parts of the 1s, 2s, and 2p atomic orbitals for hydrogen.
 Write the radial parts of the n = 1 and n = 2 Slater–type orbitals (STO).
 Check that the above five functions are normalized.
 Graph these five functions, measuring r in units of the Bohr radius.
 Graph the radial probability densities for these orbitals. Put the hydrogen orbital and the corresponding STO on the same graph so they can be compared easily.
 Adjust the zeta parameter \(\zeta\) in each case to give the best match of the radial probability density for the STO with that of the corresponding hydrogen orbital.
 Comment on the similarities and differences between the hydrogen orbitals and the STOs and the corresponding radial probability densities.
Linear Variational Method
An alternative approach to the general problem of introducing variational parameters into wavefunctions is the construction of a singleelectron wavefunction as a linear combination of other functions. For hydrogen, the radial function decays, or decreases in amplitude, exponentially as the distance from the nucleus increases. For helium and other multielectron atoms, the radial dependence of the total probability density does not fall off as a simple exponential with increasing distance from the nucleus as it does for hydrogen. More complex singleelectron functions therefore are needed in order to model the effects of electronelectron interactions on the total radial distribution function. One way to obtain more appropriate singleelectron functions is to use a sum of exponential functions in place of the hydrogenic spinorbitals.
An example of such a wavefunction created from a sum or linear combination of exponential functions is written as
\[ \varphi _{1s} (r_1) = \sum _j c_j e^{\zeta _j r_j /a_o} \label{937}\]
The linear combination permits weighting of the different exponentials through the adjustable coefficients (cj) for each term in the sum. Each exponential term has a different rate of decay through the zetaparameter \(\zeta _j\). The exponential functions in Equation \(\ref{937}\) are called basis functions. Basis functions are the functions used in linear combinations to produce the singleelectron orbitals that in turn combine to create the product multielectron wavefunctions. Originally the most popular basis functions used were the STO’s, but today STO’s are not used in most quantum chemistry calculations. However, they are often the functions to which more computationally efficient basis functions are fitted.
Physically, the \(\zeta _j\) parameters account for the effective nuclear charge (often denoted with \(Z_{eff}\). The use of several zeta values in the linear combination essentially allows the effective nuclear charge to vary with the distance of an electron from the nucleus. This variation makes sense physically. When an electron is close to the nucleus, the effective nuclear charge should be close to the actual nuclear charge. When the electron is far from the nucleus, the effective nuclear charge should be much smaller. See Slater's rules for a ruleofthumb approach to evaluate \(Z_{eff}\) values.
A term in Equation \(\ref{937}\) with a small \(\zeta\) will decay slowly with distance from the nucleus. A term with a large \(\zeta\) will decay rapidly with distance and not contribute at large distances. The need for such a linear combination of exponentials is a consequence of the electronelectron repulsion and its effect of screening the nucleus for each electron due to the presence of the other electrons.
Exercise \(\PageIndex{4}\)
Make plots of \(\varphi\) in Equation \(\ref{937}\) using three equally weighted terms with \(\zeta\) = 1.0, 2.0, and 5.0. Also plot each term separately.
Computational procedures in which an exponential parameter like \(\zeta\) is varied are more precisely called the Nonlinear Variational Method because the variational parameter is part of the wavefunction and the change in the function and energy caused by a change in the parameter is not linear. The optimum values for the zeta parameters in any particular calculation are determined by doing a variational calculation for each orbital to minimize the groundstate energy. When this calculation involves a nonlinear variational calculation for the zetas, it requires a large amount of computer time. The use of the variational method to find values for the coefficients, \(\{c_j\}\), in the linear combination given by Equation \(\ref{937}\) above is called the Linear Variational Method because the singleelectron function whose energy is to be minimized (in this case \(\varphi _{1s}\)) depends linearly on the coefficients. Although the idea is the same, it usually is much easier to implement the linear variational method in practice.
Nonlinear variational calculations are extremely costly in terms of computer time because each time a zeta parameter is changed, all of the integrals need to be recalculated. In the linear variation, where only the coefficients in a linear combination are varied, the basis functions and the integrals do not change. Consequently, an optimum set of zeta parameters were chosen from variational calculations on many small multielectron systems, and these values, which are given in Table \(\PageIndex{1}\), generally can be used in the STOs for other and larger systems.


































Exercise \(\PageIndex{5}\)
Compare the value \(\zeta _{1s}\) = 1.24 in Table \(\PageIndex{1}\) for hydrogen with the value you obtained in Exercise \(\PageIndex{3}\). and comment on possible reasons for any difference. Why are the zeta values larger for 1s than for 2s and 2p orbitals? Why do the \(\zeta _{1s}\) values increase by essentially one unit for each element from He to Ne while the increase for the \(\zeta _{2s, 2p}\) values is much smaller?
The discussion above gives us some new ideas about how to write flexible, useful singleelectron wavefunctions that can be used to construct multielectron wavefunctions for variational calculations. Singleelectron functions built from the basis function approach are flexible because they have several adjustable parameters, and useful because the adjustable parameters still have clear physical interpretations. Such functions will be needed in the HartreeFock method discussed elsewhere.
Contributors
 Adapted from "Quantum States of Atoms and Molecules" by David M. Hanson, Erica Harvey, Robert Sweeney, Theresa Julia Zielinski