微观动力学需要哪些DFT数据？微观动力学数据清单与速率常数换算讲透

说明：本文主要介绍微观动力学所需的 DFT 结构、能量、振动频率与环境数据，以及这些数据怎样用于计算速率常数、表面覆盖度和产物生成速率。

微观动力学把吸附、表面转化、扩散和脱附写成一组基元反应速率方程，计算各中间体随时间的生成与消耗。DFT 为具体原子构型提供能量，统计热力学将能量和振动信息转化为自由能，速率理论再给出基元步骤的速率常数。实验测得的产物生成速率，来自整组方程求解后的净生成通量。

反应网络中的节点对应物种与吸附位点，连线对应基元步骤。乙二醇可沿 O–H 或 C–H 断键生成不同中间体，也可发生 C–C、C–O 断裂。每个保留分支都要有热力学与动力学描述；副产物、强吸附残留物和逆反应会改变表面组成，只计算目标产物的一条路径可能漏掉竞争反应。

在固定表面位点数的平均场模型中，覆盖度满足 dθ_i/dt = Σ_jν_ijr_j。θ_i 是物种数量与参考位点数之比，ν_ij 表示步骤 j 对物种 i 的净计量变化，r_j 为按相同位点基准归一的净速率。速率常数决定转化概率，活度和覆盖度决定参与反应的物种数量，空位也参与位点守恒。

Pd(111) 上乙醇分解的能量剖面可作为 DFT 数据示例；在 600 K、乙醇分压 1 atm 的模型中，CO 和含碳中间体的覆盖度随时间变化，产物脱附速率出现不同平台。图中的能量曲线属于计算输入，覆盖度与速率曲线属于动力学输出。该时间响应模型忽略横向相互作用，其时间数值仅作定性参考。

表面即使积累了低能中间体，也要经过后续转化或脱附才能再生空位。中间体自由能、过渡态自由能和位点占用必须共同进入方程，才可计算工作条件下的覆盖度。有限时间内的速率平台可能只是暂态平台，稳态计算还需检查各物种净生成速率及初始覆盖状态的影响。

对网络内的气相分子、吸附中间体和过渡态，需要保存优化结构及电子总能，并保留清洁表面的能量。吸附构型要带有晶面、位点、化学组成和磁态信息；过渡态要连接指定初末态，并由反应方向上的单个虚频及位移模式确认。只给一个势垒数值，会丢失它对应的反应物组合。

单步电子反应能为 ΔE = E_FS − E_IS，正向电子势垒为 ΔE^‡ = E_TS − E_IS；初末态要计入计量式中的全部分子与表面。各物种采用同一元素参考、清洁表面基准和计算方法后，才能组合不同基元步骤。原子守恒的闭合反应循环应恢复相同能量；气相能量若采用实验校正，相关步骤也要一致更新。

每个稳定态及过渡态的正常振动频率用于计算零点能、热焓和振动熵。常用近似写作 G°(T) ≈ E_DFT + E_ZPE + ΔH_thermal(T) − TS°(T)。吸附 O₂ 的内部 O–O 伸缩、分子与表面的伸缩及转动模式，对金属种类的响应不同，给所有吸附物套用一组频率会改变自由能差。

过渡态的反应坐标虚频不计入稳定振动配分函数，其余实频参与热修正。低频摆动或表面迁移若明显偏离谐振近似，可采用受阻转动、受阻平动模型，并补充对应势垒。气相热力学还涉及分子质量、转动惯量、对称数和自旋简并度，可由结构或热化学数据库提供。

在 CO/O 共吸附的示例中，采用 Pd(111) 的振动贡献作为固定值，会在部分吸附能区间低估 CO 覆盖度。色标表示覆盖度误差，横纵坐标分别为 O 与 CO 的电子吸附能；即使电子能数据保留不变，热修正近似也会改变可用空位。

对满足过渡态理论假设的热活化步骤，标准态速率常数常写作 k = κ(k_BT/h)exp[−ΔG^‡/(k_BT)]，κ 为透射系数，常规近似取 1。活化熵已包含在 ΔG^‡ 中时，应避免再乘入同一组振动贡献。气体吸附还可由分子碰撞通量与黏附系数描述；无吸附势垒也仍有分子到达表面的通量限制。

正逆速率常数采用一致的无量纲活度及标准态约定时，应满足 k_f/k_r = exp[−ΔG°_rxn/(k_BT)]。同一过渡态与两侧状态自由能可同时确定正逆势垒。独立拼接两套势垒可能破坏详细平衡，使模型在热力学平衡时仍产生非零净通量。

一个中间体可能占据多个表面位点。若 θ_i 按每个参考位点上的分子数定义，位点守恒为 θ_* + Σ_iq_iθ_i = 1，其中 q_i 是占位数。位点类别、占位数及邻接关系来自结构模型；纳米颗粒中台面、边缘和角位的数量比例，还需由粒径、形貌模型或表征约束。

Pt 纳米颗粒的台面、边缘和角位具有不同近邻数量，相邻吸附物的排斥或吸引会随局域排列改变。高覆盖度下，DFT 数据应包含不同覆盖度和共吸附排列的能量差，必要时拟合成对或多体相互作用。计算结果需区分平均吸附能与新增一个分子的微分吸附能；覆盖相关模型还要处理过渡态对邻近物种的响应。

同一颗粒的 Pt(111) 与 Pt(211) 可承担不同反应步骤。水煤气变换的双位点模型中，水解离偏向台阶位，羧基分解可更多发生于台面；跨晶面迁移使两类位点共享中间体。模型因此需要位点间扩散势垒与允许迁移的邻接关系，或有依据的快速扩散近似。

双位点模型的两组 CO 覆盖度、分步净速率和温度响应，用于检验跨位点耦合的影响。温度、分压与进料组成由反应条件给定；理想气体活度为 p/p°，化学势含 k_BT ln(p/p°)。气体压力已通过活度进入速率式时，避免在标准态自由能中重复计入相同压力项。

电催化模型还需要电位标尺、pH、溶液活度及界面溶剂状态。带电吸附物与过渡态对电位和界面电场的响应可能不同，单个电位下的静态势垒通常只适用于该界面状态。电位依赖势垒可来自多个电位下的计算或经检验的近似关系；电子数可变的模型要使用与电极电子库一致的热力学势。

能量表应区分直接 DFT 结果、经验校正和标度关系估计值。大网络可利用同类断键反应的能量关系估计部分势垒，但估计式的适用材料、反应类型与残差应随数据保存。Rh(111) 合成气反应的网络示例显示，考虑 DFT 能量误差后，为达到更高置信水平，需保留更多产物分支。

80% 与 90% 置信水平对应不同候选网络，反映能量误差对机理归属的影响。对单个步骤，在 500 K、前因子固定时，活化自由能增加 0.10 eV 会使速率常数约降为十分之一；整网产物速率的变化还受覆盖度反馈影响。同一泛函对相似中间体的误差常有关联，误差传播宜保留这种关联。

速率控制度可通过小幅改变某个过渡态自由能、重新求解覆盖度与产物速率获得。扰动过渡态时保持反应平衡常数不变，正逆速率随同一过渡态共同调整。对中间体能量的扰动则会改变吸附平衡及相关势垒。补算可优先针对高敏感度且误差较大的状态，以及可能改变主要产物分支的竞争步骤。

DFT 原始数据还应保留结构、能量参考、频率、泛函、自旋态、溶剂处理和收敛信息，使势垒差能追溯到具体模型。对照实验时，产物速率、反应级数、表观活化能和覆盖度分别约束不同反应行为。将每位点速率换算为实验质量活性需要活性位数量；存在传质或温度梯度时，还需反应器模型提供局域浓度与温度。

#微观动力学 #材料计算 #第一性原理 #DFT #过渡态 #活化自由能 #振动频率 #表面覆盖度 #催化 #科研