用GGA密度泛函及其长程、色散校正方法计算各类氢键的结合能

应用广义梯度近似(GGA) (PW91和PBE)、含动能密度的广义梯度近似(meta-GGA) (M06-L)、杂化泛函(hyper-GGA)(M06-2X、X3LYP和B3LYP)及其长程校正泛函LC-DFT(CAM-B3LYP、LC-ωPBE和ωB97X)和色散校正密度泛函(DFT-D)(ωB97X-D和B97-D),用多种基函数对15种不同强度的传统氢键和非传统氢键体系的结合能进行了系统的计算与分析.并与高精度的CCSD(T)/aug-cc-pVQZ结果比较发现:在上述各类泛函中,对于氢键结合能的计算M06-2X和ωB97X-D泛函较为精确与可靠,且没有必要使用过大的基函数,6-311++G(2d,2p)或aug-cc-pVDZ水平的基组就已足够,各类泛函所计算结合能的基组重叠误差(BSSE)均较小,除ωB97X和ωB97X-D外,其它9种泛函不经BSSE校正也能得到同样甚至更准确的结果.

Calculation of the Binding Energies of Different Types of Hydrogen Bonds Using GGA Density Functional and Its Long-Range, Empirical Dispersion Correction Methods

We investigated eleven exchange-correlation energy density functionals including generalized gradient approximation (GGA) (PBE, PW91), meta-GGA (M06-L), hyper-GGA (M06-2X, B3LYP, X3LYP), LC-DFT methods (CAM-B3LYP, LC-ωPBE, ωB97X), and density functional theory with dispersion corrections (DFT-D) methods (ωB97X-D, B97-D) for their performance in describing systems with conventional and non-conventional hydrogen bonds. After comparing the results using the benchmark CCSD(T)/aug-cc-pVQZ approach we found that the M06-2X and ωB97X-D functionals provided the most accurate and reliable results for the fifteen systems studied in this work with strong, moderate, and weak hydrogen bonds. It is important to employ an appropriate basis set to predict the binding energy of hydrogen bonds for all DFT methods and we found that the basis set of 6-311++G(2d, 2p) or aug-cc-pVDZ is adequate. The effect of the basis set superposition error (BSSE) is relatively small for the DFT methods tested. All the methods except for ωB97X and ωB97X-D were found to produce equally accurate or even more accurate results without BSSE correction.