Gamess做能量分解，能量和Gauss做出来的差很多的问题

小靖靖

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

zjxitcc

我用frag_guess_wfn小程序试了下，输入文件为 AuCl3-Cl-dimer.gjf (2.25 KB, 下载次数 Times of downloads: 7)

2023-3-12 01:10 上传 Uploaded

点击下载
Click to download

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

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

复制代码

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

1. -------------
   
2. OWN BASIS SET                                       HARTREE          KCAL/MOL
   
3. -------------
   
4. ELECTROSTATIC ENERGY                  ES=         -0.023195            -14.56
   
5. EXCHANGE ENERGY                       EX=         -0.050468            -31.67
   
6. REPULSION ENERGY                     REP=          0.083485             52.39
   
7. POLARIZATION ENERGY                  POL=         -0.008302             -5.21
   
8. ELECTRON CORRELATION                CORR=         -0.013246             -8.31
   
9. TOTAL INTERACTION ENERGY               E=         -0.011726             -7.36
   
10. 
    
11. 
    
12. -------------
    
13. ALL BASIS SET                                       HARTREE          KCAL/MOL
    
14. -------------
    
15. ELECTROSTATIC ENERGY                  ES=         -0.023642            -14.84
    
16. EXCHANGE ENERGY                       EX=         -0.052073            -32.68
    
17. REPULSION ENERGY                     REP=          0.086453             54.25
    
18. POLARIZATION ENERGY                  POL=         -0.007504             -4.71
    
19. ELECTRON CORRELATION                CORR=         -0.012308             -7.72
    
20. TOTAL INTERACTION ENERGY               E=         -0.009074             -5.69

复制代码

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

小靖靖

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

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

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

小靖靖

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的话，我该加什么关键词呢？

zjxitcc

小靖靖 发表于 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

小靖靖

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

好的，谢谢老师

小靖靖

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

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

zjxitcc

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

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

小靖靖

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

好的，谢谢老师

btm

请问你这个”GAUSSIAN算出来的能量是-81.98kj/mol,能量分解算出来的总能量是-28.83kj/mol“这里的-28.83是怎么来的在输出文件中搜索什么呀？

小靖靖

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

倒数几行就是