李尚俊 談寧馨 姚 倩 李澤榮,* 李象遠(yuǎn)

(1四川大學(xué)化學(xué)學(xué)院,成都610064;2四川大學(xué)化學(xué)工程學(xué)院,成都610065)

1 引言

目前,由于燃料的使用越來(lái)越廣泛,為了提高燃燒效應(yīng)、減少排放,燃燒反應(yīng)的研究越來(lái)越重要,大量的新實(shí)驗(yàn)和傳統(tǒng)的實(shí)驗(yàn)技術(shù)應(yīng)用于碳?xì)淙剂系娜紵^(guò)程研究,與此同時(shí),隨著量子化學(xué)和反應(yīng)速率理論的發(fā)展,在高溫條件下碳?xì)淙剂系娜紵M研究已經(jīng)取得了很大成就,1,2發(fā)展建立了很多單組分碳?xì)淙剂虾筒煌旌先剂?如柴油、煤油)的替代燃燒模型,能較好地再現(xiàn)實(shí)驗(yàn)結(jié)果.但在低溫條件下(T＜900 K),由于涉及反應(yīng)類型較多,其中絕大多數(shù)反應(yīng)類型缺乏反應(yīng)速率系數(shù)的實(shí)驗(yàn)數(shù)據(jù),因而在燃燒機(jī)理模型構(gòu)建中,通常按反應(yīng)類型規(guī)則給出近似速率系數(shù),即同一反應(yīng)類型中不同反應(yīng)采用相同速率系數(shù),或由基團(tuán)加和方法給出近似速率系數(shù),目前為止所建立的燃燒模型的模擬結(jié)果很多情形下不是太令人滿意.3現(xiàn)在人們普遍認(rèn)為,在碳?xì)淙剂系母邷厝紵偷蜏厝紵龣C(jī)理中,烷基自由基與氧氣分子生成烷基過(guò)氧自由基(R+O2=ROO)的反應(yīng)是最關(guān)鍵反應(yīng),該反應(yīng)由于是強(qiáng)放熱反應(yīng),其平衡常數(shù)強(qiáng)烈地依賴于溫度,低溫下有利于ROO的生成,而高溫下有利于反應(yīng)R+O2→R′+HO2的進(jìn)行,因而該反應(yīng)及后續(xù)反應(yīng)在高溫下不考慮,而在低溫下起關(guān)鍵作用.4,5在其后續(xù)反應(yīng)中,烷基過(guò)氧自由基ROO經(jīng)分子內(nèi)氫遷移形成氫過(guò)氧烷基自由基(QOOH)的反應(yīng)為過(guò)氧烷基自由基的最主要消耗反應(yīng)通道,6,7是低溫燃燒反應(yīng)中的一類重要基元反應(yīng),因而其精確的熱力學(xué)和動(dòng)力學(xué)參數(shù)是碳?xì)浠衔锏蜏厝紵M的基礎(chǔ).該類反應(yīng)的過(guò)渡態(tài)為四元環(huán)、五元環(huán)、六元環(huán)等環(huán)狀結(jié)構(gòu).Zádor等5給出了低溫下碳?xì)浠衔锶紵婕暗年P(guān)鍵基元反應(yīng)類型,并報(bào)道了從C2到C5的烷基過(guò)氧自由基的分子內(nèi)氫遷移反應(yīng)的能壘計(jì)算的文獻(xiàn)綜述.

對(duì)該類反應(yīng)中個(gè)別反應(yīng)的計(jì)算研究,文獻(xiàn)已有大量的報(bào)道.Miller等6在近G2水平計(jì)算得到了乙基過(guò)氧自由基(CH3CH2O2)的分子內(nèi)氫遷移反應(yīng)的能壘;Taatjes3對(duì)正丙基過(guò)氧自由基(CH3CH2CH2O2)分子內(nèi)氫遷移反應(yīng)進(jìn)行了實(shí)驗(yàn)和量子化學(xué)計(jì)算研究,得到了近QCISD(T)水平的過(guò)渡態(tài)的能量,發(fā)現(xiàn)分子內(nèi)氫遷移經(jīng)六元環(huán)過(guò)渡態(tài)的能壘比經(jīng)五元環(huán)過(guò)渡態(tài)的能壘更低.Villano等7,8在CBS-QB3水平計(jì)算得到了C1-C5的烷基過(guò)氧自由基的分子內(nèi)氫遷移反應(yīng)及逆反應(yīng)的速率系數(shù);Zhang等9在CBS-QB3水平計(jì)算得到了一些C3-C5烷基過(guò)氧自由基和烯基過(guò)氧自由基分子內(nèi)氫遷移反應(yīng)的能壘.由于在燃燒機(jī)理的構(gòu)建與機(jī)理自動(dòng)生成程序中,通常將反應(yīng)按其相似性將其分為不同反應(yīng)類,如對(duì)烷烴RH的奪氫反應(yīng)根據(jù)奪氫基團(tuán)分為:RH+O=R+OH反應(yīng)類、RH+CHO=R+HCHO反應(yīng)類、RH+H=R+H2反應(yīng)類等等.且為方便起見(jiàn),在機(jī)理自動(dòng)生成程序中通常同一類分子反應(yīng)的速率常數(shù)近似地采用同一數(shù)值.反應(yīng)類的選取和劃分,即所謂的反應(yīng)規(guī)則,在燃燒機(jī)理的構(gòu)建與機(jī)理自動(dòng)生成程序中就很關(guān)鍵.但由于相似性沒(méi)有嚴(yán)格的定義或標(biāo)準(zhǔn),因而反應(yīng)類的選取和劃分往往是根據(jù)化學(xué)直覺(jué)或根據(jù)結(jié)果的好壞.對(duì)我們研究的烷基過(guò)氧自由基分子內(nèi)氫遷移反應(yīng),人們發(fā)現(xiàn)不能簡(jiǎn)化為一類.Sharma等10按氫遷移距離及遷移氫所在碳原子類型,把烷基過(guò)氧自由基分子內(nèi)氫遷移反應(yīng)進(jìn)一步劃分為許多類型,在CBS-QB3水平計(jì)算得到了從(1,3)氫遷移到(1,8)氫遷移各反應(yīng)類速率系數(shù)規(guī)則;Miyoshi11按氫遷移距離及遷移氫所在碳原子類型和過(guò)氧基所在碳原子類型,把該反應(yīng)類進(jìn)一步劃分為許多子類型,在CBS-QB3水平計(jì)算得到了從(1,4)氫遷移到(1,8)氫遷移各反應(yīng)類及逆反應(yīng)各反應(yīng)類速率系數(shù)規(guī)則.因此,為了精確描述烷基過(guò)氧自由基分子內(nèi)氫遷移反應(yīng)的反應(yīng)規(guī)則,需要將其細(xì)分為很多類型.即便如此,研究發(fā)現(xiàn)在同一類中不同反應(yīng)的速率常數(shù)相差較大,因而這種反應(yīng)規(guī)則的方法對(duì)烷基過(guò)氧自由基分子內(nèi)氫遷移反應(yīng)有較大局限性.為了構(gòu)建碳?xì)淙剂峡煽康牡蜏厝紵龣C(jī)理,最可靠的辦法是對(duì)烷基過(guò)氧自由基分子內(nèi)氫遷移反應(yīng)涉及的每一反應(yīng)分別計(jì)算出精確的速率常數(shù).但是,目前的量子化學(xué)高精度計(jì)算方法只能局限于較少碳原子數(shù)的烷基過(guò)氧自由基的分子體系.目前為止,文獻(xiàn)報(bào)道的體系為C8及以下體系.本文的目的是通過(guò)將等鍵反應(yīng)理論與過(guò)渡態(tài)理論相結(jié)合,計(jì)算得到C8以上烷基過(guò)氧自由基的分子體系分子內(nèi)氫遷移反應(yīng)的精確速率常數(shù).

2 計(jì)算方法

單分子反應(yīng)過(guò)渡態(tài)理論的反應(yīng)速率常數(shù)表達(dá)式為:

式中,κ為隧穿系數(shù);R為理想氣體常數(shù);T為溫度;σ是反應(yīng)對(duì)稱數(shù);kB為Boltzmann常數(shù);h為Planck常數(shù);Q≠為過(guò)渡態(tài)配分函數(shù);QA為反應(yīng)物A的配分函數(shù);ΔV≠為過(guò)渡態(tài)與反應(yīng)物A的包含零點(diǎn)振動(dòng)能的電子能之差,即反應(yīng)能壘.在單分子反應(yīng)過(guò)渡態(tài)理論速率常數(shù)的計(jì)算公式(1)中,振動(dòng)貢獻(xiàn)決定于過(guò)渡態(tài)和反應(yīng)物的諧振頻率,諧振頻率的計(jì)算可在較低從頭算級(jí)別計(jì)算.(1)式中關(guān)鍵項(xiàng)是勢(shì)壘項(xiàng),必須采用高精度計(jì)算方法才可得到精確數(shù)據(jù).對(duì)于小分子的基元反應(yīng)高精度的計(jì)算是可行的,但是,對(duì)于較大分子以及大分子的基元反應(yīng)高精度的量子化學(xué)計(jì)算幾乎是不可能的.

Truong等12-18以過(guò)渡態(tài)理論為基礎(chǔ),將反應(yīng)類的概念引入反應(yīng)類過(guò)渡態(tài)理論(RC-TST)中.該理論認(rèn)為同一類反應(yīng)中所有反應(yīng)的活性中心相同,因此,任意兩個(gè)反應(yīng)的速率常數(shù)的不同主要在于取代基的不同.他們通常將同一類反應(yīng)中最小反應(yīng)體系定義為參考反應(yīng),其余反應(yīng)成為目標(biāo)反應(yīng).Truong等研究發(fā)現(xiàn),參考反應(yīng)與目標(biāo)反應(yīng)的能壘差隨不同從頭算理論方法變化而變化.但是,Truong等提出的反應(yīng)類過(guò)渡態(tài)理論對(duì)反應(yīng)類并沒(méi)有給出嚴(yán)格定義,不能簡(jiǎn)單地用于烷基過(guò)氧自由基體系分子內(nèi)氫遷移反應(yīng)體系.

在我們的前期研究中,我們將等鍵反應(yīng)方法引入了反應(yīng)類的過(guò)渡態(tài)生成過(guò)程,19,20對(duì)反應(yīng)類過(guò)渡態(tài)理論提供了解釋,將等鍵反應(yīng)的定義推廣到非典型化學(xué)鍵的過(guò)渡態(tài)情形,并用于反應(yīng)類目標(biāo)反應(yīng)能壘及速率常數(shù)的精確計(jì)算.

目標(biāo)反應(yīng)能壘表示如下:

式中,P代表參考反應(yīng),T代表目標(biāo)反應(yīng),≠代表過(guò)渡態(tài),Δ?和 Δ分別為用高級(jí)別從頭算方法精確計(jì)算和低級(jí)別從頭算方法近似計(jì)算的參考反應(yīng)的反應(yīng)能壘,其差值ΔΔ為參考反應(yīng)的能壘修正值.Δ為低級(jí)別從頭算方法近似計(jì)算的目標(biāo)反應(yīng)的反應(yīng)能壘,而Δ?為對(duì)目標(biāo)反應(yīng)的反應(yīng)能壘精確計(jì)算的估計(jì).(2)、(3)式表明:目標(biāo)反應(yīng)的精確反應(yīng)勢(shì)壘可以由其低水平計(jì)算的近似結(jié)果通過(guò)參考反應(yīng)的能壘修正值修正得到.同時(shí),目標(biāo)反應(yīng)的速率常數(shù)kT可由低水平下的目標(biāo)反應(yīng)的速率常數(shù)kP及參考反應(yīng)的能壘修正值得到:

本文選取了62個(gè)過(guò)氧烷基分子內(nèi)氫遷移形成氫過(guò)氧化烷基自由基反應(yīng)作為研究對(duì)象,反應(yīng)見(jiàn)表S1(見(jiàn)Supporting Information).所有的量子化學(xué)計(jì)算都是在Gaussian 03程序21下進(jìn)行.所有反應(yīng)涉及到的物種的幾何結(jié)構(gòu)優(yōu)化、頻率計(jì)算和低級(jí)別單點(diǎn)計(jì)算均在B3LYP/6-311+G(d,p)22水平下進(jìn)行,并用內(nèi)稟反應(yīng)坐標(biāo)(IRC)理論對(duì)體系中所有反應(yīng)過(guò)渡態(tài)與對(duì)應(yīng)反應(yīng)物和產(chǎn)物的關(guān)聯(lián)性進(jìn)行了確認(rèn).本文中所有振動(dòng)均采用諧振子模型,沒(méi)有考慮非諧效應(yīng)對(duì)熱力學(xué)參數(shù)及動(dòng)力學(xué)參數(shù)的影響.熱力學(xué)配分函數(shù)頻率校正因子為0.9663,23高級(jí)別從頭算方法采用模型化學(xué)方法CBS-QB3.24之所以選擇CBS-QB3方法作為高精度的單點(diǎn)能計(jì)算,是因?yàn)锽ozzelli等25用此方法計(jì)算了一系列過(guò)氧烷基分子的熱力學(xué)參數(shù),得到了與實(shí)驗(yàn)接近的結(jié)果.

3 結(jié)果與討論

3.1 過(guò)渡態(tài)幾何中心結(jié)構(gòu)

過(guò)氧烷基分子內(nèi)氫遷移反應(yīng)的過(guò)渡態(tài)涉及碳上氫原子遷移到過(guò)氧自由基一端氧原子形成的環(huán)狀結(jié)構(gòu),(1,3)氫遷移為四元環(huán),(1,4)氫遷移為五元環(huán),(1,5)氫遷移為六元環(huán),(1,6)氫遷移為七元環(huán),等等.該類反應(yīng)中心所涉及到的原子,鍵長(zhǎng)與鍵角編號(hào)如圖 1所示,其中d1、d2、d3、d4為鍵長(zhǎng),A1、A2、A3為鍵角.

B3LYP/6-311+G(d,p)方法優(yōu)化得到的各具體反應(yīng)的反應(yīng)中心過(guò)渡態(tài)幾何結(jié)構(gòu)參數(shù)列于表S2(見(jiàn)Supporting Information).由表S2可以看出,環(huán)狀大小相同的過(guò)渡態(tài)結(jié)構(gòu)幾乎相同,所以本文對(duì)環(huán)狀大小相同的過(guò)渡態(tài)反應(yīng)中心的結(jié)構(gòu)參數(shù)進(jìn)行平均,其值列于表1.

圖1 過(guò)渡態(tài)(TS)反應(yīng)中心的幾何結(jié)構(gòu)Fig.1 Geometry of the reaction centers in the transition states(TS)

從表1可以看出,按照氫遷移距離,從(1,3)氫遷移到(1,4)氫遷移、(1,5)氫遷移、(1,6)氫遷移,其反應(yīng)中心過(guò)渡態(tài)幾何結(jié)構(gòu)參數(shù)變化明顯,但從(1,6)氫遷移到(1,7)氫遷移、(1,8)氫遷移其幾何結(jié)構(gòu)參數(shù)很接近.因此,本文根據(jù)過(guò)氧烷基分子內(nèi)氫遷移反應(yīng)的反應(yīng)中心過(guò)渡態(tài)幾何結(jié)構(gòu)參數(shù)是否守恒的特征,將(1,3)氫遷移、(1,4)氫遷移和(1,5)氫遷移各分為一類,而將(1,6)氫遷移、(1,7)氫遷移、(1,8)歸為一類反應(yīng),即為((1,n),n=6,7,8).

為了更進(jìn)一步反映這樣劃分的每一類反應(yīng)過(guò)渡態(tài)幾何結(jié)構(gòu)是否守恒,表2給出了其幾何結(jié)構(gòu)參數(shù)的最大偏差和平均偏差.從表2可以看出,各類反應(yīng)過(guò)渡態(tài)幾何結(jié)構(gòu)參數(shù)偏差很小,表明每一類反應(yīng)其過(guò)渡態(tài)幾何結(jié)構(gòu)是守恒的,因此對(duì)這樣劃分的每一類反應(yīng),可以使用等鍵反應(yīng)方法.本文將研究的62個(gè)反應(yīng)分為(1,3),(1,4),(1,5),((1,n),n=6,7,8)四類反應(yīng),分別選取每一類反應(yīng)的最簡(jiǎn)單反應(yīng),即表S1中R1、R11、R23、R57反應(yīng),作為參考反應(yīng),剩余58個(gè)反應(yīng)為目標(biāo)反應(yīng).

3.2 反應(yīng)勢(shì)壘

由于缺乏過(guò)氧烷基分子內(nèi)氫遷移反應(yīng)動(dòng)力學(xué)參數(shù)實(shí)驗(yàn)報(bào)道,本文選用CBS-QB3方法作為參考反應(yīng)的精確計(jì)算方法.同時(shí)在6-311+G(d,p)水平下,采用B3LYP方法作為參考反應(yīng)的近似計(jì)算方法.兩種方法計(jì)算得到的參考反應(yīng)的反應(yīng)勢(shì)壘和反應(yīng)焓變以及兩種方法的差值結(jié)果列于表3.

從表3可以看出,低級(jí)別的B3LYP計(jì)算得到的參考反應(yīng)勢(shì)壘與高級(jí)別的CBS-QB3結(jié)果比較,有較大偏差,其偏差范圍在2.70-11.06 kJ?mol-1.所以,低級(jí)別的B3LYP計(jì)算得到的結(jié)果有較大誤差,有必要對(duì)其修正.所以,本文對(duì)其余58個(gè)參考反應(yīng)采用等鍵反應(yīng)方法進(jìn)行了修正,修正前后的反應(yīng)勢(shì)壘見(jiàn)表S3.對(duì)較小分子體系的部分目標(biāo)反應(yīng),本文用高精度CBS-QB3方法進(jìn)行了精確計(jì)算,其結(jié)果與用等鍵反應(yīng)修正后反應(yīng)勢(shì)壘的比較列于表4.從表4可以看出,直接用密度泛函理論(DFT)方法計(jì)算的反應(yīng)勢(shì)壘的誤差范圍為-16.94--27.15 kJ?mol-1,遠(yuǎn)超化學(xué)精度,但經(jīng)等鍵反應(yīng)方法修正以后,其誤差范圍變?yōu)?1.96-2.80 kJ?mol-1,已在化學(xué)精度以內(nèi).因而,計(jì)算結(jié)果表明,過(guò)氧烷基分子內(nèi)氫遷移反應(yīng)用DFT方法經(jīng)過(guò)等鍵反應(yīng)方法修正可得到可靠的反應(yīng)勢(shì)壘.而且,從表S3可以看出,各遷移反應(yīng)類型勢(shì)壘與過(guò)渡態(tài)環(huán)狀大小有關(guān),呈現(xiàn)規(guī)律性變化,其變化規(guī)律為:(1,3)＞(1,4)＞(1,5)＞(1,6)＜(1,7)＜(1,8).

3.3 反應(yīng)速率常數(shù)

本文根據(jù)過(guò)渡態(tài)理論及等鍵反應(yīng)修正后勢(shì)壘,19,20使用chemrate程序26對(duì)所涉及的62個(gè)反應(yīng)計(jì)算了300-1500 K范圍的速率參數(shù),并按阿侖尼烏斯方程(見(jiàn)式(5)),經(jīng)最小二乘法擬合得到動(dòng)力學(xué)參數(shù)(A,Ea),結(jié)果列于表5中.

表1 反應(yīng)中心過(guò)渡態(tài)平均幾何結(jié)構(gòu)參數(shù)Table 1 Average geometrical parameters of the reaction centers for the transition states

表2 每一類反應(yīng)過(guò)渡態(tài)中心幾何結(jié)構(gòu)參數(shù)的最大絕對(duì)偏差(MAD)和平均絕對(duì)偏差(AAD)Table 2 Maximum absolute deviation(MAD)and average absolute deviation(AAD)of the geometrical parameters of the transition states for each class reaction

表3 參考反應(yīng)的反應(yīng)勢(shì)壘Table 3 Energy barriers for the reference reactions

表4 B3LYP方法修正前后勢(shì)壘與CBS-QB3結(jié)果比較Table 4 Comparison of B3LYPreaction barriers(Δ)before and after validation with the CBS-QB3 results

表4 B3LYP方法修正前后勢(shì)壘與CBS-QB3結(jié)果比較Table 4 Comparison of B3LYPreaction barriers(Δ)before and after validation with the CBS-QB3 results

acalculated at the B3LYP/6-311+G(d,p)level of theory,bdifference between density functional theory(DFT)barrier and CBS-QB3 barrier:Δ(DFT)=ΔV(B3LYP)-ΔV(CBS-QB3),cDFT results after validation by isodesmic reaction method(IRM),ddeviation between IRM value and CBS-QB3 value:Δ(IRM)=ΔV(IRM)-ΔV(CBS-QB3)

Δ(IRM)d 1.80 1.76 1.65 1.69 1.77 0.87 1.98 2.35 2.48 3.07 0.38 1.02-1.57-0.33 2.00 1.90 0.91-1.44 2.65-1.66-1.41 0.85-1.31 Reaction R2 R3 R4 R5 R6 R12 R13 R14 R15 R16 R24 R25 R26 R27 R28 R38 R39 R40 R41 R42 R49 R50 R51 Type(1,3)(1,3)(1,3)(1,3)(1,3)(1,4)(1,4)(1,4)(1,4)(1,4)(1,5)(1,5)(1,5)(1,5)(1,5)(1,6)(1,6)(1,6)(1,6)(1,6)(1,7)(1,7)(1,7)ΔVP≠/(kJ?mol-1)CBS-QB3 171.97 170.50 170.38 170.43 170.37 150.26 132.47 147.75 133.27 148.39 86.15 99.35 86.47 73.56 84.57 83.42 92.56 80.13 97.32 69.06 97.49 92.20 95.95 DFTa 181.22 179.79 179.79 179.79 179.66 139.25 153.41 138.57 153.55 124.32 90.67 103.22 92.93 78.79 87.46 84.21 94.35 84.27 97.36 73.41 101.60 94.05 99.96 Δ(DFT)b-9.25-9.29-9.41-9.37-9.29-6.78-5.66-5.30-5.16-4.58-4.51-3.87-6.46-5.22-2.89-0.79-1.79-4.14-0.04-4.36-4.11-1.85-4.01 IRMc 170.16 168.74 168.74 168.74 168.60 131.60 145.77 130.92 145.91 116.67 85.78 98.34 88.04 73.90 82.57 81.51 91.65 81.57 94.67 70.72 98.90 91.36 97.27

圖2中給出了在300-1000 K溫度范圍內(nèi)本文計(jì)算的部分反應(yīng)速率常數(shù)與文獻(xiàn)實(shí)驗(yàn)值和其他計(jì)算值的比較示意圖.對(duì)(1,4)小分子氫遷移反應(yīng)(R11):CH3CH2OO???CH2CH2OOH(見(jiàn)圖2A),Sharma等10用CBS-QB3高精度水平下計(jì)算得到500 K時(shí)速率常數(shù)為5.48×10-4s-1,而Sheng等27用CBS-Q//B3LYP/6-31G(d,p)方法所得值為1.18×10-4s-1,Huynh等28用CBS-QB3方法計(jì)算所得500 K時(shí)速率常數(shù)為2.78×10-4s-1,Baldwin等29得到的實(shí)驗(yàn)值為5.04×10-3s-1,本文通過(guò)等鍵反應(yīng)計(jì)算得到該溫度下常數(shù)為1.26×10-3s-1,與實(shí)驗(yàn)值非常接近.對(duì)(1,4)氫遷移反應(yīng)類型中(R16):(CH3)2CHCH2OO??(CH3)C?CH2OOH(見(jiàn)圖2B),Sharma等10用CBS-QB3方法計(jì)算得500 K時(shí)速率常數(shù)為4.75×10-1s-1,Baldwin等29實(shí)驗(yàn)值2.46×101s-1,而本文通過(guò)等鍵反應(yīng)計(jì)算出速率常數(shù)為2.31 s-1,與實(shí)驗(yàn)值接近.同時(shí),本文用等鍵反應(yīng)方法計(jì)算得到的反應(yīng)R24、R38、R49的速率常數(shù)也與實(shí)驗(yàn)值符合較好,因而本文為過(guò)氧烷基分子內(nèi)氫遷移這一類型反應(yīng)提供了可靠和準(zhǔn)確的動(dòng)力學(xué)數(shù)據(jù).

表5 反應(yīng)動(dòng)力學(xué)參數(shù)(A,Ea)Table 5 Reaction kinetic parameters(A,Ea)

圖2 本文計(jì)算的速率常數(shù)與文獻(xiàn)值比較Fig.2 Comparison of rate constants calculated in this study with those reported in other studies

4 結(jié)論

在B3LYP/6-311+G(d,p)方法下優(yōu)化得到了過(guò)氧烷基分子內(nèi)氫遷移反應(yīng)的過(guò)渡態(tài),提出了根據(jù)過(guò)渡態(tài)反應(yīng)中心幾何結(jié)構(gòu)是否守恒作為反應(yīng)類的判據(jù),并將過(guò)氧烷基分子內(nèi)氫遷移反應(yīng)劃分為(1,3),(1,4),(1,5),((1,n),n=6,7,8)四種反應(yīng)類型.將過(guò)渡態(tài)理論與反應(yīng)類等鍵反應(yīng)方法結(jié)合,對(duì)B3LYP/6-311+G(d,p)方法水平計(jì)算得到的過(guò)氧烷基氫遷移反應(yīng)體系中62個(gè)反應(yīng)的勢(shì)壘、焓變和速率常數(shù)進(jìn)行了修正,所得結(jié)果與高精度CBS-QB3方法計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果一致.因此,本文通過(guò)反應(yīng)類等鍵反應(yīng)方法解決了反應(yīng)類的嚴(yán)格定義和大分子體系的過(guò)氧烷基分子內(nèi)氫遷移反應(yīng)的動(dòng)力學(xué)參數(shù)的精確計(jì)算問(wèn)題,對(duì)碳?xì)浠衔锏牡蜏厝紵龣C(jī)理研究具有重要意義.

Supporting lnformation:Table S1 listed all reactions.Optimized geometrical parameters of the reaction centers of the transition states were given in Table S2.Table S3 listed all reaction barriers,reaction enthalpy changes.This information is available free of chargeviathe internet at http://www.whxb.pku.edu.cn.

(1) Battin-Leclerc,F.Prog.Energy Combust.Sci.2008,34,440.doi:10.1016/j.pecs.2007.10.002

(2) Pousse,E.;Glaude,P.A.;Fournet,R.;Battin-Leclerc,F.Combust.Flame2009,156,954.doi:10.1016/j.combustflame.2008.09.012

(3) Taatjes,C.A.J.Phys.Chem.A2006,110,4299.doi:10.1021/jp056997f

(4) Walker,R.W.;Morley,C.Basic Chemistry of Combustion.InLow-Temperature Combustion and Autoignition;Pilling,M.J.Ed.;Elsevier:Amsterdam,The Netherlands,1997;pp 1-124.

(5) Zádor,J.;Taatjes,C.A.;Fernandes,R.X.Prog.Energy Combust.Sci.2011,37,371.doi:10.1016/j.pecs.2010.06.006

(6) Miller,J.A.;Klippenstein,S.J.;Robertson,H.Proc.Combust.Inst.2000,28,1479.doi:10.1016/S0082-0784(00)80544-5

(7)Villano,S.M.;Huynh,L.K.;Carstensen,H.H.;Dean,A.M.J.Phys.Chem.A2012,116,5068.doi:10.1021/jp3023887

(8)Villano,S.M.;Huynh,L.K.;Carstensen,H.H.;Dean,A.M.J.Phys.Chem.A2011,115,13425.doi:10.1021/jp2079204

(9) Zhang,F.;Dibble,T.S.J.Phys.Chem.A2011,115,655.doi:10.1021/jp1111839

(10)Sharma,S.;Raman,S.;Green,W.H.J.Phys.Chem.A2010,114,5689.doi:10.1021/jp9098792

(11) Miyoshi,A.J.Phys.Chem.A2011,115,3301.doi:10.1021/jp112152n

(12) Truong,T.N.J.Chem.Phys.2000,113,4957.doi:10.1063/1.1287839

(13)Huynh,L.K.;Ratkiewicz,A.;Truong,T.N.J.Phys.Chem.A2006,110,473.doi:10.1021/jp051280d

(14)Muszynska,M.;Ratkiewicz,A.;Huynh,L.K.;Truong,T.N.J.Phys.Chem.A2009,113,8327.doi:10.1021/jp903762x

(15) Bankiewicz,B;Huynh,L.K.;Ratkiewicz,A.;Truong,T.N.J.Phys.Chem.A2009,113,1564.doi:10.1021/jp808874j

(16) Zhang,S.W.;Truong,T.N.J.Phys.Chem.A2003,107,1138.doi:10.1021/jp021265y

(17)Kungwan,N.;Truong,T.N.J.Phys.Chem.A2005,109,7742.doi:10.1021/jp051799+

(18) Huynh,L.K.;Truong,T.N.Theor.Chem.Account2008,120,107.doi:10.1007/s00214-007-0311-9

(19)Wang,B.Y.;Li,Z.R.;Tan,N.X.;Yao,Q.;Li,X.Y.J.Phys.Chem.A2013,117,3279.doi:10.1021/jp400924w

(20)Wang,B.Y.;Tan,L.X.;Yao,Q.;Li,Z.R.;Li,X.Y.Acta Phys.-Chim.Sin.2012,28,2824.[汪必耀,談寧馨,姚 倩,李澤榮,李象遠(yuǎn).物理化學(xué)學(xué)報(bào),2012,28,2824.]doi:10.3866/PKU.WHXB201209053

(21) Frisch,M.J.;Trucks,G.W.;Schlegel,H.B.;et al.Gaussian 03,Revision B.01;Gaussian Inc.:Pittsburgh,PA,2003.

(22) Scott,A.P.;Radom,L.J.Phys.Chem.1996,100,16502.doi:10.1021/jp960976r

(23)NIST.Computational Chemistry Comparison and Benchmark Database.NIST Standard Reference Database Number 101.2013,available at http://webbook.nist.gov/

(24) Montgomery,J.A.;Frisch,M.J.;Ochterski,J.W.;Petersson,G.A.J.Chem.Phys.1999,110,2822.doi:10.1063/1.477924

(25) Zhu,L.;Bozzelli,J.W.;Kardos,L.M.J.Phys.Chem.A2007,111,6361.doi:10.1021/jp070342s

(26) Mokrushin,V.;Tsang,W.Chemratev.1.5.8;National Institute of Standards and Technology:Gaithersburg,MD,2009.

(27)Sheng,C.Y.;Bozzelli,J.W.;Dean,A.M.;Chang,A.Y.J.Phys.Chem.A2002,106,7276.doi:10.1021/jp014540+

(28)Huynh,L.K.;Carstensen,H.H.;Dean,A.M.J.Phys.Chem.A2010,114,6594.

(29)Baldwin,R.R.;Hisham,M.W.M.;Walker,R.W.J.Chem.Soc.Faraday Trans.11982,78,1615.doi:10.1039/f19827801615