CF2ClC(O)OCH2CH3+OH的反應(yīng)機(jī)理及動力學(xué)性質(zhì)的理論研究

朱　鵬, 段雪梅, 劉靖堯

(吉林大學(xué)理論化學(xué)研究所, 長春 130021)

?

CF2ClC(O)OCH2CH3+OH的反應(yīng)機(jī)理及動力學(xué)性質(zhì)的理論研究

朱鵬, 段雪梅, 劉靖堯

(吉林大學(xué)理論化學(xué)研究所, 長春 130021)

摘要利用雙水平直接動力學(xué)方法, 在MCG3-MPWB//M06-2X/aug-cc-pVDZ水平上研究了CF2ClC(O)OCH2CH3+OH的微觀反應(yīng)機(jī)理. 得到了反應(yīng)物CF2ClC(O)OCH2CH3的5種穩(wěn)定構(gòu)象(RC1～RC5), 并對每一構(gòu)象考察了發(fā)生在—CH3和—CH2—基團(tuán)上的所有可能氫提取反應(yīng)通道. 利用改進(jìn)的變分過渡態(tài)理論(ICVT)結(jié)合小曲率隧道效應(yīng)校正(SCT)計算了各反應(yīng)通道的速率常數(shù), 分析了各構(gòu)象反應(yīng)位點選擇性. 結(jié)果表明, 對于構(gòu)象RC1和RC2, 低溫時氫提取反應(yīng)主要發(fā)生在—CH2—基團(tuán)上; 而對于構(gòu)象RC3, RC4和RC5, 發(fā)生在—CH3基團(tuán)上的氫提取反應(yīng)通道在整個溫度區(qū)間內(nèi)占絕對優(yōu)勢. 根據(jù)Boltzmann配分函數(shù)計算總包反應(yīng)速率常數(shù), 在298 K溫度下計算的體系總包反應(yīng)速率常數(shù)與實驗值相符, 進(jìn)而給出200～1000 K溫度范圍內(nèi)擬合了速率常數(shù)的三參數(shù)Arrhenius表達(dá)式: koverall=5.45×10-25T4.54exp(-685/T).

關(guān)鍵詞直接動力學(xué)方法; 速率常數(shù); 變分過渡態(tài)理論; 密度泛函理論; 反應(yīng)機(jī)理

由于氟利昂(CFCs)對大氣臭氧層的破壞以及引發(fā)溫室效應(yīng)已被禁止使用, 因此尋找環(huán)境友好型CFC替代物已成為實驗和理論研究的熱點. 近年來, 氟氯取代醚(HCFEs)[1,2]類化合物被認(rèn)為是CFC新一代替代品, 已引起廣泛關(guān)注. HCFE類化合物具有與CFC類似的物理化學(xué)性質(zhì), 且分子中的醚鍵能夠提高其在對流層中的反應(yīng)活性, 從而縮短其大氣壽命. 氟氯取代脂(CFESs)是HCFE化合物降解的主要產(chǎn)物, 同時也是大氣中重要的中間體. 由于降解產(chǎn)物CFES中含有氯原子, 釋放氯原子進(jìn)入同溫層后仍會破壞臭氧層. 在大氣中, CFES降解的主要途徑是與OH自由基的反應(yīng), 因此, 為了更好評估CFES對大氣環(huán)境的影響, 研究CFES與OH的反應(yīng)機(jī)理及動力學(xué)性質(zhì)是十分必要的.

本文主要關(guān)注CF2ClC(O)OCH2CH3與OH自由基的反應(yīng). 顯然, 由于反應(yīng)物CF2ClC(O)OCH2CH3含有較長的碳鏈, 因此通常具有不止一種穩(wěn)定構(gòu)象; 并且對于每個構(gòu)象與OH的反應(yīng), 存在從2個反應(yīng)位置(—CH3和—CH2—)發(fā)生的多條氫提取通道. 然而, 實驗上Blanco等[3]使用相對速率方法只給出了該反應(yīng)在298 K時的反應(yīng)速率常數(shù)[(5.4±1.5)×10-13cm3·molecule-1·s-1], 而對于其它溫度下的反應(yīng)速率常數(shù)及反應(yīng)位點選擇性和產(chǎn)物分支比等信息均未見報道. 因此, 從理論上研究反應(yīng)CF2ClC(O)OCH2CH3+OH的微觀機(jī)理和動力學(xué)性質(zhì)具有重要意義.

本文采用量子化學(xué)計算方法, 在MCG3-MPWB//M06-2X/aug-cc-pVDZ水平上, 得到了反應(yīng)物CF2ClC(O)OCH2CH3的5個穩(wěn)定構(gòu)象, 并確定了與之相關(guān)的共23條不同的反應(yīng)通道. 采用直接動力學(xué)方法[4～10], 應(yīng)用變分過渡態(tài)理論(VTST)[11,12]在較大溫度范圍內(nèi)計算了各條氫提取通道速率常數(shù), 進(jìn)而預(yù)測出總反應(yīng)速率常數(shù), 并分析每個構(gòu)象反應(yīng)位點的選擇性及其對總反應(yīng)的貢獻(xiàn); 此外, 在計算中利用等化學(xué)鍵反應(yīng)估算了反應(yīng)物和產(chǎn)物自由基的標(biāo)準(zhǔn)生成焓.

1計算方法

所有電子結(jié)構(gòu)計算均采用Gaussian 09[13]程序完成. 在M06-2X/aug-cc-pVDZ[14]水平上優(yōu)化各駐點的幾何結(jié)構(gòu), 并通過頻率分析確定反應(yīng)物、產(chǎn)物和中間體為局部極小點, 而過渡態(tài)有唯一虛頻. 通過內(nèi)稟反應(yīng)坐標(biāo)(IRC)計算得到各反應(yīng)通道的最小能量途徑(MEP)確認(rèn)過渡態(tài)連接了相應(yīng)的反應(yīng)物構(gòu)象和產(chǎn)物. 在MEP上選取適當(dāng)點進(jìn)行頻率分析, 獲得所選點的梯度和海森(Hessian)矩陣. 為了獲得更為準(zhǔn)確的能量, 使用多系數(shù)外推的密度泛函方法MCG3-MPWB[15]計算了穩(wěn)定點和MEP上4個非駐點的單點能量. 計算中采用0.9721作為M06-2X/aug-cc-pVDZ 的頻率校正因子. MCG3-MPWB計算在結(jié)合Gaussian 03程序[16]的MLGAUSS 2.0程序[17]下完成.

采用變分過渡態(tài)理論結(jié)合內(nèi)推單點能方法(ISPE)[18]獲得體系雙水平(MCG3-MPWB//M06-2X/aug-cc-pVDZ)勢能面. 在此基礎(chǔ)上, 利用改進(jìn)的變分過渡態(tài)理論(ICVT)[19]結(jié)合小曲率隧道效應(yīng)校正(SCT)[20,21]近似, 計算了每條反應(yīng)通道的速率常數(shù). 小曲率隧道校正方法假設(shè)最可幾隧道效應(yīng)路徑是沿反應(yīng)路徑附近區(qū)域. 值得注意的是, 該反應(yīng)體系中包含了2個重原子之間的輕粒子的遷移, 在低溫區(qū)間大曲率隧道校正(LCT, Large Curvature Tunneling)可能更加重要. 然而, 由于LCT方法需要更加全面的勢能面信息, 僅僅獲得最小能量路徑的信息是不夠的, 因此受計算條件的限制, 目前還不能完成大曲率的計算; 另外, 對于類似的氫提取反應(yīng), 已有理論研究[22～28]表明利用SCT方法得到的反應(yīng)速率常數(shù)值與實驗結(jié)果相符. 因此, 本文對該反應(yīng)體系采用SCT方法進(jìn)行速率常數(shù)的計算. 在計算OH自由基的電子配分函數(shù)時, 考慮了由自旋軌道耦合引起的能級劈裂, 劈裂能級為140 cm-1. 所有動力學(xué)計算均采用Polyrate 9.7程序包[29]完成.

2結(jié)果與討論

2.1穩(wěn)定點的性質(zhì)

由于反應(yīng)物CF2ClC(O)OCH2CH3存在2個可以產(chǎn)生不同結(jié)構(gòu)的轉(zhuǎn)動, 通過對分子內(nèi)基團(tuán)的相對扭轉(zhuǎn), 可以得到不同構(gòu)象. 以具有Cs對稱性的優(yōu)化結(jié)構(gòu)RC1(見圖1)為初始構(gòu)象, 在M06-2X/aug-cc-pVDZ水平上得到了隨二面角Cl—C—C—O和C—O—C—C變化的掃描勢能曲線[圖2(A), (B)]. 對掃描勢能曲線上對應(yīng)的極小點結(jié)構(gòu)進(jìn)行全優(yōu)化后, 得到穩(wěn)定構(gòu)象RC2和RC3(圖1). 從圖2可見, 這2個過程的轉(zhuǎn)換能壘分別約為0.96和5.15 kJ/mol. 為了確保得到所有可能的CF2ClC(O)OCH2CH3分子的穩(wěn)定構(gòu)象, 本文以構(gòu)象RC3為基礎(chǔ), 得到了扭轉(zhuǎn)二面角Cl—C—C—O的掃描勢能曲線[圖2(C)]. 從圖2(C)可見, 勢能曲線上還存在2個能量極小點, 經(jīng)全優(yōu)化后得到反應(yīng)物構(gòu)象RC4和RC5(圖1). 在M06-2X/aug-cc-pVDZ水平下, RC3與RC4和RC5之間的旋轉(zhuǎn)能壘分別約為0.46和0.71 kJ/mol. 在MCG3-MPWB//M06-2X/aug-cc-pVDZ水平下的計算結(jié)果表明, 反應(yīng)物的5個構(gòu)象中, RC2, RC4和RC5的能量接近, 其中RC4最穩(wěn)定, RC2和RC5的能量僅比RC4分別高0.79和0.38 kJ/mol, 而RC1和RC3的能量比RC4分別高5.31和6.07 kJ/mol.

Fig.1　Optimized geometries parameters of conformers RC1—RC5 at the M06-2X/aug-cc-pVDZ level　Bond lengths are in nm and bond angles are in degree.

Fig.2　Torsional potential profiles for torsion angle Cl—C—C—O for RC1(A), torsion angle 　C—O—C—C for RC1(B) and torsion angle Cl—C—C—O for RC3(C)

在與OH自由基的反應(yīng)中, H原子可以從CF2ClC(O)OCH2CH3中的—CH2—或—CH3基團(tuán)上提取. 對于具有Cs對稱性的構(gòu)象RC1, 存在平面內(nèi)氫(R1a)和平面外氫(R1b)兩類氫原子, 其中—CH3和—CH2—基團(tuán)中的兩對面外H原子是等價的, 對反應(yīng)的貢獻(xiàn)以對稱數(shù)的形式反映在總反應(yīng)速率常數(shù)中, 因此, RC1有3條可區(qū)分的氫提取反應(yīng)通道:

(R1a)

(R1b)

(R1c)

對于具有C1對稱性的構(gòu)象RC2—RC5, 由于—CH3和—CH2—基團(tuán)上的全部5個H原子均不等價, 得到與之對應(yīng)5條不同的氫提取反應(yīng)通道:

(Rna, Rnb, Rnc)

(Rnd, Rne)

因此, 對于反應(yīng)物CF2ClC(O)OCH2CH3的5個穩(wěn)定構(gòu)象, 我們共找到23條氫提取通道. 本文僅以能量最低的構(gòu)象RC4為例, 圖3給出了M06-2X/aug-cc-pVDZ水平上5條反應(yīng)通道的過渡態(tài)及中間體的幾何結(jié)構(gòu)及相應(yīng)的OH和H2O的實驗值, 其它通道上駐點的結(jié)構(gòu)及頻率見圖S1(見本文支持信息). 從圖3可見, 過渡態(tài)構(gòu)型中形成的鍵的伸長遠(yuǎn)大于相應(yīng)斷裂的鍵, 意味著過渡態(tài)都是和反應(yīng)物相似的, 即“早壘”型過渡態(tài). 如, 在TS4a中, 即將斷裂的C—H鍵長比單獨的反應(yīng)物CF2ClC(O)OCH2CH3中的C—H鍵長伸長了9%, 而即將形成的H—O鍵長比H2O中H—O鍵長伸長了39%. 反應(yīng)體系的其它氫提取過渡態(tài)具有與TS4a類似的結(jié)構(gòu)特性. 此外, 由于反應(yīng)物含有電負(fù)性較高的O原子和F原子, 因此在每條氫提取反應(yīng)通道的入口和出口處都能找到氫鍵絡(luò)合物, 這些絡(luò)合物結(jié)構(gòu)中, O…H和H…F的氫鍵距離均小于其范德華半徑之和(0.272和0.267 nm), 其它鍵長則均很接近其在相應(yīng)反應(yīng)物或產(chǎn)物結(jié)構(gòu)中的鍵長, 這表明CF3C(O)OCH2CH3+OH反應(yīng)可能要經(jīng)過一個間接的機(jī)理. 值得注意的是, 對于從—CH3和—CH2—基團(tuán)上的氫提取反應(yīng)通道, 由于OH自由基進(jìn)攻H原子的位置和方向的不同, 所得到的產(chǎn)物具有不同的構(gòu)象. 如, 從氫提取反應(yīng)通道(R4d)和(R4e)分別得到了P4b和P4c 2種不同結(jié)構(gòu)的產(chǎn)物自由基(見圖3).

Fig.3　Optimized geometry parameters of reactants, products, transition-states, and hydrogen-bond 　complexes for the H-abstraction channels of R4 at the M06-2X/aug-cc-pVDZ　Bond lengths are in nm and angles are in degree.

2.2能量計算

精確的標(biāo)準(zhǔn)生成焓在評估反應(yīng)熱動力學(xué)性質(zhì)和建立大氣模型過程中十分重要. 由于在實驗上缺少反應(yīng)物和產(chǎn)物的標(biāo)準(zhǔn)生成焓值(ΔH{0—f,298), 因此, 使用等化學(xué)鍵反應(yīng)(Ⅰ)～(Ⅲ)預(yù)測CF2ClC(O)OCH2CH3, CF2ClC(O)OCH2CH2和CF2ClC(O)OCHCH3的各構(gòu)象的ΔH{0—f,298值[30,31].

(Ⅰ)

(Ⅱ)

(Ⅲ)

在MCG3-MPWB//M06-2X/aug-cc-pVDZ水平上得到了反應(yīng)(Ⅰ)～(Ⅲ)中所有物種的幾何構(gòu)型、頻率和能量, 然后計算出各反應(yīng)的反應(yīng)焓值, 再結(jié)合已知參考化合物的標(biāo)準(zhǔn)生成焓值[32](CH4: -74.89, CHF2Cl: -481.37, CH3F: -234.18, CF3Cl: -709.65, CH3: 145.65和CH3C(O)OCH2CH3: -445.22 kJ/mol)來計算反應(yīng)物和產(chǎn)物自由基的各構(gòu)象對應(yīng)的標(biāo)準(zhǔn)生成焓值, 結(jié)果列于表1. 從表1可見, 對于構(gòu)象RC1, 以及相應(yīng)的產(chǎn)物自由基P1a, P1b和P1c, 其估算的標(biāo)準(zhǔn)生成焓分別為-838.47,-625.09, -628.02和-649.78 kJ/mol. 根據(jù)其余構(gòu)象及其產(chǎn)物自由基與RC1及產(chǎn)物自由基之間焓值的差值, 很容易得到所有反應(yīng)物構(gòu)象和產(chǎn)物自由基的標(biāo)準(zhǔn)生成焓值.

Table 1　Standard enthalpies of formation(ΔH{0—f,298) at 298 K calculated at

圖4為在MCG3-MPWB//M06-2X水平上反應(yīng)CF3C(O)OCH2CH3+ OH的勢能面示意圖, 其中能量包含了零點能校正. 以(R4)反應(yīng)的5條反應(yīng)通道為例, 可以看到,各反應(yīng)通道均為放熱反應(yīng), 其中CH3上氫提取通道(R4a)～(R4c)的反應(yīng)能為66.94 kJ/mol, 而—CH2—上氫提取通道由于產(chǎn)物構(gòu)象不同反應(yīng)能分別為-76.15和-88.70 kJ/mol. 反應(yīng)首先經(jīng)過反應(yīng)前期的氫鍵絡(luò)合物CR4a, CR4b, CR4c和CR4e, 相應(yīng)的能量分別為-9.62, -5.86, -14.64和-7.53 kJ/mol, 其中, (R4c)和(R4d) 2條反應(yīng)通道經(jīng)由相同的反應(yīng)前期絡(luò)合物CR4c. 然后, 從絡(luò)合物出發(fā), 經(jīng)過氫轉(zhuǎn)移過渡態(tài)生成產(chǎn)物絡(luò)合物, 進(jìn)而分解得到產(chǎn)物自由基CF2ClC(O)OCH2CH2(P4a)或CF2ClC(O)OCHCH3(P4b和P4c)和H2O. 絡(luò)合物CP4a, CP4b, CP4c, CP4d和CP4e分別比相應(yīng)的獨立產(chǎn)物穩(wěn)定11.72, 7.95, 10.46, 11.72和5.86 kJ/mol. 從圖4(D)可見, —CH3上的3個氫提取通道(R4a)～(R4c)和—CH2—上的2個氫提取通道(R4d)和(R4e)的反應(yīng)能壘分別為9.20, 14.64, 6.28, 12.13和8.37 kJ/mol, 其中通道(R4c)能壘最低. 因此, 僅從能量角度分析, 預(yù)測對于構(gòu)象RC4發(fā)生在—CH3基團(tuán)上的氫提取反應(yīng)更占優(yōu)勢, 是動力學(xué)上最可行的反應(yīng)通道. 從圖4(C)和(E)可見, 構(gòu)象RC3和RC5具有類似的反應(yīng)位置選擇性. 與反應(yīng)通道(R4)不同, 在圖4(A)和(D)給出的構(gòu)象RC1和RC2的各反應(yīng)通道中, —CH2—上的氫提取通道在動力學(xué)和熱力學(xué)上是最可行的, 是主要的反應(yīng)通道. 如, 通道(R1c)的反應(yīng)能壘為5.02 kJ/mol, 比—CH3上的氫提取通道(R1a)和(R1b)分別低了9.62和3.35 kJ/mol, 因此可以預(yù)測氫提取反應(yīng)主要發(fā)生在構(gòu)象RC1(或RC2)的—CH2—基團(tuán)上.

Fig.4　Schematic potential energy profiles for the RC1+OH reaction(A), the RC2+OH reaction(B), the RC3+OH reaction(C), the RC4+OH reaction(D) and the RC5 + OH reaction(E)Relative energies with ZPE atthe MCG3-MPWB//M06-2X/aug-cc-pVDZ level are in kJ/mol.

2.3速率常數(shù)計算

在MCG3-MPWB//M06-2X/aug-cc-pVDZ水平上, 使用改進(jìn)的變分過渡態(tài)理論結(jié)合小曲率隧道效應(yīng)校正(ICVT/SCT)計算了全部23條氫提取反應(yīng)通道在200～1000 K溫度范圍內(nèi)的速率常數(shù). 圖5給出了RC4的主要反應(yīng)通道(R4c)和(R4e)的TST, ICVT和ICVT/SCT速率常數(shù)曲線而其它穩(wěn)定體反應(yīng)通道的TST, ICVT和ICVT/SCT速率常數(shù)見圖S2(見本文支持信息). 從圖5(A)可見, 對于通道(R4c), 在298和500 K時,kICVT/kTST分別為0.92和0.93, 意味著變分效應(yīng)在整個溫度區(qū)間內(nèi)起著一定作用, 且反應(yīng)通道(R2c), (R3b), (R3c), (R3e), (R4a), (R4b), (R4d)和(R5c)也具有同樣的變分效果. 而對于反應(yīng)通道(R4e), 從圖5(B)可見, 在整個溫度區(qū)間范圍內(nèi), TST曲線和ICVT曲線幾乎完全重合, 即變分效應(yīng)很小或可以被忽略. 通道(R1a), (R2a), (R2e), (R5b), (R5d)和(R5e)具有類似的變分效應(yīng). 另一方面, 從圖5可見, 低溫時ICVT/SCT速率常數(shù)遠(yuǎn)大于ICVT速率常數(shù), 如, 反應(yīng)通道(R4c)的kICVT/SCT/kICVT值在200, 500和1000 K時分別為8.9, 2.1和1.03, 這表明小曲率隧道校正在低溫區(qū)間對這2個通道的速率常數(shù)計算中都起著重要作用. 其它各條氫提取反應(yīng)通道都具有類似的隧道效應(yīng), 具體速率常數(shù)結(jié)果列于表S2(見本文支持信息).

Fig.5　TST, ICVT and ICVT/SCT rate constants calculated at the MCG3-MPWB//M06-2X/aug-cc-pVDZ level 　versus 1000/T between 200 and 1000 K for path (R4c)(A) and for path (R4e)(B)

對每一穩(wěn)定構(gòu)象(RC1～RC5) 的總反應(yīng)速率常數(shù)為與其所相應(yīng)的氫提取反應(yīng)通道速率常數(shù)之和, 即

k1=(k1a+k1b)+k1c=k1(—CH3)+k1(—CH2—)

kn=(kna+knb+knc)+(knd+kne)=kn(—CH3)+kn(—CH2—),n=2—5

Fig.6　Plots of the calculated branching ratios versus 1000/T between 200 and 1000 K for the reaction RC1+OH→products(A) and the reaction RC4+OH→products(B)

圖6(A)和(B)分別給出了RC1和RC4的各氫提取反應(yīng)通道的產(chǎn)物分支比與溫度的關(guān)系, 而其余穩(wěn)定體的各氫提取反應(yīng)通道的產(chǎn)物分支比與溫度的關(guān)系在圖S2(見本文支持信息)給出. 從圖6(A)可見, 對于構(gòu)象RC1, 在200 K時, 其所對應(yīng)的的3條氫提取通道的產(chǎn)物分支比分別為0.08, 0.12和0.80, 隨著溫度的升高,k1c/k1比值下降, 而k1a/k1和k1b/k1比值上升, 當(dāng)達(dá)到500 K時三者的比值分別為0.31, 0.29和0.40. 表明通道(R1c)在低溫區(qū)間對構(gòu)象RC1做主要貢獻(xiàn), 而隨著溫度的增加通道(R1a)和(R1b)的作用越來越大. 考慮到提取—CH3上的3個氫原子在實驗上的不可分辨性, 本文給出—CH3基團(tuán)上氫提取對構(gòu)象RC1的總貢獻(xiàn)k1(—CH3)/k1, 其中k1(—CH3)=k1a+k1b. 顯然, 在低溫時, 從—CH2—基團(tuán)上的氫提取是反應(yīng)的主要通道, 而隨著溫度的升高, 從—CH3基團(tuán)上發(fā)生的氫提取通道對構(gòu)象RC1貢獻(xiàn)越來越大, 當(dāng)T>380 K時成為主要的反應(yīng)通道. 構(gòu)象RC2具有與RC1類似的氫提取位置選擇性. 圖6(B)給出的RC4的各通道速率常數(shù)分支比隨溫度變化的關(guān)系與RC1不同. 在整個溫度區(qū)間內(nèi), 發(fā)生在—CH3基團(tuán)上的氫提取反應(yīng)通道對構(gòu)象RC4的貢獻(xiàn), 即k4(—CH3)=k4a+k4b+k4c占絕對優(yōu)勢, 而發(fā)生在—CH2—基團(tuán)上的氫提取反應(yīng)通道[kn(—CH2—)=k4d+k4e]對總反應(yīng)的貢獻(xiàn)較小. 構(gòu)象RC3和RC5具有類似的反應(yīng)位置選擇性, 即—CH3基團(tuán)上的氫提取是反應(yīng)的主要通道.

反應(yīng)CF3C(O)OCH2CH3+OH→產(chǎn)物的總包反應(yīng)速率常數(shù)(koverall)由下式獲得:

koverall=ω1k1+ω2k2+ω3k3+ω4k4+ω5k5

=ω1(k1a+k1b+k1c)+ω2(k2a+k2b+k2c+k2d+k2e)+ω3(k3a+k3b+k3c+k3d+k3e)+

ω4(k4a+k4b+k4c+k4d+k4e)+ω5(k5a+k5b+k5c+k5d+k5e)

=k1T+k2T+k3T+k4T+k5T

式中:ωn(n=1～5)為利用Boltzmann配分函數(shù)計算得到的所有構(gòu)象對總包反應(yīng)貢獻(xiàn)的權(quán)重因子, 其表達(dá)式如下:

式中:R為理想氣體常數(shù), ΔEn為各構(gòu)象RCn(n=1～5)與構(gòu)象RC4的能量差. 表2給出了理論計算得出的速率常數(shù)和相應(yīng)實驗值. 在298 K時, 計算的雙水平ICVT/SCT速率值為9.12×10-13cm3·molecule-1·s-1, 與實驗值(5.4±1.5)×10-13cm3·molecule-1·s-1[3]相符. 由于實驗未提供298 K以外其它溫度的反應(yīng)速率常數(shù), 因此為了給實驗研究提供參考, 本文擬合了koverall在200～1000 K溫度范圍內(nèi)的三參數(shù)速率常數(shù)表達(dá)式, 即koverall= 5.45×10-25T4.54exp(-685/T).

Table 2Rate constants for each conformer and the overall rate constant[koverall/(cm3·molecule-1·s-1)] in the temperature range 200—1000 K at the MCG3-MPWB//M06-2X level

T/K1012k11012k21012k31012k41012k51012koverall2000.3550.4640.4710.5510.3600.4642500.4250.7560.6140.6950.5160.6452960.5291.120.8030.8770.7960.8982980.5361.140.8130.8880.8120.912(0.54±0.15)[3]3500.7051.671.110.1611.401.344000.9262.311.500.1492.351.925001.564.062.640.2366.123.806002.536.514.403.6013.77.038005.8814.110.47.3347.519.4100011.826.120.713.212143.5

Fig.7　Plot of the contributions of RC1, RC2, RC3, RC3 and RC4 to the reaction versus 1000/T between 200 and 1000 K

為了描述各構(gòu)象對總反應(yīng)的貢獻(xiàn), 圖7給出了knT/koverall隨溫度變化的曲線. 從圖7可見, 在200 K時k1T/koverall,k2T/koverall,k3T/koverall,k4T/koverall和k5T/koverall的比值分別為0.01, 0.31, 0.01, 0.45和0.22, 表明在低溫區(qū)間構(gòu)象RC4對總反應(yīng)的貢獻(xiàn)最大, 是最重要的反應(yīng)物構(gòu)象, 其次是構(gòu)象RC2和RC5, 而構(gòu)象RC1和RC3對總包反應(yīng)的貢獻(xiàn)可以忽略不計; 而隨著溫度的升高,k4T/koverall和k2T/koverall的比值下降,k5T/koverall的比值則升高, 特別是在400 K及更高的溫度范圍內(nèi),k5T/koverall的比值急速上升. 當(dāng)T>500 K時, 各構(gòu)象對總反應(yīng)速率的貢獻(xiàn)順序為k5T/koverall>k2T/koverall>k4T/koverall>k3T/koverall>k1T/koverall, 表明在高溫區(qū)間構(gòu)象RC5對總反應(yīng)速率具有決定性的貢獻(xiàn).

綜上所述, 采用直接動力學(xué)方法, 研究了CF2ClC(O)OCH2CH3與OH自由基的反應(yīng)的微觀反應(yīng)機(jī)理和動力學(xué)性質(zhì). 在M06-2X/aug-cc-pVDZ水平上確定了反應(yīng)物的5個穩(wěn)定的幾何構(gòu)象(RC1～RC5), 其中構(gòu)象RC4為最穩(wěn)定構(gòu)象, RC2和RC5的能量接近. 由于這5個構(gòu)象之間的能量差較小, 所以在動力學(xué)計算中它們對總反應(yīng)速率的貢獻(xiàn)應(yīng)考慮在內(nèi). 對于每一個構(gòu)象, 分別研究了發(fā)生在—CH2—和—CH3基團(tuán)上的氫提取反應(yīng)通道, 得到23條氫提取反應(yīng)通道. 此外, 由于缺少關(guān)于反應(yīng)物及產(chǎn)物自由基的標(biāo)準(zhǔn)生成焓數(shù)據(jù), 用等化學(xué)鍵法估算了反應(yīng)物及產(chǎn)物自由基的標(biāo)準(zhǔn)生成焓. 使用改進(jìn)的變分過渡態(tài)理論結(jié)合小曲率隧道效應(yīng)校正計算了各個氫提取通道在200～1000 K內(nèi)的速率常數(shù), 計算結(jié)果表明, 對于構(gòu)象RC1和RC2, 氫提取反應(yīng)在低溫區(qū)間內(nèi)主要發(fā)生在—CH2—基團(tuán)上, 隨著反應(yīng)溫度的升高, 發(fā)生在—CH3基團(tuán)上的氫提取反應(yīng)通道成為主要反應(yīng)通道. 而對于構(gòu)象RC3, RC4和RC5, 在整個溫度區(qū)間內(nèi), 發(fā)生在—CH3基團(tuán)上的氫提取反應(yīng)通道都占絕對優(yōu)勢. 根據(jù)Boltzmann配分函數(shù)計算200～1000 K溫度區(qū)間內(nèi)總包反應(yīng)的速率常koverall, 298 K的計算結(jié)果與實驗值相吻合, 進(jìn)而擬合了總反應(yīng)在較大溫度區(qū)間內(nèi)的三參數(shù)的Arrhenius表達(dá)式:koverall=5.45×10-25T4.54exp(-685/T). 另外, 分析了5種構(gòu)象對總包反應(yīng)的貢獻(xiàn), 表明在低溫區(qū)間構(gòu)象RC4對總反應(yīng)的貢獻(xiàn)最大; 其次為RC2和RC5; 隨著溫度升高, 構(gòu)象RC5的貢獻(xiàn)變得越來越重要, 在高溫區(qū)間對總反應(yīng)貢獻(xiàn)最大.

支持信息見http://www.cjcu.jlu.edu.cn/CN/10.7503/cjcu20150568.

參考文獻(xiàn)

[1]Oyaro N., Sellevag S. R., Nielsen C. J.,Environ.Sci.Technol., 2004, 38, 5567—5576

[2]Sulback Andersen M. P., Nielsen O. J., Wallington T. J., Hurley M. D., DeMoore G. W.,J.Phys.Chem.A, 2005, 109, 3926—3934

[3]Blanco M. B., Teruel M. A.,Chem.Phys.Lett., 2007, 441, 1—6

[4]Hu W. P., Truhlar D. G.,J.Am.Chem.Soc., 1996, 118, 860—869

[5]Truhlar D. G., Garrett B. C., Klippenstein S. J.,J.Phys.Chem., 1996, 100, 12771—12800

[6]Truhlar D. G., Gordon M. S.,Science, 1990, 249, 491—498

[7]Tucker S. C., Truhlar D. G.,NewTheoreticalConceptsforUnderstandingOrganicReaction, Kluwer, Dordrecht, 1989, 291—346

[8]Truhlar D. G., Gordon M. S., Stechler R.,Chem.Rev., 1987, 87, 217—236

[9]Lu D. H., Truong T. N., Melissas V. S., Lynch G. C., Liu Y. P., Grarrett B. C., Stechler R., Issacson A. D., Rai S. N., Hancock G. C., Lauderdale J. G., Joseph T., Truhlar D. G.,Comput.Phys.Commum., 1992, 71, 235—262

[10]Garrett B. C., Truhlar D. G.,J.Phys.Chem., 1991, 95, 10374—10379

[11]Truhlar D. G., Garrett B. C.,Acc.Chem.Res., 1980, 13, 440—448

[12]Truhlar D. G., Garrett B. C.,Annu.Rev.Phys.Chem., 1984, 35, 159—189

[13]Frisch M. J., Trucks G. W., Schlegel H. B., Scuseria G. E., Robb M. A., Cheeseman J. R., Scalmani G., Barone V., Mennucci B., Petersson G. A., Nakatsuji H., Caricato M., Li X., Hratchian H. P., Izmaylov A.F., Bloino J., Zheng G., Sonnenberg J. L., Hada M., Ehara M., Toyota K., Fukuda R., Hasegawa J., Ishida M., Nakajima T., Honda Y., Kitao O., Nakai H., Vreven T., Montgomery J. A. Jr., Peralta J. E., Ogliaro F., Bearpark M., Heyd J. J., Brothers E., Kudin K. N., Staroverov V. N., Kobayashi R., Normand J., Raghavachari K., Rendell A., Burant J. C., Iyengar S. S., Tomasi J., Cossi M., Rega N., Millam N. J., Klene M., Knox J. E., Cross J. B., Bakken V., Adamo C., Jaramillo J., Gomperts R., Stratmann R. E., Yazyev O., Austin A. J., Cammi R., Pomelli C., Ochterski J. W., Martin R. L., Morokuma K., Zakrzewski V. G., Voth G. A., Salvador P., Dannenberg J. J., Dapprich S., Daniels A. D., Farkas O., Foresman J. B., Ortiz J. V., Cioslowski J., Fox D. J.,Gaussian09,RevisionA.1, Gaussian Inc., Wallingford CT, 2009

[14]Zhao Y., Truhlar D.G.,Theor.Chem.Acc., 2008, 120, 215—241

[15]Zhao Y., Truhlar D. G., J.Phys.Chem.A, 2005, 109, 4209—4212

[16]Frisch M. J., Trucks G. W., Schlegel H. B., Scuseria G. E., Robb M. A., Cheeseman J. R., Zakrzewski V. G., Montgomery J. A. Jr., Stratmann R. E., Burant J. C., Dapprich S., Millam J. M., Daniels A. D., Kudin K. N., Strain M. C., Farkas O., Tomasi J., Barone V., Cossi M., Cammi R., Mennucci B., Pomelli C., Adamo C., Clifford S., Ochterski J., Petersson G. A., Ayala P. Y., Cui Q., Morokuma K., Malick D. K., Rabuck A. D., Raghavachari K., Foresman J. B., Cioslowski J., Ortiz J. V., Boboul A. G., Stefnov B. B., Liu G., Liaschenko A., Piskorz P., Komaromi L., Gomperts R., Martin R. L., Fox D. J., Keith T., Al-Laham M. A., Peng C. Y., Nanayakkara A., Gonzalez C., Challacombe M., Gill P. M. W., Johnson B., Chen W., Wong M. W., Andres J. L., Gonzalez C., Head-Gordon M., Replogle E. S., Pople J. A,Gaussian03,RevisionA.1, Guassian Inc., Pittsburgh PA, 2003

[17]Zhao Y., Truhlar D. G.,MLGAUSS-Version2.0, University of Minnesota, Minneapolis, 2005

[18]Chuang Y. Y., Corchado J. C., Truhlar D. G.,J.Phys.Chem.A, 1999, 103, 1140—1149

[19]Garrett B. C., Truhlar D. G., Grev R. S., Magnuson A. W.,J.Phys.Chem., 1980, 84, 1730—1748

[20]Truhlar D. G.,J.Comput.Chem., 1991, 12, 266—270

[21]Liu Y. P., Lynch G. C., Truong T. N., Lu D. H., Truhlar D. G., Garrett B. C.,J.Am.Chem.Soc., 1993, 115, 2408—2415

[22]Galano A., Alvarez-Idaboy J. R., Ruiz-Santoyo M. E., Vivier-Bunge A.,Chem.Phys.Chem, 2004, 5, 1379—1388

[23]Taghikhani M., Parsafar G. A., Sabzyan H.,J.Phys.Chem.A, 2005, 109, 8158—8167

[24]Wang Y. X., Gao H., Wang Q., Liu J. Y.,Chem.J.ChineseUniversities, 2010, 31(6), 1240—1245(王永霞, 高紅, 王欽, 劉靖堯. 高等學(xué)校化學(xué)學(xué)報, 2010, 31(6), 1240—1245)

[25]Ren H., Zhang L. L., Wang R. S., Pan X. M.,J.Phys.Chem.A, 2012, 116, 10647—10655

[26]Jin T. Y., Wang Q., Liu J. Y.,Chem.J.ChineseUniversities, 2013, 34(3), 641—649(金銅音, 王欽, 劉靖堯. 高等學(xué)?；瘜W(xué)學(xué)報, 2013, 34(3), 641—649)

[27]Zhu P., Ai L. L., Wang H., Liu J. Y.,Comput.Theor.Chem., 2014, 1029, 91—98

[28]Zhu P., Duan X. M., Liu J. Y.,J.Fluorine.Chem., 2015, 176, 61—70

[29]Corchado J. C., Chuang Y. Y., Fast P. L., Hu W. P., Liu Y. P., Lynch G. C., Nguyen K. A., Jackels C. F., Ramos A. F., Ellingson B. A., Lynch B. J., Melissas V. S., Villa J., Rossi I., Coitino E. L., Pu J., Albu T. V., Steckler R., Garrett B. C., Isaacson A. D., Truhlar D. G.,POLYRATE,Version9.7, University of Minnesota, Minneapolis, 2007

[30]Jin X. H., Hu B. C., Jia H. Q., Lu C. X.,Chem.J.ChineseUniversities, 2013, 34(7), 1685—1690(金興輝, 胡炳成, 賈歡慶, 呂春緒. 高等學(xué)?；瘜W(xué)學(xué)報, 2013, 34(7), 1685—1690)

[31]Wang H. X., Wang B. Y., Zhang J. L., Li Z. R., Li X. Y.,Chem.J.ChineseUniversities, 2011, 32(5), 1123—1128(王海霞, 汪必耀, 張俊玲, 李澤榮, 李象遠(yuǎn). 高等學(xué)校化學(xué)學(xué)報, 2011, 32(5), 1123—1128)

[32]Linstrom P. J., Mallard W. G.,NISTChemistryWebbook, http://webbook.Nist.Gov/chemistry

Mechanism and Kinetics of the Hydrogen-abstraction Reaction of

CF2ClC(O)OCH2CH3with OH Radicals?

ZHU Peng, DUAN Xuemei, LIU Jingyao*

(InstituteofTheoreticalChemistry,JilinUniversity,Changchun130021,China)

AbstractThe mechanism of the hydrogen abstraction reaction of CF2Cl(O)OCH2CH3+OH was studied theoreticallyviaa dual-level direct dynamics method at the MCG3-MPWB//M06-2X/aug-cc-pVDZ level. Five stable conformers(RC1—RC5) of the reactant CF2Cl(O)OCH2CH3were located, and for each conformer, the possible H-abstraction channels from —CH3and —CH2— groups were taken into account. The rate constants were calculated using the improved canonical variational transition-state theory(ICVT) with the small-curvature tunneling correction(SCT) and the selectivity of the reaction sites of CF2Cl(O)OCH2CH3was evaluated. The results show that for the conformers RC1 and RC2, the H-abstraction reactions mainly take place at the —CH2— group at low temperature, while for the conformers RC3, RC4 and RC5, the hydrogen abstraction from the —CH3group is the major channel in the whole considered temperature range. The overall rate constant is obtained by considering the weight factors of the five conformers calculated from the Boltzmannn distribution function. It is found, that the calculatedkoverallat 298 K is in good agreement with the available experimental data. The three parameter expression for the overall reaction within 200—1000 K is fitted tokoverall=5.45×10-25T4.54exp(-685/T).

KeywordsDirect dynamics method; Rate constant; Variational transition-state theory; Density functional theory; Reaction mechanism

(Ed.: Y, Z)

? Supported by the National Natural Science Foundation of China(No.21373098).

doi:10.7503/cjcu20150522

基金項目：國家自然科學(xué)基金(批準(zhǔn)號: 21373098)資助.

收稿日期:2015-07-20. 網(wǎng)絡(luò)出版日期: 2015-12-20.

中圖分類號O641; O643

文獻(xiàn)標(biāo)志碼A

聯(lián)系人簡介: 劉靖堯, 女, 博士, 教授, 博士生導(dǎo)師, 主要從事化學(xué)反應(yīng)動力學(xué)研究. E-mail: ljy121@jlu.edu.cn