Wavefunction Slater

McLean, A. D., and McLean, R. S. (1981). Roothaan-Hartree-Fock atomic wavefunctions. Slater basis set expansions for A=55-92. Atomic Data and Nuclear Data Tables 26, 197-401. 

QM grew out of studies of blackbody radiation and of the photoelectric effect. Besides QM, radioactivity and relativity contributed to the transition from classical to modem physics. The classical Rutherford nuclear atom, the Bohr atom, and the Schrodinger wave-mechanical atom are discussed. Hybridization, wavefunctions, Slater determinants and other basic concepts are explained. For obtaining eigenvectors and eigenvalues from the secular equations the elegant and simple matrix diagonalization method is explained and used. All the necessary mathematics is explained. 

The Slater wavefunction differs from the Hartree function not only in being composed of spin orbitals rather than just spatial orbitals, but also in the fact that it is not a simple product of one-electron functions, but rather a determinant (Section 4.3.3) whose elements are these functions. To construct a Slater wavefunction (Slater determinant) for a closed-shell species (the only kind we consider in any detail here), we use each of the occupied spatial orbitals to make two spin orbitals, by multiplying the spatial orbital by a and, separately, by jl. The spin orbitals are then filled with the available electrons. An example should make the procedure clear (Fig. 5.2). Suppose we wish to write a Slater determinant for a four-electron... 

Combination of the two AOs with the spin functions a and p produces four spin orbitals each containing one electron. This means that only one wavefunction (Slater determinant) can be written down the ground state of Be (or any other Group 2 atom) is non-degenerate. 

In the configuration interaction method, the wavefunction is expanded in a basis of Slater determinants. The Hartree-Fock determinant, noted is taken as a reference zeroth-order wavefunction. Slater determinants corresponding to excited configurations are generated by swapping occupied MOs V a with virtual (unoccupied) MOs ipr and can be classified with respect to the number of excited electrons. Singly excited Slater determinants are noted 3, doubly excited Slater determinants triply excited Slater determinants X so on. The configuration interaction wavefunction then reads... 

For atoms with more than two electrons, the wavefunctions must become more elaborate to satisfy the symmetrization principle. However, John Slater developed a general method for reliably generating many-electron spin-spatial wavefunctions, antisymmetric with respect to P21 (exchange of the electron labels 1 and 2), for any number of electrons. We call these wavefunctions Slater determinants, because they are obtained by taking the determinant of a matrix of possible one-electron wavefunctions. For example, for ground state He, there are two possible one-electron spin-spatial wavefunctions for each electron Isa and lsj8. We set up a 2 X 2 matrix in which each row corresponds to a different electron and each column to a different wavefunction ... 

The normalisation factor is assumed. It is often convenient to indicate the spin of each electron in the determinant this is done by writing a bar when the spin part is P (spin down) a function without a bar indicates an a spin (spin up). Thus, the following are all commonly used ways to write the Slater determinantal wavefunction for the beryllium atom (which has the electronic configuration ls 2s ) ... 

VVc can now see why the normalisation factor of the Slater determinantal wavefunction is I v/N . If each determinant contains N terms then the product of two Slater determinants, ldeU rminant][determinant], contains (N ) terms. However, if the spin orbitals form an oi lhonormal set then oidy products of identical terms from the determinant will be nonzero when integrated over all space. We Ccm illustrate this with the three-electron example, k ljiiiidering just the first two terms in the expansion we obtain the following ... 

Ihe one-electron orbitals are commonly called basis functions and often correspond to he atomic orbitals. We will label the basis functions with the Greek letters n, v, A and a. n the case of Equation (2.144) there are K basis functions and we should therefore xpect to derive a total of K molecular orbitals (although not all of these will necessarily 3e occupied by electrons). The smallest number of basis functions for a molecular system vill be that which can just accommodate all the electrons in the molecule. More sophisti- ated calculations use more basis functions than a minimal set. At the Hartree-Fock limit he energy of the system can be reduced no further by the addition of any more basis unctions however, it may be possible to lower the energy below the Hartree-Fock limit ay using a functional form of the wavefunction that is more extensive than the single Slater determinant. 

In accordance with the variation theorem we require the set of coefficients that gives the lowest-energy wavefunction, and some scheme for changing the coefficients to derive that wavefunction. For a given basis set and a given functional form of the wavefunction (i.e. a Slater determinant) the best set of coefficients is that for which the energy is a minimum, at which point... 

In order to calculate higher-order wavefunctions we need to establish the form of the perturbation, f. This is the difference between the real Hamiltonian and the zeroth-order Hamiltonian, Remember that the Slater determinant description, based on an orbital picture of the molecule, is only an approximation. The true Hamiltonian is equal to the sum of the nuclear attraction terms and electron repulsion terms ... 

To express, in terms of Slater determinants, the wavefunctions corresponding to each of the states in each of the levels, one proceeds as follows ... 

The single Slater determinant wavefunction (properly spin and symmetry adapted) is the starting point of the most common mean field potential. It is also the origin of the molecular orbital concept. 

Corrections to the mean-field model are needed to describe the instantaneous Coulombic interactions among the electrons. This is achieved by including more than one Slater determinant in the wavefunction. 

The first-order MPPT wavefunction can be evaluated in terms of Slater determinants that are excited relative to the SCF reference function k. Realizing again that the perturbation coupling matrix elements I>k H i> are non-zero only for doubly excited CSF s, and denoting such doubly excited i by a,b m,n the first-order... 

Such a determinantal wavefunction is called a Slater determinant, after Slater (1929), and you should appreciate that a... 

I don t mean that such a wavefunction is necessarily very accurate you saw a minute ago that the LCAO treatment of dihydrogen is rather poor. I mean that, in principle, a Slater determinant has the correct spatial and spin symmetry to represent an electronic state. It very often happens that we have to take combinations of Slater determinants in order to make progress for example, the first excited states of dihydrogen caimot be represented adequately by a single Slater determinant such as... 

As computational facilities improve, electronic wavefunctions tend to become more and more complicated. A configuration interaction (Cl) calculation on a medium-sized molecule might be a linear combination of a million Slater determinants, and it is very easy to lose sight of the chemistry and the chemical intuition , to say nothing of the visualization of the results. Such wavefunctions seem to give no simple physical picture of the electron distribution, and so we must seek to find ways of extracting information that is chemically useful. 

I will refer to the Hartree model from time to time in the text. Hartree s energies were in poor agreement with experiment. With the benefit of hindsight he should have allowed for indistinguishability and the Pauli principle. This was Fock s contribution to the field he wrote the wavefunction as what we would now recognize as a Slater determinant. Such a wavefunction automatically satisfies the Pauli principle. 

Don t confuse the state wavefunction with a molecular orbital we might well want to build the state wavefunction, which describes all the 16 electrons, from molecular orbitals each of which describe a single electron. But the two are not the same. We would have to find some suitable one-electron wavefunctions and then combine them into a slater determinant in order to take account of the Pauli principle. 

In Chapter 6, I discussed the open-shell HF-LCAO model. 1 considered the simple case where we had ti doubly occupied orbitals and 2 orbitals all singly occupied by parallel spin electrons. The ground-state wavefunction was a single Slater determinant. I explained that it was possible to derive an expression for the electronic energy... 

The total wavefunction, , is an antisymmetrized product of the one-electron functions i/q (a Slater determinant). The i/tj are called one-electron functions since they depend on the coordinates of only one electron this approximation is embedded in all MO methods. The effects that are missing when this approximation is used go under the general name of electron correlation. 

To ensure this, the-many-body wavefunction can be written as a Slater determinant of one particle wavefunctions - this is the Hartree Fock method. The drawbacks of this method are that it is computationally demanding and does not include the many-body correlation effects. 

The HPHF wavefunction for an 2n electron system, in a gronnd state of S qnantum number, even or odd, is written as a linear combination of only two DODS Slater determinants, built up with spinorbitals which minimize the total energy [1-2] ... 

Also, an alternative formulation of equation (17) can be conceived if one wants to distinguish between ground state, monoexcitations, biexcitations,. .. and so on. Such a possibility is symbolized in the following Cl wavefunction expression for n electrons, constructed as to include Slater determinants up to the p-th (pp) unoccupied ones l9klk=i,ni Then, the Cl wavefunction is written in this case as the linear combination ... 

Density functions can be obtained up to any order from the manipulation of the Slater determinant functions alone as defined in section 5.1 or from any of the linear combinations defined in section 5.2. Density functions of any order can be constructed by means of Lowdin or McWeeny descriptions [17], being the diagonal elements of the so called m-th order density matrix, as was named by Lowdin the whole set of possible density functions. For a system of n electrons the n-th order density function is constructed from the square modulus of any n-electron wavefunction attached to the n-electron system somehow. 

Using a unnormalized n-electron Slater determinant D(j) as system wavefunction, constructed as discussed in section 5.1, then one can write the n-th order density function p " (j) as ... 

The calculated perturbed wavefunctions have been rewritten in terms of a combination of normalized Slater orbitals in real form. Ref. 6 gives a detailed illustration for... 

---