Explanation of the Source of Very Large Errors in Many Exchange鈥揅orrelation Functionals for Vanadium Dimer

详细信息    查看全文

• 作者：Wenjing Zhang ; Donald G. Truhlar ; Mingsheng Tang
• 刊名：Journal of Chemical Theory and Computation
• 出版年：2014
• 出版时间：June 10, 2014
• 年：2014
• 卷：10
• 期：6
• 页码：2399-2409
• 全文大小：483K
• 年卷期：v.10,no.6(June 10, 2014)
• ISSN：1549-9626

文摘

Vanadium dimer is a notoriously difficult case for Kohn鈥揝ham (KS) density functional theory with currently available approximations to the exchange鈥揷orrelation (xc) functionals, and many approximate xc functionals yield an exceedingly large error in the calculated bond energy. In this paper, we first test the bond energies estimated by 43 xc functionals and the Hartree鈥揊ock (HF) method. The results further confirm the large errors and show that, with the experimental bond energy being 64.2 kcal/mol, the KS calculations give predictions all over the map with errors ranging from 鈭?1.5 to +60.5 kcal/mol, and the HF method performs much worse with an error of 鈭?24.4 kcal/mol! The reason for these very large errors is examined in this article by analyzing the atomic and molecular orbital energies calculated by various xc functionals. The results show that the errors in estimates of the bond energy of vanadium dimer can primarily be related to the calculated energy gap between the 4s and 3d_z2 atomic orbitals of the vanadium atom and especially to the 3d_z2 orbital energy. This interesting relation between the errors in the calculated bond energy and the magnitudes of the single-particle orbital energies provides a constructive alternative to the common but more sterile explanation that it is the static correlation energy due to multicenter left鈥搑ight correlation that makes the vanadium dimer and many other transition metal compounds so difficult for Kohn鈥揝ham calculations. One of the most important factors in determining the critical atomic orbital energy is the amount of nonlocal HF exchange that is included in the xc functional, but it is still difficult to explain why different local functionals (functionals with no HF exchange) yield quite different results. We conclude that improving calculations of orbital energies of atoms may provide a route to improving the accuracy of theoretical predictions of molecular bond energies for systems containing metal atoms.