表面限域与界面限域：原子坐标、扩散方向及层间配位的对比

说明：本文华算科技主要介绍表面限域和界面限域的原子坐标特征、扩散方向、层间配位，以及区分几何约束与化学作用的计算量。

表面/界面限域由哪些坐标定义？

表面限域保留一个开放反应半空间。吸附原子或中间体位于固体外表面，法向坐标受吸附势阱约束，面内坐标仍可沿台阶、条纹、空位列或晶格沟槽迁移。限域强弱不由模型盒子的真空层厚度命名，而由势能面上可达位点、相邻位点间能垒和脱附代价共同刻画。条纹表面可把二维随机迁移压缩为近似一维的优选方向，即使吸附物上方仍是气相或真空。

将三维概率密度 ρ(x,y,z) 投影到表面法向，可得到吸附态与夹层态的 z 坐标分布。开放表面上的吸附态在法向 z 方向形成靠近表面的概率峰，远离表面后进入气相区域；面内 x、y 方向仍存在周期位点。夹层客体的 z 分布位于两片材料之间，并随上、下侧键长出现一个或多个局域峰。概率峰占据的坐标范围可由 AIMD 轨迹统计，低温静态结构则给出其中一个能量极小点。

MnO(001)/Pt(111) 的表面形成周期条纹，Mn₃O₄(001)/Pt(111) 则呈网格状起伏。两者都属于开放表面模型，差别写在面内 O 位点排列和各向异性的局域几何里。表面限域在这类体系中表现为横跨条纹与沿条纹方向具有不同的迁移代价，并非把氢封闭在孔腔中。

界面限域增加第二个固体接触面。客体原子位于两层材料之间，法向位移会改变上、下两侧的键长；面内迁移还会改变与两层晶格的配准。单层表面上的吸附位和多层夹缝中的插层位具有不同的配位数、形成能与可达构型。计算模型应保留两侧材料、真实层间距和客体原子的相对高度。

Zn 单原子位于 1T‑MoS₂ 基面上方时，HAADF‑STEM 强度可区分替位模型与吸附模型；放入层间后，上、下两片 S 基面共同参与配位。双侧配位是界面限域的结构信号。只有第二层使能量、键长或构型分布发生可分辨变化时，双侧作用才进入限域定义；第二层远离客体后，上侧相互作用趋近于零，模型回到单侧吸附极限。

表面限域为何产生方向依赖扩散？

表面扩散可写成位点序列上的热激活跃迁。沿 x、y 两个方向的最低迁移能垒分别为 E_m,x 与 E_m,y 时，扩散系数 D_x、D_y 由各自的跃迁频率控制。E_m,x 与 E_m,y 的差值比单个吸附能更接近“方向受限”的物理含义；吸附能描述某个位点的稳态深度，迁移能垒描述相邻稳态之间的鞍点高度。

过渡态理论把方向 i 上的跃迁频率写成 k_i≈ν_iexp[−E_m,i/(k_BT)]，扩散系数还含跳跃距离与可用相邻位点数。两个方向相差 0.2–0.4 eV 的能垒在室温下会产生很大的速率比。能垒差控制指数项，振动频率 ν_i、空位密度和覆盖度则进入前因子与可用路径数。

氢在条纹 MnO 单层上由 Pt 裸露区进入氧化物岛，羟基化亮线沿 [01̅1] 方向向岛内延伸；网格状 Mn₃O₄ 上的亮区由岛边缘近似各向同性地推进。时间分辨 STM 提供迁移前沿，O 1s 谱中晶格氧与 OD 峰面积记录羟基化进程。方向与速率是两个独立观测量：条纹选择迁移方向，位点配位和 O–H 稳定程度调节绝对速率。

NEB 能量曲线把这种差别转成可计算量。MnO/Pt(111) 沿 [01̅1] 的两步能垒约为 1.01 eV，横跨条纹的步骤升至 1.47 eV 和 1.21 eV；Mn₃O₄/Pt(111) 所示路线为 1.19 eV 与 1.23 eV。最低连续能垒序列决定优选迁移方向，孤立地挑出某一步最低值会漏掉后续高鞍点。

图中的 O‑p 带中心由 −3.12 eV 移至 −2.61 eV，并伴随 H^* 吸附能由 −0.39 eV 变为 −1.13 eV。后者对应更深的 O–H 稳态，迁移时付出的能量随之增大。这里的表面限域含有晶格各向异性与局域成键强度两部分；温度、羟基覆盖度、氧空位和表面重构会改变实际跃迁频率，静态 NEB 给出的是选定构型下的势能路线。

界面限域为何改变配位与层间距？

客体原子进入层间后，上、下两片晶格会弛豫到新的平衡距离。层间距过小时，短程排斥抬高插层能；距离过大时，客体与其中一侧的键合减弱，双侧约束逐渐退化为单侧吸附。层间距是界面限域的显式结构变量，其数值应和插层浓度、层数、堆垛次序及色散修正共同报告。

Zn 插层使 1T‑MoS₂ 的 (001) 衍射峰由 14.07°移至 13.59°，对应层间距由 0.629 nm 增至 0.651 nm，增幅约 3.4%。HRTEM 晶格条纹给出局域距离分布。峰移记录长程平均周期，显微图记录纳米区域起伏；平均层间距与局域配位构型承担不同信息，二者数值无需逐点相等。

层间距变化只回答几何响应，Zn 的化学身份和近邻原子还由 XAS 限定。Zn K 边位置接近 ZnO 参照，对应以 Zn²⁺ 为主的状态；R 空间 EXAFS 在 Zn–S 距离附近出现主峰，未出现 Zn 箔中明显的 Zn–Zn 壳层。Zn–S 配位而非 Zn–Zn 聚集支持单原子处于硫基面包围的局域环境。

DFT 优化中，单层 1T‑MoS₂ 上的 Zn–S 距离超过约 2.95 Å；含上下基面的模型可形成四面体型 Zn–S 配位。几何夹持与成键配位在同一次弛豫中耦合，固定层间距的扫描可得到纯几何趋势，全弛豫模型则包含层片形变和键能变化。两类结果应采用相同泛函、色散方案、超胞尺寸和 Zn 化学势。

插层形成能可写成 E_ins=E_host+Zn−E_host−μ_Zn。若 host 与 host+Zn 采用不同层数、堆垛或应变，E_ins 会混入基体差值。固定层间距 d 扫描 E(d) 可定位排斥区、平衡距离和解离区；每个 d 上再次优化 Zn 的面内位置，可识别配准变化。化学势与结构参照必须一致，否则能量排序对应不同插层反应

如何区分几何约束与化学作用信号？

表面模型可设置沿条纹、横跨条纹和网格表面三组等覆盖度路线，分别比较 E_m、跳跃距离与尝试频率；界面模型可设置无 Zn、单侧 Zn 吸附、双层 Zn 插入以及多种层间距，比较形成能、Zn–S 配位数和 ΔG_H*。只改变一个结构变量的成对模型可把晶格方向、第二基面和客体成键的贡献分开记录。

静态电荷密度与 PDOS 取自单个弛豫构型，NEB 描述预设初态与终态之间的最小能量路线，AIMD 则记录有限温度下已经发生的位点交换；三类结果对应不同时间尺度。表面重构或层间滑移若慢于 AIMD 时长，轨迹会停留在一个亚稳构型。多起点轨迹与独立 NEB 路线可检验势能面上是否存在未采样的迁移方向。

Zn/1T‑MoS₂ 的三层模型把 H^* 放在基面 S 位点。该论文给出的 ΔG_H* 从未插层模型的 −5.40 eV 变为 Zn 插层模型的 0.00294 eV，数值对应 0 V（vs RHE）和 pH=0 的指定计算设置。自由能变化描述 H^* 吸附—脱附热力学；HER 电流还受质子迁移、界面电阻、可接近位点数量和电位下表面状态控制。

几何约束优先反映在方向分辨能垒、层间距、原子高度与取向分布；化学作用优先反映在配位数、键长、电荷分区、PDOS 和吸附自由能。AIMD 还能给出 D_x/D_y、驻留时间和配位数随时间的分布。若缩短层间距后扩散方向改变而电荷与成键特征近似保持，几何贡献占主导；若固定几何后替换壁面元素仍使 ΔG_H*、PDOS 或形成能显著变化，差异来自客体—基面相互作用。

有限温度轨迹还能生成自由能 F(q)=−k_BTlnP(q)+C。表面体系的 q 可取 x/y 位移、O–H 配位数或岛边缘距离，界面体系的 q 可取 z 坐标、上/下侧 Zn–S 键长或层间距。自由能盆地的数量与盆地间势垒分别对应可驻留构型和构型互变代价；采样时长应覆盖多次跃迁，单次短轨迹只记录其访问过的局部区域。

表面限域的计算对象是开放表面上的位点网络，界面限域的计算对象是两侧材料共同限定的局域构型。前者以 E_m,x/E_m,y 和面内轨迹检验方向选择，后者以层间距、上/下侧键长、Zn–S 配位数和 ΔG_H* 检验双侧作用；温度、覆盖度、电荷态和化学势应与所比较的原子模型保持一致。