1 Background and theoretical context

The Kohn–Sham equations arise from the problem of describing interacting electrons in atoms, molecules, and solids. In a fully quantum treatment, the electrons are correlated through their mutual repulsion and the antisymmetry of the many-body wavefunction. Direct solution becomes impractical for systems with more than a few electrons, which led to the search for reformulations that retain essential physical information while reducing computational complexity.

1.1 Many-electron Schrödinger equation

The nonrelativistic many-electron Schrödinger equation describes a system in terms of a wavefunction that depends on the coordinates of all electrons simultaneously. This wavefunction contains complete information about the state, but its dimensionality grows rapidly with particle number. As a result, the exact equation is analytically solvable only in very limited cases.

1.2 Challenges of electron correlation

Electron correlation refers to the coupled motion of electrons caused mainly by Coulomb repulsion and by the antisymmetry requirement of fermionic wavefunctions. Mean-field pictures that treat electrons as moving independently often miss important correlation effects, such as bond energies, magnetic behavior, and dispersion interactions. Capturing these effects accurately is one of the central difficulties of electronic-structure theory.

1.3 Density functional theory

Density functional theory reformulates the many-electron problem in terms of the electron density rather than the full many-body wavefunction. Since the density depends only on three spatial variables, it offers a major simplification. In principle, the exact ground-state properties of a system can be derived from the density alone.

1.3.1 Hohenberg–Kohn theorems

The Hohenberg–Kohn theorems establish that the ground-state electron density uniquely determines the external potential, up to a constant, and that the exact ground-state energy is obtained by minimizing an energy functional of the density. These results provide the conceptual basis for density functional theory and justify the use of the density as the fundamental variable.

1.3.2 Ground-state density as the basic variable

By making the density the key quantity, density functional theory avoids direct handling of the full many-electron wavefunction. This shift makes it possible to construct practical approximations that are often accurate enough for chemistry and materials science, especially for ground-state structures and energies.

1.4 Motivation for the Kohn–Sham approach

Although the exact density functional exists in principle, its explicit form is unknown. The Kohn–Sham approach introduces an auxiliary system of noninteracting electrons designed to reproduce the exact ground-state density. This strategy separates the difficult many-body effects into a remainder term, the exchange-correlation functional, while preserving a tractable single-particle framework.

2 Formulation of the Kohn–Sham equations

The Kohn–Sham equations describe electrons moving in an effective potential chosen so that the resulting density matches that of the interacting system. Each electron is represented by a Kohn–Sham orbital, and the orbitals are solved self-consistently because the effective potential depends on the density built from those same orbitals.

2.1 Auxiliary noninteracting electron system

The auxiliary system contains noninteracting electrons occupying orbitals whose combined density reproduces the exact ground-state density of the interacting system. This does not imply that the physical electrons cease to interact; rather, their interaction effects are encoded in the effective potential.

2.2 Effective single-particle potential

The effective potential is the sum of several contributions, each reflecting a different physical influence on the electrons. Together, these terms define a one-electron problem whose solution approximates or, in the exact limit, reproduces the true density.

2.2.1 External potential

The external potential represents the interaction of electrons with nuclei and any applied fields. In molecular and solid-state systems, it is typically determined by the arrangement of atomic centers and their charges.

2.2.2 Hartree potential

The Hartree potential is the classical electrostatic potential generated by the electron density itself. It accounts for the average Coulomb repulsion among electrons, but it does not include exchange or correlation effects.

2.2.3 Exchange-correlation potential

The exchange-correlation potential contains the nonclassical effects of antisymmetry and many-body correlation. It is the least known part of the theory and must usually be approximated. Its quality strongly influences the accuracy of Kohn–Sham calculations.

2.3 Kohn–Sham orbitals

Kohn–Sham orbitals are mathematical functions whose occupations determine the electron density. They are not generally interpreted as exact physical orbitals, but they provide a useful and often chemically intuitive representation of the electronic structure.

2.4 Self-consistent field equations

Because the effective potential depends on the density, and the density depends on the orbitals, the equations must be solved iteratively. An initial guess is refined through repeated updates until the input and output densities are consistent within a chosen tolerance.

3 Mathematical structure

The Kohn–Sham formalism is built around a set of coupled eigenvalue equations and an associated total-energy functional. The structure resembles a single-particle quantum problem, but the dependence of the potential on the density makes the equations nonlinear in practice.

3.1 Eigenvalue problem

Each Kohn–Sham orbital satisfies a differential equation with an effective one-electron Hamiltonian. The resulting eigenvalues are often used as indicators of electronic levels, although they are not always equal to observable excitation energies.

3.2 Electron density expression

The electron density is obtained by summing the squared magnitudes of the occupied orbitals, weighted by their occupations. For closed-shell systems, this yields a compact expression that links the many-electron state to the auxiliary single-particle picture.

3.3 Total energy functional

The total energy in Kohn–Sham theory is expressed as a functional of the density and, equivalently, of the orbitals. It includes kinetic, electrostatic, and exchange-correlation contributions, along with interaction terms involving the external potential.

3.3.1 Kinetic-energy term

The kinetic energy of the auxiliary noninteracting system is computed directly from the Kohn–Sham orbitals. This term is treated exactly within the Kohn–Sham construction, which is one reason the method is often more accurate than purely density-based approximations.

3.3.2 Coulomb interaction term

The Coulomb term accounts for the electrostatic repulsion between electrons through the Hartree energy. It is straightforward to evaluate from the density, though it must be corrected by exchange-correlation effects to avoid double counting of interaction contributions.

3.3.3 Exchange-correlation energy

The exchange-correlation energy collects the difference between the true interacting system and the auxiliary noninteracting model. It includes exchange effects from antisymmetry and correlation effects from electron-electron motion beyond the mean field.

3.4 Functional derivatives

The exchange-correlation potential is obtained as the functional derivative of the exchange-correlation energy with respect to the density. This derivative enters the Kohn–Sham equations and determines how the effective potential responds to changes in the electron distribution.

4 Exchange-correlation approximations

Because the exact exchange-correlation functional is unknown, practical calculations rely on approximations. Different forms balance accuracy, efficiency, and general applicability, and the choice of functional often determines the reliability of a calculation.

4.1 Local density approximation

The local density approximation assumes that the exchange-correlation energy at each point depends only on the density at that point, using results from the uniform electron gas as a reference. It is simple and often effective for systems with slowly varying densities.

4.2 Generalized gradient approximation

Generalized gradient approximations improve on local models by including density gradients. This additional information often yields better descriptions of molecular geometries, bond strengths, and structural trends in solids.

4.3 Meta-GGA functionals

Meta-GGA functionals incorporate further ingredients such as the kinetic-energy density or higher-order density information. They can improve accuracy in many cases while retaining a semilocal character.

4.4 Hybrid functionals

Hybrid functionals mix exchange from Kohn–Sham theory with a portion of exact exchange from Hartree–Fock theory. This often improves the description of band structures, reaction energies, and molecular properties, though at greater computational cost.

4.5 Beyond-standard approximations

More advanced approaches include double hybrids, range-separated methods, and functionals designed for specific classes of systems. These methods seek higher fidelity for challenging cases where standard approximations are insufficient.

5 Solution methods

Solving the Kohn–Sham equations in practice requires numerical algorithms, basis representations, and convergence control. The implementation is usually tailored to the geometry and periodicity of the system under study.

5.1 Self-consistent iteration

A typical calculation begins with an initial density or set of orbitals, then repeatedly constructs the effective potential and solves the Kohn–Sham equations until convergence is reached. Mixing schemes are often used to stabilize the iteration.

5.2 Basis-set choices

The orbitals must be represented in a finite basis or on a numerical grid. The choice affects accuracy, computational scaling, and ease of implementation.

5.2.1 Plane-wave basis

Plane waves are widely used for periodic systems because they naturally match crystal symmetry and provide systematic convergence. They are especially common in solid-state and surface calculations.

5.2.2 Localized atomic orbitals

Localized atomic orbitals are compact and efficient for molecules and large systems with localized bonding. They can reduce computational effort, though basis-set completeness must be carefully managed.

5.2.3 Real-space grids

Real-space grids discretize the equations directly in position space. This approach offers flexibility and can be convenient for complex geometries, finite systems, and certain numerical methods.

5.3 Convergence criteria

Convergence is typically assessed using changes in total energy, density, orbital occupations, or residual norms. Strict criteria are important when comparing energies between similar structures or when subtle properties are of interest.

5.4 Numerical implementation

Efficient implementations rely on linear algebra, fast evaluation of potentials, and strategies for handling large matrices or grids. Modern codes often include parallel computing techniques to treat systems with many atoms or large basis sets.

6 Applications

Kohn–Sham density functional theory is used across physics, chemistry, and materials science. Its combination of manageable cost and useful accuracy makes it a standard tool for studying electronic structure.

6.1 Molecular electronic structure

In molecular systems, Kohn–Sham calculations are used to predict geometries, vibrational trends, charge distributions, and relative stabilities. They also help interpret bonding patterns and compare isomer energies.

6.2 Solid-state physics

For crystalline materials, the method is central to calculating electronic bands, densities of states, and equilibrium structures. It provides a practical framework for studying semiconductors, insulators, and metals.

6.3 Surface and interface studies

Surfaces and interfaces often have electronic properties that differ from the bulk because of reduced coordination and local chemical environments. Kohn–Sham methods are widely used to analyze adsorption, reconstruction, and work-function changes.

6.4 Materials design and prediction

The method supports computational screening of new compounds and estimation of structural, electronic, and mechanical properties. It is frequently used in the early stages of materials discovery to narrow experimental searches.

6.5 Chemical reactivity analysis

Kohn–Sham orbitals and densities are commonly used to discuss reactivity, charge transfer, and frontier orbital effects. Although such interpretations are approximate, they offer practical insight into reaction pathways and intermolecular interactions.

7 Limitations and interpretive issues

Despite its wide success, Kohn–Sham density functional theory has limitations rooted mainly in approximate exchange-correlation functionals and in the interpretation of orbital energies. These issues can affect both quantitative predictions and physical interpretation.

7.1 Approximate nature of exchange-correlation functionals

Since the exact exchange-correlation functional is not known, all practical calculations involve approximations. Errors may appear in bond energies, equilibrium structures, magnetic properties, and weak interactions, depending on the functional used.

7.2 Band-gap problem

Standard Kohn–Sham calculations often underestimate band gaps in solids. This reflects limitations of approximate functionals and the distinction between Kohn–Sham eigenvalues and true excitation energies.

7.3 Self-interaction error

In many approximate functionals, an electron spuriously interacts with itself through incomplete cancellation of Coulomb terms. This self-interaction error can distort charge localization, reaction barriers, and dissociation behavior.

7.4 Degeneracy and fractional occupation

Systems with nearly degenerate states may require fractional occupations or ensemble descriptions. Proper handling of such cases is important for metallic systems, open shells, and cases with symmetry breaking.

7.5 Relationship to true quasiparticle energies

Kohn–Sham eigenvalues should not be treated automatically as experimental ionization energies or electron affinities. While some orbital energies have useful connections to physical observables, the relationship is indirect and often approximate.

Several extensions broaden the Kohn–Sham framework to include spin, time dependence, finite temperature, and orbital-dependent effects. These generalizations retain the density-based philosophy while addressing additional physical phenomena.

8.1 Spin-polarized Kohn–Sham equations

Spin-polarized formulations treat spin-up and spin-down densities separately. They are essential for magnetic materials, open-shell molecules, and systems where spin asymmetry matters.

8.2 Time-dependent density functional theory

Time-dependent density functional theory extends the density approach to time-dependent processes. It is used for electronic excitations, response properties, and the dynamics of systems under external perturbations.

8.3 Finite-temperature formulations

Finite-temperature density functional theory introduces thermal occupations and entropy contributions. It is useful for warm dense matter, metallic systems, and calculations in which electronic temperature is relevant.

8.4 Orbital-dependent density functionals

Orbital-dependent functionals depend explicitly on the Kohn–Sham orbitals rather than only on the density. They can improve accuracy in some settings, but they usually require more elaborate algorithms.

8.5 Generalized Kohn–Sham schemes

Generalized Kohn–Sham methods allow a broader class of effective operators, including nonlocal terms. They provide a flexible framework for incorporating more sophisticated approximations while retaining a single-particle structure.

9 Historical development and impact

The Kohn–Sham equations transformed density functional theory from a formal concept into a practical computational framework. Their influence has spread across multiple scientific disciplines and continues to shape electronic-structure methods.

9.1 Original 1965 formulation

Walter Kohn and Lu Jeu Sham introduced the method in 1965 as a way to construct a practical implementation of density functional theory. Their formulation connected the exact density principle to a solvable auxiliary system, making applications far more feasible.

9.2 Influence on computational chemistry

In computational chemistry, the Kohn–Sham approach became a standard tool for predicting molecular structures and energetics. Its relatively favorable cost-to-accuracy ratio made it attractive for systems too large for high-level wavefunction methods.

9.3 Role in condensed-matter theory

In condensed-matter physics, the formalism became central to first-principles modeling of electronic structure in solids. It supports the study of bonding, defects, phonon-related properties, and electronic response in complex materials.

9.4 Nobel Prize recognition

The foundational role of density functional theory and the Kohn–Sham framework was recognized in the Nobel Prize in Chemistry awarded to Walter Kohn in 1998. This distinction reflected the broad scientific impact of the approach on theoretical and computational science.