任 奎， 龍 軍， 任 強(qiáng)， 李永祥， 代振宇

(中國石化 石油化工科學(xué)研究院 催化材料與反應(yīng)工程國家重點(diǎn)實(shí)驗(yàn)室， 北京 100083)

隨著我國第六階段車用汽油標(biāo)準(zhǔn)的實(shí)施，汽油中烯烴、芳烴和苯的含量進(jìn)一步降低[1]。我國汽油池中，催化裂化汽油和催化重整汽油所占比例已超過90%，烷基化和異構(gòu)化汽油比例不足1%，這導(dǎo)致了我國汽油標(biāo)準(zhǔn)中烯烴和芳烴含量長期高于歐美等發(fā)達(dá)國家[2-5]。而烷基化汽油以異構(gòu)烷烴為主(約90%以上)，具有辛烷值高和不排放SOx、NOx等優(yōu)點(diǎn)[6-8]。因此，提高烷基化汽油的比例，是解決我國汽油質(zhì)量升級的有效途徑。

傳統(tǒng)的液體酸烷基化技術(shù)(硫酸法和氫氟酸法)雖然已工業(yè)化多年，卻一直存在酸耗高、毒性大等諸多問題[9-11]。而固體酸催化劑，尤其是分子篩催化劑，具有安全、易分離等特點(diǎn)，工業(yè)化潛力較大[12]。世界上第一套商業(yè)規(guī)模的固體酸烷基化裝置于2015年8月在中國山東淄博成功啟動，且該技術(shù)獲得了2016年“美國總統(tǒng)綠色化學(xué)挑戰(zhàn)獎”中的綠色合成路線獎[13]。

異丁烷/丁烯烷基化是一個鏈?zhǔn)椒磻?yīng)，包含了質(zhì)子化、烯烴加成、氫負(fù)離子轉(zhuǎn)移和異構(gòu)化等反應(yīng)[6,8]。其中，碳正離子與丁烯異構(gòu)體之間的加成反應(yīng)是生產(chǎn)烷基化汽油的主要步驟。在異丁烷/丁烯烷基化體系中所涉及到的所有加成反應(yīng)中，叔丁基碳正離子與異丁烯和2-丁烯間的加成反應(yīng)可以直接生成帶3個甲基支鏈的C8碳正離子，這樣的碳正離子經(jīng)過與異丁烷之間的氫負(fù)離子轉(zhuǎn)移反應(yīng)即可生成三甲基戊烷(TMP)，這是烷基化反應(yīng)的理想產(chǎn)物。因此，十分有必要對這兩個基元反應(yīng)的過程進(jìn)行詳細(xì)研究。傳統(tǒng)的叔丁基碳正離子與異丁烯/2-丁烯的加成反應(yīng)機(jī)理如圖1所示。由于在液體酸中，碳正離子是獨(dú)立存在的，故圖1可以用來解釋液體酸為催化劑下的加成反應(yīng)機(jī)理。但是，文獻(xiàn)[14-16]表明，在固體酸催化劑(如分子篩)催化的反應(yīng)體系中，碳正離子以“烷氧基團(tuán)”的構(gòu)型存在于分子篩孔道內(nèi)，即中心C原子與分子篩骨架上的O原子以共價鍵的方式相連，碳正離子不再是獨(dú)立存在的，而是限定在O原子附近的一定區(qū)域內(nèi)。因此，對于固體酸催化下的烷基化反應(yīng)，考慮催化劑性質(zhì)對反應(yīng)過程的影響，并采用相應(yīng)固體酸催化劑模型來深入研究反應(yīng)機(jī)理至關(guān)重要。

與實(shí)驗(yàn)技術(shù)相比，分子模擬技術(shù)可以從分子和原子尺度，對催化劑和反應(yīng)物建立模型，并對反應(yīng)過程的細(xì)節(jié)進(jìn)行詳細(xì)探究。研究表明，在分子篩孔道內(nèi)，叔丁基碳正離子以叔丁基烷氧基團(tuán)(t-butylalkoxide, TBA)的形式存在，該結(jié)構(gòu)中的C—O鍵鍵長大于標(biāo)準(zhǔn)C—O單鍵鍵長[14, 17]，故具有一定的反應(yīng)活性。因此，本研究擬采用分子模擬的方法，以HY分子篩為催化劑模型，對異丁烷/丁烯烷基化反應(yīng)體系中TBA分別與異丁烯和2-丁烯發(fā)生加成反應(yīng)的過程進(jìn)行深入探究，并比較二者異同。

圖1 叔丁基碳正離子與異丁烯和2-丁烯加成反應(yīng)的 反應(yīng)機(jī)理示意圖Fig.1 Mechanism schematic of olefin addition reaction of t-butyl carbonium ion and isobutene/2-butene

1 Y分子篩模型的建立和計(jì)算方法

1.1 模型的建立

本研究中所用模型為120T HY分子篩模型，該模型截取自FAU分子篩晶體庫[18]，包含由1個十二元氧環(huán)連接的2個超籠，如圖2所示。根據(jù)理論計(jì)算和實(shí)驗(yàn)測定，B酸位點(diǎn)的H原子與O1原子相連[19-20]。模型的邊界用H原子飽和，形成懸掛的Si—H鍵。為了保持分子篩的固有特征，所有的Si—H懸掛鍵與原晶體骨架中的Si—O鍵的方向一致，Si—H鍵長固定為0.147 nm。

由于120T的簇模型較大，采用全密度泛函(DFT)方法進(jìn)行計(jì)算非常耗時，而且不切實(shí)際。因此，對這樣的模型，應(yīng)該優(yōu)先考慮采用量子力學(xué)與分子力學(xué)(Quantum mechanics/molecular mechanics, QM/MM)相結(jié)合的方法[21]。該方法將模型繼續(xù)分成2個部分，內(nèi)部區(qū)域和外部區(qū)域。其中，內(nèi)部區(qū)域是反應(yīng)活性區(qū)域，包含分子篩的活性位點(diǎn)和反應(yīng)物分子；外部區(qū)域是分子篩空間環(huán)境的代表。B酸位點(diǎn)周圍的一些原子和反應(yīng)物分子為QM區(qū)域，其余部分為MM區(qū)域。前者以機(jī)械嵌入的方式嵌套在分子篩骨架中，采用DFT方法進(jìn)行計(jì)算，所用模塊為DMol3[22]；后者采用分子力學(xué)方法進(jìn)行計(jì)算，所用模塊為GULP[23]。為了保證分子篩骨架不發(fā)生畸變，模型中所有QM區(qū)域原子均處于弛域狀態(tài)，而所有的非QM區(qū)域的原子都固定在其晶體結(jié)構(gòu)位置上。

圖2 計(jì)算所用120T HY分子篩模型的結(jié)構(gòu)示意圖Fig.2 Structure schematic of the 120T HY zeolite model used in calculation (a) Front view; (b) Side view

計(jì)算所使用的軟件為Biovio公司開發(fā)的Material Studios 2016。

1.2 計(jì)算參數(shù)

在DFT計(jì)算中，采用Dmol3模塊中的密度泛函GGA-PBE方法在DNP基組水平上完成內(nèi)層活性區(qū)域原子和反應(yīng)物分子的計(jì)算。能量、梯度和位移的收斂閾值分別為0.05 kJ/mol、1.05×10-11N/mol和0.0005 nm，所有計(jì)算均采用DIIS加速收斂。在QM/MM的計(jì)算中，采用UFF力場[24]，力場中原子所帶電荷采用QEq電荷[25]。采用Additive方式計(jì)算能量時，體系總能量ETotal由下式計(jì)算[24,26]。

ETotal=EMM+EQM+EQMMM

其中，EQM和EMM分別為QM和MM計(jì)算所得到的原子能量，kJ/mol；EQMMM為處于QM和MM邊界上的原子對整個體系能量的貢獻(xiàn)，kJ/mol。

QMERA模塊中過渡態(tài)的搜索采用NEB的方法[27]。該方法通過優(yōu)化2種狀態(tài)之間的一系列中間產(chǎn)物的構(gòu)象，可以在已知的反應(yīng)物和產(chǎn)物之間找到最小能量路徑，過渡態(tài)位于最低能量路徑的鞍點(diǎn)上。

在該方法中，對QM區(qū)域進(jìn)行相關(guān)計(jì)算前，首先會對QM區(qū)域的邊界進(jìn)行氫飽和處理。因此，氫的加入會對活性位點(diǎn)的性質(zhì)產(chǎn)生影響，即邊界效應(yīng)。為了消除邊界效應(yīng)，筆者采用擴(kuò)大QM區(qū)域的方法來消除[28]。采用22T～36T不同大小的QM區(qū)域，對TBA與異丁烯之間加成反應(yīng)的活化能進(jìn)行計(jì)算，活化能隨QM區(qū)域大小的變化如圖3所示。由圖3可以看到，隨著QM區(qū)域原子的增加，反應(yīng)能壘大幅度降低。當(dāng)QM區(qū)域擴(kuò)大到32T時，活化能量值收斂，即此時的QM區(qū)域的大小足以消除邊界效應(yīng)，故本研究中QM區(qū)域選擇32T。

圖3 TBA與異丁烯之間加成反應(yīng)的活化能 隨QM區(qū)域大小變化的曲線Fig.3 The reaction barrier of addition between TBA and isobutane as a function of QM region

2 結(jié)果與討論

2.1 TBA與異丁烯的加成反應(yīng)

為了確定TBA與異丁烯的作用位點(diǎn)，采用前線軌道方法對體系進(jìn)行分析。圖4為TBA和異丁烯的前線軌道分布。由圖4可以看出，TBA的HOMO主要分布在Al原子所在的六方柱籠內(nèi)，主要是與Al相鄰近的Si、O原子上；LUMO主要分布在TBA的中心C原子和與之相連的O原子上。從圖4還可以看出，異丁烯分子的HOMO主要分布在C=C雙鍵以及甲基的C—H原子上，LUMO主要分布在C—C鍵和C—H鍵上。TBA叔碳原子的LUMO與異丁烯分子中C=C雙鍵的HOMO對稱性匹配；而HOMO軌道電子能量最高，最為活潑，給電子能力強(qiáng)，容易將電子傳遞到對稱性匹配的LUMO軌道，因此，加成反應(yīng)優(yōu)先發(fā)生在二者之間。

圖5為TBA與異丁烯發(fā)生加成反應(yīng)過程各階段的構(gòu)象示意圖。反應(yīng)過程中體系內(nèi)部分原子間的距離和電荷如表1～2所示。

圖4 TBA和異丁烯的前線軌道分布Fig.4 The HOMO and LUMO of TBA and isobutene (a) LUMO of TBA (overall structure); (b) LUMO of TBA (local structure); (c) HOMO of TBA (overall structure); (d) HOMO of TBA (local structure); (e) LUMO of isobutene (overall structure); (f) HOMO of isobutene (overall structure)

圖5 TBA與異丁烯發(fā)生加成反應(yīng)過程各階段的構(gòu)象示意圖Fig.5 Structures of various steps during addition reaction of TBA and isobutene (a) TBA+Isobutene (g); (b) TBA+Isobutene*; (c) TS-1; (d) C8 carbonium ion (1); (e) Rotate 180o; (f) C8 alkoxide (1) g——Isolated state; *——Adsorbed state

由圖5和表1、表2可以看到，當(dāng)異丁烯處于游離狀態(tài)時(見圖5(a))，異丁烯分子整體為電中性，即此時異丁烯分子與TBA之間無相互作用。隨著異丁烯分子靠近TBA，二者之間會產(chǎn)生較弱的相互作用，形成吸附態(tài)的異丁烯分子，其結(jié)構(gòu)如圖5(b)所示。此時，異丁烯分子整體所帶電荷增加為+0.02 e，說明異丁烯已經(jīng)向分子篩骨架傳遞了少量電子。由于二者之間的相互作用，使得C2—O1鍵長

表1 TBA與異丁烯發(fā)生加成反應(yīng)過程中體系內(nèi)部分原子間距離Table 1 Interatomic distance between atoms during addition reaction of TBA and isobutene

g——Isolated state; *——Adsorbed state

表2 TBA與異丁烯發(fā)生加成反應(yīng)過程中體系內(nèi)部分原子及原子團(tuán)的電荷Table 2 Charge of atoms and atomic groups during addition reaction of TBA and isobutene

g——Isolated state; *——Adsorbed state

從0.163 nm增加為0.166 nm，烷氧基團(tuán)發(fā)生反應(yīng)的活性增加，有遠(yuǎn)離活性位點(diǎn)的趨勢。在該吸附態(tài)構(gòu)象中，C5原子所帶電荷為-0.23 e，而C2原子電荷為+0.09 e，二者之間存在較強(qiáng)的靜電相互作用。導(dǎo)致C2原子遠(yuǎn)離O1原子，靠近C5原子，形成吸附態(tài)的叔丁基碳正離子過渡態(tài)結(jié)構(gòu)(TS1)，如圖5(c)所示。在該TS1中，C2—O1原子間的距離由0.166 nm增加為0.258 nm，C2—C5原子間距離由0.395 nm減小為0.301 nm。此時，異丁烯分子整體所帶電荷從+0.02 e增加為+0.10 e，叔丁基碳正離子整體所帶電荷為+0.76 e，說明在這個過程中，異丁烯與TBA之間發(fā)生了明顯的電子轉(zhuǎn)移，電子從異丁烯C=C雙鍵的HOMO傳遞到了TBA的C2原子的LUMO。在TS1中，C2、C5原子所帶電荷分別為+0.18 e 和-0.24 e，二者之間的靜電作用有較大程度的增強(qiáng)，在此作用下，叔丁基碳正離子的中心C原子與異丁烯的C5原子會繼續(xù)發(fā)生電子轉(zhuǎn)移，從而形成C—C單鍵，這樣就生成了如圖5(d)所示的C8碳正離子結(jié)構(gòu)，該結(jié)構(gòu)中，C2—C5原子間距離減小為0.171 nm，C2—O1原子間距離增加為0.353 nm，中心C原子由C2轉(zhuǎn)移為C6，C6原子所帶電荷為+0.19 e。此后，該C8碳正離子會在活性中心附近以吸附狀態(tài)的形式發(fā)生旋轉(zhuǎn)，并以C6為中心C原子吸附在活性位點(diǎn)O1上(見圖5(e))，此時，C6、O1原子所帶電荷分別為+0.22 e、-0.88 e，二者之間的靜電作用使得C6原子極易與O1原子成鍵，從而形成如圖5(f)所示的C8烷氧基團(tuán)。

圖6為TBA與異丁烯加成反應(yīng)的能量示意圖。由圖6可以看到，游離態(tài)異丁烯分子進(jìn)入分子篩孔道形成吸附態(tài)異丁烯分子的過程放熱為34.6 kJ/mol，克服14.9 kJ/mol的能壘后，形成叔丁基碳正離子的過渡態(tài)，過渡態(tài)的叔丁基碳正離子與異丁烯形成C8碳正離子，體系能量降低了50.5 kJ/mol，隨后C8碳正離子在活性位點(diǎn)旋轉(zhuǎn)180°，并以新的中心C原子吸附在活性氧原子O1上，最終生成C8烷氧基團(tuán)，從而完成加成反應(yīng)。

圖6 TBA與異丁烯加成反應(yīng)的能量示意圖Fig.6 Energy diagram for addition reaction of TBA and isobutene Energies are given relative to separated 2-butane and TBA. g——Isolated state; *——Adsorbed state

2.2 TBA與2-丁烯之間的加成反應(yīng)

為了確定TBA與2-丁烯之間發(fā)生作用的位點(diǎn)，采用前線軌道方法對體系進(jìn)行分析，圖7為TBA和2-丁烯的前線軌道分布。從圖7(a)～(d)可以看出，TBA的HOMO主要分布在Al原子所在的六方柱籠內(nèi)，主要是與Al相鄰近的Si、O原子上；LUMO主要分布在TBA的中心C原子和與之相連的O原子上。從圖7(e)～(f)可以看出，2-丁烯的HOMO主要分布在C=C雙鍵以及甲基的C—H原子上，LUMO主要分布在C—H鍵上。由此可知，TBA的C2原子的LUMO與2-丁烯分子中C=C雙鍵的HOMO對稱性匹配；而HOMO軌道電子能量最高，給電子能力強(qiáng)，容易將電子傳遞到對稱性匹配的LUMO軌道，因此，加成反應(yīng)優(yōu)先發(fā)生在二者之間。

圖8為TBA與2-丁烯發(fā)生加成反應(yīng)過程各階段的構(gòu)象示意圖。反應(yīng)過程中體系內(nèi)部分原子間的距離和電荷如表3～4所示。

圖7 TBA和2-丁烯的前線軌道分布Fig.7 The HOMO and LUMO of TBA and 2-butene (a) LUMO of TBA (overall structure); (b) LUMO of TBA (local structure); (c) HOMO of TBA (overall structure); (d) HOMO of TBA (local structure); (e) LUMO of 2-butene (overall structure); (f) HOMO of 2-butene (overall structure)

圖8 TBA與2-丁烯發(fā)生加成反應(yīng)過程各階段的構(gòu)象示意圖Fig.8 Structures of various steps during addition reaction of TBA and 2-butene (a) TBA+2-Butene (g); (b) TBA+2-Butene*; (c) TS-2; (d) Triangle-structured intermediate; (e) C8 carbonium ion (2); (f) C8 alkoxide (2) g——Isolated state; *——Adsorbed state

ItemBondlengthordistance/nmC2—O1C2—C6C2—C7C6—C7C7—O1Al—O1Al—O2TBA+2-Butene(g)0.1630.535—0.134—0.1910.171TBA+2-Butene*0.1660.3950.4120.134—0.1900.171TS20.2630.3230.3410.134—0.1780.174Triangle-structuredintermediate0.3190.2640.2710.136—0.1770.175C8carboniumion(2)—0.1680.2490.1430.3680.176—C8alkoxide(2)——0.1590.1540.1550.1880.172

g——Isolated state; *——Adsorbed state

表4 TBA與2-丁烯發(fā)生加成反應(yīng)過程中體系內(nèi)部分原子及原子團(tuán)的電荷Table 4 Charge of atoms and atomic groups during addition reaction of TBA and 2-butene

g——Isolated state; *——Adsorbed state

由圖8和表3、表4可以看到，當(dāng)2-丁烯分子處于游離狀態(tài)時(見圖8(a))，2-丁烯分子整體為電中性，即此時2-丁烯分子與TBA之間無相互作用。其中，C6、C7原子所帶電荷分別為-0.07 e和-0.08 e，叔丁基烷氧基中心C原子所帶電荷為+0.09 e，C2—O1鍵長為0.163 nm。隨著2-丁烯分子靠近TBA，二者之間會產(chǎn)生較弱的相互作用，形成吸附態(tài)的2-丁烯分子，其結(jié)構(gòu)如圖8(b)所示。此時，2-丁烯分子整體所帶電荷為+0.03 e，說明2-丁烯分子已經(jīng)向分子篩骨架傳遞了少量的電子。由于二者之間的相互作用，使得C2—O1鍵的鍵長從0.163 nm增加為0.166 nm，TBA發(fā)生反應(yīng)的活性增加，有遠(yuǎn)離活性位點(diǎn)的趨勢。在該吸附態(tài)構(gòu)象中，C6、C7原子電荷均為-0.09 e，而C2原子電荷為+0.10 e，C=C雙鍵的C原子與TBA的中心C原子間存在較強(qiáng)的靜電作用，使得C2原子逐步遠(yuǎn)離了O1原子，靠近異丁烯的C=C雙鍵C原子，形成如圖8(c)所示的叔丁基碳正離子的過渡態(tài)結(jié)構(gòu)(TS2)。此時，C2—O1原子間的距離由0.166 nm增加為0.263 nm，C2—C6、C2—C7原子間距離分別由原來的0.395 nm和0.412 nm變?yōu)?.323 nm和0.341 nm。在這個過程中，2-丁烯分子整體所帶電荷從+0.03 e增加為+0.07 e，叔丁基碳正離子整體所帶電荷為+0.80 e，說明在這個過程中，2-丁烯與TBA之間繼續(xù)發(fā)生電子轉(zhuǎn)移，即2-丁烯C=C 雙鍵的HOMO上的電子傳遞到了TBA的C2原子的LUMO。在TS2中，C6、C7原子所帶電荷均為-0.09 e，而碳正離子中心C原子C2所帶電荷增加為+0.19 e，此時，C2與C6、C7之間的作用力增強(qiáng)，會繼續(xù)靠近，但由于TBA與2-丁烯的加成反應(yīng)是一個由中心C原子為叔C的結(jié)構(gòu)生成了一個中心C原子為仲C的碳正離子結(jié)構(gòu)，即反應(yīng)后碳正離子的穩(wěn)定性降低，因此，叔丁基碳正離子不能與2-丁烯直接生成C8碳正離子，而是首先在C2、C6、C7這3個原子間形成一個“三角形結(jié)構(gòu)”中間產(chǎn)物，如圖8(d)所示。在這個中間產(chǎn)物結(jié)構(gòu)中，C2—O1原子間距離增加為0.319 nm，C2與C6、C7原子間距離分別減小為0.264 nm和0.271 nm，C6=C7雙鍵鍵長由0.134 nm增加為0.136 nm。此時，2-丁烯分子整體電荷從過渡態(tài)結(jié)構(gòu)的+0.07 e增加為+0.25 e，碳正離子所帶電荷從+0.80 e降低為+0.65 e，說明在形成“三角形結(jié)構(gòu)”中間產(chǎn)物的過程中，2-丁烯分子與碳正離子之間發(fā)生了非常明顯的電子轉(zhuǎn)移，電子由2-丁烯分子C=C 雙鍵的HOMO傳遞到TBA中C2原子的LUMO。

要完成加成反應(yīng)，該“三角形結(jié)構(gòu)”中間產(chǎn)物需要在活性位點(diǎn)旋轉(zhuǎn)180°。該過程中，由于C2原子遠(yuǎn)離了O1原子，其所帶正電荷增加，奪電子的能力增強(qiáng)，容易奪取C6=C7雙鍵的電子，從而形成如圖8(e)所示的C8碳正離子(2)，并以C7為中心C原子吸附在活性位點(diǎn)O1上。由于該C8碳正離子的中心碳原子C7所帶電荷為+0.19 e，與分子篩骨架上帶負(fù)電荷的O1原子存在較強(qiáng)的靜電作用，易形成如圖8(f)所示的C8烷氧基團(tuán)(2)，從而完成加成反應(yīng)。

圖9為TBA與2-丁烯加成反應(yīng)的能量示意圖。由圖9可以看到，游離態(tài)2-丁烯分子進(jìn)入分子篩孔道形成吸附態(tài)2-丁烯分子的過程放熱為32.9 kJ/mol，克服29.0 kJ/mol的能壘后，形成叔丁基碳正離子的過渡態(tài)；隨后過渡態(tài)的叔丁基碳正離子與2-丁烯形成“三角形結(jié)構(gòu)”中間產(chǎn)物，體系能量降低了30.5 kJ/mol，“三角形結(jié)構(gòu)”中間產(chǎn)物在活性位點(diǎn)上旋轉(zhuǎn)180°形成C8碳正離子，并以烷氧基團(tuán)的形式吸附在活性位點(diǎn)O1上從而完成加成反應(yīng)。

2.3 TBA與異丁烯和2-丁烯發(fā)生加成反應(yīng)的對比

TBA與異丁烯和2-丁烯2個加成反應(yīng)的異同之處如下：

(1)2個反應(yīng)的過渡態(tài)均為叔丁基碳正離子，但TBA的C—O鍵斷裂后，碳正離子與丁烯之間的作用不同。其中，TS1結(jié)構(gòu)中碳正離子與異丁烯分子整體所帶電荷分別為+0.76 e和+0.10 e，TS2結(jié)構(gòu)中碳正離子與2-丁烯分子整體所帶電荷分別為+0.80 e 和+0.07 e?？梢姡谶^渡態(tài)結(jié)構(gòu)時，異丁烯與碳正離子之間存在更大的電荷轉(zhuǎn)移，即二者間的相互作用力更強(qiáng)。相比之下，2-丁烯與碳正離子之間的電子轉(zhuǎn)移相對明顯，即二者間的相互作用相對較弱。因此，異丁烯與TBA發(fā)生加成反應(yīng)時，反應(yīng)能壘較低，僅為14.9 kJ/mol，小于2-丁烯與TBA發(fā)生加成反應(yīng)時的反應(yīng)能壘29.0 kJ/mol。

(2)發(fā)生反應(yīng)的丁烯異構(gòu)體的結(jié)構(gòu)不同導(dǎo)致加成反應(yīng)路徑的本質(zhì)不同。TBA與異丁烯發(fā)生加成反應(yīng)的反應(yīng)物和產(chǎn)物的中心C原子類型均為叔C，反應(yīng)前后中心C原子類型穩(wěn)定性不變；而TBA與2-丁烯發(fā)生加成反應(yīng)的反應(yīng)物中心碳原子類型為叔C，產(chǎn)物的中心C原子類型為仲C，中心C原子類型穩(wěn)定性降低。這就導(dǎo)致了TBA分別與異丁烯和 2-丁烯發(fā)生加成反應(yīng)的路徑存在本質(zhì)的區(qū)別。其中，異丁烯可直接與過渡態(tài)的碳正離子生成C8碳正離子，而 2-丁烯則首先與過渡態(tài)的碳正離子生成“三角形結(jié)構(gòu)”中間產(chǎn)物，再經(jīng)過旋轉(zhuǎn)生成C8碳正離子，進(jìn)而完成加成反應(yīng)。

圖9 TBA與2-丁烯加成反應(yīng)的能量示意圖Fig.9 Energy diagram for addition reaction of TBA and 2-butene Energies are given relative to separated 2-butane and TBA. g——Isolated state; *——Adsorbed state

3 結(jié) 論

(1)TBA與異丁烯分子發(fā)生加成反應(yīng)是通過叔丁基碳正離子過渡態(tài)進(jìn)行的，由于異丁烯分子與碳正離子間的相互作用力較強(qiáng)，故生成該碳正離子需要克服的能壘較低，僅為14.9 kJ/mol。隨后異丁烯分子C=C雙鍵中帶較多負(fù)電荷的端基C原子與碳正離子中帶正電的中心C原子形成共價鍵，生成C8碳正離子，進(jìn)而形成C8烷氧基團(tuán)。

(2)TBA與2-丁烯分子發(fā)生加成反應(yīng)也是通過叔丁基碳正離子過渡態(tài)進(jìn)行的，但是與異丁烯相比，2-丁烯分子與碳正離子間的相互作用力相對更弱，故生成該碳正離子需要克服的能壘相對較高，為29.0 kJ/mol。而由于TBA與2-丁烯分子發(fā)生加成反應(yīng)后中心C原子穩(wěn)定性降低，故需要首先生成“三角形結(jié)構(gòu)”中間產(chǎn)物，隨后中間產(chǎn)物在活性位點(diǎn)旋轉(zhuǎn)180°，形成C8碳正離子，最終形成C8烷氧基團(tuán)完成加成反應(yīng)。

[1] GB 17930-2016, 中華人民共和國國家標(biāo)準(zhǔn)[S].

[2] 朱玉琴, 司云航, 朱憶寧, 等. 我國車用汽油標(biāo)準(zhǔn)現(xiàn)狀及發(fā)展趨勢[J]. 天然氣化工(C1化學(xué)與化工), 2014, 39(6): 77-81. (ZHU Yuqin, SI Yunhang, ZHU Yining, et al. Current situation and development trend of motor gasoline standards in China[J]. Natural Gas Chemical Industry, 2014, 39(6): 77-81.)

[3] 王吉云, 張清新. 開好烷基化裝置促進(jìn)汽油質(zhì)量升級[J]. 石油石化節(jié)能與減排, 2012, 2(5): 25-28. (WANG Jiyun, ZHANG Qingxin. Energy gasoline quality upgrade by alkylation operation[J]. Conservation and Emission Reduction in Petroleum and Petrochemical Industry, 2012, 2(5):25-28.)

[4] 許友好, 屈錦華, 楊永壇, 等. MIP系列技術(shù)汽油的組成特點(diǎn)及辛烷值分析[J]. 石油煉制與化工, 2009, 40(1): 10-14. (XU Youhao, QU Jinhua, YANG Yongtan, et al. Analysis of octane number and composition characteristics of MIP naphtha[J]. Petroleum Processing Petrochemicals, 2009, 40(1): 10-14.)

[5] 馬愛增. 連續(xù)重整裝置設(shè)計(jì)參數(shù)研究[J]. 石油煉制與化工, 2013,44(6): 64-69. (MA Aizeng. Study of design parameters of CCR reforming units[J]. Petroleum Processing Petrochemicals, 2013,44(6): 64-69.)

[6] 耿英杰. 烷基化生產(chǎn)工藝與技術(shù)[M]. 北京: 中國石化出版社, 1993.

[7] 張劍秋. 降低汽油烯烴含量的催化裂化新材料探索[D]. 北京: 石油化工科學(xué)研究院, 2001.

[8] 林世雄. 石油煉制工程[M]. 第三版. 北京: 石油工業(yè)出版社, 2000: 507-515.

[9] 歐陽健, 鄭明光, 張紹良, 等. DUPONT工藝硫酸烷基化裝置的腐蝕與防護(hù)[J]. 石油化工腐蝕與防護(hù), 2012, 29(6): 31-35. (OUYANG Jian, ZHENG Mingguang, ZHANG Shaoliang, et al. Corrosion of sulfuric acid alkylation unit of DUPONT process and protection[J]. Corrosion & Protection in Petrochemical Industry, 2012, 29(6): 31-35.)

[10] 梁毅. 硫酸法烷基化處理段沖蝕腐蝕分析及處理[J]. 石油化工設(shè)備, 2014, 43(5): 104-107. (LIANG Yi. Erosion corrosion analysis and treatment of sulfuric acid alkylation processing section[J]. Petro-chemical Equipment, 2014, 43(5): 104-107.)

[11] LULY M H, POINTNER B E. Hydrogen fluoride compositions: US, 2011152432A1[P]. 2011-06-23.

[12] FELLER A, LERCHER J A. Chemistry and technology of isobutane/alkene alkylation catalyzed by liquid and solid acids [J]. Advances in Catalysis, 2004, 48: 229-295.

[13] Presidential Green Chemistry Challenge: 2016 Greener Synthetic Pathways Award[EB/OL]. https://www.epa.gov/greenchemistry/presidential-green-chemistry-challenge-2016-greener-synthetic-pathways-award.

[14] ROSENBACH N, DOS SANTOS A P A, FRANCO M, et al. Thetert-butyl cation on zeolite Y: A theoretical and experimental study[J]. Chemical Physics Letters, 2010, 485(1-3): 124-128.

[15] MULLEN G M, JANIK M J. Density functional theory study of alkane-alkoxide hydride transfer in zeolites[J]. ACS Catalysis, 2011, 1(2):105-115.

[16] 任奎, 龍軍, 任強(qiáng), 等. HY分子篩催化異丁烷/丁烯烷基化反應(yīng)中氫負(fù)離子轉(zhuǎn)移過程的模擬[J]. 石油學(xué)報(石油加工), 2017, 33(6): 703-713. (REN Kui, LONG Jun, REN Qiang, et al. A computational study on the hydride transfer process in butene alkylation with isobutane over HY zeolite[J]. Acta Petrolei Sinica (Petroleum Processing Section), 2017, 33(6): 703-713.)

[17] 魯玉瑩, 李永祥, 龍軍, 等. HY催化C4烷基化中異丁烯質(zhì)子化反應(yīng)的分子模擬[J]. 石油學(xué)報(石油加工), 2014, 30(5): 765-771. (LU Yuying, LI Yongxiang, LONG Jun, et al. Molecular simulation of isobutene protonation in C4alkylation catalyzed by HY zeolite[J]. Acta Petrolei Sinica (Petroleum Processing Section), 2014, 30(5): 765-771.)

[18] OLSON D H, DEMPSEY E. The crystal structure of the zeolite hydrogen[J]. Journal of Catalysis, 1969, 13(2): 221-231.

[19] HILL J, FREEMAN C M, DELLEY B. Bridging hydroxyl groups in faujasite: Periodic vs cluster density functional calculations[J]. The Journal of Physical Chemistry A, 1999,103(19): 3772-3777.

[21] SHERWOOD P, DE VRIES A H, GUEST M F, et al. QUASI: A general purpose implementation of the QM/MM approach and its application to problems in catalysis[J]. Journal of Molecular Structure (THEOCHEM), 2003, 632(1-3): 1-28.

[22] DELLEY B. From molecules to solids with the DMol3approach[J]. The Journal of Chemical Physics, 2000, 113(18): 7756-7764.

[23] GALE J D, ROHL A L. The general utility lattice program (GULP)[J]. Molecular Simulation, 2003, 29(5): 291-341.

[24] RAPPE A K, CASEWIT C J, COLWELL K S, et al. UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations[J]. Journal of the American Chemical Society, 1992,114(25): 10024-100035.

[25] RAPPE A K, GODDARD W A. Charge equilibration for molecular dynamics simulations[J]. The Journal of Physical Chemistry, 1991, 95(8): 3358-3363.

[26] BAKOWIES D, THIEL W. Hybrid models for combined quantum mechanical and molecular mechanical approaches[J]. The Journal of Physical Chemistry, 1996, 100(25): 10580-10594.

[27] HENKELMAN G, JNSSON H. Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points[J]. Journal of Chemical Physics, 2000, 113(22): 9978-9985.

[28] SUN Y, YANG J, ZHAO L, et al. A two-layer ONIOM study on initial reactions of catalytic cracking of 1-butene to produce propene and ethene over HZSM-5 and HFAU zeolites[J]. The Journal of Physical Chemistry C, 2010, 114(13): 5975-5984.