计算化学公社

标题: GAMESS计算EDA，极化能量（POLARIZATION ENERGY ）特别高 [打印本页]

作者
Author: yang_lzu 时间: 2017-12-19 15:51
标题: GAMESS计算EDA，极化能量（POLARIZATION ENERGY ）特别高
老师：
   您好，这是我进行能量分解分析的输入文件和输出文件，不知道为什么我计算出来的那个极化能会如此高，另外我改变了方法和和基组还是会发生同样的情况。但是奇怪的是，我利用这个输入文件计算另外一个相同体系的不同结合方式其中有一个极化能就正常，请问老师我问题是出在哪里呢？
输入文件：
$contrl scftyp=rhf dfttyp=m06-2x runtyp=eda ISPHER=1 $end
$system timlim=123456 $end
$basis gbasis=accd $end
$guess guess=huckel $end
$lmoeda matom(1)=11 12 mcharg(1)=0 0 mmult(1)=1 1 $end
$system mwords=800 $end
$data
Molecule specification
C1
C  6.000000 4.231722 -1.226748 0.291319
C  6.000000 4.354037 -0.503653 -0.896659
C  6.000000 4.303158 0.887516 -0.835250
C  6.000000 4.132816 1.491394 0.407171
C  6.000000 4.019830 0.674902 1.533389
N  7.000000 4.066637 -0.659057 1.489587
H  1.000000 4.389179 1.488297 -1.740299
H  1.000000 4.265084 -2.317510 0.277363
H  1.000000 4.485494 -1.024195 -1.843833
H  1.000000 4.081960 2.573980 0.508336
H  1.000000 3.882877 1.117993 2.521117
C  6.000000 -2.829329 1.317416 0.048227
C  6.000000 -1.464361 1.155150 -0.156006
C  6.000000 -0.896644 -0.115047 -0.206563
C  6.000000 -1.723368 -1.224060 -0.049602
C  6.000000 -3.089092 -1.071207 0.155292
C  6.000000 -3.643213 0.202159 0.204312
F  9.000000 -3.869057 -2.139165 0.301830
F  9.000000 -3.360160 2.537222 0.092529
F  9.000000 -0.709567 2.241260 -0.301750
F  9.000000 -1.220764 -2.454752 -0.093085
F  9.000000 -4.949132 0.352290 0.398208
Br  35.000000 0.943821 -0.330965 -0.477233
$END
输出文件：
-------------
OWN BASIS SET                                     HARTREE       KCAL/MOL
-------------
ELECTROSTATIC ENERGY                ES=       -0.003284          -2.06
EXCHANGE ENERGY                      EX=       -0.004069          -2.55
REPULSION ENERGY                   REP=       0.015372             9.65
POLARIZATION ENERGY                POL=    3550.050693       2227692.31
DFT DISPERSION ENERGY             DISP=       -0.011223          -7.04
TOTAL INTERACTION ENERGY HF OR DFT    E=    3550.047490       2227690.30

作者
Author: wei 时间: 2017-12-19 18:46
没有收敛

作者
Author: yang_lzu 时间: 2017-12-19 20:42

wei 发表于 2017-12-19 18:46
没有收敛

像碰到这种问题一般要怎么解决呢？需要加什么关键词吗？还是要调节结构什么的？因为这是第一次使用GAMESS，很多还不是很懂。

作者
Author: sobereva 时间: 2017-12-19 20:59

yang_lzu 发表于 2017-12-19 20:42
像碰到这种问题一般要怎么解决呢？需要加什么关键词吗？还是要调节结构什么的？因为这是第一次使用GAMESS ...

看本论坛GAMESS分类里的帖子，很多说解决SCF收敛问题的

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