DOSCAR

画DOSCAR

关于状态密度的基础知识：Density of States - Engineering LibreTexts​​​​​​

                                            电子态密度如何分析？ - 知乎 (zhihu.com)

                                            态密度图分析 - 简书 (jianshu.com)

画出准确的DOSCAR要先进行静态计算，首先复制 CONTCAR 到 POSCAR 之后继续一步静态计算 (ISTART=1; NSW=0) 。

输入——相关的控制参数有：

NBLOCK：每多少离子步计算一次XDATCAR和DOSCAR，默认值为1，一般设置为1，对计算速度影响较小，超过1000步可以设置为10-100。如果SMASS = -1,控制多少离子步算一次动能。

KBLOCK：默认值为NSW，在 KBLOCK*NBLOCK 离子步骤之后，平均对相关函数和 DOS 将写入PCDAT和DOSCAR，即更新一次。 

EMIN，EMAX：DOS能量范围

NEDOS:

ISMEAAR:

RWIGS:

LORBIT：小于10 的时候必须设置RWIGS才能计算PROCAR，

输出——输出的信息内容有：

DOSCAR中，信息为：

​

line1:离子数（包含空球），离子数，0或1（无/有分波dos），NCDIJ

line2:单胞体积（1，A的三次方），基矢长度（3，m），POTIM（s）

line3:初始温度TEBEG

line4:CAR

line5:体系名字

line6:Emax, Emin,NEDOS,Efermi,1.000

line7: energy  DOS(up)  DOS(down)  integrated DOS(up)  integrated DOS(down)

包含72个离子的体系，line6重复73次,后面的72次为单个离子的贡献，每列分别为

energy s-DOS(up) s-DOS(down) p-DOS(up) p-DOS(down) d-DOS(up) d-DOS(down)

PROCAR中，信息为：

​

line1: K点数量，能级数量，离子数量。

line2: K点标号，坐标和权重

后面是每条能级的标号，能量值和占据状态，接下来是各个离子的贡献，分成spd，对于 lm-resolved计算，也有px,py pz等投影。如果设置ISPIN = 2，能级会输出两遍，分别对应自旋向上和向下。

如果设置 LNONCOLLINEAR=.TRUE. 会在每条能级的后面增加额外的三个投影，和磁矩相关。

在vaspkit中，输入11 --->5之后选择原子或轨道即可计算相应的原子轨道或者总的DOS。