求助VASP计算变形电荷密度的方法和步骤

lora2021

请教大家：
       最近在比较两个化合物的电子结构差异，看到文献计算了两种物质的变形电荷密度，如图所示，我查了一下文献，有给出计算公式：
SCF收敛以后体系的电荷密度减去组成它的各个原子在自由状态下的球对称的电荷密度，也称为变形电荷密度
(Deformation charge density)
∆𝜌 = 𝜌(𝐴𝐵𝑠𝑒𝑙𝑓−𝑐𝑜𝑛𝑠𝑖𝑠𝑡𝑒𝑛𝑡) − 𝜌(𝐴𝐵𝑎𝑡𝑜𝑚𝑖𝑐)
这里面的第一项容易理解，就是化合物AB的电荷密度，后面的“各个原子在自由状态下的球对称的电荷密度”，这个要怎么计算？

看说明书是这样解释：新建立一个文件夹，做非自洽计算，输出原子电荷密度的叠加，有点晕
这个是要怎么做？

lora2021

也就是后面的非自洽计算，K点要怎么设置？scf计算是3 2 1，INCAR里面设置和band或都是计算一样吗？

lora2021

问题解决了，嘿嘿。

柠檬鬼lemon

lora2021 发表于 2023-10-23 19:17
问题解决了，嘿嘿。

你好，我现在也要计算变形电荷密度，请问你是怎么做的呀，在做非自洽计算是怎么改变INCAR的参数呀，将ICHARG =12吗，其他还有没有什么需要改的，那KPOINTS需要改吗，求赐教

cyx98

柠檬鬼lemon 发表于 2023-10-27 10:13
你好，我现在也要计算变形电荷密度，请问你是怎么做的呀，在做非自洽计算是，是读取自洽计算的CHGCAR然后 ...

https://blog.shishiruqi.com/2019/07/12/chgdiff/
根据描述，
（1）自洽计算优化分子 LCHARG=.TRUE.
（2）优化好的分子CONTCAR改为POSCAR，设置ICHARG = 12，NSW=0，LCHARG=.TRUE.
（3）将（1）电荷密度 -（2）电荷密度

柠檬鬼lemon

cyx98 发表于 2023-10-27 14:40
https://blog.shishiruqi.com/2019/07/12/chgdiff/
根据描述，
（1）自洽计算优化分子 LCHARG=.TRUE.

好的好的，谢谢，就是ISTART =0 ,ICHARG=12然后KPOINTS还和之前做自洽计算的时候一样就可以了是嘛

cyx98

柠檬鬼lemon 发表于 2023-10-27 16:37
好的好的，谢谢，就是ISTART =0 ,ICHARG=12然后KPOINTS还和之前做自洽计算的时候一样就可以了是嘛

是的