劉文霞,王金道,許有偉
(1 中國石油蘭州化工研究中心,甘肅蘭州 730060;2 中國石油蘭州石化公司,甘肅蘭州 730060)
目前,二氧化硅制備技術(shù)處于粒徑小、比表面積大的方向。多峰介孔二氧化硅是無定型二氧化硅重要的成員,而方石英內(nèi)部具有一定無定型二氧化硅形狀,因此選擇方石英晶體作為吸附質(zhì)來研究二氧化硅對MAO 的吸附。晶型二氧化硅中的氧分布與無定型二氧化硅中的羥基類型形成映襯關(guān)系,基于此,了解無定型二氧化硅中羥基對MAO 形成的配位。
MAO 是鋁和氧原子隨機(jī)排列且自由價(jià)被甲基取代基飽和的化合物。由于合成條件、分離手段、所用溶劑等原因,其晶體結(jié)構(gòu)的定義是非常困難的工作。已經(jīng)提出了許多結(jié)構(gòu),例如聚合物線性鏈、二維結(jié)構(gòu)、環(huán)和三維籠結(jié)構(gòu)[1-5],如圖1 所示。
分子在表面上的吸附是表面催化反應(yīng)的第一步,產(chǎn)生的吸附行為對于理解這樣的反應(yīng)至關(guān)重要。量子化學(xué)計(jì)算是模擬化學(xué)反應(yīng)細(xì)節(jié)的強(qiáng)大工具,除了其他信息外,它們還為各種構(gòu)象提供分子優(yōu)化及能量計(jì)算,有助于更好地理解界面相互作用[6-8]。
圖1 MAO 的結(jié)構(gòu)形式Fig.1 The structure of MAO
密度泛函數(shù)理論(DFT)計(jì)算已被廣泛用于模擬表面結(jié)構(gòu)的反應(yīng)性。文中采用DFT 模擬并研究MAO 在方石英表面上的吸附,通過了解不同晶面對MAO 的吸附情況,以此推斷無定型二氧化硅中什么類型的羥基對MAO 作用最強(qiáng),從而通過改性無定型二氧化硅的制備工藝,不斷提升此種類型的羥基在硅膠內(nèi)的分布,以此達(dá)到對MAO的充分利用,降低茂金屬催化劑的制備成本。
方石英晶體是立方晶系,空間群具有a=b=c=7.16?的空間特性。二氧化硅晶體結(jié)構(gòu)采用Material Studio 中的CASTEP 模塊進(jìn)行優(yōu)化[9],在進(jìn)行結(jié)構(gòu)優(yōu)化過程中,選擇由Perdew、Burke 和Ernzerhof[10]提出的包含Becke交換和相關(guān)表達(dá)式的廣義梯度近似(GGA)函數(shù)。通過將方石英晶體的表面切割成給定的結(jié)晶學(xué)(hkl)平面,在垂直于SiO2板的方向上施加10.0 埃的真空厚度以避免兩層之間的相互作用。在所有情況下,互易空間用Monkhorst-Pack 特殊K 點(diǎn)的方案表示[11],其中網(wǎng)格參數(shù)設(shè)定為3×4×1,將分?jǐn)?shù)頂部添加到1.5,并且能量切斷的模型被選擇定制,域值從260.0eV 改變?yōu)?00.0eV。文本中使用的單元格顯示樣式是Default。
從圖2 可以看出,低能量部分主要是O 2s 原子軌道,與Si 3s 和3p 混合,在0~11eV 之間,主要是價(jià)帶,主要是O 2p,與一些Si 3s 和3p 混合;導(dǎo)帶部分主要是Si軌道。
圖2 二氧化硅中Si 和O 的局部態(tài)密度(PDOS)Fig.2 Local state density of Si and O in SiO2
方石英二氧化硅形成的帶隙為4.943eV,能帶分布及形成的PDOS 如圖3 所示。
圖3 二氧化硅產(chǎn)生的能帶分布及PDOSFig.3 Energy band distribution and PDOS produced by silica
二氧化硅不同晶面及MAO 吸附在二氧化硅不同晶面產(chǎn)生的電子結(jié)構(gòu)進(jìn)行計(jì)算,結(jié)果如圖4 所示。
圖4 不同二氧化硅晶面及MAO 吸附于二氧化硅不同晶面的PDOS 分布Fig.4 PDOS distribution of different SiO2 crystal planes and MAO adsorbed on different SiO2 crystal planes
通過比較SiO2晶體的PDOS 與三個(gè)SiO2晶面的PDOS,發(fā)現(xiàn)當(dāng)形成不同的晶面時(shí),SiO2晶體的帶寬變寬。頻帶越寬,頻帶圖上的波動(dòng)越大,這表明頻帶中電子的有效質(zhì)量越小,非局域性越大,并且頻帶的原子軌道擴(kuò)展越強(qiáng)。然而,由SiO2(100)和SiO2(110)形成的帶明顯寬于由SiO2(111)形成的帶寬,這表明SiO2(111)表面上的電子定位非常強(qiáng)并且有效質(zhì)量比較大。
當(dāng)MAO 吸附于不同晶面時(shí),MAO 中的Al 原子的p軌道存在空位,與SiO2中的O 形成配位,即產(chǎn)生吸附作用。從圖4 中可以發(fā)現(xiàn),當(dāng)MAO 吸附于不同晶面時(shí),在導(dǎo)帶處均出現(xiàn)了一個(gè)明顯的p 軌道,這與Al 原子的3p 軌道的特征峰值正好吻合,Al 原子的軌道分布如圖5所示。然而MAO 在SiO2(111)形成的導(dǎo)帶峰寬明顯要窄于MAO 在其他兩面晶面,說明MAO 在SiO2(111)形成的電子局域性非常強(qiáng),有效質(zhì)量相對較大;而MAO 在SiO2(100)面形成導(dǎo)帶寬度又明顯寬于MAO 在SiO2(110)面形成的導(dǎo)帶,說明此帶中電子有效質(zhì)量越小、非局域的程度越大。而根據(jù)計(jì)算獲得,MAO 在SiO2(100)面吸附形成的帶隙為0.016eV,在Si(110)面吸附形成的帶隙為0.086eV,而在SiO2(111)面形成的帶隙為0.180eV,帶隙越寬,說明Al 和O 兩個(gè)原子的成鍵越強(qiáng),這也說明MAO 在SiO2(111)吸附越穩(wěn)定。
圖5 鋁原子的PDOS 分布Fig.5 PDOS distribution of aluminum atom
圖6 二氧化硅與MAO 作用機(jī)理Fig.6 Mechanism of action of silica and MAO
無定形二氧化硅是烯烴聚合,尤其是茂金屬催化體系最重要的載體,它是SiO2·nH2O 無定型物,表面存在有自由羥基、氫鍵羥基、締合羥基和吸附水等不同的基團(tuán)。在進(jìn)行MAO 的吸附實(shí)驗(yàn)前,需要對二氧化硅進(jìn)行熱活化處理,以使二氧化硅表面的各種基團(tuán)具有適當(dāng)?shù)谋壤⑶曳植季鶆颉?/p>
二氧化硅與MAO 作用機(jī)理如圖6 所示。
無定形二氧化硅表面羥基類型主要有四種,如圖7所示。
圖7 無定形二氧化硅表面羥基類型Fig.7 Amorphous silica surface hydroxyl type
從映襯關(guān)系可以發(fā)現(xiàn),自由羥基上的O 原子與SiO2(111)晶面上的原子分布接近,而橋連羥基上的O 原子與SiO2(110)晶面上的O 原子分布接近,硅氧烷橋鍵中的O 原子與SiO2(100)晶面上的O 原子分布接近。自由羥基對MAO 的吸附最強(qiáng),因此,在制備硅膠過程中,盡可能讓自由羥基分布于硅膠表面及孔道內(nèi)部,以達(dá)到對MAO 最大的吸附。映襯關(guān)系如圖8 所示。
圖8 二氧化硅不同晶面與無定形羥基類型映襯關(guān)系Fig.8 The relationship between different crystal faces of silica and amorphous hydroxyl type
無定形二氧化硅表面羥基數(shù)量隨著溫度的升高而不斷降低,結(jié)果見表1。
表1 二氧化硅表面羥基的總數(shù)量隨溫度關(guān)系Table 1 The total number of hydroxyl groups on the surface of silica varies with temperature
Pascal 認(rèn)為[12],硅膠表面上自由羥基和橋連羥基的相對分布強(qiáng)烈地取決于樣品的預(yù)處理溫度和孔半徑。較小的孔半徑有利于橋連羥基與游離羥基的縮合,橋連羥基的減少與縮合反應(yīng)一致,反應(yīng)如圖9 所示。
圖9 縮合反應(yīng)式Fig.9 Condensation reaction
從縮合反應(yīng)式可知,自由羥基數(shù)的減少是由橋連羥基的隨機(jī)縮合引起的。
(1)MAO 在SiO2(111)面形成的帶隙值最大,帶隙越寬,說明Al 和O 兩個(gè)原子的成鍵越強(qiáng),產(chǎn)生的吸附越穩(wěn)定。
(2)二氧化硅晶面上的O 原子分布與硅膠的羥基類型存在映襯關(guān)系,自由羥基對應(yīng)SiO2(111)晶面,橋連羥基對應(yīng)SiO2(110)晶面。
(3)橋連羥基的減少與縮合反應(yīng)一致,而自由羥基的減少是由橋連羥基隨機(jī)縮合引起的。