24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 564  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

[交流] Gaussian03的O3LYP是错的

截至目前,所有用Gaussian03作的O3LYP计算,结果都是错的。DFT专家Truhlar的小组早在04年就发现Gaussian 03的O3LYP有问题,无法重复O3LYP文献里的标准结果。见

Phys. Chem. Chem. Phys., 2004, 6, 673 - 676,
Tests of second-generation and third-generation density functionals for thermochemical kinetics
Yan Zhao, Jingzhi Pu, Benjamin J. Lynch and Donald G. Truhlar

不过文章没给出解决办法。通过分析源代码,发现O3LYP的这个BUG隐藏得很深。OPTX和Slater交换泛函最初的因子是正确的,但不知在什么位置发生了改变,因此需要用IOP指令重新覆盖OPTX和Slater的因子:

# o3lyp IOP(3/77=0813309262)

改后的O3LYP可以给出与其他程序一致的结果。测试结果见

http://bbs.zixia.net/disparticle ... ;ID=36527&pos=7
[search]O3LYP[/search]
回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

IOP(3/77=0813309262)可以用于gaussian03.b.03,c.01,c.02(估计其他的b.0x版本也都可以),但是对于d.01无效。对于gaussian03的其他版本,可以运行这个F离子的测试:

# cc-pvtz IOp(3/74=-24) IOP(3/77=0813309262)

-1 1
F

如果计算的能量在-99.6980... Hartree左右,说明指定IOP(3/77=0813309262)后的O3LYP是对的。

O3LYP的结构优化+频率计算还有个问题。opt和freq要分开计算,写成两个输入文件,或定义两步任务。要是写在一起,由于gaussian的另一个bug,导致输入的IOP(3/77=0813309262)无法传递到freq计算中,因此在频率步骤中计算的能量、梯度、力常数、偶极矩等,都是错的。

[ Last edited by beefly on 2008-10-5 at 10:57 ]
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2008-10-05 07:41:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

wxjbuilder

木虫 (知名作家)

多谢指教
能不能介绍一下O3LYP与B3LYP有啥区别?优缺点?
悠悠天地,奈何独立苍茫?
2楼2008-09-22 20:52:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

主要差别有两处。一是交换泛函不同,OPTX是一种改良的Beck 88泛函,从理论上讲要更准确(从实用角度可能未必如此)。二是HF精确交换的因子不同,几乎是1:2的关系。对主族元素构成的体系可能区别不大,对于过渡元素,特别是一些电子相关能比较强的3d金属,计算结果对HF因子有很强的依赖性,O3LYP和B3LYP的结果有可能完全不同。

其他交换、相关项以及不同的因子,造成的差别通常不大。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2008-09-22 22:32:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


不妨将此问题投到CCL上讨论
5楼2008-10-05 08:54:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本科郑州大学,一志愿华东师范大学282求调剂 +14 熊哥xtk 2026-04-07 17/850 2026-04-07 15:08 by WTUChen
[考研] 305分求调剂 +3 哈_哈_哈_哈_哈 2026-04-04 5/250 2026-04-07 14:49 by 哈_哈_哈_哈_哈
[考研] 0854求调剂 +9 亨氏番茄沙司 2026-04-06 10/500 2026-04-07 14:37 by shdgaomin
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +12 哇呼哼呼哼 2026-04-01 13/650 2026-04-07 10:02 by zhen~
[考研] 22408 调剂材料 +6 我叫ez 2026-04-06 7/350 2026-04-07 09:04 by MW_WTU
[考研] 304求调剂 +7 素年祭语 2026-04-06 10/500 2026-04-06 23:31 by 素年祭语
[考研] 285求调剂 +15 哦呦呼o 2026-04-04 17/850 2026-04-06 23:02 by chenzhimin
[考研] 华南理工0703化学,总分336求调剂 +10 手机用户 2026-04-03 11/550 2026-04-06 22:20 by michael2011
[考研] 一志愿生物与医药,296分,求调剂 +11 66鹿 2026-04-03 13/650 2026-04-06 21:45 by 52305043001
[考研] 283求调剂 +5 baiiyu 2026-04-05 6/300 2026-04-05 20:35 by 啵啵啵0119
[考研] 308求调剂 +4 maverick^_^ 2026-04-03 4/200 2026-04-05 19:08 by 蓝云思雨
[考研] 一志愿西北农林畜牧专硕336分求调剂 +3 5ourr 2026-04-03 3/150 2026-04-05 10:40 by JOKER0401
[考研] 求生物学学硕调剂——364分 +7 云朵遛弯指南 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
[考研] 286求调剂 +3 草木不言 2026-04-04 3/150 2026-04-04 22:40 by lbsjt
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 22:08 by 啵啵啵0119
[考研] 22408,264求调剂 +3 ywh729 2026-04-03 4/200 2026-04-04 11:04 by ywh729
[考研] 学硕288调剂!!! +3 小王xw123 2026-04-03 3/150 2026-04-03 21:20 by 啵啵啵0119
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 08工科求调剂290分 +5 1314捧花 2026-04-02 8/400 2026-04-02 13:16 by 乔哒哒哒
[考研] 生物学327,求调剂 +5 书上的梅子 2026-04-01 6/300 2026-04-02 06:47 by ilovexiaobin
信息提示
请填处理意见