基于HEAT协议与二维主方程的高精度理论研究：C2H4 + OH反应的动力学机制与星际介质新通道

《Physical Chemistry Chemical Physics》：A high-level theoretical study of the C2H4 + OH reactionOpen Access

　　

乙烯（C2H4）与羟基自由基（OH）的反应在碳氢化合物燃烧和大气化学中均发挥着至关重要的作用。在本理论研究中，研究人员采用高精度HEAT协议，结合实用的二维主方程（2DME）分析和半经典过渡态理论（SCTST）对该反应进行了系统研究。计算得到的HEAT反应焓与基准ATcT值表现出极佳的一致性，偏差完全落在HEAT协议对稳定分子约0.2 kcal mol?1的名义不确定度范围之内。研究确认了两种主要反应机制：OH加成和直接H夺取。在环境大气条件下，OH加成路径生成热稳定化产物HOCH2CH2为主要通道，而直接H夺取在燃烧温度下占主导地位，两种机制的转变发生在约800–850 K。通过将OH加成步骤内过渡态（TS1）的HEAT计算能量调整?0.32 kcal mol?1，几乎所有可用的室温下压力相关衰减曲线（1–1000 Torr）实验数据集，以及作为温度函数（600–2000 K）的压力相关（0.01–100 atm）速率数据均被成功复现，误差一般在20%以内。通过加成/消除机制生成的乙烯醇在高温火焰条件下被表征，其估算产率为10 ± 2%，与文献值高度吻合。最后，通过快速隧穿生成甲醛 + CH3自由基和乙醛 + H原子的通道被首次从理论上预测并表征，被认为是星际介质极低温条件下最主要——甚至可能是唯一的反应路径。

论文解读：C₂H₄ + OH反应的高精度理论综合研究

一、研究背景与科学问题

乙烯（C₂H₄）是化石燃料燃烧和大气碳氢化合物光氧化过程中的关键中间体。其与羟基自由基（OH）的反应在多种环境中充当C₂H₄的主要汇，通过多条路径生成不同的初级产物。在对流层条件下，该反应是甲醛（CH₂O）、一氧化碳（CO）和对流层臭氧（O₃）等二次污染物形成的重要途径；在内燃机（ICEs）等高温环境中，它是包括甲烷在内的碳氢化合物复杂燃烧机制中的关键中间步骤；而在星际介质（ISM）的极低温度条件下，该反应可能有助于甲醛和乙醛（CH₃CHO）的形成。

尽管该反应已被广泛研究，但以往基于不同理论水平的势能面（PES）研究存在典型0.5–1 kcal mol^?1的不确定度，且与实验吻合需要将OH加成路径关键过渡结构（TS1）的能量调整高达1.3 kcal mol^?1。近期Li和Long使用GMM(P).L方法改进了PES，特别是TS1能量，实现了对该系统约0.2 kcal mol^?1的估计精度，但未考虑燃烧相关的高温条件和星际介质的极低温条件。因此，亟需采用更高精度的理论方法对该反应进行全面、系统的重新研究，以覆盖从星际介质极低温到燃烧高温的宽温度范围。

二、研究内容与主要结论

研究人员采用高精度外推从头算热化学（HEAT）协议，结合实用的二维（E,J）分辨主方程（2DME）分析和W. H. Miller的半经典过渡态理论（SCTST），对C₂H₄ + OH反应的最低双重态势能面进行了全面映射。计算得到的HEAT反应焓与基准ATcT值表现出极佳的一致性，偏差小于0.15 kcal mol^?1。研究确认了两种主要反应机制：OH加成（路径I）和直接H夺取（路径II）。在环境大气条件下，OH加成生成热稳定化产物HOCH₂CH₂（INT1）为主要通道；直接H夺取在燃烧温度下占主导，两种机制的转变发生在约800–850 K。通过将TS1的HEAT计算能量调整?0.32 kcal mol^?1，几乎所有可用的室温压力相关衰减曲线（1–1000 Torr）实验数据以及压力相关（0.01–100 atm）速率数据（600–2000 K）均被成功复现，误差一般在20%以内。通过加成/消除机制生成的乙烯醇在高温火焰条件下被表征，估算产率为10 ± 2%。甲醛 + CH₃自由基和乙醛 + H原子的生成通道被首次从理论上预测并表征为星际介质极低温条件下最主要——甚至唯一的反应路径。

三、关键理论方法

研究人员主要运用了以下关键技术方法：第一，高精度HEAT协议，采用全电子耦合簇方法ae-CCSD(T)/cc-pCVQZ进行几何优化，对动力学关键的TS1采用更高精度的HEAT-456(Q_Λ)方法，利用(Q_Λ)变体的微扰四重激发校正以提高对过渡态部分键的描述精度；第二，二维（E,J）分辨主方程（2DME）分析，用于处理压力敏感的OH加成路径，求解能量和角动量分辨的动力学弛豫矩阵；第三，半经典过渡态理论（SCTST），基于VPT2框架内完全耦合的非简谐振动，自动包含多维隧穿修正，用于计算紧致内过渡态的振转态和；第四，微变分RRKM理论，用于表征无垒OH + C₂H₄ → PRC缔合步骤中松散的变分外过渡态TS0；第五，一维受阻内转动（1D-HIR）处理，对TS1、TS2、TS5、TS6、INT1和INT2中的大振幅扭转运动单独求解薛定谔方程获得量子态。

四、研究结果

4.1 势能面（PES）

研究人员采用HEAT-345(Q_Λ)方法构建了最低双重态势能面，对动力学关键的TS1采用HEAT-456(Q_Λ)方法。HEAT计算结果与ATcT基准值偏差小于0.15 kcal mol^?1。与文献值对比发现，Cleary等人（CBS-APNO//BHandHLYP/6-311+G(3df,2p)水平）的结果偏差最大，高估INT1和INT2的稳定性分别超过1和2 kcal mol^?1；Ali和Barker（CCSD(T)/aVTZ//BHandHLYP/aVTZ水平）的反应势垒与HEAT值吻合较好，但某些产物能量不准确；Senosiain等人（RQCISD(T)/CBS//(U)QCISD(T)/6-311G(d,p)水平）的结果总体吻合良好，但预测TS1为+0.2 kcal mol^?1的正势垒，而本研究预测为?0.68 kcal mol^?1的浸没势垒，与Cleary和Ali和Barker的结果一致。

4.2 TS1能量分解分析

研究人员将HEAT-456(Q_Λ)方法与GMM(P).L方法计算的TS1能量分解项进行了对比。总差异为?0.21 kcal mol^?1，来源于三个具体贡献：自旋轨道耦合（?0.09 kcal mol^?1）、零点能ZPE（?0.08 kcal mol^?1）和CCSDTQ(P)/CBS校正（?0.05 kcal mol^?1）。自旋轨道校正采用Van Vleck方法评估，ZPE差异源于TS1中1D-HIR的处理方式，而?0.05 kcal mol^?1的差异是两种模型的内在差异。微扰五重激发（P）被省略因其对TS1能量贡献可忽略。

4.3 动力学计算与实验对比

研究人员采用“双过渡结构”（two-TS）模型处理浅的范德华预反应复合物（PRC），有效态和通过虚拟入口TS01计算，绕过了PRC阱。通过将TS1能量调整?0.32 kcal mol^?1，室温下1–1000 Torr的压力相关衰减曲线和600–2000 K温度范围0.01–100 atm的压力相关速率数据均被成功复现，误差一般在20%以内。

4.4 产物分布与星际介质通道

乙烯醇通过加成/消除机制在高温火焰条件下生成，估算产率为10 ± 2%，与文献值吻合。甲醛 + CH₃自由基和乙醛 + H原子通过快速隧穿生成的通道被首次从理论上预测和表征，被认为是星际介质极低温条件下最主要——甚至唯一的反应路径。

五、讨论与结论

研究人员指出，主方程模拟在约1 μs时终止，以防止热稳定化的HOCH₂CH₂热再解离回反应物掩盖初级动力学。在低至中等温度下可使用汇收集碰撞稳定化的中间体，但在高温下热稳定化中间体可被碰撞再活化发生分解，因此汇不适用。能量和角动量转移概率的可分离性假设以及平均向下转动能量转移等于平均向下振动能量转移的假设，对计算衰减曲线影响适中，差异可通过适当选择〈ΔE_vib〉_d补偿。对称陀螺近似对速率常数的影响预计在约1%以内。TS1中两个镜像鞍点间的扭转势垒仅3 cm^?1，远低于该HIR模式的ZPE 24 cm^?1，因此不应作为速率倍增的对映体处理。

研究结论表明，HEAT协议结合2DME分析和SCTST能够以约0.2 kcal mol^?1的精度描述C₂H₄ + OH反应的完整动力学，在宽温度（极低温至燃烧温度）和压力范围内与实验数据良好吻合。该研究首次从理论上预测了星际介质条件下甲醛 + CH₃和乙醛 + H的生成通道，为理解星际化学中复杂有机分子的形成提供了新的理论依据，同时也为燃烧化学和大气化学中的反应动力学模型提供了高精度的速率常数数据。