| 查看: 561 | 回复: 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] |
回复此楼» 猜你喜欢
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有95人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有0人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
诚邀有志之士加入江苏大学环境学院合成生物学课题组！(长期接收学硕专硕、硕博连读)
已经有16人回复
-
一志愿南大理论与计算化学方向（070300）342分，求调剂
已经有1人回复
-
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
-
调剂招生
已经有7人回复
-
河南理工大学化工学院招收调剂生啦！
已经有1人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13046.7
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.1小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
3楼2008-09-22 22:32:13
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
- 专业: 理论和计算化学
- 管辖: 计算模拟
|
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
