| Hartree–Fock | |
|---|---|
| Name | Hartree–Fock method |
| Caption | Self-consistent field schematic |
| Developer | Douglas Hartree, Vladimir Fock |
| Introduced | 1928–1930s |
| Field | Quantum mechanics; Quantum chemistry |
| Application | Electronic structure theory; atomic and molecular calculations |
Hartree–Fock
Hartree–Fock is an approximate method for determining the wave function and energy of a quantum many-body system in a stationary state. It reduces the many-electron Schrödinger equation to a set of self-consistent single-particle equations by approximating the full wave function as a single antisymmetric Slater determinant of one-electron orbitals. Hartree–Fock underpins much of modern quantum chemistry and provides a baseline for more accurate correlated methods.
The Hartree–Fock approach emerged from independent developments by Douglas Hartree (self-consistent field method) and Vladimir Fock (incorporating antisymmetry via determinants) in the late 1920s and 1930s. It synthesized ideas from Erwin Schrödinger wave mechanics and Wolfgang Pauli's antisymmetry requirement. Early applications addressed atoms and small molecules, influencing computational programs developed at institutions such as Harvard University, Cambridge University, and later national laboratories including Lawrence Berkeley National Laboratory and Argonne National Laboratory. The method's historical significance lies in enabling tractable electronic structure calculations before the advent of modern high-performance computing and in establishing terminology still used in methods like density functional theory and post-Hartree–Fock techniques.
Hartree–Fock approximates an N-electron wave function by a single Slater determinant constructed from N orthonormal spin-orbitals. The variational principle applied to this ansatz yields the Hartree–Fock equations: a set of coupled integro-differential equations for the spin-orbitals. In second-quantized form the method minimizes the expectation value of the non-relativistic electronic Hamiltonian (kinetic energy, electron–nucleus attraction, electron–electron repulsion) subject to orthonormality constraints, often expressed using the Fock operator F̂. The Fock operator contains a one-electron term (core Hamiltonian) and mean-field potentials: the Coulomb operator J and the non-local exchange operator K arising from antisymmetry. Solutions are obtained iteratively to self-consistency, called the self-consistent field (SCF) procedure. Key theoretical concepts include Koopmans' theorem, which relates orbital energies to ionization potentials in the frozen-orbital approximation, and Brillouin's theorem governing single excitations from the HF ground state.
Practical Hartree–Fock calculations use finite basis sets (e.g., Gaussian or Slater-type orbital) to convert integro-differential equations to matrix form, yielding the Roothaan–Hall equations for closed-shell systems and the Pople–Nesbet equations for open-shell cases. Efficient implementations rely on algorithms for two-electron integral evaluation (direct integral, density fitting/auxiliary basis), SCF convergence accelerators (DIIS, level shifting), and symmetry exploitation (point groups from molecular symmetry). Major quantum chemistry packages implementing Hartree–Fock include Gaussian, GAMESS, NWChem, Psi4, ORCA, and TURBOMOLE. High-performance implementations leverage parallel computing on systems provided by national supercomputing centers and use libraries for linear algebra (e.g., BLAS, LAPACK).
Choice of basis set strongly influences Hartree–Fock accuracy. Common families include minimal bases (STO-nG), split-valence (e.g., 6-31G), correlation-consistent sets (Dunning's cc-pVnZ), and polarized/diffuse augmented variants (e.g., aug-cc-pVTZ). Approximations such as Restricted HF (RHF), Unrestricted HF (UHF), and Restricted Open-shell HF (ROHF) handle spin constraints and open-shell species. Extensions incorporate relativistic effects via scalar-relativistic Hamiltonians (Douglas–Kroll–Hess) or four-component Dirac–Hartree–Fock for heavy elements, and account for environmental effects through embedding schemes like QM/MM used with packages such as CHARMM or AMBER. Techniques like density fitting (RI) and Cholesky decomposition accelerate integral computations.
Hartree–Fock provides molecular geometries, vibrational frequencies (after Hessian evaluation), and qualitative electronic structures for atoms, molecules, and solids. It serves as the reference state for correlated methods including Møller–Plesset perturbation theory (MP2), configuration interaction (CI), and coupled cluster (CC) methods, and supplies orbitals for time-dependent Hartree–Fock (TDHF) and response theories to compute excitation spectra. In solid-state physics, Hartree–Fock variants inform band structure calculations and hybrid methods combining HF exchange with density functional theory (e.g., B3LYP, PBE0). Applications span spectroscopy, reaction mechanism studies, and benchmarks for force-field parameterization used in molecular dynamics.
Hartree–Fock neglects dynamic and static electron correlation beyond mean-field exchange, leading to characteristic errors: incorrect dissociation limits, overestimation of total energies, and missing dispersion interactions. Remedies involve post-Hartree–Fock correlation methods: perturbative MP2, multi-reference approaches (CASSCF), Configuration Interaction (CI), and Coupled Cluster with Single, Double, and perturbative Triple excitations (CCSD(T)), often regarded as a "gold standard". Alternative strategies employ Density Functional Theory for more favorable cost–accuracy tradeoffs or explicitly correlated methods (R12/F12) to accelerate basis set convergence. The choice among these depends on system size, desired accuracy, and available computational resources at facilities such as National Institutes of Health or university computing clusters.
Category:Quantum chemistry Category:Computational chemistry