计算化学公社

标题: 能量分解输出文件没有色散项，而且总能量为正 [打印本页]

作者
Author: gyl 时间: 2020-8-29 09:40
标题: 能量分解输出文件没有色散项，而且总能量为正
刚接触能量分解，有些问题想请教一下大神！我想用m062x/def2tzvp做能量分解，针对弱作用，看过一些帖子，说用高斯可以生成fort7的文件，然后加一下关键词，就能保存gamess的输入文件，虽然能正常结束，但是分解出来没有色散项，而且总能量为正的，不知道怎么回事，麻烦各位大神能指点一下

，下面是输入文件和输出文件的截图，和生成fort7的输入文件截图 (, 下载次数 Times of downloads: 24) (, 下载次数 Times of downloads: 16)

作者
Author: zjxitcc 时间: 2020-8-29 10:44
本帖最后由 zjxitcc 于 2020-8-29 13:32 编辑

原理我是不太懂，不过操作上，可能你要算DFT或>=MP2级别，输出内容才会有dispersion项

提醒你注意：高斯输入文件内不用写方法名称，写了也没用，因为fort.7里面的关键词很简陋，需要自己修改的。你的inp文件里面就是基于RHF的能量分解，而不是基于M062X的，不知道你是没意识到、还是故意贴个RHF的出来。基于RHF的能量分解结果就不用看了，用点像样的计算方法。

作者
Author: sobereva 时间: 2020-8-29 14:07
原始的morokuma形式的EDA只能用于HF，没法用于DFT
如果你做更普适的LMOEDA，并需要获得色散作用项，需要在后HF下做，LMOEDA会把后HF相关能的变化视为色散作用能，但这样做并非总合理。J. Phys. Chem. Lett. 2014, 5, 1380-1385这篇文章指出，某些情况下，这部分能量当中其实掺杂着少许静电作用，不纯粹是色散作用。（若用LMOEDA结合普通泛函，会把复合物的DFT相关能和单体的DFT相关能求差来获得色散能，而这种色散作用能的定义显然是错误的）
用修改版GAMESS-US做GKS-EDA的话，会把DFT-D色散校正能作为色散作用能，这么考虑是否合理取决于用什么泛函。下文里有相关讨论
DFT-D色散校正的使用
http://sobereva.com/210

作者
Author: gyl 时间: 2020-8-29 17:11

zjxitcc 发表于 2020-8-29 10:44
原理我是不太懂，不过操作上，可能你要算DFT或>=MP2级别，输出内容才会有dispersion项

提醒你注意 ...

好的，辛苦了，非常感谢你吆

！还有一个小问题，如果用MOKIT程序中的fch2inp小程序从高斯fchk文件生成GAMESS .inp文件，算出来是不是也没有色散项，因为按照sobereva老师所说，我这个方法不属于HF方法。

作者
Author: zjxitcc 时间: 2020-8-29 17:20

gyl 发表于 2020-8-29 17:11
好的，辛苦了，非常感谢你吆！还有一个小问题，如果用MOKIT程序中的fch2inp小程序从高斯fchk文件生 ...

恩，没有。fch2inp小程序只是产生坐标、基组和轨道而已。计算方法还是得你自己设。是不属于后HF方法，你少了“后”字。

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