表面羟基化怎么建模：终止方式、OH覆盖、质子排布与水的化学势约束

说明：本文华算科技主要介绍表面羟基化模型中的原子组成、表面终止、羟基覆盖与质子排布，以及水的化学势、溶液酸碱条件和实验信号对结构选择的约束。

表面羟基化指固体表面形成与基底成键的 O–H 基团。在许多氧化物上，水的氧与配位不饱和金属位点结合，一根 O–H 键断开，质子转移给邻近晶格氧。弛豫后的表面含有水来源的羟基和质子化的晶格氧，二者都带 O–H 键，却有不同的氧原子来源和金属配位。

若清洁基底记为 S，添加一分子中性水后，整个模型只新增一个 O 和两个 H。水在金属位与晶格氧之间解离时，可形成两个表面羟基；若直接在两个位置各添加一个 OH，则新增两个 O，已经改变了氧计量。ZnO 的裸表面与解离水结构提供了这种原子计数的实例。

羟基的端位、双配位与三配位通常按其氧原子连接的金属数区分。水来源的 OH 也可能同时连接两个金属，质子化晶格氧则可能保留多个金属邻居。图中 ZnO 的水来源 OH 处于双金属配位环境，另一个 H 与原有表面氧成键；用统一竖直的 OH 替代两种几何会改变局部氢键。

仅给原有表面氧添加 H，描述的是另一种组成变化。中性 H 原子、溶液质子和质子加电子具有不同电子计数；中性 OH 自由基与 OH⁻也有电子数差别。对中性水解离，晶胞总电荷可保持不变，局域电荷仍会重排。计算中的总电子数、净电荷和可能的自旋态，应与所写反应的原子和电荷守恒一致。

表面上还可保留完整配位水与未质子化晶格氧。羟基化结构允许这些物种共存；把每个暴露金属都设成金属–OH，会预先排除分子吸附和部分解离。对有氧空位的表面，水的氧还可能填入空位，模型需记录填位前后的缺陷数及氢所在位置。

晶面指数只规定截面朝向，最外层的元素组成与重构决定实际暴露的位点。极性氧化物的任意体相截断可能产生人工电场，需比较满足计量和电荷补偿要求的终止。单面羟基化薄片具有法向偶极时，应检查真空静电势与偶极修正；增加薄片厚度和横向尺寸后，目标吸附能差应趋于收敛。

α-Al₂O₃(0001) 的一类饱和模型采用两侧羟基化结构，表面 OH 的氧与两个 Al 配位。同一套计算还比较增加 H 或 OH 的反应自由能：在 pH = 0、相对标准氢电极的 0–1.23 V 区间，所选单侧晶胞的 12 个羟基构型保持饱和，再增加第 13 个 H 或 OH 均不利。这个数目对应指定晶胞与终止。

跨材料比较可报告羟基面密度 ρ_OH = N_OH/A，单位为 OH·nm⁻²，其中 A 与羟基计数取同一侧表面。覆盖度 θ = N/N_site 还需指定参考位点；在水解离成对产物中，总羟基数可为解离水数的两倍。端位 OH、质子化晶格氧和配位水分别计数，可避免不同分母造成的“满覆盖”歧义。

含晶界的氧化铝表面出现连接一个、两个或三个 Al 的羟基，局部酸碱反应能随配位变化。同样数量的氢可以占据不同氧原子，O–H 取向又改变氢键供受体排列。模型集合应包含相邻与分离的质子、不同配位位点，以及羟基和配位水的共存构型。

两种晶界结构的侧视和俯视图保留了不同的金属–氧骨架，配位标号可定位加氢反应发生在哪里。对每一种骨架，在相同组成和净电荷下比较多种质子排列，才是在比较构型能量。若小晶胞强制所有 OH 同向、等间距排列，所得能量还包含人为有序及周期副本相互作用。

几何优化会在局部势能面内调整键长，却可能保留初始质子排列。对能量接近的候选结构，可改变氢键朝向、放大晶胞复算，并检查受约束原子是否妨碍表面弛豫。金属–氧重构、质子转移和水脱附若导致终态物种变化，应按优化后的组成重新分类，不能仍沿用初始“羟基覆盖率”。

不同水分子数的表面交换了不同数量的原子，电子总能不能直接排序。若所有候选结构具有同一基底，只与水库交换中性 H₂O，可用Δγ_n = [G_n − G₀ − nμ_w]/A_act比较。n 是新增水分子数，μ_w为水的化学势，A_act是参与交换的表面积；两侧等价且一起羟基化时取 2A，固定背面的单侧差分取 A。

气相水的化学势含温度、分压与熵贡献，理想气体近似为 μ(T,p) = μ°(T) + k_BT ln(p/p°)。表面自由能可由 DFT 电子能、零点能和热修正估计。ZnO 的相图比较裸露、吸氢和吸水结构，水分支采用水的化学势；氢分支对应氢气库，不能将两类吸附物的化学势当成同一个物理量。

对于图中的 ZnO 结构 A，一个解离水分子的构型在所给水化学势约 −1.85 至 −1.15 eV区间占优，水化学势升高后，三分子水的候选结构能量更低。新增水数改变曲线斜率，交点给出这些候选态之间的平衡转换。若终止之间还交换金属或额外氧，应增加相应化学势项，并限制在基底相允许存在的条件内。

平衡相图没有记录水解离要经过的中间构型。TiO₂ 的动态模拟区分直接向晶格氧供质子和经邻近水分子的两步转移：金红石(001)与锐钛矿(101)呈现不同偏好。邻近水参与反应时，两种羟基产物可相隔更远，静态建模只保留最近邻羟基对会漏掉这类构型。

氢环境投影图和状态转换统计来自DFT 数据训练的机器学习势模拟，记录了质子在不同环境间的变化。对具体材料，解离与复合的自由能差决定平衡比例，势垒影响到达平衡的时间。短时间轨迹若始终停在初始状态，只能报告该时间内的驻留；覆盖比例需用足够的解离、复合事件或自由能采样确定。

若多个羟基排布的自由能差接近 k_BT，有限温度下可能共同占据，只保留最低自由能构型会忽略不同排列的构型熵。模型可通过构型统计计算这些排布的权重。浸液表面采用液态水化学势，不能将理想水蒸气公式无限外推到凝结后的高压区；溶液酸碱变化则另由质子交换描述。

在液态水中，表面羟基与配位水持续交换质子。显式水层可描述氢键重排与水参与的质子转移，连续介质主要近似溶剂的平均极化，难以单独给出这些局部反应。统计羟基时应根据 O–H 配位随时间更新物种身份，并检查水层厚度、采样时长与不同初始排列对平均组成的影响。

锐钛矿 TiO₂(101) 的中性溶液模拟中，两类解离产物平均占据比例接近；碱性条件增加金属位吸附 OH⁻，酸性条件下部分端位 OH 接受质子形成配位水。图中两类曲线分别计数吸附质子和吸附羟基，酸碱条件使二者出现数量差，表面电荷也随之改变。

把几个 H 从表面删除并指定“高 pH”，尚未定义溶液的质子库。pH 通过质子化学势与表面酸碱平衡影响组成；显式模型还需设置溶液离子、补偿电荷及其位置。电极体系涉及电子交换时，需给定电位标尺并采用一致的电子化学势。带净电荷薄片的均匀补偿背景会产生有限尺寸误差，其作用与溶液反离子的局域屏蔽不同。

实验约束应对应模型的晶面、预处理和水暴露状态。ZnO 在法向与掠射出射条件下的 O 1s 谱具有不同表面权重；室温暴露于 D₂O 后，高结合能侧出现新增组分。插图记录残余气体暴露过程中峰面积随时间增加，提示转移和等待时间也能改变测量时的表面组成。

该 ZnO 表面在额外通水之前已有解离水形成的羟基，部分信号与其他表面氧重叠。因而新增高结合能峰的面积主要记录指定暴露过程的光谱变化，不能直接换算成全部羟基数量。候选结构的位点分辨芯能级位移、实验出射方向和吸附水贡献，应共同用于模拟 O 1s 谱。

对可获得振动光谱的样品，孤立羟基、氢键羟基和配位水具有不同振动环境；氢键通常使 O–H 伸缩频率降低并产生分布。H/D 替换后对应振动移向低频，可检验含氢物种的峰归属。模型对照应保留实际含水量和羟基邻居，比较计算频率、同位素位移及峰随脱水处理的变化。