24小时热门版块排行榜    

查看: 777  |  回复: 4
当前主题已经存档。

beefly

专家顾问 (职业作家)

地沟油冶炼专家

[交流] 【原创】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]牛肉一般地
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

★
小木虫(金币+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 ]
Simpler, stronger, more tolerant
2楼2010-01-05 06:46:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

★
小木虫(金币+0.5):给个红包,谢谢回帖交流
两位专家的讨论让人受益匪浅啊。
前途光明,出路难觅!
3楼2010-01-05 08:17:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

受教了!
4楼2010-01-05 10:07:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

引用回帖:
Originally posted by neweroica at 2010-1-5 06:46:


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

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

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

# hf/3-21g n ...

谢谢!按照你的解释,我试了试其它程序的NBO接口(pc-gamess,q-chem,orca),发现都有这个问题。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
5楼2010-01-05 12:04:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 beefly 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 WqtG0Aem8Ufl 2026-10-24 6/300 2026-10-26 08:26 by 13760108691
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 GVGRDgNfv7P4 2026-10-23 8/400 2026-10-26 08:07 by zjoO7LQNO3k6
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:18 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:12 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:05 by cLOojpOxkX0G
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 5/250 2026-10-26 01:46 by cLOojpOxkX0G
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 DI1qFHFK0KbE 2026-10-24 3/150 2026-10-26 01:37 by cLOojpOxkX0G
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 7xHjyB3bAiwg 2026-10-24 5/250 2026-10-26 01:20 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:54 by cLOojpOxkX0G
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:19 by cLOojpOxkX0G
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 GVGRDgNfv7P4 2026-10-23 5/250 2026-10-26 00:08 by cLOojpOxkX0G
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 3/150 2026-10-25 20:11 by cLOojpOxkX0G
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 3/150 2026-10-25 14:23 by XiYP04It3Evb
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 12:59 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:51 by XiYP04It3Evb
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 00:15 by XcOYG53d3x4p
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 18:23 by WqtG0Aem8Ufl
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 15:41 by Equinoxhua
[文学芳草园] 亲密关系 +7 宋昙 2026-10-21 13/650 2026-10-23 17:07 by 小牧11
[考博] 郑州大学申博 +5 stars~ 2026-10-19 8/400 2026-10-22 15:06 by stars~
信息提示
请填处理意见