计算化学公社

标题: VASP结果分析：如何画沿z轴方向的平均差分电荷密度图 [打印本页]

作者
Author: 启明 时间: 2016-11-19 13:09
标题: VASP结果分析：如何画沿z轴方向的平均差分电荷密度图
本帖最后由启明于 2016-11-19 13:06 编辑

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

(The Averaged Electron Difference Density Plot)

Wang’s Group at NCU

2016年11月17日星期四

基本目的

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

(, 下载次数 Times of downloads: 70)

学习基本内容

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

详细内容见PDF附件：
(, 下载次数 Times of downloads: 1243)

【参考文献：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】

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

作者
Author: hakuna 时间: 2017-10-17 18:07

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

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

作者
Author: 123654 时间: 2023-10-15 20:51
您好，由于网址改版，您引用的第一个有关DFT的PDF已经失效，如果方便的话，您可以发到我邮箱中吗？944981075@qq.com 不胜感激

作者
Author: 123wo 时间: 2024-1-25 10:26
这里的x,y,z方向，是相互垂直的吗，请问。还是沿着晶格矢量方向（这个可能是倾斜的）。

作者
Author: PLwang 时间: 2025-1-9 11:21
【转载他人】https://blog.csdn.net/weixin_43948934/article/details/135538495。一个更好用的方法，不需要matlab，直接使用vaspkit生成。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)

Powered by Discuz! X3.3