Introduction to Density Functional Theory (DFT)

Density Functional Theory (DFT) is based on the principle that all ground-state properties of an electronic system are fully determined by its electron density $ \rho(r) $. In other words, the behavior of a system of electrons and atomic nuclei can be described entirely through this scalar function of position.

Hohenberg–Kohn Theorems

The foundation of DFT lies in two theorems established by Hohenberg and Kohn:

Theorem 1:

The external potential $ V_{\text{ext}}(r) $ — and consequently, the total energy — is a unique functional of the electron density $ \rho(r) $. The total energy functional can be written as:

$$ E[\rho(r)] = F[\rho(r)] + \int V_{\text{ext}}(r) \rho(r) \, d^3r $$

$ E[\rho(r)] $: total energy as a functional of the density
$ F[\rho(r)] $: a universal functional including kinetic and electron–electron interaction energies
$ V_{\text{ext}}(r) $: external potential due to nuclei

Theorem 2:

The ground-state electron density $ \rho_0(r) $ that minimizes $ E[\rho(r)] $ is the exact ground-state density. Thus, all properties of the system can be found by variationally minimizing the energy functional with respect to $ \rho(r) $.

The Kohn–Sham Equations

In practice, DFT is implemented through the Kohn–Sham (KS) equations, which map the original system of strongly interacting electrons onto a fictitious system of non-interacting electrons moving in an effective one-particle potential:

$$ \left( -\frac{\hbar^2}{2m} \nabla^2 + v_{\text{eff}}(r) \right) \varphi_i(r) = \epsilon_i \varphi_i(r) $$

The effective potential $ v_{\text{eff}}(r) $ includes three components:

External potential due to nuclei: $ V_{\text{ext}}(r) $
Coulomb (Hartree) potential from electron–electron repulsion:

$$ v_H(r) = e^2 \int \frac{\rho(r')}{|r - r'|} \, d^3r' $$

Exchange–correlation potential:

$$ v_{xc}(r) = \frac{\delta E_{xc}[\rho(r)]}{\delta \rho(r)} $$

The exchange–correlation term $ E_{xc}[\rho] $ encapsulates all many-body effects beyond the Hartree term and is not known exactly. Much of DFT research focuses on finding accurate approximations to this functional. One widely used approximation is the B3LYP functional.

By solving the Kohn–Sham equations self-consistently, one obtains the ground-state electron density and total energy of the system, with significantly reduced computational cost compared to wavefunction-based methods.

DFT Caveats & Common Errors

Error modes: self‑interaction, delocalization, static correlation, and the derivative discontinuity (band gaps, charge transfer).
Janak & piecewise linearity: only with the exact functional does $-\varepsilon_\text{HOMO}$ equal the first ionization potential.
Dispersion: include D3/D4/VV10 (or use functionals with built‑in dispersion).

Ensemble Interpretation of Density Functional Theory

Although DFT is often presented as a ground-state theory, it can also be understood through the lens of ensemble theory. In fact, the original Hohenberg–Kohn theorems allow for ensemble densities in the presence of degeneracy or fractional occupation. This connects DFT to statistical mechanics by recognizing that the electron density $ \rho(r) $ represents a statistical average over quantum microstates.

More generally, in thermal DFT or finite-temperature extensions, the system is explicitly treated as an ensemble of quantum states with weights given by the Fermi–Dirac distribution. The energy functional then incorporates entropy contributions, and the minimizing density corresponds to a thermal ensemble average.

From the viewpoint of microstate engineering, this highlights that DFT solutions describe not just a single configuration, but a weighted average over accessible electronic structures — and by controlling potential energy surfaces or coupling to environments, one modulates this ensemble in a physically meaningful way.