| 查看: 562 | 回复: 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] |
回复此楼» 猜你喜欢
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
汕头大学化学招调剂，相关内容如下。
已经有4人回复
-
物理化学论文润色/翻译怎么收费?
已经有71人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
-
分子动力学模拟合作需求
已经有3人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
5楼2008-10-05 08:54:20
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
