Abstract:We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn-Sham Density Functional Theory functionals on simulation of small to medium-size proteins in vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO-LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the Kohn-Sham (KS) self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver [Phys. Chem. Chem. Phys. 26, 6557 (2024)] to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree-Fock and KS solutions; the charge delocalization can occur not only between charged fragments (such as in zwitterionic polypeptides) but also may involve neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of {\em some} range-separated hybrid functionals, but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping "frontier" orbitals.

What problem does this paper attempt to address?

This paper mainly investigates the nonphysical charge delocalization/fractionalization problem that arises in simulating small to medium-sized proteins using the Kohn-Sham density functional theory (KS DFT). This problem leads to self-interaction errors, resulting in incorrect electron density distribution, distorted total energy values, and a decrease or disappearance of the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), thereby causing difficulties in solving the Kohn-Sham self-consistent field (SCF) problem. The authors use a robust local convergence SCF solver to address these issues and reveal the detailed structure of charge delocalization through natural deformation orbital analysis. The study found that charge delocalization occurs not only between charged fragments, but also potentially involves neutral fragments. While certain range-separated hybrid functionals can partially alleviate this problem, they may not effectively solve it in all cases. The paper also points out that even for semi-local KS functionals prone to such defects, charge-localized solutions without delocalization can be found, but these solutions exhibit non-Aufbau characteristics and instability in the frontier orbitals. In the paper, the authors demonstrate how local solvers can find solutions with properly positioned charges in some cases through the analysis of 17 peptides, even though these solutions are not energy-minimized. Additionally, they scan combinations of 20 natural amino acids with various functionals to determine which ones are likely to trigger such issues. The research results indicate that range-separated functionals are helpful to a large extent, but cannot completely eliminate the nonphysical charge delocalization and SCF solving problems. Finally, the paper emphasizes the need to exercise caution and suggests further exploration of solutions and the use of more advanced KS functionals when encountering unexpected small or zero HOMO-LUMO gaps and abnormal SCF convergence patterns in Kohn-Sham DFT simulations.

Revisiting Artifacts of Kohn-Sham Density Functionals for Biosimulation

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

Delocalization error poisons the density-functional many-body expansion

Symmetry Dilemma of Doubly Hybrid Density Functionals for Equilibrium Molecular Property Calculations.

Stretching bonds in Density Functional Theory without artificial symmetry breaking

How Semilocal Are Semilocal Density Functional Approximations? -Tackling Self-Interaction Error in One-Electron Systems

Range separation and local hybridization in density functional theory

Scaling Correction Approaches for Reducing Delocalization Error in Density Functional Approximations

Subsystem Density Functional Theory Augmented by a Delta Learning Approach to Achieve Kohn–Sham Accuracy

Multiconfiguration Pair-Density Functional Theory is Free from Delocalization Error

Orbital Anatomy of Self-Interaction in Kohn-Sham Density Functional Theory

Many-electron self-interaction and spin polarization errors in local hybrid density functionals

The role of electron localization in density functionals

Requirements for very high temperature Kohn-Sham density functional simulations and how to bypass them

On the Analysis of the Discretized Kohn-Sham Density Functional Theory

Density Functionals with Quantum Chemical Accuracy: From Machine Learning to Molecular Dynamics

Consistent analytical second derivatives of the Kohn–Sham DFT energy in the framework of the Conductor-like screening Model (COSMO) through Gaussian charge distributions

Understanding density driven errors for reaction barrier heights

Delocalization Error Of Density-Functional Approximations: A Distinct Manifestation In Hydrogen Molecular Chains