Effective local potentials for density and density-matrix functional approximations with non-negative screening density.

A way to improve the accuracy of the spectral properties in density functional theory (DFT) is to impose constraints on the effective, Kohn-Sham (KS), local potential [J. Chem. Phys. 136, 224109 (2012)]. As illustrated, a convenient variational quantity in that approach is the "screening" or "electron repulsion" density, ρrep, corresponding to the local, KS Hartree, exchange and correlation potential through Poisson's equation. Two constraints, applied to this minimization, largely remove self-interaction errors from the effective potential: (i) ρrep integrates to N - 1, where N is the number of electrons, and (ii) ρrep ≥ 0 everywhere. In this work, we introduce an effective "screening" amplitude, f, as the variational quantity, with the screening density being ρrep = f2. In this way, the positivity condition for ρrep is automatically satisfied, and the minimization problem becomes more efficient and robust. We apply this technique to molecular calculations, employing several approximations in DFT and in reduced density matrix functional theory. We find that the proposed development is an accurate, yet robust, variant of the constrained effective potential method.

Density inversion method for local basis sets without potential auxiliary functions: inverting densities from RDMFT.

A density inversion method is presented, to obtain the constrained, optimal, local potential that has a prescribed asymptotic behaviour and reproduces optimally any given ground-state electronic density. This work builds upon the method of [Callow et al., J. Chem. Phys., 2020, 152, 164114.] and differs in the expansion of the screening density in orbital basis element products instead of basis functions of an additional auxiliary set. We demonstrated the method by applying it to densities from DFT, Hartree-Fock, CAS-SCF and RDMFT calculations. For RDMFT, we demonstrate that density inversion offers a viable single-particle description by comparing the ionization potentials for atomic and molecular systems to the corresponding experimental values. Finally, we show that with the present method, accurate correlation potentials can be obtained from the inversion of accurate densities.

Density functionals with spin-density accuracy for open shells.

Electrons in zero external magnetic field can be studied with the Kohn-Sham (KS) scheme of either density functional theory (DFT) or spin-DFT (SDFT). The latter is normally used for open-shell systems because its approximations appear to model better the exchange and correlation (xc) functional, but also because, so far the application of DFT implied a closed-shell-like approximation. In the first part of this Communication, we show that correcting this error for open shells allows the approximate DFT xc functionals to become as accurate as those in SDFT. In the second part, we consider the behavior of SDFT for zero magnetic field. We show that the KS equations of SDFT emerge as the generalized KS equations of DFT in this limit, thus establishing a so far unknown link between the two theories.

Improving the exchange and correlation potential in density-functional approximations through constraints.

We review and expand on our work to impose constraints on the effective Kohn-Sham (KS) potential of local and semi-local density-functional approximations. Constraining the minimisation of the approximate total energy density-functional invariably leads to an optimised effective potential (OEP) equation, the solution of which yields the KS potential. We review briefly our previous work on this and demonstrate with numerous examples that despite the well-known mathematical issues of the OEP with finite basis sets, our OEP equations are numerically robust. We demonstrate that appropriately constraining the 'screening charge' which corresponds to the Hartree, exchange and correlation potential not only corrects its asymptotic behaviour but also allows the exchange and correlation potential to exhibit a non-zero derivative discontinuity, a feature of the exact KS potential that is necessary for the accurate prediction of band-gaps in solids but very hard to capture with semi-local approximations.

Density-inversion method for the Kohn-Sham potential: Role of the screening density.

We present a method to invert a given density and find the Kohn-Sham (KS) potential in Density Functional Theory (DFT) that shares the density. Our method employs the concept of screening density, which is naturally constrained by the inversion procedure and thus ensures that the density being inverted leads to a smooth KS potential with correct asymptotic behavior. We demonstrate the applicability of our method by inverting both local and non-local (Hartree-Fock and coupled cluster) densities; we also show how the method can be used to mitigate the effects of self-interactions in common DFT potentials with appropriate constraints on the screening density.

Orbitals from local RDMFT: Are they Kohn-Sham or natural orbitals?

Recently, an approximate theoretical framework was introduced, called local reduced density matrix functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1-RDM) are minimized under the additional condition that the optimal orbitals satisfy a single electron Schrödinger equation with a local potential. In the present work, we focus on the character of these optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained with the full minimization (without the extra condition) by contrasting them against the exact NOs and orbitals from a density functional calculation using the local density approximation (LDA). We find that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of the description of strong static correlation, this finding opens the way to a mixed density/density matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation at a cost only minimally higher than a density functional calculation.

A correction for the Hartree-Fock density of states for jellium without screening.

We revisit the Hartree-Fock (HF) calculation for the uniform electron gas, or jellium model, whose predictions--divergent derivative of the energy dispersion relation and vanishing density of states (DOS) at the Fermi level--are in qualitative disagreement with experimental evidence for simple metals. Currently, this qualitative failure is attributed to the lack of screening in the HF equations. Employing Slater's hyper-Hartree-Fock (HHF) equations, derived variationally, to study the ground state and the excited states of jellium, we find that the divergent derivative of the energy dispersion relation and the zero in the DOS are still present, but shifted from the Fermi wavevector and energy of jellium to the boundary between the set of variationally optimised and unoptimised HHF orbitals. The location of this boundary is not fixed, but it can be chosen to lie at arbitrarily high values of wavevector and energy, well clear from the Fermi level of jellium. We conclude that, rather than the lack of screening in the HF equations, the well-known qualitative failure of the ground-state HF approximation is an artifact of its nonlocal exchange operator. Other similar artifacts of the HF nonlocal exchange operator, not associated with the lack of electronic correlation, are known in the literature.

Quasi-particle energy spectra in local reduced density matrix functional theory.

Recently, we introduced [N. N. Lathiotakis, N. Helbig, A. Rubio, and N. I. Gidopoulos, Phys. Rev. A 90, 032511 (2014)] local reduced density matrix functional theory (local RDMFT), a theoretical scheme capable of incorporating static correlation effects in Kohn-Sham equations. Here, we apply local RDMFT to molecular systems of relatively large size, as a demonstration of its computational efficiency and its accuracy in predicting single-electron properties from the eigenvalue spectrum of the single-particle Hamiltonian with a local effective potential. We present encouraging results on the photoelectron spectrum of molecular systems and the relative stability of C20 isotopes. In addition, we propose a modelling of the fractional occupancies as functions of the orbital energies that further improves the efficiency of the method useful in applications to large systems and solids.

Electronic non-adiabatic states: towards a density functional theory beyond the Born-Oppenheimer approximation.

A novel treatment of non-adiabatic couplings is proposed. The derivation is based on a theorem by Hunter stating that the wave function of the complete system of electrons and nuclei can be written, without approximation, as a Born-Oppenheimer (BO)-type product of a nuclear wave function, X(R), and an electronic one, ΦR(r), which depends parametrically on the nuclear configuration R. From the variational principle, we deduce formally exact equations for ΦR(r) and X(R). The algebraic structure of the exact nuclear equation coincides with the corresponding one in the adiabatic approximation. The electronic equation, however, contains terms not appearing in the adiabatic case, which couple the electronic and the nuclear wave functions and account for the electron-nuclear correlation beyond the BO level. It is proposed that these terms can be incorporated using an optimized local effective potential.

Parametric representation of open quantum systems and cross-over from quantum to classical environment.

The behavior of most physical systems is affected by their natural surroundings. A quantum system with an environment is referred to as open, and its study varies according to the classical or quantum description adopted for the environment. We propose an approach to open quantum systems that allows us to follow the cross-over from quantum to classical environments; to achieve this, we devise an exact parametric representation of the principal system, based on generalized coherent states for the environment. The method is applied to the s = 1/2 Heisenberg star with frustration, where the quantum character of the environment varies with the couplings entering the Hamiltonian H. We find that when the star is in an eigenstate of H, the central spin behaves as if it were in an effective magnetic field, pointing in the direction set by the environmental coherent-state angle variables (θ, ϕ), and broadened according to their quantum probability distribution. Such distribution is independent of ϕ, whereas as a function of θ is seen to get narrower as the quantum character of the environment is reduced, collapsing into a Dirac-δ function in the classical limit. In such limit, because ϕ is left undetermined, the Von Neumann entropy of the central spin remains finite; in fact, it is equal to the entanglement of the original fully quantum model, a result that establishes a relation between this latter quantity and the Berry phase characterizing the dynamics of the central spin in the effective magnetic field.

Constraining density functional approximations to yield self-interaction free potentials.

Self-interactions (SIs) are a major problem in density functional approximations and the source of serious divergence from experimental results. Here, we propose to optimize density functional total energies in terms of the effective local potential, under constraints for the effective potential that guarantee it is free from SI errors and consequently asymptotically correct. More specifically, we constrain the Hartree, exchange and correlation potential to be the electrostatic potential of a non-negative effective repulsive density of N - 1 electrons. In this way, the optimal effective potentials exhibit the correct asymptotic decay, resulting in significantly improved one-electron properties.