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

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

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

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

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

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

1 計(jì)算方法

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

2 結(jié)果與討論

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

其中,反應(yīng)物用R來(lái)表示,反應(yīng)中間體用IM來(lái)表示,過(guò)渡態(tài)用TS來(lái)表示,產(chǎn)物則用P來(lái)表示。圖1和圖2分別列出了反應(yīng)物、產(chǎn)物、中間體和過(guò)渡態(tài)的幾何構(gòu)型和參數(shù)。表1列出了反應(yīng)物、產(chǎn)物、中間體和過(guò)渡態(tài)的總能量(total energies,TE)、零點(diǎn)振動(dòng)能(zeropoint vibration energies,ZPVE)及相對(duì)能量(relative energies,RE)。其中,相對(duì)能量RE=相對(duì)反應(yīng)物的總能量－相對(duì)反應(yīng)物的零點(diǎn)振動(dòng)能。根據(jù)各駐點(diǎn)的相對(duì)能量繪制噻吩結(jié)構(gòu)與O2的反應(yīng)勢(shì)能面剖面圖,如圖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水平下計(jì)算反應(yīng)物和產(chǎn)物的幾何構(gòu)型(鍵長(zhǎng)單位: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程序計(jì)算結(jié)果可知,噻吩結(jié)構(gòu)與O2的結(jié)合方式主要有3種。第1種:兩個(gè)O原子進(jìn)攻C4H4S中S原子連接的兩個(gè)C原子,形成不穩(wěn)定的過(guò)渡態(tài)TS1,這一過(guò)程無(wú)需克服化學(xué)勢(shì)壘,反而為進(jìn)一步的異構(gòu)化和分解反應(yīng)提供了21.17 kJ/mol(圖3)的能量,因此預(yù)測(cè)此種結(jié)合是此反應(yīng)的主要反應(yīng)入口;第2種:兩個(gè)O原子的電子分別向C4H4S上S—C4鍵上S原子和C4原子靠近,氧分子的雙鍵打開(kāi)變成單鍵,與S原子和C4原子形成環(huán)狀結(jié)構(gòu),經(jīng)過(guò)過(guò)渡態(tài)TS7形成穩(wěn)定的中間體IM6,此過(guò)程僅需克服15.40 kJ/mol的勢(shì)壘,IM6的能量比反應(yīng)物低149.70 kJ/mol,因此,這種結(jié)合方式也是很容易形成的;第3種:兩個(gè)O原子共同進(jìn)攻C4H4S的S原子,形成不穩(wěn)定的過(guò)渡態(tài)TS12,經(jīng)過(guò)過(guò)渡態(tài)TS12形成穩(wěn)定的中間體IM8,這個(gè)過(guò)程需要克服的勢(shì)壘為315.99 kJ/mol,很顯然這種結(jié)合方式在反應(yīng)路徑的競(jìng)爭(zhēng)中處于劣勢(shì)。

圖2 在B3LYP/6－311G水平下計(jì)算中間體和過(guò)渡態(tài)的幾何構(gòu)型(鍵長(zhǎng)單位: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)物、中間體和過(guò)渡態(tài)的總能量(TE)、零點(diǎn)振動(dòng)能(ZPVE)以及相對(duì)能量(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)勢(shì)能面剖面圖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)過(guò)過(guò)渡態(tài)TS1生成穩(wěn)定的中間體IM1,此過(guò)程具體分析見(jiàn)2.1節(jié)。從圖2可以看出,中間體IM1的S—C4鍵拉長(zhǎng),經(jīng)過(guò)過(guò)渡態(tài)TS2,然后S—C4鍵斷開(kāi)生成中間體IM2,中間體IM1的S—C4鍵的鍵長(zhǎng)為1.939×10－10m,而過(guò)渡態(tài)TS2的S—C4鍵長(zhǎng)達(dá)到2.878×10－10m,IM1→IM2需要克服的勢(shì)壘為72.39 kJ/mol(圖3)。中間體IM2上與C1原子相連的H原子發(fā)生C1—S鄰位遷移,同時(shí)O1—O2鍵作用力減弱并斷開(kāi),經(jīng)過(guò)過(guò)渡態(tài)TS3形成中間體IM3,這一過(guò)程需要克服的勢(shì)壘為72.95 kJ/mol。中間體IM3的O2原子的電子向C1原子靠近,同時(shí)S—C1鍵變?nèi)?H—S基從C1原子向C2原子遷移,經(jīng)過(guò)過(guò)渡態(tài)TS4生成中間體IM4,O2原子與C1原子連接使中間體IM4成為環(huán)狀結(jié)構(gòu),IM3→IM4需要克服的勢(shì)壘為90.80 kJ/mol。中間體IM4經(jīng)過(guò)旋轉(zhuǎn)過(guò)渡態(tài)TS5異構(gòu)化為中間體IM5,這一過(guò)程需要克服的勢(shì)壘為8.68 kJ/mol。中間體IM5的C1—O2鍵和C1—C2鍵同時(shí)拉長(zhǎng),經(jīng)過(guò)過(guò)渡態(tài)TS6分解為產(chǎn)物P1(HSC2H2CHO+CO),由IM5→P1需要克服很高的勢(shì)壘(195.61 kJ/mol),此步在反應(yīng)路徑的競(jìng)爭(zhēng)中處于劣勢(shì)。這條路徑的速控步驟為過(guò)渡態(tài)TS6(勢(shì)壘為195.61 kJ/mol)。

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

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

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

Path 2和Path 3生成相同的產(chǎn)物P2(C4H4O+ SO)。Path2的速控步驟為T(mén)S11(勢(shì)壘為117.06 kJ/mol);Path 3的速控步驟為T(mén)S12(勢(shì)壘為315.99 kJ/mol);Path 2的速控步驟反應(yīng)能壘比Path 3的速控步驟反應(yīng)能壘低198.93 kJ/mol,顯然,Path 2在反應(yīng)路徑的競(jìng)爭(zhē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)路徑的競(jìng)爭(zhēng)中也優(yōu)于Path 1。

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

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

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

Path 4和Path 5生成相同的產(chǎn)物P3(2C2H2+ SO2)。Path 4的第1速控步驟為過(guò)渡態(tài)TS13(勢(shì)壘為421.61 kJ/mol),第2速控步驟為過(guò)渡態(tài)TS14(勢(shì)壘為244.22 kJ/mol);Path 5的第1速控步驟與Path 4相同,第2速控步驟為過(guò)渡態(tài)TS12(勢(shì)壘為315.99 kJ/mol);從兩條反應(yīng)路徑的第1速控步驟可以看出,它們?cè)诜磻?yīng)路徑的競(jìng)爭(zhēng)中同樣處于劣勢(shì),更無(wú)力與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)過(guò)過(guò)渡態(tài)TS15異構(gòu)化為中間體IM11,IM6→IM11需要克服的勢(shì)壘為123.98 kJ/mol。中間體IM11經(jīng)過(guò)一個(gè)旋轉(zhuǎn)過(guò)渡態(tài)TS16異構(gòu)化為中間體IM12,O1原子以O(shè)2原子為原點(diǎn)繞其旋轉(zhuǎn),異構(gòu)化前二面角D(OC4C1C3)為0.03°,而異構(gòu)化后中間體IM12的二面角D (OC4C1C3)為179.99°,這個(gè)過(guò)程需要克服的勢(shì)壘為150.93 kJ/mol。中間體IM12上與C4原子相連的H原子的電子向O1原子靠近,經(jīng)過(guò)過(guò)渡態(tài)TS17后,這個(gè)H原子與C4原子斷開(kāi),與O1原子連接,在H原子向O1原子遷移的過(guò)程中,S原子與C4原子的作用力加強(qiáng)形成單鍵,形成更穩(wěn)定的中間體IM13,IM12→IM13需要克服的勢(shì)壘為117.89 kJ/mol。中間體IM13經(jīng)過(guò)H—O1基從O2原子遷移到C3原子,形成更穩(wěn)定的中間體IM14,IM13→IM14經(jīng)過(guò)過(guò)渡態(tài)TS18,這個(gè)過(guò)程需要克服的勢(shì)壘為82.24 kJ/mol。中間體IM14的O2—C4鍵是一條雙鍵,O2原子和C4原子相互作用力很強(qiáng),又因?yàn)樘幵诓粩噌尫拍芰康沫h(huán)境中,S—C4和C3—C4鍵拉長(zhǎng),經(jīng)過(guò)過(guò)渡態(tài)TS19,兩條鍵斷裂分解為產(chǎn)物P4,IM14→P4需要克服的勢(shì)壘為124.17 kJ/mol。

Path 6生成產(chǎn)物P4(HOC2H2CHS+CO),這條路徑的速控步驟為過(guò)渡態(tài)TS16(勢(shì)壘為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)路徑中的競(jìng)爭(zhēng)力強(qiáng)于Path 1,Path 2的競(jìng)爭(zhēng)力優(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)計(jì)算可知,噻吩結(jié)構(gòu)與O2反應(yīng)的6條反應(yīng)路徑中,路徑Path 2(R→TS7→IM6→TS8→IM7→TS9→IM8→TS10→IM9→TS11→P2)的反應(yīng)勢(shì)壘最低,為主反應(yīng)路徑,其次就是路徑Path 6(R→TS7→IM6→TS15→IM11→TS16→IM12→TS17→IM13→TS18→IM14→TS19→P4)。

(2)計(jì)算結(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)無(wú)需從外界吸收熱量,與氧氣接觸即可被氧化。

[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.煤中有機(jī)硫形態(tài)結(jié)構(gòu)和熱解過(guò)程硫變遷特性的研究[J].燃料化學(xué)學(xué)報(bào),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ò)程中硫氮分配及遷移規(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é)報(bào),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)對(duì)氧分子的物理吸附機(jī)理[J].煤炭學(xué)報(bào),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] 鄧存寶,王繼仁,葉 兵,等.煤表面對(duì)單氧分子的物理吸附機(jī)理[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào),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對(duì)非晶態(tài)合金Ni－B/薄水鋁石催化劑上噻吩加氫脫硫性能的影響[J].催化學(xué)報(bào),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].無(wú)機(jī)材料學(xué)報(bào),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ī)理的理論研究[J].煤炭學(xué)報(bào),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é)任編輯:畢永華

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

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