详谈使用CP2K产生给Multiwfn用的molden格式的波函数文件

UW_0728.

2026.2版（预计六月底或七月初发布）CP2K将直接支持在生成的Molden下写入晶胞和有效核电荷数信息，分别为[Cell]和[Pseudo]板块（前者写入晶格矢量，后者与ORCA中[Pseudo]的格式一致：“[元素名]  [序号]  [有效核电荷数]”；考虑过[Nval]，但鉴于个别计算中可能对同种元素的不同原子使用不同有效核电荷数的赝势，使用[Pseudo]处理起来更方便，相关讨论可以请看https://github.com/cp2k/cp2k/issues/4353#issuecomment-4265218000）；相应关键词分别为“WRITE_CELL”和“WRITE_PSEUDO”，需手动开启。

此外[Cell]中晶胞参数的单位与[Atom]一致，由UNIT关键词统一处理，同时支持Angstrom（`[Cell] Angs`）和Bohr（`[Cell] AU`）。

输入示例如下：

1. &MO_MOLDEN
   
2.   NDIGITS 9
   
3.   UNIT ANGSTROM
   
4.   WRITE_CELL T
   
5.   WRITE_PSEUDO T
   
6. &END MO_MOLDEN

复制代码

Molden开头示例：

1. [Molden Format]
   
2. [Cell] Angs
   
3.      3.574965       0.000000       0.000000
   
4.      0.000000       3.574965       0.000000
   
5.      0.000000       0.000000       3.574965
   
6. [Atoms] Angs
   
7. C      1     6       0.000000      -0.000000       0.000000
   
8. C      2     6       2.681237       0.893728       2.681237
   
9. C      3     6      -0.000000       1.787483       1.787483
   
10. C      4     6       2.681237       2.681237       0.893728
    
11. C      5     6       1.787483       0.000000       1.787483
    
12. C      6     6       0.893728       0.893728       0.893728
    
13. C      7     6       1.787483       1.787483      -0.000000
    
14. C      8     6       0.893728       2.681237       2.681237
    
15. [Pseudo]
    
16. C      1     4
    
17. C      2     4
    
18. C      3     4
    
19. C      4     4
    
20. C      5     4
    
21. C      6     4
    
22. C      7     4
    
23. C      8     4

复制代码

我试着根据CP2K中molden文件输出的更新，做了一版新的fileIO.f90：

fileIO.f90 (371.56 KB, 下载次数 Times of downloads: 4)

2026-6-4 01:23 上传 Uploaded

点击下载
Click to download

改动就两处：一是Molden那块读取[Cell]的时候不再直接Bohr转Ang（但考虑到向后兼容性，当molden中[Cell]后没有跟单位说明时依然按Angstrom处理）；二是把[Pseudo]读取时强制要求必须是ORCA的molden文件的限制去掉了（同时把位置挪到了[Nval]读取的下面）。我用卢老师的博文《使用Multiwfn结合CP2K计算晶体中原子的氧化态》（http://bbs.keinsci.com/thread-45876-1-1.html）里的例子测试了一下，结果显示没什么问题（跟改动之前算得的结果一样）。如果卢老师有空的话可以烦请看一看。

sobereva

UW_0728. 发表于 2026-6-3 19:52
2026.2版（预计六月底或七月初发布）CP2K将直接支持在生成的Molden下写入晶胞和有效核电荷数信息，分别为[C ...

等CP2K新版本正式版出来后我会跟进