差分电荷密度 matlab,关注：电荷密度计算及图形化处理，电声耦合常数计算

1) 电荷密度计算及图形化处理；

2) 电声耦合常数计算；

3) Raman光谱计算与分析

题记：好久没碰过计算了，一篇稿子搁置了一年，是时候重新拿起来了，希望两周内能完成；这段时间希望在朋友们的帮助下，着力解决上述三个问题

参考资料：

(一)若干问题

1、请问价电子密度等值面图与电荷密度图的区别？

电子密度等值面图是电荷密度图的一种表达方式，但通常来说喜欢用差分电荷密度表示

电荷密度的分布是三维空间上的，所以通常是表征某一个平面上的电荷密度或是差分电荷密度的分布

这样图像就会呈现x,y为想表征的平面坐标，用等值面表征电荷分布

2、VASP怎么做出价电子密度等值面图？

可以用VESTA沿想要观测的表面表征

3、分析电学性质原理一般有哪些方法？

通常是态密度，差分电荷密度，能带结构。

差分电荷密度通常是用来分析在体系得失电子的问题，比如说做掺杂，或是其他点缺陷，通常对比形成缺陷前后的差分电荷密度讨论缺陷对电荷转移的影响，进而讨论对体系稳定性影响

至于颜色就是表示了电荷聚集或偏离的程度

(二)、VASP计算PartialPartial charge density

Partial charge density计算或称为Band decomposed charge density计算，即计算特定的某个(或某些)k点和本征值(这些k点和本征值是相互对应的)所对应的本征波函数的平方(也就是电荷密度)。特别是用在 STM的计算中，以及分析特定能量范围内或能量点的化学键特征(或atomic characteristic)。电荷密度 定义

下面以计算金刚石结构Si(晶格常数为5.40 angstrom)的VBM(价带顶)为例：

一、先自洽计算得到收敛的电荷密度(CHGCAR和CHG)和波函数(W***ECAR).

%%%%%%%%%%%%%%%%%%

INCAR输入文件如下：

SYSTEM = Si

ENCUT = 400

ISTART = 0

ICHARG = 2

ISMEAR = -5

EDIFF = 1