4.2 Theoretical Background

4.2.2 Hartree-Fock Theory

(June 30, 2021)

As with much of the theory underlying modern quantum chemistry, the HF approximation was developed shortly after publication of the Schrödinger equation, but remained a qualitative theory until the advent of the computer. Although the HF approximation tends to yield qualitative chemical accuracy, rather than quantitative information, and is generally inferior to many of the DFT approaches available, it remains as a useful tool in the quantum chemist’s toolkit. In particular, for organic chemistry, HF predictions of molecular structure are very useful.

Consider once more the Roothaan-Hall equations, Eq. (4.11), or the Pople-Nesbet equations, Eq. (4.13), which can be traced back to Eq. (4.8), in which the effective potential υeff depends on the SCF methodology. In a restricted HF (RHF) formalism, the effective potential can be written as

υeff=aN/2[2J^a(1)-K^a(1)]-A=1MZAr1A (4.15)

where the Coulomb and exchange operators are defined as

J^a(1)=ψa(2)1r12ψa(2)𝑑𝐫2 (4.16)

and

K^a(1)ψi(1)=[ψa(2)1r12ψi(2)𝑑𝐫2]ψa(1) (4.17)

respectively. By introducing an atomic orbital basis, we obtain Fock matrix elements

Fμν=Hμνcore+Jμν-Kμν (4.18)

where the core Hamiltonian matrix elements

Hμνcore=Tμν+Vμν (4.19)

consist of kinetic energy elements

Tμν=ϕμ(𝐫)(-12^2)ϕν(𝐫)𝑑𝐫 (4.20)

and nuclear attraction elements

Vμν=ϕμ(𝐫)(-AZA|𝐑A-𝐫|)ϕν(𝐫)𝑑𝐫 (4.21)

The Coulomb and exchange elements are given by

Jμν=λσPλσ(μν|λσ) (4.22)

and

Kμν=12λσPλσ(μλ|νσ) (4.23)

respectively, where the density matrix elements are

Pμν=2a=1N/2CμaCνa (4.24)

and the two electron integrals are

(μν|λσ)=ϕμ(𝐫1)ϕν(𝐫1)(1r12)ϕλ(𝐫2)ϕσ(𝐫𝟐)𝑑𝐫1𝑑𝐫2. (4.25)

Note:  The formation and utilization of two-electron integrals is a topic central to the overall performance of SCF methodologies. The performance of the SCF methods in new quantum chemistry software programs can be quickly estimated simply by considering the quality of their atomic orbital integrals packages. See Appendix B for details of Q-Chem’s AOInts package.

Substituting the matrix element in Eq. (4.18) back into the Roothaan-Hall equations, Eq. (4.11), and iterating until self-consistency is achieved will yield the RHF energy and wave function. Alternatively, one could have adopted the unrestricted form of the wave function by defining separate α and β density matrices:

Pμνα=a=1nαCμaαCνaαPμνβ=a=1nβCμaβCνaβ (4.26)

The total electron density matrix 𝐏=𝐏α+𝐏β. The unrestricted α Fock matrix,

Fμνα=Hμνcore+Jμν-Kμνα, (4.27)

differs from the restricted one only in the exchange contributions, where the α exchange matrix elements are given by

Kμνα=λNσNPλσα(μλ|νσ) (4.28)