VASP结果分析：如何画沿z轴方向的平均差分电荷密度图

启明

VASP结果分析：如何画沿z轴方向的平均差分电荷密度图

(The Averaged Electron Difference Density Plot)

Wang’s Group at NCU

2016年11月17日星期四

• 基本目的
  

如何画沿z轴方向的平均差分电荷密度图，如下图：

• 学习基本内容
  
  

•    第一用VESTA软件画三维差分密度图
•    第二搞懂电荷密度函数，差分电荷密度函数，平均电荷密度函数的物理意义
•    第三搞懂VASP中CHGCAR文件里面数据的物理含义
•    第四安装MATLAB软件
•    第五学习MATLAB软件基本操作，以及实现平均电子密度图
  

详细内容见PDF附件：
VASP结果分析：如何画沿z轴方向的平均电荷密度图.pdf (735.76 KB, 下载次数 Times of downloads: 1243)

2016-11-19 13:06 上传 Uploaded

点击下载
Click to download

【参考文献：1. http://www.uam.es/personal_pdi/ciencias/jcuevas/Talks/JC-Cuevas-DFT.pdf

2.  http://pubs.rsc.org/en/Content/ArticleLanding/2015/CP/C5CP00011D#!divAbstract

3. http://wenku.baidu.com/view/5e6f080b27284b73f24250b9.html】

tuqingyun

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

hakuna

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，数据和图应该不是一回事

Zxin_up

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

123654

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

123wo

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

PLwang

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