噻吩結(jié)構(gòu)與O2反應(yīng)機理的理論研究

戴鳳威,鄧存寶,鄧漢忠,王雪峰,高 飛,張 勛

(1.遼寧工程技術(shù)大學(xué)安全科學(xué)與工程學(xué)院,遼寧阜新 123000;2.遼寧工程技術(shù)大學(xué)材料科學(xué)與工程學(xué)院,遼寧阜新 123000;3.遼寧工程技術(shù)大學(xué)礦業(yè)學(xué)院,遼寧阜新 123000)

噻吩結(jié)構(gòu)與O2反應(yīng)機理的理論研究

戴鳳威1,鄧存寶1,鄧漢忠2,王雪峰1,高 飛1,張 勛3

(1.遼寧工程技術(shù)大學(xué)安全科學(xué)與工程學(xué)院,遼寧阜新 123000;2.遼寧工程技術(shù)大學(xué)材料科學(xué)與工程學(xué)院,遼寧阜新 123000;3.遼寧工程技術(shù)大學(xué)礦業(yè)學(xué)院,遼寧阜新 123000)

為研究煤自燃的反應(yīng)機理,利用Gaussian 03程序,采用密度泛函理論(DFT)方法,在B3LYP/6－311G水平下研究煤結(jié)構(gòu)中噻吩型有機硫與O2反應(yīng)機理。由計算結(jié)果可知,反應(yīng)物的初始連接方式有3種,其中Path 1的反應(yīng)物經(jīng)過過渡態(tài)TS1異構(gòu)化為中間體IM1是無需克服勢壘的放熱過程,這說明煤中的噻吩結(jié)構(gòu)無需從外界吸收熱量,與氧氣接觸即可被氧化。此反應(yīng)主要有6條反應(yīng)路徑,其中,Path 2是反應(yīng)的主要反應(yīng)路徑,其速控步驟為過渡態(tài)TS6,需要克服的勢壘為117.06 kJ/mol,其產(chǎn)物P2(C4H4O+SO)是反應(yīng)的主要產(chǎn)物。

噻吩;反應(yīng)機理;Gaussian 03程序;密度泛函理論

煤是一種具有芳香性的有機大分子,其有機質(zhì)主要是由碳、氫、氧、氮和硫等元素組成。煤中的有機硫賦存形式主要有3種,分別為硫醇、硫醚和噻吩[1－2]。其中,噻吩型有機硫的含量非常多,在煙煤中的體現(xiàn)尤為突出[3－4]。隨著科學(xué)技術(shù)的發(fā)展,科研工作者對煤炭自燃的研究越來越深入,將量子化學(xué)理論應(yīng)用于煤炭研究已經(jīng)成為了一種新的研究方法。謝克昌院士等采用量子化學(xué)原理對煤中含硫化合物的脫除進(jìn)行了研究[5－6]。王繼仁等采用量子化學(xué)原理,從不同角度研究了煤的活性基團(tuán)對氧分子的物理吸附和化學(xué)吸附,研究結(jié)果表明:煤表面?zhèn)孺溁鶊F(tuán)對氧分子的吸附能大于苯環(huán)對氧分子的吸附能[7－9],煤發(fā)生氧化自燃主要是由于含N,S等活性基團(tuán)被氧化引起。且噻吩的硫原子的兩對孤電子中的一對與噻吩的兩個雙鍵共軛,從而形成離域π鍵,因此噻吩具有芳香性,其芳香性僅略弱于苯。噻吩的C—S鍵在光催化作用下易發(fā)生斷裂[10－11],易于發(fā)生氧化或加氫反應(yīng)。因此,煤結(jié)構(gòu)中含硫結(jié)構(gòu)與O2的反應(yīng)機理已成為亟待深入研究的問題。本文采用Gaussian 03程序針對煤中噻吩型有機硫與氧氣的反應(yīng)機理進(jìn)行了研究,為探索預(yù)防煤炭自燃技術(shù)奠定理論基礎(chǔ)。

1 計算方法

全部計算利用Gaussian 03程序進(jìn)行[12],采用密度泛函理論(density functional theory,DFT),首先在B3LYP/6－311G水平下對反應(yīng)物、產(chǎn)物、中間體和過渡態(tài)分子進(jìn)行幾何優(yōu)化,對每一個駐點進(jìn)行振動頻率分析,其中穩(wěn)定構(gòu)型的全部頻率為正,而過渡態(tài)構(gòu)型有且只有一個虛頻[13－14]。然后,在同一水平下進(jìn)行內(nèi)稟反應(yīng)坐標(biāo)(intrinsic reaction coordinate,IRC)計算,確認(rèn)反應(yīng)物、中間體、過渡態(tài)和產(chǎn)物的相關(guān)性,證明過渡態(tài)的正確連接[15－16]。在以下的討論中,將反應(yīng)物(C4H4S+O2)的總能量設(shè)為能量零點。

2 結(jié)果與討論

對于噻吩結(jié)構(gòu)與O2的反應(yīng),計算所得的各條反應(yīng)路徑如下(Path 1～Path 6)所示。

其中,反應(yīng)物用R來表示,反應(yīng)中間體用IM來表示,過渡態(tài)用TS來表示,產(chǎn)物則用P來表示。圖1和圖2分別列出了反應(yīng)物、產(chǎn)物、中間體和過渡態(tài)的幾何構(gòu)型和參數(shù)。表1列出了反應(yīng)物、產(chǎn)物、中間體和過渡態(tài)的總能量(total energies,TE)、零點振動能(zeropoint vibration energies,ZPVE)及相對能量(relative energies,RE)。其中,相對能量RE=相對反應(yīng)物的總能量－相對反應(yīng)物的零點振動能。根據(jù)各駐點的相對能量繪制噻吩結(jié)構(gòu)與O2的反應(yīng)勢能面剖面圖,如圖3所示。從圖3可看出,4種產(chǎn)物中P1,P2和P4的能量均比反應(yīng)物低,熱力學(xué)是可行的;而P3的能量比反應(yīng)物高,熱力學(xué)上是不利的產(chǎn)物。

2.1 反應(yīng)物的初始連接

圖1 在B3LYP/6－311G水平下計算反應(yīng)物和產(chǎn)物的幾何構(gòu)型(鍵長單位:10－10m,鍵角單位:(°))Fig.1 Geometries of reactants and products at the B3LYP/6－311G level(Bond lengths:10－10meters,bond angles:degrees)

噻吩結(jié)構(gòu)和O2反應(yīng)的歷程是復(fù)雜的多步反應(yīng),經(jīng)歷最初的結(jié)合以及隨后的異構(gòu)和離解,最終生成多種產(chǎn)物。根據(jù)Gaussian 03程序計算結(jié)果可知,噻吩結(jié)構(gòu)與O2的結(jié)合方式主要有3種。第1種:兩個O原子進(jìn)攻C4H4S中S原子連接的兩個C原子,形成不穩(wěn)定的過渡態(tài)TS1,這一過程無需克服化學(xué)勢壘,反而為進(jìn)一步的異構(gòu)化和分解反應(yīng)提供了21.17 kJ/mol(圖3)的能量,因此預(yù)測此種結(jié)合是此反應(yīng)的主要反應(yīng)入口;第2種:兩個O原子的電子分別向C4H4S上S—C4鍵上S原子和C4原子靠近,氧分子的雙鍵打開變成單鍵,與S原子和C4原子形成環(huán)狀結(jié)構(gòu),經(jīng)過過渡態(tài)TS7形成穩(wěn)定的中間體IM6,此過程僅需克服15.40 kJ/mol的勢壘,IM6的能量比反應(yīng)物低149.70 kJ/mol,因此,這種結(jié)合方式也是很容易形成的;第3種:兩個O原子共同進(jìn)攻C4H4S的S原子,形成不穩(wěn)定的過渡態(tài)TS12,經(jīng)過過渡態(tài)TS12形成穩(wěn)定的中間體IM8,這個過程需要克服的勢壘為315.99 kJ/mol,很顯然這種結(jié)合方式在反應(yīng)路徑的競爭中處于劣勢。

圖2 在B3LYP/6－311G水平下計算中間體和過渡態(tài)的幾何構(gòu)型(鍵長單位:10－10m,鍵角單位:(°))Fig.2 Geometries of intermediate and transition states at the B3LYP/6－311G level(Bond lengths:10－10meters,bond angles:degrees)

表1 噻吩與O2反應(yīng)的反應(yīng)物、產(chǎn)物、中間體和過渡態(tài)的總能量(TE)、零點振動能(ZPVE)以及相對能量(RE)Table 1 Total energies(TE),zero-point vibration energies(ZPVE),and relative energies(RE)ofreactants, products,intermediate,and transition states on the reaction for thiophene with O2kJ/mol

圖3 在B3LYP/6－311G水平下噻吩結(jié)構(gòu)與O2反應(yīng)的反應(yīng)勢能面剖面圖Fig.3 Schematic profile of the potential energy surface for thiophene structure and O2at the B3LYP/6－311G level

2.2 反應(yīng)路徑分析

(1)噻吩結(jié)構(gòu)與O2反應(yīng)生成產(chǎn)物P1。

Path 1:從圖3可以看出,反應(yīng)物R經(jīng)過過渡態(tài)TS1生成穩(wěn)定的中間體IM1,此過程具體分析見2.1節(jié)。從圖2可以看出,中間體IM1的S—C4鍵拉長,經(jīng)過過渡態(tài)TS2,然后S—C4鍵斷開生成中間體IM2,中間體IM1的S—C4鍵的鍵長為1.939×10－10m,而過渡態(tài)TS2的S—C4鍵長達(dá)到2.878×10－10m,IM1→IM2需要克服的勢壘為72.39 kJ/mol(圖3)。中間體IM2上與C1原子相連的H原子發(fā)生C1—S鄰位遷移,同時O1—O2鍵作用力減弱并斷開,經(jīng)過過渡態(tài)TS3形成中間體IM3,這一過程需要克服的勢壘為72.95 kJ/mol。中間體IM3的O2原子的電子向C1原子靠近,同時S—C1鍵變?nèi)?H—S基從C1原子向C2原子遷移,經(jīng)過過渡態(tài)TS4生成中間體IM4,O2原子與C1原子連接使中間體IM4成為環(huán)狀結(jié)構(gòu),IM3→IM4需要克服的勢壘為90.80 kJ/mol。中間體IM4經(jīng)過旋轉(zhuǎn)過渡態(tài)TS5異構(gòu)化為中間體IM5,這一過程需要克服的勢壘為8.68 kJ/mol。中間體IM5的C1—O2鍵和C1—C2鍵同時拉長,經(jīng)過過渡態(tài)TS6分解為產(chǎn)物P1(HSC2H2CHO+CO),由IM5→P1需要克服很高的勢壘(195.61 kJ/mol),此步在反應(yīng)路徑的競爭中處于劣勢。這條路徑的速控步驟為過渡態(tài)TS6(勢壘為195.61 kJ/mol)。

(2)噻吩結(jié)構(gòu)與O2反應(yīng)生成產(chǎn)物P2。

Path 2:反應(yīng)物R經(jīng)過過渡態(tài)TS7生成穩(wěn)定的中間體IM6,此過程具體分析見2.1節(jié)反應(yīng)物的初始連接。中間體IM6的O1—O2鍵拉長,經(jīng)過過渡態(tài)TS8后斷開,異構(gòu)化為穩(wěn)定的中間體IM7,中間體IM6的O1—O2鍵長為1.483×10－10m,過渡態(tài)TS8的O1—O2鍵長為1.550×10－10m,IM6→IM7過程需要克服的勢壘為6.41 kJ/mol。中間體IM7的O1—S鍵以S—C1鍵為軸旋轉(zhuǎn),經(jīng)過過渡態(tài)TS9異構(gòu)化為中間體IM8,中間體IM7的二面角D(O1SC1C2)為0.25°,過渡態(tài)TS9的二面角D(O1SC1C2)為61.20°,中間體IM8的二面角D(O1SC1C2)為180.00°,由IM7→IM8需要克服的勢壘為60.64 kJ/mol。中間體IM8的C1原子自旋經(jīng)過過渡態(tài)TS10,形成中間體IM9,由IM8→IM9需要克服的勢壘為57.02 kJ/mol。中間體IM9的O2原子開始向C1原子靠近,同時S—C1鍵略微拉長,通過過渡態(tài)TS11后,O2原子和C1原子連接成單鍵,而S—C1鍵斷裂,分解為穩(wěn)定的產(chǎn)物P2,中間體IM9的S—C1鍵長為1.710×10－10m,過渡態(tài)TS11的S—C1鍵長為1.846×10－10m,IM9→P2需要克服的勢壘為117.06 kJ/mol。

Path 3和Path 2的不同點在于反應(yīng)物R經(jīng)過不同的過渡態(tài)形成中間體IM8,在Path 3中,R→IM8經(jīng)過過渡態(tài)TS12,這一過程需要克服勢壘為315.99 kJ/mol,具體分析見2.1節(jié)反應(yīng)物的初始連接。由IM8→P2的過程與Path 2的過程相同,在此不再贅述,具體分析請見上文。

Path 2和Path 3生成相同的產(chǎn)物P2(C4H4O+ SO)。Path2的速控步驟為TS11(勢壘為117.06 kJ/mol);Path 3的速控步驟為TS12(勢壘為315.99 kJ/mol);Path 2的速控步驟反應(yīng)能壘比Path 3的速控步驟反應(yīng)能壘低198.93 kJ/mol,顯然,Path 2在反應(yīng)路徑的競爭中優(yōu)于Path 3;Path 2與Path 1相比,Path 2的速控步驟反應(yīng)能壘比Path 1的速控步驟反應(yīng)能壘低78.55 kJ/mol,因此,Path 2在反應(yīng)路徑的競爭中也優(yōu)于Path 1。

(3)噻吩結(jié)構(gòu)與O2反應(yīng)生成產(chǎn)物P3。

Path 4的反應(yīng)物R異構(gòu)化為中間體IM8的過程與Path 2相同,具體分析請見Path 2。在Path 4中,中間體IM8的S原子與C4原子作用力加強并形成單鍵,同時O2原子從C4原子向S原子遷移,經(jīng)過過渡態(tài)TS13形成中間體IM10,過渡態(tài)TS13的S—C4鍵長為2.035×10－10m,中間體IM10的S—C4鍵長為1.908×10－10m,這一過程需要克服的勢壘相當(dāng)?shù)拇?421.61 kJ/mol),在反應(yīng)路徑的競爭中必處劣勢。中間體IM10的S原子與兩個O原子相連接的作用力很強,而環(huán)狀結(jié)構(gòu)中C1—C2鍵和C3—C4鍵由于共軛鍵的作用作用力也很強,因此S—C1鍵、S—C2鍵和C2—C3鍵同時拉長,經(jīng)過過渡態(tài)TS14分解為產(chǎn)物P3,IM10→P3需要克服的勢壘為244.22 kJ/mol。

Path 5中反應(yīng)物R到中間體IM8的過程和Path 3相同,而Path 5中中間體IM8到產(chǎn)物P3的過程與Path 4相同,在此不再分析,具體分析請見上文。

Path 4和Path 5生成相同的產(chǎn)物P3(2C2H2+ SO2)。Path 4的第1速控步驟為過渡態(tài)TS13(勢壘為421.61 kJ/mol),第2速控步驟為過渡態(tài)TS14(勢壘為244.22 kJ/mol);Path 5的第1速控步驟與Path 4相同,第2速控步驟為過渡態(tài)TS12(勢壘為315.99 kJ/mol);從兩條反應(yīng)路徑的第1速控步驟可以看出,它們在反應(yīng)路徑的競爭中同樣處于劣勢,更無力與Path 2和Path 1相比;Path 4的第2速控步驟反應(yīng)能壘比Path 5的第2速控步驟反應(yīng)能壘低71.77 kJ/mol,顯然,Path 4優(yōu)于Path 5。

(4)噻吩結(jié)構(gòu)與O2反應(yīng)生成產(chǎn)物P4。

Path 6的初始連接方式與Path 2相同,分析同上。Path 6中中間體IM6的S—O1鍵斷裂,經(jīng)過過渡態(tài)TS15異構(gòu)化為中間體IM11,IM6→IM11需要克服的勢壘為123.98 kJ/mol。中間體IM11經(jīng)過一個旋轉(zhuǎn)過渡態(tài)TS16異構(gòu)化為中間體IM12,O1原子以O(shè)2原子為原點繞其旋轉(zhuǎn),異構(gòu)化前二面角D(OC4C1C3)為0.03°,而異構(gòu)化后中間體IM12的二面角D (OC4C1C3)為179.99°,這個過程需要克服的勢壘為150.93 kJ/mol。中間體IM12上與C4原子相連的H原子的電子向O1原子靠近,經(jīng)過過渡態(tài)TS17后,這個H原子與C4原子斷開,與O1原子連接,在H原子向O1原子遷移的過程中,S原子與C4原子的作用力加強形成單鍵,形成更穩(wěn)定的中間體IM13,IM12→IM13需要克服的勢壘為117.89 kJ/mol。中間體IM13經(jīng)過H—O1基從O2原子遷移到C3原子,形成更穩(wěn)定的中間體IM14,IM13→IM14經(jīng)過過渡態(tài)TS18,這個過程需要克服的勢壘為82.24 kJ/mol。中間體IM14的O2—C4鍵是一條雙鍵,O2原子和C4原子相互作用力很強,又因為處在不斷釋放能量的環(huán)境中,S—C4和C3—C4鍵拉長,經(jīng)過過渡態(tài)TS19,兩條鍵斷裂分解為產(chǎn)物P4,IM14→P4需要克服的勢壘為124.17 kJ/mol。

Path 6生成產(chǎn)物P4(HOC2H2CHS+CO),這條路徑的速控步驟為過渡態(tài)TS16(勢壘為150.93 kJ/ mol);與Path 1相比,Path 6速控步驟的反應(yīng)能壘比Path 1速控步驟的低44.68 kJ/mol;而比生成P2的Path 2高33.87 kJ/mol;因此,Path 6在反應(yīng)路徑中的競爭力強于Path 1,Path 2的競爭力優(yōu)于Path 6。

綜上所述,噻吩結(jié)構(gòu)和O2反應(yīng)的6條反應(yīng)路徑的難易程度為:Path 2＞Path 6＞Path 1＞Path 3＞Path 4＞Path 5。

3 結(jié) 論

(1)經(jīng)計算可知,噻吩結(jié)構(gòu)與O2反應(yīng)的6條反應(yīng)路徑中,路徑Path 2(R→TS7→IM6→TS8→IM7→TS9→IM8→TS10→IM9→TS11→P2)的反應(yīng)勢壘最低,為主反應(yīng)路徑,其次就是路徑Path 6(R→TS7→IM6→TS15→IM11→TS16→IM12→TS17→IM13→TS18→IM14→TS19→P4)。

(2)計算結(jié)果表明,噻吩結(jié)構(gòu)與O2反應(yīng)的4種產(chǎn)物中,產(chǎn)物P2(C4H4O+SO)是反應(yīng)的主產(chǎn)物, P3(2C2H2+SO2)是最難生成的產(chǎn)物。

(3)煤結(jié)構(gòu)中噻吩結(jié)構(gòu)無需從外界吸收熱量,與氧氣接觸即可被氧化。

[1] Gryglewicz G,Jasieńko S.Sulfur groups in the cokes obtained from coals of different ranks[J].Fuel Process Technol.,1988,19(1): 51－59.

[2] Miura K,Mae K,Shimada M,et al.Analysis of formation rates of sulfur-containing gases during the pyrolysis of various coals[J].Energ.Fuel,2001,15(3):629－636.

[3] 孫成功,李保慶,Snape C E.煤中有機硫形態(tài)結(jié)構(gòu)和熱解過程硫變遷特性的研究[J].燃料化學(xué)學(xué)報,1997,25(4):358－362.

Sun Chenggong,Li Baoqing,Snape C E.Characterization of organic sulfur forms in some chinese coals by high pressure TPR and sulphur transfer during hydropyrolysis[J].Journal of Fuel Chemistry and Technology,1997,25(4):358－362.

[4] 謝建軍,楊學(xué)民,呂雪松,等.煤熱解過程中硫氮分配及遷移規(guī)律研究進(jìn)展[J].化工進(jìn)展,2004,23(11):1214－1218.

Xie Jianjun,Yang Xuemin,Lü Xuesong,et al.Progress on transformation behavior of sulfur and nitrogen during coal pyrolysis[J].Chemical Industry and Engineering Progress,2004,23(11):1214－1218.

[5] Ling Lixia,Zhang Riguang,Wang Baojun,et al.Density functional theory study on the pyrolysis mechanism of thiophene in coal[J].Journal of Molecular Structure:THEOCHEM,2009,905(1):8－12.[6] Ling Lixia,Zhang Riguang,Wang Baojun,et al.DFT study on the sulfur migration during benzenethiol pyrolysis in coal[J].Journal of Molecular Structure:THEOCHEM,2010,952(1):31－35.

[7] 鄧存寶,王繼仁,鄧漢忠,等.氧在煤表面—CH2—NH2基團(tuán)上的化學(xué)吸附[J].煤炭學(xué)報,2009,34(9):1234－1238.

Deng Cunbao,Wang Jiren,Deng Hanzhong,et al.Chemical adsorption of O2adsorbed in the coal surface—CH2—NH2group[J].Jounal of China Coal Society,2009,34(9):1234－1238.

[8] 鄧存寶,王雪峰,王繼仁,等.煤表面含S側(cè)鏈基團(tuán)對氧分子的物理吸附機理[J].煤炭學(xué)報,2008,33(5):556－560.

Deng Cunbao,Wang Xuefeng,Wang Jiren,et al.Physical adsorption mechanism of coal surface containing sulfur group adsorption to more oxygen molecule[J].Jounal of China Coal Society,2008,33(5): 556－560.

[9] 鄧存寶,王繼仁,葉 兵,等.煤表面對單氧分子的物理吸附機理[J].中國礦業(yè)大學(xué)學(xué)報,2008,37(2):171－175.

Deng Cunbao,Wang Jiren,Ye Bing,et al.Physical mechanism of a single oxygen molecule adsorbs to the coal surface[J].Journal of China University of Mining&Technology,2008,37(2):171－175.

[10] 亓 雪,石秋杰,諶偉慶,等.Mo對非晶態(tài)合金Ni－B/薄水鋁石催化劑上噻吩加氫脫硫性能的影響[J].催化學(xué)報,2012,33 (3):543－549.

Qi Xue,Shi Qiujie,Chen Weiqing,et al.Effect of Mo on performance of Ni-B/boehmite amorphous alloy catalyst for thiophene hydrodesulfurization[J].Chinese Journal of Catalysis,2012,33 (3):543－549.

[11] 高曉明,付 峰,武玉飛,等.Co－BiVO4光催化劑的制備及其用于光催化氧化噻吩[J].無機材料學(xué)報,2012,27(10):1073－1078.

Gao Xiaoming,Fu Feng,Wu Yufei,et al.Preparation of Co-BiVO4photocatalyst and its application in the photocatalytic oxidative thiophene[J].Journal of Inorganic Materials,2012,27(10):1073－1078.

[12] Frisch M J,Trucks G W,Schlegel H B,et al.Gaussian 03w,version6.0,Gaussian[M].Inc.,Wallingford,CT,2004.

[13] Govind N,Petersen M,Fitzgerald G,et al.A generalized synchronous transit method for transition state location[J].Computational Materials Science,2003,28(2):250－258.

[14] Malick D K,Petersson G A,Montgomery Jr J A.Transition states for chemical reactions I.Geometry and classical barrier height[J].The Journal of Chemical Physics,1998,108(14):5704－5713.

[15] Bakken V,Helgaker T.The efficient optimization of molecular geometries using redundant internal coordinates[J].The Journal of Chemical Physics,2002,117(20):9160－9174.

[16] Ayala P Y,Schlegel H B.A combined method for determining reaction paths,minima,and transition state geometries[J].The Journal of Chemical Physics,1997,107(2):375－384.

Theoretical study on the mechanism of the thiophene structure with O2reaction

DAI Feng-wei1,DENG Cun-bao1,DENG Han-zhong2,WANG Xue-feng1,GAO Fei1,ZHANG Xun3
(1.College of Safety Science and Engineering,Liaoning Technical University,Fuxin 123000,China;2.College of Materials Science and Engineering,Liaoning Technical University,Fuxin 123000,China;3.College of Mining Engineering,Liaoning Technical University,Fuxin 123000,China)

In order to study the reaction mechanism of coal spontaneous combustion,the calculations were performed with the Gaussian 03 program.At the B3LYP/6－311G level,the reaction mechanism of thiophene type organic sulfur in coal structure with O2was studied by using density functional theory(DFT)method.The calculations show that initial connection of the reactant has three forms,the reactant of path 1 in it is readily isomerized to the intermediate IM1through the transition state TS1,this process is exothermic process without overcome potential barrier,it show that thiophene type organic sulfur in coal contact with oxygen to be oxidized without absorbing heat from outside;this reaction has six reaction paths,the path 2 is the main reaction path,its rate-controlling step is the transition state TS6,its reaction energy barrier is 117.06 kJ/mol,P2(C4H4O+SO)is the main product.

thiophene;reaction mechanism;Gaussian 03 program;density functional theory

TQ530

A

0253－9993(2014)04－0699－06

戴鳳威,鄧存寶,鄧漢忠,等.噻吩結(jié)構(gòu)與O2反應(yīng)機理的理論研究[J].煤炭學(xué)報,2014,39(4):699－704.

10.13225/j.cnki.jccs.2013.0562

Dai Fengwei,Deng Cunbao,Deng Hanzhong,et al.Theoretical study on the mechanism of the thiophene structure with O2reaction[J].Journal of China Coal Society,2014,39(4):699－704.doi:10.13225/j.cnki.jccs.2013.0562

2013－05－02 責(zé)任編輯:畢永華

國家自然科學(xué)基金面上資助項目(51274112);國家自然科學(xué)基金資助項目(51174108)

戴鳳威(1987—),女,遼寧遼陽人,博士研究生。Tel:0418－3351702,E－mail:daifengwei214@163.com