| 查看: 566 | 回复: 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] |
回复此楼» 猜你喜欢
河南理工大学化工学院招收调剂生啦！
已经有1人回复
-
国家级人才课题组招研究生 （催化方向）调剂
已经有13人回复
-
物理化学论文润色/翻译怎么收费?
已经有225人回复
-
[天津]中国民航大学 招收材料类研究生
已经有0人回复
-
国家大人才课题组招收化学、材料化工类2026年调剂生
已经有0人回复
-
浙师大化与材学院2026年有大量硕士研究生调剂名额 （线上 调剂）
已经有9人回复
-
吉林工程技术师范学院 材料与化工 （0856） 招收调剂
已经有5人回复
-
国家大人才课题组招收化学、材料化工类2026年调剂生
已经有0人回复
-
国家大人才课题组招收化学、材料化工类2026年调剂生
已经有0人回复
-
吉林工程技术师范学院 交叉学科研究院（08/07） 招收调剂
已经有3人回复
-
国家大人才课题组招收化学、材料化工类2026年调剂生
已经有0人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13046.7
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.1小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
6楼2008-10-05 10:19:23
wxjbuilder
木虫 (知名作家)
- 应助: 1 (幼儿园)
- 贵宾: 0.05
- 金币: 2986.4
- 散金: 1761
- 红花: 2
- 沙发: 1
- 帖子: 5392
- 在线: 273.2小时
- 虫号: 442356
- 注册: 2007-10-27
- 性别: GG
- 专业: 理论和计算化学
2楼2008-09-22 20:52:38
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13046.7
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.1小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
3楼2008-09-22 22:32:13
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13046.7
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.1小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
|
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 ] |
4楼2008-10-05 07:41:17
