| Thomas–Fermi model | |
|---|---|
| Name | Thomas–Fermi model |
| Caption | Schematic of electronic density in an atom as approximated by a statistical model |
| Field | Quantum physics |
| Introduced | 1927 |
| Authors | Llewellyn Thomas; Enrico Fermi |
| Notable authors | Paul Dirac |
| Related | Density functional theory; Hartree–Fock method; Kohn–Sham equations |
Thomas–Fermi model
The Thomas–Fermi model is a statistical approximation for the distribution of electrons in atoms, molecules, and solids that replaces the many-electron wavefunction by an electron density. Developed independently by Llewellyn Thomas and Enrico Fermi in 1927, it was among the first density-based approaches to electronic structure and laid conceptual groundwork for modern Density functional theory and computational methods in Quantum physics.
The model originated amid early efforts to make quantum mechanics applicable to multi-electron systems without solving the full Schrödinger equation for each electron. Llewellyn Thomas and Enrico Fermi proposed treating electrons as a degenerate Fermi gas subject to the electrostatic potential of nuclei and the mean field of other electrons. Its statistical nature linked it to earlier work on the Fermi–Dirac statistics and to contemporaneous methods such as the Hartree and Hartree–Fock method. The model influenced later advances by Paul Dirac (exchange corrections) and motivated the rigorous development of Density functional theory by Walter Kohn and others.
The core assumption is that, locally, electrons can be approximated as a homogeneous Fermi gas characterized by a local chemical potential determined by the external (nuclear) potential and the self-consistent electrostatic potential. This semiclassical viewpoint connects to the Fermi energy and to phase-space counting methods used in statistical mechanics. The derivation invokes the Pauli exclusion principle through occupation of single-particle states up to the local Fermi momentum and neglects explicit antisymmetrized multi-electron wavefunctions beyond mean-field exchange approximations. The approach is consistent with the thermodynamic limit and with semiclassical expansions such as the WKB approximation and the Thomas–Fermi screening concept in plasmas and metals.
The Thomas–Fermi model expresses the electron number density n(r) in terms of the local potential φ(r) via a power-law relation derived from the density of states of a free electron gas: n(r) = (1/3π^2)[2m(μ - eφ(r))/ħ^2]^{3/2} for regions where μ > eφ(r), where μ is the chemical potential. The potential satisfies Poisson's equation ∇^2φ(r) = −4πe[n(r) − ρ_nuc(r)] with nuclear charge density ρ_nuc(r). For atoms, this reduces to the dimensionless Thomas–Fermi differential equation for the screening function χ(x): d^2χ/dx^2 = χ^{3/2}/√x with appropriate boundary conditions. The model yields scaling laws such as total energy ∝ Z^{7/3} for atomic number Z, matching asymptotic results from semiclassical analysis and rigorous bounds by mathematical physicists.
Corrections to the original model address exchange and inhomogeneity. Paul Dirac introduced an exchange correction (Thomas–Fermi–Dirac, TFD) adding a local exchange energy density proportional to n(r)^{4/3}. Gradient corrections, such as the Weizsäcker correction, incorporate dependence on ∇n to partially recover shell structure and kinetic energy of inhomogeneous densities. Higher-order semiclassical expansions (e.g., the Kirzhnits expansion) and generalized gradient approximations link the Thomas–Fermi picture to modern exchange-correlation functional construction in Density functional theory. Relativistic extensions (Thomas–Fermi–Dirac–Weizsäcker and relativistic Thomas–Fermi models) address heavy elements and high-density regimes encountered in astrophysical contexts like white dwarf modeling.
Despite its simplicity, the model provides useful scaling insights for atomic energies, radii, and ionization trends across the periodic table, and serves as a baseline for more refined calculations. In solid-state physics it underpins Thomas–Fermi screening theory of metals and doped semiconductors, influencing understanding of screening lengths and plasmon behavior. The model has been applied in plasma physics, condensed matter, and astrophysics for approximate equations of state in high-density environments, and as an initial guess for self-consistent field calculations in computational chemistry and materials modeling codes (e.g., as starting densities for Hartree–Fock or Kohn–Sham equations solvers).
The Thomas–Fermi model neglects shell structure, discrete orbital effects, and detailed exchange–correlation beyond local approximations, making it inaccurate for small Z atoms, molecules, and chemical bonding where orbital-specific features matter. It is semiclassical and fails to capture quantum oscillations described by the Lindhard function or Friedel oscillations at atomic scales without further corrections. Its validity improves in the high-Z or high-density limit, where smooth, averaged properties dominate and asymptotic expansions become accurate. Rigorous mathematical results (e.g., by Elliott Lieb and collaborators) have established bounds and the precise limit in which Thomas–Fermi theory becomes exact.
Numerical solution typically proceeds by discretizing the radial Thomas–Fermi equation for atoms or by iterative self-consistent field methods coupling Poisson's equation with the density relation on spatial grids. Finite-difference, finite-element, and spectral methods are used; acceleration techniques include Newton–Raphson, multigrid solvers, and mixing schemes familiar from electronic structure codes. Modern implementations often appear inside larger packages for density functional theory as initializers or low-cost models for large systems, and in plasma/astrophysical EOS codes where computational efficiency is essential. Benchmarks against Hartree–Fock method and Kohn–Sham calculations quantify systematic errors and guide hybrid approaches that combine Thomas–Fermi ideas with orbital-based corrections.
Category:Quantum mechanics Category:Density functional theory