© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:0711.2996 · Physical Review B

Generalized Kohn-Sham system in one-matrix functional theory

Ryan Requist, Oleg Pankratov

Published Jun 25, 2008• 37 citations

Abstract

A system of electrons in a local or nonlocal external potential can be described with one-matrix functional theory (1MFT), which is similar to density-functional theory (DFT) but takes the one-particle reduced density matrix (one-matrix) instead of the density as its basic variable. Within 1MFT, Gilbert derived [Phys. Rev. B 12, 2111 (1975)] effective single-particle equations analogous to the Kohn-Sham (KS) equations in DFT. The self-consistent solution of these 1MFT-KS equations reproduces not only the density of the original electron system but also its one-matrix. While in DFT it is usually possible to reproduce the density using KS orbitals with integer (0 or 1) occupancy, in 1MFT reproducing the one-matrix requires in general fractional occupancies. The variational principle implies that the KS eigenvalues of all fractionally occupied orbitals must collapse at self-consistency to a single level. We show that as a consequence of the degeneracy, the iteration of the KS equations is intrinsically divergent. Fortunately, the level-shifting method, commonly introduced in Hartree-Fock calculations, is always able to force convergence. We introduce an alternative derivation of the 1MFT-KS equations that allows control of the eigenvalue collapse by constraining the occupancies. As an explicit example, we apply the 1MFT-KS scheme to calculate the ground state one-matrix of an exactly solvable two-site Hubbard model.

DOI: 10.1103/physrevb.77.235121

Related papers

Generalization of the Kohn-Sham system enabling representing arbitary one electron density matrices

van Dam, J J Hubertus

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.

arXiv:1512.06120•Dec 7, 2015•0 citations

Density-matrix functional theory for the N-particle ground state

Gil Zumbach, Klaus Maschke

We discuss the one-particle density-matrix functional theory for the N-particle ground-state problem. Using the variational principle we obtain a set of self-consistent equations which could be useful for practical calculations. The derived relation between the eigenvalues of the density matrix and the chemical potential agrees with that of Gilbert. The present theory is compared with the Kohn–Sham and the Hartree–Fock approach to the ground state problem. It is found that both approaches can be obtained from the present theory under certain assumptions.

The Journal of Chemical Physics•Jun 15, 1985•72 citations

Accurate real-time evolution of electron densities and ground-state properties from generalized Kohn-Sham theory

M. J. P. Hodgson, Jack Wetherell

The exact static and time-dependent Kohn-Sham (KS) exchange-correlation potential is extremely challenging to approximate as it is a local multiplicative potential that depends on the electron density everywhere in the system. The KS approach can be generalized by allowing part of the potential to be spatially nonlocal. We take this nonlocal part to be that of unrestricted Hartree-Fock theory. The additional local correlation potential in principle ensures that the single-particle density exactly equals the many-body density. In our case, the local correlation potential is predominantly nearsighted in its dependence on the density and hence an (adiabatic) local-density approximation to this potential yields accurate ground-state properties and real-time densities for one-dimensional test systems.

arXiv:1909.02510•Physical Review A•Mar 2, 2020•3 citations

Variational principle for the ground-state energy as a functional of the one-particle density matrix: Beyond Hartree-Fock theory

Abraham Klein, Reiner M. Dreizler

In parallel with standard density-functional theory, we study the energy of the ground state of a finite many-body system as a functional of the one-particle density matrix. We show that the formulation of a variational principle that is valid beyond the Hartree-Fock limit requires that two-body correlations be included not only in the ground-state energy but also in the constraints. As an illustration, we apply a linear-response argument to derive formulas for first-order corrections to the Hartree-Fock density matrix. Further analysis suggests an approach in terms of the density matrix of an independent-particle system, which can be introduced by the application of an alternative variational principle. This approach is reminiscent of Kohn-Sham theory, but the effective external potential is not required to be local. This variational method can be implemented in a systematic fashion by means of the linked-cluster expansion. In an appendix we study a variant of the Hohenberg-Kohn theorem for nonlocal potentials.

Physical Review A•Apr 1, 1998•11 citations

Exposing minimal composition of Kohn–Sham theory and its extendability

Hitoshi Nakada

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.

arXiv:2108.01232•Physica Scripta•Aug 29, 2023•3 citations

Excited-state-specific Kohn–Sham formalism for the asymmetric Hubbard dimer

Pierre‐François Loos, Sara Giarrusso

Building on our recent study [Giarrusso and Loos, J. Phys. Chem. Lett. 14, 8780 (2023)], we explore the generalization of the ground-state Kohn-Sham (KS) formalism of density-functional theory (DFT) to the (singlet) excited states of the asymmetric Hubbard dimer at half-filling. While we found that the KS-DFT framework can be straightforwardly generalized to the highest-lying doubly excited state, the treatment of the first excited state presents significant challenges. In particular, using a density-fixed adiabatic connection, we show that the density of the first excited state lacks non-interacting v-representability. However, by employing an analytic continuation of the adiabatic path, we demonstrate that the density of the first excited state can be generated by a complex-valued external potential in the non-interacting case. More practically, by performing state-specific KS calculations with exact and approximate correlation functionals-each state possessing a distinct correlation functional-we observe that spurious stationary solutions of the KS equations may arise due to the approximate nature of the functional.

arXiv:2412.14945•The Journal of Chemical Physics•Apr 8, 2025•6 citations

A system of electrons in a local or nonlocal external potential can be described with one-matrix functional theory (1MFT), which is similar to density-functional theory (DFT) but takes the one-particle reduced density matrix (one-matrix) instead of the density as its basic variable. Within 1MFT, Gilbert derived [Phys. Rev. B 12, 2111 (1975)] effective single-particle equations analogous to the Kohn-Sham (KS) equations in DFT. The self-consistent solution of these 1MFT-KS equations reproduces not only the density of the original electron system but also its one-matrix. While in DFT it is usually possible to reproduce the density using KS orbitals with integer (0 or 1) occupancy, in 1MFT reproducing the one-matrix requires in general fractional occupancies. The variational principle implies that the KS eigenvalues of all fractionally occupied orbitals must collapse at self-consistency to a single level. We show that as a consequence of the degeneracy, the iteration of the KS equations is intrinsically divergent. Fortunately, the level-shifting method, commonly introduced in Hartree-Fock calculations, is always able to force convergence. We introduce an alternative derivation of the 1MFT-KS equations that allows control of the eigenvalue collapse by constraining the occupancies. As an explicit example, we apply the 1MFT-KS scheme to calculate the ground state one-matrix of an exactly solvable two-site Hubbard model.

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.

We discuss the one-particle density-matrix functional theory for the N-particle ground-state problem. Using the variational principle we obtain a set of self-consistent equations which could be useful for practical calculations. The derived relation between the eigenvalues of the density matrix and the chemical potential agrees with that of Gilbert. The present theory is compared with the Kohn–Sham and the Hartree–Fock approach to the ground state problem. It is found that both approaches can be obtained from the present theory under certain assumptions.

The exact static and time-dependent Kohn-Sham (KS) exchange-correlation potential is extremely challenging to approximate as it is a local multiplicative potential that depends on the electron density everywhere in the system. The KS approach can be generalized by allowing part of the potential to be spatially nonlocal. We take this nonlocal part to be that of unrestricted Hartree-Fock theory. The additional local correlation potential in principle ensures that the single-particle density exactly equals the many-body density. In our case, the local correlation potential is predominantly nearsighted in its dependence on the density and hence an (adiabatic) local-density approximation to this potential yields accurate ground-state properties and real-time densities for one-dimensional test systems.

In parallel with standard density-functional theory, we study the energy of the ground state of a finite many-body system as a functional of the one-particle density matrix. We show that the formulation of a variational principle that is valid beyond the Hartree-Fock limit requires that two-body correlations be included not only in the ground-state energy but also in the constraints. As an illustration, we apply a linear-response argument to derive formulas for first-order corrections to the Hartree-Fock density matrix. Further analysis suggests an approach in terms of the density matrix of an independent-particle system, which can be introduced by the application of an alternative variational principle. This approach is reminiscent of Kohn-Sham theory, but the effective external potential is not required to be local. This variational method can be implemented in a systematic fashion by means of the linked-cluster expansion. In an appendix we study a variant of the Hohenberg-Kohn theorem for nonlocal potentials.

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.

Building on our recent study [Giarrusso and Loos, J. Phys. Chem. Lett. 14, 8780 (2023)], we explore the generalization of the ground-state Kohn-Sham (KS) formalism of density-functional theory (DFT) to the (singlet) excited states of the asymmetric Hubbard dimer at half-filling. While we found that the KS-DFT framework can be straightforwardly generalized to the highest-lying doubly excited state, the treatment of the first excited state presents significant challenges. In particular, using a density-fixed adiabatic connection, we show that the density of the first excited state lacks non-interacting v-representability. However, by employing an analytic continuation of the adiabatic path, we demonstrate that the density of the first excited state can be generated by a complex-valued external potential in the non-interacting case. More practically, by performing state-specific KS calculations with exact and approximate correlation functionals-each state possessing a distinct correlation functional-we observe that spurious stationary solutions of the KS equations may arise due to the approximate nature of the functional.