24小时热门版块排行榜

返回列表

当前主题已经存档。

beefly

专家顾问 (职业作家)

地沟油冶炼专家

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

[交流] 【原创】Gaussian的NBO有BUG，会导致错误的MO bond order

做这样一个计算
# hf/3-21g nosymm pop=nboread

test

0 1
f
h 1 0.96

$nbo bndidx $end

得到MO bond order=1.1382

去掉nosymm，重复上面的计算，会得到错误结果MO bond order=-0.4741

这是由对称性错误导致的，使部分矩阵元的正负号发生了不应有的改变。

此问题出现在gaussian的各个版本中，包括产生的*.47文件。

回复此楼

» 猜你喜欢

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

1楼2010-01-05 05:02:27

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

neweroica

木虫 (著名写手)

应助: 53 (初中生)
金币: 4997.5
散金: 301
红花: 31
帖子: 2471
在线: 216.6小时
虫号: 36775
注册: 2004-02-11
专业: 理论和计算化学

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

引用回帖:

Originally posted by beefly at 2010-1-5 05:02:
做这样一个计算
# hf/3-21g nosymm pop=nboread

test

0 1
f
h 1 0.96

$nbo bndidx $end

得到MO bond order=1.1382

去掉nosymm，重复上面的计算，会得到错误结果MO bond order=-0.4741

这是由 ...

不是对称性的原因，而是计算时分子取向的原因。

即使给了NOSYMM照样可能得出错误结果（因为取决于你给定的分子取向）

试着计算下面两个输入文件：

# hf/3-21g nosymm pop=nboread

test

0 1
F 0. 0. 0.
H 0. 0. 0.96

$nbo bndidx $end

和

# hf/3-21g nosymm pop=nboread

test

0 1
F 0. 0. 0.
H 0. 0. -0.96

$nbo bndidx $end

第一个给出正确的MO键级，第二个则给出那个错误值。这两个输入文件的分子轴正好反向。

回到lz的例子，不加nosymm时，gaussian会在计算前先自动调整分子取向为“Standard orientation”，这里就是：
F 0.000000 0.000000 0.096000
H 0.000000 0.000000 -0.864000

这是因为，gaussian的标准取向是这样选取坐标原点的：(核电荷数*原子到坐标原点的距离)对所有原子的代数和为零。

因此，gaussian的标准取向正好是给出错误MO键级的取向，所以不加nosymm时，总会得到错误数值。

不管怎么样，这的确是个bug，值得注意。尤其是：即使加了nosymm也不能保证不会出错。

[ Last edited by neweroica on 2010-1-5 at 06:48 ]