计算化学公社

 找回密码 Forget password
 注册 Register
Views: 166|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[Multiwfn使用咨询] Multiwfn结合CP2K进行NAdO分析时出错

[复制链接 Copy URL]

88

帖子

0

威望

416

eV
积分
504

Level 4 (黑子)

本帖最后由 2429783189 于 2025-12-25 10:43 编辑

各位老师,我在服务器上使用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",想请教一下这种情况应该如何解决。


Fe.inp

7.86 KB, 下载次数 Times of downloads: 2

NADO-FE.txt

91 Bytes, 下载次数 Times of downloads: 2

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

2#
发表于 Post on 2025-12-26 04:23:47 | 只看该作者 Only view this author
按照屏幕上的提示输入。你的截图里没体现你输入了什么,没法回答
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

88

帖子

0

威望

416

eV
积分
504

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-12-26 10:53:42 | 只看该作者 Only view this author
本帖最后由 2429783189 于 2025-12-26 17:39 编辑
sobereva 发表于 2025-12-26 04:23
按照屏幕上的提示输入。你的截图里没体现你输入了什么,没法回答

好的老师,我载入molden文件后输入0查看轨道,输出如下

一开始我按照屏幕上的提示输入:1242和5835,结果算出来的NAdO轨道本征值都是0,

是否是molden文件没有记录轨道信息的原因呢,我计算使用的是OT,是否需要按照http://bbs.keinsci.com/thread-31827-1-1.html这帖的方法从out文件里读轨道,然后再按照博文教程做NAdO分析呢?

88

帖子

0

威望

416

eV
积分
504

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2025-12-26 17:45:05 | 只看该作者 Only view this author
本帖最后由 2429783189 于 2025-12-26 18:19 编辑
sobereva 发表于 2025-12-26 04:23
按照屏幕上的提示输入。你的截图里没体现你输入了什么,没法回答

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




6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

5#
发表于 Post on 2025-12-27 02:14:58 | 只看该作者 Only view this author
2429783189 发表于 2025-12-26 10:53
好的老师,我载入molden文件后输入0查看轨道,输出如下

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

用对角化产生的molden文件,如果还有问题,把molden文件和csr文件压缩后传网盘我看看
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

88

帖子

0

威望

416

eV
积分
504

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2025-12-27 14:15:28 | 只看该作者 Only view this author
本帖最后由 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,还是只能硬着头皮解决对角化的收敛呢?

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

7#
发表于 Post on 2025-12-28 14:32:36 | 只看该作者 Only view this author
2429783189 发表于 2025-12-27 14:15
老师,这是个反铁磁性的d族金属表面slab模型,不用OT用对角化很难收敛,设置正确的电子组态初猜,读取OT ...

本来涉及金属表面的体系就应该对角化+smearing。OT就算收敛了,都有可能收敛到不稳定波函数
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

88

帖子

0

威望

416

eV
积分
504

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2026-1-5 17:44:44 | 只看该作者 Only view this author
sobereva 发表于 2025-12-28 14:32
本来涉及金属表面的体系就应该对角化+smearing。OT就算收敛了,都有可能收敛到不稳定波函数

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

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

9#
发表于 Post on 2026-1-6 00:30:48 | 只看该作者 Only view this author
2429783189 发表于 2026-1-5 17:44
谢谢老师,换成对角化算出来的molden之后没有任何问题了,但是仍然有一个疑惑,我这个体系是正磁矩和负磁 ...

同时讨论
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

88

帖子

0

威望

416

eV
积分
504

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2026-1-8 01:12:36 | 只看该作者 Only view this author

好的,感谢老师的回复

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-23 23:52 , Processed in 0.282552 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list