The Hartree-Fock (HF) method approximates the many-electron wavefunction as a single Slater determinant -- an antisymmetrized product of one-electron orbitals that automatically satisfies the Pauli exclusion principle. Each electron moves in the mean field of all other electrons, and the orbitals are optimized iteratively: guess orbitals, compute the mean field (Fock operator), solve for new orbitals, repeat until self-consistency. The variational principle guarantees the HF energy is an upper bound to the true energy. Basis sets (STO-3G, 6-31G*, cc-pVDZ, etc.) expand each molecular orbital in a finite set of known functions, and basis set size controls accuracy versus cost. The fundamental limitation is that HF neglects electron correlation -- the instantaneous electron-electron interactions beyond the mean-field approximation -- which typically accounts for ~1% of total energy but can be chemically decisive for bond energies and reaction barriers.
Run HF calculations on small molecules (H2, H2O, HF) using computational chemistry software with progressively larger basis sets. Compare the computed bond lengths and energies to experimental values and to correlated methods, seeing how the correlation energy gap persists regardless of basis set completeness.
The central challenge of quantum chemistry is the many-electron problem. The Schrödinger equation for a molecule with N electrons contains interaction terms between every pair of electrons, making an exact analytical solution impossible for N > 1. The Hartree-Fock method attacks this problem with an elegant simplification: replace the instantaneous electron-electron repulsion with an average, or mean, field. Each electron is treated as if it moves independently in the combined field of the nuclei and the averaged repulsion from all other electrons.
The wavefunction built from this approximation is not just a product of one-electron functions — that would violate the quantum mechanical requirement that the wavefunction change sign when any two electrons are swapped (the antisymmetry principle, which enforces the Pauli exclusion principle). The Slater determinant solves this: it is constructed so that swapping any two rows (i.e., two electrons) changes the determinant's sign, and setting two rows equal makes it zero (two electrons cannot occupy the same state). So the Slater determinant is a compact, elegant way to build antisymmetry into a product of one-electron orbitals.
The orbitals themselves are not known in advance. This leads to the self-consistent field procedure: start with a reasonable guess for the orbitals, compute the mean field (encoded in the Fock operator) that each electron experiences, solve the resulting eigenvalue equations for a new set of orbitals, then use those new orbitals to recompute the Fock operator, and repeat. You keep cycling until the orbitals from one iteration match the orbitals used to build the Fock operator — that is, until the field is self-consistent. The variational principle guarantees that the converged HF energy is an upper bound to the true ground-state energy.
In practice, molecular orbitals are expanded in a basis set — a finite collection of known mathematical functions (typically Gaussian-type orbitals centered on atoms). The choice of basis set determines both accuracy and computational cost. Small basis sets like STO-3G are fast but inaccurate; larger ones like cc-pVTZ are more accurate but expensive. Importantly, making the basis set larger improves accuracy toward the Hartree-Fock limit but cannot recover correlation energy — that is a fundamental limitation of the single-determinant approximation, not a basis set problem.
The correlation energy is the gap between the Hartree-Fock limit energy and the true energy. It represents the instantaneous electron-electron repulsions that the mean-field picture ignores — electrons actually avoid each other in real time, not just on average. This typically accounts for under 1% of total energy but can be chemically decisive for bond dissociation, reaction barriers, and dispersion interactions. Post-HF methods (MP2, coupled cluster, configuration interaction) address this by mixing in excited Slater determinants, and density functional theory takes a different route entirely — but HF remains the conceptual and practical foundation for most of these approaches.