计算化学公社

标题: vasp怎样计算 Mulliken 电荷分布？ [打印本页]

作者
Author: wmg166 时间: 2021-11-24 15:54
标题: vasp怎样计算 Mulliken 电荷分布？
怎样计算布局分析(population analysis)？用 vaspkit怎样生成计算Mulliken 电荷的INCAR？ vaspkit 里 “ BD) Bader Charge Analysis ” 就是 Mulliken 电荷分布？吗

作者
Author: sobereva 时间: 2021-11-24 16:31
Bader电荷正经名字叫AIM电荷，不同原子电荷的差异在《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）里说得很清楚，仔细看此文搞清楚关系

VASP不清楚。CP2K可以直接算或结合Multiwfn算Mulliken电荷

作者
Author: wmg166 时间: 2021-11-24 17:15

sobereva 发表于 2021-11-24 16:31
Bader电荷正经名字叫AIM电荷，不同原子电荷的差异在《原子电荷计算方法的对比》（http://www.whxb.pku.edu. ...

3Q，万分感谢

作者
Author: granvia 时间: 2021-11-24 20:23
Mulliken电荷是使用原子轨道基组时采用的布居分析方法。VASP是用平面波基组来展开波函数的，并没有以原子为中心的原子轨道，所以不可能存在Mulliken电荷

作者
Author: 卡开发发 时间: 2021-11-24 20:52

granvia 发表于 2021-11-24 20:23
Mulliken电荷是使用原子轨道基组时采用的布居分析方法。VASP是用平面波基组来展开波函数的，并没有以原子为 ...

理论上可以将平面波投影到局域轨道（例如原子轨道）上然后再进行Mulliken电荷分析，比如CASTEP就可以做。不过没人去给vasp写Mulliken罢了，实际上之前PeriodicNBO也有上述过程，理论上魔改下也能做。

作者
Author: granvia 时间: 2021-11-24 22:42

卡开发发发表于 2021-11-24 20:52
理论上可以将平面波投影到局域轨道（例如原子轨道）上然后再进行Mulliken电荷分析，比如CASTEP就可以做。 ...

VASP自带的算法是投影到以各原子为中心的球谐函数上，还考虑Wigner–Seitz半径吧，问题是这种投影就不是严格的波函数展开了，因为这里的原子轨道不构成完备集了吧？

作者
Author: 卡开发发 时间: 2021-11-25 00:35

granvia 发表于 2021-11-24 22:42
VASP自带的算法是投影到以各原子为中心的球谐函数上，还考虑Wigner–Seitz半径吧，问题是这种投影就不是 ...

真要做确实不能用那个自带的。有专门的算法用spillage来衡量这种投影的质量，比如Plato这些程序（10.1103/PhysRevB.62.4899）就是采用这种方式来构造数值基组，肯定达不到绝对完备。

作者
Author: wmg166 时间: 2021-11-26 11:38

卡开发发发表于 2021-11-25 00:35
真要做确实不能用那个自带的。有专门的算法用spillage来衡量这种投影的质量，比如Plato这些程序（10.1103 ...

大神一席话，在下醍醐灌顶，恍然大悟

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