计算化学公社

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

[GAMESS-US] Gaussian 与Gamess 的能量相差10倍

[复制链接 Copy URL]

87

帖子

0

威望

380

eV
积分
467

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
大家好,我是用Gaussian09进行优化,然后用Gamess进行能量分解,但是却发现高斯的相互作用能,要比Gamess的高10倍,我用的是Mp2方法,aug-cc-pVDZ基组,构型优化和能量分解都没有问题,但是就是能量二者相差10倍,不知道是什么原因,希望大家能给予解答,万分感谢!!

446

帖子

0

威望

7797

eV
积分
8243

Level 6 (一方通行)

2#
发表于 Post on 2015-9-29 09:39:30 | 只看该作者 Only view this author
检查:1.计算能量的分子构型,单体重新优化了?抑或用双体里相应的构型?
2.BSSE
3.full和frozencore是否一致
4.自洽场收敛情况

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

87

帖子

0

威望

380

eV
积分
467

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2015-9-29 09:44:05 | 只看该作者 Only view this author
yflchx 发表于 2015-9-29 09:39
检查:1.计算能量的分子构型,单体重新优化了?抑或用双体里相应的构型?
2.BSSE
3.full和frozencore是否 ...

1.单体的能量都是用二体里的构型计算的;
2.在优化二体的时候,都加了counter=2;
3.因为是新手,不知道怎么看。。。
4.单体、二体优化的时候都是正常结束,没有虚频

446

帖子

0

威望

7797

eV
积分
8243

Level 6 (一方通行)

4#
发表于 Post on 2015-9-29 09:52:19 | 只看该作者 Only view this author
xiaobogaga 发表于 2015-9-29 09:44
1.单体的能量都是用二体里的构型计算的;
2.在优化二体的时候,都加了counter=2;
3.因为是新手,不知 ...

3.是指是不是计算能量的时候一个软件用了MP2(fc),另一个软件用了MP2(full)?

4.看看SCF自洽场能量迭代的情况,对比一下两个程序。

87

帖子

0

威望

380

eV
积分
467

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2015-9-29 17:01:10 | 只看该作者 Only view this author
yflchx 发表于 2015-9-29 09:52
3.是指是不是计算能量的时候一个软件用了MP2(fc),另一个软件用了MP2(full)?

4.看看SCF自洽场能量迭 ...

3. 计算单体和二体的时候命令行只写了 MP2/基组 怎么看是(fc) 还是(full) 呢

446

帖子

0

威望

7797

eV
积分
8243

Level 6 (一方通行)

6#
发表于 Post on 2015-9-29 17:55:38 | 只看该作者 Only view this author
xiaobogaga 发表于 2015-9-29 17:01
3. 计算单体和二体的时候命令行只写了 MP2/基组 怎么看是(fc) 还是(full) 呢

默认是FC,frozencore。这个问题可以跳过了。

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

7#
发表于 Post on 2015-9-29 23:28:40 | 只看该作者 Only view this author
碰见问题先把问题尽可能简化,一点一点对比来分析原因。
比如考察相互作用能的计算差异,你就先分析单点能的计算差异,别把其它乱七八糟的因素引入,否则肯定弄不清楚。
你就用相同的结构,在MP2/aug-cc-pVDZ下用高斯和GAMESS-US各算个单点。注意都先别用冻核(两个程序在这点上有差异),看看HF能量和MP2相关能是否都能对上。假设没对上,然后你再看比如是否是实际计算用的基函数数目不对(比如高斯会自动去除线性依赖的基函数)。
慢慢一点一点对比总能搞清楚。
北京科音自然科学研究中心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

87

帖子

0

威望

380

eV
积分
467

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2015-9-30 13:52:18 | 只看该作者 Only view this author
sobereva 发表于 2015-9-29 23:28
碰见问题先把问题尽可能简化,一点一点对比来分析原因。
比如考察相互作用能的计算差异,你就先分析单点能 ...

非常感谢sob老师的耐心解答和指导,您的话也为我以后的科研指点了方向,我再试试,谢谢您

87

帖子

0

威望

380

eV
积分
467

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2015-9-30 16:46:25 | 只看该作者 Only view this author
sobereva 发表于 2015-9-29 23:28
碰见问题先把问题尽可能简化,一点一点对比来分析原因。
比如考察相互作用能的计算差异,你就先分析单点能 ...

老师你好,我用Gassian和Gamess 分别计算了能量为什么会差这么多,附件中是输入和输出文件

MP2-ptCl2(NH3)2-ClF-sp.gjf

5.76 KB, 下载次数 Times of downloads: 13

高斯输入文件

MP2-ptCl2(NH3)2-ClF-sp.out

82.61 KB, 下载次数 Times of downloads: 5

高斯输出文件

1.inp

6.2 KB, 下载次数 Times of downloads: 22

Gamesss输入文件

1.out

663.8 KB, 下载次数 Times of downloads: 10

Gamess输出文件

44

帖子

0

威望

663

eV
积分
707

Level 4 (黑子)

10#
发表于 Post on 2015-10-1 02:44:04 | 只看该作者 Only view this author
基组输入错误。两个计算的基函数数目不同。

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

11#
发表于 Post on 2015-10-1 04:25:35 | 只看该作者 Only view this author
xiaobogaga 发表于 2015-9-30 16:46
老师你好,我用Gassian和Gamess 分别计算了能量为什么会差这么多,附件中是输入和输出文件

完全没有可比性嘛
1 一个是单点一个是counterpoise任务
2 一个对配体用cc-pVDZ一个用aug-cc-pVDZ
3 一个用冻核一个没用冻核
北京科音自然科学研究中心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

87

帖子

0

威望

380

eV
积分
467

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2015-10-1 18:20:58 | 只看该作者 Only view this author
jjjspring 发表于 2015-10-1 02:44
基组输入错误。两个计算的基函数数目不同。

你好,我是新手,请问正确的输入应该怎么写,比如我用Gamess 计算的时候用MP2(full),在Gamess的输入文件里应该怎么写关键词呢?还有在Gamess里accd 表示 aug-cc-pvdz 机组吗?

87

帖子

0

威望

380

eV
积分
467

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2015-10-1 18:22:26 | 只看该作者 Only view this author
sobereva 发表于 2015-10-1 04:25
完全没有可比性嘛
1 一个是单点一个是counterpoise任务
2 一个对配体用cc-pVDZ一个用aug-cc-pVDZ

那sob 老师,如果gaussian 命令行里只写方法/基组 是不是就是计算单点能了?

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

14#
发表于 Post on 2015-10-1 22:31:33 | 只看该作者 Only view this author
xiaobogaga 发表于 2015-10-1 18:22
那sob 老师,如果gaussian 命令行里只写方法/基组 是不是就是计算单点能了?

北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

15#
发表于 Post on 2015-10-1 22:42:59 | 只看该作者 Only view this author
xiaobogaga 发表于 2015-10-1 18:20
你好,我是新手,请问正确的输入应该怎么写,比如我用Gamess 计算的时候用MP2(full),在Gamess的输入文件 ...

ACCD就是aug-cc-pVDZ
不想冻核就写
$MP2 NACORE=0 $END
北京科音自然科学研究中心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, 2025-8-14 07:45 , Processed in 0.171118 second(s), 26 queries , Gzip On.

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