高熵氧化物的表面不是一个“平均位点”。同一个晶面上，中心金属可能不同，第一近邻金属也不同，旁边还可能有氧空位、羟基或局域畸变。吸附物真正感受到的是一个很小的局域环境，而不是化学式里写出来的整体比例。局域位点筛选的目标，是找出哪些原子环境真正改变了吸附强弱。
这篇不再讨论高熵氧化物表面如何从零搭建，而是专门解决筛选顺序：先看中心金属，还是先看近邻元素？氧空位要不要单独分组？局域畸变是结果，还是也能作为筛选指标？把这几个问题拆开，后面算 *OH、*O、*OOH、*COOH 或 *H 才不会变成随机枚举。
中心金属先决定直接成键对象吗？
吸附物通常先和表面中心金属成键。比如 *OH、*O 常通过 O 端接近金属位，*COOH 可能同时受到金属位和邻近氧的影响。中心金属的 d 电子占据、价态偏好和 M—O 键强度，会直接影响吸附能，所以它仍然是第一层分组。
但中心金属相同，不代表吸附能相同。高熵氧化物里，同一个 Co、Ni、Fe 或 Mn 位点周围可能接着完全不同的近邻金属。中心金属回答“谁和吸附物直接成键”，近邻环境回答“这个成键位点被怎样调节”。
图1 高熵氧化物表面结构与氧空位相关表征示例，说明表面位点需要结合局域环境理解。DOI：10.1021/acsami.4c22352
为什么同一金属会出现不同吸附能？
同一中心金属周围若有更多亲氧元素，吸附含氧中间体时局域电子补偿会不同；若邻近的是半径较大的金属，表面 M—O 键长和吸附角度也可能改变。吸附能的差异，有时不是中心金属本身变了，而是它所在的配位场变了。
筛选时可以先按中心金属分组，再在每组内部按第一近邻元素数量排序。比如统计目标位点周围 3–6 Å 内 Ni、Co、Fe、Mn、Cr 的数量，同时记录表面配位数和最近氧空位距离。这样得到的是一组可解释的局域样本，而不是一堆孤立能量点。
如果后续要做机器学习或描述符拟合，这一步也很重要。中心金属、近邻元素数、空位距离和局域键长可以转成特征；而只给一个“高熵氧化物位点编号”，模型无法知道不同样本之间真正差在哪里。
第一近邻主要调电子还是调几何？
第一近邻元素对吸附的影响通常有两条路。一条是电子路：邻近金属改变中心金属的电荷、磁矩和 d 态位置，让中心位点更容易或更不容易给吸附物提供电子。另一条是几何路：不同半径和价态让表面产生局域拉伸或压缩，改变吸附键长。
这两条路需要分开看。若 Bader 电荷、磁矩和 PDOS 有明显变化，而键长差异不大，说明近邻效应更偏电子调节；若电荷变化不明显，但 M—O、M—ads 键长系统变化，说明几何畸变可能更关键。近邻元素不是背景成分，它可能通过电子和几何两条路改变吸附。
图2 高熵氧化物表面结构与反应响应示例，可用于理解不同局域组成对活性位的调节。DOI：10.1021/acsami.4c22352
第二近邻还值得管吗？
如果体系是金属、合金或高度离域电子结构，第二近邻常常也会影响吸附趋势。高熵氧化物里第二近邻的影响通常弱一些，但在强畸变、氧空位富集或小超胞模型中，它仍可能通过应变和电荷补偿进入结果。
实操上不建议一开始就把第二近邻所有组合都枚举完。可以先算中心金属和第一近邻覆盖较全的一批位点，再检查残差：如果同一中心金属/第一近邻组合仍分散很大，再加入第二近邻、次表层金属或局域畸变指标。
氧空位是活性位还是调节因素？
氧空位可能直接成为吸附位，也可能只是在旁边调节金属位。前者通常表现为吸附物落入空位附近，和低配位金属形成更强相互作用；后者则表现为吸附物仍在金属顶位或桥位，但空位改变了附近金属价态和键长。
dVO−ads = |rVO − rads|
式(1)不是能量公式，而是一个筛选指标：用氧空位到吸附中心的距离，把位点分成空位邻近、空位次邻近和远离空位三类。氧空位是否影响吸附，首先要看它和吸附物是不是处在同一个局域环境里。
图3 氧空位形成能与表面位点描述符示例，说明缺陷位置和局域环境会共同影响缺陷代价。DOI：10.1021/acsomega.6c03286
空位形成能和吸附能谁先算？
若目标是找真实可存在的活性位，应先算空位形成能，再算吸附；若目标是解释某个已知缺陷样品的反应趋势，可以先围绕实验指认的空位附近位点做吸附。两种顺序都可以，但结论口径不同。
比较时至少保留两列能量：Evac 和 Eads/ΔGads。Evac 很高而吸附能很好，只能说明该位点可能很活泼，但缺陷浓度未必高；Evac 低而吸附过强，则可能带来位点占据或表面重构。
如果研究对象是反应中的工作表面，还可以增加补氧或补羟基对照。一个空位在真空模型里很活泼，但在含水或含氧环境中很快被填补，就不应继续按裸空位来解释催化循环。
图4 高熵氧化物表面在缺陷调节下的局域响应示例，提示吸附判断需要结合结构演化。DOI：10.1021/acsami.4c22352
样本数量怎样控制？
第一轮不追求把所有随机结构算完，而是覆盖主要局域类型。可以选择 3–5 种中心金属，每种中心金属下挑 2–4 类第一近邻组合，再对有无邻近氧空位做对照。这样样本数量大约在几十个以内，既能看趋势，也不会让计算量失控。
每个样本都要统一 slab 厚度、固定层数、真空层、k 点密度、磁性设置和吸附物覆盖度。若一个位点优化后吸附物迁移到相邻金属，最终分类要按优化后结构重命名，而不是继续按初始放置位点解释。
SYSTEM = HEO_site_screen_ads_relax
ENCUT = 520
EDIFF = 1E-5
EDIFFG = -0.03
ISPIN = 2
IBRION = 2
NSW = 180
ISIF = 2
LDIPOL = .TRUE.
IDIPOL = 3
# 每个样本记录：中心金属、第一近邻、d_VO-ads、最终吸附位点
最终归因看哪几个指标？
建议把结果分成三层。第一层是几何：最终吸附位点、M—ads 键长、M—O 键长和表面位移。第二层是能量：Evac、Eads、ΔGads 或关键反应台阶。第三层是电子：Bader 电荷、磁矩、PDOS 和差分电荷密度。
如果中心金属改变后能量整体平移，说明中心元素主导；如果同一中心金属内部分散很大，近邻元素或氧空位更关键；如果同一组元素的构型差异仍很大，局域畸变就不能忽略。谁更影响吸附，不靠直觉排序，而靠同一批样本中的分组差异来判断。
图5 多主元氧化物组成与电子性质变化示例，说明多元素引入会改变整体电子结构背景。DOI：10.1038/s41598-023-29477-0
要点一：高熵氧化物筛位点时，中心金属、近邻元素、氧空位和局域畸变要分层记录。
要点二：中心金属决定直接成键对象，近邻元素通过电子和几何两条路改变吸附。
要点三：氧空位既可能是活性位，也可能只是调节因素，需要按距离和形成能分组。
要点四：最终归因应来自分组比较，而不是只看一个最低吸附能构型。
