Density Functional Theory (DFT)

Why DFT matters

• ~80% of computational chemistry papers use DFT. The dominant electronic-structure method since the late 1990s. B3LYP alone appears in roughly half of computational papers; PBE dominates solid-state work; M06-2X and ωB97X-D are competitive for non-covalent interactions.
• Cost scaling makes 1000-atom systems tractable. Pure GGA functionals scale as O(N³) where N is basis size. Hybrid functionals scale as O(N⁴) because of exact exchange. Compare CCSD(T) at O(N⁷) — capped at ~30 atoms — and DFT at > 1000 atoms is routine.
• Hohenberg-Kohn 1964: a half-page proof that changed everything. Pierre Hohenberg and Walter Kohn published two theorems showing the ground-state energy is determined entirely by the electron density. The proofs are short and rigorous; in principle no information is lost relative to wavefunction methods.
• Kohn-Sham 1965 made it practical. Walter Kohn and Lu Sham introduced a fictitious system of non-interacting electrons that produces the same density as the real one. The resulting Kohn-Sham equations are one-electron Schrödinger-like equations solvable by SCF iteration — the same machinery as Hartree-Fock.
• Materials Genome Initiative runs on DFT. The 2011 US initiative cataloged the electronic structure of essentially every known stoichiometric crystal in databases like Materials Project, OQMD, and AFLOW — tens of thousands of compounds, all computed with PBE and a fixed pseudopotential library.
• Drug discovery uses DFT for binding-pose energetics. Mechanism-of-action studies for kinase inhibitors, oxidative metabolism prediction (cytochrome P450 site selectivity), and conformational energy ranking all rely on DFT — usually B3LYP or ωB97X-D — for thousands of poses.
• The Walter Kohn / John Pople 1998 Nobel. Kohn for DFT itself; Pople for the development of computational methods including the Gaussian program, which made these methods routinely usable. They split the prize evenly.

Common misconceptions

• DFT is exact in principle. The Hohenberg-Kohn theorem is exact in principle. But the universal functional F[ρ] is not known explicitly, so all practical DFT uses approximate exchange-correlation functionals. The functional is the source of essentially all DFT errors.
• Kohn-Sham orbitals are real. No. They are the orbitals of a fictitious non-interacting system constructed to give the same density as the real interacting system. KS orbital energies are not strictly ionization energies (only the HOMO energy maps to −IP for an exact functional, by Janak's theorem), and KS orbital shapes only approximately represent the real wavefunction.
• DFT is variational without caveats. The energy E[ρ] is variational with respect to the exact functional. Approximate functionals can give energies below the true ground state, especially for stretched bonds and transition states.
• Hybrid functionals are always better. Hybrids improve thermochemistry of small organic molecules. They are often worse for transition-metal complexes, where the Hartree-Fock exchange admixture overlocalizes d-electrons. PBE and SCAN often outperform B3LYP for crystals and metals.
• DFT handles dispersion forces. Pure local and semilocal functionals miss London dispersion entirely. Grimme's -D3 (2010) and -D4 (2019) corrections add an empirical pairwise dispersion term −C₆/r⁶ to the energy. Without it B3LYP underbinds the benzene dimer by ~10 kJ/mol.
• Better functional always means better results. Functional rankings depend on the property and the system. M06-2X is excellent for non-covalent interactions; PBE excels for solids; SCAN is balanced; B3LYP is a standard reference for organics. There is no universal best functional.

From wavefunctions to densities and back

The wavefunction ψ(r₁,...,r_N) for an N-electron system depends on 3N spatial coordinates. Even storing it on a grid for benzene (N = 42) is impossible — 1000⁴² points exceeds the number of atoms in the observable universe. The density ρ(r) = N ∫ |ψ|²dr₂...dr_N depends on just 3 coordinates regardless of N. Hohenberg and Kohn's first theorem says the ground-state ρ uniquely determines the external potential V_ext(r) up to a constant, which determines the Hamiltonian, which determines ψ and every other property. Their second theorem says the true ground-state ρ minimizes a universal functional E[ρ] = ∫V_ext(r)ρ(r)dr + F[ρ]. So in principle solving for ρ (3D problem) gives all ground-state information without ever computing ψ (3N-D problem).

The Kohn-Sham construction makes this practical. Introduce a fictitious system of N non-interacting electrons with a Slater-determinant wavefunction built from one-electron orbitals φ_i. The density is ρ(r) = Σ|φ_i(r)|². Choose the effective potential v_eff(r) such that this fictitious density equals the real interacting density. The total energy decomposes into E = T_s[ρ] + ∫V_ext(r)ρ(r)dr + J[ρ] + E_xc[ρ], where T_s is the non-interacting kinetic energy, J is the classical Hartree (Coulomb) repulsion, and E_xc is the exchange-correlation functional that absorbs all the residual quantum effects. The Kohn-Sham one-electron equations resemble Hartree-Fock; they're solved by SCF iteration with cost O(N³) for the dominant matrix operations.

The catch is that the exact E_xc[ρ] is unknown. All practical DFT approximates it. Local Density Approximation (LDA) takes E_xc[ρ] = ∫ε_xc^HEG(ρ(r))ρ(r)dr from the homogeneous electron gas energy density. Generalized Gradient Approximation (GGA) adds dependence on ∇ρ; PBE (Perdew-Burke-Ernzerhof, 1996) is the standard GGA. Hybrid functionals mix a fraction of exact (Hartree-Fock) exchange; B3LYP uses 20% HF exchange combined with Becke's B88 exchange (1988) and the LYP correlation functional. Range-separated hybrids (ωB97X) use long-range HF and short-range DFT; double-hybrids add a fraction of MP2 correlation. Each rung adds cost and improves typical errors but introduces functional-specific failure modes.

DFT functionals — Jacob's ladder, with cost and accuracy

Rung	Functional class	Examples	Cost scaling	Typical thermochemistry error	Best use cases
1	LDA	SVWN, PW92	O(N3)	~70 kJ/mol	Solids, metals (overbinding)
2	GGA	PBE, BLYP, PW91	O(N3)	~30 kJ/mol	Solid-state, materials, fast screening
3	meta-GGA	TPSS, M06-L, SCAN	O(N3)	~20 kJ/mol	Versatile; SCAN balances solids and molecules
4	Hybrid (global)	B3LYP, PBE0	O(N4)	~15-20 kJ/mol	Organic thermochemistry, ground-state
4'	Range-separated hybrid	ωB97X-D, CAM-B3LYP, LC-BLYP	O(N4)	~10-15 kJ/mol	Charge-transfer states, large π-systems
5	Double-hybrid	B2PLYP, DSD-PBEP86	O(N5)	~5-10 kJ/mol	High-accuracy benchmark; expensive

HF vs MP2 vs CCSD(T) vs DFT — when to use which

Method	Cost scaling	Tractable system size	Accuracy on bond energies	Strength	Weakness
Hartree-Fock	O(N4)	~500 atoms	50-100 kJ/mol error	Reference for post-HF; reproducible	No correlation, no dispersion
MP2	O(N5)	~100 atoms	~30 kJ/mol error	Cheap correlation; good for noncovalent	Diverges for stretched bonds
CCSD	O(N6)	~50 atoms	~10 kJ/mol error	Reliable for closed-shell	Misses triple excitations
CCSD(T)	O(N7)	~30 atoms	~4 kJ/mol error	"Gold standard" benchmark	Costly; multireference fails
DFT (PBE)	O(N3)	> 1000 atoms	~30 kJ/mol error	Materials/solids; fast	Underestimates band gaps
DFT (B3LYP)	O(N4)	~500 atoms	~15-30 kJ/mol error	Organic chemistry workhorse	Misses dispersion; -D3 patch

Applications and examples

• ~80% of computational chemistry papers use DFT. Surveys of leading journals (J. Chem. Phys., JCTC, J. Phys. Chem.) consistently show B3LYP, PBE, ωB97X-D, and M06-2X dominating the methods sections. CCSD(T) appears mainly as a benchmark or for reactions where accurate barriers are essential.
• Materials Project and OQMD. Open databases of DFT-computed properties for ~140,000 inorganic crystals (PBE-level, GGA pseudopotentials). They power high-throughput materials screening for batteries, photocatalysts, thermoelectrics, and topological materials.
• Catalyst design. Heterogeneous catalysis on metal surfaces — N₂ dissociation on Fe (Haber), CO₂ reduction on Cu, oxygen reduction on Pt — is computed with PBE plus van der Waals corrections. Activation barriers within 20 kJ/mol of experiment guide alloy and facet choice.
• Drug discovery. Cytochrome-P450 metabolism site prediction, kinase inhibitor binding poses, organic photocatalyst screening (Doyle, Knowles labs) all use B3LYP or M06-2X with -D3. Throughput of 10³-10⁵ compounds per study is routine.
• The 1998 Nobel Prize in Chemistry. Walter Kohn for DFT itself (Hohenberg-Kohn 1964, Kohn-Sham 1965); John Pople for "the development of computational methods in quantum chemistry," principally the Gaussian program. They split the prize evenly. Becke's 1988 functional and the 1993 B3LYP publication had been transformational by then.