氧空位形成能怎么算？不同位置怎么比较？

说明：本文华算科技主要介绍氧空位形成能的能量定义、氧原子参照、不同氧位的比较条件，以及带电氧空位的费米能级项。

氧空位形成能的定义与能量来源

晶格里原本由氧占据的位点失去一个氧原子，形成氧空位，记作 V_O。氧位编号、晶体表面或体相位置、电荷态共同规定一个缺陷计算对象。在 La_1.85Sr_0.15CuO₄ 中，平面氧位和顶位氧位的配位不同；从这两个位置移走氧，得到的是两个独立构型。

计算中准备原子数为 N 的完整超胞，再从指定氧位移走一个氧，得到含 N－1 个原子的缺陷超胞。两套结构各自弛豫，分别取得完整晶体总能 E_host 和缺陷总能 E_def。仅做 E_def－E_host 会漏掉移走的氧：移出的氧进入外部氧原子库，能量差需计入该原子库的化学势。

对于电中性的单个氧空位，缺氧反应可记为“完整晶体 → 含 V_O⁰ 的晶体＋氧库中的一个 O”。按这个方向定义，E_f(V_O⁰)＝E_def(V_O⁰)－E_host＋μ_O，单位通常为每个空位的 eV。氧从晶格移出，所以 μ_O 前面是正号；若采用相反的原子数符号约定，通式外观会改变，实际缺氧反应仍应给出相同能量差。

若氧库选气态 O₂，零温电子能近似下可取 μ_O＝½E_tot(O₂)，于是公式中的第三项就是半个氧分子的总能。一个 O₂ 分子含两个氧原子，系数 ½ 来自原子数守恒。计算 O₂ 时，其三重态自旋和充分隔离的分子模型会影响参照能。

若同一超胞一次移走两个氧，能量式中的氧库项变为 2μ_O；除以二才得到每个空位的平均形成能。两个空位靠近时，局域弛豫和电子占据可能互相影响，平均值包含空位间相互作用，不能用它代替孤立单空位的结果。负形成能则提示给定参照下稀释缺陷模型可能失去适用性，需检查非化学计量相或结构重构。

氧化学势为何改变形成能？

μ_O 表示外界接收一个氧原子时的能量参照。氧气丰富时，氧原子库的 μ_O 较高；同一缺陷总能差加上较高的 μ_O，氧空位形成能随之升高。贫氧环境降低 μ_O，同一氧位的形成能随之下降。这种平移来自外部氧库，晶体中的缺氧构型并未因更换参照而自动改变。

富氧上限常以 O₂ 为参照；贫氧下限受到宿主与竞争相的稳定条件约束。以 Y₂Si₂O₇ 为例，固定宿主相时，2μ_Y＋2μ_Si＋7μ_O 对应宿主相每化学式单元在该温度下的自由能，Y₂SiO₅、SiO₂ 等竞争相进一步限制可选的化学势。贫氧值受共存相的热力学约束；越过稳定区域，原先指定的宿主相可能不再是平衡相。

研究温度和氧分压时，氧气参照还含热运动与熵贡献。理想气体近似下，μ_O(T,p)＝½μ(O₂)(T,p°)＋½k_BT ln(p/p°)。升高氧分压会提高 μ_O；温度效应则需连同氧气标准化学势的温度变化计算。这里的 μ(O₂) 必须与所用 DFT 总能和热化学修正采用同一参照。

有些还原气氛把氧送入 H₂O，氧库可写为 μ_O＝μ(H₂O)－μ(H₂)。这与半个 O₂ 分子的参照定义不同，同一缺陷会得到不同的绝对形成能。比较两篇论文的 eV 数值前，氧库、温度、分压及 O₂ 能量修正必须一致。气相 O₂ 的泛函误差也可能给绝对值引入系统偏移。

不同氧位和超胞的结果怎样比较？

在同一完整超胞、同一电荷态和同一 μ_O 下，比较氧位 A 与 B 的单空位，完整晶体总能和氧库项相互抵消：E_f(A)－E_f(B)＝E_def(A)－E_def(B)。这使不同氧位的相对排序主要取决于缺氧结构的弛豫总能。若一个模型是体相超胞，另一个是表面 slab，完整模型的能量和氧位环境均已改变，不能套用这条抵消关系。

Cu₂O/Cu 界面提供了具体例子：台阶界面、平整界面和界面下方氧位的邻近配位各异。移走同一种氧，金属与氧化物的局域弛豫代价仍会不同；该论文计算中，界面台阶位的氧空位形成能低于平整界面和界面下方位点。图中还给出 Cu₂O 晶格常数随缺氧比例的变化，位点差和缺氧浓度应分别记录。

周期性计算把超胞内的氧空位复制成空间中等间距排列的镜像，缺陷在模型中并非真正孤立。超胞太小，空位之间的弹性畸变、局域电子和静电响应可能穿过超胞边缘与周期镜像耦合；同一缺陷在小超胞中还代表更高的缺氧浓度。超胞尺寸影响的是 E_def－E_host 本身，换一套氧分子参照无法消除这种有限尺寸误差。

ABO₃ 钙钛矿数据库曾比较 9 原子和 79 原子模型得到的氧空位形成能。部分化合物两种尺寸的能量差明显，另一些相近。9 原子和 79 原子模型的柱高差随化合物改变；同为一个氧空位，超胞尺寸不同也对应不同的缺氧浓度和周期镜像间距。

一个可信的位点排序还取决于缺陷构型是否充分弛豫。过渡金属氧化物失氧后，多余电子可局域在邻近金属，也可出现另一种磁有序或局域结构。只从一个初始磁矩和一种电子局域方案出发，可能停在较高能的局部极小值；这会直接抬高所报的形成能。

表面模型还涉及 slab 厚度、真空层、表面终止和单面或双面造空位。完整与缺陷模型若采用不同泛函、U 值或自旋设置，两项总能便不再处于同一计算参照；k 点和截断能未收敛也会改变能量差。形成能衡量缺氧反应的热力学代价；氧在晶格里移动的速率由迁移势垒控制，低形成能本身不给出迁移速度。

带电氧空位为何还要计入费米能级？

单斜 ZrO₂ 的氧位分为三配位和四配位两类：前者与三个 Zr 邻接，后者与四个 Zr 邻接。从这两类氧位分别移走一个电中性氧原子，缺氧结构的弛豫总能可能不同。若缺陷超胞相对中性状态再失去电子，形成 V_O⁺ 或 V_O²⁺。氧位与电荷态分别规定缺陷模型；同一氧位的中性与带电空位还需按电子库化学势比较。

以价带顶 E_VBM 为能量零点，费米能级 E_F 从价带顶向上计，常用式为 E_f(V_O^q)＝E_def(q)－E_host＋μ_O＋q(E_VBM＋E_F)＋E_corr。q 为缺陷超胞的净电荷数，正值表示移走电子；E_corr 汇集带电周期镜像和电势对齐等修正，具体项随修正方案改变。

固定氧化学势后，形成能对 E_F 的斜率等于 q：中性态为水平线，二价正电态的斜率为 ＋2。ZrO₂ 图把中性三配位氧空位设为相对零点，原子氧库项在不同电荷态的差值中抵消。相对形成能图只比较该论文内部的电荷态，纵轴零点不能当成从完美晶体造空位的绝对能量。

不同电荷态的稳定区间由各自形成能的最低包络给出。同一种缺陷的两条稳定电荷态线相交，可定义热力学电荷转变能级；不同缺陷种类的曲线交点没有这个含义。实际平衡费米能级受电荷中性条件约束，施主、受主及自由载流子共同决定可取位置，不能随意挑选某个 E_F 给缺陷排名。

Y₂Si₂O₇ 和 Yb₂Si₂O₇ 的图同时改变氧库与电子库：富氧、贫氧两组曲线的纵向位置不同；V_O⁰ 的线近乎水平，V_O²⁺ 随 E_F 上升。氧间隙还出现在同一张图里，空位与间隙的原子数项符号相反，不可把后者曲线误认成另一种空位电荷态。

E_corr 的数值对高电荷态尤其敏感：带电缺陷与周期镜像发生长程库仑作用，势能零点也可能与完整晶体不同。介电响应、超胞形状和修正公式应与报告值对应，电势对齐若已含在所选修正方法中，不能再次累加。泛函引起的带隙误差还会改变可用的 E_F 范围。

静态 DFT 总能描述的是指定缺陷构型的 0 K 能量。高温平衡浓度还受振动、构型和气相熵影响，相应贡献进入形成自由能；淬火样品的氧空位数还保留退火时的生成与迁移历史。室温测得的缺氧量未必等于将室温 E_f 直接代入平衡浓度公式得到的数值。