计算化学公社

 找回密码 Forget password
 注册 Register
Views: 9165|回复 Reply: 3
打印 Print 上一主题 Last thread 下一主题 Next thread

[GAMESS-US] GAMESS计算EDA,极化能量(POLARIZATION ENERGY )特别高

[复制链接 Copy URL]

15

帖子

0

威望

765

eV
积分
780

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
老师:
     您好,这是我进行能量分解分析的输入文件和输出文件,不知道为什么我计算出来的那个极化能会如此高,另外我改变了方法和和基组还是会发生同样的情况。但是奇怪的是,我利用这个输入文件计算另外一个相同体系的不同结合方式其中有一个极化能就正常,请问老师我问题是出在哪里呢?
输入文件:
$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

93

帖子

0

威望

1788

eV
积分
1881

Level 5 (御坂)

2#
发表于 Post on 2017-12-19 18:46:10 | 只看该作者 Only view this author
没有收敛

15

帖子

0

威望

765

eV
积分
780

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2017-12-19 20:42:28 | 只看该作者 Only view this author

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

6万

帖子

99

威望

6万

eV
积分
125151

管理员

公社社长

4#
发表于 Post on 2017-12-19 20:59:41 | 只看该作者 Only view this author
yang_lzu 发表于 2017-12-19 20:42
像碰到这种问题一般要怎么解决呢?需要加什么关键词吗?还是要调节结构什么的?因为这是第一次使用GAMESS ...

看本论坛GAMESS分类里的帖子,很多说解决SCF收敛问题的
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 08:03 , Processed in 0.168283 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list