说明：本文华算科技主要介绍同一个分子在 DFT 吸附计算中为什么要测试多个吸附构型，以及初始位点、取向、局域极小值、构型初筛和全优化结果怎样影响吸附能与电子结构判断。

吸附构型是什么？

在表面吸附计算里，一个分子放到 slab 上并不是只有一种“摆法”。同一个 CO2、H2O、O2、1-hexene 或反应中间体，可以靠近 top 位、bridge 位、hollow 位、缺陷位、掺杂位和边缘位，也可以用不同分子轴、不同倾角、不同吸附原子朝向表面。计算输入文件里的初始坐标，只是势能面上的一个起点；结构优化会沿着受力方向下滑到附近的局域极小值。

这里的“多个构型”至少包含三类变量。第一类是 吸附位点，也就是分子靠近哪一个表面原子或哪一组表面原子；第二类是 分子取向，例如 C=C 键平行或倾斜于 Fe(110) 表面，羰基 O 原子朝向金属还是背离金属；第三类是 初始距离和覆盖度，分子离表面太远会接近物理吸附，离得过近又可能引入不真实的排斥或结构畸变。

真实吸附态由总能、原子受力和电子结构共同约束，不能只靠初始摆放直觉来判断。若只把分子放在一个看起来顺眼的位置，优化后的 Eads 可能只是某个局域极小值，而不是该模型下更低能的吸附态。对于金属表面、氧化物缺陷、单原子位点和掺杂碳材料，同一分子在不同构型里的成键原子、键长、吸附高度和电荷重排区域都可能改变，表面局域配位环境差别会进入能量排序。

图1. 多构型吸附计算流程示意。

一个完整的吸附构型测试流程，通常会把分子在非等价位点上平移、旋转，再做局部弛豫或初筛。流程中 构型生成不是装饰步骤，它位于能量排序、最低能结构筛选和电子结构分析之前。若起点集合太窄，后面的 PDOS、差分电荷密度、Bader/Hirshfeld 电荷和反应路径都可能围绕一个偏高能构型展开。

多个构型测试的核心作用，是避免把单一初始坐标误当成分子在该表面上的代表性吸附态。吸附能常写作 Eads=Eslab+molecule-Eslab-Emolecule。这个式子看起来只比较三个能量，实际要求三套结构状态可对应：slab 要与吸附模型里的表面相同，孤立分子要用同一泛函和自旋设置，吸附体系要完成足够充分的弛豫。

吸附位点如何影响分子相互作用？

表面不是平整的数学平面。金属晶面上有不同配位数的原子，氧化物表面有阳离子、桥氧、端位氧和氧空位，碳材料上可能有 N、B、S、P 掺杂位或 Stone-Wales 缺陷。分子靠近哪个局域环境，会改变前线轨道与表面态的重叠方式，也会改变短程 Pauli 排斥、范德华作用和电荷重新分布。

以 Fe(110) 表面的高对称位点为例，top、bridge 和 hollow 位对应不同数量的近邻 Fe 原子。一个含 C=C 键的 1-hexene 分子放在这些位置时，双键 π 轨道可能与一个 Fe 原子相互作用，也可能跨在两个 Fe 原子之间形成更分散的金属-碳相互作用。配位数、吸附高度和分子轴方向发生变化后，优化过程会进入不同能量盆地。

图2. Fe(110) 表面的非等价吸附位点。

在 Si 掺杂重构 C(001) 这类缺陷或掺杂表面上，位点差异会放大到局域应变和电子态变化中。Si 原子附近既有掺杂原子引起的局域应变，又有相邻 C 原子的电子态变化。这类体系不适合只拿一个“看起来靠近活性中心”的初始构型做机制判断；分子靠近 Si 原子、相邻 C 原子或远离掺杂区域，电荷密度重排的空间区域不同，吸附能排序也可能换位。

图3. Si 掺杂重构 C(001) 表面的候选吸附位点。

分子本身的对称性、柔性和官能团分布，也会把同一个吸附物拆成多个构型族。线性小分子主要改变 端基朝向和倾角，柔性有机分子还会改变碳链构象、官能团接近表面的顺序和分子内扭转角。1-hexene 中 C=C 键、烷基链和端位氢原子对应不同相互作用区域，同一分子轴旋转几十度，就可能让双键从主吸附片段变成侧向接触片段。

图4. 1-hexene 分子的原子编号与双键位置。

局域极小值并不等价于计算错误，优化后得到多个构型，说明势能面上存在多个可停留的吸附态。真正要检查的是这些吸附态之间的能量差是否超过热扰动量级，关键键长是否合理，是否出现分子断键、表面原子异常拔出或 slab 整体漂移。若两个构型能量只差几十 meV，还应谨慎比较泛函、vdW 修正、覆盖度和偶极修正造成的误差范围。

构型初筛和全优化为什么都要做？

候选构型数量一多，直接对每个模型做完整优化会消耗大量机时。一个 4×4 金属表面，如果选 8 个非等价位点、12 个旋转角、2 个吸附高度，已经得到 192 个初始结构；再把分子换成柔性含氧有机物或含多个吸附原子的中间体，搜索空间会更大。构型初筛的作用是用较低成本排除明显高能或不合理的起点。

初筛用于缩小候选集合，全优化用于得到可分析结构，两者回答的是不同问题。初筛通常固定部分自由度或采用较短弛豫，用来判断哪些构型值得继续算；全优化则要让分子和可动表面原子充分响应，得到可用于 Eads、键长、PDOS、差分电荷密度和反应路径的结构。只用初筛能量做最终排序，容易低估表面重排、分子弯曲和吸附键形成带来的能量变化。

图5. 1-hexene 在 Fe(110) 上优化后的多组吸附构型与吸附能。

多构型优化后的结果会形成一个能量排序表，并把势能面上的相邻能量盆地暴露出来。最低能构型通常作为后续电子结构分析的主要模型，但偏高能构型不宜随手删掉。这些信息比单个最低能数字更能反映吸附势能面的形状：分子是否总会回到同一吸附姿态，哪些起点会转化为同一个终态，哪些构型只比最低能高 0.05-0.15 eV，哪些构型属于明显不稳定或受边界条件影响的状态。

构型筛选的输出应是一组可以复查的模型状态，不是一堆图片。全优化后的结构要看力收敛，更要看化学图像。若分子从 bridge 位滑到 hollow 位，说明输入坐标并没有保持为独立终态；若表面金属原子被拉出过多，可能需要扩大 slab 或检查覆盖度；若吸附分子发生断键，后续应区分吸附态和反应态。

对 *CO、*OOH 和 *CHO 这类催化反应中间体，多个吸附构型会直接改变可连接的反应路径。*CO 的 C 端吸附和 O 端吸附对应不同的成键方式，*OOH 的 O-O 轴向和表面金属位点关系会影响 O-O 键长，*CHO 的 C 端靠近金属还是 O 端靠近金属，会影响后续加氢位点。NEB 起点和终点必须来自同一反应通道内的相邻状态，前序吸附态没有筛干净时，能垒比较很容易被构型差异干扰。

吸附分析为何需匹配对应构型？

最低能结构可以作为主线模型，但机制分析不能只停在 Eads 排序。要看分子与表面的相互作用到底来自哪里：吸附键长、吸附高度、关键键角、差分电荷密度等值面、PDOS 轨道峰位和 Bader/Hirshfeld 电荷应围绕同一个优化构型展开。若电子结构图来自构型 A，而反应能垒来自构型 B，文章里的机制解释就会错位。

图6. 1-hexene 在 Fe(110) 上的最低能吸附构型侧视图。

侧视结构更适合检查吸附高度、分子弯曲和表面原子位移。对于范德华吸附，分子可能整体保持较远距离，Eads 主要受色散相互作用影响；对于化学吸附，某些键长会明显缩短，表面原子可能向上松弛。同一个能量值背后可能是完全不同的结构来源，这也是多构型比较必须结合几何参数的原因。

图7. 1-hexene 在 Fe(110) 上的最低能吸附构型俯视图。

俯视结构用来确认吸附位点和分子轴方向。一个构型在侧视图里看起来差别不大，俯视后可能已经从 top 位滑到 bridge 位，或者从某个金属原子上方转向相邻空隙。位点归属不能只靠初始命名，应以优化后最近邻原子、键长和分子中心位置来判定。

多构型数据进入机制分析时，建议保留三个层面的记录：第一，初始构型怎样生成，包括位点、旋转角和初始高度；第二，优化后是否汇聚到同一类终态；第三，最低能结构与次低能结构的能差是否足够大。若能差很小，机制解释要写成“在该模型下若选取最低能构型，电子结构信号表现为……”，而不是把该构型写成所有条件下唯一吸附方式。

对有电位、溶剂或覆盖度影响的界面体系，多个构型测试还要与模型边界配套。真空 slab 中最低能的平躺构型，在显式水层或电场作用下可能被氢键网络、双电层取向和离子分布改变；低覆盖度下稳定的单分子吸附，在高覆盖度下可能受分子-分子排斥或偶极相互作用影响。构型筛选要与覆盖度、溶剂模型、外加电场和参考态共同标注，后续图谱解释才不会越过模型范围。

同一个分子测试多个吸附构型，是在同一套 slab、覆盖度、泛函、vdW 修正和弛豫条件下采样势能面。吸附构型进入机制分析的前提，是把最低能构型、次低能构型和异常构型对应到后续判断：Eads 反映哪个结构的能量，PDOS 分析的是哪个成键片段，差分电荷密度对应哪个吸附高度，NEB 路径连接的是哪两个终态。