Hartree–Fock Theory: 100 Self-Contained Facts

A comprehensive collection of facts about Hartree–Fock theory, organized by topic.

I. Foundations of the Electronic Structure Problem

  1. Hartree–Fock theory is an approximate method for solving the time-independent electronic Schrödinger equation for atoms and molecules.
  2. Hartree–Fock assumes the Born–Oppenheimer approximation, meaning nuclei are treated as fixed while solving for electronic motion.
  3. The exact electronic wavefunction for N electrons depends on 3N spatial coordinates and N spin coordinates.
  4. The electronic Hamiltonian contains electron kinetic energy, electron–nucleus attraction, and electron–electron Coulomb repulsion.
  5. The electron–electron repulsion term 1/r_ij couples all electrons and prevents exact separation of variables.
  6. The presence of electron–electron repulsion makes the many-electron Schrödinger equation non-separable.
  7. Solving the exact electronic problem scales exponentially with the number of electrons.
  8. Hartree–Fock reduces the complexity of the many-electron problem through approximation.
  9. The goal of Hartree–Fock is to find the best single-determinant approximation to the ground state.
  10. Hartree–Fock becomes exact for one-electron systems because electron–electron repulsion is absent.

II. Slater Determinant and Antisymmetry

  1. Hartree–Fock approximates the many-electron wavefunction using a single Slater determinant.
  2. A Slater determinant is constructed from one-electron spin orbitals arranged in determinant form.
  3. A spin orbital depends on spatial coordinates and a spin coordinate.
  4. The determinant structure guarantees antisymmetry under exchange of any two electrons.
  5. Antisymmetry means that exchanging two electrons changes the sign of the wavefunction.
  6. The antisymmetry property enforces the Pauli exclusion principle.
  7. If two electrons occupy the same spin orbital, the Slater determinant equals zero.
  8. A normalized Slater determinant requires orthonormal spin orbitals.
  9. The determinant contains N! terms for N electrons.
  10. The use of a single determinant neglects configurations beyond that determinant.

III. Variational Principle and Energy Minimization

  1. Hartree–Fock theory is based on the variational principle of quantum mechanics.
  2. The variational principle states that any approximate normalized wavefunction gives an energy greater than or equal to the exact ground-state energy.
  3. Hartree–Fock minimizes the total electronic energy with respect to variations in the spin orbitals.
  4. The minimization is performed under the constraint that orbitals remain orthonormal.
  5. Orthonormal orbitals are mutually orthogonal and individually normalized to one.
  6. Lagrange multipliers are introduced to enforce orbital orthonormality during minimization.
  7. Applying the stationary condition yields the Hartree–Fock equations.
  8. The Hartree–Fock equations resemble eigenvalue equations.
  9. The Hartree–Fock equations are nonlinear because the operator depends on the orbitals being solved for.
  10. The solutions of the Hartree–Fock equations are called canonical orbitals.

IV. The Fock Operator

  1. The effective one-electron operator in Hartree–Fock theory is called the Fock operator.
  2. The Fock operator contains a one-electron term plus Coulomb and exchange contributions.
  3. The one-electron term includes electron kinetic energy and nuclear attraction.
  4. The Coulomb operator represents classical electrostatic repulsion between electrons.
  5. The Coulomb operator depends on the average electron density.
  6. The exchange operator arises from antisymmetry of the Slater determinant.
  7. The exchange operator has no classical analogue.
  8. The exchange operator is nonlocal because it depends on orbital values at multiple spatial points.
  9. Exchange interactions occur only between electrons of the same spin.
  10. Opposite-spin electrons do not contribute exchange terms.

V. Self-Consistent Field (SCF) Procedure

  1. The Fock operator depends on the occupied orbitals.
  2. The occupied orbitals are obtained by solving the Fock operator eigenvalue equation.
  3. This circular dependence makes Hartree–Fock a self-consistent theory.
  4. Hartree–Fock equations are solved iteratively in the self-consistent field (SCF) method.
  5. The SCF procedure begins with an initial guess for orbitals or electron density.
  6. From this guess, a Fock operator is constructed.
  7. The Fock operator is diagonalized to obtain new orbitals.
  8. A new electron density is constructed from the updated orbitals.
  9. The SCF process repeats until convergence is achieved.
  10. Convergence means successive iterations produce negligible changes in energy and density.

VI. Matrix Formulation (Roothaan Equations)

  1. In practical calculations, orbitals are expanded in a finite basis set of known functions.
  2. Gaussian-type orbitals are commonly used as basis functions in molecular Hartree–Fock calculations.
  3. Expanding orbitals in a basis converts differential equations into matrix equations.
  4. The Hartree–Fock equations in a basis are called the Roothaan equations.
  5. The Roothaan equations have the matrix form F C = S C ε.
  6. F is the Fock matrix constructed from one- and two-electron integrals.
  7. S is the overlap matrix containing overlaps between basis functions.
  8. C is the matrix of orbital expansion coefficients.
  9. ε is a diagonal matrix containing orbital energies.
  10. The Roothaan equations form a generalized eigenvalue problem.

VII. Energy Expression and Orbital Energies

  1. The Hartree–Fock total energy is not equal to the simple sum of orbital energies.
  2. Orbital energies double-count electron–electron interactions.
  3. The correct Hartree–Fock energy subtracts half of the Coulomb and exchange contributions.
  4. Only occupied orbitals contribute directly to the Hartree–Fock energy.
  5. Virtual orbitals do not contribute directly to the Hartree–Fock total energy.
  6. The Hartree–Fock energy is variational within a given basis set.
  7. Increasing the basis set size lowers or maintains the Hartree–Fock energy.
  8. The difference between exact energy and Hartree–Fock energy is called correlation energy.
  9. Correlation energy arises from electron motions that are correlated beyond mean-field exchange.
  10. Hartree–Fock includes exchange exactly within a single determinant but neglects dynamic correlation.

VIII. Physical Interpretation

  1. Hartree–Fock treats each electron as moving in an average field created by all other electrons.
  2. The average field smooths out instantaneous electron–electron interactions.
  3. Hartree–Fock orbitals often resemble chemically intuitive bonding and antibonding orbitals.
  4. Koopmans' theorem states that the negative of an occupied orbital energy approximates the ionization energy.
  5. Koopmans' theorem assumes that orbitals do not relax after electron removal.
  6. Hartree–Fock can predict spin polarization in open-shell systems.
  7. Restricted Hartree–Fock (RHF) uses the same spatial orbital for paired electrons of opposite spin.
  8. Unrestricted Hartree–Fock (UHF) allows different spatial orbitals for alpha and beta electrons.
  9. UHF may suffer from spin contamination, meaning the wavefunction is not an eigenfunction of total spin.
  10. Restricted open-shell Hartree–Fock (ROHF) is designed for open-shell systems while preserving some spin symmetry.

IX. Computational Scaling and Practical Aspects

  1. The computational cost of Hartree–Fock scales approximately as N^4 with respect to basis size.
  2. The N^4 scaling arises from evaluating two-electron integrals involving four basis functions.
  3. Two-electron integrals are the most computationally expensive part of Hartree–Fock.
  4. Integral screening techniques reduce the number of integrals that must be evaluated.
  5. Density fitting approximates electron repulsion integrals to accelerate calculations.
  6. Direct SCF methods compute integrals on the fly instead of storing them.
  7. DIIS (Direct Inversion in the Iterative Subspace) is commonly used to accelerate SCF convergence.
  8. Convergence difficulties may arise for near-degenerate or strongly correlated systems.
  9. Hartree–Fock typically converges reliably for closed-shell molecules.
  10. Modern Hartree–Fock implementations rely on optimized linear algebra libraries.

X. Relation to Correlated Methods

  1. Hartree–Fock provides a reference wavefunction for post-Hartree–Fock correlation methods.
  2. Configuration Interaction (CI) improves upon Hartree–Fock by including multiple determinants.
  3. Full Configuration Interaction (FCI) includes all possible determinants within a basis and is exact within that basis.
  4. Møller–Plesset perturbation theory treats electron correlation as a perturbation to the Hartree–Fock reference.
  5. Coupled-cluster theory builds correlated wavefunctions using exponential operators applied to the Hartree–Fock determinant.
  6. Density Functional Theory (DFT) replaces the wavefunction with the electron density as the central variable.
  7. Hybrid DFT functionals mix Hartree–Fock exchange with density functional exchange.
  8. Hartree–Fock becomes exact in the absence of electron–electron interaction.
  9. Even in a complete basis set, a single-determinant Hartree–Fock wavefunction cannot describe strong correlation.
  10. Hartree–Fock is the foundational mean-field approximation underlying most modern electronic structure methods.