态密度计算
态密度：表示单位能量范围内所允许的电子数，也就是说电子在某一能量范围的分布情况。

因为原子轨道主要是以能量的高低去划分的，所以态密度图能反映出电子在各个轨道的分布情况，反映出原子与原子之间的相互作用情况，并且还可以揭示化学键的信息。

态密度有分波态密度（PDOS）和总态密度（TDOS）形式。

原则上讲，态密度可以作为能带结构的一个可视化结果。

很多分析和能带的分析结果可以一一对应，很多术语也和能带分析相通。

但是因为它更直观，因此在结果讨论中用得比能带分析更广泛一些。

计算过程：主要分成三步：一、结构优化；二、静态自洽计算；三、非自洽计算。

1，结构优化：原子弛豫，确定体系内每个原子位置。

常用INCAR。

2，静态自洽计算：（得到自洽的电荷密度CHG、CHGCAR和E-fermi，提供给下一步非自洽计算用）
INCAR设置注意，ICHARG = 2
3，非自洽计算（准确计算电荷分布）
INCAR设置：ISTART=1（若存在WAVECAR文件时取1）；ICHARG=11（表示从CHGCAR中读入电荷分布，并且在计算中保持不变）；RWIGS (或LORBIT=11（或10），这时可不设RWIGS)；
计算完成时，生成DOSCAR，采用spit_dos.dl小程序把dos分开（注意vp.dl要拷到同目录下），会生成N+1个文件，DOS0为总态密度，DOS1到DOSN为N个原子的分态密度。

每个分态密度共7列分布为
—能量→Sup→Sdown→Pup→Pdown→Dup→Ddown
不知道从态密度能否定性分析出来，因为态密度越尖，则电子的局域性越强，
修正版的splitdos有三个文件：vp、sumdos和split_dos.ksh
INCAR设置：
ISTART = 1；ICHARG = 11
LORBIT = 10
【对于PAW势,可设置LORBIT = 10，此时可不用设置RWIGS参数】或者设置RWIGS参照POTCAR。