计算化学公社

标题: 求助：差分电荷密度能否识别转移/极化的电荷属于哪个轨道（s, p, or d) [打印本页]

作者
Author: liheng 时间: 2023-9-10 22:21
标题: 求助：差分电荷密度能否识别转移/极化的电荷属于哪个轨道（s, p, or d)
本帖最后由 liheng 于 2023-9-10 22:23 编辑

最近投稿的文章，审稿人其中一个问题是：如何区分不同轨道电子对总电荷密度变化的贡献？但是我不知道怎么回。因为做差分电荷密度图时用的CHGCAR没有相关信息啊。是否有其他办法可以判断不同轨道电子对电荷转移的贡献呢

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

作者
Author: liheng 时间: 2023-9-12 03:05

sobereva 发表于 2023-9-10 22:42
说清楚算的什么体系、用的什么程序
一般来说，这种问题可以对电子转移之前和之后的状态计算原子轨道的布居 ...

sob老师，我用的VASP的CHGCAR来计算的差分电荷密度。如果要用cp2k也行，只是您说的做Mulliken布居分析我还不太了解，不过问题应该不大。再简单描述一下我的问题吧。我计算的是铯原子在钼金属上的吸附体系，其中要分析功函数的变化。我在一篇老的PRL（DOI:https://doi.org/10.1103/PhysRevLett.48.1128）里面看到作者做了差分电荷密度图（见下），就牵扯到电荷的转移/极化什么的。然后他说铯在表面吸附时价电子（6s）会向表面极化，而5p电子会向真空极化（也就是和6s的极化相反），然后我在文中就引用了这个结论。再然后审稿人就问通过差分电荷密度图是如何判断5p和6s的极化方向的。我但是我在那篇PRL里面也没有看出作者是如何判断的。

作者
Author: sobereva 时间: 2023-9-12 16:31

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

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

作者
Author: liheng 时间: 2023-9-12 17:32

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

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

作者
Author: sobereva 时间: 2023-9-12 19:45

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）里有完整的分析例子
(, 下载次数 Times of downloads: 3)

作者
Author: liheng 时间: 2023-9-12 19:49

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

好的，谢谢sob老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)