Cl(2P3/2) + C2H5I反应多通道路径的高精度从头算（ab initio）表征

《Physical Chemistry Chemical Physics》：High-level ab initio characterization of the multichannel Cl(2P3/2) + C2H5I reaction

　　

Cl(²P_3/2) + C₂H₅I反应的势能面（PES）通过高精度电子结构计算进行描述，涵盖氢抽提和碘抽提路径，以及几条经由瓦尔登翻转（Walden inversion）或正面进攻（front-side attack）进行的取代路径，包括原子（H、I）和基团交换（CH₂I、CH₃）机理。所有稳定点的几何结构和简谐振动频率在MP2/aug-cc-pVDZ和CCSD(T)-F12b/aug-cc-pVDZ理论水平下确定，单点能在最准确几何结构下使用耦合簇方法结合aug-cc-pVTZ和aug-cc-pVQZ基组进一步精修。为达到化学精度，五项额外能量修正——考虑核相关、标量相对论、自旋-轨道和后CCSD(T)效应——被纳入CCSD(T)-F12b/aug-cc-pVQZ单点能中。所得基准数据实现了反应路径的详细绘制，包括过渡态以及反应前和反应后极小点的识别，这些物质引导系统从反应物向势能面上的各个产物通道演进。速率系数使用过渡态理论（包括Wigner隧穿修正）确定，并与文献理论值进行比较。

海洋大型藻类和浮游植物的代谢活动释放多种含碘有机化合物进入沿海和开阔海洋大气，其中包括碘乙烷（C₂H₅I）等小分子碘化挥发性物质；这些碘化物种通过两条不同路径间接影响臭氧损耗：一条是近紫外吸收导致高量子产率碘原子释放的光解离路径，另一条是由氯原子等多种大气氧化剂引发的氧化降解路径。现有关于Cl + C₂H₅I反应的研究主要集中于动力学视角，对其底层能量学或稳定结构信息涉及甚少。该研究通过高阶从头算方法全面绘制了反应势能面，为理解这一重要大气化学反应提供了精确的电子结构基础。

前人动力学研究表明，Cl + C₂H₅I反应存在两条主要路径：Cl与C₂H₅I结合形成分子加合物（C₂H₅I?Cl）的通道，以及生成HCl和C₂H₄I·自由基的H-抽提通道。低温时加合物形成被有效稳定，而高温时快速预解离使H-抽提成为主导路径。然而，除加合物表征外，反应路径的理论工作十分有限，且既往研究仅报告了连接反应物与产物的少数稳定点。鉴于现有研究在能量学和稳定结构信息方面的不足，研究人员开展了这项高精度从头算研究，旨在提供基于基准质量几何结构和相对能的准确描述。

该研究发表于《Physical Chemistry Chemical Physics》期刊，其核心意义在于：首次全面绘制了Cl + C₂H₅I反应的多通道势能面，涵盖抽提和取代两大类反应路径，为后续动力学建模和轨迹模拟研究奠定了坚实基础。

研究人员采用以下关键技术方法：几何优化与频率计算在MP2/aug-cc-pVDZ和CCSD(T)-F12b/aug-cc-pVDZ水平进行；单点能计算使用CCSD(T)-F12b方法结合aug-cc-pVTZ和aug-cc-pVQZ基组；为达到化学精度，系统考虑了五项能量修正——核相关修正（Δ_core）、标量相对论修正（Δ_rel）、自旋-轨道耦合修正（Δ_SO）、CCSDT修正（δ[T]）和CCSDT(Q)修正（δ[(Q)]）；自旋-轨道修正采用多参考构型相互作用方法（MRCI）结合Davidson校正（+Q）使用Breit-Pauli算符计算；基准经典能量通过将上述修正加入CCSD(T)-F12b/aug-cc-pVQZ单点能获得；基准绝热能量则进一步加入零点振动能（ΔZPE）修正；动力学计算采用标准过渡态理论（TST）并包含Wigner隧穿修正；所有计算使用MOLPRO程序包（2015.1版本及2023.2版本）和MRCC程序完成。

**入口通道特征**

研究人员对Cl + C₂H₅I反应的入口通道进行了九种不同取向的一维势能曲线扫描，采用MP2和CCSD(T)-F12b方法结合aug-cc-pVDZ基组。仅两种取向表现出深于1.5 kcal mol^-1的明显势阱，对入口通道动力学有显著影响。其一为Cl原子从碘原子侧垂直于C–C键接近C₂H₅I分子，形成约7 kcal mol^-1（MP2水平）或近11 kcal mol^-1（耦合簇水平）的极小点，C–Cl距离约4 ?；其二为Cl原子沿C–C键轴向碘原子连接的碳原子接近，类似深度约为7 kcal mol^-1（MP2）或约10 kcal mol^-1（CCSD(T)-F12b），C–Cl距离约3 ?。除一种取向外，MP2和耦合簇方法的势阱位置基本一致；在两种最有利接近方式中，耦合簇方法所得的势阱深度比MP2深约4 kcal mol^-1。

**抽提路径**

碘抽提（IABS）路径中，研究人员通过一维扫描发现了与明显势阱相关的构型，其中反应物侧极小点（IABS preMIN）经几何优化和频率分析确认为真正的势能面极小点，基准经典（绝热）相对能为-12.85（-12.27）kcal mol^-1。该路径上定位了两个过渡态：IABS TS1为Cl原子沿C–C键轴上方接近的直线构型，IABS TS2为从侧面指向碘原子的外部进攻构型。较低能垒对应IABS TS2，经典（绝热）相对能为11.47（10.11）kcal mol^-1，而IABS TS1显著较高，为20.60（19.88）kcal mol^-1。产物侧存在两个分离极小点IABS postMIN1和IABS postMIN2，后者仅略高于前者0.44（0.35）kcal mol^-1。最终产物ICl和乙基自由基（C₂H₅·）的能量比反应物高11.32（8.01）kcal mol^-1，表明该过程为吸热反应。实验测得的IABS preMIN形成焓与计算结果吻合良好，特别是与第三定律实验值-12.67 ± 1 kcal mol^-1最为接近。

氢抽提（HABS）可在三个不同位点发生。抽象碘原子直接相连碳上的氢（HABS TS3）能垒最低，经典（绝热）相对能为1.46（-2.65）kcal mol^-1；从与碘取代基同侧的甲基抽氢（HABS TS1）的能垒略高，为2.03（-2.27）kcal mol^-1；反式取向抽氢（HABS TS2）能垒最高，为5.36（0.61）kcal mol^-1。考虑零点振动能后，三个通道中的两个已无能垒。产物侧存在对应的三个后反应极小点，其能量排序与过渡态不同。根据抽氢位点不同，产物包括CH₃CHI和CH₂CH₂I两种类型（均伴随HCl）。甲基抽氢产生的CH₂CH₂I存在两种构象：碘原子位于C–C键上方的桥式（Bridged）构型，以及碘原子向外的反式（Anti）构型。Bridged-CH₂CH₂I + HCl为最低能量产物态，经典（绝热）相对能为-8.25（-12.56）kcal mol^-1；Anti-CH₂CH₂I + HCl为-1.75（-6.18）kcal mol^-1；CH₃CHI + HCl最高但仍为放热，为-0.30（-5.10）kcal mol^-1。CH₂CH₂I分子可能解离生成C₂H₄ + I + HCl三碎片产物，相对能为-7.13（-12.12）kcal mol^-1。

**取代路径**

四种取代路径的特征如下：碘取代（I-substitution）是唯一热力学上有利的取代过程，产物CH₃CH₂Cl + I比反应物低28.71（27.98）kcal mol^-1，是所有研究反应中最稳定的产物路径。该路径同时具有最低取代能垒：瓦尔登翻转过渡态（IS W TS）为13.59（13.13）kcal mol^-1，正面进攻过渡态（IS FS TS）为20.59（20.00）kcal mol^-1；正面进攻路径还存在全局最低的后反应极小点IS FS postMIN，为-30.35（-29.gon kcal mol^-1。

CH₂I取代为最深吸热取代通道，产物CH₂I + CH₃Cl位于反应物上方9.42（5.73）kcal mol^-1；两个过渡态分别为瓦尔登翻转（CH₂I S W TS，34.29(32.60) kcal mol^-1）和正面进攻（CH₂I S FS TS，57.73(55.92) kcal mol^-1），瓦尔登路径上还存在后反应络合物CH₂I S W postMIN（7.79(4.50) kcal mol^-1）。

CH₃取代略比CH₂I取代吸热，产物在反应物上方14.93（10.32）kcal mol^-1；瓦尔登翻转（CH₃ S W TS）和正面进攻（CH₃ S FS TS）过渡态分别为43.42（41.11）和54.64（52.48）kcal mol^-1，对应后反应络合物能量为13.26（9.19）和13.84（9.59）kcal mol^-1。

H取代是最吸热的通道，产物为三种C₂H₄ClI型物质伴随H原子。反式-CH₂ClCH₂I + H最低，为21.79（16.53）kcal mol^-1；顺式异构体（cis-CH₂ClCH₂I）略高，为23.45（18.15）kcal mol^-1；两个卤素连于同一碳原子的CClHICH₃最吸热。能量顺序可归因于空间位阻效应：大体积卤素取代基靠近时产生增强排斥。反式路径上两个过渡态（HS W TS1为43.19（39.13），HS FS TS1为57.00（52.66）kcal mol^-1）和一个后反应极小点HS FS postMIN1（21.48(16.57) kcal mol^-1）；顺式路径类似，但后反应极小点不遵循过渡态的能量顺序。

**能量数据综合分析**

所有稳定点的能量数据汇总于表1中。MP2/aug-cc-pVDZ与CCSD(T)-F12b/aug-cc-pVDZ能量存在显著差异，最大偏差见于CH₂I S FS TS（8.32 kcal mol^-1）。aug-cc-pVDZ至aug-cc-pVQZ基组的耦合簇能量差异通常在1 kcal mol^-1内，三ζ至四ζ转变引入的变化常低于0.1 kcal mol^-1，表明显相关方法具有优良的基组收敛性。五项修正项的纳入对达到亚化学精度不可或缺：自旋-轨道效应最为显著（因含两个较重卤素原子，Δ_SO通常为0.8–1.0 kcal mol^-1，产物通道中可达6.32 kcal mol^-1，源于产物中自由碘原子的强自旋-轨道耦合）；其余四项修正通常低于1 kcal mol^-1，但δ[T]对IABS TS1和IABS TS2达约-1.1 kcal mol^-1，Δ_core对若干点超过1 kcal mol^-1。

与活性热化学表（ATcT）依据赫斯定律推导的反应焓比较显示：C₂H₄ + I + HCl通道的计算值-12.12 kcal mol^-1与ATcT值-12.49 ± 0.12 kcal mol^-1吻合在0.4 kcal mol^-1内；C₂H₅Cl + I通道的-27.98 kcal mol^-1与-28.22 ± 0.13 kcal mol^-1也符合良好；仅C₂H₅ + ICl通道差异较大，可能与ATcT数据本身的问题有关。

**动力学结果**

TST速率系数与Jia等人的ICVT/SCT结果比较显示显著偏差，特别是在低温区，随温度升高渐趋一致。传统TST处理所得较大速率系数可能部分源于隧穿效应的简化描述：Wigner修正仅基于过渡态虚频率，而SCT形式更为严格；此外，入口通道低能极小点（IABS preMIN）的存在使系统在相遇后可能花费非可忽视时间，结构重排可能使系统从碘抽提转向氢抽提通道，这些动力学效应未在传统TST中充分考虑；最后，TST忽略再穿越效应也是差异来源。尽管本研究电子结构水平更高，但鉴于ICVT/SCT对整个最小能量路径和隧穿效应的更严谨处理，文献值可能提供更可靠的动力学描述，尤其在低温区。

**研究结论**

研究人员通过高精度从头算方法全面绘制了Cl(²P_3/2) + C₂H₅I反应的势能面，系统研究了抽提和取代两类反应路径，包括H-和I-抽提通道，以及经由瓦尔登翻转或正面进攻进行的碘取代、CH₂I取代、CH₃取代和H取代过程。能量轮廓基于CCSD(T)-F12b/aug-cc-pVDZ优化几何结构和CCSD(T)-F12b/aug-cc-pVQZ单点能构建，并加入标量相对论、核相关、自旋-轨道和后CCSD(T)贡献修正。热力学上，碘取代和H-抽提通道最为有利，其中碘取代更深（-28.71/-27.98 kcal mol^-1）。尽管碘取代产物更稳定，H-抽提过程在动力学上更受青睐，其活化能垒更低。最低碘取代能垒（瓦尔登翻转机制）为13.59（13.13）kcal mol^-1，而最高H-抽提能垒（甲基反式抽氢）仍低8.23（12.51）kcal mol^-1。最低的H-抽提能垒位于碘取代碳上（1.46/-2.65 kcal mol^-1），甲基同侧抽氢次之（2.03/-2.27 kcal mol^-1）。碘抽提及CH₃-、CH₂I-、H-取代通道在热力学和动力学上均不利。入口通道浅势阱可能影响整体反应动力学，一维扫描鉴定了两个具有明显势阱的构型，其中之一优化为稳定的反应前络合物（-12.85/-12.27 kcal mol^-1）。基于最准确参数估计的H-抽提速率系数与传统ICVT/SCT结果存在分歧，有待未来研究解决。

该研究为Cl + C₂H₅I势能面的全面高精度表征，揭示了控制竞争反应路径的热力学和动力学特征，为后续动力学建模研究和势能面开发奠定了基础，将促进对该反应动力学的定量理解。