从CP2K中的的WANNIER_CUBES输出的是什么？

赵小壮

我在CP2K中想要输出一个电荷密度网格，要求这个电荷网格基于的基的localize方法与 计算与wannier function center的相同。我在手册中找到了这个东西

https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/XAS/LOCALIZE/PRINT/WANNIER_CUBES.html

然而这个只说了一句介绍：

Controls the printing of the wannier functions

我不知道这个输出的具体是什么东西，然后我自己测试了一下一个纯水体系输出的wannier_cube文件，我发现格点上的值取负，再乘上格子体积dx*dy*dz后，
加和为4.1164810168439532 |e|  非零因此应该不是全电子+核的电荷密度网格，而电荷加和分配到每个水分子上又是0.几，说明也不是价电子的电荷密度网格，不然应该是8。

所以请问这个到底输出的是什么东西啊？谢谢

sobereva

不知道你的实际目的是什么
Multiwfn基于CP2K的molden文件可以做轨道定域化得到最大定域化Wannier函数，可以给出质心位置以及cub文件

详谈使用CP2K产生给Multiwfn用的molden格式的波函数文件
http://sobereva.com/651（http://bbs.keinsci.com/thread-31827-1-1.html）

Multiwfn的轨道定域化功能的使用以及与NBO、AdNDP分析的对比
http://sobereva.com/380（http://bbs.keinsci.com/thread-6053-1-1.html）

赵小壮

sobereva 发表于 2023-11-17 03:31
不知道你的实际目的是什么
Multiwfn基于CP2K的molden文件可以做轨道定域化得到最大定域化Wannier函数，可 ...

谢谢sob老师，我的目的是对比价电子密度等值面在localization之前和之后的差异，

并要求这个localization是localization wanniner function的方法

我现在有CP2K局域化之前的cube文件，记录的是价电子密度，用的命令是FORCE_EVAL/DFT/PRINT/E_DENSITY_CUBE

我现在没有的是在localization后输出的cube文件，我只找到了WANNIER_CUBES这个命令，但输出的经检验不是价电子的电荷密度网格，难道是波函数网格？

最后我理解sob老师您的回复的意思是,Multiwfn是可以基于CP2K的轨道文件做相同的localization并给出的cub文件，并且这个cub文件是价电子密度网格的文件吗（在CP2K给出的轨道文件时用的时赝势的情况下）

sobereva

赵小壮 发表于 2023-11-20 10:32
谢谢sob老师，我的目的是对比价电子密度等值面在localization之前和之后的差异，

并要求这个localizat ...

我没用过那个，可能给你的是localization wanniner function的轨道波函数的格点数据，凭借图形判断

Multiwfn做完轨道定域化之后，内存里的轨道就变成了定域化轨道，你可以基于这些轨道计算和导出各种格点数据，比如轨道波函数值的cube、轨道概率密度的cube、特定一些轨道贡献的电子密度的cube等等，极其灵活。