计算化学公社

标题: Gaussian 与Gamess 的能量相差10倍 [打印本页]

作者
Author: xiaobogaga 时间: 2015-9-29 09:12
标题: Gaussian 与Gamess 的能量相差10倍
大家好，我是用Gaussian09进行优化，然后用Gamess进行能量分解，但是却发现高斯的相互作用能，要比Gamess的高10倍，我用的是Mp2方法，aug-cc-pVDZ基组，构型优化和能量分解都没有问题，但是就是能量二者相差10倍，不知道是什么原因，希望大家能给予解答，万分感谢！！

作者
Author: yflchx 时间: 2015-9-29 09:39
检查：1.计算能量的分子构型，单体重新优化了？抑或用双体里相应的构型？
2.BSSE
3.full和frozencore是否一致
4.自洽场收敛情况

作者
Author: xiaobogaga 时间: 2015-9-29 09:44

yflchx 发表于 2015-9-29 09:39
检查：1.计算能量的分子构型，单体重新优化了？抑或用双体里相应的构型？
2.BSSE
3.full和frozencore是否 ...

1.单体的能量都是用二体里的构型计算的；
2.在优化二体的时候，都加了counter=2；
3.因为是新手，不知道怎么看。。。
4.单体、二体优化的时候都是正常结束，没有虚频

作者
Author: yflchx 时间: 2015-9-29 09:52

xiaobogaga 发表于 2015-9-29 09:44
1.单体的能量都是用二体里的构型计算的；
2.在优化二体的时候，都加了counter=2；
3.因为是新手，不知 ...

3.是指是不是计算能量的时候一个软件用了MP2(fc)，另一个软件用了MP2(full)？

4.看看SCF自洽场能量迭代的情况，对比一下两个程序。

作者
Author: xiaobogaga 时间: 2015-9-29 17:01

yflchx 发表于 2015-9-29 09:52
3.是指是不是计算能量的时候一个软件用了MP2(fc)，另一个软件用了MP2(full)？

4.看看SCF自洽场能量迭 ...

3. 计算单体和二体的时候命令行只写了 MP2/基组怎么看是(fc) 还是(full) 呢

作者
Author: yflchx 时间: 2015-9-29 17:55

xiaobogaga 发表于 2015-9-29 17:01
3. 计算单体和二体的时候命令行只写了 MP2/基组怎么看是(fc) 还是(full) 呢

默认是FC，frozencore。这个问题可以跳过了。

作者
Author: sobereva 时间: 2015-9-29 23:28
碰见问题先把问题尽可能简化，一点一点对比来分析原因。
比如考察相互作用能的计算差异，你就先分析单点能的计算差异，别把其它乱七八糟的因素引入，否则肯定弄不清楚。
你就用相同的结构，在MP2/aug-cc-pVDZ下用高斯和GAMESS-US各算个单点。注意都先别用冻核（两个程序在这点上有差异），看看HF能量和MP2相关能是否都能对上。假设没对上，然后你再看比如是否是实际计算用的基函数数目不对（比如高斯会自动去除线性依赖的基函数）。
慢慢一点一点对比总能搞清楚。

作者
Author: xiaobogaga 时间: 2015-9-30 13:52

sobereva 发表于 2015-9-29 23:28
碰见问题先把问题尽可能简化，一点一点对比来分析原因。
比如考察相互作用能的计算差异，你就先分析单点能 ...

非常感谢sob老师的耐心解答和指导，您的话也为我以后的科研指点了方向，我再试试，谢谢您

作者
Author: xiaobogaga 时间: 2015-9-30 16:46

sobereva 发表于 2015-9-29 23:28
碰见问题先把问题尽可能简化，一点一点对比来分析原因。
比如考察相互作用能的计算差异，你就先分析单点能 ...

老师你好，我用Gassian和Gamess 分别计算了能量为什么会差这么多，附件中是输入和输出文件

作者
Author: jjjspring 时间: 2015-10-1 02:44
基组输入错误。两个计算的基函数数目不同。

作者
Author: sobereva 时间: 2015-10-1 04:25

xiaobogaga 发表于 2015-9-30 16:46
老师你好，我用Gassian和Gamess 分别计算了能量为什么会差这么多，附件中是输入和输出文件

完全没有可比性嘛
1 一个是单点一个是counterpoise任务
2 一个对配体用cc-pVDZ一个用aug-cc-pVDZ
3 一个用冻核一个没用冻核

作者
Author: xiaobogaga 时间: 2015-10-1 18:20

jjjspring 发表于 2015-10-1 02:44
基组输入错误。两个计算的基函数数目不同。

你好，我是新手，请问正确的输入应该怎么写，比如我用Gamess 计算的时候用MP2(full)，在Gamess的输入文件里应该怎么写关键词呢？还有在Gamess里accd 表示 aug-cc-pvdz 机组吗？

作者
Author: xiaobogaga 时间: 2015-10-1 18:22

sobereva 发表于 2015-10-1 04:25
完全没有可比性嘛
1 一个是单点一个是counterpoise任务
2 一个对配体用cc-pVDZ一个用aug-cc-pVDZ

那sob 老师，如果gaussian 命令行里只写方法/基组是不是就是计算单点能了？

作者
Author: sobereva 时间: 2015-10-1 22:31

xiaobogaga 发表于 2015-10-1 18:22
那sob 老师，如果gaussian 命令行里只写方法/基组是不是就是计算单点能了？

是

作者
Author: sobereva 时间: 2015-10-1 22:42

xiaobogaga 发表于 2015-10-1 18:20
你好，我是新手，请问正确的输入应该怎么写，比如我用Gamess 计算的时候用MP2(full)，在Gamess的输入文件 ...

ACCD就是aug-cc-pVDZ
不想冻核就写
$MP2 NACORE=0 $END

作者
Author: xiaobogaga 时间: 2015-10-2 10:42

sobereva 发表于 2015-10-1 22:42
ACCD就是aug-cc-pVDZ
不想冻核就写
$MP2 NACORE=0 $END

非常感谢sob老师的耐心指导！！

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