态密度（Density of States, DOS）是第一性原理电子结构分析的核心工具，直观呈现电子在不同能级的分布规律，是判断材料金属/半导体属性、解析轨道成键作用、定位催化活性中心的必备手段。

一、VASP 计算标准流程

态密度计算在VASP中需要两步：第一步是常规的静态自洽计算，产生收敛的电荷密度；第二步是读取该电荷密度进行非自洽计算，在一条致密的K点网格上重新计算本征值。

这样两步走的方法平衡了精度与效率——自洽部分用常规K点即可收敛电荷密度，非自洽部分则用更密的K点来解析能带能量分布。

1.静态自洽计算（SCF）

以优化后的 CONTCAR 为初始结构，固定原子位置（NSW=0），开展高精度电子自洽迭代，生成收敛的电荷密度文件（CHGCAR）。这一步建议收紧电子收敛阈值至EDIFF=1E-6，保证电荷密度精度，为后续 DOS 计算提供可靠基础。

2.非自洽 DOS 计算

设置ICHARG=11读取上一步的电荷密度，在不更新电荷的前提下专门求解本征值计算态密度。这一步可单独加密 k 点、提升采样点数。

通常要求K点间距达到0.01–0.02 Å⁻¹（对应约1000个不可约K点以上）。对于块体硅的常规计算，8×8×8自洽搭配16×16×16甚至24×24×24非自洽是常见操作。k点过疏会出现大量虚假尖峰、带隙计算偏差大，需在算力范围内适当加密。

二、数据处理与可视化

VASP输出的DOSCAR文件结构固定：头部信息后是总DOS（能量、DOS、积分DOS），随后依次为每个原子的投影态密度。手工提取效率低，推荐用vaspkit功能21或自行编写解析脚本。

绘图前必须将能量对齐到费米能级。费米能级可从OUTCAR中 grep "E-fermi" 获取，或在DOSCAR头部第6行读出。将所有能量列减去E_Fermi，使费米能级位于0 eV。半导体的带隙直接读取价带顶到导带底的能量差即可，注意四面体方法下带边可能稍有噪声，可结合能带计算确认。

与能带计算的区别。能带计算需要沿高对称k点路径（Line-mode KPOINTS），而态密度需要均匀致密网格。两者从同一CHGCAR出发但KPOINTS完全不同。若需画能带+DOS联合图，分别进行两步非自洽计算，不要混用KPOINTS。

三、常见错误与避坑

1.忘记第一步自洽计算就直接做DOS：非自洽需要CHGCAR，没有则读取初始电荷密度，结果不可靠。

2.非自洽K点沿用自洽的稀疏网格：DOS毛刺多，峰值位置偏移，这是最频繁出现的问题。

3.金属用ISMEAR=-5导致负DOS：换用ISMEAR=0并给合理SIGMA，或尝试-15。

4.展宽值未在图中标注：审稿人常要求说明展宽方法和数值，养成在图中或图注标明的习惯。

5.未对齐费米能级：直接绘图会得到偏移的能量坐标，失去物理意义。

【科学计算 优选至简】 

至简科技提供专业的第一性原理、分子动力学、量子化学、机器学习、有限元仿真、相图计算等代算服务。 

✅50➕全职工程师 500➕合作高校企业
✅部分合作项目己发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊