24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 558  |  回复: 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的回帖

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的回帖

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的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


不妨将此问题投到CCL上讨论
5楼2008-10-05 08:54:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

我以前的单位买了gaussian的site license,我有帐号,能用,但不知道是否合法。现在的单位是上了gaussian的黑名单的。而且测试的gaussian03版本太多了。去ccl讨论,弄不好会给自己惹麻烦。

[ Last edited by beefly on 2008-10-5 at 10:23 ]
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
6楼2008-10-05 10:19:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 beefly 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂 +11 A child 2026-04-04 11/550 2026-04-06 23:05 by chenzhimin
[考研] 考研调剂一志愿北矿安全科学与工程学硕316 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
[考研] 求助071001调剂!!! +4 黄守松 2026-04-05 5/250 2026-04-06 10:55 by 1028907439
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 22:23 by hemengdong
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 考研生物学考A区211,初试322,科目生化和生物综合,求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
[考研] 324求调剂 +9 想上学求调 2026-04-03 9/450 2026-04-04 23:57 by 果冻大王
[考研] 292分,材料与化工,申请调剂 +22 程晴之 2026-04-01 26/1300 2026-04-04 22:03 by hemengdong
[考研] 数二英二348求调剂 +4 hxdzj1 2026-04-03 5/250 2026-04-03 21:25 by zhq0425
[考研] 372分材料与化工(085600)一志愿湖南大学求调剂 +3 蓝笺片 2026-04-03 4/200 2026-04-03 17:58 by Jimmyandyou
[考研] 285求调剂 +7 AZMK 2026-04-02 9/450 2026-04-03 11:12 by wanwan00
[考研] 330求调剂 +3 白神呜呼呼 2026-04-02 3/150 2026-04-03 10:15 by 蓝云思雨
[考研] 一志愿大工学硕,求调剂 +4 yub0811 2026-04-02 4/200 2026-04-02 21:36 by 百灵童888
[考研] 298求B区调剂 +4 zzz,,r 2026-04-02 5/250 2026-04-02 12:17 by 土木硕士招生
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +9 双马尾痞老板2 2026-04-01 9/450 2026-04-02 12:14 by oooqiao
[考研] 0710生物学求调剂 +9 manman511 2026-04-01 9/450 2026-04-02 10:00 by zxl830724
[考研] 【求调剂】新能源材料本科,一志愿211,初试321 +6 求调剂学校, 2026-04-02 6/300 2026-04-02 09:41 by 晴空210210
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 一志愿 南京航空航天大学 ,080500材料科学与工程学硕 +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
信息提示
请填处理意见