计算化学公社

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

作者
Author: gaolei0915 时间: 2015-4-25 17:33
标题: Gamess做能量分解，能量和Gauss做出来的差很多的问题
以下是我的输入文件，第一次用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，为什么差了这么多呢？

作者
Author: brothers 时间: 2015-4-25 20:42
你用gaussian怎么算的？

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

作者
Author: gaolei0915 时间: 2015-4-26 12:25

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

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

作者
Author: gaolei0915 时间: 2015-4-26 12:28

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

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

作者
Author: brothers 时间: 2015-4-26 14:21
本帖最后由 brothers 于 2015-4-27 10:44 编辑

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啊，不知道比例是不是有点高。我做的也不多，没见过这么高的

作者
Author: sobereva 时间: 2015-4-26 18:27

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

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

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

作者
Author: gaolei0915 时间: 2015-4-27 08:29

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

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

作者
Author: gaolei0915 时间: 2015-4-27 08:29

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

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

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