如何计算差分电荷密度并可视化

在文献中有时候能够看到差分电荷密度图，它显示了材料中电荷的转移情况。一般是存在于两种不同材料的交界面处。

比如一个异质结AB（如果是一个器件，我们只需要取它的中心区），我们首先需要将它转换成一个晶体材料，这在Device Studio里面可以实现，需要注意的是，一旦转换成晶体材料，我们只想要得到的是异质结交界面处的电荷转移情况，所以两边缘处也就是异质结的两端需要加上足够的真空层（一般两端各加上8~10埃就够了），以确保在转换成晶体后两端不会受到周期的影响而产生作用。

需要注意的是，再加上真空层之前，你需要记住材料的晶格常数，以及加了多少真空层和加了真空层之后的晶格常数，因为在后面我们计算完成之后需要还原。

加上真空层，生成POSCAR文件，就可以进行正常的计算了，首先就是对结构进行优化，你当然可以选择只对交界面处进行某个方向的部分优化，因为我们只关心交界面处的电荷转移情况，那么你就需要将ISIF设置成2，不改变原胞的形状和体积，然后再POSCAR里将相应的原子的相应方向设置成T，其他都是F即可（如果你采用的是SelecFive dynamics）。KPOINTS只需选取周期方向，其他都为1即可，例如x输运y周期z真空，那么k点选取1 4 1即可。

优化完成之后，比较关键，进行整体和两部分的自洽计算，整体不用说，需要强调的是两边的计算。我们不能在DS的视图里面直接删去原子然后再生成POSCAR，而最好是在整体优化后的CONTCAR里面删去对应的原子，将其修改成对应部分的POSCAR，注意不要改变晶格常数大小，然后进行自洽计算。

OK，完成以上步骤，我们已经得到了AB和A和B的CHGCAR三个电荷密度，接下来我们打开VESTA。

在VESTA里打开我们的AB电荷密度，点击Edit-Edit Data-Volumetric Data

在Isosurfaces里import导入你需要减去的A的CHGCAR

如图，选择第二项，Subtract from current data从当前数据中减去，点击OK，再对B的CHGCAR进行同样的操作，最后检查框里是否正确

没问题就可以点击OK，生成差分电荷密度图了。

这时候你看到的图可能会比较乱，那可能是由于分辨度太大或太小的问题

如图，选择等值面，可以调节level和color，以及其他设置美化图形。

由于我们之前已经记录了加上真空层前原始的晶格常数，你当然可以对它进行还原。

如图，点击File-Export Raster Image还可以保存图像到本地文件夹。