化学反应模拟：从ReaxFF反应力场到Kinetic Monte Carlo的多尺度建模路径

化学反应模拟在催化、燃烧和材料降解领域面对的是多时间尺度难题。化学键的断裂和生成发生在飞秒级，但反应器层面的转化率变化在秒到小时级。跨越10¹⁵的时间尺度不能用一把尺子丈量——需要DFT（飞秒→皮秒）→ReaxFF分子动力学（皮秒→纳秒）→动力学蒙特卡洛KMC（纳秒→秒/小时）的多尺度串联。每往上一层，牺牲一些原子细节，换取时间尺度的跨越。

DFT：势垒的基准精度

反应模拟的起点是基元反应势垒的DFT计算。过渡态搜索用CI-NEB（Climbing-Image Nudged Elastic Band）或Dimer方法，在PBE-D3(BJ)级别上得到势垒Ea和反应热ΔH。团队在一个CH₄催化部分氧化（CPOx）的项目里，用VASP计算了甲烷在Rh(111)表面C-H键解离的活化能：第一个C-H（CH₄→CH₃+H）Ea≈0.72 eV，第二个C-H（CH₃→CH₂+H）Ea≈0.52 eV，反应热分别是+0.12 eV和-0.26 eV。

DFT计算给了基元步骤的势垒和反应能，这是后面所有模拟的输入数据。精度要求不低：如果势垒误差>0.1 eV，反应速率在300°C下的Arrhenius外推就会偏差5倍以上——足以让KMC预测的反应转化率完全偏离实验。

ReaxFF：反应网络的发现

DFT给出的是假设好的基元反应路径上的势垒。但如果反应网络本身不完全——因为实际催化表面上的反应路径远比直觉设想的复杂——ReaxFF反应分子动力学可以用来”发现”未被假设的反应通道。

ReaxFF力场的核心是键级依赖（bond-order dependent）的势能表达式，允许化学键在MD过程中自然断裂和生成。对Rh(111)表面CH₄氧化的ReaxFF模拟（1000 K、覆盖度0.25 ML、500 ps），观察到的反应事件里，CH₃+O→CH₃O（甲氧基形成）和CH₂+O→CH₂O（甲醛形成）两条路径同时存在——这在初期的DFT网络假设中只包括了CH₃O路径。ReaxFF跑完后再补做了CH₂O路径的DFT过渡态计算，验证了它在热力学和动力学上的可行性。

这就是ReaxFF的多尺度角色：不是替代DFT算势垒（精度差一个量级），而是在DFT提供的势垒网络上”跑一遍”来发现遗漏的反应通道。发现通道→用DFT精算势垒→更新KMC的反应网络，这才是ReaxFF的正确用法。

KMC：宏观动力学的桥接

当反应网络包含30个以上基元步骤时（CH₄氧化在Rh上的网络约36个基元反应），用微观动力学方程（microkinetic modeling，MKM）的ODE做数值积分已经是常规做法。但MKM假设所有表面位点性质相同（平均场近似），对于催化表面上的空间非均匀性（比如台阶位和平台位的不同反应性）无法捕捉。

KMC不需要平均场假设——它在真实或近真实的表面上模拟每个基元反应的随机发生，反应速率由DFT给的速率常数控制。团队在Rh(111)用MKM和KMC分别计算了CH₄转化率-pH₂O分压曲线，MKM预测转化率在pH₂O>0.3 atm后明显高于KMC——差异来自MKM忽略了H₂O对Rh台阶位（高活性位点）的优先占据作用。KMC捕捉到了这个空间位阻效应，预测结果与实验流反应器数据一致性明显更好。

多尺度耦合的精度传递

化学反应模拟的精度从DFT势垒→ReaxFF反应网络→KMC宏观动力学是逐级传递的。DFT的势垒误差（~0.1 eV）在经过Arrhenius指数传播后在KMC层面放大约5-10倍——这是多尺度模拟的”误差放大器”效应。应对策略是在DFT阶段用高精度泛函（HSE06或RPA）做关键路径的势垒校正、在KMC阶段做DFT势垒±0.05 eV的敏感性分析来评估宏观预测的不确定区间。

更多内容请访问 https://www.keyanxueshu.com/