| 查看: 2492 | 回复: 10 | |||
[求助]
Multiwfn读取NBO的.40文件是对应MO吗？
|
|
如题，manual里说NBO的.40文件是NLMO，我听说这种轨道的性质和MO是相同的，请问这种说法是否正确？ 另外，我试着做了.40的电子密度格点文件，又对同一结构(相同坐标)再次用高斯计算能量，得到MO轨道的fch文件，同样做电子密度格点。两个格点文件的质量选择相同，二者相减会得到差值，也就是说.40和MO电子密度格点不一样。这是正常的吗？ |
回复此楼» 猜你喜欢
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
求助立方晶体In(OH)3的CIF卡片
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有177人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有3人回复
-
阴离子交换膜电解二氧化碳还原
已经有0人回复
-
九江学院2026年最新高层次人才招聘公告
已经有0人回复
-
深圳信息职业技术学院2026年教师招聘公告（最新）
已经有0人回复
-
研究生招生
已经有0人回复
-
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 建议multiwfn开发GUI界面 文件打开什么的效率太低了
已经有10人回复
- 如何用Multiwfn软件做出P DOS?
已经有4人回复
- 使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响
已经有17人回复
- Multiwfn的初步学习，越详细越好。
已经有6人回复
- MultiWFN想要生成用于NBO分析的cub时不能载入.37文件
已经有3人回复
- 使用Multiwfn绘制NBO及相关轨道
已经有30人回复
- 用Multiwfn、Gaussview、Molekel、VMD观看龙虾、盆景、骨盆、大脑
已经有16人回复
- 使用Multiwfn做拓扑分析以及计算孤对电子角度
已经有18人回复
- 回答一些关于Multiwfn的疑问以及未来Multiwfn的发展打算
已经有37人回复
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
2楼2013-06-13 18:10:14
3楼2013-06-14 09:59:40
4楼2013-06-14 10:49:50
5楼2013-06-14 10:51:59
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
gmy1990: 金币+2 2013-06-14 16:53:14
klaus1987: 金币+5, ★★★★★最佳答案 2013-06-16 15:28:25
gmy1990: 金币+2 2013-06-14 16:53:14
klaus1987: 金币+5, ★★★★★最佳答案 2013-06-16 15:28:25
|
你可以去比较.40和fch当中特定的分子轨道的密度，但是不要去比较它们计算的总密度。 用.40作为输入后，选主功能6再选3，可以看到所有轨道都是占据的，也就是说，以.40为输入时无法正确判断轨道占据数，因此计算的总密度毫无意义。这是因为.40里面本身就没有提供轨道占据信息。计算分子性质的时候，不要用NBO plot文件，而应当直接用fch/wfn/wfx/molden作为输入。 只有以.31(NLMO)作为输入时，计算出的分子性质才和实际一致，这是手册中这句话的含义"Notice that among all types of the orbitals generated by NBO program, only using NLMO to calculate real space functions is meaningful!” |
6楼2013-06-14 14:11:24
7楼2013-06-16 15:29:56
8楼2013-06-16 17:18:12
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
★ ★ ★
gkf高: 金币+3 2013-06-17 11:17:08
gkf高: 金币+3 2013-06-17 11:17:08
|
这些.xx文件之间也不能相互比较计算出来的总的密度或其它实空间函数 原理上说，PNAO/NAO/PNHO/NHO/PNBO/NBO/PNLMO/NLMO/MO这些不管用哪种类型轨道，实空间函数（如总密度）的计算结果都是完全一致的，所以实空间函数的计算只要用MO就行了（.fch/.molden/.wfx/.wfn作为输入文件），基于其它的轨道来计算由于结果和基于MO的完全一致，所以没必要用它们。 而之所以我提到不要用NBO程序给出的轨道来计算（NLMO除外），一方面是因为这些NBO plot文件里往往没有正确记录或者根本没记录轨道占据数；另一方面是由于以某些类型轨道为基时（如NBO），整体的密度矩阵不是对角矩阵，此时以Multiwfn目前基于MO计算密度的方式来计算NBO表达的密度结果会和实际有一定偏差。 |
9楼2013-06-16 18:57:29
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
★ ★ ★
gkf高: 金币+3 2013-06-17 11:17:18
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朵）
klaus198710楼2013-06-16 19:01:50
50