缺陷形成能全解析：从热力学定义、带电超胞修正到缺陷浓度计算

说明：本文华算科技主要介绍缺陷形成能的热力学定义、原子与电子参照、带电超胞修正，以及有限温度下由形成自由能推算缺陷浓度的方法。

缺陷形成能是什么？

在完整晶体中拿走一个氧原子，计算程序会给出缺陷超胞的总能，但这个数值尚未构成可比较的热力学量。氧原子转移到哪种环境、晶体与哪类电子态交换电荷、超胞周期复制引入多大伪相互作用，都要记入同一份能量账本。缺陷形成能衡量指定外界条件下生成某个缺陷态的能量代价，数值随原子化学势、费米能级和电荷态变化。

形成一个点缺陷可以写成平衡反应：完整宿主转化为含缺陷宿主，并与外部原子库、电子库交换粒子。反应前后必须满足元素数目守恒与总电荷守恒。这一定义让空位、间隙原子、替位和反位缺陷进入同一套计量规则，也限定了每一项能量的物理来源。

常用符号下，缺陷 D^q 的 0 K 形成能写为：

E_f(D^q) = E_tot(D^q) − E_tot(host) − Σ_in_iμ_i + q(E_F + E_VBM) + E_corr

这里规定：向缺陷超胞加入 i 类原子时 n_i > 0，从超胞移走时 n_i n_X = −1，原子库项便成为 +μ_X。有些文献采用相反的 n_i 号约定，公式与粒子计数定义必须成对读取，否则同一空位会出现相反的化学势趋势。q > 0 表示缺陷超胞向电子库放出电子，E_F 通常从价带顶 E_VBM 起算。

同一计数法可直接处理其他缺陷：加入一个 X 间隙原子时 n_X = +1，原子库项为 −μ_X；用 B 替换 A 时，n_B = +1、n_A = −1，对应 −μ_B + μ_A。化学势前的正负号由原子流向决定，不应按“空位”或“间隙”名称凭经验填写。复合缺陷只需把每种原子的净交换数列入求和。

E_tot(host) 与 E_tot(D^q) 要使用相容的泛函、赝势、截断能和 k 点精度。缺陷会驱动邻近原子重排、成键甚至改变自旋构型，因而缺陷侧总能应取目标电荷态充分弛豫后的构型。直接删除原子后得到的单点能量仍混有未释放的局域应变。

两个总能之差往往远小于各自的绝对值，宿主与缺陷超胞的数值设置一旦不一致，基组、布里渊区采样或磁态误差就会进入差值。体积是否固定也应服从研究对象：稀释体缺陷通常保持宿主晶格矢量，只放松离子；高缺陷浓度或应力控制问题可能允许晶胞改变。能量比较必须共享同一热力学约束，否则弛豫自由度差异会伪装成形成能差异。

CaZrO₃ 的 Ca、Zr 与 O 空位在弛豫前后呈现不同位移范围。图中的黑圈标出优化后原子位置，黄色球标出空位。结构弛豫能属于形成反应的一部分，它会随缺陷种类、电荷态和局域配位改变，统一套用固定修正值缺少物理依据。

原子化学势为什么受宿主稳定域约束？

μ_i 表示晶体与外界交换 i 类原子时的能量参照。它由实验环境和共存相限定，不能在数轴上任意选择。以 A_aB_bO_c 为例，宿主平衡要求 aμ_A + bμ_B + cμ_O = μ_host；每个竞争相又给出一条析出限制。全部等式与不等式的公共区域才是宿主稳定域，富氧、贫氧对应其中不同极限。

若以单质为零点写 μ_i = μ_i⁰ + Δμ_i，通常有 Δμ_i ≤ 0；取正值会使相应单质析出。气相元素还要把温度与分压写入 μ(T, p)。例如氧库常由 O₂ 给出，富氧条件对应较高的 μ_O，并非在缺陷超胞中多放几个氧原子。

单质固体参照可取稳定相每原子总能，气体参照还涉及分子能、零点能和温压修正。半局域泛函对 O₂ 等分子的键能偏差会整体移动含氧缺陷形成能，处理时可采用经验证的热化学修正。一组缺陷必须沿用同一套原子库零点，混用不同数据库或不同泛函的 μ_i 会破坏能量闭合。

化学势项让同一种结构缺陷随制备气氛改变能量排序。移走一个氧原子时 n_O = −1，E_f 中出现 +μ_O；μ_O 降低，氧空位形成能随之下降。相反，移走 Ca 或 Zr 的空位会响应各自原子库。报告形成能时应标明稳定域坐标或富、贫条件，缺少 μ_i 取值的结果没有确定的原子库参照。

图4给出 CaZrO₃ 稳定域内 A、B、X、C、D 五组条件下的空位形成能。氧空位从富氧端走向贫氧端明显降低，Ca、Zr 空位则呈另一趋势。个别带电氧空位出现负值时，当前假设的原子库与电子库会推动该缺陷生成，随后成分、费米能级或共存相发生响应。负形成能对应所设参照状态缺乏自洽平衡，而非孤立超胞可以无限释放能量。

带电缺陷为什么需要有限尺寸修正？

平面波 DFT 以周期性条件复制超胞。中性局域缺陷也会与周期镜像发生弹性和波函数重叠；带电缺陷还会产生长程库仑相互作用。程序通常加入均匀补偿背景以维持周期体系电中性，补偿背景是数值构造，真实材料中的电荷补偿仍来自载流子、其他缺陷或掺杂原子。

三维各向同性近似下，主导误差常按 q²/(εL) 衰减，其中 ε 是介电响应，L 是超胞特征长度。图5比较点电荷与有限宽度电荷在周期阵列中的能量变化，两条曲线都随 1/L 改变。增大超胞只能逐渐削弱镜像作用，一次小超胞计算无法自动得到无限稀释极限。

E_corr 的职责包括镜像电荷校正、势能零点对齐，并在相应体系中处理带填充或各向异性介电响应。修正方案要匹配体相、表面或二维模型的周期性设定；给含真空层的薄膜直接套三维立方公式，可能留下显著尺寸依赖。介电常数、超胞形状、缺陷电荷局域范围都应进入收敛检查。

势能对齐与镜像电荷项在部分方案中已经耦合，重复相加会产生过修正。可靠做法是对若干超胞尺寸比较原始形成能、各修正分量和修正后结果，并检查远离缺陷处的电势与电荷密度。中性缺陷没有 q² 项，仍可能保留弹性镜像和缺陷波函数重叠。超胞收敛应按具体缺陷逐项验证，一个体系的尺寸经验无法直接移植到另一种局域态。

费米能级项记录电子交换：形成能对 E_F 的斜率为电荷数 q。实际 E_F 由自由载流子、掺杂物和全部带电缺陷共同满足电中性方程。形成能与费米能级需要自洽求解，外加一个任意 E_F 后直接宣称缺陷占优，可能忽略补偿缺陷引起的钉扎。有限尺寸修正也不会消除泛函带隙偏差、错误局域构型或浅能级带填充误差。

0 K 形成能怎样进入有限温缺陷浓度？

常规静态 DFT 公式给出 0 K 电子能主项。实验温度下应转向形成自由能：G_f(T, p) 可在 E_f 基础上加入振动、取向、自旋、电子激发和体积功贡献。固体中的 pΔV 常较小，振动熵与多个低能构型的统计权重可能达到可观幅度，尤其适用于软晶格和高温过程。

CdTe 中 Te_i⁺¹ 有两个相近局域极小值，价差仅 18 meV；V_Te⁺² 的亚稳构型则高出约 1.8 eV。前者在室温可以访问多个构型，后者主要停留在基态。热可及构型要通过配分函数共同贡献自由能，仅保留最低点会漏掉结构熵。

稀缺陷近似下，平衡浓度可写为 c ≈ N_sitegexp[−G_f/(k_BT)]。N_site 是可用位点数，g 包含简并度。指数关系会放大几十分之一电子伏特的自由能差。图7中的 Te_i⁺¹ 在 840 K 加入各类熵贡献后，预测浓度比简化模型高约两个数量级。浓度计算应使用形成自由能，0 K 形成能适合作为起点。

该指数式采用稀释、近似独立的缺陷。浓度升高后，缺陷之间会通过弹性场、静电作用或成键发生吸引与排斥，空位团簇的自由能也不再等于单空位形成能的整数倍。此时可引入结合能、格点模型或团簇展开处理相互作用。单缺陷形成能对应稀释极限，把它外推到高度非化学计量材料时需另建统计模型。

平衡浓度还假设缺陷能在实验时间内生成、迁移和湮灭。淬火材料可能把高温缺陷群冻结到低温，低形成能而高迁移势垒的缺陷也可能长期滞留。原子缺陷数目冻结后，电子俘获与释放仍可能在测量温度重新分配电荷态。形成自由能控制平衡占比，迁移势垒控制到达平衡的时间尺度；解读实验缺陷浓度时，应明确合成温度、气氛、冷却路径和测量温度。