一：计算前的准备工作

1，得到一个稳定的弛豫过后的POSCR结构

2，对结构进行一次静态计算

二：计算态密度的参数准备

对于态密度计算INCAR参数的一些思考

1. ISTAT=1 这里我们选择1，但可以选择的还有0，在INCAR中不填，系统应该会默认为1 。他的目的是继承上一步的静态计算当中得到的WAVECAR，以减少计算的耗时。
2. ICHARG=11 这一项还可以选择 0，1，2和12 选择用11进行计算的话就会读取我们在上一步静态计算过程中得到的CHGCAR并进行态密度计算这样，可以节约时间，以便平衡加大K点等其他参数所所带来的耗时增加。
3. IBRION=-1 这里选择 -1 我们在结构优化的时候得到的体系已经是最稳定的这里选择-1就是把所有原子固定住。（思考： 为什么要固定？如果不固定是不是也能算呢？固定是为了节约计算量以节约时间吗？如果不固定也能算，那不固定对态密度峰值精度的影响大吗？）
4. NSW=0 因为不需要调整结构所以不需要弛豫步数
5. LCHCAR=.F. 数据有点大，可以输出，但一般不选择输出节约空间
6. LWAVE=.F. 数据有点大，可以输出，但一般不选择输出节约空间
7. ISMEAR=0 （疑问：我都是用了0，但我也不知道为什么）
8. SIGMA=0.1 这个参数对能带或态密度中间的间隙影响很大，越接近0能带间隙就会越大，但相应计算的耗时也会增加（所以有时也不经怀疑自己，这弄计算的纯粹就是在凑）
9. LORBIT=11 这个参数可以将每一个原子对应的轨道的贡献输出到OUTCAR里
10. ISPIN=2 这个参数用于计算磁性，慎用，开启后会态密度图会成上下分布，对于磁性体系会出现上下分裂的情况
11. NEDOS=3000 加大态密度计算的撒点数，这样可以让态密度的峰值更平顺
12. EMIN=-20 决定态密度中x轴（能量）的最小取值范围
13. EMAX=20 决定态密度中x轴（能量）的最大取值范围
    其余的值没有写在INCAR里了，系统会自动使用默认值

对于态密度计算KPOINTS取值的一些思考

取值尽量大一点

三：计算态密度需要的文件

INCAR
KPOINTS
POSCAR （上一步静态得到的CONTCAR）
POTCAR
WAVCAR （上一步静态得到的）
CHGCAR （上一步静态得到的）

四：计算态密度

将任务投入计算列队

五：处理数据

将计算得到的vasprun.xml使用p4vasp.exe处理，导出数据
后使用matlab或origin绘图

最后：留疑

静态计算的数据也可以绘制态密度图，但经过态密度加强过的峰值会更平缓，