24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

专家经验: +459
QC强帖: 12
应助: 88 (初中生)
金币: 13381.5
散金: 27307
红花: 302
帖子: 4685
在线: 929小时
虫号: 408372
注册: 2007-06-21
性别: GG
专业: 理论和计算化学
管辖: 计算模拟

[交流] 【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人回复
固相反应中颗粒形貌变化计算模拟学习求助已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有185人回复
大湾区大学刘天辉课题组招收2027级博士生（与中山大学物理学院联合培养）已经有63人回复
电催化氧化技术处理高盐难降解废水具备独特优势已经有0人回复
电化学水处理四种主流阳极材料性能与工况适配性对比探讨已经有0人回复
深圳理工大学-深圳大学联合培养项目2027年接收推荐免试研究生已经有2人回复
铁离子电池求助已经有0人回复
想问为什么自制半电池循环只有50圈？已经有0人回复
硕士研究生全国统考英语二真题合集已经有0人回复

beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地

1楼 2010-01-09 04:21:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

fooo

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 6828
散金: 15
红花: 1
帖子: 318
在线: 872.6小时
虫号: 710228
注册: 2009-02-27
专业: 理论和计算化学
管辖: 量子化学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

厉害，佩服！

引用回帖:

Originally posted by beefly at 2010-1-9 04:21:
如果基组含有g函数，那么Gaussian-NBO产生的*.31和*.47文件存在错误。以Ne原子的HF/cc-pvqz计算为例。计算完成后，打开NBO程序产生的*.47文件。在$CONTRACT部分：

第一个错误。没有对应于g函数的归一化收缩因子 ...

回复此楼