从限域效应到限域催化：自由能、过渡态与产物选择的计算解析

说明：本文华算科技主要介绍限域空间怎样改写催化自由能、局部覆盖度、过渡态、产物选择与周转条件，以及材料计算怎样辨认其中的几何、电子与传质贡献。

限域效应与限域催化的基本概念怎样进入自由能

反应物进入孔道、纳米笼、碳管内腔或壳层包围的空腔后，平移、转动和构象采样都受到空间约束。材料计算中的限域效应可定义为同一化学状态 X 在受限空间与开放表面之间的自由能差：ΔG_conf(X)=G_X^conf−G_X^open。当这种差值改变反应物富集、活化自由能、产物离开速率或路径分支时，体系呈现限域催化特征。

ΔG_conf可按物理来源写成孔壁相互作用、分子或骨架形变、溶剂化与局域电势变化、构象熵损失等项。孔壁色散与静电作用可能降低某个吸附态的能量，狭窄空间引起的键角弯曲则增加形变能；可用构型数量减少会使 −TΔS 项升高。每一项都依赖同一温度、化学势、覆盖度和电荷状态，参考条件改变后，数值不可横向拼接。

CoPc 位于碳纳米管外表面或内腔时，金属中心没有换成另一种元素，空间位置却改变了 CO 的驻留时间、CoPc 的几何形变和 ^*CO 的后续反应机会。这类内外对照把活性中心化学组成尽量固定，使 ΔG_conf指向碳管曲率、空腔尺度与局部浓度。若内外模型还伴随不同缺陷、质子化状态或分子聚集，能量差中会混入新的化学变量。

催化自由能还要与可观测速率对应。单分子实验可分别拟合吸附平衡常数 K_AR、表观反应常数 k_eff和孔长依赖传质速率。2.2 nm 孔道中，外壳加长使 K_AR下降而 k_eff升高；壳层继续增厚后，整体反应速率受扩散拖慢。吸附、化学转化与传质在同一限域结构中可以朝不同方向变化。

孔内覆盖度和过渡态为何会一起变化

线性纳米孔对反应物接近位点时的姿态产生筛选，使平面分子的长轴偏向孔轴方向，平躺吸附所需的构象采样空间随之减少。取向分布 P(θ)、吸附高度、关键二面角和分子—孔壁最近距离可量化这种预组织。姿态靠近反应坐标时，反应物进入过渡态所需的几何重排减小；姿态受到强迫扭曲时，初态或鞍点会积累形变能。

局部覆盖度由吸附自由能、孔内化学势和迁移速度共同控制。周期 DFT 的单分子构型只给出低覆盖度切片；覆盖度序列、从头算热力学或恒化学势采样可描述位点占据随压力、电位和温度的变化。平均覆盖度与瞬时邻近分布承担不同角色：前者决定空位比例，后者决定两个中间体能否在同一时刻靠近并发生双分子反应。

CoPc 进入不同内径的碳纳米管后，N—Co—N 夹角偏离平面构型，活性中心本身随受限空间发生形变；Co K 边 1s→4p_z跃迁强度随内径出现对应变化。DFT 得到的 ^*CO 吸附能也呈非单调趋势。孔径不是单调活性旋钮：极窄内腔可能阻止分子进入，适中曲率可调节配位几何，较宽内腔的形变逐渐减弱。

活化自由能的限域贡献写成 ΔΔG^‡=[G_TS^conf−G_IS^conf]−[G_TS^open−G_IS^open]。孔壁若对过渡态的稳定幅度超过初态，ΔΔG^‡为负；初态吸附获得更强稳定而过渡态收益较小，能垒反而升高。初态与过渡态必须采用同一覆盖度和参考态，否则差值会夹带额外分子或化学势贡献。

限域怎样改变产物选择与周转条件

中间体在孔内停留时间增加后，原本容易脱附的物种可继续反应。CoPc 在开放表面主要把 CO₂还原为 CO；碳管内腔提高 CO 的局部浓度，并使弯曲的 CoPc 加强 ^*CO 相互作用。原位 XANES 的 1s→3d 跃迁和 ATR-SEIRAS 的约 1920 cm⁻¹信号记录了电位下的 ^*CO 覆盖。驻留时间与吸附能共同设定中间体库存，两者会改变深度还原分支获得反应物种的概率。

在 U=−1.0 V_RHE、pH=6.8 的计算条件下，^*CO 继续质子—电子转移可通向 ^*CHO、^*CH₂O、^*CH₂OH 和 CH₃OH。CH₄路径中的 C—O 断裂具有较高能量代价。选择性来自竞争路径的相对自由能，孔内富集改变的是各路径实际经历的覆盖度和化学势，单条 0 K 总能曲线无法承载这两项环境变量。

蛋黄—壳层 Au@Cu₂O 纳米反应器给出了覆盖度改写路径的另一种情形。壳层过厚时，CO₂到达 Au 核受到抑制；中等壳厚兼顾供料与 CO 驻留，^*CO 氢化更有利；薄壳与小腔形成较高 ^*CO 覆盖，C—C 偶联比例上升。Cu(111) 上 ^*CHO 与 ^*OCCOH 形成能随覆盖度呈相反的非线性变化，3/16 与 6/16 ML对应不同的 CH₄/C₂H₄自由能排序。

周转频率还受产物移除控制。强孔壁作用可延长关键中间体寿命，也可能让产物占据位点；长孔道提高碰撞与再吸附概率，扩散时间则按孔长、曲折度和物种尺寸增长。微观动力学模型应把表面反应速率常数、吸附平衡、覆盖度和扩散通量写入同一反应网络。速率峰值常位于供料、活化和移除三类时间尺度相近的区间，而非孔径最小或吸附最强的一端。

怎样确认催化差异来自限域本身

受限样品与开放样品常伴随位点配位、负载量、表面羟基和可接近面积变化。Co 单原子体系中，介孔壳层既改变空间位置，又把主要 Co—N 配位从 Co—N₃移向 Co—N₄。PMS 吸附电荷重排、吸附能、原位 FT-IR 和生成 ¹O₂的自由能路径都随样品变化。几何约束与位点化学必须分别计量，否则“限域贡献”会包含配位数和表面基团的作用。

计算归因可设置两组互补对照。原子组成、活性中心、覆盖度和电荷状态保持一致，只改变孔壁距离或曲率，可得到纯几何约束项；几何坐标保持接近，只替换孔壁官能团或表面电荷，可读取孔壁化学项。每组模型对吸附态、过渡态和产物态采用相同超胞、泛函、色散修正与溶剂处理，能量差才对应被改变的变量。

静态 DFT 给出构型能量与电子结构，AIMD 提供孔内溶剂、氢键和短时构象涨落，增强采样得到沿孔轴或反应坐标的平均力势，经典或机器学习势分子动力学可延伸扩散时间尺度。四类量对应四个物理对象：ΔG_conf描述状态稳定，ΔΔG^‡描述化学步骤，驻留时间分布描述局部库存，扩散系数与平均首达时间描述供料和产物移除。

实验对照还应固定活性位数量与工作态组成，并报告按位点归一的速率。孔壁阻挡探针分子时，化学吸附测得的位点数可能偏低；高负载使纳米反应器彼此接触后，粒子间传质会改写单腔行为。在相同电位、pH、覆盖度与表面构型下，孔内外过渡态自由能差、^*CO 驻留分布和产物扩散势垒分别限定路径选择、局部库存与位点再生速率。