计算化学公社

标题: 高斯中B2PLYP色散校正的写法疑问 [打印本页]
作者
Author: WaterMirror    时间: 2022-4-30 07:15
标题: 高斯中B2PLYP色散校正的写法疑问
尊敬的各位老师:

   高斯的说明书中指出: B2PLYPD3这个关键词代表B2PLYP双杂化泛函加D3(BJ)校正。请问各位老师下面的三种写法（使用高斯16A03版本）:
   1. # opt freq b2plypd3 jun-cc-pvtz
   2. # opt freq b2plyp jun-cc-pvtz em=gd3bj
   3. # opt freq b2plyp jun-cc-pvtz em=gd3

  是否1和2完全等价呢？我自己实测了几个例子 发现1和2输出的能量都是一样的，和3的输出能量不同，但是仍然有些不放心....

  以及，高斯内置了b2plypd3这个关键词之后，对于b2plyp双杂化泛函，还能用em=gd3,em=gd3(bj)的方式调用D3,D3(BJ)校正吗？对于3这种写法，D3(zero)校正有没有加上去呢.....?

  谢谢各位老师！



作者
Author: sobereva    时间: 2022-4-30 07:40
对于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)。
作者
Author: WaterMirror    时间: 2022-4-30 07:47
sobereva 发表于 2022-4-30 07:40
对于G16，1和2等价。不要用3

B2PLYP的D3参数问题详见此文

谢谢老师！



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3