Paula Mori-Sánchez

Revealing Noncovalent Interactions

Molecular structure does not easily identify the intricate noncovalent interactions that govern many areas of biology and chemistry, including design of new materials and drugs. We develop an approach to detect noncovalent interactions in real space, based on the electron density and its derivatives. Our approach reveals the underlying chemistry that compliments the covalent structure. It provides a rich representation of van der Waals interactions, hydrogen bonds, and steric repulsion in small molecules, molecular complexes, and solids. Most importantly, the method, requiring only knowledge of the atomic coordinates, is efficient and applicable to large systems, such as proteins or DNA. Across these applications, a view of nonbonded interactions emerges as continuous surfaces rather than close contacts between atom pairs, offering rich insight into the design of new and improved ligands.

Insights into Current Limitations of Density Functional Theory

Density functional theory of electronic structure is widely and successfully applied in simulations throughout engineering and sciences. However, for many predicted properties, there are spectacular failures that can be traced to the delocalization error and static correlation error of commonly used approximations. These errors can be characterized and understood through the perspective of fractional charges and fractional spins introduced recently. Reducing these errors will open new frontiers for applications of density functional theory.

Localization and Delocalization Errors in Density Functional Theory and Implications for Band-Gap Prediction

The band-gap problem and other systematic failures of approximate exchange-correlation functionals are explained from an analysis of total energy for fractional charges. The deviation from the correct intrinsic linear behavior in finite systems leads to delocalization and localization errors in large and bulk systems. Functionals whose energy is convex for fractional charges such as the local density approximation display an incorrect apparent linearity in the bulk limit, due to the delocalization error. Concave functionals also have an incorrect apparent linearity in the bulk calculation, due to the localization error and imposed symmetry. This resolves an apparent paradox and identifies the physical nature of the error to be addressed to obtain accurate band gaps from density functional theory.

Many-electron self-interaction error in approximate density functionals.

One of the most important challenges in density functional theory (DFT) is the proper description of fractional charge systems relating to the self-interaction error (SIE). Traditionally, the SIE has been formulated as a one-electron problem, which has been addressed in several recent functionals. However, these recent one-electron SIE-free functionals, while greatly improving the description of thermochemistry and reaction barriers in general, still exhibit many of the difficulties associated with SIE. Thus we emphasize the need to surpass this limit and shed light on the many-electron SIE. After identifying the sufficient condition for functionals to be free from SIE, we focus on the symptoms and investigate the performance of most popular functionals. We show that these functionals suffer from many-electron SIE. Finally, we give a SIE classification of density functionals.

Fractional charge perspective on the band gap in density-functional theory

The calculation of the band gap by density-functional theory (DFT) is examined by considering the behavior of the energy as a function of number of electrons. It is explained that the incorrect band-gap prediction with most approximate functionals originates mainly from errors in describing systems with fractional charges. Formulas for the energy derivatives with respect to number of electrons are derived, which clarify the role of optimized effective potentials in prediction of the band gap. Calculations with a recent functional that has much improved behavior for fractional charges give a good prediction of the energy gap and also

Development of exchange-correlation functionals with minimal many-electron self-interaction error

{\ensuremath{\epsilon}}_{\mathrm{HOMO}}\ensuremath{\simeq}\ensuremath{-}I

Self-interaction-free exchange-correlation functional for thermochemistry and kinetics

for finite systems. Our results indicate that it is possible, within DFT, to have a functional whose eigenvalues or derivatives accurately predict the band gap.

Fractional spins and static correlation error in density functional theory

New exchange-correlation functionals that address the important issue of many-electron self-interaction are developed. This is carried out by considering the performance of the functional on systems with fractional numbers of electrons at the same time as more standard thermochemical tests. The inclusion of Coulomb-attenuated exchange in the functional is facilitated by use of the adiabatic connection coupled with a short-range and long-range splittings. The new functionals have a good performance on thermochemistry and a much improved description of the total energy versus number of electrons and henceforth a much smaller many-electron self-interaction error.

Discontinuous Nature of the Exchange-Correlation Functional in Strongly Correlated Systems

We develop a self-interaction-free exchange-correlation functional which is very accurate for thermochemistry and kinetics. This is achieved by theoretical construction of the functional form and nonlinear fitting. We define a simple interpolation of the adiabatic connection that uses exact exchange, generalized gradient approximation (GGA) and meta-GGA functionals. The performance is optimized by fitting a small number of empirical parameters. Overall the new functional improves significantly upon hybrids and meta-GGAs while correctly describing one-electron systems. The mean absolute error on a large set of reaction barriers is reduced to 1.99 kcal/mol.

Delocalization errors in density functionals and implications for main-group thermochemistry

Electronic states with fractional spins arise in systems with large static correlation (strongly correlated systems). Such fractional-spin states are shown to be ensembles of degenerate ground states with normal spins. It is proven here that the energy of the exact functional for fractional-spin states is a constant, equal to the energy of the comprising degenerate pure-spin states. Dramatic deviations from this exact constancy condition exist with all approximate functionals, leading to large static correlation errors for strongly correlated systems, such as chemical bond dissociation and band structure of Mott insulators. This is demonstrated with numerical calculations for several molecular systems. Approximating the constancy behavior for fractional spins should be a major aim in functional constructions and should open the frontier for density functional theory to describe strongly correlated systems. The key results are also shown to apply in reduced density-matrix functional theory.