| 查看: 1233 | 回复: 7 | |||||
| 当前主题已经存档。 | |||||
beefly专家顾问 (职业作家)
地沟油冶炼专家
|
[交流]
【beefly个人文集】总结：Gaussian结果无法用其它程序重复的原因 已有1人参与
|
||||
|
周期体系的量子化学计算结果，通常无法用其他程序重复，因为倒空间积分、库仑势等的算法存在差别。但是对于孤立体系的计算（原子、分子、团簇），不同的程序中的同样方法原则上应该得到相同的结果（在计算精度内），当然前提是分子结构相同。但也有一些例外。对于Gaussian而言，可能的原因有： 1. 使用的基组不同。即使是同名基组，因发表年代不同也会略有差别，因此要把基组打印出来，仔细检查基组是否完全相同。 2. 直角基函数和球谐基函数。不同程序的默认设置不同，结果也是不一样的。Gaussian中的默认设置随基组类型而定（见基组说明部分），或者在输入中直接定义。 3. 对于DFT计算，需要考虑积分格点的差别。不过在一般情况下，对总能量的影响不会非常大。 4. 对于杂化泛函，还有个更重要的原因：相名不同义。 最著名的是B3LYP：在局域关联部分，Gaussian用VWN-3而不是原始文献建议的VWN-5，幸好大多数密度泛函程序同时提供两个版本的B3LYP。 最近发现的是O3LYP，在文献中，OPTX项的因子定义不明确导致混乱：在Gaussian中采用的因子是0.8133*1.43169（据称计算结果和O3LYP原始文献一致；但也有反例，见PCCP，2004，6，673），而PC-GAMESS和ORCA程序采用0.8133（最新版本的PC-GAMESS也支持前一种定义）。此外，Jaguar和PQS程序的O3LYP结果都不同于以上两者，技术细节不清楚。 5. Post-HF计算也存在相名不同义的问题。如Gaussian的MR-MP2（CASSCF+MP2）方法，和其它程序的CASPT2、MR-MP2等方法的结果是不一样的。 6. Douglas-Kroll标量相对论计算。大多数程序使用点核模型，Gaussian使用更精确的Gaussian有限核模型。对于重元素体系的计算，总能量有显著差别。 7. 对于Post-HF计算，冻结轨道设置不同。Gaussian默认冻结全部芯轨道。如果HF参考能量相同，但Post-HF能量不同，可能就是这个原因。 8. 指定的或程序猜测的占据轨道不同。 9. 其它原因，如做了近似（如密度拟合，局域化轨道），程序BUG，输入错误，等。 [ Last edited by yjcmwgk on 2010-6-16 at 21:37 ] |
回复此楼» 猜你喜欢
材料科学与工程 或 材料与化工（专硕）调剂 要求：初试考408
已经有0人回复
-
298分 070300求调剂
已经有14人回复
-
物理化学论文润色/翻译怎么收费?
已经有158人回复
-
硕士研究生调剂
已经有0人回复
-
湖南大学张世国团队招26年第二批博士—离子液体、高分子与凝胶材料
已经有0人回复
-
湖南农业大学能源材料创新团队招收2025化学专业研究生（学硕）调剂
已经有0人回复
-
2026年西华大学材料学院-材料近净成形与表面技术研究团队-招收调剂研究生
已经有8人回复
-
调剂信息
已经有0人回复
-
河北农大生物与医药专业招收调剂，名额充足
已经有0人回复
-
齐鲁工业大学 国家重点实验室-新能源、锂离子电池、材料、化学方向
已经有0人回复
-
杰青/长江团队招收调剂
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
quantumor
金虫 (著名写手)
快乐兔子
- 应助: 13 (小学生)
- 金币: 8506.3
- 散金: 10
- 红花: 14
- 帖子: 1898
- 在线: 215.9小时
- 虫号: 630331
- 注册: 2008-10-19
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
学习了，不错。 2楼2009-11-13 09:56:53
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42654.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
3楼2009-11-13 14:15:44
4楼2009-11-13 21:17:19
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13043.7
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.1小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
5楼2009-11-14 00:22:30
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
|
比较了一下手上的程序包，Turbomole和Gaussian是一致的, molpro是不同的。 Turbomole的说明是（和Gaussian对应的文献JCP, 110,6158 (1999)相同) ------------------ density functional ------------------ PBE0 hybrid functional exchange: 3/4 (LDA+PBE) + 1/4 HF correlation: LDA (PW) + PBE spherical integration : Lebedev's spherical grid spherical gridsize : 5 i.e. gridpoints : 590 value for diffuse not defined radial integration : Chebyshev 2nd kind (scaling 3) radial gridsize : 8 integration cells : 1 partition function : becke partition sharpness : 3 ..... molpro手册的说明是PBE0 = 0.25d0 EXACT + 0.75d0 PBEX + PW91C 可能是相关泛函不同，不知道molpro为什么这么定义PBE0........ |
6楼2009-11-14 03:08:05
7楼2009-11-15 07:15:20
bluewhale
铁杆木虫 (正式写手)
- 应助: 58 (初中生)
- 金币: 8438.3
- 散金: 50
- 红花: 10
- 帖子: 890
- 在线: 362.2小时
- 虫号: 278289
- 注册: 2006-09-10
- 专业: 凝聚态物性I:结构、力学和
8楼2010-04-17 23:06:25
