g09D内置def2赝势与EMSL上下载的有差异

captain

请教大神，计算中涉及到了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
请问二者是否有真正的差别？还是只是输出显示的问题？
谢大神指导！

sobereva

等价的。因为在最下面有
F - G
                                                      2        3.1232691      -28.88430917    0.00000000

captain

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                     相比，正负号正好相反啊？

beefly

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

captain

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

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