简单Multiwfn使用问题和波函数分析问题答疑专贴

Zhengmaxmell · 发表于 Post on 2019-3-16 08:17:15

sobereva 发表于 2019-3-16 07:26
gview打开输出文件，保存成新的输入文件，0 1改成1 2就完了

好的，非常感谢您的帮助！

Zhengmaxmell · 发表于 Post on 2019-3-16 10:27:21

sobereva 发表于 2019-3-16 07:26
gview打开输出文件，保存成新的输入文件，0 1改成1 2就完了

sob老师，您好。我按照您的指点将natural state的输出文件另存为+1state的输入文件并更改电荷与自旋多重度，然后去掉opt关键词来获得对应wfn文件，但是gaussian运行结束后，没有在预设的文件夹中获得对应的wfn文件，反而在scratch文件中获得了1 2 1.5 6 1.5 9 1.0.的文件，请问这是什么情况？

usedgo · 发表于 Post on 2019-3-16 10:39:59

sobereva 发表于 2019-3-16 07:30
得不到是操作问题。你先把图中2c、3c、4c都找到了，然后直接让Multiwfn搜索36中心的，必定直接就找到了那 ...

感谢sob老师的回复，另外，在搜索36c的时候，搜索了好长时间都没反应。是因为36中心的计算量太大吗？

sobereva · 发表于 Post on 2019-3-16 11:08:17

Zhengmaxmell 发表于 2019-3-16 10:27
sob老师，您好。我按照您的指点将natural state的输出文件另存为+1state的输入文件并更改电荷与自旋多重 ...

把输入文件里连接关系部分去掉

sobereva · 发表于 Post on 2019-3-16 11:09:50

usedgo 发表于 2019-3-16 10:39
感谢sob老师的回复，另外，在搜索36c的时候，搜索了好长时间都没反应。是因为36中心的计算量太大吗？

如果你已经把1~4c的轨道都挑出来了，之后选择3 Set the number of centers in the next exhaustive search，输入36，然后选选项2搜索36c，应当很快就出来才对。
如果不是这样，把你选择搜索36c选项时的截图发出来，或者把Multiwfn中从一开始到最后的所有输出拷到一个文本文件里传上来

Zhengmaxmell · 发表于 Post on 2019-3-16 13:01:23

sobereva 发表于 2019-3-16 11:08
把输入文件里连接关系部分去掉

谢谢sob老师，问题已经解决了！

moritaichi · 发表于 Post on 2019-3-16 23:43:24

老师您好，我想询问一下我在使用Multiwfn的一个小小的疑惑。
读取例子里自带的水的fch文件，使用SCPA方法计算得到了水的TDOS图。
提取数据后用origin画图，负无穷到费米能级的积分只有0.3，按理来积TDOS应该是总电子数10才对呀。
不知道我的理解哪里出了差错，有时间的话还请指点一下，不胜感激。

sobereva · 发表于 Post on 2019-3-17 06:20:03

moritaichi 发表于 2019-3-16 23:43
老师您好，我想询问一下我在使用Multiwfn的一个小小的疑惑。
读取例子里自带的水的fch文件，使用SCPA方法 ...

有诸多因素：
1 你的横坐标范围是默认值，而不是所有占据轨道所处的能级范围
2 单位问题
3 你对费米能级和TDOS图理解有问题。Multiwfn并不给出费米能级，孤立体系也没有确切的费米能级这个概念。如果你只是从低到高一直积分到HOMO轨道能量，那么由于HOMO的竖线展宽后在其右侧也有分布，将之忽略掉也是不对的

moritaichi · 发表于 Post on 2019-3-17 11:47:22

sobereva 发表于 2019-3-17 06:20
有诸多因素：
1 你的横坐标范围是默认值，而不是所有占据轨道所处的能级范围
2 单位问题

老师您好，是这样的，我的横坐标是参考了我用VASP做出的冰晶体态密度范围的，单位也在multiwfn里换成了eV，是-25eV到15eV，按理说不会差太多；
单位统统换成了eV，而且我尝试过，无论单位是用eV还是a.u.，由于是态密度，横纵坐标会一起变化，态密度的积分两个单位都是0.3.
费米能级是我按照multiwfn里的显示费米能级所得。
我也尝试过从负无穷积分到正无穷，这样包括了右侧分布，但积分值仍然只有0.5。

sobereva · 发表于 Post on 2019-3-17 11:59:37

moritaichi 发表于 2019-3-17 11:47
老师您好，是这样的，我的横坐标是参考了我用VASP做出的冰晶体态密度范围的，单位也在multiwfn里换成了eV ...

举个例子，H2分子，这是TDOS图

把导出的DOS_curve.txt

DOS_curve.txt (70.31 KB, 下载次数 Times of downloads: 0) 弄到origin里积分

积分结果为0.2。你在屏幕上可以看到“4 Set scale ratio for DOS curve, current: 0.10000”，即是说默认是对DOS乘了0.1，因此将积分值除以0.1，就是2，正好对应图中的两个能级

moritaichi · 发表于 Post on 2019-3-17 12:27:14

sobereva 发表于 2019-3-17 11:59
举个例子，H2分子，这是TDOS图

把导出的DOS_curve.txt弄到origin里积分

嗯嗯，就是这个4ratio，感谢老师的耐心讲解，真的很感谢。

498746012 · 发表于 Post on 2019-3-18 09:28:36

sob老师您好！multiwfn的pipek定域化功能可以看出某个键是共价键还是离子键吗？判断共价键和离子键的最好方法是什么？是NRT中的bond order吗里边有离子键的成分和共价键的成分

sobereva · 发表于 Post on 2019-3-18 09:36:42

498746012 发表于 2019-3-18 09:28
sob老师您好！multiwfn的pipek定域化功能可以看出某个键是共价键还是离子键吗？判断共价键和离子键的最好方 ...

这里写得非常详细
Multiwfn支持的分析化学键的方法一览
http://sobereva.com/471

Zhengmaxmell · 发表于 Post on 2019-3-18 10:24:53

sob老师，您好。之前向您请教过如何利用Multiwfn计算福井函数。现在我能够将graph导出，并且能够分辨f-的分布，我想请教如何能够导出或显示各原子f-的具体数值？

498746012 · 发表于 Post on 2019-3-18 20:08:14

sobereva 发表于 2019-3-18 09:36
这里写得非常详细
Multiwfn支持的分析化学键的方法一览
http://sobereva.com/471

谢谢sob老师

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Multiwfn使用咨询] 简单Multiwfn使用问题和波函数分析问题答疑专贴