Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review A

Axiomatic formulations of the Hohenberg-Kohn functional

Paul W. Ayers

Published Jan 31, 2006• 57 citations
Publisher

Abstract

Properties of the Hohenberg-Kohn functional are considered. In particular, the Hohenberg-Kohn functional should (a) give correct results in the variational principle and should be (b) continuous, (c) convex, and (d) size consistent. All of these properties are satisfied by the Legendre-transform functional (equivalently, the density matrix constrained search functional) and, moreover, this is the only functional that possesses all these properties. Not only that, but the Legendre-transform functional is determined uniquely by requiring (a), (b), and either (c) or (d). This shows how an ``axiomatic'' approach to constructing the Hohenberg-Kohn functional leads naturally to the Legendre-transform functional. Among all functionals consistent with the variational principle, the Legendre-transform functional is the smallest. One corollary to this approach is a simple proof of the equivalence of the Legendre-transform and density-matrix constrained search functionals. For completeness, the Appendix shows that ensemble-$v$-representable densities lie dense in the set of $N$-representable densities.

DOI
10.1103/physreva.73.012513

Related papers

Generalizations of the Hohenberg-Kohn theorem: I. Legendre Transform Constructions of Variational Principles for Density Matrices and Electron Distribution Functions

Paul W. Ayers, Sidney Golden, Mel Levy

Given a general, N-particle Hamiltonian operator, analogs of the Hohenberg-Kohn theorem are derived for functions that are more general than the particle density, including density matrices and the diagonal elements thereof. The generalization of Lieb's Legendre transform ansatz to the generalized Hohenberg-Kohn functional not only solves the upsilon-representability problem for these entities, but, more importantly, also solves the N-representability problem. Restricting the range of operators explored by the Legendre transform leads to a lower bound on the true functional. If all the operators of interest are incorporated in the restricted maximization, however, the variational principle dictates that exact results are obtained for the systems of interest. This might have important implications for practical work not only for density matrices but also for density functionals. A follow-up paper will present a useful alternative approach to the upsilon- and N-representability problems based on the constrained search formalism.

The Journal of Chemical Physics•Jan 27, 2006•92 citations

Density-functional theory on graphs

Markus Penz, Robert van Leeuwen

The principles of density-functional theory are studied for finite lattice systems represented by graphs. Surprisingly, the fundamental Hohenberg-Kohn theorem is found void, in general, while many insights into the topological structure of the density-potential mapping can be won. We give precise conditions for a ground state to be uniquely v-representable and are able to prove that this property holds for almost all densities. A set of examples illustrates the theory and demonstrates the non-convexity of the pure-state constrained-search functional.

arXiv:2106.15370•The Journal of Chemical Physics•Dec 7, 2021•35 citations

Exact density-functional theory as parallel ensemble variational hierarchies: from Lieb's formulation to Kohn-Sham theory

Nan Sheng

Exact density-functional theory is recast here as two parallel exact ensemble variational hierarchies: an interacting hierarchy rooted in Lieb's ensemble formulation and a noninteracting hierarchy rooted in exact noninteracting ensemble theory. In optimization terms, $N$-representability is primal feasibility, Legendre-Fenchel duality equates the primal and dual values, $v$-representability is dual attainment, and the Hohenberg-Kohn theorem gives uniqueness, modulo constants, of an attained local potential. The Kohn-Sham construction couples the interacting density-space optimality condition to a compatible noninteracting dual realization on their common $N$-representable density domain. State-class restrictions yield the Levy-Lieb and single-determinant branches, while fractional particle number and fractional occupations lead naturally to piecewise linearity, one-sided chemical potentials, Janak-type relations, and the derivative discontinuity. This organization locates the exactness of Kohn-Sham theory in the preservation of the interacting density-space optimization together with its compatible noninteracting realization, without implying a general many-body spectral interpretation of Kohn-Sham eigenvalues.

arXiv:2603.23399•Mar 24, 2026•0 citations

The Hohenberg–Kohn theorem

Carsten A. Ullrich, Adam Wasserman

Abstract In this chapter, the Hohenberg-Kohn theorem of DFT is formulated and proved. This theorem says that there is a unique one-to-one correspondence between potentials and ground-state densities in electronic many-body systems. The original proof of Hohenberg and Kohn is presented and discussed, which proceeds in two steps: first, the mapping from potentials to wave functions is shown to be unique, followed by the mapping from wave functions to densities. Two corollaries of the theorem are then discussed: the variational principle, whereby the exact ground-state energy can be found by variation of the energy functional, and then the universality of parts of this energy functional. The so-called constrained search formalism is also briefly discussed.

•Aug 20, 2026•0 citations

Elementary properties of an energy functional of the first-order reduced density matrix

Robert A. Donnelly, Robert G. Parr

The Hohenberg–Kohn theorem implies the existence of an energy functional based solely on the first-order reduced density matrix of the ground state of an atomic or molecular system. Application of the variational principle for the functional EvM[γ] generates a set of coupled Euler equations for the representation coefficients and spin orbitals of a rank-M approximation to the exact ground-state density matrix. Defining the (assumed Hermitian) kernel FM[γ;x′,x]≡δEvM/δγ (x,x′), the equations in an arbitrary representation for the approximate density matrix, γ (x,x′) =ΣijMψ1(x) γijψj* (x′) are the following: ℱdξFM[γ;x′,ξ]ψi(ξ) =μψi(x′); i=1, 2,...,M; ℱdξFM[γ;x′,ξ] ΣMiψi(ξ) γ ij=ΣMiψi(x′) λij; j=1, 2, ...,M. The quantity μ is the chemical potential of the system of interest and the λij are a set of M2 Lagrange multipliers constraining the orthonormality of the spin orbital basis {ψ}. The coefficients γij must be chosen such that FM has the degenerate eigenvalue spectrum, FiiM=μ, i=1, 2,...,M, for all partially occupied orbitals. The spin orbitals, for a fixed set of coefficients, must be simultaneously determined so that λ is Hermitian. Provided all occupation numbers lie on the open interval (0,1), the stationary density matrix itself obeys the eigenvalue equation ℱdξFM[γ;x,ξ]γ (ξ,x′) =μγ (x,x′), and for any stationary density matrix the following commutation rules are valid: [FM[γ], γ]=0; [λ, γ]=0. The matrices γ and λ consequently can be brought simultaneously to diagonal form, and the canonical representation of the energy functional provides an eigenvalue equation determining the natural spin orbitals. The value of the chemical potential is μ=εi/ni, i=1, 2, ...,M; ε=diag (λ); from which follows the distribution function governing the occupation numbers of a stationary density matrix. Two limiting forms of the variational principle are examined, the exact and Hartree–Fock functionals, and related previous work by Gilbert is discussed. The physical content of the equations is illuminated by identification of the chemical potential as the negative of the electronegativity; Sanderson’s Principle of Electronegativity Equalization follows.

The Journal of Chemical Physics•Nov 15, 1978•285 citations

ON THE N-REPRESENTABILITY AND UNIVERSALITY OF F[ρ] IN THE HOHENBERG-KOHN-SHAM VERSION OF DENSITY FUNCTIONAL THEORY

Eduardo V. Ludeña, Francesc Illas, Alejandro Ramírez-Solís

We discuss two basic problems in the Hohenberg-Kohn-Sham version of density functional theory, HKS-DFT: the first, the N-representability of the functional F[ρ] and, the second, the universality of F[ρ]. In relation to the first, we show that F[ρ] must satisfy N-representability conditions that follow from those on the 2-matrix D2(r1, r2; r′, r′2). In the case of the second, we provide arguments based on the equivalence between ab initio DFT and HKS-DFT to show that the functional F[ρ] is not universal.

International Journal of Modern Physics B•Oct 20, 2008•12 citations
Properties of the Hohenberg-Kohn functional are considered. In particular, the Hohenberg-Kohn functional should (a) give correct results in the variational principle and should be (b) continuous, (c) convex, and (d) size consistent. All of these properties are satisfied by the Legendre-transform functional (equivalently, the density matrix constrained search functional) and, moreover, this is the only functional that possesses all these properties. Not only that, but the Legendre-transform functional is determined uniquely by requiring (a), (b), and either (c) or (d). This shows how an ``axiomatic'' approach to constructing the Hohenberg-Kohn functional leads naturally to the Legendre-transform functional. Among all functionals consistent with the variational principle, the Legendre-transform functional is the smallest. One corollary to this approach is a simple proof of the equivalence of the Legendre-transform and density-matrix constrained search functionals. For completeness, the Appendix shows that ensemble-vvv-representable densities lie dense in the set of NNN-representable densities.
Given a general, N-particle Hamiltonian operator, analogs of the Hohenberg-Kohn theorem are derived for functions that are more general than the particle density, including density matrices and the diagonal elements thereof. The generalization of Lieb's Legendre transform ansatz to the generalized Hohenberg-Kohn functional not only solves the upsilon-representability problem for these entities, but, more importantly, also solves the N-representability problem. Restricting the range of operators explored by the Legendre transform leads to a lower bound on the true functional. If all the operators of interest are incorporated in the restricted maximization, however, the variational principle dictates that exact results are obtained for the systems of interest. This might have important implications for practical work not only for density matrices but also for density functionals. A follow-up paper will present a useful alternative approach to the upsilon- and N-representability problems based on the constrained search formalism.
The principles of density-functional theory are studied for finite lattice systems represented by graphs. Surprisingly, the fundamental Hohenberg-Kohn theorem is found void, in general, while many insights into the topological structure of the density-potential mapping can be won. We give precise conditions for a ground state to be uniquely v-representable and are able to prove that this property holds for almost all densities. A set of examples illustrates the theory and demonstrates the non-convexity of the pure-state constrained-search functional.
Exact density-functional theory is recast here as two parallel exact ensemble variational hierarchies: an interacting hierarchy rooted in Lieb's ensemble formulation and a noninteracting hierarchy rooted in exact noninteracting ensemble theory. In optimization terms, NNN-representability is primal feasibility, Legendre-Fenchel duality equates the primal and dual values, vvv-representability is dual attainment, and the Hohenberg-Kohn theorem gives uniqueness, modulo constants, of an attained local potential. The Kohn-Sham construction couples the interacting density-space optimality condition to a compatible noninteracting dual realization on their common NNN-representable density domain. State-class restrictions yield the Levy-Lieb and single-determinant branches, while fractional particle number and fractional occupations lead naturally to piecewise linearity, one-sided chemical potentials, Janak-type relations, and the derivative discontinuity. This organization locates the exactness of Kohn-Sham theory in the preservation of the interacting density-space optimization together with its compatible noninteracting realization, without implying a general many-body spectral interpretation of Kohn-Sham eigenvalues.
Abstract In this chapter, the Hohenberg-Kohn theorem of DFT is formulated and proved. This theorem says that there is a unique one-to-one correspondence between potentials and ground-state densities in electronic many-body systems. The original proof of Hohenberg and Kohn is presented and discussed, which proceeds in two steps: first, the mapping from potentials to wave functions is shown to be unique, followed by the mapping from wave functions to densities. Two corollaries of the theorem are then discussed: the variational principle, whereby the exact ground-state energy can be found by variation of the energy functional, and then the universality of parts of this energy functional. The so-called constrained search formalism is also briefly discussed.
The Hohenberg–Kohn theorem implies the existence of an energy functional based solely on the first-order reduced density matrix of the ground state of an atomic or molecular system. Application of the variational principle for the functional EvM[γ] generates a set of coupled Euler equations for the representation coefficients and spin orbitals of a rank-M approximation to the exact ground-state density matrix. Defining the (assumed Hermitian) kernel FM[γ;x′,x]≡δEvM/δγ (x,x′), the equations in an arbitrary representation for the approximate density matrix, γ (x,x′) =ΣijMψ1(x) γijψj* (x′) are the following: ℱdξFM[γ;x′,ξ]ψi(ξ) =μψi(x′); i=1, 2,...,M; ℱdξFM[γ;x′,ξ] ΣMiψi(ξ) γ ij=ΣMiψi(x′) λij; j=1, 2, ...,M. The quantity μ is the chemical potential of the system of interest and the λij are a set of M2 Lagrange multipliers constraining the orthonormality of the spin orbital basis {ψ}. The coefficients γij must be chosen such that FM has the degenerate eigenvalue spectrum, FiiM=μ, i=1, 2,...,M, for all partially occupied orbitals. The spin orbitals, for a fixed set of coefficients, must be simultaneously determined so that λ is Hermitian. Provided all occupation numbers lie on the open interval (0,1), the stationary density matrix itself obeys the eigenvalue equation ℱdξFM[γ;x,ξ]γ (ξ,x′) =μγ (x,x′), and for any stationary density matrix the following commutation rules are valid: [FM[γ], γ]=0; [λ, γ]=0. The matrices γ and λ consequently can be brought simultaneously to diagonal form, and the canonical representation of the energy functional provides an eigenvalue equation determining the natural spin orbitals. The value of the chemical potential is μ=εi/ni, i=1, 2, ...,M; ε=diag (λ); from which follows the distribution function governing the occupation numbers of a stationary density matrix. Two limiting forms of the variational principle are examined, the exact and Hartree–Fock functionals, and related previous work by Gilbert is discussed. The physical content of the equations is illuminated by identification of the chemical potential as the negative of the electronegativity; Sanderson’s Principle of Electronegativity Equalization follows.
We discuss two basic problems in the Hohenberg-Kohn-Sham version of density functional theory, HKS-DFT: the first, the N-representability of the functional F[ρ] and, the second, the universality of F[ρ]. In relation to the first, we show that F[ρ] must satisfy N-representability conditions that follow from those on the 2-matrix D2(r1, r2; r′, r′2). In the case of the second, we provide arguments based on the equivalence between ab initio DFT and HKS-DFT to show that the functional F[ρ] is not universal.