Skip to content

Density Functional Theory: Conceptual Overview

Theorem 1: The ground-state electron density n(r)n(\mathbf{r}) uniquely determines the external potential Vext(r)V_{\text{ext}(\mathbf{r})} (up to an additive constant), and hence the full many-body Hamiltonian and all ground-state properties.

Theorem 2: The ground-state energy is a functional of the density: E[n]=FHK[n]+Vext(r)n(r)d3rE[n] = F_{\text{HK}[n] + \int V_{\text{ext}(\mathbf{r})n(\mathbf{r})\,d^3r}} And the variational principle applies: E0E[n]E_0 \leq E[n] for any trial density n(r)n(\mathbf{r}).

The interacting system is mapped to a fictitious system of non-interacting electrons in an effective potential:

\left[-\frac{\hbar^2}{2m}\nabla^2 + V_{\text{eff}[n](\mathbf{r})\right]\psi_i(\mathbf{r}) = \varepsilon_i\psi_i(\mathbf{r})}

n(r)=i=1Nψi(r)2(summing over occupied states)n(\mathbf{r}) = \sum_{i=1}^{N}|\psi_i(\mathbf{r})|^2 \quad \text{(summing over occupied states)}

Veff=Vext+VH[n]+Vxc[n]V_{\text{eff} = V_{\text{ext} + V_H[n] + V_{\text{xc}[n]}}}

VH[n](r)=e2n(r")rrd3r(Hartree potential)V_H[n](\mathbf{r}) = e^2\int\frac{n(\mathbf{r}")}{|\mathbf{r} - \mathbf{r}'|}\,d^3r' \quad \text{(Hartree potential)}

The exchange-correlation functional Vxc[n]V_{\text{xc}[n]} contains all many-body effects beyond the classical Hartree approximation.

The Hartree potential includes the interaction of each electron with itself. This self-interaction error is not cancelled by the local density approximation (LDA) for VxcV_{\text{xc}}. Consequences:

  • Wrong asymptotic behaviour: Veff(r)e2/rV_{\text{eff}(r \to \infty) \to -e^2/r} (correct) vs. Veff0V_{\text{eff} \to 0} (LDA, wrong)
  • Underestimation of band gaps by 30—50%
  • Incorrect description of charge transfer excitations

Hybrid functionals (e.g., B3LYP, HSE06) and range-separated functionals partially correct this.

Worked Example 14.1: Thomas--Fermi Theory

The simplest density functional theory: the Thomas—Fermi model treats the kinetic energy as a local functional of the density:

TTF[n]=3210m(3π2)2/3n5/3(r)d3r=CTFn5/3d3rT_{\text{TF}[n] = \frac{3\hbar^2}{10m}(3\pi^2)^{2/3}\int n^{5/3}(\mathbf{r})\,d^3r = C_{\text{TF}\int n^{5/3}\,d^3r}}

For an atom with nuclear charge ZeZeMinimising E[n]=TTF[n]Ze2n(r)/rd3r+12e2n(r)n(r)/rrd3rd3rE[n] = T_{\text{TF}[n] - Ze^2\int n(\mathbf{r})/r\,d^3r + \frac{1}{2}e^2\iint n(\mathbf{r})n(\mathbf{r}')/|\mathbf{r}-\mathbf{r}'|\,d^3rd^3r'}:

The variational equation gives:

CTFn2/3=Ze2re2n(r)rrd3rC_{\text{TF}\,n^{2/3} = \frac{Ze^2}{r} - e^2\int\frac{n(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\,d^3r'}

This integral equation can be solved by scaling: n(r)=(Z/a03)g(r/a0Z1/3)n(r) = (Z/a_0^3)\,g(r/a_0 Z^{-1/3}) where gg is a universal function.

The Thomas—Fermi energy: ETF=37(9π/2)2/3Z7/3e22a0=20.8Z7/3E_{\text{TF} = -\frac{3}{7}(9\pi/2)^{2/3}\frac{Z^{7/3}e^2}{2a_0} = -20.8\,Z^{7/3}} eV.

This gives reasonable total energies for heavy atoms but fails qualitatively for light atoms (no shell structure, no chemical bonding).

flowchart TD
A[15_Density Functional Theory Conceptual Overview] --> B[Key Concepts]
A --> C[Core Principles]
A --> D[Practical Applications]
B --> E[Fundamental definitions]
C --> F[Design patterns]
D --> G[Real-world usage]

Problem. Find the energy levels and normalised wave functions for a particle in a 1D infinite square well of width aa.

Solution. ψn(x)=2asin(nπxa)\psi_n(x) = \sqrt{\frac{2}{a}}\sin\left(\frac{n\pi x}{a}\right), En=n2π222ma2E_n = \frac{n^2\pi^2\hbar^2}{2ma^2}, n=1,2,3,n = 1, 2, 3, \ldots

\blacksquare

Problem. An electron is confined to a region of width 0.1nm0.1 \mathrm{ nm}. Estimate the minimum uncertainty in its momentum.

Solution. Δx=0.1×109m\Delta x = 0.1 \times 10^{-9} \mathrm{ m}. Δp2Δx=1.055×10342×1010=5.28×1025kgm/s{\Delta p \geq \frac{\hbar}{2\Delta x} = \frac{1.055 \times 10^{-34}}{2 \times 10^{-10}} = 5.28 \times 10^{-25} \mathrm{ kg\cdot} m/s}.

\blacksquare

  • Confusing the wave function and probability density. ψ(x)2|\psi(x)|^2 is the probability density; ψ\psi itself is complex and not directly observable. Fix: P(x)dx=ψ(x)2dxP(x)\, dx = |\psi(x)|^2\, dx is the probability of finding the particle in [x,x+dx][x, x + dx].
  • Wrong commutator interpretation. If [A^,B^]=0[\hat{A}, \hat{B}] = 0, the observables share eigenstates and can be simultaneously measured. If not, they obey the uncertainty principle. Fix: [x^,p^]=i[\hat{x}, \hat{p}] = i\hbar implies ΔxΔp/2\Delta x \cdot \Delta p \geq \hbar/2.
  • Confusing time-dependent and time-independent Schrödinger equations. Time-dependent: iψt=H^ψi\hbar\frac{\partial\psi}{\partial t} = \hat{H}\psi. Time-independent: H^ϕ=Eϕ\hat{H}\phi = E\phi. Fix: Time-independent gives stationary states (energy eigenvalues); time-dependent describes evolution.
  • Postulates: state vector ψ|\psi\rangle, observables as Hermitian operators, measurement gives eigenvalues.
  • Schrödinger equation: iψ/t=H^ψi\hbar\partial\psi/\partial t = \hat{H}\psi; stationary states: H^ϕn=Enϕn\hat{H}\phi_n = E_n\phi_n.
  • Commutators and uncertainty: [A^,B^]0ΔAΔB12[A^,B^][\hat{A}, \hat{B}] \neq 0 \Rightarrow \Delta A \cdot \Delta B \geq \frac{1}{2}|\langle[\hat{A}, \hat{B}]\rangle|.
  • Key systems: infinite square well, harmonic oscillator, hydrogen atom.
TopicSiteLink
[Quantum Physics]A-LevelView
[Quantum Physics]IBView
[Quantum Physics]UniversityView

Density functional theory replaces the complicated many-body wavefunction with the much simpler electron density, which depends on only three spatial coordinates instead of three times the number of electrons. The Hohenberg-Kohn theorem guarantees that the ground state energy is a functional of the density, and the Kohn-Sham equations reformulate the problem as independent particles moving in an effective potential. This makes electronic structure calculations tractable for systems with hundreds of atoms. The exchange-correlation functional captures all the complicated many-body effects, and improving it is the central challenge of the field.

  • Hohenberg—Kohn mapping. The external potential Vext(r)V_{\text{ext}}(\mathbf{r}) is a unique functional of the ground-state density n(r)n(\mathbf{r}) (up to a constant), so all ground-state properties are functionals of nn alone.
  • Kohn—Sham self-consistency. The effective potential Veff[n]V_{\text{eff}}[n] depends on n(r)n(\mathbf{r}), which is built from the Kohn—Sham orbitals that are themselves solutions of the equations with Veff[n]V_{\text{eff}}[n]. This requires iterative self-consistent solution.

This section provides detailed coverage of advanced concepts, including full derivations, proofs, and extended examples.

Complete mathematical derivations and proofs are provided where appropriate. Each step is explained to ensure understanding of the underlying reasoning.

Advanced examples demonstrate the application of concepts to complex problems. These examples go beyond standard exam questions to develop deeper understanding.

This material connects to current research and advanced applications in the field. Understanding these connections provides context for the study material.

Ensure you have mastered the prerequisite material before attempting this advanced content.

This section provides detailed coverage of advanced concepts, including full derivations, proofs, and extended examples.

Complete mathematical derivations and proofs are provided where appropriate. Each step is explained to ensure understanding of the underlying reasoning.

Advanced examples demonstrate the application of concepts to complex problems. These examples go beyond standard exam questions to develop deeper understanding.

This material connects to current research and advanced applications in the field. Understanding these connections provides context for the study material.

Ensure you have mastered the prerequisite material before attempting this advanced content.