sobereva 发表于 2023-9-12 19:45
先产生molden文件
详谈使用CP2K产生给Multiwfn用的molden格式的波函数文件
http://sobereva.com/651（h ...

好的，谢谢sob老师

liheng 发表于 2023-9-12 17:32
好的sob老师，我马上做cp2k计算，但是这个布居分析有没有详细点介绍或者推荐啊？

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

然后给Multiwfn，主功能7的子功能5里面选1，得到各个基函数的Mulliken布居数，按下文了解基函数与原子轨道的对应关系便知道原子轨道布居数
利用布居分析判断基函数与原子轨道的对应关系
http://sobereva.com/418（http://bbs.keinsci.com/thread-10094-1-1.html）


另外，北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里有完整的分析例子

sobereva 发表于 2023-9-12 16:31
这种图太抽象，只能根据原子轨道分布的大致范围、密度差图中的电子密度增/减的区域来猜，远没有布居分析 ...

好的sob老师，我马上做cp2k计算，但是这个布居分析有没有详细点介绍或者推荐啊？

liheng 发表于 2023-9-12 03:05
sob老师，我用的VASP的CHGCAR来计算的差分电荷密度。如果要用cp2k也行，只是您说的做Mulliken布居分析我 ...

这种图太抽象，只能根据原子轨道分布的大致范围、密度差图中的电子密度增/减的区域来猜，远没有布居分析有说服力

说清楚算的什么体系、用的什么程序
一般来说，这种问题可以对电子转移之前和之后的状态计算原子轨道的布居数，彼此求差，就定量知道各个原子轨道的电子密度分布变化的贡献了
如果用的VASP，具体实现我不清楚。若用CP2K，可以用基于其molden文件用Multiwfn做Mulliken布居分析给出这样的信息。