请问cp2k-2026.2怎样实现使用OT时轨道能量和系数等波函数信息写入MOLDEN，并获得PDOS

reid · 发表于 Post on 3 day ago

本帖最后由 reid 于 2026-8-15 16:28 编辑

如题，尝试了多种办法，始终不能实现这一目的，请指教，谢谢。

8月15日：结合楼下几位热心朋友的提示，已解决：不要用cp2k产生的[Pseudo]，删掉，手写[Nval]字段，就可以了。

liuxp · 发表于 Post on 前天 07:42

https://manual.cp2k.org/trunk/me ... cular_orbitals.html 我看这个写了还没有试一试。

reid · 发表于 Post on 前天 15:29

liuxp 发表于 2026-8-15 07:42
https://manual.cp2k.org/trunk/methods/electronic_structure/molecular_orbitals.html 我看这个写了还没 ...

感谢回答。前天我尝试过了，获得的molden文件没法用Multiwfn获得PDOS。

UW_0728. · 发表于 Post on 前天 15:56

reid 发表于 2026-8-15 15:29
感谢回答。前天我尝试过了，获得的molden文件没法用Multiwfn获得PDOS。

具体是什么症状？

需要提示一点，2026.2版CP2K的molden文件与文章《详谈使用CP2K产生给Multiwfn用的molden格式的波函数文件》（http://bbs.keinsci.com/thread-31827-1-1.html）里面推荐修改的molden文件的组成不一样：文章里记录有效核电荷数用的是[Nval]，但CP2K自己的实现采用了跟ORCA一样的[Pseudo]。当前版本（2026.7.15）Multiwfn源代码里仍然只允许读明确被识别为ORCA输出的molden文件里的[Pseudo]字段（也就是其他程序产生的molden文件如果带有[Pseudo]是无法被Multiwfn读取的）。因此，如果原封不动用当前版本（2026.7.15）Multiwfn读取2026.2版CP2K生成的molden文件，仍然无法正确读取有效核电荷数，因此不能直接做波函数分析。请等下一个Multiwfn版本跟进，或者仍然继续手写[Nval]字段。

liuxp · 发表于 Post on 前天 16:05

UW_0728. 发表于 2026-8-15 15:56
具体是什么症状？

需要提示一点，2026.2版CP2K的molden文件与文章《详谈使用CP2K产生给Multiwfn用的mo ...

我上午测试的时候也是发现了这个问题，而且cube文件好像并没有输出这些

UW_0728. · 发表于 Post on 前天 16:14

liuxp 发表于 2026-8-15 16:05
我上午测试的时候也是发现了这个问题，而且cube文件好像并没有输出这些

cube的输出例程跟molden不一样（毕竟完全是两种不同的文件类型），但我印象中cube文件是正确记录晶胞和有效电荷信息的

reid · 发表于 Post on 前天 16:21

UW_0728. 发表于 2026-8-15 15:56
具体是什么症状？

需要提示一点，2026.2版CP2K的molden文件与文章《详谈使用CP2K产生给Multiwfn用的mo ...

确实是这样。一定要手写[Nval]字段

liuxp · 发表于 Post on 前天 16:39

UW_0728. 发表于 2026-8-15 16:14
cube的输出例程跟molden不一样（毕竟完全是两种不同的文件类型），但我印象中cube文件是正确记录晶胞和有 ...

也是2026.2办更新的吗？我没找到这个

UW_0728. · 发表于 Post on 前天 16:54

liuxp 发表于 2026-8-15 16:39
也是2026.2办更新的吗？我没找到这个

不是

Uus/pMeC6H4-/キ · 发表于 Post on yesterday 18:36

读取molden只认[Nval]不认[Pseudo]纯粹是后处理软件版本没跟上，不是CP2K的问题。[Pseudo]不仅是ORCA的既有惯例，而且能处理[Nval]没法恰当考虑的同一元素不同赝势的情形（在更新molden的读取代码的提议里早就说了这点），还可以把鬼原子用0来标记。

我印象中cube文件是正确记录晶胞和有效电荷信息

cube文件里有各个方向的格点数与格点间距，但是从http://bbs.keinsci.com/thread-58707-1-1.html来看这些信息不会自动转化成晶胞参数。往cube文件里写入有效核电荷的修改https://github.com/cp2k/cp2k/pull/4584在2026.1版和2026.2版里并没有完全覆盖所有种类的cube文件（但对http://bbs.keinsci.com/thread-46269-1-1.html所述用途已生效），剩余部分需要等https://github.com/cp2k/cp2k/pull/5718合并。

UW_0728. · 发表于 Post on yesterday 20:38

还有，你可以查看https://manual.cp2k.org/trunk/me ... _structure/dos.html，2026.2版本的CP2K支持生成可作图的PDOS曲线，可以试一下

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

[CP2K] 请问cp2k-2026.2怎样实现使用OT时轨道能量和系数等波函数信息写入MOLDEN，并获得PDOS