The Fock Operator

Introduction

Hartree–Fock turns a many-electron variational problem into coupled one-electron-like equations. The operator in those equations is the Fock operator. It contains the electron's kinetic energy and attraction to fixed nuclei, plus terms generated by occupied electrons. Separating direct Coulomb from exchange reveals what the mean field captures and why the equations must be solved self-consistently.

Core explanation

At fixed nuclear geometry, the one-electron core Hamiltonian h includes electronic kinetic energy and attraction to each nucleus. For electron coordinate r, one may write h = −(ħ²/2m e)∇² − Σ A Z Ae²/(4πε₀ r−R A ) in SI-style notation for point nuclei. Electron–electron repulsion is absent from h; it enters the effective Fock operator through occupied orbitals. Nucleus–nucleus repulsion is a constant for a fixed geometry and is added to the total molecular energy rather than acting as an electron-coordinate operator.

In spin-orbital language, a schematic Fock operator is f = h + Σ j^occupied(J j−K j), with the understanding that exchange vanishes between orthogonal opposite-spin functions and conventions differ in whether self terms are explicitly included. J j describes direct Coulomb repulsion from the average density of occupied spin orbital j. Acting on a test orbital, it resembles multiplication by an electrostatic potential obtained by integrating that density against 1/r 12. It is a mean-field term; the instantaneous position of electron j is not explicitly correlated with the test electron's position.

K j is the exchange operator. It arises from cross terms in the antisymmetrised determinant and depends on the test function over space, so it is nonlocal rather than an ordinary classical potential at one point. Its nonzero contribution requires compatible spin. For an electron interacting with its own spin orbital, the corresponding direct and exchange contributions cancel, preventing unphysical self-repulsion in the exact Hartree–Fock operator construction. Opposite-spin electrons still repel through J even though their exchange integral vanishes.

For a closed-shell restricted spatial-orbital formulation, each occupied spatial orbital j contains one α and one β electron. The Fock operator for a spatial orbital is commonly written F = h + Σ j^occupied spatial(2J j−K j). The factor two in front of J counts both spin occupancies, while one exchange term reflects same-spin exchange. This formula is not contradictory to the spin-orbital form; it packages spin sums differently. A student should state which convention is being used before comparing factors.

The occupied orbitals determine J and K, but the Fock operator determines the orbitals. Thus F is not a fixed operator known in advance from nuclear positions alone. Starting orbitals produce an initial F, whose eigenvectors yield updated orbitals and a new density. Iteration continues until the occupied density and operator agree within a convergence tolerance. Multiple self-consistent solutions can sometimes exist, and the lowest stable solution under a chosen ansatz must be distinguished from a merely converged one.

Fock eigenvalues are orbital energies within a mean-field model. The total energy cannot be obtained by simply adding all occupied eigenvalues because each electron pair's interaction appears in the one-electron equations more than once. A separate energy expression subtracts the appropriate double counting and includes nuclear repulsion. This matters for comparing structures, reaction energies and computational outputs.

The Fock operator depends on the chosen determinant and basis. A restricted closed-shell determinant, an unrestricted open-shell determinant or a different orbital occupancy leads to a different mean field. Changing the basis changes how accurately orbital shapes can respond. Exchange is treated exactly within a single determinant, but dynamical and static correlation beyond one determinant remain outside ordinary Hartree–Fock.

Step-by-step reasoning

Write the core h first and identify which physical terms it contains. Add direct Coulomb J from occupied density and exchange K from antisymmetry, stating a spin-orbital or closed-shell spatial convention. Check that opposite-spin repulsion remains and same-spin self-interaction cancels. Then explain that F must be rebuilt from the orbitals it produces, leading to self-consistency.

Visual explanation

Draw a test electron in the field of nuclei and a diffuse cloud of occupied electron density. Label the nuclear attraction and kinetic contribution h, the cloud's direct repulsion J and a separate arrow for exchange K tied to same-spin orbital overlap. A loop from occupied orbitals to J/K and back to new orbitals depicts self-consistency.

Real-world analogy

A person navigating through a crowd responds to an average crowd pattern, but that pattern depends on how everyone chooses to move. Iterating predictions until the crowd pattern and individual paths agree resembles self-consistency. Exchange has no classical crowd analogue; it comes from indistinguishable electron antisymmetry.

Real-world example

In an atom or molecule with a closed-shell reference, software constructs Coulomb and exchange matrices from the current density and solves for improved orbitals. The relative orbital shapes can help interpret bonding, while the total energy is computed with the proper pair-interaction correction. A failure to converge may require a better initial guess or algorithm, not necessarily a change in physical chemistry.

Why?

Why is exchange nonlocal? Its action on a test function involves an integral containing that function at another electron coordinate, multiplied by an occupied orbital. It cannot generally be represented as a simple number depending only on the test electron's position. This mathematical form follows from exchanging electron labels in the determinant.

Common misconception

Exchange is not a classical force that acts only because electrons have like electrical charge; all electrons have the same charge, while exchange requires matching spin structure. Another error is to interpret 2J−K as two physical repulsions plus one negative particle. It is the algebra of closed-shell occupancy and antisymmetry within a mean-field operator.

Worked example

Consider a one-electron system in a spin-orbital Hartree–Fock formulation. If the occupied orbital itself appears in the formal sum, its direct self term J self is matched by its exchange self term K self when acting on that same orbital, so J self−K self cancels. The remaining orbital equation contains the core h, as it should for one electron. A method leaving a net Coulomb self-repulsion would be physically inconsistent for this simple case.

Quick check

1. Does h include electron–electron repulsion? Answer: No. It contains electron kinetic energy and attraction to fixed nuclei; occupied-electron effects enter J and K. 2. Why is the Fock operator rebuilt during iteration? Answer: Its Coulomb and exchange terms depend on the occupied orbitals, which themselves are obtained from the operator.

Exam focus

State whether you use spin-orbital J−K or closed-shell spatial 2J−K notation. Explain direct versus exchange and why total energy is not an uncorrected orbital-eigenvalue sum. Mention that nuclear repulsion is included in the final fixed-geometry energy, not in the one-electron core operator acting on r.

Advanced insight

Exchange is exact for the chosen single determinant, yet Hartree–Fock's mean-field form can still give an inaccurate electron-pair distribution. Correlated methods alter the many-electron state beyond one determinant. Thus the phrase exact exchange does not mean exact total electronic energy or absence of electron correlation error.

Summary

The Fock operator combines a one-electron core Hamiltonian with occupied-orbital Coulomb and exchange contributions. Coulomb represents average electrostatic repulsion; exchange arises from antisymmetry and is nonlocal and spin selective. Because the operator depends on its own occupied solutions, Hartree–Fock requires self-consistent iteration and a separate total-energy expression.

Practice questions

1. What does the factor 2 in closed-shell 2J j−K j count? Answer: Two opposite-spin electrons occupy each spatial orbital j and both contribute direct Coulomb density, while exchange applies only between same-spin components. 2. Why is it wrong to say opposite-spin electrons do not interact in Hartree–Fock? Answer: Their exchange contribution vanishes due to spin orthogonality, but their direct Coulomb repulsion remains. 3. If one sums all occupied Fock eigenvalues to compare two molecular geometries, what problem occurs? Answer: Electron–electron interactions are double counted and nuclear repulsion may be missing. The proper Hartree–Fock total-energy expression must be used.