VASP计算中PDOS加和不等于总DOS 核心原因与误差溯源方法

说明：本文华算科技主要介绍总态密度与原子投影态密度的数学关系，区分投影空间缺失、局域基组误差和绘图口径差异，并说明怎样根据差额随能量的分布判断一组 PDOS 是否适合支撑轨道成分分析。

同一次电子结构计算得到的总 DOS 与 PDOS 共用一组本征能量，但两者统计的对象并不相同。总 DOS 记录每个能级在能量轴上的权重，PDOS 还要询问波函数有多少成分落入指定原子或轨道的投影子空间。本征能量相同，不会自动保证投影权重之和等于一，差额正是从这里产生。

把第 n 条能带、k 点权重和展宽函数写开，总态密度可记作 D(E) = Σ_nkw_kg(E − ε_nk)。这里每个归一化本征态携带一个完整的态权重。若用投影算符 P_A 标记原子 A，对应 PDOS 还乘上 ⟨ψ_nk|P_A|ψ_nk⟩。PDOS 纵坐标是能级权重经过投影后的分量，并非另一套独立能级。

若所有原子投影满足 Σ_AP_A = I，且各曲线使用相同的 k 点权重、能量网格、展宽与自旋定义，那么每个本征态的投影权重之和都是 1，逐能量相加便能复原总 DOS。等式成立的核心条件是投影算符在目标空间内构成单位算符，不是“把元素列表写全”这一项表面条件。

计算输出里的“原子 PDOS”常由有限半径的原子球、赝原子分波或一组局域轨道生成。这些对象能标记局域成分，却未必覆盖平面波表示的全部希尔伯特空间。对某个本征态定义 q_nk = Σ_A⟨ψ|P_A|ψ⟩，当 q_nk 小于 1 时，缺少的部分仍属于该本征态，只是没有归入已选投影。若局域函数互相重叠而处理方式不当，q_nk 也可能大于 1。

用同一能量网格定义差额 R(E) = D(E) − Σ_AD_A(E)，它对应各本征态的 1 − q_nk 经展宽后的叠加。差额曲线不是新的电子态，它记录当前投影定义没有复原的谱权重；出现负值时，则要检查投影重叠、双重计数或归一化。

原子球投影只积分球内或增强区内的波函数成分，原子球之间的间隙区并未自然归属于某个原子。共价键中心、层间区域、表面真空侧和较离域的导带态都可能在球外保留可观权重。此时各原子 PDOS 的和通常低于总 DOS，差额较大处对应球外权重较集中的能量范围，而不是少读了某个元素文件。

角动量分量也有截断。只输出 s、p、d 分量时，较高 l 分量、半芯态或额外径向函数可能留在求和之外；局域基若按照孤立原子的占据轨道选取，还可能漏掉固体成键需要的极化函数。“所有原子都参与求和”与“每个原子的局域基都足够”是两项独立检查，前一项通过后仍可留下能量选择性的尖锐差额。

β-Be 的自由原子基组仅含 1s 和 2s 函数，但体相金属成键含有 2p 成分。原论文中，缺少 2p 时绝对 charge spilling 约为 19%，加入 2p 后降到 1.73%，相应密度差的最大值缩小约两个数量级。这里的缺口来自局域基不能重构原波函数，其能量分布会跟缺失轨道参与成键的位置相连。

局域轨道彼此非正交时，普通投影、对偶基投影和正交化后的投影给出的权重并不完全一致。某个波函数成分可能同时落入多个重叠轨道，直接把分量当成互斥概率相加会出现过计数。ΣPDOS 高于总 DOS 并非“间隙区贡献为负”，更常指向局域函数重叠、权重定义或数据后处理中的重复累计。

charge spilling 衡量原波函数投到局域基后损失的占据态权重。图中数据集约 99% 的化合物低于 5%，仍有少量体系远高于该范围。百分数必须连同所用基组、能区与材料化学一起读取；5% 是这项数据验证所采用的分界，并非适用于所有投影方法的普适合格线。

总 DOS 常按整个晶胞给出，元素 PDOS 却可能按单个原子平均。含四个等价 O 原子的晶胞若只把一条“O 原子平均 PDOS”加一次，结果自然少四倍；反过来，把已经按元素求和的曲线再乘原子数会造成过计数。单位中的 states·eV⁻¹·cell⁻¹ 与 states·eV⁻¹·atom⁻¹ 必须对应同一计数对象，元素标签相同不足以确认口径。

自旋极化图常把 spin-down 画在横轴下方，这个负号只服务于视觉分离。若把上下自旋的有符号纵坐标直接相加，得到的是两自旋之差对应的形状，不是总态数。总 DOS 求和使用两支自旋的正权重；向下绘图的分量应取绝对值后再按同一归一化累计。

能量轴不一致时，同一个本征能级也会在逐点相减中留下差值。总 DOS 与 PDOS 若采用不同展宽宽度、四面体积分与高斯展宽混用，或分别插值到不同网格，同一个能级会变成宽度和峰高不同的曲线。比较前应使用同一零点、网格与展宽，费米能平移要施加到全部数据，不能只移动其中一组。

曲线还可能经历峰值归一化、面积归一化或按化学计量数缩放。这样的图适合比较峰位与形状，却不再保留原始态数。图5中的两套金刚石 PDOS 在统一能区内比较带心、带宽、轨道百分比和指纹相似度，每项指标都要求相同的能区与归一化定义，否则看似接近的峰形不具备同一数值含义。

化学计量、自旋、能量网格和缩放产生的差异，可以在投影本身完备时制造不相等。只有排除这些口径项，差额曲线才指向投影空间。若整条 ΣPDOS 只差一个恒定倍数，优先核对原子数和自旋；若峰位略有错开，优先核对能量零点与插值。

把原始总 DOS 与未经峰值归一化的各原子 PDOS 放到同一网格后，差额的形状比一个总积分数字包含更多信息。平滑且跨越宽能区的正差额常与原子球外的离域权重相符；只在某组窄峰处突增的差额更像缺少特定角动量或径向函数；与某元素 PDOS 成固定比例的偏差则应核对等价原子倍数。

从最低占据能级积分到 E_F，可检查占据态谱权重损失是否与 charge spilling 或程序给出的投影质量指标相呼应；对未占据区还应单独观察，因为只针对占据态定义的指标不会约束高能导带。一个较小的占据态 spilling 不能替高能区 PDOS 作担保，研究吸附激发态、XAS 末态或高能导带时尤其如此。

可接受的误差形式取决于 PDOS 将要回答的科研任务。研究费米能级附近由哪类轨道主导时，应关注该窄能窗内的残差和轨道排序；计算电子数或轨道占据时，积分归一化与投影完备程度更严格；比较构型前后的峰位移动时，两种构型必须沿用相同投影半径、局域基和能量处理，共同存在的球外权重才不会伪装成结构效应。

大规模数据比较也表明，单看 ΣPDOS 是否贴合总 DOS 过于单薄。带心与带宽检查谱权重的位置和离散程度，指纹相似度检查整段曲线形状，spilling 约束波函数投影损失。这些量分别检验不同误差，组合起来才能区分恒定缩放、局部缺峰和广泛的基组失配。

YBa₂Cu₃O_6.5 的计算把电荷投影到 muffin-tin 球内的位点对称局域轨道。论文明确指出，原子球对总体积的覆盖有限，绝对空穴数无法据此分配；在投影定义保持不变时，结构位移前后的相对占据变化仍可估计。PDOS 求和不闭合限制的是绝对归属，并不会抹去同口径下可重复的峰位移动和谱权重转移。

当 ΣPDOS 与总 DOS 不重合，记录差额出现的能区、正负号和随构型的变化，再核对投影空间、基组、原子倍数、自旋与展宽。能量局域的缺口适合追查遗漏轨道，宽能区偏低适合检查球外权重，恒定比例偏差适合检查计数口径；只有这些信号都与研究能窗相容，原子与轨道归属才有明确的物理含义。