李對春，邢 斌，劉紅艷，劉 平，王寶俊，李瑞豐

(1.太原理工大學 化學化工學院，山西 太原 030024;2.大同大學 化學與環(huán)境工程學院，山西 大同 037009;3.中國科學院 山西煤炭化學研究所,山西 太原 030001;4.太原理工大學 煤科學與技術教育部和山西省重點實驗室,山西 太原 030024;5.太原理工大學 能源化工催化研究中心，山西 太原 030024)

0 引 言

低碳烯烴中乙烯、丙烯是現(xiàn)代化學工業(yè)的基本有機原料，隨著化工產(chǎn)業(yè)的飛速發(fā)展，其需求量將越來越大。制取乙烯、丙烯的傳統(tǒng)路線是通過石腦油裂解生產(chǎn)，但由于石油是不可再生資源，儲量有限，價格起伏較大，尋求替代傳統(tǒng)的非石油路線生產(chǎn)低碳烯烴勢在必行[1-4]。甲醇制烯烴(methanol-to-olefins，MTO)是替代傳統(tǒng)路線的新工藝，是環(huán)境催化領域的研究熱點[2]。目前，甲醇的制備路線已成熟完善，不僅可通過生物質(zhì)發(fā)酵制得，還可以通過煤、天然氣經(jīng)由合成氣(CO+H2)制得。而MTO是該替代路線的關鍵環(huán)節(jié)。

目前，共提出20多種甲醇制烯烴反應機理，但至今沒有得出與反應完全符合的機理[2]，其中，最具代表性的為直接反應機理和間接反應機理。根據(jù)反應過程中產(chǎn)生的特殊中間體，直接反應機理又被分為卡賓機理[1,5]、氧鎓葉立德機理[6-7]、碳正離子機理[8]、自由基機理[9]等。Dahl等[10-12]提出了間接反應機理(烴池機理)，即限域在分子篩孔道內(nèi)的烴池物種與分子篩作為共催化劑，通過與甲醇發(fā)生連續(xù)的甲基化，使碳鏈增長，經(jīng)低碳烯烴消除反應生成產(chǎn)物烯烴(乙烯和丙烯)，同時，烴池物種和分子篩的B酸性位點再生，完成一個催化循環(huán)，烴池機理得到了研究者們的普遍認可。以烯烴為基礎的循環(huán)起源于30年前由Dessau[13-14]提出的烯烴甲基化和裂化路線。王傳明等[15]通過密度泛函理論證明了SAPO-34分子篩上發(fā)生的MTO反應除了聚甲基苯外，烯烴本身也是主要的烴池物種，即在MTO催化反應中同時存在芳烴路線和烯烴路線。

在MTO反應中，SAPO-34分子篩表現(xiàn)出優(yōu)異的催化性能。SAPO-34分子篩具有8元環(huán)構(gòu)成的橢球形籠和三維孔道結(jié)構(gòu)，這種分子篩的小孔結(jié)構(gòu)、中等酸性強度及良好的水熱穩(wěn)定性使其具有良好的MTO催化性能[16-17]。不同的分子篩具有不同的骨架拓撲結(jié)構(gòu)及酸性質(zhì)，包括酸性位點的分布、酸性位點的密度及強度等，這些因素決定了其催化活性和產(chǎn)物的選擇性。其中，分子篩催化劑酸性位點的強度是影響MTO反應催化性能的重要因素[18-20]，但通過試驗手段無法明確單個因素對MTO反應催化性能的影響。因此，本文采用色散力矯正的密度泛函理論方法，構(gòu)建SAPO-34分子篩模型，以2，3-二甲基-2-丁烯(iso-C6)為烴池物種[21]，選擇不同金屬雜原子Fe、Co和Ni同晶取代Al原子對SAPO-34分子篩催化劑進行改性，考察了不同金屬原子Fe、Co和Ni改性的SAPO-34 分子篩中B酸位點強度的變化，并分析改性SAPO-34分子篩對MTO反應的催化活性和產(chǎn)物選擇性的影響。

1 計算方法與模型

1.1 計算方法

采用VASP軟件包的密度泛函理論方法(DFT)進行結(jié)構(gòu)優(yōu)化和總能量計算[22-24]，電子的交換關聯(lián)勢采用廣義梯度近似(GGA)下的PBE泛函來描述[25]，用PAW方法描述離子實，價電子用平面波基組擴展[26-27]，同時考慮了庫倫相互作用的色散效應[28]。設定截斷能為400 eV，通過CI-NEB方法尋找過渡態(tài)[29]，并通過頻率分析確認過渡態(tài)，在計算過程中，還考慮了零點能(ZPE)校正[30-31]。對于SAPO-34分子篩及Fe、Co和Ni改性的SAPO-34分子篩，布里淵帶使用1×1×1的Monkhost-Pack篩狀的k點空間取樣[32]。計算過程中基于自旋極化方法考慮了Fe、Co和Ni原子的磁性，結(jié)構(gòu)優(yōu)化時晶胞參數(shù)保持不變，直至每個原子上的力小于0.02 eV/Angst。

1.2 計算模型

SAPO-34分子篩具有菱沸石(CHA)骨架結(jié)構(gòu)單元[33]，屬于小孔分子篩[34]。CHA晶胞單元中含有36個T原子(1個Si原子，18個Al原子，17個P原子)和72個O原子。由于Si、P原子價態(tài)不同，需要用額外的一個質(zhì)子(H+)補償晶胞中產(chǎn)生的一個負電荷來保持整體結(jié)構(gòu)的電中性，該位點(Si-O(H)-Al)即為Br?nsted酸性中心[35]。SAPO-34的周期性結(jié)構(gòu)如圖1所示，優(yōu)化的晶胞參數(shù)為a=b=13.85×10-10m,c=15.08×10-10m，與國際分子篩協(xié)會數(shù)據(jù)庫中的晶胞參數(shù)(a=b=13.68×10-10m,c=14.77×10-10m)基本相符[36]。

金屬雜原子Fe、Co和Ni通過同晶取代Al原子的方式進入分子篩骨架，得到MeAPO-34(Me=Fe、Co、Ni)催化劑。由于Fe原子與Al原子的價態(tài)相同，所以Fe(+3)原子替換Al(+3)原子時不會產(chǎn)生電荷差，而Co(+2)原子和Ni(+2)原子與Al(+3)原子的價態(tài)不同，所以替換Al原子時，同樣用質(zhì)子(H+)補償晶胞中產(chǎn)生一個負電荷。FeAPO-34，CoAPO-34和NiAPO-34分子篩優(yōu)化后的穩(wěn)定結(jié)構(gòu)如圖2所示，計算模型中所有的原子保持弛豫。

2 結(jié)果與討論

2.1 金屬雜原子引入SAPO-34分子篩后對酸性強度的影響

分子篩催化劑的B酸位點的強度是影響MTO反應催化性能的重要因素[19-20,37]。B酸定義為任何化學物種(分子或離子)能夠失去或貢獻出氫離子(質(zhì)子)。去質(zhì)子化能是評價酸性強度的標準之一，根據(jù)B酸的定義，本文采用去質(zhì)子化能EDPE(從Si-O(H)-Al位點上失去質(zhì)子(H+)所需要的能)描述B酸性位點的酸性強度[38-40]。計算公式[39-40]為

圖2 3種金屬原子Fe、Co、Ni改性的SAPO-34結(jié)構(gòu)示意Fig.2 Configuration of FeAPO-34,CoAPO-34 and NiAPO-34

EDPE=EMeAPO-34-+EH+-EMeAPO-34

其中，EMeAPO-34為MeAPO-34分子篩的結(jié)構(gòu)能量;EMeAPO-34-為MeAPO-34結(jié)構(gòu)中失去一個質(zhì)子H+后的結(jié)構(gòu)能量;EH+為一個質(zhì)子H+的能量。EDPE越小，對應的B酸位點越容易貢獻出H+質(zhì)子,B酸位點的酸性越強。EDPE的計算結(jié)果見表1，其與Jones.等[41]的計算結(jié)果一致。B酸位點周圍有4個O位點，為了比較不同MeAPO-34結(jié)構(gòu)中B酸的酸性強度，分別計算4個O位點的EDPE值和EDPE平均值(Eav-DPE)(表1)。

Eav-DPE=(EDPE-1+EDPE-2+EDPE-3+EDPE-4)/4

從表1可知，當3種金屬雜原子Fe、Co和Ni引入SAPO-34結(jié)構(gòu)，NiAPO-34結(jié)構(gòu)的Eav-DPE值比未改性SAPO-34結(jié)構(gòu)的Eav-DPE值大，說明NiAPO-34結(jié)構(gòu)中的B酸位點的強度較弱。但FeAPO-34和CoAPO-34結(jié)構(gòu)的Eav-DPE值均比未改性SAPO-34結(jié)構(gòu)的Eav-DPE值小，說明FeAPO-34和CoAPO-34結(jié)構(gòu)中的B酸位點的強度較強。因此，改性前后SAPO-34結(jié)構(gòu)的B酸位點的酸性強度順序為FeAPO-34>CoAPO-34>SAPO-34>NiAPO-34，表明金屬雜原子引入SAPO-34結(jié)構(gòu)中能改變分子篩的B酸位點的酸性強度。

表1SAPO-34結(jié)構(gòu)及MeAPO-34結(jié)構(gòu)中的EDPE和Eav-DPE
Table1CalculatedEDPEandEav-DPEofSAPO-34andMeAPO-34structure

結(jié)構(gòu)位點nEDPE/(kJ·mol-1)Eav-DPE/(kJ·mol-1)11 173.5SAPO-34-Si-OHn-Al21 166.71 168.131 158.541 173.711 037.4FeAPO-34-Si-OHn-Al21 009.01 017.431 006.341 016.911 092.4CoAPO-34-Si-OHn-Al21 082.41 073.431 039.141 079.811 191.9NiAPO-34-Si-OHn-Al21 159.81 174.031 174.841 169.6

2.2 金屬雜原子引入SAPO-34分子篩后對MTO反應性能的影響

在MeAPO-34穩(wěn)定結(jié)構(gòu)的基礎上，基于有機活性中心為iso-C6分子烯烴循環(huán)的烴池機理，即iso-C6分子與分子篩作為共催化劑，與甲醇發(fā)生連續(xù)的甲基化，使碳鏈增長，然后發(fā)生低碳烯烴消除反應，生成產(chǎn)物乙烯和丙烯，同時，烴池物種iso-C6分子和分子篩的B酸位點再生，形成一個完整的催化循環(huán)。反應路線如圖3所示。

反應步驟可分為:

1)甲基化反應

(M1)

(CH3)3C+—C(CH3)(CH2CH3)+H2O+Z-

(M2)

(M3)

圖3 MTO反應的烯烴循環(huán)路徑Fig.3 Reaction pathway of olefin-based for MTO conversion

2)過渡反應

(CH3)3C—C(CH3)(CH2)+H2O+HZ

(D1)

(CH3)3C—C(CH3)(CHCH3)+H2O+HZ

(D2)

H2O+Z-

(C1)

C+(CH3)2(C3H7)+H2O+Z-

(C2)

3)烯烴脫除反應

(E1)

(E2)

各個反應步驟的活化能(Ea)和反應熱(ΔH)見表2。本文選擇甲基化反應步驟(M1、M2、M3)和低碳烯烴(乙烯和丙烯)消除步驟(E1、E2)來評價分子篩的酸性強度對MTO反應催化性能的影響。

表2MeAPO-34結(jié)構(gòu)中不同反應步驟的活化能
(Ea)和反應熱(ΔH)
Table2CalculatedEaandΔHofdifferentreactionstepsoverMeAPO-34structurekJ/mol

MTO反應的第1步是甲基化反應(M1)，其反應物、產(chǎn)物和過渡態(tài)(S1、S2、TS1)結(jié)構(gòu)如圖4所示。在FeAPO-34，CoAPO-34和NiAPO-34結(jié)構(gòu)上發(fā)生M1反應的反應熱分別為26.1、-13.5和-28.9 kJ/mol，需要克服的能壘分別為104.2、99.4和112.9 kJ/mol。在NiAPO-34結(jié)構(gòu)中，M1反應放熱最多，需要克服的能壘最高，說明Ni的存在不利于CH3OH的活化;但在FeAPO-34和CoAPO-34結(jié)構(gòu)中，M1反應需要克服的能壘較接近，說明FeAPO-34和CoAPO-34對CH3OH活化的催化活性相近。

第2步甲基化反應(M2)的反應物、產(chǎn)物和過渡態(tài)(S4、S5、TS3)結(jié)構(gòu)與M1反應的S1、S2、TS1結(jié)構(gòu)相似，其鍵長見表3。在FeAPO-34、CoAPO-34和NiAPO-34結(jié)構(gòu)上發(fā)生M2反應的反應熱分別為-25.1、34.7和-39.6 kJ/mol，需要克服的能壘分別為98.4、100.3和148.6 kJ/mol。NiAPO-34結(jié)構(gòu)仍不利于CH3OH的活化，F(xiàn)eAPO-34和CoAPO-34結(jié)構(gòu)對CH3OH活化的催化活性相近。

乙烯消除反應步驟(E1)的反應物、產(chǎn)物和過渡態(tài)(S6，S7，TS5)結(jié)構(gòu)如圖5所示。在FeAPO-34、CoAPO-34和NiAPO-34結(jié)構(gòu)上發(fā)生乙烯消除反應的反應熱分別為-12.5、5.8和-54.0 kJ/mol，需要克服的能壘分別為86.8、130.3和83.0 kJ/mol。結(jié)果表明:FeAPO-34和NiAPO-34表現(xiàn)出較好的催化活性，而CoAPO-34的催化活性較差，與Kang等[42]通過快速水熱晶化法合成的Fe、Co、Ni三種金屬改性的SAPO-34分子篩上對乙烯選擇性的順序一致。

FeAPO-34、CoAPO-34和NiAPO-34的生成產(chǎn)物丙烯的反應路徑中，甲基化反應(M3)的反應熱分別為32.8、-6.8和-56.0 kJ/mol，需要克服的能壘分別為104.2、93.6和138.0 kJ/mol，與M1、M2反應的計算結(jié)果相似，即NiAPO-34結(jié)構(gòu)對甲基化反應的催化活性最差，F(xiàn)eAPO-34和CoAPO-34結(jié)構(gòu)對甲基化反應的催化活性相近。M3甲基化反應步驟中的反應物、產(chǎn)物和過渡態(tài)(S9，S10，TS7)結(jié)構(gòu)的鍵長見表3。

丙烯消除反應(E2)的反應物、產(chǎn)物和過渡態(tài)(S11，S12，TS9)結(jié)構(gòu)與E1反應中的S6，S7，TS5結(jié)構(gòu)相似，其鍵長見表4。在FeAPO-34，CoAPO-34和NiAPO-34結(jié)構(gòu)上發(fā)生丙烯消除反應的反應熱分別為-9.6、-34.7和46.3 kJ/mol，需要克服的能壘分別為38.6、41.5和131.2 kJ/mol，可知，對于丙烯消除反應E2，F(xiàn)eAPO-34和CoAPO-34結(jié)構(gòu)催化活性明顯比NiAPO-34結(jié)構(gòu)的催化活性好。

圖4 不同MeAPO-34結(jié)構(gòu)中M1的反應物(S1)、產(chǎn)物(S2)和過渡態(tài)(TS1)結(jié)構(gòu)示意Fig.4 Configurations of reactant (S1),product (S2) and transition state (TS1) of M1 reaction step over different MeAPO-34 structures

反應步驟MeAPO-34結(jié)構(gòu)結(jié)構(gòu)鍵長/10-10 mC1—C2C2—O1H1—O1H1—O2FeAPO-343.361.451.431.06CoAPO-34反應物(S4)3.401.461.271.17NiAPO-343.621.451.601.03FeAPO-342.242.041.021.40M2CoAPO-34過渡態(tài)(TS3)2.251.991.011.76NiAPO-342.251.991.011.76FeAPO-341.463.101.041.62CoAPO-34產(chǎn)物(S5)1.583.970.981.97NiAPO-341.583.970.981.97FeAPO-343.531.461.421.08CoAPO-34反應物(S9)3.781.451.341.12NiAPO-343.261.462.241.04FeAPO-342.322.041.021.61M3CoAPO-34過渡態(tài)(TS7)2.351.991.011.71NiAPO-342.192.191.011.45FeAPO-341.604.010.983.00CoAPO-34產(chǎn)物(S10)1.623.200.991.88NiAPO-341.594.371.002.09

注:C1為有機分子的C原子;C2為CH3OH分子中的C原子;O1為CH3OH分子中的O原子;H1為分子篩骨架上B酸位點的H原子;O2為分子篩骨架上B酸位點的O原子。

表4不同MeAPO-34結(jié)構(gòu)上丙烯消除反應E2的反應物，產(chǎn)物和過渡態(tài)結(jié)構(gòu)中的鍵長
Table4Correspondingbondlengthofreactant,productandtransitionstateofE2reactionstepsoverdifferentMeAPO-34structures

MeAPO-34結(jié)構(gòu)結(jié)構(gòu)鍵長/10-10 mC1—C2'C3—H1'H1'—O1'H2—O1'H2—O2FeAPO-341.581.102.770.991.87CoAPO-34反應物(S11)1.571.102.510.982.04NiAPO-341.561.102.520.992.25FeAPO-342.891.291.381.051.52CoAPO-34過渡態(tài)(TS9)3.211.341.341.011.73NiAPO-342.901.261.431.012.06FeAPO-343.233.020.981.551.03CoAPO-34產(chǎn)物(S12)3.251.951.011.261.18NiAPO-343.591.951.011.261.18

注:C2′為丙烯分子中的C原子;C3為丙烯分子中的C原子;H1′為丙烯分子中的H原子;O1′為H2O分子中的O原子;H2為H2O分子中的H原子。

圖5 不同MeAPO-34結(jié)構(gòu)中乙烯消除反應E1的反應物(S6)、產(chǎn)物(S7)和過渡態(tài)(TS5)結(jié)構(gòu)示意Fig.5 Configurations of reactant (S6),product (S7) and transition state (TS5) of E1 reaction step over different MeAPO-34 structures

由表4可知，在3種不同的MeAPO-34結(jié)構(gòu)中，F(xiàn)eAPO-34和CoAPO-34上丙烯消除反應E2需要克服的反應能壘比乙烯消除反應E1要低，而NiAPO-34上丙烯消除反應E2需要克服的反應能壘比乙烯消除反應E1高，說明在本文的烯烴循環(huán)反應路徑中，F(xiàn)eAPO-34和CoAPO-34分子篩對丙烯的選擇性好，NiAPO-34分子篩對乙烯的選擇性好。

3 結(jié) 論

1)通過比較不同MeAPO-34結(jié)構(gòu)中Si-O(H)-Al酸性位點的去質(zhì)子化能可知，MeAPO-34分子篩的B酸位的酸性強度順序為FeAPO-34>CoAPO-34>SAPO-34>NiAPO-34，不同金屬雜原子的引入對SAPO-34分子篩B酸位點的酸性強度起調(diào)節(jié)作用。

2)以iso-C6分子為烴池物種，研究烯烴循環(huán)反應路徑對MTO反應中的甲基化反應步驟(M1、M2、M3)和烯烴消除反應(E1、E2)的催化活性，結(jié)果表明，F(xiàn)e、Co和Ni改性的SAPO-34分子篩、 M1、M2、M3反應的催化活性順序為FeAPO-34≈CoAPO-34>NiAPO-34;對乙烯生成反應(E1)的催化活性順序為FeAPO-34≈NiAPO-34>CoAPO-34;對丙烯生成反應(E2)的催化活性順序為FeAPO-34≈CoAPO-34>NiAPO-34。

3)NiAPO-34結(jié)構(gòu)中B酸位的酸性強度比FeAPO-34和CoAPO-34結(jié)構(gòu)中的弱，表明強酸性位點對甲基化反應(M1、M2、M3)和丙烯消除反應E2的催化活性好，而弱酸性位點對乙烯消除反應E1的催化活性好。

4)對于烯烴循環(huán)反應路徑而言，F(xiàn)eAPO-34和CoAPO-34分子篩對丙烯的選擇性好，NiAPO-34分子篩對乙烯的選擇性好，即強酸位點有利于丙烯產(chǎn)物的生成，而弱酸位點有利于乙烯產(chǎn)物的生成。計算結(jié)果為MTO催化反應中的分子篩催化劑改性提供理論依據(jù)，從而制備出產(chǎn)物選擇性和催化活性更好的分子篩催化劑。