计算化学公社

标题: Gamess 赝势只能用于spdfg轨道？ —— cc-pvtz-pp基组报错 [打印本页]

作者
Author: brothers 时间: 2015-1-23 22:02
标题: Gamess 赝势只能用于spdfg轨道？ —— cc-pvtz-pp基组报错
本帖最后由 brothers 于 2015-1-24 00:11 编辑

因为在gaussian里计算用的是cc基组，所以希望gamess能量分解也在cc基组下完成。

但是在gamess中使用cc-pvtz基组，并为Os指定相应赝势，结果计算停止“THE ECP INTEGRALS ARE LIMITED TO SPDFG CORE POTENTIALS.”
这是因为cc-pvtz-pp的赝势部分包含有h轨道成分？为什么Gamess有这种错误？如何解决？

本问题在小木虫上问过，但没人回答。谢谢各位~~现在gamess的使用就卡在这里了......

作者
Author: sobereva 时间: 2015-1-23 23:43
这是因为cc-pVTZ-PP对于Os的赝势定义部分有h角动量，GAMESS-US目前对于ECP积分只能处理到g角动量。

非要用GAMESS-US结合这个基组的话，只能做个近似处理，把h角动量砍掉，此帖第5节有说明
赝势的函数形式以及在量子化学程序中定义的方式
http://sobereva.com/188

作者
Author: jjjspring 时间: 2015-1-24 00:15
如果这个赝势对大于等于h的角动量定义是零，那么砍掉h也没有影响了？
h-ul potential
1
2 1.00000000 0.00000000

作者
Author: brothers 时间: 2015-1-24 00:16

sobereva 发表于 2015-1-23 23:43
这是因为cc-pVTZ-PP对于Os的赝势定义部分有h角动量，GAMESS-US目前对于ECP积分只能处理到g角动量。

非要 ...

非常感谢！！
刚发现原帖内容有错，已经更正
“这是因为cc-pvtz-pp的赝势部分包含有g轨道成分？”更正为“这是因为cc-pvtz-pp的赝势部分包含有h轨道成分？”

麻烦您再帮我看看另一个帖子吧，最近纠结死了.....
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=670

作者
Author: brothers 时间: 2015-1-24 01:05
成功了。按照http://sobereva.com/188 的讲解，
1、将cc基组的赝势L设为4（仅包含spdfg轨道），并删除h轨道定义部分
2、将g轨道部分的定义剪切下来，放在最前面，像原来的h一样设为0(第二条是1.25日新修正的，最初并未设为0。为什么要设为0，可参考http://bbs.keinsci.com/forum.php ... a=page%3D1&mobile=2)
3、OK~~~
总结下来(并根据sob老师的讲解作了修改)以供以后的同学参考。
ps：
前些天意识到是h轨道的问题时，曾经找了L=4的赝势及相应基组作了计算，比如gamess input中的sbkjc基组，成功了。于是虎虎地就去把cc基组的pp部分的h部分直接删除了，算出来结果有问题....原来是角动量的顺序还有要求

~~~受教了~~非常感谢~~

作者
Author: brothers 时间: 2015-1-24 01:14

jjjspring 发表于 2015-1-24 00:15
如果这个赝势对大于等于h的角动量定义是零，那么砍掉h也没有影响了？
h-ul potential
1

应该是的吧，数学上为0，删掉理应没有影响~
gamess里算的话，可以找个L<=4并且最高角动量d值为0的试试，比如Cu的赝势就只有spdfg，其中g轨道d值为0.抽空算一个看看，然后再来补充。
回去睡觉了。

还有一个问题亟待解决~~各位大神走过路过不要错过

http://bbs.keinsci.com/forum.php?mod=viewthread&tid=670

作者
Author: sobereva 时间: 2015-1-24 02:36

jjjspring 发表于 2015-1-24 00:15
如果这个赝势对大于等于h的角动量定义是零，那么砍掉h也没有影响了？
h-ul potential
1

有影响。因为去掉h的话，就得让换g写到前头去，并定义为>=g为0。但是g原先不是0，所以等于g角动量部分没有感受到赝势。但由于角动量很高，造成的影响可能不大。

作者
Author: brothers 时间: 2015-1-24 13:15

sobereva 发表于 2015-1-24 02:36
有影响。因为去掉h的话，就得让换g写到前头去，并定义为>=g为0。但是g原先不是0，所以等于g角动量部分没 ...

sob老师的这条回复让我又产生了一个疑问，就是：
Stuttgart赝势难道必须定义 >=某个L的势为0吗？这有没有数学或物理上的要求？
我昨天的做法是将g的定义直接挪到最上面，也就是对>=g的部分全部使用和g相同的势（不为0），即下面这种形式，计算也能成功。
看了sob老师的回复感觉应该确认一下是否有问题~~谢谢老师~~
  PT-ECP GEN    60    4
  2    ----- g-ul potential    -----
-11.65174900    2    5.45202000
-14.37552500    2    5.41258500
  2    ----- s-ul potential    -----
429.64608700    2    14.60450000
   73.15688400    2    7.21828700
  4    ----- p-ul potential    -----
   88.02291700    2    11.57716200
175.99819600    2    10.88384300
   13.68227400    2    6.42440300
   27.41465100    2    5.22419800
  4    ----- d-ul potential    -----
   43.55785200    2    7.69961000
   65.36910800    2    7.55080800
   7.01859600    2    3.96116400
   11.39173300    2    3.87277700
  2    ----- f-ul potential    -----
   10.71022000    2    3.37986900
   14.27812500    2    3.32625500
  $END

作者
Author: brothers 时间: 2015-1-24 13:19

brothers 发表于 2015-1-24 13:15
sob老师的这条回复让我又产生了一个疑问，就是：
Stuttgart赝势难道必须定义 >=某个L的势为0吗？这有没 ...

赝势原始形式如下
$ECP
  PT-ECP GEN    60    5
  1    ----- h-ul potential    -----
   0.00000000    2    1.00000000
  2    ----- s-ul potential    -----
429.64608700    2    14.60450000
   73.15688400    2    7.21828700
  4    ----- p-ul potential    -----
   88.02291700    2    11.57716200
175.99819600    2    10.88384300
   13.68227400    2    6.42440300
   27.41465100    2    5.22419800
  4    ----- d-ul potential    -----
   43.55785200    2    7.69961000
   65.36910800    2    7.55080800
   7.01859600    2    3.96116400
   11.39173300    2    3.87277700
  2    ----- f-ul potential    -----
   10.71022000    2    3.37986900
   14.27812500    2    3.32625500
  2    ----- g-ul potential    -----
-11.65174900    2    5.45202000
-14.37552500    2    5.41258500
  $END

作者
Author: zjxitcc 时间: 2022-4-13 23:57
本帖最后由 zjxitcc 于 2022-4-13 23:59 编辑

发现GAMESS更新了：GAMESS 2021开始已经支持赝势H角动量，无需再做此类修改，直接算就行了。

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