计算化学公社

标题: g09D内置def2赝势与EMSL上下载的有差异 [打印本页]

作者
Author: captain 时间: 2020-3-13 10:19
标题: g09D内置def2赝势与EMSL上下载的有差异
请教大神，计算中涉及到了Cs元素，用了def2基组及赝势，
从EMSL上下载的赝势如下
CS    0
CS-ECP    3    46
f potential
  1
2    3.1232690          -28.8843090
s-f potential
  3
2    4.0797500          84.5477300
2    2.4174060          16.6541730
2    3.1232690          28.8843090
p-f potential
  3
2    5.5140800          157.0490590
2    2.1603160          26.4233070
2    3.1232690          28.8843090
d-f potential
  3
2    1.8074100          13.1727530
2    0.8581820             3.3428330
2    3.1232690          28.8843090

直接在g09中使用其内置的def2基组（def2TZVP GFPrint），输出文件中显示如下
                                    G and up
                                                   2       1.0000000       0.00000000 0.00000000
                                    S - G
                                                   2       4.0797507    84.54773062 0.00000000
                                                   2       2.4174064    16.65417395 0.00000000
                                    P - G
                                                   2       5.5140802    157.04905939 0.00000000
                                                   2       2.1603168    26.42330704 0.00000000
                                    D - G
                                                   2       1.8074100    13.17275381 0.00000000
                                                   2       0.8581822       3.34283398 0.00000000
                                    F - G
                                                   2       3.1232691    -28.88430917 0.00000000

发现s、p、d对>=L势的修正都少了这么一项：
2    3.1232690          28.8843090
请问二者是否有真正的差别？还是只是输出显示的问题？
谢大神指导！

作者
Author: sobereva 时间: 2020-3-13 10:53
等价的。因为在最下面有
F - G
2 3.1232691 -28.88430917 0.00000000

作者
Author: captain 时间: 2020-3-13 11:03

sobereva 发表于 2020-3-13 10:53
等价的。因为在最下面有
F - G
2 3.123 ...

大神，
这一项 2 3.1232691 -28.88430917 0.00000000 不是 f对>=L势的修正吗？
而且与这一项 2 3.1232690 28.8843090 相比，正负号正好相反啊？

作者
Author: beefly 时间: 2020-3-13 15:04
def2是turbomole作者开发的，赝势来自sdd，但是当时的turbomole不支持g、h赝势函数，所以就降到了f，同时为了补偿带来的误差，额外加入了s-f、p-f、d-f投影项。他们自称造成的误差很小，但是似乎没人对比过。所以像gaussian等程序为了避免可能存在的问题，仍然坚持用原始的sdd赝势。

作者
Author: captain 时间: 2020-3-13 15:54
本帖最后由 captain 于 2020-3-13 16:04 编辑

beefly 发表于 2020-3-13 15:04
def2是turbomole作者开发的，赝势来自sdd，但是当时的turbomole不支持g、h赝势函数，所以就降到了f，同时为 ...

明白了，谢谢指导！
另外，对于这个例子，二者都是支持到 f 吧，g09中 g 及其以上（2 1.0000000 0.00000000 0.00000000）是没有作用的。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)