| Hartree–Fock method | |
|---|---|
| Name | Hartree–Fock method |
| Field | Quantum mechanics |
| Developer | Douglas Hartree; Vladimir Fock |
| Introduced | 1928–1930 |
| Related | Configuration interaction, Density functional theory, Møller–Plesset perturbation theory |
Hartree–Fock method
The Hartree–Fock method is an approximate method for solving the many-body problem in quantum mechanics by reducing an interacting system of fermions to an effective single-particle description using antisymmetric trial wavefunctions. It yields a set of self-consistent integro-differential equations for molecular or atomic orbitals and forms the foundation for many post‑Hartree–Fock electronic structure methods in quantum chemistry and condensed matter physics.
The Hartree–Fock (HF) approach was developed to incorporate the effects of the Pauli exclusion principle for systems of electrons while treating electron–electron interactions at a mean-field level. Early formulations by Douglas Hartree (self-consistent field) and Vladimir Fock introduced the use of antisymmetrized product states (Slater determinant) to enforce fermionic antisymmetry. HF replaces the true many-electron wave function with a single Slater determinant built from one-electron spin orbitals, producing an average or "mean" potential experienced by each electron called the Hartree–Fock potential. The method captures exchange effects exactly for a given single-determinant ansatz but neglects dynamic electron correlation beyond exchange, motivating subsequent correlated methods such as Configuration interaction and Coupled cluster theory.
In the HF formalism the electronic Hamiltonian for N electrons, H = Σ_i h(i) + 1/2 Σ_{i≠j} g(i,j), is projected onto a trial Slater determinant |Φ⟩. Applying the variational principle yields the Hartree–Fock equations: F[ρ] φ_p = ε_p φ_p, where F is the Fock operator depending on the one‑particle density matrix ρ, φ_p are canonical spin orbitals and ε_p are orbital energies. The Fock operator comprises the one-electron operator h (kinetic + nuclear attraction) plus the Coulomb operator J and the nonlocal exchange operator K: F = h + J − K. Practical implementations discretize orbitals in a finite basis set, commonly Gaussian-type orbitals (GTOs) or Slater-type orbitals (STOs), converting the integro-differential equations into a matrix eigenvalue problem (Roothaan or Pople equations) solved as a generalized eigenvalue problem FC = SCε where S is the overlap matrix.
HF solutions are found via the self-consistent field (SCF) iterative procedure. Starting from an initial guess for the density (e.g., superposition of atomic densities or Hückel method-like guesses), one constructs the Fock matrix, solves the Roothaan equations for molecular orbital coefficients, updates the density matrix, and repeats until convergence of energy and density. Convergence techniques widely used include density mixing, level shifting, Direct inversion in the iterative subspace (DIIS), and damping. SCF may encounter issues such as variational collapse, multiple local minima, or slow convergence near degeneracies; remedies include orbital localization, symmetry constraints, and use of more robust initial guesses from Hartree or Kohn–Sham calculations.
HF has multiple spin and symmetry adaptations: Restricted Hartree–Fock (RHF) enforces identical spatial orbitals for paired electrons in closed-shell systems; Unrestricted Hartree–Fock (UHF) allows different spatial orbitals for different spin components and can describe open-shell systems and spin polarization but may suffer from spin contamination. Restricted open-shell HF (ROHF) provides a compromise for certain radicals. Post‑HF extensions include Møller–Plesset perturbation theory (MPn) and Coupled cluster (CC) methods that recover correlation energy missing in HF. Hybrid methods combine HF exchange with Density functional theory (DFT) exchange–correlation functionals (e.g., B3LYP, PBE0), creating HF-DFT hybrids that improve thermochemistry and barrier heights by mixing exact exchange with approximate functionals. Other variants include time‑dependent HF (TDHF) for excited states and generalized Hartree–Fock (GHF) that relaxes spin restrictions.
HF is a standard starting point for electronic structure calculations of atoms, molecules, and periodic systems in quantum chemistry and solid state physics. It provides molecular geometries, orbital energies, and qualitative descriptions of bonding used in studies at institutions such as Harvard University, Massachusetts Institute of Technology, University of Cambridge, and national laboratories (e.g., Lawrence Berkeley National Laboratory). HF underpins basis-set development (e.g., 6-31G basis families) and serves as reference wavefunctions for correlated methods in computational packages like Gaussian, GAMESS, NWChem, Quantum ESPRESSO (for periodic HF), and Molpro. In condensed matter, HF approximations are related to mean-field theories such as the Hartree–Fock approximation to the Hubbard model and inform studies of ferromagnetism, Wigner crystals, and Fermi liquid behavior.
While HF captures exchange exactly within a single determinant, it neglects dynamic correlation arising from instantaneous electron–electron interactions, leading to errors in total energies, bond dissociation, dispersion forces, and reaction barriers. Remedies include perturbative corrections (MP2), variational configuration interaction (CISD, CISDT), coupled-cluster (CCSD, CCSD(T)), and multi-reference methods (CASSCF) for strong correlation. HF's lack of electron correlation also impacts properties such as excitation spectra and response functions, often necessitating time‑dependent DFT (TDDFT) or equation‑of‑motion coupled-cluster approaches for accurate predictions. Conceptually, HF remains valuable as a variational, computationally efficient mean-field baseline and as the zeroth‑order reference in many-body perturbation theory and Green's function approaches like GW approximation.
Category:Quantum chemistry Category:Computational chemistry