Ground-state energy as a simple sum of orbital energies in Kohn-Sham theory: a shift in perspective through a shift in potential

Extended N-centered ensemble density functional theory of double electronic excitations

A recent work (arXiv:2401.04685) has merged N-centered ensembles of neutral and charged electronic ground states with ensembles of neutral ground and excited states, thus providing a general and in-principle exact (so-called extended N-centered) ensemble density functional theory of neutral and charged electronic excitations. This formalism made it possible to revisit the concept of density-functional derivative discontinuity, in the particular case of single excitations from the highest occupied Kohn-Sham (KS) molecular orbital, without invoking the usual "asymptotic behavior of the density" argument. In this work, we address a broader class of excitations, with a particular focus on double excitations. An exact implementation of the theory is presented for the two-electron Hubbard dimer model. A thorough comparison of the true physical ground- and excited-state electronic structures with that of the fictitious ensemble density-functional KS system is also presented. Depending on the choice of the density-functional ensemble as well as the asymmetry of the dimer and the correlation strength, an inversion of states can be observed. In some other cases, the strong mixture of KS states within the true physical system makes the assignment "single excitation" or "double excitation" irrelevant.

Effective homogeneity of Fermi-Amaldi-containing exchange-correlation functionals

Parr and Ghosh [Phys. Rev. A. 51 3564 (1995)] demonstrated that when near-exact electron densities and potentials are used, the exchange-correlation energies of first- and second-row atoms are well-described by a combination of the Fermi-Amaldi functional with a functional that is homogeneous of degree one under density scaling. Insight into this observation is provided by considering their work from the perspective of the effective homogeneity of the overall exchange-correlation functional. By considering a general form that combines the Fermi-Amaldi functional with a functional that is homogeneous of degree k, it is shown that for these atoms, the functional of Parr and Ghosh (k = 1) exhibits essentially optimal effective homogeneities on the electron-deficient side of the integer. Percentage errors in effective homogeneities are close to percentage errors in exchange-correlation energies.

A Density Functional Theory for the Average Electron Energy

A formally exact density functional theory (DFT) determination of the average electron energy is presented. Our theory, which is based on a different accounting of energy functional terms, partially solves one well-known downside of conventional Kohn-Sham (KS) DFT: that electronic energies have but tenuous connections to physical quantities. Calculated average electron energies are close to experimental ionization potentials (IPs) in one-electron systems, demonstrating a surprisingly small effect of self-interaction and other exchange-correlation errors in established DFT methods. Remarkable agreement with ab initio quantum mechanical calculations of multielectron systems is demonstrated using several flavors of DFT, and we argue for the use of the average electron energy as a design criterion for density functional approximations.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Incorporation of the Fermi-Amaldi Term into Direct Energy Kohn-Sham Calculations

In direct energy Kohn-Sham (DEKS) theory, the density functional theory electronic energy equals the sum of occupied orbital energies, obtained from Kohn-Sham-like orbital equations involving a shifted Hartree exchange-correlation potential, which must be approximated. In the present study, the Fermi-Amaldi term is incorporated into approximate DEKS calculations, introducing the required -1/r contribution to the exchange-correlation component of the shifted potential in asymptotic regions. It also provides a mechanism for eliminating one-electron self-interaction error, and it introduces a nonzero exchange-correlation component of the shift in the potential that is of appropriate magnitude. The resulting electronic energies are very sensitive to the methodologies considered, whereas the highest occupied molecular orbital energies and exchange-correlation potentials are much less sensitive and are similar to those obtained from DEKS calculations using a conventional exchange-correlation functional.

Ensemble Density Functional Theory of Neutral and Charged Excitations : Exact Formulations, Standard Approximations, and Open Questions

Recent progress in the field of (time-independent) ensemble density-functional theory (DFT) for excited states are reviewed. Both Gross-Oliveira-Kohn (GOK) and N-centered ensemble formalisms, which are mathematically very similar and allow for an in-principle-exact description of neutral and charged electronic excitations, respectively, are discussed. Key exact results, for example, the equivalence between the infamous derivative discontinuity problem and the description of weight dependencies in the ensemble exchange-correlation density functional, are highlighted. The variational evaluation of orbital-dependent ensemble Hartree-exchange (Hx) energies is discussed in detail. We show in passing that state-averaging individual exact Hx energies can lead to severe (although solvable) v-representability issues. Finally, we explore the possibility of using the concept of density-driven correlation, which has been introduced recently and does not exist in regular ground-state DFT, for improving state-of-the-art correlation density-functional approximations for ensembles. The present review reflects the efforts of a growing community to turn ensemble DFT into a rigorous and reliable low-cost computational method for excited states. We hope that, in the near future, this contribution will stimulate new formal and practical developments in the field.

Relationships between Orbital Energies, Optical and Fundamental Gaps, and Exciton Shifts in Approximate Density Functional Theory and Quasiparticle Theory

The relationships between Kohn-Sham (KS) and generalized KS (GKS) density functional orbital energies and fundamental gaps or optical gaps raise many interesting questions including the physical meanings of KS and GKS orbital energies when computed with currently available approximate density functionals (ADFs). In this work, by examining three diverse databases with various ADFs, we examine such relations from the point of view of the exciton shift of quasiparticle theory. We start by calculating a large number of excitation energies by time-dependent density functional theory (TDDFT) with a large number of ADFs. To relate the exciton shift implicit in TDDFT to the exciton shift that is explicit in Green's function theory, we define the exciton shift in TDDFT as the difference of the response shift and the quasiparticle shift. We found a strong correlation between the response shift and the amount of Hartree-Fock exchange included in the density functional, with the response shift varying between -1 and 5 eV. This range is an order of magnitude larger than the mean errors of the TDDFT excitation energies. This result suggests that, with currently available functionals, the KS or GKS orbital energies should be treated as intermediate mathematical variables in the calculation of excitation energies rather than as the energies of independent-particle reference states for quasiparticle theory.

Screened range-separated hybrid by balancing the compact and slowly varying density regimes: Satisfaction of local density linear response

Due to their quantitative accuracy and ability to solve several difficulties, screened range-separated hybrid exchange-correlation functionals are now a standard approach for ab initio simulation of condensed matter systems. However, the screened range-separated hybrid functionals proposed so far are biased either toward compact or slowly varying densities. In this paper, we propose a screened range-separated hybrid functional, named HSEint, which can well describe these density regimes, achieving good accuracy for both molecular and solid-state systems. The semilocal part of the proposed functional is based on the PBEint generalized gradient approximation [E. Fabiano et al., Phys. Rev. B 82, 113104 (2010)], constructed for hybrid interfaces. To improve the functional performance, we employ exact or nearly exact constraints in the construction of range-separated hybrid functional, such as recovering of the local density linear response and semiclassical atom linear response.

Kinetic Correlation Functionals from the Entropic Regularization of the Strictly Correlated Electrons Problem

In this work, we study the entropic regularization of the strictly correlated electrons formalism, discussing the implications for density functional theory and establishing a link with earlier works on quantum kinetic energy and classical entropy. We carry out a very preliminary investigation (using simplified models) on the use of the solution of the entropic regularized problem to build approximations for the kinetic correlation functional at large coupling strengths. We also analyze lower and upper bounds to the Hohenberg-Kohn functional using the entropic regularized strictly correlated electrons problem.

Nine questions on energy decomposition analysis

The paper collects the answers of the authors to the following questions: Is the lack of precision in the definition of many chemical concepts one of the reasons for the coexistence of many partition schemes? Does the adoption of a given partition scheme imply a set of more precise definitions of the underlying chemical concepts? How can one use the results of a partition scheme to improve the clarity of definitions of concepts? Are partition schemes subject to scientific Darwinism? If so, what is the influence of a community's sociological pressure in the "natural selection" process? To what extent does/can/should investigated systems influence the choice of a particular partition scheme? Do we need more focused chemical validation of Energy Decomposition Analysis (EDA) methodology and descriptors/terms in general? Is there any interest in developing common benchmarks and test sets for cross-validation of methods? Is it possible to contemplate a unified partition scheme (let us call it the "standard model" of partitioning), that is proper for all applications in chemistry, in the foreseeable future or even in principle? In the end, science is about experiments and the real world. Can one, therefore, use any experiment or experimental data be used to favor one partition scheme over another? © 2019 Wiley Periodicals, Inc.

Strain Stabilization of Superionicity in Copper and Lithium Selenides

Superionic (SI) phases have utility as solid electrolytes for next generation battery technology, but these phases are typically not stable at room temperature. Our density functional theory calculations demonstrate that compressive lattice strain can stabilize SI phases of Cu₂Se and Li₂Se, two potential solid electrolytes. Electronic and bonding insights into this effect are obtained. In the ordered, non-SI phase, cations are localized primarily in tetrahedral (T) interstices with little access to the higher-energy octahedral (O) sites, but 1-2% compressive strain promotes attractive stabilization of the O cations with 6-fold coordination to Se anions, at the expense of the stability of 4-fold-coordinated T cations. In such compressed lattices, cations can access both T and O sites, resulting in a cation-disordered, SI phase. Thus, lattice strain is demonstrated as a handle for controlling ionic structure and transport and accomplishing ambient temperature superionicity.

Approximating the Shifted Hartree-Exchange-Correlation Potential in Direct Energy Kohn-Sham Theory

Levy and Zahariev [Phys. Rev. Lett. 113 113002 (2014)] have proposed a new approach for performing density functional theory calculations, termed direct energy Kohn-Sham (DEKS) theory. In this approach, the electronic energy equals the sum of orbital energies, obtained from Kohn-Sham-like orbital equations involving a shifted Hartree-exchange-correlation potential, which must be approximated. In the present study, density scaling homogeneity considerations are used to facilitate DEKS calculations on a series of atoms and molecules, leading to three nonlocal approximations to the shifted potential. The first two rely on preliminary Kohn-Sham calculations using a standard generalized gradient approximation (GGA) exchange-correlation functional and the results illustrate the benefit of describing the dominant Hartree component of the shift exactly. A uniform electron gas analysis is used to eliminate the need for these preliminary Kohn-Sham calculations, leading to a potential with an unconventional form that yields encouraging results, providing strong motivation for further research in DEKS theory.

Derivative discontinuity with localized Hartree-Fock potential

The localized Hartree-Fock potential has proven to be a computationally efficient alternative to the optimized effective potential, preserving the numerical accuracy of the latter and respecting the exact properties of being self-interaction free and having the correct -1/r asymptotics. In this paper we extend the localized Hartree-Fock potential to fractional particle numbers and observe that it yields derivative discontinuities in the energy as required by the exact theory. The discontinuities are numerically close to those of the computationally more demanding Hartree-Fock method. Our potential enjoys a "direct-energy" property, whereby the energy of the system is given by the sum of the single-particle eigenvalues multiplied by the corresponding occupation numbers. The discontinuities c↑ and c↓ of the spin-components of the potential at integer particle numbers N↑ and N↓ satisfy the condition c↑N↑ + c↓N↓ = 0. Thus, joining the family of effective potentials which support a derivative discontinuity, but being considerably easier to implement, the localized Hartree-Fock potential becomes a powerful tool in the broad area of applications in which the fundamental gap is an issue.