能量分解输出文件没有色散项，而且总能量为正

gyl · 发表于 Post on 2020-8-29 09:40:36

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

，下面是输入文件和输出文件的截图，和生成fort7的输入文件截图

zjxitcc · 发表于 Post on 2020-8-29 10:44:48

本帖最后由 zjxitcc 于 2020-8-29 13:32 编辑

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

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

sobereva · 发表于 Post on 2020-8-29 14:07:43

原始的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

gyl · 发表于 Post on 2020-8-29 17:11:11

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

提醒你注意 ...

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

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

zjxitcc · 发表于 Post on 2020-8-29 17:20:29

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[GAMESS-US] 能量分解输出文件没有色散项，而且总能量为正