说明：本文华算科技主要介绍DFT中氧空位的结构性质、形成能、电子结构、迁移能垒和吸附能。
氧空位模型是什么？
氧空位在计算模型里不是一个抽象符号，而是晶格中某个 O 原子被移除后的缺陷位点。这个位点可以写作 VO，也可以按电荷态写成 VOq。真正进入 DFT 计算时，研究者必须先指定它位于体相、表面还是界面，位于桥氧、顶位氧还是晶格内部氧，再给出超胞大小、表面层数、真空层和缺陷浓度。同样叫氧空位，位点不同，近邻金属配位、局域应变和电子占据都会变成另一套结果。
图 1 里的两个构型提醒我们，DFT 能计算的并不是“材料里有没有空位”这句话本身，而是某个构型在给定计算条件下的能量、几何和电子态。轴向 VO 与平面 VO 对 TiO6 八面体的破坏方式不同，缺氧后留下的电子也会分配到不同的 Ti 近邻或空位附近态。写氧空位性质时，缺陷位点、缺陷浓度、参考结构必须和后面的能量、DOS、吸附结果一起出现。
形成能是最常被放在第一位的量。中性氧空位的基本比较可以写成缺陷超胞能量减去完整超胞能量，再加回一个氧化学势；带电缺陷还会出现 qEF、价带顶参照和有限尺寸修正。它回答的是在氧富、氧贫、不同费米能级条件下，某类 VO 的生成代价有多高。Eform低，只说明热力学生成更容易；它不自动给出迁移速度、吸附强弱或催化活性。这些量要从后续结构、电子和反应模型里分别计算。
结构性质能给出哪些局域畸变信息？
移走 O2- 后，原来由氧桥连的 M-O-M 键断开，金属原子的配位数下降，邻近金属可能向空位内侧靠近，也可能因为电子还原和离子半径变化向外松弛。DFT 优化后的坐标可以给出 M-O 键长、M-O-M 键角、金属-金属距离、八面体倾斜角和表面台阶重构。这些几何量会改变后面 DOS 峰位、局域磁矩和吸附位点。
CeO2(110) 的还原表面给了一个很直观的例子：VO 附近会出现 nO、nCe、nnCe 等不同近邻位置，图中 Ce-O 距离已经不是完整表面的均一配位。对于氧化物来说，Ce4+/Ce3+、Ti4+/Ti3+、Fe3+/Fe2+ 这类还原对会改变离子半径，结构弛豫便和电荷局域绑在一起。只报“出现氧空位”太粗；更有用的写法是列出 VO 一近邻和二近邻的键长变化，说明哪一圈原子承担了主要畸变。
M-O-M 夹角变小后，d-p-d 轨道重叠减弱，带宽可能收窄；表面 O-VO-M 构型让吸附分子弯折时，分子的反键轨道占据可能增加。晶格畸变还可以用 原子位移、局部体积、八面体畸变参数或径向分布来量化。键长、键角、位移和配位数合在一起，才构成氧空位的结构性质。
图 3 把结构信息推进到迁移路径附近的配位壳层。随着半径 R 扩大，被统计的 X-O* 配对数快速增加，而 ICOHP 总和在一定范围后趋于平缓。这类图可以帮助判断结构分析需要多大的邻域：只看最近邻，容易漏掉迁移氧和二近邻金属之间的键合贡献；邻域取到能量量基本收敛后，能垒来源可以对应到具体键合项。这里的重点不是把 ICOHP 当成万能指标，而是用它确认“哪一圈键”正在参与 VO 相关过程。
图 4 中 CO2 靠近 VO 后不再保持线性，O-C-O 夹角从气相接近 180° 弯到 136.9° 或 140.4°，C-O 键长也出现一长一短的分化。这些几何量比单独的吸附能更能说明吸附状态：分子是否被活化、哪个 O 与空位附近 Ce 配位、电子是否进入 CO2 反键轨道，都要从键长和夹角里读。Eads给出吸附前后能量差，结构参数给出吸附后分子处在什么化学状态。
电子结构性质怎样区分缺陷态/电荷重排
氧空位移走的是一个带负电的晶格氧，计算里常会留下两个电子。它们可能局域在相邻金属 d 或 f 轨道上，形成 Ti3+、Ce3+ 这样的还原中心；也可能在空位附近形成较分散的缺陷态。DOS、PDOS、能带结构、差分电荷密度、Bader 电荷、局域磁矩和 ELF 都可以用来追踪电荷重新分布，但每一种图回答的对象不同。DOS 是能量分布，能带是 k 空间色散，电荷密度是空间分布，Bader 是分区电子数。
图 5 把氧空位的电子态分成三种读法。能带图中被标出的低能带显示缺陷态在布里渊区里的色散；右侧 DOS 把 VO、近邻 Ti 和远邻 Ti 的贡献分开；下方实空间图显示空位相关轨道集中在缺陷周围。若只给一张总 DOS，读者只能看到费米能级附近多了峰；加入投影和实空间图，才能判断这个峰来自空位腔、金属 d 轨道，还是两者杂化。
在 CeO2 这类强局域体系中，氧空位常和 Ce f 态有关。图 6 的自旋分辨 PEDOS 可以读出费米能级附近 Ce f 态的位置，以及 O p 态在价带区的分布。这里常见的适用范围是：DOS 峰只能说明某类轨道在某个能区有贡献，它不能替代空间电荷图来说明电子具体集中在哪个 Ce 原子附近。如果要确认 Ce3+ 位点，还要检查自旋密度、局域磁矩或 Bader 电荷。
图 7 里蓝色等值面集中在两个 nCe 位点附近，说明还原电子并非均匀散在整个表面。写电荷密度时要区分三类图：总电荷密度显示电子云背景，差分电荷密度显示缺陷或吸附前后的电子增减，过剩自旋密度显示未成对电子位置。Bader 电荷给出每个原子分区内电子数的变化，但分区方法会影响绝对值，更适合比较同一套模型中的相对变化。
在电子结构图里，ELF 更接近局域电子对和成键区域的图像。图 8 中 CO2 周围的高 ELF 区域和电荷转移量，可以辅助判断吸附后分子是否获得电子、C-O 键是否削弱。它和差分电荷密度的读法不同：差分图看电子增减，ELF 看局域化程度。若讨论“氧空位活化 CO2”，较扎实的证据组合通常是弯曲构型、C-O 键长变化、Eads、PDOS 轨道重叠、电荷转移和 ELF 局域图一起出现。
能量性质怎样连接形成、迁移和吸附？
氧空位的能量性质至少分三类。Eform 衡量从完整晶格产生 VO 的代价，受氧化学势、费米能级、电荷态和缺陷浓度影响；Emig 衡量氧原子或空位在相邻位点之间跳迁的势垒，通常来自 NEB 或类似最小能量路径；Eads 衡量分子吸附到缺陷表面前后的能量差，参考态要写清楚是气相分子、洁净缺陷表面还是完整表面。三者的参考态不同，数值不能混用成一个“氧空位活性”。
图 9 的横轴是反应坐标，纵轴是相对能量。初态到最高点的能量差就是迁移势垒 Emig。如果 VO 形成能低但迁移势垒高，材料可能容易生成缺陷，却不容易让缺陷快速扩散；若形成能和迁移能垒都低，缺陷浓度和扩散都可能更活跃。讨论离子导体、还原气氛下的氧迁移、催化循环中的晶格氧补位时，Emig 比 DOS 峰更接近动力学问题。
把 CO2、O2、H2O、OH*、O* 这类物种放到完整表面和含 VO 表面上比较，可以判断空位是否稳定某个中间体。吸附能更负通常表示吸附后能量降低更多，但仍要检查构型是否合理、分子是否过度束缚、后续步骤是否受阻。氧空位的催化影响往往来自形成能、吸附能和迁移能垒之间的平衡。
实际阅读一篇氧空位 DFT 文章，可以按对象分表：结构栏写 VO 位点、键长、键角、配位数和畸变；热力学栏写 Eform 与化学势条件；电子栏写 DOS、能带、差分电荷密度、Bader、磁矩和 ELF；动力学栏写 Emig、路径和过渡态；表面反应栏写 Eads 与吸附构型。这样读，氧空位就不再只是论文里的一个符号，而是一组可以彼此校验的结构、电子和能量结果。
