高熵氧化物表面怎么建模？

引言

高熵氧化物做催化模型时，难点不只是“元素多”。同一个表面上，金属种类、近邻配位、氧空位位置和局域畸变会同时变化；一个 *OH 或 *OOH 看到的并不是平均组成，而是几埃范围内的具体原子环境。高熵氧化物表面建模的核心，是把多金属位点、缺陷位置和局域结构变化放到同一个局域模型里。

如果只把几个金属随机替换进一个 slab，再挑一个最低能吸附构型，得到的结果很容易带有偶然性。更合适的做法，是先确定材料的母相和价态补偿方式，再围绕表面金属位、近邻金属组合、氧空位和畸变程度做分组。这样后面讨论 OER、ORR 或 CO2RR 时，才知道活性差异来自元素本身，还是来自缺陷和局部结构。

高熵氧化物为什么先看局域环境？

多金属位点会改变什么？

高熵氧化物通常含有四种或更多金属元素，这些金属可能有不同半径、价态偏好和 M—O 键强度。即使整体晶格仍保持钙钛矿、尖晶石、萤石或岩盐骨架，表面上的一个金属位也会受到邻近金属组成影响。对吸附物来说，真正起作用的是中心金属和邻近 O/M 原子的组合，而不是化学式里的平均比例。

例如一个表面 B 位金属附近如果有更多亲氧元素，*OH 或 *O 可能更容易稳定；如果邻近金属更容易变价，吸附过程中电荷补偿也会不同。多金属位点带来的不是简单的元素叠加，而是局域配位场和电荷补偿方式的改变。

图1 高熵钙钛矿氧化物表面转化与氧空位相关表征示例，提示多金属氧化物的反应表面会随缺陷和环境变化。DOI：10.1021/acsami.4c22352

为什么不能只用平均结构？

平均结构能描述晶格常数和整体相，但不能代表每个活性位。高熵体系里，同一个名义晶面上会出现多种不等价局域环境：中心金属不同、第一近邻不同、次表层组成不同、氧空位距离不同。把这些环境平均掉，会让吸附能和电子结构的差异被抹平。

在 VASP 里，平均结构适合做起点，不适合直接作为活性结论。更有用的模型是若干个局域代表：金属 A 富集位、金属 B 富集位、混合近邻位、邻近氧空位位、远离氧空位位。每个模型只改变一个主要因素，结果才有可比性。

氧空位应该放在哪里？

空位靠近活性位会带来什么？

氧空位会降低邻近金属配位数，也会改变局域电荷。若空位紧挨吸附位，反应物可能直接落入空位附近，形成更强的 M—O、M—C 或 M—N 相互作用；若空位在次表层，它可能通过电子结构和应变间接影响表面吸附。两种情况不能混为一谈。

Evac= Eslab−O+ μO− Eslab

式(1)用于比较同一 slab 中不同氧位点的空位形成代价。μO可以采用氧富/氧贫边界，或用同一参照做相对比较。氧空位形成能回答空位是否容易出现，吸附自由能回答这个空位附近的中间体是否合适。

图2 氧空位形成能预测框架示意，可用于理解不同表面氧位点的缺陷代价差异。DOI：10.1021/acsomega.6c03286

表面空位和次表层空位能直接比吗？

不能简单直接比。表面空位会改变吸附物落点，次表层空位更多影响金属价态、局部应变和电子态。若研究的是 OER 中 *OH、*O、*OOH 的吸附，表面空位可能直接参与成键；若研究的是导电性或整体氧迁移，次表层空位也可能更重要。

实际建模可以先把空位分成三类：目标吸附位第一近邻、第二近邻、远离吸附位。每类至少保留 2–3 个不同金属近邻组合。若不同空位类型优化后都变成类似结构，说明表面有明显自发重排；若能量和吸附构型分散很大，就要保留多种候选结构，而不是只取一个最低能点。

图3 高熵氧化物表面氧空位与局域结构变化的显微证据，说明缺陷位置会影响反应表面。DOI：10.1021/acsami.4c22352

局域畸变要怎样进入模型？

畸变为什么不是小修小补？

高熵氧化物里，不同金属半径和价态偏好会让 M—O 键长分布变宽。某些金属周围的氧八面体可能拉伸，某些位置可能压缩，表面再叠加低配位和吸附物作用，局域畸变会进一步放大。吸附能差异有时不是中心元素单独决定，而是中心位点的实际键长和配位角决定。

所以优化后不能只看总能。需要记录目标金属—O 键长、吸附物—金属键长、邻近金属位移、表层起伏和空位附近配位数。局域畸变把“同一个元素位点”拆成了多个不同的真实活性环境。

图4 多主元氧化物结构与组分示例，展示多金属引入后局域结构和电子性质可能发生变化。DOI：10.1038/s41598-023-29477-0

哪些结构指标最值得保留？

第一类是几何指标，包括 M—O 平均键长、最大/最小键长、表面金属配位数、吸附物键长和吸附角度。第二类是电子指标，包括 Bader 电荷、局域磁矩、PDOS 中金属 d 态与吸附物 p 态的重叠、差分电荷密度。第三类是缺陷指标，包括 Evac、补氧能和空位迁移倾向。

这些指标最好围绕同一个局域位点记录，而不是分散在不同模型里。若某个位点 Evac低、*OH 适中、优化后 M—O 键长保持合理，就说明它有成为有效反应位的条件；若吸附能理想但结构大幅塌陷，后续就应按重构表面重新命名。

VASP 里怎样安排计算？

局域模型怎样分组？

建议先固定一个母相晶面，例如 (001)、(011) 或实验明确暴露的晶面。然后建立若干表面模型：无空位表面、表面氧空位、次表层氧空位、不同金属近邻组合。所有模型使用相同 slab 厚度、真空层、固定层数、k 点密度和泛函设置。

每个模型至少要经历三步：干净表面优化、缺陷表面优化、关键中间体吸附优化。若体系含过渡金属，要设置合理初始磁矩，并检查优化后磁矩是否异常翻转。高熵氧化物比较的关键，是让每个局域模型只改变一个主要变量。

SYSTEM = HEO_surface_VO_ads_relax
ENCUT  = 520
EDIFF  = 1E-5
EDIFFG = -0.03
ISPIN  = 2
IBRION = 2
NSW    = 180
ISIF  = 2
LDIPOL = .TRUE.
IDIPOL = 3

# 同一晶面下分别准备：clean、VO_surface、VO_subsurface、ads_on_VO

吸附和缺陷结果怎样连起来？

先用 Evac找出哪些氧位点更容易缺失，再在这些局域环境上放 *OH、*O、*OOH、*H 或 *COOH 等关键中间体。若空位形成代价很高，但吸附自由能很好，需要说明这种位点可能浓度有限；若空位容易形成但吸附过强，也不能把它直接当作最佳位。

Eads= Eslab+X− Eslab− EX

式(2)中的 Eslab必须和吸附态对应同一个局域缺陷模型。不能用无空位 slab 去扣空位吸附态，也不能把不同金属近邻的 slab 混在一起。这样算出的 Eads才真正代表“某个空位附近某个中间体”的结合强度。

图5 高熵氧化物表面在氧空位调控下的结构演化证据，提示缺陷模型需要结合最终反应表面理解。DOI：10.1021/acsami.4c22352

本章要点总结

要点一：高熵氧化物表面不能只用平均组成代表，局域金属组合决定真实吸附环境。

要点二：氧空位既要看形成代价，也要看它与关键中间体的空间距离和成键方式。

要点三：局域畸变会改变 M—O 键长、配位数和电子结构，是吸附能差异的重要来源。

要点四：VASP 比较时应统一晶面、slab 设置和参照态，按局域模型逐组分析。

#华算科技朱老师 #VASP吸附与催化计算教程 #高熵氧化物 #氧空位 #局域畸变 #表面建模 #OER