Modern computational quantum chemistry has developed largely from applications of the Hartree–Fock method to atoms and molecules [1–3]. A simple representation of a manyelectron atom is given by a Slater determinant constructed from occupied spinorbitals: , where spinorbitals are products of the form . A Slater determinant automatically satisfies the Pauli exclusion principle, requiring that spinorbitals be, at most, singly occupied. The Hamiltonian for an electron atom, in hartree atomic units, is given by: . The corresponding approximation to the total energy is then given by , where , and are, respectively, the core, coulomb and exchange integrals. This Demonstration carries out a simplified Hartree–Fock computation on the ground states of the atoms He to Ne, to 10. Assume a single Slater determinant, thus possibly shortchanging openshell configurations. The computation involves only , and orbitals, approximated by modified Slatertype orbitals: , , . These are normalized and mutually orthogonal, which simplifies the computation. The orbital parameters are chosen so as to minimize the total energy, in accordance with the variational principle. You can vary these using the sliders. Actually, since the possible values of vary over inconveniently large ranges, the sliders vary internal parameters that determine them. The values of are shown in the graphic. The blue bar shows the exact energy of the atomic ground state (actually, the exact nonrelativistic value) in hartrees. The red bar shows the calculated energy, given the displayed values of . You should try to adjust the three parameters to get the lowest possible energy. You can never reach the exact energy, however, since a single Hartree–Fock determinant fails to account for correlation energy, which involves instantaneous electronelectron interactions.
Following are the optimized values of the energies and orbital parameters, found using the builtin Mathematica function FindMinimum: Z = 2, {2.84766, {α 1.6875}} Z = 3, {7.42617, {α 2.69716, β 0.721027}} Z = 4, {14.5413, {α 3.71867, β 1.10264}} Z = 5, {24.4442, {α 4.72562, β 1.50386, γ 1.2137}} Z = 6, {37.4687, {α 5.72987, β 1.90153, γ 1.53991}} Z = 7, {53.9194, {α 6.73255, β 2.29803, γ 1.86317}} Z = 8, {74.1007, {α 7.73422, β 2.69406, γ 2.18487}} Z = 9, {98.3167, {α 8.73521, β 3.08988, γ 2.50566}} Z = 10, {126.872, {α 9.73572, β 3.48558, γ 2.82585}} These can be slightly improved by hand. [1] S. M. Blinder, Introduction to Quantum Mechanics, Amsterdam: Elsevier, 2004 section 9.7. [2] A. Szabo and N. S. Ostlund, Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory, Mineola, NY: Dover, 1996. [3] S. M. Blinder, "Introduction to the Hartree–Fock Method," in Mathematical Physics in Theoretical Chemistry (S. M. Blinder and J. E. House, eds.), Elsevier, 2018, in press.
