We establish one-body reduced density matrix (1RDM) functional theory for the canonical ensemble in a finite basis set at an elevated temperature. Including temperature guarantees the differentiability of the universal functional by occupying all states and additionally not fully occupying the states in a fermionic system. We use the convexity of the universal functional and invertibility of the potential-to-1RDM map to show that the subgradient contains only one element which is equivalent to differentiability. This allows us to show that all 1RDMs with a purely fractional occupation number spectrum ($0<{n}_{i}<1\phantom{\rule{0.28em}{0ex}}{\ensuremath{\forall}}_{i}$) are uniquely $v$ -representable up to a constant.We establish one-body reduced density matrix (1RDM) functional theory for the canonical ensemble in a finite basis set at an elevated temperature. Including temperature guarantees the differentiability of the universal functional by occupying all states and additionally not fully occupying the states in a fermionic system. We use the convexity of the universal functional and invertibility of the potential-to-1RDM map to show that the subgradient contains only one element which is equivalent to differentiability. This allows us to show that all 1RDMs with a purely fractional occupation number spectrum ($0<{n}_{i}<1\phantom{\rule{0.28em}{0ex}}{\ensuremath{\forall}}_{i}$) are uniquely $v$-representable up to a constant.

arXiv:2209.11663Physical Review AFeb 10, 20239 citations

Density functional theory at finite temperatures

Mario V. Stoitsov, Ivan Zh. Petkov

Annals of PhysicsMay 1, 198856 citations

Approximate energy functionals for one-body reduced density matrix functional theory from many-body perturbation theory

Klaas J. H. Giesbertz, Anna-Maija Uimonen, Robert van Leeuwen

We develop a systematic approach to construct energy functionals of the one-particle reduced density matrix (1RDM) for equilibrium systems at finite temperature. The starting point of our formulation is the grand potential Ω [ G ] regarded as variational functional of the Green’s function G based on diagrammatic many-body perturbation theory and for which we consider either the Klein or Luttinger–Ward form. By restricting the input Green’s function to be one-to-one related to a set on one-particle reduced density matrices (1RDM) this functional becomes a functional of the 1RDM. To establish the one-to-one mapping we use that, at any finite temperature and for a given 1RDM γ in a finite basis, there exists a non-interacting system with a spatially non-local potential v [ γ ] which reproduces the given 1RDM. The corresponding set of non-interacting Green’s functions defines the variational domain of the functional Ω . In the zero temperature limit we obtain an energy functional E [ γ ] which by minimisation yields an approximate ground state 1RDM and energy. As an application of the formalism we use the Klein and Luttinger–Ward functionals in the GW-approximation to compute the binding curve of a model hydrogen molecule using an extended Hubbard Hamiltonian. We compare further to the case in which we evaluate the functionals on a Hartree–Fock and a Kohn–Sham Green’s function. We find that the Luttinger–Ward version of the functionals performs the best and is able to reproduce energies close to the GW energy which corresponds to the stationary point.We develop a systematic approach to construct energy functionals of the one-particle reduced density matrix (1RDM) for equilibrium systems at finite temperature. The starting point of our formulation is the grand potential Ω [ G ] regarded as variational functional of the Green’s function G based on diagrammatic many-body perturbation theory and for which we consider either the Klein or Luttinger–Ward form. By restricting the input Green’s function to be one-to-one related to a set on one-particle reduced density matrices (1RDM) this functional becomes a functional of the 1RDM. To establish the one-to-one mapping we use that, at any finite temperature and for a given 1RDM γ in a finite basis, there exists a non-interacting system with a spatially non-local potential v [ γ ] which reproduces the given 1RDM. The corresponding set of non-interacting Green’s functions defines the variational domain of the functional Ω . In the zero temperature limit we obtain an energy functional E [ γ ] which by minimisation yields an approximate ground state 1RDM and energy. As an application of the formalism we use the Klein and Luttinger–Ward functionals in the GW-approximation to compute the binding curve of a model hydrogen molecule using an extended Hubbard Hamiltonian. We compare further to the case in which we evaluate the functionals on a Hartree–Fock and a Kohn–Sham Green’s function. We find that the Luttinger–Ward version of the functionals performs the best and is able to reproduce energies close to the GW energy which corresponds to the stationary point.

arXiv:1804.09921The European Physical Journal BNov 1, 201815 citations

Static correlated functionals for reduced density matrix functional theory

Carlos L. Benavides-Riveros, Miguel A. L. Marques

arXiv:1810.09039The European Physical Journal BJun 1, 201814 citations

The Interface Between Reduced Density Matrices and Density Functional Theory

John E. Harriman

Springer eBooksJan 1, 198710 citations

v

-representable up to a constant.