| 查看: 554 | 回复: 4 | |||||
| 当前主题已经存档。 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
beefly专家顾问 (职业作家)
地沟油冶炼专家
|
[交流]
【beefly个人文集】Gaussian-NBO产生*.31、*.47文件的一个BUG
|
||||
|
如果基组含有g函数，那么Gaussian-NBO产生的*.31和*.47文件存在错误。以Ne原子的HF/cc-pvqz计算为例。计算完成后，打开NBO程序产生的*.47文件。在$CONTRACT部分： 第一个错误。没有对应于g函数的归一化收缩因子（CG）。仔细检查，发现打印到f函数里去了（位于CF）。 第二个错误。再找到NCOMP这行数组的最后一项，正确的话应当是9或15，分别对应球g函数或直角g函数的个数。若用原始的NBO 3.1，没有问题。如果用从源代码编译的NBO 5，这个数是7或10，真的变成f函数了。 同样，打开*.31文件。如果用NBO 3.1，那么最后一部分数据对应g函数的归一化收缩因子，你会发现都是0，g的因子错误地打印到f函数区（倒数第二部分）。如果用NBO 5，则根本没有g函数区。 因此，如果有g函数，*.31、*.47文件需要做些修改，才能用于GENNBO或用来画图。 我用pc-games和nwchem做上面的计算，一切正常。 最近读了一下NBO的源代码，发现NBO 3的$CONTRACT不支持g函数，但其它数据区可以。因此，在有g基函数的情况下，Gaussian做一般的NBO分析没问题。但是在产生.47或.31文件的时候，Gaussian的接口程序做了点小手脚，把g当作f函数输出给NBO。 接上NBO 5以后，虽然支持g函数了，但Gaussian的这段接口代码在编译的时候没改过来，所以仍然是按照f函数处理。 [ Last edited by yjcmwgk on 2010-6-16 at 21:35 ] |
回复此楼» 猜你喜欢
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有292人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有0人回复
-
诚邀有志之士加入江苏大学环境学院合成生物学课题组！(长期接收学硕专硕、硕博连读)
已经有16人回复
-
一志愿南大理论与计算化学方向（070300）342分，求调剂
已经有1人回复
-
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
-
调剂招生
已经有5人回复
-
河南理工大学化工学院招收调剂生啦！
已经有1人回复
lihb734
铁杆木虫 (职业作家)
站在计算化学入门的门槛上
- QC强帖: 2
- 应助: 132 (高中生)
- 金币: 7495.4
- 散金: 9378
- 红花: 74
- 帖子: 3107
- 在线: 1177.3小时
- 虫号: 448403
- 注册: 2007-11-01
- 性别: GG
- 专业: 理论和计算化学
3楼2010-01-09 11:18:55
tangsw911
铁杆木虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 10053.8
- 散金: 5
- 红花: 2
- 帖子: 715
- 在线: 585.7小时
- 虫号: 294160
- 注册: 2006-11-05
- 性别: GG
- 专业: 结构化学
2楼2010-01-09 09:42:25
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13043.7
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.1小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
4楼2010-01-10 14:34:07
5楼2010-01-11 03:33:06
