计算化学公社

标题: 使用cp2k计算bader电荷，得到的ACF.dat中Si的电荷全为0 [打印本页]

作者
Author: 。。。 时间: 2024-11-12 14:20
标题: 使用cp2k计算bader电荷，得到的ACF.dat中Si的电荷全为0
本帖最后由。。。于 2024-11-12 14:21 编辑

各位老师好，我想证明加入了Fe会对MFI_Pt6中Pt上的电荷密度产生影响，所以用CP2K计算了bader电荷，但是我使用的是结构优化后的结构做单点计算，但是得到的ACF.dat中Si的CHARGE、MIN DIST和ATOMIC VOL均为0，并且ACF.dat中的原子的XYZ坐标和我提供的结构的坐标完全对不上，想请教各位老师是哪里的问题导致的这样的结果。以及我想证明Fe使得Pt富电子，可以看.out中Mulliken Population Analysis部分的电荷吗

作者
Author: sobereva 时间: 2024-11-13 04:47
AIM电荷是很差的选择
有好得多的原子电荷可用，仔细看
使用Multiwfn对周期性体系计算Hirshfeld(-I)、CM5和MBIS原子电荷
http://sobereva.com/712（http://bbs.keinsci.com/thread-46269-1-1.html）
一篇深入浅出、完整全面介绍原子电荷的综述文章已发表！
http://sobereva.com/714（http://bbs.keinsci.com/thread-46067-1-1.html）

作者
Author: 卡开发发 时间: 2024-11-13 08:22
（Bader-）AIM电荷确实不是很好的选择，或许也可以试试Chargemol的DDEC6电荷。对静电势和偶极矩复现的角度，CM5确实可能会更好。

若就你这里Bader没办法正确划分的原因讨论，是应当使用全电子计算得到的电子密度来划分原子盆然后用赝势得到的电子密度进行积分来规避数值问题，对Henckelman组的bader程序你可以使用-ref参数来指定专门用于划分原子盆的参考数据，你可以把全电子得到的密度网格文件名放在-ref参数后，这类似于VASP中由AECCAR0和2所合成的CHGCAR_sum。

作者
Author: funok 时间: 2025-4-4 10:37
坐标对不上应该是Bohr和Angstrom之间的转换，ACF.dat中的XYZ坐标乘上0.5291772083，就可以得到你xyz文件中的坐标。

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