缺陷电荷态是 VASP 里很容易混淆的一类问题。很多人看到 q = +1 或 q = −1，就以为是在给某个原子手动指定价态；实际不是。带电缺陷计算改变的是整个超胞的电子数，体系再用均匀背景电荷保持周期性边界下的电中性。缺陷电荷态不是局域价态按钮，而是总电子数、费米能级和电荷修正共同定义的热力学状态。
这类计算更适合半导体、绝缘体、氧化物缺陷和界面陷阱；对金属 slab 或强金属性催化表面，要格外谨慎，因为费米能级不在带隙中移动，背景电荷还会强烈影响真空电势。处理前要先确认体系是否真的需要带电缺陷口径。
NELECT 改的是超胞总电子数
在 VASP 中，缺陷电荷态通常通过 NELECT 调整总电子数。若从中性缺陷中拿走一个电子，常记为 q = +1；若加入一个电子，常记为 q = −1。这个符号不是某个原子价态的直接读数，而是缺陷超胞相对于中性参照多了或少了电子。
周期性计算不能让整个超胞无限带净电，程序会引入均匀背景电荷，使总体系在数学上可计算。这一步会带来长程库仑相互作用和势能零点问题，所以带电缺陷不能只改 NELECT 就结束。NELECT 只是建立电荷态模型的入口，后面还要处理背景电荷、势能对齐和有限尺寸修正。
图1 带电缺陷形成能计算流程示意，展示电荷态、修正项和形成能之间的关系。DOI：10.1038/s41598-026-47176-4
电荷态和 Bader 电荷有什么区别？
电荷态 q 是模型总电子数的定义；Bader 电荷是计算完成后按电子密度分区得到的局域电荷估计。q = +1 不代表某个原子的 Bader 电荷正好少 1e，因为电荷可能分散在缺陷附近、带边态或整个晶格中。
因此，写结果时可以说“计算了 q = 0、+1、−1 电荷态”，再用 Bader、PDOS、缺陷态分布说明额外电子或空穴主要落在哪里。不要把 Bader 数值当成电荷态定义，也不要用电荷态符号直接推断某个原子的氧化态。
带电缺陷会和电子库交换电子
带电缺陷形成能通常写成中性形成能、元素化学势项、电荷项和修正项的组合。常见形式是：
Eform(Dq) = Etot(Dq) − Etot(bulk) − Σniμi + q(EF + EVBM) + Ecorr 式(1)
这里 EF 表示相对于价带顶的费米能级位置。q > 0 的缺陷相当于向电子库交出电子，形成能随 EF 升高而升高；q F 升高而降低。不同电荷态的斜率不同，所以形成能线会相交。
图2 缺陷形成能随费米能级变化及费米能级钉扎的示意图，说明不同电荷态具有不同斜率。DOI：10.1038/srep16977
为什么费米能级不能随便取？
对半导体和绝缘体，EF 通常在 VBM 到 CBM 之间扫描。不同生长条件、掺杂条件或外部环境会让费米能级靠近价带或导带，从而改变哪一种电荷态更稳定。若只在一个固定 EF 下报形成能，很可能漏掉电荷态转变。
这也是形成能图常画成一组折线的原因。每条线代表一个电荷态，最低包络线代表在该费米能级下最稳定的缺陷状态。带电缺陷形成能不是单个数，而是随费米能级位置变化的一组热力学关系。
图3 NaCl 中空位缺陷形成能随费米能级变化的计算示例，可用于理解最低包络线和电荷补偿。DOI：10.1038/srep16977
周期性镜像会相互作用
带电超胞在周期性边界下会和自己的镜像发生库仑相互作用。超胞越小，这个假相互作用越明显；介电常数越低，电荷越局域，修正也越重要。若不做修正，不同电荷态之间的能量差可能被系统性偏移。
常见处理包括增加超胞尺寸、做势能对齐、使用 Freysoldt、Makov-Payne 或相关有限尺寸修正。具体方法取决于体系维度、介电响应和缺陷电荷局域程度。二维 slab 和界面模型还要额外注意真空层、电偶极和面外电势平台。
图4 磁铁矿中不同 Fe 空位电荷态的形成能示例，显示电荷态和局域环境会同时影响缺陷稳定性。DOI：10.1021/acs.jpcc.5c02852
金属 slab 为什么要特别谨慎？
金属体系没有明确带隙，费米能级被金属态钉住。用 NELECT 给金属 slab 加电荷时，更像是在模拟表面充电或电场效应，而不是半导体缺陷的离散电荷态。此时形成能随 EF 扫描的意义会变弱，背景电荷还会改变真空电势和功函数。
如果目标是电化学表面带电，更合适的说法是带电 slab 或固定电子数模型；如果目标是半导体缺陷热力学，才更适合 q = +1、0、−1 电荷态和电荷转变能级。带电缺陷和带电表面不是同一个问题，不能只因为都改 NELECT 就共用解释。
交点代表稳定电荷态切换
当两个电荷态形成能相等时，对应的 EF 位置就是电荷转变能级。例如 q = 0 和 q = −1 的交点说明费米能级高于该位置后，负电荷态更稳定；低于该位置时，中性态更稳定。这个量可以和带隙中的陷阱能级、载流子俘获或界面态联系起来。
ε(q/q’) = [Eform(Dq; EF=0) − Eform(Dq’; EF=0)] / (q’ − q) 式(2)
转变能级不是 DOS 里某个峰的位置。DOS 缺陷态展示电子能级分布，热力学转变能级来自不同电荷态总能差。两者可以互相参考，但不能直接画等号。
图5 SiO2/Si 界面氢相关缺陷的电荷转变能级示例，显示不同缺陷位点可产生不同陷阱位置。DOI：10.1038/s41598-026-47176-4
至少交代四件事
第一，说明 q 的定义和 NELECT 改变量；第二，说明元素化学势边界，比如 O-rich/O-poor、S-rich/S-poor 或金属富集条件；第三，说明 EF 扫描范围和 VBM/CBM 的来源；第四，说明是否做了电荷修正、势能对齐和超胞尺寸检查。
如果只写“q = −1 更稳定”，但没有给出费米能级位置，这句话是不完整的。更好的表达是：“在 EF 靠近 CBM 时，D−1 形成能最低；在 EF 靠近 VBM 时，D0 或 D+1 更稳定。”这样才能把电荷态和材料环境联系起来。
催化体系里什么时候需要它？
若催化剂是半导体氧化物、硫化物、光催化材料或绝缘载体，带电缺陷可能影响载流子、表面电势和吸附位点。若体系是金属电极或高导电合金，通常不建议直接套用半导体带电缺陷形成能图，而应转向表面充电、恒电位近似或功函数/电势分析。
还要注意，带电缺陷最好和中性缺陷一起报告。中性缺陷提供结构和化学势基准，带电缺陷提供费米能级依赖；如果只给某一个 q 态，很难判断它是在整个带隙内都稳定，还是只在靠近 VBM 或 CBM 的窄区间稳定。对于含过渡金属的氧化物，还应同时检查局域磁矩是否随 q 改变，否则额外电子可能落到金属 d 态而不是目标缺陷态。
报告结果时，也建议把费米能级范围写清楚。横坐标通常从 VBM 到 CBM，而不是任意取一个电子能量零点；若使用实验带隙、杂化泛函带隙或校正带隙，需要说明采用哪一种带隙边界。这样别人看到斜率、交点和最低能电荷态时，才知道这些线对应的是同一套能量参照。
要点一：缺陷电荷态通过超胞总电子数定义，不等于某个原子的手动价态。
要点二：带电缺陷形成能随费米能级变化，不应只报一个无条件的数值。
要点三：背景电荷、周期性镜像和势能零点会影响能量，需要修正或至少说明边界。
要点四：半导体带电缺陷和金属带电 slab 是两类问题，不能因为都改 NELECT 就共用结论。
