The Kohn-Sham Equations and Exchange-Correlation Functionals in DFT

Density Functional Theory (DFT) is a powerful quantum mechanical method used to investigate the electronic structure of materials. At its heart lies the Kohn-Sham (KS) scheme, which transforms the intractable many-body problem into a solvable problem of non-interacting electrons moving in an effective potential. This section delves into the core of this transformation: the Kohn-Sham equations and the crucial role of exchange-correlation functionals.

The Kohn-Sham Ansatz

The fundamental theorem of DFT states that the ground-state properties of a system are uniquely determined by its ground-state electron density, $\rho(r)$ . However, directly solving the Schrödinger equation for a many-electron system is computationally prohibitive. The Kohn-Sham ansatz provides a clever workaround by introducing a fictitious system of non-interacting electrons that has the same ground-state density as the real, interacting system. This allows us to use single-particle equations, similar to the Hartree-Fock method, but with an effective potential that accounts for all the complex many-body interactions.

The Kohn-Sham equations map a complex interacting system to a simpler non-interacting one with the same electron density.

The Kohn-Sham scheme introduces a fictitious system of non-interacting electrons that mimics the density of the real, interacting system. This allows us to use single-particle equations.

The Kohn-Sham (KS) scheme proposes that there exists a set of non-interacting electrons moving in an effective potential, $v_{KS}(r)$ , such that their ground-state density, $\rho_{KS}(r)$ , is identical to the ground-state density of the real, interacting system, $\rho_{real}(r)$ . The kinetic energy of this non-interacting system, $T_s[\rho]$ , is then used as an approximation for the kinetic energy of the real system. The total energy functional is written as $E[\rho] = T_s[\rho] + E_{ext}[\rho] + E_H[\rho] + E_{xc}[\rho]$ , where $E_{ext}$ is the external potential energy, $E_H$ is the Hartree (classical electrostatic repulsion) energy, and $E_{xc}$ is the exchange-correlation energy.

The Kohn-Sham Equations

By applying the variational principle to the Kohn-Sham energy functional, we arrive at a set of single-particle Schrödinger-like equations for the fictitious non-interacting electrons. These are the Kohn-Sham equations:

The Kohn-Sham equations are: $\left(-\frac{\hbar^2}{2m}\nabla^2 + v_{KS}(r)\right)\phi_i(r) = \epsilon_i\phi_i(r)$ . Here, $\phi_i(r)$ are the Kohn-Sham orbitals, and $\epsilon_i$ are their corresponding orbital energies. The effective potential $v_{KS}(r)$ is given by $v_{KS}(r) = v_{ext}(r) + v_H(r) + v_{xc}(r)$ . The term $v_{ext}(r)$ is the external potential (e.g., from atomic nuclei), $v_H(r)$ is the Hartree potential (classical electron-electron repulsion), and $v_{xc}(r)$ is the exchange-correlation potential, which is the functional derivative of the exchange-correlation energy with respect to the density: $v_{xc}(r) = \frac{\delta E_{xc}[\rho]}{\delta \rho(r)}$ . The electron density is then constructed from the occupied KS orbitals: $\rho(r) = \sum_{i=1}^{N} |\phi_i(r)|^2$ , where N is the number of electrons.

📚

Text-based content

Library pages focus on text content

Exchange-Correlation Functionals: The Heart of DFT Accuracy

The accuracy of DFT calculations hinges entirely on the approximation used for the exchange-correlation energy functional, $E_{xc}[\rho]$ . This term encapsulates all the complex quantum mechanical effects of electron exchange and correlation that are not captured by the kinetic energy of the non-interacting system or the Hartree term. Since the exact form of $E_{xc}[\rho]$ is unknown for systems with more than one electron, various approximations have been developed, forming a hierarchy often referred to as 'Jacob's Ladder'.

Functional Type	Key Idea	Examples
Local Density Approximation (LDA)	Assumes the exchange-correlation energy at a point depends only on the density at that point, as if it were a homogeneous electron gas.	SVWN, PZ81
Generalized Gradient Approximation (GGA)	Includes the density gradient in addition to the density itself, accounting for inhomogeneity.	PBE, BLYP, PW91
Meta-GGA	Includes the kinetic energy density or the Laplacian of the density.	TPSS, SCAN
Hybrid Functionals	Mixes a fraction of exact (Hartree-Fock) exchange with GGA or meta-GGA exchange and correlation.	B3LYP, PBE0, HSE06

The choice of exchange-correlation functional is critical. LDA is simple but often underestimates bond lengths and overestimates binding energies. GGAs generally improve upon LDA, while meta-GGAs and hybrids offer higher accuracy but at increased computational cost.

Self-Consistency Loop

Solving the Kohn-Sham equations is an iterative process. Since the effective potential $v_{KS}(r)$ depends on the electron density $\rho(r)$ , which in turn depends on the KS orbitals $\phi_i(r)$ , the equations must be solved self-consistently. This involves starting with an initial guess for the density, calculating $v_{KS}(r)$ , solving the KS equations to obtain new orbitals, constructing a new density from these orbitals, and repeating until the density (and energy) converges.

Loading diagram...

Practical Considerations

In practice, DFT calculations are performed using specialized software packages. Understanding the underlying KS equations and the nature of exchange-correlation functionals is crucial for interpreting results, choosing appropriate approximations, and troubleshooting calculations. The trade-off between accuracy and computational cost is a constant consideration in selecting an exchange-correlation functional for a given problem.

Learning Resources

Density Functional Theory (DFT) - Wikipedia(wikipedia)

Provides a broad overview of DFT, its history, fundamental theorems, and common approximations, including Kohn-Sham theory.

Introduction to Density Functional Theory - Materials Science(tutorial)

An interactive tutorial covering the basics of DFT, including the Kohn-Sham equations and the role of functionals, with visual aids.

Kohn-Sham Equations - Quantum ESPRESSO Documentation(documentation)

Details the Kohn-Sham equations as implemented in the Quantum ESPRESSO software package, offering practical context.

Exchange-Correlation Functionals in DFT - A Review(paper)

A comprehensive review article discussing the various types of exchange-correlation functionals and their performance in DFT calculations.

Understanding DFT: The Kohn-Sham Equations(video)

A video explanation of the Kohn-Sham equations, breaking down the concepts for easier comprehension.

The 'Jacob's Ladder' of Density Functionals(paper)

An article that categorizes and explains the hierarchy of DFT exchange-correlation functionals, often referred to as 'Jacob's Ladder'.

Introduction to Density Functional Theory - Lecture Notes(documentation)

Lecture notes providing a solid theoretical foundation for DFT, including detailed derivations of the Kohn-Sham equations.

Practical DFT Calculations - A Guide for Chemists(blog)

A blog post offering practical advice on performing DFT calculations, including selecting appropriate functionals and interpreting results.

The Exchange-Correlation Problem in DFT(video)

A video focusing specifically on the challenges and approximations related to the exchange-correlation functional in DFT.

Computational Materials Science with DFT - Overview(documentation)

A presentation slide deck offering a concise overview of DFT basics, including the Kohn-Sham approach, relevant for computational materials science.