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

计算化学公社

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

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

[复制链接]

2万

帖子

25

威望

2万

eV
积分
44998

管理员

公社社长

发表于 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


评分

参与人数 4eV +9 收起 理由
Reid_Chan + 4
阿锋001 + 1 这就有点坑啦
liyuanhe211 + 1 WTF?!
Novice + 3 谢谢分享

查看全部评分

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

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

501

帖子

0

威望

921

eV
积分
1422

Level 4 (黑子)

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

本版积分规则

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

GMT+8, 2019-8-26 18:24 , Processed in 0.157075 second(s), 24 queries .

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