sobereva 发表于 2022-4-30 07:40
对于G16，1和2等价。不要用3

B2PLYP的D3参数问题详见此文

谢谢老师！

对于G16，1和2等价。不要用3

B2PLYP的D3参数问题详见此文
谈谈不同量子化学程序计算结果的差异问题
http://sobereva.com/573（http://bbs.keinsci.com/thread-20122-1-1.html）

要注意不同程序，哪怕相同的程序的不同版本，用的D3参数也有可能不同。参数比较混乱的典型是B2PLYP，这个泛函总共有4套D3参数：
(1)在Grimme的2010年DFT-D3的原文的表4中，B2PLYP的DFT-D3零阻尼参数是S6=0.5，Sr,6=1.332，S8=1.000
(2)在Grimme的PWPB95-D3的原文J. Chem. Theor. Comput., 7, 291 (2011)中，作者重新拟合了B2PLYP的DFT-D3零阻尼参数，由表3可见参数为S6=0.64，Sr,6=1.427，S8=1.022
(3)在Grimme的专门对比零阻尼和BJ阻尼形式的文章J. Comput. Chem., 32, 1456 (2011)的表2中，B2PLYP的DFT-D3(BJ)参数为S6=0.5，S8=1.0860，a1=0.3451，a2=4.7735
(4)在Grimme的GMTKN30测试集文章Phys. Chem. Chem. Phys.,13, 6670 (2011)的表S1中，B2PLYP的DFT-D3(BJ)参数为S6=0.64，S8=0.9147，a1=0.3065，a2=5.0570
可见哪怕同一种阻尼函数都有两套参数，发表年份还这么接近，不导致不同程序出现分歧才怪，所以互相对比时要看清楚输出信息中有无提示具体参数、看看手册里的参数引的哪篇文章。在Gaussian中，B2PLYP加DFT-D3(BJ)校正有两种写法，一个是B2PLYPD3，一个是B2PLYP em=GD3BJ，然而在G09中这两种写法的色散校正值结果不同，估计就是因为用的D3(BJ)参数不同所致，但手册里却没提。而在G16中这两种写法的结果则变得统一了，都等于G09中B2PLYPD3的写法，此时用的参数对应的是上面的(4)。