缺陷形成能的有限尺寸修正怎么处理？

引言

带电缺陷形成能常见的偏差来自周期性超胞。VASP 用周期边界求解库仑相互作用，缺陷电荷会和自己的周期镜像、均匀补偿背景以及能带边参考同时发生耦合。有限尺寸修正的核心任务，是把不同超胞和不同电荷态放回同一能量标尺。若只把一个 q 态的总能代入形成能公式，结果看似完整，实际仍可能含有明显的超胞尺寸误差。

这个问题在半导体和绝缘体缺陷里尤其明显。氧空位、碘空位、金属空位、间隙原子和替位掺杂都可能带电；q 的绝对值越大、介电常数越低、超胞越小，镜像电荷项越需要处理。本文以 ZnO 缺陷修正算例为主线，说明形成能、转变能级和 VASP 后处理记录怎样衔接。

为什么带电缺陷会有超胞误差？

带电缺陷的超胞误差首先体现在同一 q 态随模型尺寸变化的形成能漂移。图1以 ZnO 中的氧空位和八面体间隙锌为例，左右横轴分别为 L⁻¹ 和 L⁻³，L 表示平均缺陷间距；实心、空心符号分别表示 Madelung 修正前、后的结果。修正后仍可能保留近似 L⁻³ 的尺寸依赖，需要继续外推到无限大超胞极限。形成能是否随超胞尺寸收敛，是判断修正是否足够的直接入口。不同缺陷和电荷态各有自己的稀释极限，不能把曲线接近解释为所有缺陷趋于相同形成能。

图1 ZnO 氧空位和八面体间隙锌的形成能尺寸依赖。左右横轴分别为 L⁻¹ 和 L⁻³，实心与空心符号分别表示 Madelung 修正前后；条件为贫氧极限、费米能级位于价带顶。DOI: 10.1103/PhysRevB.77.245202。

形成能写作式(1)。E_tot(D^q) 和 E_tot(bulk) 来自 VASP 总能；n_iμ_i 处理缺失或加入的原子；q(E_F+E_VBM+ΔV) 把费米能级与电势对齐接入；E_corr 处理有限尺寸误差。E_corr 不是额外调参，而是形成能公式中和带电边界条件对应的一项。式(1)把电势对齐与镜像电荷修正分列；若后处理程序输出的 FNV 总修正已包含类对齐项，应将其整体计入 E_corr，不再额外叠加同一来源的 qΔV。

E_f(D^q) = E_tot(D^q) – E_tot(bulk) – Σn_iμ_i + q(E_F + E_VBM + ΔV) + E_corr　　式(1)

中性缺陷 q=0 时，长程静电修正项通常不进入；带电缺陷 q=+1、+2、-1、-2 时，修正项会直接改变折线截距。若比较 V_O²⁺ 与 V_Zn^2-，两者不仅缺陷类型不同，q 的符号、局域弛豫、介电响应和补偿背景也不同。所有 q 态应使用同一套体相参照、赝势、泛函、超胞矩阵和带边定位方式。

VASP 计算阶段不能把修正留到最后才想。带电超胞需要在 INCAR 中通过 NELECT 改变电子数，并保留相同晶格、相同 k 点、相同 POTCAR 与相同能量读取口径。二维 slab 或表面缺陷还涉及真空方向和偶极修正，体相 FNV 修正公式不能原样替代低维模型的电势处理。

电势对齐和介电常数怎样进入修正？

缺陷附近的原子弛豫会使电势出现明显局域起伏，电势对齐需要在远离缺陷的区域取样。图2比较立方 BN 中硼空位及其他缺陷的平面平均电势和原子位点电势，说明应先扣除模型点电荷的长程势，再检查剩余电势差的远场行为。电势对齐 ΔV 用于校准缺陷超胞和完美超胞的静电势参考。对弛豫后平面平均电势起伏较大的体系，可用远场原子位点电势估计对齐常数，同时记录取样范围和离散程度。

图2 缺陷诱导电势、模型点电荷势及其差值的比较。立方 BN、MgO 和 Si 算例展示平面平均与原子位点取样的区别；远场取样区域用于估计剩余电势差。DOI: 10.1103/PhysRevB.89.195205。

修正常用两类信息：介电常数 ε 和远离缺陷区域的电势差。ε 可以来自文献、实验或 DFPT 计算；各向异性材料应说明是否采用张量或有效平均值。对于极性半导体，电子介电响应和离子介电响应的选择要与缺陷弛豫条件一致。只修正电子极化或包含离子松弛，给出的 E_corr 会对应不同物理边界。

式(2)写出常见的 Madelung 型主项。α 是与超胞形状有关的常数，L 是超胞特征长度。这个表达式能说明尺寸误差为何随 q² 增大，但实际 FNV 或 sxdefectalign 等流程还会处理模型电荷势、短程对齐和各向异性屏蔽。

E_MP ≈ q²α / (2εL)　　式(2)

局域结构弛豫会改变缺陷电荷的空间分布。图3比较完美 ZnO、不同电荷态氧空位和中性间隙锌的能带，并用下方波函数模平方展示箭头所指电子态的空间分布；氧空位的局域态与间隙锌的延展态具有不同特征。修正前要先确认缺陷态是否局域在目标缺陷附近。若额外电子占据导带延展态，还要检查带填充及尺寸依赖，不能仅凭 DOS 峰的位置判断局域性。

图3 ZnO 完美晶体、氧空位和间隙锌的能带及选定电子态的波函数模平方。HSE 混合泛函的交换比例 a=0.375，下方空间分布用于区分局域态与延展态。DOI: 10.1103/PhysRevB.77.245202。

VASP 后处理时，LOCPOT 可用于提取平面平均或球平均静电势，vasprun.xml 和 OUTCAR 用于记录总能、费米能级和最终电荷状态。Bader 电荷、PARCHG 和局域磁矩不是修正公式的替代项，但能判断多余电子或空穴是否留在缺陷附近。若电子被远端表面态或另一个缺陷吸收，形成能曲线应回到电子态设置重新检查。

转变能级怎样受修正项影响？

热力学转变能级来自稳定形成能折线的交点，交点位置取决于两个 q 态的截距和斜率。图4以 β-Li₂TiO₃ 和六方 BN 为例，比较不同超胞尺寸、形状与屏蔽模型下的相对形成能，说明各向同性近似与各向异性修正可能得到不同的能量截距。不同 q 态的 E_corr 会改变折线交点，因此会改变转变能级位置。图4比较的是同一缺陷电荷态的尺寸收敛；计算转变能级时，还需分别取得两个电荷态在同一修正约定下的收敛形成能。

图4 β-Li₂TiO₃ 中钛空位与六方 BN 中反位缺陷的各向异性修正。结构、电势及形成能尺寸收敛相互对照，能量零点取最大超胞的各向异性 FNV 修正结果。DOI: 10.1103/PhysRevB.89.195205。

转变能级可由式(3)计算。q 和 q′ 是两个电荷态，E_f(D^q; E_F=0) 表示费米能级取 VBM 时的形成能截距。若两个 q 态使用不同超胞尺寸、不同电势对齐方式或不同带隙修正，ε(q/q′) 就不再是同一能量标尺上的交点。

ε(q/q′) = [E_f(D^q;0) – E_f(D^q′;0)] / (q′ – q)　　式(3)

以 V_O 为例，+2、+1、0 三个电荷态可能对应不同局域弛豫。若 V_O²⁺ 的结构来自带电优化，而 V_O⁰ 沿用了中性初猜但未充分弛豫，形成能交点会同时混入修正误差和结构误差。对浅能级缺陷，带填充修正、带边修正和超胞尺寸收敛要一起记录；对深能级缺陷，PARCHG 可帮助确认态的空间局域性。

若使用 HSE06 或 GGA+U 校正带隙，还要分清两种计算层级。PBE 大超胞加修正常用于筛选趋势，HSE06 小超胞用于带边和缺陷态校准；两者合并时应说明能量对齐方式。带隙修正不能自动替代带电缺陷有限尺寸修正，因为一个处理能带位置，一个处理周期性电荷相互作用。

VASP 任务和结果记录怎样安排？

推荐把每个缺陷按 defect/q/supercell 三个维度分开。例如 V_O/q+2/2x2x2、V_O/q+2/3x3x2、V_O/q0/2x2x2 各自保存端点优化、静态单点、LOCPOT 和后处理结果。NELECT、总能、VBM、E_F、ΔV、E_corr、ε、超胞尺寸和最终局域结构都要记录在同一矩阵中。

SYSTEM = charged_defect_static

ENCUT = 520

EDIFF = 1E-6

ISMEAR = 0

SIGMA = 0.05

ISPIN = 2

LASPH = .TRUE.

LVTOT = .TRUE.

LVHAR = .TRUE.

NELECT = neutral_supercell_electrons – q

ISYM = 0

上面只是静态计算参考示例。ENCUT、k 点、超胞大小和 ε 的来源都需按材料测试；含 d/f 电子体系还要同步说明 U 值、MAGMOM 和 POTCAR 版本。若比较 V_O、V_Zn 与 Zn_i，各缺陷的原子化学势边界也要保持同一相稳定区，不能把 O-rich 与 Zn-rich 结果混在同一默认结论里。

比较多种 ZnO 缺陷时，电荷态、费米能级范围和化学势条件需要同时明确。图5上方分别为贫氧和富氧极限下的形成能，下方汇总稳定电荷态之间的转变能级；图中还包含氢杂质，不能把所有曲线都称为本征缺陷。有限尺寸修正后的形成能图才适合继续讨论缺陷浓度、施主受主行为和迁移路径。若后续还要算迁移势垒，同一 q 态的端点和 NEB 图像应继承这套修正口径。

图5 ZnO 缺陷和氢杂质的形成能与热力学转变能级。上方左、右面板分别为贫氧和富氧极限，下方为稳定电荷态之间的转变能级；结果包含有限尺寸修正。DOI: 10.1103/PhysRevB.77.245202。

最终记录不宜只保留修正后数值。原始形成能、ΔV、E_corr、所用介电常数、修正程序版本、超胞尺寸序列和带边来源都应保留。若修正前后排序改变，应同时报告改变发生在哪些 q 态之间，并说明该排序是否通过更大超胞复算支撑。

跨超胞核对时，可以把同一缺陷同一 q 态在 2×2×2、3×3×2 与更大超胞中的修正前后形成能同时画出。若修正后斜率仍然明显，说明介电常数、电势对齐区域或缺陷态局域性仍需检查；若修正后不同尺寸结果进入很窄区间，再讨论转变能级和热平衡浓度才有数值基础。

实际项目里还会遇到两类边界。第一类是浅施主或浅受主，额外电子或空穴可能进入带边延展态，此时需要结合带填充修正和带边对齐，不能只看局域缺陷中心。第二类是表面或界面缺陷，真空方向屏蔽和面内屏蔽不同，体相各向同性 ε 只能作为近似参照。当模型维度、缺陷局域性或带边占据发生变化时，修正方案也要随模型重新说明。

如果同一缺陷在中性态和带电态下局域构型差别很大，修正项只处理静电边界，不能消除结构采样不足。此时应把每个 q 态分别优化到局域极小值，再用相同后处理口径计算形成能；对多稳态缺陷，还要把不同局域构型的能量差纳入同一 q 态内部比较。

本章要点总结

要点：带电缺陷形成能包含周期镜像电荷、补偿背景和电势零点带来的有限尺寸误差。

要点：E_corr、ΔV、介电常数和带边位置应与同一套超胞和电荷态计算对应。

要点：转变能级来自不同 q 态形成能折线交点，修正项会直接影响交点位置。

要点：形成能图用于缺陷浓度、施主受主行为和迁移路径前，应保留原始能量、修正项和超胞收敛记录。