【转载他人】https://blog.csdn.net/weixin_43948934/article/details/135538495。一个更好用的方法，不需要matlab，直接使用vaspkit生成。

这里的x,y,z方向，是相互垂直的吗，请问。还是沿着晶格矢量方向（这个可能是倾斜的）。

您好，由于网址改版，您引用的第一个有关DFT的PDF已经失效， 如果方便的话，您可以发到我邮箱中吗？944981075@qq.com 不胜感激

如果对比不同体系，比如晶格常数不同，或者划分网格点数不同的体系（即网格点划分密度不同）差分电荷密度时，不应该直接相加来表示积分平面积分值。
1.应该按照网格的实际三维xyz数据（Bohr或者Å）和相应的电荷数据（e-/Bohr或e-/Å)进行相乘，再相加的值表示积分。
2.平面平均除以的面积也应该是实际的平面面积，而不是80*80。

tuqingyun 发表于 2017-10-17 14:32
您好，感谢您的好贴，我有几个问题想请教您：1、按照这个方法，作图的时候横坐标一开始是0~400，怎么下一张 ...

楼主这里的CHGCAR和结构似乎不是一回事
CHGCAR里电荷密度部分是一个80x80x128矩阵，也就是沿slab的x,y,z三各方向分别划分了80、80个128个格点，也就是说，CHGCAR提供了128个垂直于z轴的平面电荷密度数据，做x-y面平均电荷密度均值随z的变化，只需要每个面的电荷密度求平均即可，可以有128个值。
画图时，z轴分成了128个格点，z轴长度是已知的，你可以将128个格点，转化成128个长度刻度，和前面的128个平均电荷密度对画就ok了，楼主将z轴单位取成Bohr，那就再除以0.529。
楼主给的slab的z轴只有11.xxx，换为Bohr也就22.xxx，所以楼主画的图前面有350多个格点，后面有30几个Bohr，数据和图应该不是一回事

您好，感谢您的好贴，我有几个问题想请教您：1、按照这个方法，作图的时候横坐标一开始是0~400，怎么下一张图就变成0~30？  2、文献上的横坐标单位一般是埃，纵坐标单位一般是（e/埃），不知是如何转换的？