计算化学公社

标题: Gamess做能量分解，能量和Gauss做出来的差很多的问题 [打印本页]

作者
Author: 小靖靖 时间: 2023-3-11 20:39
标题: Gamess做能量分解，能量和Gauss做出来的差很多的问题
各位老师，我用GAMESS做gks-EDA,做出来的能量和我用GAUSSIAN算出来的差很多，用的是相同的方法和基组(PBE0/def2-tzvp)，GAUSSIAN算的时候加了BSSE矫正，并且用的是二体减去两个单体，单体的坐标是从优化好的二体里面拿坐标直接单点的，GAUSSIAN算出来的能量是-81.98kj/mol,能量分解算出来的总能量是-28.83kj/mol,然后，我用GAMESS优化了单体和二体，算出来的能量是-37.58kj/mol,问题不知道出在哪里了，请各位老师指点指点。下面是GAMESS做能量分解的输入文件，和输出文件的截图，以及用GAMESS优化的输入文件和输出文件。

作者
Author: zjxitcc 时间: 2023-3-12 01:10
本帖最后由 zjxitcc 于 2023-3-12 01:19 编辑

我用frag_guess_wfn小程序试了下，输入文件为 (, 下载次数 Times of downloads: 7)

frag_guess_wfn会自动调用高斯做DFT计算、产生GKS-EDA输入文件，过程中会打印片段能量、总相互作用能结果等信息

i= 1, frags(i)%e = -2299.670950960, frags(i)%ssquare= 0.00
i= 2, frags(i)%e = -2299.670950960, frags(i)%ssquare= 0.00
i= 3, frags(i)%e = -2299.672252700, frags(i)%ssquare= 0.00
i= 4, frags(i)%e = -2299.672252700, frags(i)%ssquare= 0.00
i= 5, frags(i)%e = -4599.353574560, frags(i)%ssquare= 0.00
If you are performing a GKS-EDA calculation using DFT or a LMO-EDA calculation
using HF, then the following energy should be close to 'TOTAL INTERACTION ENERGY'
of 'OWN BASIS SET' section in GAMESS output:
-7.32 kcal/mol
and the following energy should be close to 'TOTAL INTERACTION ENERGY' of
'ALL BASIS SET' section in GAMESS output:
-5.69 kcal/mol

复制代码

产生AuCl3-Cl-dimer.inp文件，提交给XEDA，很快算完，结果如下

-------------
OWN BASIS SET HARTREE KCAL/MOL
-------------
ELECTROSTATIC ENERGY ES= -0.023195 -14.56
EXCHANGE ENERGY EX= -0.050468 -31.67
REPULSION ENERGY REP= 0.083485 52.39
POLARIZATION ENERGY POL= -0.008302 -5.21
ELECTRON CORRELATION CORR= -0.013246 -8.31
TOTAL INTERACTION ENERGY E= -0.011726 -7.36
-------------
ALL BASIS SET HARTREE KCAL/MOL
-------------
ELECTROSTATIC ENERGY ES= -0.023642 -14.84
EXCHANGE ENERGY EX= -0.052073 -32.68
REPULSION ENERGY REP= 0.086453 54.25
POLARIZATION ENERGY POL= -0.007504 -4.71
ELECTRON CORRELATION CORR= -0.012308 -7.72
TOTAL INTERACTION ENERGY E= -0.009074 -5.69

复制代码

可见两个程序的总相互作用能可以严格对应上。所以你的高斯计算应该有问题，可以对照frag_guess_wfn的输出结果逐一对比。（这也侧面反映出使用frag_guess_wfn的好处，不仅节省计算时间，两个程序还能严格对应、不会因人工操作出错）

作者
Author: 小靖靖 时间: 2023-3-12 09:56

zjxitcc 发表于 2023-3-12 01:10
我用frag_guess_wfn小程序试了下，输入文件为

frag_guess_wfn会自动调用高斯做DFT计算、产生GKS-EDA输入 ...

好的，谢谢老师，我在好好检查一下

作者
Author: 小靖靖 时间: 2023-3-14 15:21

zjxitcc 发表于 2023-3-12 01:10
我用frag_guess_wfn小程序试了下，输入文件为

frag_guess_wfn会自动调用高斯做DFT计算、产生GKS-EDA输入 ...

老师，您好，我想请教一个问题，我GAUSSIAN算的时候，用的是PBE0-D3/def2-tzvp,这样的话，用GAMESS做gks-eda的时候，是不是也应该加D3呀？如果想在GAMESS加D3的话，我该加什么关键词呢？

作者
Author: zjxitcc 时间: 2023-3-14 16:16

小靖靖发表于 2023-3-14 15:21
老师，您好，我想请教一个问题，我GAUSSIAN算的时候，用的是PBE0-D3/def2-tzvp,这样的话，用GAMESS做gks- ...

如果你是使用frag_guess_wfn小程序产生的GKS-EDA输入文件，那么文件里D3关键词自动写好，省事方便。如果你要手动书写GKS-EDA输入文件，那你得自己写D3关键词，见GAMESS文档gamess/docs-input.txt的IDCVER

作者
Author: 小靖靖 时间: 2023-3-14 18:19

zjxitcc 发表于 2023-3-14 16:16
如果你是使用frag_guess_wfn小程序产生的GKS-EDA输入文件，那么文件里D3关键词自动写好，省事方便。如果 ...

好的，谢谢老师

作者
Author: 小靖靖 时间: 2023-3-15 08:44

zjxitcc 发表于 2023-3-14 16:16
如果你是使用frag_guess_wfn小程序产生的GKS-EDA输入文件，那么文件里D3关键词自动写好，省事方便。如果 ...

好的，非常感谢老师，关于wfn小程序，实不相瞒，不是不想用，是不会用，我和我同学一起讨论也不太会，可能是我太笨

作者
Author: zjxitcc 时间: 2023-3-15 14:15

小靖靖发表于 2023-3-15 08:44
好的，非常感谢老师，关于wfn小程序，实不相瞒，不是不想用，是不会用，我和我同学一起讨论也不太会，可 ...

是frag_guess_wfn，不是wfn。关于frag_guess_wfn的最简易使用方式
https://jeanwsr.gitlab.io/mokit- ... want-frag_guess_wfn

作者
Author: 小靖靖 时间: 2023-3-15 19:36

zjxitcc 发表于 2023-3-15 14:15
是frag_guess_wfn，不是wfn。关于frag_guess_wfn的最简易使用方式
https://jeanwsr.gitlab.io/mokit-doc ...

好的，谢谢老师

作者
Author: btm 时间: 2023-3-20 19:15
请问你这个”GAUSSIAN算出来的能量是-81.98kj/mol,能量分解算出来的总能量是-28.83kj/mol“这里的-28.83是怎么来的在输出文件中搜索什么呀？

作者
Author: 小靖靖 时间: 2023-3-20 20:28

btm 发表于 2023-3-20 19:15
请问你这个”GAUSSIAN算出来的能量是-81.98kj/mol,能量分解算出来的总能量是-28.83kj/mol“这里的-28.83是 ...

倒数几行就是

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

Powered by Discuz! X3.3