DFT 如何研究催化剂表面重构？

说明：本文华算科技主要介绍DFT如何定义催化剂表面重构、生成候选原子结构、比较反应条件下的表面自由能，并用能垒与有限温度模拟判断重构能否在实验时间内发生。

催化表面模型通常从体相晶体切出薄层，再加入真空层并周期复制。刚切出的体相终止面只是初始坐标；原子弛豫后，表层键长和层间距可以改变。若二维平移周期、表面原子数和成键拓扑保持原样，这类变化通常称为表面弛豫。

表面重构涉及二维周期、表面化学计量、堆垛次序或成键网络的改变，例如从（1×1）变为（2×1）、台阶原子增多、金属原子迁出晶格位点，或吸附物形成新的有序覆盖层。

DFT实际处理的是给定晶胞、原子数、电子数和外部约束下的一组局部能量极小构型。建模时需写明晶面指数、表面超胞、薄层厚度、固定层数、真空层以及是否采用对称薄层。相同晶面可以拥有多种表面周期和表层组成，从单一理想终止面启动的优化局限在该初始构型附近的势能盆。

RuO₂(110)的四组候选结构具有相同基底，重排集中在表面约3–5 Å的区域。RuO₄–(2×1)与RuO₂–(1×1)位于热力学稳定包络线上，Ru₄O₉–(1×1)和Ru₈O₁₇–(1×2)接近凸包。这里的化学计量和表面周期是结构身份的一部分，仅标注“RuO₂(110)”不足以区分这四类模型。

清洁表面的自发重排、气氛诱导的缺位或增原子结构、吸附层有序化、电极溶出与再沉积、表面氧化物或羟基层生成，都可进入重构研究。每种情形使用的粒子库不同。吸附物身份、覆盖度、表面电荷和溶剂环境必须写入模型定义，真空中的中性裸表面能量无法代表所有工作状态。

某个重构态能量较低，表示选定参照和约束下存在热力学驱动力。原子从初始表面迁移到该构型仍需跨越势垒，反应时间内能否完成取决于扩散路径、温度和吸附物交换速率。终态排序与重构速率属于两类计算任务，后者需要过渡态或有限温度动力学信息。

候选集合可由不同表面超胞、终止层、表层原子数、空位位置、台阶密度、吸附位点和覆盖度组合而成。已知实验周期时，可以人工枚举少量结构；未知周期或成分时，组合数会迅速增长。一次几何优化只沿初始构型所在的势能盆下降，不会自动搜索另一种周期、另一组原子数或远离初态的成键方式。

进化算法、遗传算法、盆地跳跃、极小值跳跃和随机结构搜索通过变异、交叉、增删原子或改变吸附层排布产生新候选。低成本势函数或机器学习势可承担早期采样，低能结构再交给DFT复算。搜索数据库还需按对称性、原子邻接或结构指纹去重，避免大量重复弛豫。结构生成规则决定了搜索能够访问的构型空间，未允许改变原子数的搜索无法发现变组成重构。

计算收敛应观察低能包络是否在多轮搜索后保持稳定，并用不同随机种子重复搜索。能量接近的构型应作为低能集合保留，因为有限温度下可能同时占据，且泛函、溶剂和振动自由能的修正足以交换其排序。最低能候选附近的结构密度能反映是否存在大量亚稳态，单独报告一个最优坐标会丢失这部分信息。

Cu表面的研究在每个晶面上采样约千个CO覆盖结构，允许CO与给定化学势的粒子库交换，再以恒电位DFT重算全局最低能态和低能亚稳集合。PBE、RPBE和随机相位近似对CO位点偏好给出不同结果，研究者据此修正全局搜索中的位点能量。泛函误差会改变吸附位点与平衡覆盖度，并继续传递到表面稳定性排序。

薄层化学计量改变后，DFT总能包含不同数量原子的能量贡献，比较时需引入化学势。常用量可写成γ=(G_slab−ΣN_iμ_i)/A；双面等价薄层还需按两个表面归一。面积A、原子数N_i和化学势μ_i须采用同一参照；截断能、k点密度、薄层厚度、磁态、色散修正与溶剂模型也应保持可比。

反应气氛中的表面会与气相、电解液和电子库交换粒子。每个候选结构的表面自由能或表面巨势可概括为Ω_i=[G_i−ΣN_jμ_j−qU]/A，具体形式随电荷控制方式和参照电极而变。稳定表面由给定T、p、U和pH下最低的Ω_i确定，而非真空DFT总能最低的结构。

气相分子的μ(T,p)包含温度和分压项，因而氧化性、还原性气氛会改变缺氧或富氧终止面的斜率。计算氢电极近似常写为μ(H⁺+e⁻)=½G(H₂)−eU_SHE−k_BTln10·pH。电位标尺、pH、气体分压和温度应随相图同时给出；改用RHE标尺时，pH项会进入电位换算。

裸Cu表面在零电荷电位以及−1.5 V vs SHE下仍以致密的Cu(111)和Cu(100)较稳定。加入各晶面的平衡CO覆盖后，台阶与扭折表面因CO结合更强而降低自由能，稳定顺序发生反转。吸附覆盖带来的自由能收益可补偿缺陷表面的形成代价，这给出平整表面生成台阶和扭折位点的热力学驱动力。

H和CO覆盖会分别改变Cu(111)、Cu(100)与Cu(110)的巨势界面能，Wulff构形据各晶面能重新分配外露面积。原研究给出的电位标尺误差约为0.3 V，说明相界位置受模型近似影响。泛函、隐式溶剂、双电层电容和振动熵都会移动相界；接近相界的结构排序应附能量差和误差来源。

相图中的每个区域表示平衡假设下的最低自由能构型。若候选库遗漏溶剂共吸附、表面溶出物、不同超胞或非化学计量结构，得到的区域只能在现有候选之间成立。相图范围等于候选结构、粒子库与自由能修正共同限定的范围，计算记录中应保留所有近稳结构及其参照。

自由能排序给出起态与终态，原子迁移过程还需单独求解。路径明确时，可用爬山弹性带或dimer方法定位鞍点；表面原子交换、空位扩散和吸附物辅助迁移往往含多个中间态。活化自由能ΔG^‡进入速率常数k≈(k_BT/h)exp(−ΔG^‡/k_BT)，相差0.2 eV的势垒在室温下可造成数个数量级的速率差。

AIMD按飞秒时间步推进原子轨迹，常见计算长度为皮秒到数十皮秒，适合观察快速键交换、溶剂响应和高温重排。高势垒事件在短轨迹中可能一次也不出现。增强采样可沿选定集体变量跨越势垒，动力学蒙特卡洛则把已知事件及速率推进到秒级以上。采样时长必须与待观察事件的时间尺度匹配，完整的短程轨迹只覆盖所模拟的皮秒区间，秒级迁移需由势垒或长时间模型检验。

Pd_n/CeO₂(111)的多尺度模型用DFT提供团簇能量、CO吸附能与扩散势垒，再训练簇展开哈密顿量。遗传算法搜索不同尺寸的团簇结构，动力学蒙特卡洛处理原子迁移、聚并与CO吸附事件。DFT在这里负责标定局部能量和事件势垒，长时间结构演化由统计动力学方法承担。

从分散Pd原子出发，CO覆盖和团簇层数比例在早期迅速变化，随后进入缓慢演化区间；实线与阴影分别记录三次模拟的均值和95%置信区间。时间序列可区分快速局部重排与动力学冻结的亚稳集合，而最低能团簇本身没有这项时间信息。

结构与实验信号的比较应从同一工作条件生成可观测量：候选原子模型可用于模拟STM、振动频率、X射线散射或吸附能分布，动力学集合则按占比平均。候选结构收敛、条件相关自由能、转变势垒和时间分布分别回答“有哪些构型”“哪种受热力学偏爱”“能否跨越路径”和“实验时段内保留哪些结构”。四类结果的温度、压力、电位、覆盖度与模型组成需保持一致。