VASP计算态密度 — 从总态密度到分轨道投影的电子结构解析

态密度（Density of States, DOS）是DFT计算中最基础也最有用的电子结构信息。能带结构告诉你电子在k空间的行为，态密度告诉你电子在能量空间的分布——两者互补。而VASP的分波态密度（PDOS）更进一步：它可以告诉你某个能量区间的主要轨道成分，从而直接揭示化学键的电子结构本质。但态密度计算看似简单（改几个INCAR标签就行），实际从参数设置到数据解读，每个环节都藏着容易犯错的细节。

一、VASP态密度计算的关键参数设置

VASP计算态密度的核心INCAR参数就几个，但组合不当会让结果严重失真：

ISMEAR的选择。这是最关键的参数。对半导体和绝缘体，ISMEAR=-5（四面体方法+Blöchl修正）是不二之选——它基于k空间的四面体插值，不引入任何人工展宽，得到的DOS最接近真实态密度。对金属体系，四面体方法同样适用，但需要确保k点足够密（金属的能带在Fermi面附近变化剧烈）。ISMEAR=0（Gaussian展宽）在需要平滑DOS曲线时可用，但要记住SIGMA越小越接近真实，典型的金属体系取SIGMA=0.05-0.1 eV。

LORBIT。LORBIT=11输出原子分辨的球谐投影——每个原子的s、p_y、p_z、p_x、d_xy、d_yz、d_z2、d_xz、d_x2分量分别输出。这是做PDOS分析的必需参数。注意LORBIT=11和LORBIT=10的区别：11额外输出lm分解的投影，适合需要精确分析轨道杂化细节的场景；10只输s/p/d总投影，对大多数用途已足够。

NEDOS。态密度的能量网格点数，默认301对定性分析够用，但对精确的d带中心积分或Fermi能级附近的精细结构分析，建议2000-3000。能量分辨率 = (E_max – E_min) / NEDOS，默认约0.05 eV/点，NEDOS=3000时约0.005 eV/点——这对d带中心的积分精度有0.01-0.02 eV的影响。

EMIN/EMAX。DOS的能量范围。默认范围通常足够（覆盖所有能带），但如果你只关心Fermi能级附近的行为（如带隙分析），缩小范围可以节省DOSCAR文件大小。

二、两步法：为什么DOS计算要独立跑

VASP计算态密度的标准做法是”两步法”：

第一步：静态自洽计算，用常规k点网格，获得收敛的电荷密度CHGCAR。 第二步：读入CHGCAR做非自洽DOS计算，使用更密的k点网格（通常加密2-3倍），设置ICHARG=11。

两步法的意义在于：自洽计算需要优化电荷密度，k点密度影响迭代收敛行为，一般不需要特别高的k点；但DOS计算直接依赖能带在k空间的插值密度，k点越密，DOS越精确。分开做可以在不增加自洽迭代成本的前提下获得高精度DOS。

注意：ICHARG=11表示从CHGCAR读入电荷密度但不再自洽更新，这是非自洽计算模式。如果设成ICHARG=1（从CHGCAR读入电荷密度后继续自洽），DOS计算会额外耗十几小时的迭代时间。

三、态密度的物理解读：从总DOS到轨道投影

拿到DOSCAR后，第一步是画出总DOS图。总DOS的结构特征可以提供丰富的物理信息：

Fermi能级附近的高DOS：表明该能量处有大量可用的电子态。对于金属催化剂，Fermi能级处的DOS值与催化活性直接相关——更多的电子态意味着更强的电子给受能力。

带隙区域的零DOS：半导体和绝缘体的特征。PBE给出的带隙通常低于实验值，这需要心里有数。

尖锐的DOS峰（van Hove奇点）：来源于能带在特定k点的平带（色散很小）。在低维材料（如石墨烯）中，van Hove奇点非常明显，且与光学吸收峰的位置直接对应。

但总DOS只是全貌，真正的洞察来自分波态密度（PDOS）。以TiO₂为例：价带顶（-6到0 eV）主要由O 2p轨道贡献，导带底（2到5 eV）主要由Ti 3d轨道贡献。这个简单的轨道分析直接揭示了TiO₂的光催化激发本质是O 2p→Ti 3d的电荷转移跃迁。

四、轨道杂化分析：PDOS揭示的化学键信息

VASP计算态密度的PDOS分析是判断化学键类型的有力工具。判断规则很简单：如果两个原子的轨道在同一个能量区间出现明显的PDOS峰重叠，说明这两个轨道之间存在杂化——即化学键的电子结构基础。

以钙钛矿MAPbI₃为例：Pb 6s和I 5p在-5到0 eV区间有显著的PDOS重叠，这是Pb-I键的共价成分来源；Pb 6s的孤对电子在价带顶附近形成窄峰（约-1到0 eV），这个孤对电子的立体化学活性是钙钛矿结构软化和缺陷容忍的电子根源。

再以CO在Pt(111)表面吸附的经典体系为例：CO的5σ和2π分子轨道与Pt的d带在-8到-5 eV和Fermi能级附近分别发生杂化。5σ-Pt_d的成键态被完全占据（在-8到-7 eV形成宽峰），2π-Pt_d的反键态部分被推到Fermi能级以上。PDOS图上这两个特征峰的位置和填充程度，就是Blyholder模型在DFT层面的定量呈现。

五、态密度与实验的可比性分析

DOS是电子结构的理论输出，但如何和实验对比？最常见的是与光电子能谱（XPS/UPS）和X射线吸收谱（XAS）的对应。价带XPS/UPS测量的是占据态的态密度（加上光电离截面权重和仪器展宽），XAS测量的是从芯能级到非占据态的跃迁——反映的是非占据态的PDOS（加上芯能级选择规则）。

对比时需要做的处理：将DFT的PDOS按对应轨道的原子光电离截面加权（这个截面数据可从文献中获取），再加Gaussian展宽模拟仪器的能量分辨率（通常0.3-0.5 eV），最后将能量轴对齐（DFT的Fermi能级=实验的零点）。

一个常见的错误是直接把裸PDOS和实验XPS谱对比——没有加截面权重和展宽，峰值相对强度完全不对，导致错误的轨道归属结论。

六、VASP态密度计算的常见翻车与修复

一、ISMEAR=-5的负DOS。四面体方法在k点不够密时会在Fermi能级附近产生负的DOS值——这是插值方法的数值伪影，不是物理信号。修复方法：加密k点（让四面体更小）或换ISMEAR=0配合理SIGMA。

二、DOSCAR中PDOS的排列顺序。不同LORBIT设置下PDOS的排列不同。LORBIT=10输出s/p/d三个总量；LORBIT=11输出s, p_y, p_z, p_x, d_xy, d_yz, d_z2, d_xz, d_x2九个分量。要注意p_x/p_y/p_z的顺序和d轨道的标识——在四方晶系和六方晶系中，d轨道的”标准命名”会因坐标系选择而不同。

三、自旋极化的DOS。ISPIN=2时，DOSCAR输出的能量坐标对自旋上和自旋下是共用的，但态密度数据分开。做积分时确保两个自旋通道都包括了。

四、Fermi能级的定义。DOSCAR的能量零点就是计算中的Fermi能级。对半导体，这个Fermi能级通常在带隙中央（这是VASP的默认设置，而非掺杂后的真实Fermi能级位置）——所以DOS图上的带隙可能和”价带顶到导带底”对应不上，需要手动确认。

配图建议：

1. TiO₂总态密度与分波态密度图（ALT：”VASP计算态密度TiO2总DOS与Ti 3d O 2p分波投影态密度带隙分析”）：双层图，上为总DOS标注带隙，下为Ti 3d和O 2p的PDOS分色显示。
2. CO/Pt(111)吸附体系PDOS图（ALT：”VASP态密度计算CO在Pt(111)表面吸附的分子轨道与Pt d带投影PDOS杂化分析”）：展示气相CO的分子轨道PDOS和吸附后CO与Pt d的杂化PDOS对比。
3. 计算DOS与实验XPS对比图（ALT：”VASP态密度PDOS加截面权重展宽后与实验价带XPS谱对比”）：上层为DFT PDOS（加权重+展宽），下层为实验XPS谱。