Search

Standard approximations for the exchange-correlation (XC) functional in Kohn-Sham density functional theory (KS-DFT) typically lead to unacceptably large errors when applied to strongly-correlated electronic systems. Partition-DFT (PDFT) is a formally exact reformulation of KS-DFT in which the ground-state density and energy of a system are obtained through self-consistent calculations on isolated fragments, with a partition energy representing the \textit{inter}-fragment interactions. Here we show how typical errors of the local density approximation (LDA) in KS-DFT can be largely suppressed through a simple approximation, the generalized overlap approximation (GOA), for the partition energy in PDFT. Our method is illustrated on simple models of one-dimensional strongly-correlated linear hydrogen chains. The GOA, when used in combination with the LDA for the fragments, improves the LDA dissociation curves of hydrogen chains and produces results that are comparable to those of spin-unrestricted LDA, but without breaking the spin symmetry. GOA also induces a correction to the LDA electron density that partially captures the correct density dimerization in strongly-correlated hydrogen chains. Moreover, with an additional correction to the partition energy, the approximation is shown to produce dissociation energies in quantitative agreement to calculations based on the Density Matrix Renormalization Group method.

Accurate first-principles calculations for the energies, charge distributions, and spin symmetries of many-electron systems are essential to understand and predict the electronic and structural properties of molecules and materials. Kohn-Sham density functional theory (KS-DFT) stands out among electronic-structure methods due to its balance of accuracy and computational efficiency. It is now extensively used in fields ranging from materials engineering to rational drug design. However, to achieve chemically accurate energies, standard density functional approximations in KS-DFT often need to break underlying symmetries, a long-standing "symmetry dilemma". By employing fragment spin densities as the main variables in calculations (rather than total molecular densities as in KS-DFT), we present an embedding framework in which this symmetry dilemma is resolved for the case of stretched molecules. The spatial overlap between fragment densities is used as the main ingredient to construct a simple, physically-motivated approximation to a universal functional of the fragment densities. This 'overlap approximation' is shown to significantly improve semi-local KS-DFT binding energies of molecules without artificial symmetry breaking.

We explore a variety of unsolved problems in density functional theory, where mathematicians might prove useful. We give the background and context of the different problems, and why progress toward resolving them would help those doing computations using density functional theory. Subjects covered include the magnitude of the kinetic energy in Hartree-Fock calculations, the shape of adiabatic connection curves, using the constrained search with input densities, densities of states, the semiclassical expansion of energies, the tightness of Lieb-Oxford bounds, and how we decide the accuracy of an approximate density., Comment: 21 pages, 6 figures

Partition Density Functional Theory (P-DFT) is a density embedding method that partitions a molecule into fragments by minimizing the sum of fragment energies subject to a local density constraint and a global electron-number constraint. To perform this minimization, we study a two-stage procedure in which the sum of fragment energies is lowered when electrons flow from fragments of lower electronegativity to fragments of higher electronegativity. The global minimum is reached when all electronegativities are equal. The non-integral fragment populations are dealt with in two different ways: (1) An ensemble approach (ENS) that involves averaging over calculations with different numbers of electrons (always integers); and (2) A simpler approach that involves fractionally occupying orbitals (FOO). We compare and contrast these two approaches and examine their performance in some of the simplest systems where one can transparently apply both, including simple models of heteronuclear diatomic molecules and actual diatomic molecules with 2 and 4 electrons. We find that, although both ENS and FOO methods lead to the same total energy and density, the ENS fragment densities are less distorted than those of FOO when compared to their isolated counterparts, and they tend to retain integer numbers of electrons. We establish the conditions under which the ENS populations can become fractional and observe that, even in those cases, the total charge transferred is always lower in ENS than in FOO. Similarly, the FOO fragment dipole moments provide an upper bound to the ENS dipoles. We explain why, and discuss implications., Comment: 11 pages, 12 figures