2026.1版CP2K已经可以在OT计算时往molden文件里输出轨道能量了，相应地对本文进行了修改

sobereva 发表于 2026-1-6 22:40
Multiwfn绘制DOS界面里定义片段时，把恰当的基函数加进去片段就完了

Multiwfn手册里给了球谐型基函数 ...

太感谢您了！

如果各位看到一个多小时前新鲜出炉的CP2K 2026.1的链接https://github.com/cp2k/cp2k/releases/tag/v2026.1，可能会在New Features里面看到一行

Molden output with cell/ECP/ghost/orbital(OT) information (#4364)

很可惜这恐怕是个假消息，开发者直接复制了一个特性请求https://github.com/cp2k/cp2k/issues/4353的标题上来，但是没有留意到相应的代码提交https://github.com/cp2k/cp2k/pull/4364并未完全实现所有的功能。这里面只有OT算法的占据轨道能量能按社长此文第4节所说的输出到molden里，而晶胞、ECP等都还没自动输出，仍然需要按此文所说的另外补充。

编辑：已经修改了

Better availability of occupied orbital eigenvalues from OT (#4364)

sobereva 发表于 2025-4-4 16:39
按下文基于CP2K的molden文件绘制DOS，其中说了怎么定义片段
使用Multiwfn绘制态密度(DOS)图考察电子结构 ...

你好sob老师，我用CP2K计算了有机分子吸附在金基底表面（金原子数目远多于分子原子数），在PDOS中把金基底对应的s,p,d轨道部分投影出来，比如我想具体看dz2轨道的成分投影，请教您该怎么做呢？

lycheeho 发表于 2025-4-3 15:30
你好Sob老师，我想在CP2K计算出的PDOS中提取比如Au元素的s轨道和dz2轨道成分，请问该怎么做呢？还有想请教 ...

按下文基于CP2K的molden文件绘制DOS，其中说了怎么定义片段
使用Multiwfn绘制态密度(DOS)图考察电子结构
http://sobereva.com/482（http://bbs.keinsci.com/thread-13242-1-1.html）

如果你搞不清楚基函数和原子轨道的对应关系，参考下文
利用布居分析判断基函数与原子轨道的对应关系
http://sobereva.com/418（http://bbs.keinsci.com/thread-10094-1-1.html）

计算级别可以

你好Sob老师，我想在CP2K计算出的PDOS中提取比如Au元素的s轨道和dz2轨道成分，请问该怎么做呢？还有想请教一下，我想结合PDOS分析有机分子吸附在金表面后分子轨道的排布，用PBE+D3(BJ)结合DZVP-MOLOPT-SR-GTH基组在精度上是否够用？

sobereva 发表于 2024-7-10 15:05
没完整输入文件没法说

老师，输入文件是文件夹中bulk_TDDFT.inp, 问题是无法生成.molden文件。另一输入文件BaBS4_TDDFT.inp在cp2k中计算时提示其中的Eu元素所用的TZVP-GTH-q17不存在于GTH_Basis_Sets中(Eu来自于BaB2S4.cif结构中的Ba做了一个Eu原子替换)，请问应该换成哪个基组能对其进行TDDFT计算呢？
链接：https://pan.baidu.com/s/1UnhXfj5tRwmfMzxsNBTmlQ?pwd=fvtn
提取码：fvtn

Hlin 发表于 2024-7-9 21:04
老师们，请问为什么用中的示例文件产生cp2k输入文件bulk_TDDFT.inp后，修改了其中的&DFT后的字段，但cp2k运 ...

没完整输入文件没法说

老师们，请问为什么用中的示例文件产生cp2k输入文件bulk_TDDFT.inp后，修改了其中的&DFT后的字段，但cp2k运行完后没有得到molden文件呢？

sobereva 发表于 2024-2-1 02:55
帖子里说了

因此没法得到PDOS，计算空轨道的PDOS需要空轨道的波函数

好的，明白了。谢谢卢老师。

ball2006 发表于 2024-1-31 11:46
卢老师，我用了Multiwfn 3.8Dev (2024-1-31)版本是可以导入之前给您上传的molden文件(空轨道我增加了500 ...

帖子里说了

但注意此时用OT产生的molden文件里仍然只有占据轨道的波函数

因此没法得到PDOS，计算空轨道的PDOS需要空轨道的波函数

sobereva 发表于 2024-1-31 07:53
问题在于你的输出文件里求解空轨道时没完全收敛，导致输出了额外信息，妨碍了Multiwfn读取数据：
  Lowe ...

卢老师，我用了Multiwfn 3.8Dev (2024-1-31)版本是可以导入之前给您上传的molden文件(空轨道我增加了500个)，但在画DOS图设置片段输出PDOS，无法输出PDOS的曲线。报错如下：

1. Fragment is saved
   
2. 
   
3. 
   
4. Fragment    1, number of basis functions:  2592
   
5. Fragment    2, number of basis functions:   360
   
6. Fragment    3, number of basis functions:   360
   
7. Fragment    4, number of basis functions:   360
   
8. Fragment    5, number of basis functions:   144
   
9. Fragment    6, number of basis functions:   864
   
10. Fragment    7, has not been defined
    
11. Fragment    8, has not been defined
    
12. Fragment    9, has not been defined
    
13. Fragment   10, has not been defined

复制代码

1. Calculating alpha orbital composition, please wait...
   
2.     1249 orbitals will be calculated
   
3. Progress: [##################################################]   100.0 %     \
   
4. Calculating beta orbital composition, please wait...
   
5.     1250 orbitals will be calculated
   
6. Progress: [##################################################]   100.0 %     /
   
7. Calculation of orbital composition took up wall clock time         0 s
   
8. Note: There are orbitals with zero energy, they are automatically ignored
   
9. 
   
10. Center of TDOS:    3.039651 eV
    
11. Center of PDOS  1:         NaN eV
    
12. Center of PDOS  2:         NaN eV
    
13. Center of PDOS  3:         NaN eV
    
14. Center of PDOS  4:         NaN eV
    
15. Center of PDOS  5:         NaN eV
    
16. Center of PDOS  6:         NaN eV
    
17. 
    
18. Note: The vertical dash line corresponds to HOMO level at   5.80484 eV

复制代码

sobereva 发表于 2024-1-31 07:53
问题在于你的输出文件里求解空轨道时没完全收敛，导致输出了额外信息，妨碍了Multiwfn读取数据：
  Lowe ...

明白了，谢谢卢老师。