第16届北京科音初级量子化学培训班将于10月3~6日于北京举行,是量子化学初学者迅速、正确上手计算的最佳机会!详情请点击此链接查看,欢迎参加!

计算化学公社

 找回密码
 现在注册!
查看: 6980|回复: 1

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

[复制链接]

3万

帖子

99

威望

3万

eV
积分
72834

管理员

公社社长+计算化学玩家

发表于 2019-6-25 11:01:25 | 显示全部楼层 |阅读模式
这个帖子目的是给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


评分

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

查看全部评分

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395(已满),2号:466017436(已满),3号:764390338(可加)。合计8000人,讨论范畴相同
思想家公社的门口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!

852

帖子

0

威望

2488

eV
积分
3340

Level 5 (御坂)

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

本版积分规则

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

GMT+8, 2021-9-25 23:03 , Processed in 0.775225 second(s), 25 queries .

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