请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 394|回复: 11

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

[复制链接]

85

帖子

0

威望

739

eV
积分
824

Level 4 (黑子)

发表于 2018-6-26 15:47:58 | 显示全部楼层 |阅读模式
本帖最后由 柒月小鱼 于 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

103

帖子

0

威望

1110

eV
积分
1213

Level 4 (黑子)

发表于 2018-6-26 18:15:43 | 显示全部楼层
两个软件用的泛函和基组都不一样当然没有可比性了

85

帖子

0

威望

739

eV
积分
824

Level 4 (黑子)

 楼主| 发表于 2018-6-26 19:30:17 | 显示全部楼层
Accelerator 发表于 2018-6-26 18:15
两个软件用的泛函和基组都不一样当然没有可比性了

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

57

帖子

0

威望

911

eV
积分
968

Level 4 (黑子)

石墨

发表于 2018-6-26 22:56:38 | 显示全部楼层
柒月小鱼 发表于 2018-6-26 19:30
是的
可是我纠结的是1.低级别的计算数据与实验更相近
2.在想我的ORCA是不是还是有些地方没有弄好

泛函之间没有固定的谁高谁低, 要看是否适合体系
比较泛函的适用性时应该固定基组, 比较基组的精度时也应该固定泛函
长衫我亦何为者,也在游人笑语中.

1万

帖子

25

威望

1万

eV
积分
32364

管理员

公社社长

发表于 2018-6-26 23:06:15 | 显示全部楼层
问这个问题说明你对量化计算的基本思想还很糊涂

要么对比程序实现的差异,要么对比理论方法的差异,要么对比基组的差异,要么考察影响数值精度的设定的差异,要么对比理论方法对实验结果重现精度的差异...比其中一个就应当令其它因素保持不变。
你当前问题把乱七八糟因素全搅合在一起,根本什么都没法讨论
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

85

帖子

0

威望

739

eV
积分
824

Level 4 (黑子)

 楼主| 发表于 2018-6-27 09:08:09 | 显示全部楼层
Graphite 发表于 2018-6-26 22:56
泛函之间没有固定的谁高谁低, 要看是否适合体系
比较泛函的适用性时应该固定基组, 比较基组的精度时也应 ...

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

85

帖子

0

威望

739

eV
积分
824

Level 4 (黑子)

 楼主| 发表于 2018-6-27 09:12:49 | 显示全部楼层
sobereva 发表于 2018-6-26 23:06
问这个问题说明你对量化计算的基本思想还很糊涂

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

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

1万

帖子

25

威望

1万

eV
积分
32364

管理员

公社社长

发表于 2018-6-27 14:19:17 | 显示全部楼层
柒月小鱼 发表于 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程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

85

帖子

0

威望

739

eV
积分
824

Level 4 (黑子)

 楼主| 发表于 2018-6-29 09:21:47 | 显示全部楼层
sobereva 发表于 2018-6-27 14:19
那你就没必要去和Gaussian去对比,都用ORCA就完了
根本没必要用SDD,def2-TZVP本身对Ir就有定义,何况 ...

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

ppy.rar

528.46 KB, 下载次数: 4

输出文件

85

帖子

0

威望

739

eV
积分
824

Level 4 (黑子)

 楼主| 发表于 2018-7-3 16:37:33 | 显示全部楼层
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真的很快
测试了几天
还请老师指教

1万

帖子

25

威望

1万

eV
积分
32364

管理员

公社社长

发表于 2018-7-4 00:59:42 | 显示全部楼层
柒月小鱼 发表于 2018-7-3 16:37
老师我在def2-QZVPP级别下计算了单点 PWPB95-D3BJ和B2PLYP-D3BJ
两者对比试验的效果确实不如用ORCA/B3LY ...

既然结果是这样,就用这个吧
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

85

帖子

0

威望

739

eV
积分
824

Level 4 (黑子)

 楼主| 发表于 2018-7-4 09:42:46 | 显示全部楼层
sobereva 发表于 2018-7-4 00:59
既然结果是这样,就用这个吧

好的 谢谢老师耐心指导
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-7-20 07:25 , Processed in 0.212296 second(s), 26 queries .

快速回复 返回顶部 返回列表