二维材料加应变后，如何分析催化活性？

引言

二维材料的应变调控看起来很直观：把晶格拉长或压缩，活性位的键长、轨道重叠和吸附自由能都会改变。但在 VASP 结果里，晶格拉伸、PDOS 位移和 ΔG_H* 并不会自动形成一条简单直线。应变是结构变量，催化活性是吸附态、电子结构和反应路径共同给出的结果。

应变模型先限定哪类自由度？

二维材料施加应变时，第一步不是直接放 H，而是确定应变定义。常用双轴应变可写为 ε = (a – a₀)/a₀ × 100%，其中 a₀ 是零应变优化后的面内晶格常数，a 是施加应变后的面内晶格常数。对单层 MoS₂、石墨烯负载 MoS₂ 或其他二维催化剂，应变通常只作用在面内方向，真空层高度保持足够大。

VASP 实操中，可以先完成零应变结构优化，再按 −5% 到 +5% 等范围缩放面内晶格，随后固定面内晶胞并优化离子位置。对二维 slab，常用 ISIF = 2 或等价设置保持面内应变；如果释放晶胞，体系可能回到低应变状态，最后得到的就不是指定应变下的吸附能。应变计算的可比性来自相同晶格约束，而不是每个模型都追求最低总能。

图1 MoS₂/graphene 不同缺陷位点在双轴应变下的 H 吸附构型。图中同时展示了顶视图、侧视图和 H 吸附位置。DOI: 10.1039/d2ra07363c。

图1 给出了多个缺陷 MoS₂/graphene 位点在应变下的 H 吸附构型。这个图适合说明一个关键边界：应变不是孤立改变一个键长，而是同时改变缺陷邻近原子、载体耦合和吸附位点几何。若 H 在拉伸后换到另一个位点，ΔG_H* 的变化就不再只代表原位点键长变化。

在应变序列里，先看同一位点是否还能被追踪：同一 V_S2 位点在 −3%、0%、+3% 下都保持 H 吸附在相同 Mo 邻近位点，才适合讨论应变趋势；如果 +3% 下 H 转移到边缘或缺陷孔另一侧，应先把它列为构型切换，再讨论活性变化。相同应变序列里，吸附位点身份必须先保持可追踪。

ΔGH* 随应变为什么可能非单调？

HER 中常用的氢吸附自由能可写为：

式(1) ΔG_H* = E_slab+H – E_slab – 1/2E_H₂ + ΔZPE – TΔS

这个式子说明 ΔG_H* 同时包含吸附态总能、裸表面应变能、H₂ 参照和热力学校正。应变改变的不只是 E_slab+H，也会改变裸表面局域结构，因此比较时必须让每个应变下的 E_slab 与对应吸附态成对出现。

图2 展示了不同缺陷 MoS₂/graphene 在双轴应变下的 ΔG_H* 变化。可以看到，不同缺陷的曲线方向并不相同：有的拉伸让 H 吸附更接近热中性，有的则使吸附偏强或偏弱。应变对 ΔG_H* 的影响具有位点依赖性，不能把拉伸或压缩直接等同于活性提升。

图2 不同缺陷 MoS₂/graphene 催化剂的 ΔG_H* 与双轴应变关系。图中展示了缺陷类型和应变方向共同影响氢吸附自由能。DOI: 10.1039/d2ra07363c。

在 VASP 任务安排上，建议每个应变点都包含三类计算：应变后裸二维材料优化、同一位点 H 吸附优化、吸附态频率或热力学校正估算。若只在零应变下做 ZPE/熵，而其他应变点直接套同一修正，需要说明这是近似处理；若 H 振动频率随应变变化明显，则应对关键应变点重新计算频率。

换到缺陷横向比较时，V_S、V_S2、V_MoS3 的 ΔG_H* 可能在同一应变下接近，但这并不代表它们的电子结构相同。若一个位点来自缺硫造成的 Mo 暴露，另一个来自 Mo-S 空位复合缺陷，它们的局域配位、载体耦合和磁性都可能不同。因此，跨缺陷比较应把构型图、ΔG_H* 和 PDOS 放在同一组证据里。

PDOS位移能否直接解释吸附变化？

PDOS 常被用来解释应变调控，因为拉伸或压缩会改变金属 d 态和吸附物轨道的重叠。图3 给出缺陷 MoS₂/graphene 中 Mo 4d 态以及应变下 d 态中心和 H 化学键合能的关系。它说明电子结构描述符确实能提供线索，但线索要回到具体位点和吸附态中验证。

图3 缺陷 MoS₂/graphene 的 Mo 4d PDOS、d 态中心和 H 化学键合能关系。图中比较了缺陷类型及双轴应变对电子态位置的影响。DOI: 10.1039/d2ra07363c。

PDOS 位移和 ΔG_H* 之间不能简单写成一一对应。第一，PDOS 峰位变化可能来自局域 Mo 配位改变，也可能来自载体耦合改变；第二，H 吸附后还会引入 H 1s 与 Mo d 态杂化；第三，结构重构能可能抵消一部分电子结构收益。PDOS 能说明轨道能级和杂化趋势，不能单独替代吸附自由能计算。

更可靠的分析顺序是先确认几何位点，再比较 ΔG_H*，最后用 PDOS、差分电荷密度或 Bader 电荷解释趋势来源。若 PDOS 显示 d 态中心上移，但 ΔG_H* 反而远离热中性，应检查是否出现 H-M 键过强、吸附位点切换、局域畸变能增加或载体界面电子转移改变。

图4 的自由能台阶进一步说明，同一类缺陷在不同应变下可能从吸附偏弱转向吸附偏强，也可能在热中性附近出现最优区间。对教程类计算而言，最有价值的不是挑一个最低点，而是给出一套可以复现的判断：相同位点、相同应变约束、相同参照态和相同后处理口径。

图4 缺陷 MoS₂/graphene 在不同双轴应变下的 HER 自由能台阶。图中给出多个缺陷模型中 ΔG_H* 对应变的响应差异。DOI: 10.1039/d2ra07363c。

VASP结果怎样形成闭环判断？

闭环判断至少包含四个层次。第一是结构层：零应变和各应变点的晶格常数、真空层、固定层或面内约束要一致。第二是能量层：每个应变下都用对应的裸表面与吸附态相减，不把零应变裸表面拿来做所有应变点参照。第三是电子结构层：PDOS 静态计算使用相同投影方式、相同能量零点和相同展宽。第四是反应层：当 ΔG_H* 接近热中性时，再看 Volmer、Tafel 或 Heyrovsky 势垒是否同步改善。

应变能本身也应记录，可写为：
式(2) E_strain(ε) = E_slab(ε) – E_slab(0)。
如果某个应变点虽然给出较好的 ΔG_H*，但 E_strain 很高，或者结构优化后出现明显断键、翘曲和虚频，该点更适合作为机理边界，而不是直接作为可实现材料状态。

图5 给出了 V_S2-MoS₂/graphene 上 Volmer、Tafel 和 Heyrovsky 反应的最小能量路径。这个图把吸附自由能进一步连接到动力学势垒：ΔG_H* 接近零只能说明 H* 热力学吸附较合适，不能自动说明水解离、H-H 偶联或电化学脱附也快。活性判断的最后一步应把自由能描述符和关键基元步骤能垒分开验证。

图5 V_S2-MoS₂/graphene 上 Volmer、Tafel 和 Heyrovsky 反应的最小能量路径。图中给出初态、过渡态、终态和对应势垒。DOI: 10.1039/d2ra07363c。

实际交付一组应变计算时，可以把结果组织成三张图：ΔG_H*-应变曲线、代表性位点的 PDOS 对比、关键路径势垒图。若三者方向一致，可以说应变通过局域电子态调节 H 吸附并影响后续步骤；若三者不一致，就应把差异说成结构位点切换、电子态响应或动力学步骤分离，而不是强行用单一描述符解释全部现象。

还要给应变范围设置物理边界。二维材料在小应变范围内可能保持晶格连续，超过一定范围后会出现局域断键、褶皱、层间滑移或载体脱耦。VASP 结果中若伴随异常大的残余力、明显结构重排或声子虚频，应把该应变点从常规趋势线中单独拿出讨论。可实现的应变区间应同时满足结构稳定、位点身份清楚和能量比较口径一致。

本章要点总结

要点：二维材料应变要先定义面内晶格约束，不能让优化过程自动释放目标应变。

要点：ΔG_H* 需要用同一应变下的裸表面和吸附态成对相减。

要点：PDOS 位移能解释电子结构趋势，但不能单独替代吸附自由能和构型检查。

要点：当 ΔG_H* 接近热中性时，还应继续检查关键基元步骤势垒和结构稳定性。