计算化学公社

标题: gamess能量分解时赝势出错 [打印本页]

作者
Author: 放飞心情 时间: 2017-3-2 10:59
标题: gamess能量分解时赝势出错
各位老师好，我在用gamess做能量分解时，需要对Au加赝势，把g以上的角动量砍掉了，计算出现下面的错误，是不是哪里没设置对，真心求教各位大神
ORTHOGONALIZING FOUND A REDUNDANT SET OF VECTORS ONLY  133 INDEPENDENT VECTORS FOUND
BUT THERE SHOULD BE  136 TOTAL MOS
EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT Thu Mar  2 10:36:52 2017

作者
Author: sobereva 时间: 2017-3-2 13:58
Au后面不该是19.，应当是79.
改完如果还不行，换基组再试。

作者
Author: 放飞心情 时间: 2017-3-2 17:03
谢谢sob老师，现在能算了。之前用的坐标是gaussian加了赝势生成的orb文件的坐标，当时把orb里的坐标和基组信息都拷到gamess里直接算，得到的相互作用能为正值，与gaussian的-54.56kcal/mol相差太大。考虑到gamess的赝势需要自己定义，我只把坐标保留了，基组自己去emsl上下的，没注意到前面坐标的问题。现在算出来的相互作用能为-50.57kcal/mol，跟gaussian的结果差了4kcal/mol

作者
Author: 放飞心情 时间: 2017-3-2 17:33

放飞心情发表于 2017-3-2 17:03
谢谢sob老师，现在能算了。之前用的坐标是gaussian加了赝势生成的orb文件的坐标，当时把orb里的坐标和基组 ...

我又把gaussian计算相互作用能时的Au的赝势按同样的截断方式处理了，发现对原来的能量没有影响，不知道gamess算出来的4kcal/mol的误差来源是什么。。

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