Elimination of Spurious Fractional Charges in Dissociating Molecules by Correcting the Shape of Approximate Kohn–Sham Potentials

Outstanding improvement in removing the delocalization error by global natural orbital functional

This work assesses the performance of the recently proposed global natural orbital functional (GNOF) against the charge delocalization error. GNOF provides a good balance between static and dynamic electronic correlations leading to accurate total energies while preserving spin, even for systems with a highly multi-configurational character. Several analyses were applied to the functional, namely, (i) how the charge is distributed in super-systems of two fragments, (ii) the stability of ionization potentials while increasing the system size, and (iii) potential energy curves of a neutral and charged diatomic system. GNOF was found to practically eliminate the charge delocalization error in many of the studied systems or greatly improve the results obtained previously with PNOF7.

Non-empirical, low-cost recovery of exact conditions with model-Hamiltonian inspired expressions in jmDFT

Density functional theory (DFT) is widely applied to both molecules and materials, but well known energetic delocalization and static correlation errors in practical exchange-correlation approximations limit quantitative accuracy. Common methods that correct energetic delocalization errors, such as the Hubbard U correction in DFT+U or Hartree-Fock exchange in global hybrids, do so at the cost of worsening static correlation errors. We recently introduced an alternate approach [Bajaj et al., J. Chem. Phys. 147, 191101 (2017)] known as judiciously modified DFT (jmDFT), wherein the deviation from exact behavior of semilocal functionals over both fractional spin and charge, i.e., the so-called flat plane, was used to motivate functional forms of second order analytic corrections. In this work, we introduce fully nonempirical expressions for all four coefficients in a DFT+U+J-inspired form of jmDFT, where all coefficients are obtained only from energies and eigenvalues of the integer-electron systems. We show good agreement for U and J coefficients obtained nonempirically as compared with the results of numerical fitting in a jmDFT U+J/J' correction. Incorporating the fully nonempirical jmDFT correction reduces and even eliminates the fractional spin error at the same time as eliminating the energetic delocalization error. We show that this approach extends beyond s-electron systems to higher angular momentum cases including p- and d-electrons. Finally, we diagnose some shortcomings of the current jmDFT approach that limit its ability to improve upon DFT results for cases such as weakly bound anions due to poor underlying semilocal functional behavior.

Visualizing atomic sizes and molecular shapes with the classical turning surface of the Kohn-Sham potential

The Kohn-Sham potential [Formula: see text] is the effective multiplicative operator in a noninteracting Schrödinger equation that reproduces the ground-state density of a real (interacting) system. The sizes and shapes of atoms, molecules, and solids can be defined in terms of Kohn-Sham potentials in a nonarbitrary way that accords with chemical intuition and can be implemented efficiently, permitting a natural pictorial representation for chemistry and condensed-matter physics. Let [Formula: see text] be the maximum occupied orbital energy of the noninteracting electrons. Then the equation [Formula: see text] defines the surface at which classical electrons with energy [Formula: see text] would be turned back and thus determines the surface of any electronic object. Atomic and ionic radii defined in this manner agree well with empirical estimates, show regular chemical trends, and allow one to identify the type of chemical bonding between two given atoms by comparing the actual internuclear distance to the sum of atomic radii. The molecular surfaces can be fused (for a covalent bond), seamed (ionic bond), necked (hydrogen bond), or divided (van der Waals bond). This contribution extends the pioneering work of Z.-Z. Yang et al. [Yang ZZ, Davidson ER (1997) Int J Quantum Chem 62:47-53; Zhao DX, et al. (2018) Mol Phys 116:969-977] by our consideration of the Kohn-Sham potential, protomolecules, doubly negative atomic ions, a bond-type parameter, seamed and necked molecular surfaces, and a more extensive table of atomic and ionic radii that are fully consistent with expected periodic trends.

How Interatomic Steps in the Exact Kohn-Sham Potential Relate to Derivative Discontinuities of the Energy

Accurate density functional calculations hinge on reliable approximations to the unknown exchange-correlation (xc) potential. The most popular approximations usually lack features of the exact xc potential that are important for an accurate prediction of the fundamental gap and the distribution of charge in complex systems. Two principal features in this regard are the spatially uniform shift in the potential, as the number of electrons infinitesimally surpasses an integer, and the spatial steps that form, for example, between the atoms of stretched molecules. Although both aforementioned concepts are well known, the exact relationship between them remained unclear. Here we establish this relationship via an analytical derivation. We support our result by numerically solving the many-electron Schrödinger equation to extract the exact Kohn-Sham potential and directly observe its features. Spatial steps in the exact xc potential of a full configuration-interaction (FCI) calculation of a molecule are presented in three dimensions.