How Do Self-Interaction Errors Associated with Stretched Bonds Affect Barrier Height Predictions?

Abstract

Density functional theory (DFT) suffers from selfinteraction errors (SIEs) that generally result in the underestimation of chemical reaction barrier heights. This is commonly attributed to the tendency of density functional approximations to overstabilize delocalized densities that typically occur in the stretched bonds of transition state structures. The Perdew−Zunger self-interaction correction (PZSIC) and locally scaled selfinteraction correction (LSIC) improve the prediction of barrier heights of chemical reactions, with LSIC giving better accuracy than PZSIC on average. These methods employ an orbital-byorbital correction scheme to remove the one-electron SIE. In the context of barrier heights, this allows an analysis of how the selfinteraction correction (SIC) for each orbital contributes to the calculated barriers using Fermi−Loẅdin orbitals (FLOs). We hypothesize that the SIC contribution to the reaction barrier comes mainly from a limited number of orbitals that are directly involved in bond-breaking and bond-making in the reaction transition state. We call these participant orbitals (POs), in contrast to spectator orbitals (SOs) which are not directly involved in changes to the bonding. Our hypothesis is that ΔE Total SIC ≈ ΔE PO SIC , where ΔE Total SIC is the difference in the SIC corrections for the reactants or products and the transition state. We test this hypothesis for the reaction barriers of the BH76 benchmark set of reactions. We find that the stretched-bond orbitals indeed make the largest individual SIC contributions to the barriers. These contributions increase the barrier heights relative to LSDA, which underpredicts the barrier. However, the full stretched-bond hypothesis does not hold in all cases for either PZSIC or LSIC. There are many cases where the total SIC contribution from the SOs is significant and cannot be ignored. The size of the SIC contribution to the barrier height is a key indicator. A large SIC correction is correlated to a large LSDA error in the barrier, showing that PZSIC properly gives larger corrections when corrections are needed most. A comparison of the performance of PZSIC and LSIC shows that the two methods have similar accuracy for reactions with large LSDA errors, but LSIC is clearly better for reactions with small errors. We trace this to an improved description of reaction energies in LSIC.

How Does HF-DFT Achieve Chemical Accuracy for Water Clusters?

“…Most of the BH76 transition states are stretched radicals, and none is like stretched H 2 , so the large negative FEs of the BH76 barrier heights are qualitatively explained. An independent analysis 56 of the BH76 reactions found (as did ref 7) that errors made by DFAs were in large part due to delocalized orbitals in the stretched bonds of transition states. Residual errors due to the less-chemically active orbitals were also significant in certain cases.…”

Analytic Gradients for the Electrostatic Embedding QM/MM Model in Periodic Boundary Conditions Using Particle-Mesh Ewald Sums and Electrostatic Potential Fitted Charge Operators

“…The profile computed using DFT, instead, recovers the relative energies of the cis and trans isomers correctly but largely underestimates the height of the barrier. This behavior can be explained by the tendency of DFT, due to self-interaction error, to overstabilize delocalized densities that usually characterize the distorted geometries of transition-state structures. , …”

“…This behavior can be explained by the tendency of DFT, due to selfinteraction error, to overstabilize delocalized densities that usually characterize the distorted geometries of transition-state structures. 116,117 4.2.2. AAPF Peptide Isomerization in Water.…”

Unconventional Error Cancellation Explains the Success of Hartree–Fock Density Functional Theory for Barrier Heights

“…Clearly, improved functionals will need the right amount of fully nonlocal density dependence in both the exchange-correlation energy and the exchange-correlation potential. Self-interaction corrections , to DFAs, while needing improvement for some properties, appear for barrier heights to get the right answer for the right reason, by significantly reducing both functional- and density-driven errors.…”