Unconstrained Variational Determination of the Kohn-Sham Potential

Density Functional Theory's Kohn-Sham (KS) potential emerges as the minimizing effective potential in an unconstrained variational scheme that does not involve fixing the unknown single-electron density. The physical content behind the virtual KS system, that of a non-interacting electronic system representing optimally the interacting one, is brought to light

Local reduced-density-matrix-functional theory: Incorporating static correlation effects in Kohn-Sham equations

We propose a novel scheme to bring reduced density matrix functional theory (RDMFT) into the realm of density functional theory (DFT) that preserves the accurate density functional description at equilibrium, while incorporating accurately static and left-right correlation effects in molecules and keeping the good computational performance of DFT-based schemes. The key ingredient is to relax the requirement that the local potential is the functional derivative of the energy with respect to the density. Instead, we propose to restrict the search for the approximate natural orbitals within a domain where these orbitals are eigenfunctions of a single-particle hamiltonian with a local effective potential. In this way, fractional natural occupation numbers are accommodated into Kohn-Sham equations allowing for the description of molecular dissociation without breaking spin symmetry. Additionally, our scheme provides a natural way to connect an energy eigenvalue spectrum to the approximate natural orbitals and this spectrum is found to represent accurately the ionization potentials of atoms and small molecules

Nonuniqueness in spin-density-functional theory on lattices

In electronic many-particle systems, the mapping between densities and spin magnetizations, {n(r), m(r)}, and potentials and magnetic fields, {v(r), B(r)}, is known to be nonunique, which has fundamental and practical implications for spin-density-functional theory (SDFT). This paper studies the nonuniqueness (NU) in SDFT on arbitrary lattices. Two new, non-trivial cases are discovered, here called local saturation and global noncollinear NU, and their properties are discussed and illustrated. In the continuum limit, only some well-known special cases of NU survive.Comment: 4 pages, 1 figur

Local reduced-density-matrix-functional theory: Incorporating static correlation effects in Kohn-Sham equations

We propose a scheme to bring reduced-density-matrix-functional theory into the realm of density functional theory (DFT) that preserves the accurate density functional description at equilibrium, while incorporating accurately static and left-right correlation effects in molecules and keeping the good computational performance of DFT-based schemes. The key ingredient is to relax the requirement that the local potential is the functional derivative of the energy with respect to the density. Instead, we propose to restrict the search for the approximate natural orbitals within a domain where these orbitals are eigenfunctions of a single-particle Hamiltonian with a local effective potential. In this way, fractional natural occupation numbers are accommodated into Kohn-Sham equations allowing for the description of molecular dissociation without breaking spin symmetry. Additionally, our scheme provides a natural way to connect an energy eigenvalue spectrum to the approximate natural orbitals and this spectrum is found to represent accurately the ionization potentials of atoms and small molecules.N.N.L. acknowledges financial support from the GSRT, Greece, Polynano-Kripis project (Grant No. 447963); N.H. from an Emmy-Noether grant from Deutsche Forschungsgemeinschaft; and A.R. from the European Community’s FP7 through the CRONOS project, grant agreement no. 280879, the European Research Council Advanced Grant DYNamo (ERC-2010-AdG-267374), and Grupos Consolidados UPV/EHU del Gobierno Vasco (Grant: IT578-13).Peer Reviewe

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. A90, 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 solidsN.N.L. acknowledges financial support from the GSRT action KPHΠIΣ, project “New multifunctional Nanostructured Materials and Devices – POLYNANO” No. 447963, N.H. from a DFG Emmy-Noether grant, and A.R. from the European Research Council Advanced Grant No. ERC-2010-AdG-267374, Spanish Grant No. FIS2010-21282-C02-01, Grupo Consolidado UPV/EHU (IT578-13), and European Commission Project No. CRONOS(280879-2).Peer Reviewe

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

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