24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2519  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

klaus1987

铜虫 (小有名气)

[求助] Multiwfn读取NBO的.40文件是对应MO吗?

如题,manual里说NBO的.40文件是NLMO,我听说这种轨道的性质和MO是相同的,请问这种说法是否正确?
另外,我试着做了.40的电子密度格点文件,又对同一结构(相同坐标)再次用高斯计算能量,得到MO轨道的fch文件,同样做电子密度格点。两个格点文件的质量选择相同,二者相减会得到差值,也就是说.40和MO电子密度格点不一样。这是正常的吗?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


★ ★ ★
gkf高: 金币+3 2013-06-17 11:17:18
以下是我开发Multiwfn时针对这个问题写的笔记,涉及到一些细节,可以参考

“在计算密度等属性时,按照现有方法,将NBO占据数乘上NBO轨道波函数平方,再加和,得到的并不是真正的总密度。因为以NBO为基时密度矩阵并不是对角的,因为获得NBO时并不是直接对整体密度矩阵对角化,而是依次将各个块密度矩阵对角化。为了得到正确的密度,必须输出以NBO为基的密度矩阵,将每个矩阵元乘上对应的两个NBO在此处的波函数值,然后整体加和,也就是必须把非对角元考虑进去。对于NAO也是如此。不过,按Multiwfn现有方法计算,对于NBO、NAO、NHO,密度积分倒是和期望的一致。因为它们都是相互正交的。使用NBO算的实空间函数和真实值偏差不大。
NLMO和普通方式获得的定域化轨道一样,以其为基时整体密度矩阵也是对角化的,因此可以直接将密度乘上波函数平方加和得到真实密度。正如NBCP文章里式子所述。
至于前头带P的就更不能这么算了,连正交都不满足,因此积分值都不对。”

» 本帖已获得的红花(最新10朵)

10楼2013-06-16 19:01:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+3 2013-06-13 21:05:02
klaus1987: 金币+5, ★★★很有帮助 2013-06-14 09:54:20
.40是MO,.39才是NLMO(手册写错了)
NLMO和MO截然不同,前者可认为是对后者酉变换产生的定域化轨道
.40和fch里的MO应当是一样的,存在差值可能是微小的数值精度问题所致,或者什么其它地方没有完全对应上。最好直接通过在主功能0观看等值面来比较二者。
2楼2013-06-13 18:10:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

klaus1987

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by sobereva at 2013-06-13 18:10:14
.40是MO,.39才是NLMO(手册写错了)
NLMO和MO截然不同,前者可认为是对后者酉变换产生的定域化轨道
.40和fch里的MO应当是一样的,存在差值可能是微小的数值精度问题所致,或者什么其它地方没有完全对应上。最好直 ...

我用格点文件作差得到的差值还是比较明显的啊,待我再去重复一次比较40和MO
3楼2013-06-14 09:59:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

klaus1987

铜虫 (小有名气)

我又详细的比对了两次计算的数据,算nbo.40和mo的原子坐标是完全一致的,格点密度都是中,结果是二者的电子密度确实完全不一样。mo的格点正值乘差别元(differential element)为161.94,而nbo.40的确是383.92。
4楼2013-06-14 10:49:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286分调剂 +3 Faune 2026-03-30 3/150 2026-03-30 12:42 by wxiongid
[考研] 317分 一志愿南理工材料工程 本科湖工大 求调剂 +11 芋泥小铃铛 2026-03-28 11/550 2026-03-30 12:23 by Wang200018
[考研] 291求调剂 +7 Y-cap 2026-03-29 8/400 2026-03-29 19:53 by klasasda
[考研] 环境工程 085701,267求调剂 +6 minht 2026-03-29 6/300 2026-03-29 16:21 by 学员8dgXkO
[考研] 070305高分子化学与物理 304分求调剂 +12 c297914 2026-03-28 12/600 2026-03-29 16:21 by Serene1974
[考研] 求收留 +5 1943443204 2026-03-28 5/250 2026-03-29 12:11 by 无际的草原
[考研] 340求调剂 +6 Amber00 2026-03-26 6/300 2026-03-29 12:06 by 无际的草原
[考研] 315求调剂 +4 akie... 2026-03-28 5/250 2026-03-28 21:05 by zhq0425
[考研] 0856,材料与化工321分求调剂 +12 大馋小子 2026-03-27 13/650 2026-03-28 10:56 by self2008
[考研] 085701环境工程求调剂 +9 多久上课 2026-03-27 9/450 2026-03-28 03:58 by fmesaito
[考研] 266求调剂 +11 阳阳哇塞 2026-03-27 12/600 2026-03-27 17:56 by yu221
[考研] 308求调剂 +7 墨墨漠 2026-03-25 7/350 2026-03-27 14:47 by 狂炫麦当当
[考研] 08开头275求调剂 +4 拉谁不重要 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
[考研] 284求调剂 +11 junqihahaha 2026-03-26 12/600 2026-03-27 04:37 by wxiongid
[考研] 342求调剂 +3 加油a李zs 2026-03-26 3/150 2026-03-27 00:29 by wxiongid
[考研] 321求调剂 +6 wasdssaa 2026-03-26 6/300 2026-03-26 20:57 by sanrepian
[考研] 中国科学院深圳先进技术研究院-光纤传感课题组招生-中国科学院大学、深圳理工大学联培 +5 YangTyu1 2026-03-26 5/250 2026-03-26 18:27 by 猫咪猫咪呀
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
信息提示
请填处理意见