This treatment of the electrons and their orbitals is referred to as the hartree-level of theory. As stated above, when screened hydrogenic AOs are used to approximate the o, and ox orbitals, the resultant E, values do not produce accurate predictions. For example, the negative of E, should approximate the ionization energy for removal of an electron from the Ao. Such ionization potentials(IP s)can be measured, and the measured values do not agree well with the theoretical values when a crude screening approximation is made for the ao s 2. The Laco-Expansion To improve upon the use of screened hydrogenic AOs, it is most common to approximate each of the Hartree AOs(oki as a linear combination of so-called basis AOs x 中=Cmx using what is termed the linear-combination-of-atomic-orbitals (LCAO)expansion. In this equation, the expansion coefficients C) are the variables that are to be determined by solving the Schrodinger equation 中=E1中 After substituting the LCAO expansion for o, into this Schrodinger equaiton, multiplying 6
6 This treatment of the electrons and their orbitals is referred to as the Hartree-level of theory. As stated above, when screened hydrogenic AOs are used to approximate the fJ and fK orbitals, the resultant eJ values do not produce accurate predictions. For example, the negative of eJ should approximate the ionization energy for removal of an electron from the AO fJ . Such ionization potentials (IP s) can be measured, and the measured values do not agree well with the theoretical values when a crude screening approximation is made for the AO s. 2. The LACO-Expansion To improve upon the use of screened hydrogenic AOs, it is most common to approximate each of the Hartree AOs {fK} as a linear combination of so-called basis AOs {cm}: fJ = SmCJ,m cm . using what is termed the linear-combination-of-atomic-orbitals (LCAO) expansion. In this equation, the expansion coefficients {CJ,m} are the variables that are to be determined by solving the Schrödinger equation he fJ = eJ fJ . After substituting the LCAO expansion for fJ into this Schrödinger equaiton, multiplying
on the left by one of the basis AOs %y, and then integrating over the coordinates of the electron in o), one obtains Σ<xlx>Cl=E1E<xJx>C This is a matrix eigenvalue equation in which the E, and (Clu appear as eigenvalues and eigenvectors. The matrices <xl hel x and xl x> are called the Hamiltonian and overlap matrices, respectively. An explicit expression for the former is obtained by introducing the earlier definition of he <xlhlx>=<x-t/2m V-Ix>+<x Ze/r u> +2kCknCKx<xr)xn(r)errD|x(r)x:(r) An important thing to notice about the form of the matrix Hartree equations is that to compute the Hamiltonian matrix, one must know the LCAO coefficients CK, of the orbitals which the electrons occupy On the other hand these lcao coefficients are supposed to be found by solving the Hartree matrix eigenvalue equations. This paradox leads to the need to solve these equations iteratively in a so-called self-consistent field (SCF)technique. In the SCF process, one inputs an initial approximation to the ( CKy) coefficients This then allows one to form the hamiltonian matrix defined above The Hartree matrix equations Eu <xul x> Cu=E, Eu xl x>Cu are then solved for new(CKy,) coefficients and for the orbital energies(Ex). The new LCAO coefficients
7 on the left by one of the basis AOs cn , and then integrating over the coordinates of the electron in fJ , one obtains Sm <cn | he | cm> CJ,m = eJ Sm <cn | cm> CJ,m . This is a matrix eigenvalue equation in which the eJ and {CJ,m} appear as eigenvalues and eigenvectors. The matrices <cn | he | cm> and <cn | cm> are called the Hamiltonian and overlap matrices, respectively. An explicit expression for the former is obtained by introducing the earlier definition of he : <cn | he | cm> = <cn | – h2 /2m Ñ 2 |cm> + <cn | -Ze2 /|r |cm > + SK CK,h CK,g <cn (r) ch (r’) |(e2 /|r-r’|) | cm (r) cg (r’)>. An important thing to notice about the form of the matrix Hartree equations is that to compute the Hamiltonian matrix, one must know the LCAO coefficients {CK,g} of the orbitals which the electrons occupy. On the other hand, these LCAO coefficients are supposed to be found by solving the Hartree matrix eigenvalue equations. This paradox leads to the need to solve these equations iteratively in a so-called self-consistent field (SCF) technique. In the SCF process, one inputs an initial approximation to the {CK,g} coefficients. This then allows one to form the Hamiltonian matrix defined above. The Hartree matrix equations Sm <cn | he | cm> CJ,m = eJ Sm <cn | cm> CJ,m are then solved for “new” {CK,g} coefficients and for the orbital energies {eK}. The new LCAO coefficients
of those orbitals that are occupied are then used to form a"new"Hamiltonian matrix after which the Hartree equations are again solved for another generation of LCAO coefficients and orbital energies. This process is continued until the orbital energies and LCAO coefficients obtained in successive iterations do not differ appreciably. Upon such convergence, one says that a self-consistent field has been realized because the CK,i coefficients are used to form a Coulomb field potential that details the electron-electron Interactions 3.Ao Basis sets STOs and gtos As noted above, it is possible to use the screened hydrogenic orbitals as the xu) However, much effort has been expended at developing alternative sets of functions to use as basis orbitals. The result of this effort has been to produce two kinds of functions hat currently are widely used The basis orbitals commonly used in the lCao process fall into two primary Slater-type orbitals(STOs)Xn, 1,m(r, 0. Q)=Nn. 1, m, Y1m(0. )rn-Ie-sr are haracterized by quantum numbers n, I, and m and exponents(which characterize the orbitals radial'size')s The symbol Nn, l, m, denotes the normalization constant 2. Cartesian Gaussian-type orbitals(GTOs)Xab.c(r, 0,D)=Nab.c a xayb ze exp(-ar4) are characterized by quantum numbers a, b, and c, which detail the angular shape and direction of the orbital, and exponents a which govern the radial ' size For both types of AOs, the coordinates r, 0, and o refer to the position of the 8
8 of those orbitals that are occupied are then used to form a “new” Hamiltonian matrix, after which the Hartree equations are again solved for another generation of LCAO coefficients and orbital energies. This process is continued until the orbital energies and LCAO coefficients obtained in successive iterations do not differ appreciably. Upon such convergence, one says that a self-consistent field has been realized because the {CK,g} coefficients are used to form a Coulomb field potential that details the electron-electron interactions. 3. AO Basis Sets a. STOs and GTOs As noted above, it is possible to use the screened hydrogenic orbitals as the {cm}. However, much effort has been expended at developing alternative sets of functions to use as basis orbitals. The result of this effort has been to produce two kinds of functions that currently are widely used. The basis orbitals commonly used in the LCAO process fall into two primary classes: 1. Slater-type orbitals (STOs) cn,l,m (r,q,f) = Nn,l,m,z Yl,m (q,f) rn-1 e-zr are characterized by quantum numbers n, l, and m and exponents (which characterize the orbital’s radial 'size' ) z. The symbol Nn,l,m,z denotes the normalization constant. 2. Cartesian Gaussian-type orbitals (GTOs) ca,b,c (r,q,f) = N'a,b,c,a xa yb zc exp(-ar2), are characterized by quantum numbers a, b, and c, which detail the angular shape and direction of the orbital, and exponents a which govern the radial 'size'. For both types of AOs, the coordinates r, q, and f refer to the position of the
electron relative to a set of axes attached to the nucleus on which the basis orbital is located. Note that Slater-type orbitals(STOS) are similar to hydrogenic orbitals in the region close to the nucleus. Specifically, they have a non-zero slope near the nucleus (i.e d/dr(exp(-Cr))=0=-s). In contrast, GTOS, have zero slope near r=0 because d/dr(exp(-ar))20=0. We say that STOs display a cusp"at r=0 that is characteristic of the hydrogenic solutions, whereas GtOs do not Although STOs have the proper'cusp' behavior near nuclei, they are used primarily for atomic and linear-molecule calculations because the multi-center integrals <X(1)X(2)e/rr-r2llx(1)x( 2)> which arise in polyatomic-molecule calculations(we will discuss these intergrals later in this Chapter) can not efficiently be evalusted when STOs are employed. In contrast, such integrals can routinely be computed when GTOs are used. This fundamental advantage of gtos has lead to the dominance of these functions in molecular quantum chemistry To overcome the primary weakness of GTO functions(i.e, their radial derivatives vanish at the nucleus), it is common to combine two, three or more gtOs, with combination coefficients which are fixed and not treated as lcao parameters, into new functions called contracted GTOs or CGTOs. Typically, a series of radially tight medium, and loose GTOs are multiplied by contraction coefficients and summed to produce a CGto which approximates the proper'cusp' at the nuclear center(although no such combination of GTOs can exactly produce such a cusp because each Gto has zero slope at r=0) Although most calculations on molecules are now performed using Gaussian orbitals, it should be noted that other basis sets can be used as long as they span enough 9
9 electron relative to a set of axes attached to the nucleus on which the basis orbital is located. Note that Slater-type orbitals (STO's) are similar to hydrogenic orbitals in the region close to the nucleus. Specifically, they have a non-zero slope near the nucleus (i.e., d/dr(exp(-zr))r=0 = -z). In contrast, GTOs, have zero slope near r=0 because d/dr(exp(-ar 2 ))r=0 = 0. We say that STOs display a “cusp” at r=0 that is characteristic of the hydrogenic solutions, whereas GTOs do not. Although STOs have the proper 'cusp' behavior near nuclei, they are used primarily for atomic and linear-molecule calculations because the multi-center integrals <cm (1) ck (2)|e2 /|r1 -r2 || cn (1) cg (2)> which arise in polyatomic-molecule calculations (we will discuss these intergrals later in this Chapter) can not efficiently be evalusted when STOs are employed. In contrast, such integrals can routinely be computed when GTOs are used. This fundamental advantage of GTOs has lead to the dominance of these functions in molecular quantum chemistry. To overcome the primary weakness of GTO functions (i.e., their radial derivatives vanish at the nucleus), it is common to combine two, three, or more GTOs, with combination coefficients which are fixed and not treated as LCAO parameters, into new functions called contracted GTOs or CGTOs. Typically, a series of radially tight, medium, and loose GTOs are multiplied by contraction coefficients and summed to produce a CGTO which approximates the proper 'cusp' at the nuclear center (although no such combination of GTOs can exactly produce such a cusp because each GTO has zero slope at r = 0). Although most calculations on molecules are now performed using Gaussian orbitals, it should be noted that other basis sets can be used as long as they span enough
of the regions of space(radial and angular where significant electron density resides. In fact, it is possible to use plane wave orbitals of the form x(r, 0, d)=N exp[i(kx r sine cosd +k, r sine sind+k, r cose), where N is a normalization constant and k ky, and k, are quantum numbers detailing the momenta of the orbital along the x, y, and z Cartesian directions. The advantage to using such "simple" orbitals is that the integrals one must perform are much easier to handle with such functions. The disadvantage is that one must use many such functions to accurately describe sharply peaked charge distributions of for example, inner-shell core orbitals Much effort has been devoted to developing and tabulating in widely available locations sets of STO or GtO basis orbitals for main-group elements and transition metals. This ongoing effort is aimed at providing standard basis set libraries which 1. Yield predictable chemical accuracy in the resultant energies 2. Are cost effective to use in practical calculations 3. Are relatively transferable so that a given atom s basis is flexible enough to be used fo that atom in various bonding environments(e.g, hybridization and degree of ionization) b. The fundamental core and valence basis In constructing an atomic orbital basis, one can choose from among several classes of functions. First, the size and nature of the primary core and valence basis must be specified. Within this category, the following choices are common a minimal basis in which the number of cto orbitals is equal to the number of core and valence atomic orbitals in the atom 2. A double-zeta(dz)basis in which twice as many CGTOs are used as there are core
10 of the regions of space (radial and angular) where significant electron density resides. In fact, it is possible to use plane wave orbitals of the form c (r,q,f) = N exp[i(kx r sinq cosf + ky r sinq sinf + kz r cosq)], where N is a normalization constant and kx , ky , and kz are quantum numbers detailing the momenta of the orbital along the x, y, and z Cartesian directions. The advantage to using such “simple” orbitals is that the integrals one must perform are much easier to handle with such functions. The disadvantage is that one must use many such functions to accurately describe sharply peaked charge distributions of, for example, inner-shell core orbitals. Much effort has been devoted to developing and tabulating in widely available locations sets of STO or GTO basis orbitals for main-group elements and transition metals. This ongoing effort is aimed at providing standard basis set libraries which: 1. Yield predictable chemical accuracy in the resultant energies. 2. Are cost effective to use in practical calculations. 3. Are relatively transferable so that a given atom's basis is flexible enough to be used for that atom in various bonding environments (e.g., hybridization and degree of ionization). b. The Fundamental Core and Valence Basis In constructing an atomic orbital basis, one can choose from among several classes of functions. First, the size and nature of the primary core and valence basis must be specified. Within this category, the following choices are common: 1. A minimal basis in which the number of CGTO orbitals is equal to the number of core and valence atomic orbitals in the atom. 2. A double-zeta (DZ) basis in which twice as many CGTOs are used as there are core