{
"cells": [
{
"cell_type": "markdown",
"id": "6a6c2169",
"metadata": {},
"source": [
"# Hohenberg–Kohn theorems\n",
"\n",
"To illustrate the basic idea of DFT, let us consider an isolated molecular system with $N$ electrons. The anti-symmetric $N$-electron wave function describing the electronic ground state of this system, $\\Psi_0(\\mathbf{r}_1, \\mathbf{r}_2,.., \\mathbf{r}_N) $, depends on the spatial coordinates of the electrons, $\\{\\mathbf{r}_i\\}$, and it is a solution to the time-independent Schrödinger equation\n",
"\n",
"\\begin{equation*}\n",
"\\label{eq:sch}\n",
"\\hat{H} \\Psi_0 = E_0 \\Psi_0 \n",
"\\end{equation*}\n",
"\n",
"where $E_0$ is the associated electronic energy and $\\hat{H}$ is the electronic Hamiltonian that is composed of the kinetic energy, electron–electron repulsion, and external potential operators\n",
"\n",
"\\begin{equation*}\n",
"\\hat{H} = \\hat{T} + \\hat{U} + \\hat{V}\n",
"\\end{equation*}\n",
"\n",
"with\n",
"\n",
"\\begin{equation*}\n",
"\\hat{T} = - \\frac{\\hbar^2}{2 m_\\mathrm{e}} \\sum_i \\nabla_i^2 ; \\qquad\n",
"%\n",
"\\hat{U} = \\sum_{i The external potential $\\hat{v}(\\mathbf{r})$ of a given system is determined to within a trivial additive constant by the $v$-representable electron density of the system.\n",
"\n",
"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 $n_0(\\mathbf{r})$ is sufficient to determine the ground-state energy (and other ground-state properties) of the molecular system.\n",
"\n",
"This allows for the separation of the energy into separate the energy functional into two terms {cite}`Hohenberg1964`\n",
"\n",
"$$\n",
"E_0 =\n",
"E[n_0(\\mathbf{r})] =\n",
"F_\\mathrm{HK}[n_0(\\mathbf{r})] + V[n_0(\\mathbf{r})]\n",
"$$\n",
"\n",
"where\n",
"\n",
"$$\n",
"V[n_0(\\mathbf{r})] =\n",
"\\int v(\\mathbf{r}) n_0(\\mathbf{r}) \\, \\mathrm{d}^3\\mathbf{r}\n",
"$$\n",
"\n",
"and the Hohenberg–Kohn functional is introduced as the sum of kinetic and electron repulsion energies\n",
"\n",
"$$\n",
"F_\\mathrm{HK}[n_0(\\mathbf{r})] = T[n_0(\\mathbf{r})] + U[n_0(\\mathbf{r})]\n",
"$$ \n",
"\n",
"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 $v$-representable.\n",
"\n",
"## HK theorem II\n",
"\n",
"> From the $v$-representable trial densities $\\tilde{n}_0(\\mathbf{r})$ fulfilling\n",
">\n",
"> \\begin{equation*}\n",
"\\int \\tilde{n}_0(\\mathbf{r}) \\, \\mathrm{d}^3\\mathbf{r} = N ;\n",
"\\qquad\n",
"\\tilde{n}_0(\\mathbf{r}) \\geq 0\n",
"\\end{equation*}\n",
">\n",
"> the ground state energy $E_0$ of a molecular system can be detemined from the relation\n",
">\n",
"> \\begin{equation*}\n",
" E_0 \\leq E[\\tilde{n}_0(\\mathbf{r})]\n",
"\\end{equation*}\n",
"\n",
"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 $v$-representable. \n",
"\n",
"(N-representability)=\n",
"## $N$-representability\n",
"\n",
"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 $v$-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 $N$-representable densities.\n",
"\n",
"```{note}\n",
"An electron density is $N$-representable if it can be obtained from some anti-symmetric wave function. Such a density satisfies the following conditions \n",
"\\begin{equation*}\n",
"\\int n(\\mathbf{r}) \\, d^3\\mathbf{r} = N ; \\quad \n",
"n(\\mathbf{r}) \\geq 0; \\quad \n",
"\\int \\left| \\nabla n(\\mathbf{r})^{1/2} \\right|^{2} \\, d^3\\mathbf{r} < \\infty\n",
"\\end{equation*}\n",
"```\n",
"\n",
"For a given $N$-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 {cite}`Levy1979`\n",
"\n",
"\\begin{equation*}\n",
"F[n(\\mathbf{r})] = \\min_{\\Psi \\to n(\\mathbf{r})} \\langle \\Psi | \\hat{T} +\\hat{U}| \\Psi \\rangle \n",
"\\end{equation*}\n",
"\n",
"When applied to $v$-representable densities, it was shown that $F$ becomes equal to the HK functional such that\n",
"\n",
"$$\n",
"F[n_0(\\mathbf{r})] + V[n_0(\\mathbf{r})] = E_0\n",
"$$\n",
"\n",
"and, for $N$-representable densities, it was shown that\n",
"\n",
"$$\n",
"F[n(\\mathbf{r})] + V[n(\\mathbf{r})] \\geq E_0\n",
"$$\n",
"\n",
"This implies a variational principle with respect to $N$-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\n",
"\n",
"\\begin{equation*}\n",
"\\delta \\Big[ E[n(\\mathbf{r})] + \\mu \\Big( N - \\int n(\\mathbf{r}) \\, d^3\\mathbf{r} \\Big) \\Big] = 0 \n",
"\\end{equation*}\n",
"\n",
"where the Lagrange multiplier $\\mu$ is the chemical potential of the molecular system, i.e., $\\mu = \\mathrm{d}E/\\mathrm{d}N$. This stationary condition can alternatively be written\n",
"\n",
"\\begin{equation*}\n",
"\\mu = v(\\mathbf{r}) + \\frac{\\delta F[n(\\mathbf{r})]}{\\delta n(\\mathbf{r})} \n",
"\\end{equation*}\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "d6ade0d2-6321-4580-920a-1430c3411f4b",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.9.12"
}
},
"nbformat": 4,
"nbformat_minor": 5
}