| 查看: 2523 | 回复: 10 | |||
klaus1987兑换贵宾
|
[求助]
Multiwfn读取NBO的.40文件是对应MO吗？
|
|
如题，manual里说NBO的.40文件是NLMO，我听说这种轨道的性质和MO是相同的，请问这种说法是否正确？ 另外，我试着做了.40的电子密度格点文件，又对同一结构(相同坐标)再次用高斯计算能量，得到MO轨道的fch文件，同样做电子密度格点。两个格点文件的质量选择相同，二者相减会得到差值，也就是说.40和MO电子密度格点不一样。这是正常的吗？ |
回复此楼» 猜你喜欢
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有73人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
-
求ENDNOTE破解软件安装包
已经有4人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
-
26博士申请
已经有0人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 建议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
klaus1987
超级版主
- 应助: 14 (小学生)
- 金币: 196.7
- 红花: 3
- 帖子: 122
- 在线: 117.5小时
- 虫号: 1561644
- 注册: 2012-01-03
- 性别: GG
- 专业: 高分子材料结构与性能
3楼2013-06-14 09:59:40
klaus1987
兑换贵宾
- 应助: 14 (小学生)
- 金币: 196.7
- 红花: 3
- 帖子: 122
- 在线: 117.5小时
- 虫号: 1561644
- 注册: 2012-01-03
- 性别: GG
- 专业: 高分子材料结构与性能
4楼2013-06-14 10:49:50
klaus1987
主管区长
- 应助: 14 (小学生)
- 金币: 196.7
- 红花: 3
- 帖子: 122
- 在线: 117.5小时
- 虫号: 1561644
- 注册: 2012-01-03
- 性别: GG
- 专业: 高分子材料结构与性能
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
klaus1987
实习版主
- 应助: 14 (小学生)
- 金币: 196.7
- 红花: 3
- 帖子: 122
- 在线: 117.5小时
- 虫号: 1561644
- 注册: 2012-01-03
- 性别: GG
- 专业: 高分子材料结构与性能
7楼2013-06-16 15:29:56
klaus1987
超级版主
- 应助: 14 (小学生)
- 金币: 196.7
- 红花: 3
- 帖子: 122
- 在线: 117.5小时
- 虫号: 1561644
- 注册: 2012-01-03
- 性别: GG
- 专业: 高分子材料结构与性能
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
