24小时热门版块排行榜    

查看: 2408  |  回复: 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的回帖
信息提示
请填处理意见