转自：http://muchong.com/bbs/viewthread.php?tid=9046086&fpage=1
本帖子的内容来自于上面帖子中的附件。为了方便大家比较直接的看到，而不用下载附件，故重复发帖。@Jiansh

用VESTA 做VASP差分电荷（C-A+B）图方法：
1、导入 C 的CHGCAR（拖拽到vesta）面板框中
2、Edit——>Edit data——>Volumetric Data————>Import (选择A 的CHGCAR)——>Subtract
from current data——>OK
3、Edit——>Edit data——>Volumetric Data————>Import (选择B 的CHGCAR)——>Subtract
from current data——>OK

注意：A,B,C三个文件的格点要一一对应。
（后来发现，http://muchong.com/bbs/viewthread.php?tid=4287522&fpage=1 这个帖子也不错，很详细 @沙风）
===============
转自：http://muchong.com/bbs/viewthread.php?tid=7392808&fpage=1中的 “电荷密度差的计算.doc” @gulangyu

本文中提到电荷密度差，有时也称为差分电荷密度。
这里主要是讨论在VASP中如何得到用来分析成键前后电荷转移的电荷密度差。

此时电荷密度差定义为：delta_RHO = RHO_sc - RHO_atom
其中 RHO_sc 为自洽的面电荷密度，而 RHO_atom 为相应的非自洽的面电荷密度，是由理想的原子周围电荷分布堆彻得到的，即为原子电荷密度的叠加(a superposition of atomic charge densities)。

不同晶面的 RHO_sc 可由自洽计算的CHG或CHGCAR得到；
而计算 RHO_atom 所需的 CHG 或 CHGCAR 可由下述非自洽计算得到:
仍使用原来自洽计算时的四个输入文件，但INCAR中需要设置 ICHARG=12 和 NELM=0，其他设置不变。

需要特别注意的，应保持前后两次计算（自洽和非自洽）中的 FFT-mesh 一致。因为，只有维数一样，我们才能对两个RHO作相应的矩阵相减。不过，只要按上一段提到的设置方法做就行了，无须特别增加别的设置。

获得了RHO_sc和RHO_atom的CHG或者CHGCAR文件之后，可以利用开始提到的办法，用VESTA很容易做出电荷密度差。 返回小木虫查看更多