24小时热门版块排行榜    

查看: 563  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +3 Liwangman 2026-03-15 3/150 2026-03-15 18:16 by JourneyLucky
[基金申请] 国自科面上基金字体 +4 iwuli 2026-03-12 5/250 2026-03-15 17:07 by 风云无泪
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 085600材料与化工 求调剂 +11 enenenhui 2026-03-13 12/600 2026-03-14 21:48 by laoshidan
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 纺织、生物、化学、材料相关专业招生了 +4 耶耶业 2026-03-09 7/350 2026-03-12 19:05 by Equinoxhua
[考研] 0857 资源与环境 285分 +6 未名考生 2026-03-09 6/300 2026-03-11 21:08 by 30660438
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
信息提示
请填处理意见