本文转载于微信公众号，VASP学习交流
1. 基本介绍：电荷密度分布需要得到CHGCAR, CHGCAR是基于赝势计算得到的价电子电荷密度, 包含电子密度信息的格点文件，对于自旋极化体系, 计算还包含自旋电子密度,可以使用VESTA程序打开。
2. 这里我们用体相NaCl，Si作为例子来画电荷密度分布，比较二者的电荷密度的不同，从MS库导入即可（File——Import——Structures 分别找到Cu和Si），
因为是3D材料，故而三个方向都要设K点，具体参考以前的稿子写的建议，
POTCAR用服务器里的；
3. 体系优化一下，拿优化好的CONTCAR接着做静态计算，建议的INCAR参数可以看以前的稿子；
4. 静态计算：所以INCAR设置参数需要做如下修改：
NSW=0
LCHARG=.TRUE.
#其余参数不变
5. 静态计算结束，得到CHGCAR
5.1 前16行就是POSCAR的内容，也就是坐标信息，17行是空行
5.2 第18行是x, y, z三个方向上的格点密度，分别对应于NGXF，NGYF，NGZF
5.3 19行以后就是电荷密度信息了，一共有NGXF*NGYF*NGZF个值
6. 可视化CHGCAR，把得到的CHGCAR下载下来，直接拖到VESTA作图
6.1 三维电荷密度图：原子间电荷密度轮廓反应化学键信息，共价键：电荷集中在原子间；离子键：电荷集中在一侧原子上，其实电荷密度三维的一般看的比较少。倾向于看二维电荷密度图
6.2 二维电荷密度图，需要VESTA进一步处理
Utilities——2D Data Display——Slice——修改米勒指数，也可以事先选中三个原子，三个原子确定一个平面——选完面后点击Calculate。。。——点击OK，这个方法可以显示colorbar, 但是不能显示原子
原子间电荷密度梯度反应化学键信息
共价键：梯度不明显，比如图中Si的
离子键：梯度明显，比如图中NaCl的
