使用VASP计算，很多时候都逃不掉DOS，能带计算的相关问题，尤其是对于计算材料的童鞋们，更是家常便饭一般。群里很多人，很多新手们都时常在讨论DOS的计算。这里我们通过VASP官网的说明，解释一下算DOS的具体步骤。前面我们学会了如何拟合或者优化稳定的晶胞结构。在此基础上，我们可以计算一下相关的DOS信息。

1.1 K点数目

与结构优化相比，算DOS的时候，需要用到更多的K点数目，这是因为K点越多，画出来的DOS图质量越高。

引用官网的话：

1.2 K点数目的选取

K点数目越多越好，我们该如何设置K点数目呢？

还记的前面我们讲到的K点选择的经验规则吗？那一个规则可以认为是我们平时计算时K点选择的标配。对于DOS计算，我们就需要把配置提高一个档次了。

一般来说，K * a = 45左右之间完全可以满足你的要求，大伙可以根据这个经验来选择K点。

NEDOS这个参数在DOS图的质量上面也有着很重要的作用。比如我们的DOS能量区间范围（DOS图的横坐标）为：[-10 eV，10eV]，VASP默认的将这个能量范围分成301点，然后作图。301也就是默认的NEDOS的取值。如果我们设置的NEDOS值够大，那么DOS区间就会被区分地越精确。NEDOS的取值一般来说：

NEDOS = 3000左右就足够好了。太大也没什么意义；

NEDOS越大，VASP输出的DOSCAR，vasprun.xml文件也就越大，占用存储空间。

经常有人抱怨说自己的DOS图有很多尖锐的峰，可以尝试着通过增加NEDOS这个办法来解决。

更多的信息，自己参考一下：https://cms.mpi.univie.ac.at/wiki/index.php/NEDOS

首先我们看官网的话： https://cms.mpi.univie.ac.at/wiki/index.php/ISMEAR

也就是说 ISMEAR = -5 的时候（Blöchl修正的四面体方法），我们可以得到一个非常平滑的DOS图。

注意：

3.1 K点数目

设置ISMEAR = -5 的时候，如果K点数目K点的数目小于等于4 ， 计算会出错，得到如下的错误结果：

VERY BAD NEWS! internal error in subroutineIBZKPT:

Tetrahedron method fails for NKPT