计算化学公社

标题: GAMESS算能量分解输入及输出文件分析求助 [打印本页]

作者
Author: Finesnow 时间: 2019-5-21 19:17
标题: GAMESS算能量分解输入及输出文件分析求助
各位老师好，我用Gaussian09优化构型并计算的单点能，然后用GAMESS算的能量分解，遇到的问题如下：
1.高斯里用的B3LYP-D3(BJ)方法，在GAMESS里输入IDCVER=4，报错显示ERROR !! IDCVER MUST BE SMALLER THAN 3，但是以我对输入文件的（浅显的）理解是D3(BJ)对应IDCVER=4，D3对应的才是IDCVER=3；
2.IDCVER改为3后可以开始运行也可以正常结束，但计算结果与高斯差着将近一个数量级，Gaussian计算出的相互作用能为-0.0269 HARTREE，如下为GAMESS计算结果：（Gaussian利用的scrf=smd，Gamess里加的关键词为$SMD  solvnt=water $end，是这里出错了么）
-------------
ALL BASIS SET                                     HARTREE             KCAL/MOL
-------------
ELECTROSTATIC ENERGY       ES=       -0.193221          -121.25
EXCHANGE ENERGY                EX=       -0.090653          -56.89
REPULSION ENERGY             REP=       0.173074          108.61
POLARIZATION ENERGY       POL=       -0.079837          -50.10
GRIMME DISP CORRECTION    DC=       -0.005448          -3.42
DFT CORRELATION          CORR=       -0.004171          -2.62
TOTAL INTERACTION ENERGY HF OR DFT E=  -0.200257          -125.66

诚心求助各位老师能帮我解答一下疑惑，不胜感激

作者
Author: zjxitcc 时间: 2019-5-21 20:44
看看GAMESS输出文件里每处的SCF是否都正常收敛了。

作者
Author: Finesnow 时间: 2019-5-21 22:22

zjxitcc 发表于 2019-5-21 20:44
看看GAMESS输出文件里每处的SCF是否都正常收敛了。

老师您好，每次提问都能得到您的及时答复，我太感动太感谢啦

惭愧的是，我对GAMESS知之甚少。您说的SCF收敛我理解的是每步显示density converged就可以么，我的输出文件五步里都显示了density converged

作者
Author: sobereva 时间: 2019-5-22 02:59
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素（仔细看http://bbs.keinsci.com/thread-9348-1-1.html），我已把你的标题“能量分解输入及输出文件分析 ”改了，以后务必注意

作者
Author: Finesnow 时间: 2019-5-22 08:59

sobereva 发表于 2019-5-22 02:59
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素（仔 ...

好的，谢谢您的提醒，下次一定记得啦

作者
Author: zjxitcc 时间: 2019-5-22 15:20
我不知道你如何办到的，至少在我用的2017版本里，LMOEDA不支持DFT-D，所以B3LYP-D3(BJ)做不了LMOEDA。
光做简单的单点计算的话，B3LYP-D3(BJ)是可以用的。

作者
Author: Finesnow 时间: 2019-5-22 16:15

zjxitcc 发表于 2019-5-22 15:20
我不知道你如何办到的，至少在我用的2017版本里，LMOEDA不支持DFT-D，所以B3LYP-D3(BJ)做不了LMOEDA。
光 ...

原来如此，我是2012版的，所以更不能行吧。我刚刚尝试了一下Gaussian和GAMESS都去掉SMD关键词，这样一来计算结果几乎是一样的。只是对于Gaussian来说，在溶剂化模型下优化的构型，在真空下计算相互作用能也是可以的吧？

作者
Author: zjxitcc 时间: 2019-5-23 10:48

Finesnow 发表于 2019-5-22 16:15
原来如此，我是2012版的，所以更不能行吧。我刚刚尝试了一下Gaussian和GAMESS都去掉SMD关键词，这样一来 ...

没见过这么做的，只见过反过来的。
你这么做虽然程序步骤上允许，但是缺乏实际意义。

作者
Author: Finesnow 时间: 2019-5-23 15:24

zjxitcc 发表于 2019-5-23 10:48
没见过这么做的，只见过反过来的。
你这么做虽然程序步骤上允许，但是缺乏实际意义。

嗯嗯，但我选的体系在真空下优化总是得不到预期构型，所以选择了在溶剂模型下优化。然后用GAMESS做能量分解的时候即使加了溶剂关键词计算出来的总能量却约为高斯溶剂下得到的相互作用能的十倍（事实上几乎等于溶剂下优化、真空下算得的相互作用能。不知道是不是GAMESS里我的溶剂关键词加错了，接下来我再学习尝试一下这里）。这种情况如果都用真空下的相互作用能做一下定性分析是可以接受的吧？

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