计算化学公社

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

[GAMESS-US] Gamess能量分解结果不对,数值特别大,请帮我看一下哪里出错了

[复制链接 Copy URL]

30

帖子

0

威望

227

eV
积分
257

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
Gamess能量分解可以算了,但是结果不对,数值特别大,请老师帮我看一下哪里出错了吧

下面是out文件的结果部分
OWN BASIS SET                                       HARTREE          KCAL/MOL
-------------
ELECTROSTATIC ENERGY                  ES=         -1.098519           -689.33
EXCHANGE ENERGY                       EX=         33.763408          21186.88
REPULSION ENERGY                     REP=         -1.356356           -851.13
POLARIZATION ENERGY                  POL=      -3163.870635       -1985360.46
DFT DISPERSION ENERGY               DISP=          1.536075            963.90
TOTAL INTERACTION ENERGY HF OR DFT     E=      -3131.026028       -1964750.14


-------------
ALL BASIS SET                                       HARTREE          KCAL/MOL
-------------
ELECTROSTATIC ENERGY                  ES=         -1.997901          -1253.70
EXCHANGE ENERGY                       EX=         -0.000000             -0.00
REPULSION ENERGY                     REP=          3.672155           2304.31
POLARIZATION ENERGY                  POL=         -1.530527           -960.42
DFT DISPERSION ENERGY               DISP=          0.000000              0.00
TOTAL INTERACTION ENERGY HF OR DFT     E=          0.143727             90.19

-------------------------------------
IN SUPER MOL (FOR SPECIAL COMPARISON)               HARTREE          KCAL/MOL
-------------------------------------
ELECTROSTATIC ENERGY                  ES=       -148.588892         -93241.02
EXCHANGE ENERGY                       EX=         -0.000000             -0.00
DFT CORRELATION                     CORR=          0.000000              0.00
TOTAL INTERACTION ENERGY HF OR DFT     E=       -263.140732        -165123.44


AuCl-BNH.inp

5.69 KB, 下载次数 Times of downloads: 8

334

帖子

1

威望

6155

eV
积分
6509

Level 6 (一方通行)

2#
发表于 Post on 2021-1-13 08:55:19 | 只看该作者 Only view this author
检查一下两个monomer 和整体的SCF是不是都收敛了。 GAMESS 在计算的能量分解的时候即使SCF不收敛也不会报错 而是继续算下去。尤其是Au Cl 那个部分 总觉得怕是不好收敛。

30

帖子

0

威望

227

eV
积分
257

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-1-13 09:39:24 | 只看该作者 Only view this author
请问怎么Gamess怎么看收敛了么,什么关键词呀?和高斯里一样么

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

4#
发表于 Post on 2021-1-13 10:03:06 | 只看该作者 Only view this author
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,避免有任何歧义,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“Gamess能量分解”改了,以后务必注意,下次将扣分处理
北京科音自然科学研究中心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

1245

帖子

1

威望

3860

eV
积分
5125

Level 6 (一方通行)

5#
发表于 Post on 2021-1-13 10:40:56 | 只看该作者 Only view this author
前面scf的收敛可能出问题了

4289

帖子

4

威望

9546

eV
积分
13915

Level 6 (一方通行)

MOKIT开发者

6#
发表于 Post on 2021-1-13 11:50:40 | 只看该作者 Only view this author
xiaoshan 发表于 2021-1-13 09:39
请问怎么Gamess怎么看收敛了么,什么关键词呀?和高斯里一样么

看收敛是指打开输出文件用肉眼看。关键词又不是肉眼,不能看收敛性。在输出文件里搜ITER EX即可
自动做多参考态计算的程序MOKIT

30

帖子

0

威望

227

eV
积分
257

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2021-1-13 14:45:33 | 只看该作者 Only view this author
sobereva 发表于 2021-1-13 10:03
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,避免 ...

好的老师,下次注意

30

帖子

0

威望

227

eV
积分
257

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2021-1-14 15:23:26 | 只看该作者 Only view this author
paramecium86 发表于 2021-1-13 08:55
检查一下两个monomer 和整体的SCF是不是都收敛了。 GAMESS 在计算的能量分解的时候即使SCF不收敛也不会报错 ...

请问是看用Gaussian计算时看两个单体和整体的scf有没有收敛么?单体和整体的都收敛了。
我是分别优化单体,再优化二体,二体进行BSSE校正,做能量分解时复制校正后的原子坐标,这样做对么?

30

帖子

0

威望

227

eV
积分
257

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2021-1-14 15:25:52 | 只看该作者 Only view this author
zjxitcc 发表于 2021-1-13 11:50
看收敛是指打开输出文件用肉眼看。关键词又不是肉眼,不能看收敛性。在输出文件里搜ITER EX即可

您好,在输出文件里搜ITER EX有很多数据,老师我想问,从这些数据里能看出什么呢?能不能麻烦您帮我看一下我的输入文件有没有问题呀

4289

帖子

4

威望

9546

eV
积分
13915

Level 6 (一方通行)

MOKIT开发者

10#
发表于 Post on 2021-1-14 15:54:13 | 只看该作者 Only view this author
xiaoshan 发表于 2021-1-14 15:25
您好,在输出文件里搜ITER EX有很多数据,老师我想问,从这些数据里能看出什么呢?能不能麻烦您帮我看一 ...

每处ITER EX代表算一个单点,能量分解要算好几个单点,你直接看最后一次出现的ITER EX。从后往前搜即可
自动做多参考态计算的程序MOKIT

334

帖子

1

威望

6155

eV
积分
6509

Level 6 (一方通行)

11#
发表于 Post on 2021-1-14 19:42:54 | 只看该作者 Only view this author
本帖最后由 paramecium86 于 2021-1-14 19:44 编辑
xiaoshan 发表于 2021-1-14 15:23
请问是看用Gaussian计算时看两个单体和整体的scf有没有收敛么?单体和整体的都收敛了。
我是分别优化单 ...

高斯里收敛不代表在GAMESS里可以收敛。你的问题出在GAMESS的计算自然是在GAMESS的输出文件里查看。SCF的部分在GAMESS里都是一大段一大段的如果你耐心从头看到尾是肯定不会错过的。你如果从来都没有从头到尾仔细看过一遍GAMESS我是建议你好好的把GAMESS的输出文件好好看一遍(这种能量分解这种任务一般的输出文件不会太长也。几千行而已)。如果真的要搜索某一个关键词 如果里面真的有不收敛的地方 那可以搜索    SCF IS UNCONVERGED

30

帖子

0

威望

227

eV
积分
257

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2021-1-15 08:47:46 | 只看该作者 Only view this author
paramecium86 发表于 2021-1-14 19:42
高斯里收敛不代表在GAMESS里可以收敛。你的问题出在GAMESS的计算自然是在GAMESS的输出文件里查看。SCF的 ...

老师,我的确实没有收敛,遇到这种情况应该怎么办呢?

30

帖子

0

威望

227

eV
积分
257

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2021-1-15 08:48:15 | 只看该作者 Only view this author
zjxitcc 发表于 2021-1-14 15:54
每处ITER EX代表算一个单点,能量分解要算好几个单点,你直接看最后一次出现的ITER EX。从后往前搜即可

好的,谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 01:43 , Processed in 0.173780 second(s), 23 queries , Gzip On.

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