Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Hohenberg–Kohn theorems

To illustrate the basic idea of DFT, let us consider an isolated molecular system with NN electrons. The anti-symmetric NN-electron wave function describing the electronic ground state of this system, Ψ0(r1,r2,..,rN)\Psi_0(\mathbf{r}_1, \mathbf{r}_2,.., \mathbf{r}_N) , depends on the spatial coordinates of the electrons, {ri}\{\mathbf{r}_i\}, and it is a solution to the time-independent Schrödinger equation

H^Ψ0=E0Ψ0\hat{H} \Psi_0 = E_0 \Psi_0
(1)

where E0E_0 is the associated electronic energy and H^\hat{H} is the electronic Hamiltonian that is composed of the kinetic energy, electron–electron repulsion, and external potential operators

H^=T^+U^+V^\hat{H} = \hat{T} + \hat{U} + \hat{V}

with

T^=−ℏ22me∑i∇i2;U^=∑i<je24πε0∣ri−rj∣;V^=∑iv^(i)\hat{T} = - \frac{\hbar^2}{2 m_\mathrm{e}} \sum_i \nabla_i^2 ; \qquad % \hat{U} = \sum_{i<j}\frac{e^2}{4\pi\varepsilon_0 |\mathbf{r}_i -\mathbf{r}_j|} ; \qquad % \hat{V} = \sum_i \hat{v}(i)
(2)

For an isolated system, the latter reduces to the electron–nuclear attraction operator with a coordinate representation

v(r)=−∑KZKe24πε0∣r−RK∣v(\mathbf{r}) = - \sum_K \frac{Z_K e^2}{4\pi\varepsilon_0 |\mathbf{r} - \mathbf{R}_K|}
(3)

where ZKZ_K is the proton number of atom KK in the system.

Solving the Schrödinger equation for molecular systems is, however, a non-trivial task and requires the introduction of approximations to the NN-electron wave function. In contrast, the focus in DFT is set on the one-electron density

n0(r)=N∫∣Ψ0(r,r2,..,rN)∣2d3r2⋯d3rNn_0(\mathbf{r}) = N \int |\Psi_0(\mathbf{r}, \mathbf{r}_2,.., \mathbf{r}_N)|^2 \mathrm{d}^3\mathbf{r}_2 \cdots \mathrm{d}^3\mathbf{r}_N
(4)

which only depends on three variables namely the Cartesian components of the position vector.

The first attempts to describe ground states of atoms in terms of electron density were made in the 1920s by Thomas and Fermi. But the use of the density as the primary variable to describe the electronic ground state was not legitimized until 1964, when Hohenberg and Kohn (HK) introduced their theorems for vv-representable electron densities Hohenberg & Kohn (1964).

Considering all possible variations in the external potential (system variations), we understand that there is an infinite number of vv-representable densities associated with NN-electron systems. We refer to a general density in this set as n~0(r)\tilde{n}_0(\mathbf{r}).

HK theorem I

The external potential v^(r)\hat{v}(\mathbf{r}) of a given system is determined to within a trivial additive constant by the vv-representable electron density of the system.

In other words, there is a unique one-to-one mapping between a ground-state wave function and its one-electron density, and the knowledge of n0(r)n_0(\mathbf{r}) is sufficient to determine the ground-state energy (and other ground-state properties) of the molecular system.

This allows for the separation of the energy into separate the energy functional into two terms Hohenberg & Kohn (1964)

E0=E[n0(r)]=FHK[n0(r)]+V[n0(r)]E_0 = E[n_0(\mathbf{r})] = F_\mathrm{HK}[n_0(\mathbf{r})] + V[n_0(\mathbf{r})]
(5)

where

V[n0(r)]=∫v(r)n0(r) d3rV[n_0(\mathbf{r})] = \int v(\mathbf{r}) n_0(\mathbf{r}) \, \mathrm{d}^3\mathbf{r}
(6)

and the Hohenberg–Kohn functional is introduced as the sum of kinetic and electron repulsion energies

FHK[n0(r)]=T[n0(r)]+U[n0(r)]F_\mathrm{HK}[n_0(\mathbf{r})] = T[n_0(\mathbf{r})] + U[n_0(\mathbf{r})]
(7)

The HK functional is universal in the sense that it does not depend on the system under study as such dependencies are isolated to the external potential, and it is left undefined for densities that are not vv-representable.

HK theorem II

From the vv-representable trial densities n~0(r)\tilde{n}_0(\mathbf{r}) fulfilling

∫n~0(r) d3r=N;n~0(r)≥0\int \tilde{n}_0(\mathbf{r}) \, \mathrm{d}^3\mathbf{r} = N ; \qquad \tilde{n}_0(\mathbf{r}) \geq 0

the ground state energy E0E_0 of a molecular system can be detemined from the relation

E0≤E[n~0(r)]E_0 \leq E[\tilde{n}_0(\mathbf{r})]

We recognize this relation as the variational principle in wave function theory. Any practical application of this relation in density functional theory is, however, severely hampered by the fact that it is prohibitively difficult to ensure that density variations remain vv-representable.

NN-representability

The Hohenberg–Kohn theorems provide a theoretical foundation of DFT. Still, they do not give a recipe for the practical implementation of a computational scheme due to the strict requirement of vv-representability of densities. Fortunately, the theory can be reformulated on the grounds of the Hohenberg–Kohn theorems but for the wider class of so-called NN-representable densities.

For a given NN-representable density, Levy demonstrated that there exists a universal variational functional that delivers the associated sum of kinetic and repulsion energies. This functional can be determined by means of a constrained search over the set of wave functions that yield this density Levy (1979)

F[n(r)]=min⁡Ψ→n(r)⟨Ψ∣T^+U^∣Ψ⟩F[n(\mathbf{r})] = \min_{\Psi \to n(\mathbf{r})} \langle \Psi | \hat{T} +\hat{U}| \Psi \rangle
(8)

When applied to vv-representable densities, it was shown that FF becomes equal to the HK functional such that

F[n0(r)]+V[n0(r)]=E0F[n_0(\mathbf{r})] + V[n_0(\mathbf{r})] = E_0
(9)

and, for NN-representable densities, it was shown that

F[n(r)]+V[n(r)]≥E0F[n(\mathbf{r})] + V[n(\mathbf{r})] \geq E_0
(10)

This implies a variational principle with respect to NN-representable densities and the minimization condition for the energy functional can be written in terms of the variation of a Lagrangian that preserves the number of electrons

δ[E[n(r)]+μ(N−∫n(r) d3r)]=0\delta \Big[ E[n(\mathbf{r})] + \mu \Big( N - \int n(\mathbf{r}) \, d^3\mathbf{r} \Big) \Big] = 0
(11)

where the Lagrange multiplier μ\mu is the chemical potential of the molecular system, i.e., μ=dE/dN\mu = \mathrm{d}E/\mathrm{d}N. This stationary condition can alternatively be written

μ=v(r)+δF[n(r)]δn(r)\mu = v(\mathbf{r}) + \frac{\delta F[n(\mathbf{r})]}{\delta n(\mathbf{r})}
(12)
References
  1. Hohenberg, P., & Kohn, W. (1964). Inhomogeneous Electron Gas. Phys. Rev., 136, B864–B871. 10.1103/PhysRev.136.B864
  2. Levy, M. (1979). Universal variational functionals of electron densities, first-order density matrices, and natural spin-orbitals and solution of the v-representability problem. Proc. Natl. Acad. Sci. U.S.A., 76, 6062–6065. 10.1073/pnas.76.12.606
Electronic ground states
Density functional theory
Electronic ground states
Kohn–Sham method