2. 缺陷形成能的热力学定义

缺陷形成能是点缺陷热力学研究的核心量，它回答一个非常具体的问题：在给定温度、压强和化学环境下，把一个缺陷放进完美晶体需要付出多少能量代价。这个量之所以重要，是因为它直接决定缺陷浓度、掺杂可行性和材料在器件中的稳定性。要真正理解它，必须从热力学定义出发，把化学势、参考态和电荷态三个概念串起来。

一、从自由能差到形成能定义

设一个含 N 个原子的完美晶体，其总吉布斯自由能为 G_perfect。引入一个缺陷后，体系原子数变为 N + n_i（n_i 可正可负，表示相对完美晶体增加或减少的 i 原子数），总自由能为 G_defect。缺陷形成能 E^f 定义为：

E^f = G_defect − G_perfect − Σ_i n_i μ_i + q(E_F + E_V) + E_corr

其中 μ_i 是第 i 种原子的化学势，q 是缺陷电荷态，E_F 是相对于价带顶 E_V 的费米能级，E_corr 是有限尺寸修正项。这个公式的每一项都有明确的物理含义，下面逐项拆解。

二、化学势：环境变量的核心

化学势 μ_i 描述的是“从外部 reservoir 取走或放入一个 i 原子时体系自由能的变化”。在缺陷计算中，它不是一个可以随意取值的自由参数，而是由材料制备条件约束的。以 GaAs 中的 Ga 空位 V_Ga 为例，n_Ga = −1，形成能中会出现 −(−1)μ_Ga = +μ_Ga 项。μ_Ga 越大，V_Ga 越难形成，因为环境越“富镓”，越不倾向于丢掉镓原子。

化学势的取值受热力学稳定条件限制。对于 GaAs，必须满足 μ_Ga + μ_As = ΔH_f(GaAs)，其中 ΔH_f 是 GaAs 的形成焓。同时，μ_Ga ≤ μ_Ga(bulk)，μ_As ≤ μ_As(bulk)，否则会析出单质相。这两个约束把化学势限制在一个区间内：富 Ga 条件下 μ_Ga 取最大值 μ_Ga(bulk)，此时 μ_As = ΔH_f − μ_Ga(bulk)；富 As 条件下则相反。实际计算中，通常先算出各单质相的化学势，再根据目标生长条件选择端点或中间值。

三、参考态：为什么必须明确

形成能公式中的 G_perfect 和 μ_i 都依赖于参考态的选择。常见的参考态有三种：

• 单质标准态：以 0 K、标准压强下最稳定单质的能量为零点。这是第一性原理计算中最常用的参考，因为总能量可以直接从 DFT 得到。
• 气相参考态：对于 O、N 等元素，常以 O₂、N₂ 分子的能量为参考。需要注意 DFT 对 O₂ 结合能的低估，通常需要加修正。
• 化学势参考态：直接以目标化合物中该元素的化学势为参考，适用于已知生长窗口的体系。

参考态不一致会导致形成能数值相差几个电子伏特。例如，用 O₂ 分子作参考时，ZnO 中氧空位的形成能在富氧条件下可能比用 ZnO 块体参考低 1–2 eV。因此，在报告或比较形成能时，必须同时说明参考态和化学势条件。

四、不同电荷态下的能量表达

带电缺陷的形成能随费米能级线性变化，这是缺陷物理中最重要的关系之一。对于电荷态 q，形成能写为：

E^f(q, E_F) = E^f(q, E_V) + q(E_F − E_V)

其中 E^f(q, E_V) 是费米能级钉在价带顶时的形成能。q > 0 时斜率为正，费米能级越高，正电荷缺陷越难形成；q < 0 时斜率为负，费米能级越高，负电荷缺陷越容易形成。中性缺陷 q = 0，形成能不随费米能级变化。

把不同电荷态的直线画在同一张图上，就得到缺陷形成能图。两条线的交点对应电荷态转变的费米能级，称为热力学转变能级 ε(q/q′)。当 E_F 低于该能级时，q 态更稳定；高于该能级时，q′ 态更稳定。这个能级是判断缺陷在器件工作条件下以何种电荷态存在、是否贡献载流子的关键依据。

五、有限尺寸修正与计算实践

实际第一性原理计算中，带电缺陷的形成能还需要加上有限尺寸修正 E_corr。这是因为周期性超胞中带电缺陷之间存在镜像电荷相互作用，会人为地抬高或降低总能量。常用的修正方案包括 Makov-Payne 修正和 Freysoldt 修正。以 64 原子 GaAs 超胞为例，q = −1 的 V_Ga 在 2×2×2 超胞中的修正量约为 0.1–0.3 eV，不可忽略。

计算流程通常如下：先优化完美超胞得到 G_perfect；再分别构建中性、+1、−1、−2 等电荷态的缺陷超胞，优化结构并计算总能量；然后代入形成能公式，扫描 E_F 从 0 到带隙宽度，得到各电荷态的直线；最后取最低包络线，确定稳定电荷态和转变能级。需要注意的是，DFT 计算得到的带隙通常小于实验值，因此转变能级在带隙中的相对位置比绝对值更有参考意义。

六、小结与判断依据

缺陷形成能的热力学定义可以概括为：以完美晶体和外部原子 reservoir 为参考，把缺陷引入体系所引起的自由能变化，再加上电荷态与费米能级的耦合项。化学势决定“原子从哪来、到哪去”，参考态决定“零点在哪”，电荷态决定“电子从哪来、到哪去”。三者缺一不可。

在实际判断中，如果一条形成能曲线在目标费米能级范围内始终最低，该电荷态就是稳定态；如果两条曲线相交，交点就是转变能级，器件工作点跨过该能级时缺陷电荷态会发生改变，可能引起阈值电压漂移或漏电通道。掌握这套表达方式，才能把 DFT 计算出的能量数字转化为对材料行为的可靠判断。