24小时热门版块排行榜    

查看: 561  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 东南大学364求调剂 +3 JasonYuiui 2026-03-15 3/150 2026-03-15 18:57 by 无际的草原
[考研] 321求调剂 +3 大米饭! 2026-03-15 3/150 2026-03-15 17:48 by 哈哈哈哈嘿嘿嘿
[考研] 290求调剂 +3 孔志浩 2026-03-12 8/400 2026-03-15 15:30 by 孔志浩
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 311求调剂 +3 26研0 2026-03-15 3/150 2026-03-15 09:12 by JourneyLucky
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +3 Losir 2026-03-12 3/150 2026-03-14 12:11 by 热情沙漠
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 求调剂,一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 材料专硕274一志愿陕西师范大学求调剂 +4 薛云鹏 2026-03-13 4/200 2026-03-13 10:40 by 学员8dgXkO
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 0857 资源与环境 285分 +6 未名考生 2026-03-09 6/300 2026-03-11 21:08 by 30660438
信息提示
请填处理意见