C-H 活化计算经常从一个看似简单的问题开始：把 CH₄、CH₃OH 或芳香分子放到表面上，选一根 C-H 键断开，再用 NEB 找能垒。实际执行时，困难集中在三处：哪根 C-H 键被活化，H 转移到哪个表面原子，终态是否和初态属于同一个反应事件。C-H 活化过渡态计算的核心，是把被断键、受氢位点和原子映射先固定下来，然后再讨论能垒高低。这个界定会直接影响后面的图像检查和 TS 验证。
吸附构型怎样转成反应端点？
对 CH₄ 来说，初态通常是分子靠近表面酸碱对、金属-氧对或低配位金属位；对苯或取代芳香分子来说，初态应来自已筛过的低能吸附构型。若分子还没有稳定吸附态，直接拉长 C-H 键会把吸附、旋转、平移和断键混在一条路径里。图1给出 γ-Al₂O₃ 上极性路径和自由基路径的 C-H 活化端点，可用于区分 H 转移到 O 位还是生成表面羟基样终态。
端点选择要从几何身份开始。初态中被活化 H 与 C 的距离、H 到受体原子的距离、C 原子是否靠近活性位，都应记录；终态中 C-H 键已断开，H 与表面 O、S 或金属位形成新的键合关系。初态和终态必须对应同一根 C-H 键的连续变化，否则 NEB 能垒会同时包含分子换位或位点迁移贡献。端点几何先闭合，插值路径才有明确方向。
图1 γ-Al₂O₃(100) 表面 CH₄ 极性路径与自由基路径的 C-H 活化结构。该图用于区分受氢位点、过渡态和终态的几何身份。DOI: 10.1021/acsomega.8b02554。
受氢位点怎样枚举？
C-H 活化并不只有一条路径。氧化物表面常比较晶格 O、羟基 O、缺陷邻位 O；金属硫化物表面可能比较桥位 S、端位 S 和金属-S 对；负载金属单原子体系还要比较 M 位直接受 H 与邻近非金属受 H。不同受氢位点会改变 H 转移距离、C 原子配位和电子转移方式。
在 VASP 建模时，可以先为每个吸附构型设置两到三个受氢位点，而不是一次性枚举所有表面原子。筛选依据包括初态 H···受体距离、终态键长是否合理、终态是否为局域极小值、以及反应是否与研究问题相关。受氢位点枚举决定了能垒比较的化学含义：O 位受 H 更像酸碱协同活化，金属位参与则更强调金属-碳或金属-氢中间态。这个分类也会影响电子结构归因。
初态和终态要先独立优化吗？
初态 IS 和终态 FS 都需要先独立优化，并且使用同一 slab、同一吸附物数量、同一电荷/自旋设置和同一固定层。若初态来自 CH₄ 弱吸附，终态来自 *CH₃ + H*，二者总原子数相同，但局域键合方式不同；这正是 NEB 应描述的反应坐标。图2中的能量剖面显示，同一 γ-Al₂O₃ 不同晶面和酸碱位点会给出不同势垒，端点定义会直接进入比较结果。
终态优化尤其要注意 C 片段是否跑到另一类位点。如果 FS 中 *CH₃ 迁移到远端金属位，路径就包含断键和迁移两个过程；此时需要保留局域终态，或者把反应拆成 C-H 断裂和片段迁移两条路径分段处理。IS/FS 的物理身份一致，是 CI-NEB 图像可信的前提。反向插值能进一步核对两端是否仍属于同一事件。
图2 γ-Al₂O₃ 不同晶面和位点上的 CH₄ C-H 活化能量剖面。该图用于比较端点、晶面和酸碱位点对能垒的影响。DOI: 10.1021/acsomega.8b02554。
插值图像怎样检查？
插值后要逐张检查 POSCAR：C-H 距离应逐渐增大，H-受体距离应逐渐缩短，分子主体不应突然翻转，表面骨架不应穿插。对 CH₄ 活化，5 到 7 个中间图像通常可以作为起点；对芳香分子 C-H 活化，分子转动自由度更多，图像数和约束检查要更充分。
IDPP 插值适合减少不合理近距离接触，但它不能替代化学判断。若中间图像出现 H 穿过表面层、C 原子突然改变吸附位、或 H 被两个受体同时过近约束，应重新构造端点或加入中间态。插值质量不是美观问题，而是能否沿同一反应坐标收敛的问题。图像预检的重点是键长连续、位点连续和表面骨架连续。
活化能和反应能怎样分开？
NEB 输出的最高点能量相对初态的差值可作为电子能活化能，终态相对初态的差值对应反应能。若要进入自由能，需要对 IS、TS 近似结构和 FS 做振动校正，并明确低频模式处理口径。电子能剖面先锁定同一反应路径内的动力学差异。
Ea = ETS – EIS，ΔE = EFS – EIS 式(1)
式(1)中 ETS 通常取 CI-NEB 最高图像的能量，后续可用 dimer 或 VTST 优化过渡态并做频率确认。活化能回答跨过势垒需要多少能量，反应能回答终态相对初态是否热力学稳定。一个反应能较负的 C-H 活化路径仍可能有高势垒，一个势垒较低的路径也可能生成不稳定终态。
对比两个表面位点时，Ea 和 ΔE 要并排出现。如果位点 A 的势垒较低但终态较高，位点 B 的势垒较高但终态稳定，二者对应的是不同动力学和热力学组合，反应条件会决定哪条路径更可能占主导。
复算初期只比较电子能剖面，把所有图像相对初态平移到同一零点。路径和过渡态确认之后，再加入振动校正、温度项和外场条件，数据层级就能和结构层级对应起来。
单虚频怎样确认？
过渡态确认不能只看 NEB 剖面最高点。把最高图像提取出来，用 dimer 或有限差分频率计算检查虚频模式，理想情况下应有一个主要虚频，并且振动方向对应 C-H 拉伸与 H 向受体移动。若虚频主要是分子整体转动或表面原子滑移，说明最高点未必是目标 C-H 活化 TS。
图3把 TiO₂ 表面 CH₄ 活化的 IS、TS、FS 和自旋密度放在一起，适合说明过渡态不只是能量最高构型，还应有沿断键坐标的局域电子重排。实际计算中可同时检查 C-H 距离、H-受体距离、C 片段吸附位和虚频振动矢量。TS 验证要同时看能量最高点、几何坐标和虚频方向。这些证据能互相约束。
图3 TiO₂(110) 表面光生空穴辅助 CH₄ 活化的 IS、TS、FS 和自旋密度。该图用于连接过渡态结构与局域电子重排。DOI: 10.1021/jacsau.1c00466。
PDOS 能解释哪一段？
PDOS 可以用来分析过渡态附近 C-H σ 轨道、受体 O 2p 或金属 d 态的耦合。图4把 TiO₂ 体系中过渡态轨道、O 2p/Ti 3d 态和 C-H 活化能垒描述符联系起来，说明电子结构证据应服务于能垒差异，而不是单独代替路径计算。
若某表面的 O 2p 态更接近反应相关轨道，可能更容易接受 H 或形成过渡态相互作用；若金属 d 态和 C 片段耦合增强，终态 *CH₃ 可能更稳定。PDOS 的作用是解释 IS、TS 和 FS 之间的键合变化，它的边界停留在轨道归因，路径可信度仍由端点优化和能垒搜索决定。
图4 TiO₂ 表面 CH₄ 活化过渡态轨道、PDOS 与能垒描述符。该图用于说明电子结构怎样辅助解释 C-H 活化势垒。DOI: 10.1021/jacsau.1c00466。
不同表面位点怎样比较？
图5给出 FeS₂ 表面多个位点上 CH₄ 初始活化的反应物、过渡态和产物结构。读这类图时，应先确认每条路径的受氢位点是否同类，再比较势垒；若一个路径是 Fe-S 对协同活化，另一个路径是 S-S 位点活化，二者的能垒差异同时包含位点酸碱性、几何开放度和电子结构差异。
比较 Fe-S 对、S-S 位点和氧化物酸碱对时，受氢原子和断键坐标先限定同一类反应。同样是 CH₄ 初始活化，一个路径可能是酸碱对协同断键，另一个可能是表面硫位抽氢；即使能垒接近，它们对应的 TS 几何、H 转移方向和 C 片段终态也不同。路径身份变了，讨论对象也随之改变，能垒只能作为各自通道内部的动力学量来读。
图5 FeS₂ 表面不同位点上 CH₄ C-H 活化的反应物、过渡态和产物结构。该图用于比较活性位和受氢机制对能垒的影响。DOI: 10.1073/pnas.2012666118。
输入文件怎样设置？
第一轮先优化清洁表面和分子吸附态，ISMEAR、ENCUT、KPOINTS、固定层和自旋设置保持一致。第二轮优化 *CH₃ + H* 或芳香片段 + H* 的终态，并确认它是局域终态。第三轮用 VTST、nebmake.pl 或 ASE 生成 00、01、02 到 0N 的图像目录；每个目录保留同一原子顺序，POTCAR、KPOINTS、INCAR 和初末端一致。
爬山计算可以先从温和的力收敛标准起步，待路径形成稳定单峰后再收紧。图像数取决于 H 转移距离和分子自由度，CH₄ 通常比芳香分子更容易保持局域路径；芳香分子还要额外检查环平面转动和取代基摆动。
输入文件中需要同步固定图像数、弹簧常数、爬山开关和力收敛阈值。参数只是起点，最终由体系收敛和力标准决定。参数设置要服务于路径收敛，不能弥补端点定义不清。
结果文件怎样复查？
收敛后读取 neb.dat 或各图像 OUTCAR 的能量，画出相对能量曲线；再检查每个 CONTCAR 的 C-H 距离、H-受体距离和表面键长变化。若最高图像两侧能量不平滑，先看是否存在图像跳位；若最高点不在预期 C-H 断裂附近，重新检查受氢位点和插值路径。
结论部分只放和反应事件直接相关的量：IS/FS 结构、TS 结构、Ea、ΔE、关键键长、虚频模式和电子结构解释。C-H 活化结果可靠与否，取决于端点、路径、能垒和 TS 验证是否指向同一个化学事件。当这四类证据一致时，电场、光生空穴、缺陷或溶剂环境对势垒的调控才有清晰边界。
要点：C-H 活化要先确定被断 C-H 键和受氢位点，再构建 IS/FS 与 NEB 图像。
要点：活化能和反应能分别回答动力学跨越与终态稳定性，二者不能互相替代。
要点：TS 确认需要能量最高点、关键键长和单虚频方向共同支持。
要点：PDOS、自旋密度和电荷分析用于解释势垒差异，应回到同一条反应路径中使用。
