| 查看: 2518 | 回复: 10 | |||
[求助]
Multiwfn读取NBO的.40文件是对应MO吗？
|
|
如题，manual里说NBO的.40文件是NLMO，我听说这种轨道的性质和MO是相同的，请问这种说法是否正确？ 另外，我试着做了.40的电子密度格点文件，又对同一结构(相同坐标)再次用高斯计算能量，得到MO轨道的fch文件，同样做电子密度格点。两个格点文件的质量选择相同，二者相减会得到差值，也就是说.40和MO电子密度格点不一样。这是正常的吗？ |
回复此楼» 猜你喜欢
武汉纺织大学全国重点实验室杨应奎教授团队拟招17-19位硕博生！
已经有0人回复
-
材料与化工专硕接受跨专业调剂（重庆市属高校）
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有219人回复
-
武汉纺织大学全国重点实验室杨应奎教授团队拟招17-19位硕博生！
已经有0人回复
-
上海工程技术大学生态与功能纺织品开发研究团队招收学硕专硕研究生
已经有0人回复
-
青岛大学化学化工学院2026年招收调剂硕士研究生
已经有0人回复
-
新能源材料方向考研调剂招生（长江师范学院）
已经有0人回复
-
0854电子信息272求调剂！
已经有0人回复
-
化学-化工-材料-食品-生物-环境2026年研究生调剂（可以不考数学）
已经有0人回复
-
化学-化工-材料-食品-生物-环境-农业 2026年研究生调剂（可以不考数学）
已经有0人回复
-
070300一志愿211，312分求调剂院校
已经有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
