The formalism of Kohn and Sham uses a specific (model) hamiltonian which highly simplifies the many-electron problem to that of noninteracting fermions. The theorem of Hohenberg and Kohn tells us that, for a given ground state density, this hamiltonian is unique. In principle, this density can be chosen as that of the real, interacting system. To obtain the energy, or other properties of the real system, approximations are needed. Working with non interacting fermions is an important simplification, but it may be easier to produce approximations with different choices of the model hamiltonian. The feature that the exact density is (ideally) reproduced can be kept in the newly defined fictitious systems. Using model hamiltonians having the same form as the physical one, that is, being built of one- and two-body operators, allows to approach the physical hamiltonian arbitrarily close, and thus a systematic reduction of the approximations.
Density functional theory is currently the most widely applied method in electronic structure theory. The Kohn-Sham method, based on a fictitious system of non-interacting particles, is the work horse of the theory. The particular form of the Kohn-Sham wavefunction admits only idem-potent one electron density matrices whereas wavefunctions of correlated electrons in post-Hartree-Fock methods invariably have fractional occupation numbers. Here we show that by generalizing the orbital concept, and introducing a suitable dot-product as well as a probability density a non-interacting system can be chosen that can represent the one-electron density matrix of any system, even one with fractional occupation numbers. This fictitious system ensures that the exact electron density is accessible within density functional theory. It can also serve as the basis for reduced density matrix functional theory. Moreover, to aid the analysis of the results the orbitals may be assigned energies from a mean-field Hamiltonian. This produces energy levels that are akin to Hartree-Fock orbital energies such that conventional analyses based on Koopmans theorem are available. Finally, this system is convenient in formalisms that depend on creation and annihilation operators as they are trivially applied to single determinant wavefunctions.
The Hohenberg-Kohn theorem and the Kohn-Sham equations, which are at the basis of the density-functional theory, are reformulated in terms of a particular many-body density, which is translational and Galilean invariant and therefore is relevant for self-bound systems. In a similar way that there is a unique relation between the one-body density and the external potential that gives rise to it, we demonstrate that there is a unique relation between that particular many-body density and a definite many-body potential. The energy is then a functional of this density, and its minimization leads to the ground-state energy of the system. As a proof of principle, the analogous of the Kohn-Sham equation is solved in the specific case of
4He atomic clusters, to put in evidence the advantages of this formulation in terms of physical insights.
Abstract Reducing the many-fermion problem to a set of single-particle (s.p.) equations, the Kohn–Sham (KS) theory has provided a practical tool to implement ab initio calculations of ground-state energies and densities in many-electron systems. There have been attempts to extend the KS theory so that it could describe other physical quantities, or it could be applied to other many-fermion systems. By generalizing and reformulating the KS theory in terms of the 1-body density matrix, we expose the minimal composition of the theory that enables the reduction of the many-fermion problem to the s.p. equations. Based on the reformulation, several basic issues are reconsidered. The v- and N-representabilities for the KS theory are distinguished from those for the Hohenberg-Kohn theorem. Criteria for the extendability of the KS theory are addressed.
The analysis in this paper shows that the Hohenberg-Kohn theorem is the constellation of two statements: (i) the mathematically rigorous Hohenberg-Kohn lemma, which demonstrates that the same ground-state density cannot correspond to two different potentials of an external field, and (ii) the hypothesis of the existence of the universal density functional. Based on the obtained explicit expression for the nonrel-ativistic particle energy in a local external field, we prove that the energy of the system of more than two non-interacting electrons cannot be a functional of the inhomogeneous density. This result is generalized to the system of interacting electrons. It means that the Hohenberg-Kohn lemma cannot provide justification of the universal density functional for fermions. At the same time, statements of the density functional theory remain valid when considering any number of noninteracting ground-state bosons due to the Bose condensation effect. In the framework of the density matrix functional theory, the hypothesis of the existence of the universal density matrix functional corresponds to the cases of noninteracting particles and to interaction in the Hartree-Fock approximation.
By utilizing the knowledge that a Hamiltonian is a unique functional of its ground-state density, the following fundamental connections between densities and Hamiltonians are revealed: Given that ${\ensuremath{\rho}}_{\ensuremath{\alpha}}, {\ensuremath{\rho}}_{\ensuremath{\beta}},\dots{},{\ensuremath{\rho}}_{\ensuremath{\omega}}$ are ground-level densities for interacting or noninteracting Hamiltonians
H1,H2,…,HM (
M arbitrarily large) with local potentials
v1,
v2,
…,
vM, but given that we do not know which $\ensuremath{\rho}$ belongs with which
H, the correct mapping is possible and is obtained by minimizing $\ensuremath{\int}d\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}} [{v}_{1}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}){\ensuremath{\rho}}_{\ensuremath{\alpha}}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})+{v}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}){\ensuremath{\rho}}_{\ensuremath{\beta}}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})+\ensuremath{\cdots}{v}_{M}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}){\ensuremath{\rho}}_{\ensuremath{\omega}}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})]$ with respect to optimum permutations of the $\ensuremath{\rho}$'s among the
v's. A tight rigorous bound connects a density to its interacting ground-state energy via the one-body potential of the interacting system and the Kohn-Sham effective one-body potential of the auxiliary noninteracting system. A modified Kohn-Sham effective potential is defined such that its sum of lowest orbital energies equals the true interacting ground-state energy. Moreover, of all those effective potentials which differ by additive constants and which yield the true interacting ground-state density, this modified effective potential is the most invariant with respect to changes in the one-body potential of the true Hamiltonian. With the exception of the occurrence of certain linear dependencies,
a density will not generally be associated with any ground-state wave function (is not wave function
v representable) if that density can be generated by a special linear combination of three or more densities that arise from a common set of degenerate ground-state wave functions. Applicability of the "constrained search" approach to density-functional theory is emphasized for non-
v-representable as well as for
v-representable densities. In fact, a particular constrained ensemble search is revealed which provides a general sufficient condition for non-
v representability by a wave function. The possible appearance of noninteger occupation numbers is discussed in connection with the existence of non-
v representability for some Kohn-Sham noninteracting systems.
Based on the work of Görling and that of Levy and Nagy, a density-functional formalism for many fermionic excited states is explored through a careful and rigorous analysis of the excited-state density to external potential mapping. It is shown that knowledge of the ground-state density is a must to fix the mapping from an excited-state density to an external potential. This is the excited-state counterpart of the Hohenberg–Kohn theorem, where instead of the ground-state density the density of the excited state gives the true many-body wavefunctions of the system. Further, the excited-state Kohn–Sham system is defined by comparing its non-interacting kinetic energy with the true kinetic energy. The theory is demonstrated by studying a large number of atomic systems.