| 查看: 636 | 回复: 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] |
回复此楼» 猜你喜欢
安徽大学物质科学与信息技术研究院罗其全教授组2026年化学博士招生（计算催化组）
已经有1人回复
-
现招2026北京航空航天大学联培博士1名（北航学位）-材料AI方向
已经有3人回复
-
物理化学论文润色/翻译怎么收费?
已经有162人回复
-
普瓦捷环境与材料化学研究所 (IC2MP)博后招聘
已经有0人回复
-
香港科技大学(广州)方超课题组招聘全奖博士
已经有0人回复
-
宁夏大学化学工程与技术专业招收博士研究生
已经有1人回复
-
分子动力学 Plumed教程求助
已经有0人回复
-
南京大学软物质计算模拟方向硕博连读/直博生招生
已经有0人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
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
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13249.5
- 散金: 27307
- 红花: 301
- 帖子: 4678
- 在线: 928.4小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
3楼2008-09-22 22:32:13
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13249.5
- 散金: 27307
- 红花: 301
- 帖子: 4678
- 在线: 928.4小时
- 虫号: 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
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
5楼2008-10-05 08:54:20
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13249.5
- 散金: 27307
- 红花: 301
- 帖子: 4678
- 在线: 928.4小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
6楼2008-10-05 10:19:23
100