Multiwfn已经可以载入CP2K直接产生的含有[Pseudo]和[Cell] AU的molden文件

sobereva

CP2K 2026.2在产生molden文件方面做了改进，&MO_MOLDEN里加上WRITE_CELL T可以在molden文件里直接写入[Cell]字段记录记录晶胞信息（CP2K开发者借鉴了我在http://sobereva.com/651里对molden格式的扩展。CP2K默认以原子单位记录晶胞信息，因此字段是[Cell] AU），加上WRITE_PSEUDO T可以在molden文件里直接写入[Pseudo]字段记录各个原子的价电子数，这点和ORCA一样。

相应地，我在今日更新的Multiwfn 2026.8.19里做了改进：

• 在《使用Multiwfn非常便利地创建CP2K程序的输入文件》（http://sobereva.com/587）介绍的Multiwfn产生CP2K的输入文件的界面里，要求产生molden文件后，产生的输入文件里就会包含以上两个设置（非周期性的情况不含WRITE_CELL T）
• 可以正常载入CP2K产生的含有[Pseudo]和[Cell] AU的molden文件
  

因此，Multiwfn用户对CP2K产生的周期性的、赝势下计算的波函数文件进行分析，就不需要每次手动添加[Cell]和[Nval]信息了，方便多了。

PoorChaos

太好了

Uus/pMeC6H4-/キ

CP2K产生的molden文件的晶格尺寸和原子坐标不一定使用原子单位的[Cell] AU和[Atoms] AU，如果&MO_MOLDEN/UNIT关键词设置成非默认的ANGSTROM，就会输出标记为[Cell] Angs和[Atoms] Angs的以埃为单位的数据。（原始的molden文件规范已经有了[Atoms] (Angs|AU)的标记来支持切换长度单位。）最合适的逻辑不是直接假设用了什么单位，而是按[Cell]或[Atoms]后面的标签切换缩放因子。

sobereva

Uus/pMeC6H4-/キ 发表于 2026-8-19 14:31
CP2K产生的molden文件的晶格尺寸和原子坐标不一定使用原子单位的[Cell] AU和[Atoms] AU，如果&MO_MOLDEN/UN ...

Multiwfn发现[Cell]那行带AU或Au才会视为Bohr，否则默认当成埃，因此两种单位都能兼容，考虑得非常周全。不需要特意考虑是否[Cell]后面有Angs