计算化学公社

 找回密码 Forget password
 注册 Register
Views: 18323|回复 Reply: 3

[Gaussian/gview] Gaussian16里用PBE0关键词计算的实际上是PBE0-DH双杂化泛函

[复制链接 Copy URL]

4万

帖子

99

威望

4万

eV
积分
89975

管理员

公社社长+计算化学玩家

发表于 Post on 2019-6-25 11:01:25 | 显示全部楼层 Show all |阅读模式 Reading model
这个帖子目的是给Gaussian初学者们提个醒。

Gaussian里,知名的PBE0泛函应当写为PBE1PBE,这点手册里已经写明了。在G09里,直接写PBE0程序会报错,因为不认识
在Gaussian里有一个规则是只要不引起歧义,关键词可以最大程度简化,比如empiricaldispersion=GD3BJ可以简写为em=GD3BJ
从G16开始,支持了一个精度很差的双杂化泛函PBE0-DH。于是乎,根据关键词可简化规则,在G16里面写PBE0的话程序并不会报错,而给你算的是PBE0-DH。倘若你读取的是这个级别对应的电子能量那还好(比如读的是输出文件末尾MP2=后面的值,见《谈谈该从Gaussian输出文件中的什么地方读电子能量》http://sobereva.com/488),虽然精度不怎么样但还算勉强能用,而如果读取的是SCF Done或者末尾的HF=后面的值,那得到的根本什么都不是(实际上对应于PBE0-DH里杂化泛函部分,并不等于PBE0泛函的结果),毫无实际意义,完全不能用来发表。

所以,如果你之前用G16,取的是写了PBE0关键词得到的数据,并且当成是PBE0泛函的结果,需要把所有数据重算;如果数据已经发表了,对数据有责任心的话,应当写个correction。

Gaussian的这个关键词可简化规则在PBE0上很坑爹,Gaussian开发者应该将这种情况作为一个特例考虑,不允许只写PBE0。

为了让初学者更清楚关系,这里总结一下,对于G16(至少截止到B.01)而言:
PBE关键词=报错
PBEPBE关键词=PBE
PBE1PBE关键词=PBE0
PBE0关键词=PBE0-DH


评分 Rate

参与人数
Participants 15
eV +45 收起 理由
Reason
lurensan + 3 赞!
梵墨 + 1 谢谢
lao7 + 5 谢谢
rendong + 5 赞!
玉米猫 + 4 谢谢
量化新人 + 2 谢谢
w筱翊 + 1 谢谢
Winnie-001 + 3
逮捕小知奈 + 4 谢谢
水果皮儿 + 5 谢谢
cyx98 + 3 精品内容
Reid_Chan + 4
阿锋001 + 1 这就有点坑啦
liyuanhe211 + 1 WTF?!
Novice + 3 谢谢分享

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

1020

帖子

0

威望

3510

eV
积分
4530

Level 6 (一方通行)

发表于 Post on 2019-6-25 15:50:19 来自手机 | 显示全部楼层 Show all
G16有时SCF没收敛的计算也给Normal termination,曾被坑过,当时处理数据还以为有重大发现了。还有个bug就是,开壳层体系算NBO,如果输出chk文件,那么输出的只有alpha轨道信息

719

帖子

0

威望

1614

eV
积分
2333

Level 5 (御坂)

发表于 Post on 2022-4-17 02:32:44 | 显示全部楼层 Show all
刚刚犯错,特来看看!谢谢!

64

帖子

0

威望

778

eV
积分
842

Level 4 (黑子)

发表于 Post on 2022-5-1 10:05:14 | 显示全部楼层 Show all
同楼上,刚从坑里爬出来

本版积分规则 Credits rule

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

GMT+8, 2023-2-7 02:43 , Processed in 0.185526 second(s), 22 queries .

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