| 查看: 1231 | 回复: 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 ] |
回复此楼» 猜你喜欢
国家级人才课题组招收2026年化学，材料化工 调剂生
已经有0人回复
-
河南理工大学招收材料，化工，化学调剂生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有105人回复
-
国家级人才课题组招收2026年化学，材料化工 调剂生
已经有0人回复
-
国家级人才团体课题组招收2026届化学类学硕，专硕。
已经有1人回复
-
河南理工大学绿色催化课题组招收材料，化工，化学等方向研究生
已经有0人回复
-
河南理工大学化工学院招收调剂生啦！
已经有1人回复
-
国家级人才课题组招研究生 （催化方向）调剂
已经有13人回复
-
[天津]中国民航大学 招收材料类研究生
已经有0人回复
-
国家大人才课题组招收化学、材料化工类2026年调剂生
已经有0人回复
-
浙师大化与材学院2026年有大量硕士研究生调剂名额 （线上 调剂）
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
7楼2009-11-15 07:15:20
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
