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

gaolei0915

以下是我的输入文件，第一次用GAmess做能量分解，不知道问题出在哪里，求老师们指教

$CONTRL SCFTYP=RHF RUNTYP=eda ICHARG=0 ispher=1 mplevl=2
  MULT=1 $end
!$guess guess=hcore $end
$BASIS GBASIS=acct $END
$LMOEDA MATOM(1)=8 4 MCHARG(1)=0 0 MMULT(1)=1 1 $END
!$SCF DIIS=.t. SOSCF=.f. DIRSCF=.T. FDIFF=.F. NPUNCH=1 $END
$scf maxdii=100 $end
$SYSTEM TIMLIM=99999999 MWORDS=2000 memddi=2000 $END
$DATA
optimized by mp2/augt
C1
Ga        31.0     0.53340000    0.87176400    0.00000000
Ga        31.0    -0.21670200   -1.63795300    0.00000000
H         1.0      1.06457300    1.24150300    1.40397800  
H         1.0      1.06457300    1.24150300   -1.40397800
H         1.0     -0.59928900   -2.22881500    1.37283900
H         1.0     -0.59928900   -2.22881500   -1.37283900
H         1.0     -0.89178500   -0.05597900    0.00000000
H         1.0      1.24599900   -0.93717400    0.00000000
N         7.0     -1.00713700    2.58269800    0.00000000
H         1.0     -1.60566800    2.55207900   -0.81787100
H         1.0     -0.54175400    3.48387500    0.00000000
H         1.0     -1.60566700    2.55207900    0.81787100
$END

输出的能量部分


-------------
ALL BASIS SET                                       HARTREE          KCAL/MOL
-------------
ELECTROSTATIC ENERGY                  ES=         -0.069641            -43.70
EXCHANGE ENERGY                       EX=         -0.085443            -53.62
REPULSION ENERGY                     REP=          0.163149            102.38
POLARIZATION ENERGY                  POL=         -0.027247            -17.10
MP2 DISPERSION ENERGY               DISP=         -0.005576             -3.50
TOTAL INTERACTION ENERGY HF OR DFT     E=         -0.019182            -12.04
TOTAL INTERACTION ENERGY MP2           E=         -0.024759            -15.54





而我用gauss计算出的相互作用能是-7.04kcal，为什么差了这么多呢？

brothers

你用gaussian怎么算的？

sobereva

在Gaussian里也应当在MP2/aug-cc-pVTZ下计算相互作用能。而且应该先直接比较Gaussian和GAMESS-US在MP2/aug-cc-pVTZ级别下计算的相互作用能，看看是否有明显差异，若没什么差异再去和EDA的结果比较。
凡是碰到复杂搞不清楚的问题，先把问题简化。

gaolei0915

sobereva 发表于 2015-4-25 22:51
在Gaussian里也应当在MP2/aug-cc-pVTZ下计算相互作用能。而且应该先直接比较Gaussian和GAMESS-US在MP2/aug- ...

谢谢，我在gauss里也是同样地方法基组计算的，然后用优化好的坐标直接用Gamess做的能量分解，您的意思是我应该用Gamess也做一下优化，看优化的结果是否和Gauss的一致是吗？

gaolei0915

brothers 发表于 2015-4-25 20:42
你用gaussian怎么算的？

同样地方法基组优化的，优化好以后用优化好的坐标做的能量分解，但是GAmess这里给出的Tatal interaction energy 和Gauss计算出来用复合物减去两个单体的能量对不上

brothers

gaolei0915 发表于 2015-4-26 12:28
同样地方法基组优化的，优化好以后用优化好的坐标做的能量分解，但是GAmess这里给出的Tatal interaction  ...

gamess能量分解中的interaction energy不包含A与B从各自的稳态弛豫到AB复合物中相应结构这部分能量（E_preparation，大于0）。在gaussian里，分别优化A、B，然后优化AB，然后计算AB-A-B能量，得到的是interaction energy+E_preparation。不知道你有没有计算E_preparation？
不过，看你的结果，如果是没有计算E_preparation的话，好像E_preparation会有8.5啊，不知道比例是不是有点高。我做的也不多，没见过这么高的

sobereva

gaolei0915 发表于 2015-4-26 12:25
谢谢，我在gauss里也是同样地方法基组计算的，然后用优化好的坐标直接用Gamess做的能量分解，您的意思是 ...

理论计算相互作用能的时候一般不考虑片段的变形能。也就是说，应当只优化复合物，不要优化单体，计算单体能量所用的结构直接从优化后的复合物结构中提取。

先用Gaussian在MP2/aug-cc-pVTZ下优化复合物。然后基于此结构，分别在Gaussian和GAMESS-US下计算复合物和片段的能量，并求差值，看看两个程序给出来的差多少。如果相差明显，再比较两个程序给出的复合物或片段的能量，如果也有差异，说明光是单点能计算就已经对不上了，再去考虑是否是什么其它因素引起的。一点一点把问题拆解分析，总能搞清楚原因。

gaolei0915

brothers 发表于 2015-4-26 14:21
gamess能量分解中的interaction energy不包含A与B从各自的稳态弛豫到AB复合物中相应结构这部分能量（E_pr ...

谢谢，对我很有帮助，我试一下

gaolei0915

sobereva 发表于 2015-4-26 18:27
理论计算相互作用能的时候一般不考虑片段的变形能。也就是说，应当只优化复合物，不要优化单体，计算单 ...

谢谢您，对我很有帮助，我计算一下，十分感谢