| 查看: 684 | 回复: 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] |
回复此楼» 猜你喜欢
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Postdoc position, materials simulations/machine learning, U. of South Carolina
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有157人回复
-
祈福青基
已经有34人回复
-
牛来！米来！面来！
已经有8人回复
-
2027申博咨询
已经有11人回复
-
寻求CrC2Al的cif文件，感谢
已经有0人回复
-
ORCA 纯新手向入门指南0.7版
已经有0人回复
-
哈尔滨工业大学（深圳）绿色合成化学团队招收27年春季或秋季入学的博士生
已经有5人回复
-
【接测试】TPD-TPR、TPO、脉冲吸附、原位红外，电化学原位红外，站内联系
已经有0人回复
-
吉林大学化学学院理论化学研究所刘靖尧教授课题组招收2027级申请考核制博士研究生
已经有0人回复
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
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13383
- 散金: 27307
- 红花: 302
- 帖子: 4685
- 在线: 929小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
3楼2008-09-22 22:32:13
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13383
- 散金: 27307
- 红花: 302
- 帖子: 4685
- 在线: 929小时
- 虫号: 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
10