gamess单点能计算与gaussian结果不匹配求助

m494081672

因为感觉做的LMOEDA有问题，所以看了一下论坛的建议，从单点能开始算想检查一下，用的是gaussian的nosymm 10f 6d punch=gamess int=nobasistransform导出的基组，但是有一个单体的gamess 的能量却比gaussian高了100多个hartree，其他单体也有差异。自己感觉输入文件没有写错啊，不知道是什么原因。请各位帮忙看看。gaussian和gamess的输入输出文件都在附件中。

m494081672

迟了一两天回复，谢谢各位大神，确实是这个核电荷数目的原因，改了之后，和高斯的能量差只有10^-6hartree的差距，后面的LMOEDA也正常了，太感谢了！
之前我一直以为只要是用高斯命令导入的基组一定不会错来着，看来是我太年轻了。

zjxitcc

高斯的punch=gamess关键词存在2个已知的bug，一个是使用赝势时，输出的核电荷数目不符合GAMESS标准（GAMESS仍采取等于原子序数的核电荷数）；另一个是角动量超过G时出现问号？输出，不懂得用H, I代替（此会导致GAMESS直接报错，因此容易发现）。

我曾经向高斯官方发邮件反映过，他们也说会向开发者反映这个问题，但至少在G16 B01里这个问题还是存在。

yflchx

m494081672 发表于 2020-4-15 17:15
感谢您的回答：
1、我用gamess的B3LYPV1R，也算过了，结果和gamess的B3LYP，只差0.0001hatree，对这个体 ...

你的GAMESS的输入文件里面：

Co       17.       -4.7638975148        0.5869923980        1.1116030469


仔细看，问题是不是在这里。

m494081672

yflchx 发表于 2020-4-15 16:43
1. 直接写B3LYP，Gaussian和GAMESS是不一样的，VWN用的不一样。

2. 检查确保结构一致，基组一致，积分格 ...

感谢您的回答：
1、我用gamess的B3LYPV1R，也算过了，结果和gamess的B3LYP，只差0.0001hatree，对这个体系基本没有什么影响。不过，严格一些应该是要写B3LYPV1R的，这个文件中是忘了改。
2、结构、基组、我是用gaussian的命令转化的，应该不会有问题。
至于积分格点和收敛到的电子态，我不太明白怎么看，您能告诉我吗？
3、因为我分片段的时候有断开化学键，想要化学键和弱相互作用一起分析，所以我觉得可能相比SAPT，LMOEDA好一些。至于色散项，我参考了帖子http://bbs.keinsci.com/thread-682-1-1.html中的说法， 想要用DFT-D来额外校正。

yflchx

1. 直接写B3LYP，Gaussian和GAMESS是不一样的，VWN用的不一样。

2. 检查确保结构一致，基组一致，积分格点一致，收敛到的电子态一致。

3. 如果有可能的话，建议用PSI4做SAPT能量分解。体系大的话尝试：SAPT2+/aug-cc-pVDZ和sSAPT0/jun-cc-pVDZ。

B3LYP直接做LMO-EDA会有DISP项，用GKS-EDA才合适。