计算化学公社

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

[综合交流] Gaussian-ORCA能量对比疑惑

[复制链接 Copy URL]

239

帖子

0

威望

2333

eV
积分
2572

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 柒月小鱼 于 2018-6-26 15:49 编辑

各位老师好!
我在计算一个绝热能的时候(T1-S0)开始用的Gaussian优化计算频率,而后用ORCA计算单点,发现二者差距达到5kcal,
而且Gaussian的计算结果明显更与实验值接近,是我的ORCA输入有问题么(我计算的是金属配合物)
还请各位老师指教一二
ORCA关键词
!RI-PWPB95 RIJK D3BJ def2-TZVP def2/JK Grid5 tightSCF nopop
%basis NewECP Ir "SDD" end  
Gaussian关键词
#p freq nosymm b3lyp/genecp scf=(maxcyc=400,conver=8) EmpiricalDispersion=GD3BJ
同一原子坐标计算

result.png (7.79 KB, 下载次数 Times of downloads: 81)

result.png

482

帖子

10

威望

6934

eV
积分
7616

Level 6 (一方通行)

BSJ Institute

2#
发表于 Post on 2018-6-26 18:15:43 | 只看该作者 Only view this author
两个软件用的泛函和基组都不一样当然没有可比性了

239

帖子

0

威望

2333

eV
积分
2572

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2018-6-26 19:30:17 | 只看该作者 Only view this author
Accelerator 发表于 2018-6-26 18:15
两个软件用的泛函和基组都不一样当然没有可比性了

是的
可是我纠结的是1.低级别的计算数据与实验更相近
2.在想我的ORCA是不是还是有些地方没有弄好
我重复文献的计算,用的是我Excel第三栏的基组,得到的结果与我高斯计算结果几乎一致

433

帖子

8

威望

4317

eV
积分
4910

Level 6 (一方通行)

石墨

4#
发表于 Post on 2018-6-26 22:56:38 | 只看该作者 Only view this author
柒月小鱼 发表于 2018-6-26 19:30
是的
可是我纠结的是1.低级别的计算数据与实验更相近
2.在想我的ORCA是不是还是有些地方没有弄好

泛函之间没有固定的谁高谁低, 要看是否适合体系
比较泛函的适用性时应该固定基组, 比较基组的精度时也应该固定泛函
自在飞花轻似梦,无边丝雨细如愁。

全自动反应动力学(ReaxFF、AIMD、NEP等)后处理工具网页版:http://cc-portal.xyz/reax_tools

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

5#
发表于 Post on 2018-6-26 23:06:15 | 只看该作者 Only view this author
问这个问题说明你对量化计算的基本思想还很糊涂

要么对比程序实现的差异,要么对比理论方法的差异,要么对比基组的差异,要么考察影响数值精度的设定的差异,要么对比理论方法对实验结果重现精度的差异...比其中一个就应当令其它因素保持不变。
你当前问题把乱七八糟因素全搅合在一起,根本什么都没法讨论
北京科音自然科学研究中心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

239

帖子

0

威望

2333

eV
积分
2572

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2018-6-27 09:08:09 | 只看该作者 Only view this author
Graphite 发表于 2018-6-26 22:56
泛函之间没有固定的谁高谁低, 要看是否适合体系
比较泛函的适用性时应该固定基组, 比较基组的精度时也应 ...

是的,泛函之间没有固定的高低短长,这也是入门就该了解的
可以最近看到各种说双杂化来做单点计算的效果很好,所以想来试试
结果高级别的泛函基组与实验相差较远(这是做方法的一个分子,实验结果应该是没有问题的,多次被验证了)
当然文献中的能很好的重复实验结果,怀疑我的ORCA使用ing

239

帖子

0

威望

2333

eV
积分
2572

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2018-6-27 09:12:49 | 只看该作者 Only view this author
sobereva 发表于 2018-6-26 23:06
问这个问题说明你对量化计算的基本思想还很糊涂

要么对比程序实现的差异,要么对比理论方法的差异,要么 ...

我是这样考虑的,
做高级别的计算,按理是比低级别的计算更理想,(尤其是实验已经非常成功的东西)
我没有去做程序对比也就是这个原因,考虑PWPB95是目前的金方法
\(^o^)/~

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

8#
发表于 Post on 2018-6-27 14:19:17 | 只看该作者 Only view this author
柒月小鱼 发表于 2018-6-27 09:12
我是这样考虑的,
做高级别的计算,按理是比低级别的计算更理想,(尤其是实验已经非常成功的东西)
我 ...


那你就没必要去和Gaussian去对比,都用ORCA就完了
根本没必要用SDD,def2-TZVP本身对Ir就有定义,何况SDD本身比def2-TZVP还小,而过渡金属在整个体系里又占有最重要地位,刻意用SDD明显有害无益

倘若计算流程都合理,数据有可比性,B3LYP/6-31G**真的结果比PWPB95-D3BJ/def2-TZVP还好,那9成可能也只是对当前体系的巧合。
北京科音自然科学研究中心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

239

帖子

0

威望

2333

eV
积分
2572

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2018-6-29 09:21:47 | 只看该作者 Only view this author
sobereva 发表于 2018-6-27 14:19
那你就没必要去和Gaussian去对比,都用ORCA就完了
根本没必要用SDD,def2-TZVP本身对Ir就有定义,何况 ...

老师早上好!
我现在将我的输出文件打包传上来了,希望您可以抽空看看,嘿嘿
其中ORCA使用的nohup命令,只有一个.out。第一个单点是单重态的,第二个单点是三重态的
麻烦老师了

ppy.rar

528.46 KB, 下载次数 Times of downloads: 17

输出文件

239

帖子

0

威望

2333

eV
积分
2572

Level 5 (御坂)

10#
 楼主 Author| 发表于 Post on 2018-7-3 16:37:33 | 只看该作者 Only view this author
sobereva 发表于 2018-6-27 14:19
那你就没必要去和Gaussian去对比,都用ORCA就完了
根本没必要用SDD,def2-TZVP本身对Ir就有定义,何况 ...

老师我在def2-QZVPP级别下计算了单点 PWPB95-D3BJ和B2PLYP-D3BJ
两者对比试验的效果确实不如用ORCA/B3LYP-D3BJ/def2-TZVP
那么我的单点计算是否使用高斯在B3LYP-D3BJ/def2-TZVP计算比较好,免得别人问在同种泛函还要换软件,虽然ORCA真的很快
测试了几天
还请老师指教

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

11#
发表于 Post on 2018-7-4 00:59:42 | 只看该作者 Only view this author
柒月小鱼 发表于 2018-7-3 16:37
老师我在def2-QZVPP级别下计算了单点 PWPB95-D3BJ和B2PLYP-D3BJ
两者对比试验的效果确实不如用ORCA/B3LY ...

既然结果是这样,就用这个吧
北京科音自然科学研究中心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

239

帖子

0

威望

2333

eV
积分
2572

Level 5 (御坂)

12#
 楼主 Author| 发表于 Post on 2018-7-4 09:42:46 | 只看该作者 Only view this author
sobereva 发表于 2018-7-4 00:59
既然结果是这样,就用这个吧

好的 谢谢老师耐心指导

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 22:21 , Processed in 0.201015 second(s), 23 queries , Gzip On.

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