计算化学公社

标题: Multiwfn结合CP2K进行NAdO分析时出错 [打印本页]

作者
Author: 2429783189 时间: 2025-12-24 17:03
标题: Multiwfn结合CP2K进行NAdO分析时出错
本帖最后由 2429783189 于 2025-12-25 10:43 编辑

(, 下载次数 Times of downloads: 0) 各位老师，我在服务器上使用Multiwfn_noGUI版本结合CP2K计算产生的molden文件，想参考http://bbs.keinsci.com/thread-47176-1-1.html此帖做片段间的NAdO分析，但是进行到上图输入20做NAdO分析时，发生了如上图及下文所示的提示：
Warning: Energies of all orbitals are zero, so you need to manually input HOMO index
Input index of HOMO of alpha spin, e.g. 300
Input index of HOMO of beta spin, e.g. 299然后就因为进行不下去而报错了
我产生molden文件所用的inp文件已上传为附件“Fe.inp”，已正确在molden文件开头插入定义晶胞信息和当前赝势基组描述的各元素的价电子数信息，做IGMH分析时没有问题。在inp文件中要求将Kohn-Sham矩阵导出到.csr文件，但按照博文，分析片段间NAdO时没用到这个文件。我丢给Multiwfn_noGUI执行指令的txt文件已上传为附件“NADO-FE.txt"，想请教一下这种情况应该如何解决。

作者
Author: sobereva 时间: 2025-12-26 04:23
按照屏幕上的提示输入。你的截图里没体现你输入了什么，没法回答

作者
Author: 2429783189 时间: 2025-12-26 10:53
本帖最后由 2429783189 于 2025-12-26 17:39 编辑

sobereva 发表于 2025-12-26 04:23
按照屏幕上的提示输入。你的截图里没体现你输入了什么，没法回答

好的老师，我载入molden文件后输入0查看轨道，输出如下
(, 下载次数 Times of downloads: 0)
一开始我按照屏幕上的提示输入：1242和5835，结果算出来的NAdO轨道本征值都是0，
(, 下载次数 Times of downloads: 0)
是否是molden文件没有记录轨道信息的原因呢，我计算使用的是OT，是否需要按照http://bbs.keinsci.com/thread-31827-1-1.html这帖的方法从out文件里读轨道，然后再按照博文教程做NAdO分析呢？

作者
Author: 2429783189 时间: 2025-12-26 17:45
本帖最后由 2429783189 于 2025-12-26 18:19 编辑

sobereva 发表于 2025-12-26 04:23
按照屏幕上的提示输入。你的截图里没体现你输入了什么，没法回答

(, 下载次数 Times of downloads: 0) 额外补充一个信息，当我尝试按照上图对两个原子之间的相互作用进行NAdO分析时，导入生成的csr文件，倒是最终可以成功，不过只有alpha NAdOs不为0，beta NAdOs为0，如下图，不知道这样是否正常。随后，当输入y让Multiwfn载入当前目录下的NAdOs.mwfn时，发生闪退。我这是个自旋多重度为1的反铁磁耦合体系。
(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)

作者
Author: sobereva 时间: 2025-12-27 02:14

2429783189 发表于 2025-12-26 10:53
好的老师，我载入molden文件后输入0查看轨道，输出如下

一开始我按照屏幕上的提示输入：1242和5835， ...

用对角化产生的molden文件，如果还有问题，把molden文件和csr文件压缩后传网盘我看看

作者
Author: 2429783189 时间: 2025-12-27 14:15
本帖最后由 2429783189 于 2025-12-27 16:41 编辑

sobereva 发表于 2025-12-27 02:14
用对角化产生的molden文件，如果还有问题，把molden文件和csr文件压缩后传网盘我看看

老师，这是个反铁磁性的d族金属表面slab模型，不用OT用对角化很难收敛，设置正确的电子组态初猜，读取OT的波函数同样不好收敛，http://sobereva.com/665此帖的方法已经全部尝试，均无法解决。根据先前的尝试，OT算出molden导入csr文件时可以成功输出NAdO的本征值，但轨道的能量都是0，生成的NAdOs.mwfn导入到Multiwfn会发生闪退。请问老师，是否有方法能用OT算的molden做NAdO，还是只能硬着头皮解决对角化的收敛呢？

作者
Author: sobereva 时间: 2025-12-28 14:32

2429783189 发表于 2025-12-27 14:15
老师，这是个反铁磁性的d族金属表面slab模型，不用OT用对角化很难收敛，设置正确的电子组态初猜，读取OT ...

本来涉及金属表面的体系就应该对角化+smearing。OT就算收敛了，都有可能收敛到不稳定波函数

作者
Author: 2429783189 时间: 2026-1-5 17:44

sobereva 发表于 2025-12-28 14:32
本来涉及金属表面的体系就应该对角化+smearing。OT就算收敛了，都有可能收敛到不稳定波函数

谢谢老师，换成对角化算出来的molden之后没有任何问题了，但是仍然有一个疑惑，我这个体系是正磁矩和负磁矩交替排列的，但总的自旋多重度为0，NAdO分析结果中有Eigenvalues of alpha NAdOs和Eigenvalues of beta NAdOs，这两个应该讨论哪个呢？

作者
Author: sobereva 时间: 2026-1-6 00:30

2429783189 发表于 2026-1-5 17:44
谢谢老师，换成对角化算出来的molden之后没有任何问题了，但是仍然有一个疑惑，我这个体系是正磁矩和负磁 ...

同时讨论

作者
Author: 2429783189 时间: 2026-1-8 01:12

sobereva 发表于 2026-1-6 00:30
同时讨论

好的，感谢老师的回复

欢迎光临计算化学公社 (http://bbs.keinsci.com/)