For properties of interacting electron systems, Kohn-Sham (KS) theory is often favored over many-body perturbation theory (MBPT), owing to its low computational cost. However, the exact KS potential can be challenging to approximate, for example in the presence of localized subsystems where the exact potential is known to exhibit pathological features such as spatial steps. By modeling two electrons, each localized in a distinct potential well, we illustrate that the step feature has no counterpart in MBPTs (including Hartree-Fock and $GW)$ or hybrid methods involving Fock exchange because the spatial nonlocality of the self-energy renders such pathological behavior unnecessary. We present a quantitative illustration of the orbital-dependent nature of the nonlocal potential, and a numerical demonstration of Kohn's concept of the nearsightedness for self-energies, when two distant subsystems are combined, in contrast to the KS potential. These properties emphasize the value of self-energy-based approximations in developing future approaches within KS-like theories.
For properties of interacting electron systems, Kohn-Sham (KS) theory is often favored over many-body perturbation theory (MBPT), owing to its low computational cost. However, the exact KS potential can be challenging to approximate, for example in the presence of localized subsystems where the exact potential is known to exhibit pathological features such as spatial steps. By modeling two electrons, each localized in a distinct potential well, we illustrate that the step feature has no counterpart in MBPTs (including Hartree-Fock and GW) or hybrid methods involving Fock exchange because the spatial nonlocality of the self-energy renders such pathological behavior unnecessary. We present a quantitative illustration of the orbital-dependent nature of the nonlocal potential, and a numerical demonstration of Kohn's concept of the nearsightedness for self-energies, when two distant subsystems are combined, in contrast to the KS potential. These properties emphasize the value of self-energy-based approximations in developing future approaches within KS-like theories.
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.
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.
Aiming to combine density functional theory (DFT) and wave-function theory, we study a mapping from the many-body interacting system to an effectively interacting Kohn-Sham system instead of a noninteracting Kohn-Sham system. Because a ground state of effectively interacting systems requires having a solution for the correlated many-body wave functions, this provides a natural framework to many-body wave-function theories such as the configuration interaction and the coupled-cluster method in the formal theoretical framework of DFT. Employing simple one-dimensional two-electron systems---namely, the one-dimensional helium atom, the hydrogen molecule, and the heteronuclear diatomic molecule---we investigate properties of many-body wave functions and exact exchange-correlation potentials of effectively interacting Kohn-Sham systems. As a result, we find that the asymptotic behavior of the exact exchange-correlation potential can be controlled by optimizing that of the effective interaction. Furthermore, the typical features of the exact noninteracting Kohn-Sham system, namely, a spiky feature and a step feature in the exchange-correlation potential for the molecular dissociation limit, can be suppressed by a proper choice of the effective interaction. These findings open a possibility to construct numerically robust and efficient exchange-correlation potentials and functionals based on the effectively interacting Kohn-Sham scheme.
We use the exact strong-interaction limit of the Hohenberg-Kohn energy density functional to approximate the exchange-correlation energy of the restricted Kohn-Sham scheme. Our approximation corresponds to a highly nonlocal density functional whose functional derivative can be easily constructed, thus transforming exactly, in a physically transparent way, an important part of the electron-electron interaction into an effective local one-body potential. We test our approach on quasi-one-dimensional systems, showing that it captures essential features of strong correlation that restricted Kohn-Sham calculations using the currently available approximations cannot describe.
Many problems in physical chemistry involve systems that are coupled to an environment, such as a molecule interacting with a nearby surface. In these cases, electron density may flow between the system and its surroundings, leading to metastable states and other open-system phenomena that cannot be captured by conventional Hermitian quantum mechanics. Non-Hermitian quantum mechanics (NHQM) [1] provides a natural framework for the description of open systems. Hartree-Fock (HF) and Kohn-Sham (KS) methods can be extended to the non-Hermitian domain [2–7]. Here, we show that even for nominally isolated systems, within HF and KS, each electron can be viewed as coupled to the remaining electrons, which serve as an environment. This viewpoint reveals that an electron can leave or rejoin the system, producing self-consistent complex potentials. These unconventional solutions are absent in the conventional closed-system HF and KS approach. We show that these additional solutions can be physically meaningful. Our results highlight the importance of incorporating non-Hermitian perspectives into electronic structure calculations, offering new insights and extending the range of problems HF and KS can be applied to. [1] Moiseyev, N. Non-Hermitian Quantum Mechanics, Cambridge University Press (2011). [2] Zhou, Y., Ernzerhof, M. J. Chem. Phys. 136, 094105 (2012). [3] Ernzerhof, M. J. Chem. Phys. 125, 124104 (2006). [4] Zhou, Y., Ernzerhof, M. J. Phys. Chem. Lett. 3, 1916–1920 (2012). [5] Burton, H. G. A., Thom, A. J. W. J. Chem. Theory Comput. 12, 167–173 (2016). [6] Burton, H. G. A., Gross, M., Thom, A. J. W. J. Chem. Theory Comput. 14, 607–618 (2018). [7] Burton, H. G. A., Thom, A. J. W., Loos, P.-F. J. Chem. Theory Comput. 15, 4374–4385 (2019).
Strongly correlated electrons have been the subject of substantial theoretical interest for many years. Most work has focused on obtaining the energy in a low-cost fashion. Here, we show that even methods with good energies can yield significant "delocalization errors" that affect the orbitals and density, leading to large errors in predicting other important properties such as dipole moments. We illustrate this point by comparing existing state-of-art approaches with an accurate exchange correlation functional based on a generalised valence bond ansatz, in which orbitals and fractional occupations are treated as variational parameters via a common optimized effective potential (OEP). We show that the OEP exhibits step and peak features which, similar to the exact Kohn-Sham (KS) potential of DFT, are crucial to prevent charge delocalization. We further show that the step is missing in common approximations within reduced density matrix functional theory resulting in delocalization errors comparable to those found in DFT approximations. Finally, we explain the delocalisation error as coming from an artificial mixing of the ground state with a charge-transfer excited state which is avoided if occupation numbers exhibit discontinuities.