Functional derivative of the universal density functional in Fock space

Abstract

Within the framework of zero-temperature Fock-space density-functional theory (DFT), we prove that the G\^ateaux functional derivative of the universal density functional, ${\ensuremath{\delta}{F}^{\ensuremath{\lambda}}[\ensuremath{\rho}]∕\ensuremath{\delta}\ensuremath{\rho}(\mathbit{r})\ensuremath{\mid}}_{\ensuremath{\rho}={\ensuremath{\rho}}_{0}}$, at ground-state densities with arbitrary normalizations $(⟨{\ensuremath{\rho}}_{0}(\mathbit{r})⟩=n∊{\mathcal{R}}_{+})$ and an electron-electron interaction strength $\ensuremath{\lambda}$, is uniquely defined, but is discontinuous when the number of electrons $n$ becomes an integer, thus providing a mathematically rigorous confirmation for the ``derivative discontinuity'' initially discovered by Perdew et al. [Phys. Rev. Lett. 49, 1691 (1982)]. However, the functional derivative of the exchange-correlation functional is continuous with respect to the number of electrons in Fock space; i.e., there is no ``derivative discontinuity'' for the exchange-correlation functional at an integer electron number. For a ground-state density ${\ensuremath{\rho}}_{0,n}^{v,\ensuremath{\lambda}}(\mathbit{r})$ of an external potential $v(\mathbit{r})$, we show that ${\ensuremath{\delta}{F}^{\ensuremath{\lambda}}[\ensuremath{\rho}]∕\ensuremath{\delta}\ensuremath{\rho}(\mathbit{r})\ensuremath{\mid}}_{\ensuremath{\rho}={\ensuremath{\rho}}_{0,n}^{v,\ensuremath{\lambda}}}={\ensuremath{\mu}}_{\mathrm{SM}}^{n}\ensuremath{-}v(\mathbit{r})$, where the constant ${\ensuremath{\mu}}_{\mathrm{SM}}^{n}$ is given by the following chain of dependences: ${\ensuremath{\rho}}_{0,n}^{v,\ensuremath{\lambda}}(\mathbit{r})\ensuremath{\mapsto}[v]\ensuremath{\mapsto}{E}_{0}^{v,\ensuremath{\lambda}}(n)\ensuremath{\mapsto}{\ensuremath{\mu}}_{\mathrm{SM}}^{n}={\ensuremath{\partial}{E}_{0}^{v,\ensuremath{\lambda}}(k)∕\ensuremath{\partial}k\ensuremath{\mid}}_{k=n}$. Here $[v]$ is the class of the external potential $v(\mathbit{r})$ up to a real constant, and ${\ensuremath{\mu}}_{\mathrm{SM}}^{n}$ is the chemical potential defined according to statistical mechanics. At an integer electron number $N$, we find that there is no freedom of adding an arbitrary constant to the value of the chemical potential ${\ensuremath{\mu}}_{\mathrm{SM}}^{N}$, whose exact value is generally not the popular preference of the negative of Mulliken's electronegativity, $\ensuremath{-}\frac{1}{2}(I+A)$, where $I$ and $A$ are the first ionization potential and the first electron affinity, respectively. In addition, for any external potential converging to the same constant at infinity in all directions, we resolve that ${\ensuremath{\mu}}_{\mathrm{SM}}^{N}=\ensuremath{-}I$. Finally, the equality ${\ensuremath{\mu}}_{\mathrm{DFT}}={\ensuremath{\mu}}_{\mathrm{SM}}^{n}$ is rigorously derived via an alternative route, where ${\ensuremath{\mu}}_{\mathrm{DFT}}$ is the Lagrangian multiplier used to constrain the normalization of the density in the traditional DFT approach.