宮亞玲，龐 可，任瑞鵬

(太原理工大學(xué) 煤科學(xué)與技術(shù)教育部與山西省共建重點(diǎn)實(shí)驗(yàn)室，太原 030024)

甲醇裂解反應(yīng)是其高效利用不可避免的重要步驟，該反應(yīng)過(guò)程中催化劑扮演了重要的角色，但是催化劑的積碳和CO中毒問(wèn)題制約了甲醇高效利用。ZSM-5分子篩被認(rèn)為是甲醇裂解反應(yīng)的有效催化劑，但存在石墨碳覆蓋活性位點(diǎn)和多環(huán)芳烴堵塞孔道等積碳問(wèn)題，影響了分子篩催化劑的工業(yè)化應(yīng)用[1-3]。為了提高分子篩催化劑在甲醇裂解過(guò)程中的催化性能，研究者采用金屬對(duì)分子篩催化劑進(jìn)行改性，金屬改性后多環(huán)芳烴的形成明顯減少[4]，但是仍無(wú)法避免酸性位被覆蓋。SAPO-34分子篩也是甲醇裂解的有效催化劑，目前影響其工業(yè)化應(yīng)用的原因主要是積碳失活，積碳主要來(lái)源于反應(yīng)過(guò)程中碳的聚集[5-6]。此外貴金屬鉑也是甲醇裂解中出色的單金屬催化劑，甲醇在金屬鉑上裂解的主路徑為CH3OH→CH3O→CH2O→CHO→ CO，反應(yīng)過(guò)程中不會(huì)造成碳的積聚，但鉑的活性催化位點(diǎn)容易被CO覆蓋造成中毒[7-9]。綜上所述，可以發(fā)現(xiàn)：催化劑積碳形成的原因主要為甲醇裂解過(guò)程中生成的碳聚集覆蓋在催化活性位點(diǎn)；CO中毒的原因則為CO因其吸附能過(guò)大難以從催化劑表面解吸，沉積在表面引起中毒。因此，甲醇高效利用的前提是開發(fā)出活性中心與反應(yīng)介質(zhì)隔離且CO容易脫附的穩(wěn)定催化劑。

研究表明，限制在碳納米管通道內(nèi)的金屬納米顆粒與沉積在碳納米管管壁的相同金屬具有不同的催化活性以及穩(wěn)定性[10-14]。CUI et al[15]發(fā)現(xiàn)，將碳化鉬納米顆粒封裝在單壁碳納米管中用于析氫反應(yīng)，封裝于碳納米管內(nèi)的碳化鉬比分散在碳納米管外壁上的碳化鉬納米顆粒表現(xiàn)出更高的催化活性和穩(wěn)定性。ZHANG et al[16]將Pt簇封裝在碳納米管通道中用于甲苯氧化反應(yīng)，研究發(fā)現(xiàn)：碳納米管的封裝作用可以保護(hù)Pt簇免受氧氣的影響；與沉積在碳納米管外壁和炭黑載體表面的Pt簇相比，碳納米管封裝的Pt簇具有更高的活性和穩(wěn)定性；碳納米管缺乏強(qiáng)酸位點(diǎn)，避免了反應(yīng)過(guò)程中的碳沉積和孔道堵塞。顯然，碳納米管封裝納米金屬顆粒(M@CNT)做催化劑可以將活性中心與反應(yīng)介質(zhì)隔離，提高金屬催化劑的活性和穩(wěn)定性?；谔技{米管封裝金屬的這一特性，M@CNTs催化劑可能是甲醇裂解的優(yōu)選催化劑。

有關(guān)甲醇裂解的金屬催化劑種類有很多，其中Cu基催化劑被認(rèn)為是甲醇裂解的有效催化劑[17-19]。MEHMOOD et al[20]探索了甲醇在Cu4簇上的裂解，發(fā)現(xiàn)甲醇在Cu4簇上首先通過(guò)C—H鍵斷裂生成CH2OH，隨后通過(guò)C—H鍵斷裂最終生成CO，CO的吸附能為-1.71 eV，其易從催化劑表面解吸，反應(yīng)過(guò)程中無(wú)CO中毒問(wèn)題。Garcia-Muelas et al[17]研究了甲醇在Cu(111)、Ru(0001)、Pt(111)和Pd(111)表面的分解，其中Cu(111)優(yōu)先生成CH2O，而Pt(111)、Pd(111)和Ru(0001)最初生成大量的CO，可能會(huì)毒化表面。CO在這些金屬催化劑表面的吸附能分別為-0.25 eV、-2.78 eV、-3.21 eV和-2.04 eV，可見(jiàn)與其他金屬相比，Cu基催化劑具有優(yōu)異的CO耐受性。CO在Cu表面的低吸附能是催化劑抗中毒的主要原因。

本文運(yùn)用密度泛函理論研究了甲醇在Cu@CNT上的裂解機(jī)理，并探索了該催化劑的抗積碳性及CO耐受性。

1 計(jì)算模型及方法

1.1 計(jì)算模型

關(guān)于碳納米管的選擇，XIAO et al[21]發(fā)現(xiàn)封裝金屬的碳納米管管徑對(duì)催化活性的影響呈火山型曲線關(guān)系，內(nèi)徑約1.0 nm的碳納米管應(yīng)為最佳碳納米管。本文構(gòu)建了直徑為0.93 nm的(12，0)碳納米管，并將Cu原子封裝在碳納米管管道內(nèi)。優(yōu)化后的構(gòu)型如圖1所示，封裝的Cu原子與碳納米管管壁的C原子成鍵，C—Cu鍵長(zhǎng)為0.208 7 nm.催化劑表面有四種不同的吸附位點(diǎn)，分別為與Cu原子成鍵的C原子的頂位(T1)、T1鄰位C原子的頂位(T2)、橋位B和穴位H。

圖1 Cu@CNT催化劑最優(yōu)結(jié)構(gòu)的側(cè)視圖和俯視圖Fig.1 Side view and top view of the optimal structure of Cu@CNT catalyst

1.2 計(jì)算方法

本研究中所涉及到的所有計(jì)算均采用基于密度泛函理論的VASP軟件包(Vienna abinitio simulation package)進(jìn)行[22]。采用綴加投影波函數(shù)(PAW)來(lái)描述電子-離子的相互作用[23]，并使用廣義梯度函數(shù)(GGA)PEB計(jì)算交換相關(guān)能[24]。平面波截?cái)嗄芎筒祭餃Y區(qū)K點(diǎn)分別設(shè)置為270 eV和3×2×2.通過(guò)計(jì)算甲醇及其裂解過(guò)程中涉及到的中間體的吸附能來(lái)確定各物種的最佳吸附位點(diǎn)。吸附能Eads使用以下定義式計(jì)算：

Eads=ECu@CNT+adsorbate-ECu@CNT-Eadsorbate.

(1)

式中：ECu@CNT+adsorbate是催化劑Cu@CNT和吸附物的總能量，ECu@CNT是Cu@CNT催化劑的能量，Eadsorbate是吸附物的能量。通常Eads為負(fù)值時(shí)表示該吸附過(guò)程為放熱吸附或穩(wěn)定吸附，且Eads絕對(duì)值越大表示被吸附物與催化劑Cu@CNT表面之間的相互作用越強(qiáng)，Eads為正值時(shí)則相反[25]。采用NEB法搜索甲醇在Cu@CNT表面裂解過(guò)程中所有基元反應(yīng)的過(guò)渡態(tài)(TS)，并通過(guò)虛頻計(jì)算進(jìn)行識(shí)別[26-27]?；罨軌?Ea)的計(jì)算公式為：

Ea=ETS-EIS.

(2)

式中：ETS、EIS分別為各基元反應(yīng)中過(guò)渡態(tài)(TS)和初始態(tài)(IS)的能量。速率k可通過(guò)傳統(tǒng)過(guò)渡態(tài)理論計(jì)算，k值與指前因子(A)和活化能(Ea)有關(guān)[28]，如式(3).R和T為氣體常數(shù)和溫度(本文中為298 K)，kB和h是普朗克常數(shù)和玻爾茲曼常數(shù)，QIS和QTS分別為初始狀態(tài)和過(guò)渡狀態(tài)的分配函數(shù)。

(3)

2 結(jié)果與討論

2.1 甲醇及其他中間體的吸附

通過(guò)研究甲醇及其他中間體在幾個(gè)不同吸附位點(diǎn)的吸附情況，確定了各物種的最優(yōu)吸附位點(diǎn)。吸附能和吸附構(gòu)型的相關(guān)結(jié)構(gòu)參數(shù)列于表1中，構(gòu)型圖如圖2所示。

表1 Cu@CNT表面甲醇裂解過(guò)程中各物種的吸附能及結(jié)構(gòu)參數(shù)Table 1 Adsorption energies and structural parameters of various species during methanol decomposition on Cu@CNT

圖2 Cu@CNT表面甲醇裂解過(guò)程中各物種的穩(wěn)定吸附構(gòu)型Fig.2 Most stable adsorption configurations of various species during methanol decomposition on Cu@CNT

2.1.1CHxOH(x=2，3)

計(jì)算結(jié)果表明，甲醇的最穩(wěn)定吸附構(gòu)型是通過(guò)O原子穩(wěn)定吸附在Cu@CNT面的T1位點(diǎn)，其吸附能為-0.28 eV，甲醇在Cu(110)[29]、Cu(100)[30]和Cu(111)[31]表面的吸附能分別為-0.49 eV、-0.44 eV、-0.36 eV。此外，氣相中甲醇的O—H鍵、C—O鍵和C—H鍵的鍵長(zhǎng)分別為0.096 1 nm、0.142 5 nm和0.111 0 nm，吸附后甲醇的O—H鍵和C—O鍵伸長(zhǎng)到0.096 6 nm和0.145 6 nm，而C—H鍵縮短到1.105 nm，可見(jiàn)吸附后甲醇變化較小，表明甲醇在Cu@CNT表面的吸附較弱。羥甲基(CH2OH)是甲醇裂解過(guò)程中的重要中間體。CH2OH更傾向于通過(guò)C原子吸附在Cu@CNT面的T2位點(diǎn)，吸附能為-2.74 eV，而CH2OH在Cu(110)[29]和Cu(100)[30]表面的吸附能分別為-1.97 eV和-1.83 eV.C—O鍵和C—C鍵的鍵長(zhǎng)分別為1.454 nm和1.572 nm.

2.1.2CHxO(x=0,1,2,3)

如圖2所示，CH3OH經(jīng)O—H鍵斷裂形成的甲氧基(CH3O)優(yōu)選通過(guò)O原子吸附在Cu@CNT面T2位點(diǎn)，吸附能為-2.44 eV.CH3O中的O原子與Cu@CNT表面的C原子形成的C—O鍵與管軸方向垂直，鍵長(zhǎng)為0.149 9 nm.此外，吸附物CH3O中C—O鍵鍵長(zhǎng)(0.145 8 nm)比自由基CH3O中的C—O鍵(0.139 6 nm)長(zhǎng)0.062 nm.CH2O由CH3O中C—H鍵斷裂形成，其最穩(wěn)定的吸附構(gòu)型是通過(guò)O原子吸附在Cu@CNT表面的T1位點(diǎn)，吸附能為-0.38 eV.此外，CH2O中C—O鍵平行于催化劑管軸方向。CHO的最優(yōu)吸附位為T2位，吸附能為-1.52 eV，CHO在Cu(110)[29]面和Cu(100)[30]面的的吸附能分別為-1.52 eV和-2.23 eV.C—O鍵和C—H鍵的長(zhǎng)度分別為0.121 1 nm和0.112 9 nm.此外C—O—H所形成的鍵角為122°.CO優(yōu)選以通過(guò)C原子吸附在Cu@CNT面的T2位點(diǎn)，吸附能為-0.47 eV，而CO在Cu(110)面[29]、Cu(100)面[30]和Cu(111)面[31]的吸附能分別為-1.21 eV、-1.04 eV和-0.98 eV.

如表1所示，CH3OH、CO及CH2O的吸附能遠(yuǎn)低于其他物種的吸附能，產(chǎn)生此結(jié)果的原因是自由基中間體中具有未成對(duì)電子，與Cu@CNT的相互作用強(qiáng)于分子中間體。此外，CO的吸附能為-0.47 eV，有助于其從催化劑表面脫附，提高催化劑對(duì)CO中毒的耐受性。Cu@CNT催化劑中C—Cu鍵基本穩(wěn)定在0.208 7 nm左右，表明該催化劑構(gòu)型是穩(wěn)定的。

2.2 甲醇裂解反應(yīng)

本部分討論了甲醇在Cu@CNT催化劑作用下的裂解機(jī)理。甲醇的裂解是從其以氣態(tài)形式吸附到Cu@CNT催化劑表面開始的。甲醇有三種可能的裂解路徑，如圖3所示：路徑1，甲醇通過(guò)O—H鍵斷裂生成CH3O中間體，最終生成CO或C；路徑2，甲醇通過(guò)C—H鍵斷裂生成CH2OH中間體，最終生成CO或C；路徑3則為甲醇通過(guò)C—O鍵斷裂生成CH3中間體，隨后通過(guò)C—H鍵斷裂產(chǎn)生C.各基元反應(yīng)的反應(yīng)熱、活化能壘及速率如表2.主反應(yīng)路徑中各基元反應(yīng)的構(gòu)型圖如圖4所示。

圖3 Cu@CNT表面甲醇裂解的可能路徑Fig.3 All possible paths of methanol decomposition on Cu@CNT

表2 Cu@CNT表面甲醇裂解過(guò)程中各基元反應(yīng)的反應(yīng)熱ΔH、活化能壘Ea和在298 K下的速率kTable 2 Reaction energy ΔH, activation energy barrier Ea, and rate constant k at 298 K for the elementary steps involved in methanol decomposition on Cu@CNT

圖4 Cu@CNT表面甲醇裂解主路徑反應(yīng)的IS，TS和FS的優(yōu)化構(gòu)型Fig.4 Optimized configuration of IS, TS, and FS for methanol decomposition reaction on Cu@CNT

2.2.1CH3OH的裂解

甲醇初始裂解有以下三種不同的路徑：分別為O—H、C—H和C—O鍵斷裂。路徑1始于吸附在T1位的CH3OH，終于通過(guò)O原子吸附在T2位的CH3O和吸附于T1鄰位的H.在TS1中，O—H鍵長(zhǎng)從IS中的0.096 6 nm增加到0.119 7 nm.該過(guò)程的活化能為1.01 eV，反應(yīng)熱為0.54 eV.而在Cu(100)表面CH3OH中C—H鍵斷裂的活化能為1.60 eV[30].在路徑2中甲醇通過(guò)C—H鍵斷裂生成CH2OH和H.通過(guò)O原子吸附于T1位的CH3OH作為IS，吸附于T2位點(diǎn)CH2OH和吸附于T1位的H原子作為FS.C—H鍵長(zhǎng)從IS的0.110 8 nm伸長(zhǎng)到TS中的0.164 2 nm.計(jì)算結(jié)果表明該基元反應(yīng)的活化能為1.40 eV，反應(yīng)熱為0.28 eV.CH3OH還可通過(guò)C—O鍵斷裂生成CH3和OH，CH3OH中C—O鍵被活化生成通過(guò)C原子吸附于T2位點(diǎn)的CH3和垂直吸附于T2對(duì)位的OH.該基元反應(yīng)的活化能為1.99 eV.

甲醇初步裂解的勢(shì)能圖如圖5所示。比較甲醇初始裂解的三種不同的路徑發(fā)現(xiàn)，CH3OH通過(guò)O—H鍵斷裂形成CH3O和H，比通過(guò)C—H鍵斷裂形成CH2OH和H和通過(guò)C—O鍵斷裂形成CH3和OH具有更低的活化能，因此甲醇初步裂解優(yōu)選通過(guò)O—H鍵斷裂生成CH3O.此外，O—H鍵斷裂(1.30×10-4s-1)的速率遠(yuǎn)大于C—O鍵(1.42×10-19s-1)和C—H鍵(2.12×10-9s-1)的斷裂速率。因此后續(xù)只討論路徑1.

圖5 Cu@CNT表面甲醇初始裂解的勢(shì)能圖Fig.5 Potential energy profile of methanol initial decomposition on Cu@CNT

2.2.2CH3O的裂解

CH3O的后續(xù)裂解存在兩種路徑。一種是CH3O脫氫形成CH2O，另一種則是通過(guò)C—O鍵斷裂形成CH3和O，勢(shì)能圖如圖6所示。

對(duì)于CH3O脫氫，吸附于T2的CH3O作為IS.在FS中，CH2O從表面上升，并且H原子移動(dòng)至T1鄰位。C—O鍵逐漸向Cu@CNT表面傾斜(角度α從IS的28°增大到TS的69°)，使得C—H鍵更加靠近催化劑表面以促進(jìn)其活化。該過(guò)程的反應(yīng)熱和活化能分別為0.45 eV和0.73 eV.而CH3O在Cu(100)表面和Cu (111)表面C—H鍵斷裂的活化能分別為3.08 eV[30]和1.85 eV[31].對(duì)于C—O鍵的斷裂，活化能為2.08 eV，遠(yuǎn)高于C—H斷裂的活化能，表明對(duì)于CH3O難以發(fā)生C—O鍵的斷裂。此外，在298 K時(shí)C—H鍵和C—O鍵斷裂的速率分別為4.72×104s-1和2.42×10-22s-1.

2.2.3CH2O的裂解

如圖6所示，在Cu@CNT表面，CH2O的不同分解路徑會(huì)生成兩種不同的中間體，分別為CH2和CHO.CH2O通過(guò)脫氫反應(yīng)可生成CHO和H原子，在過(guò)渡態(tài)TS6中C—H鍵長(zhǎng)從IS中的0.112 7 nm延長(zhǎng)到0.128 7 nm，F(xiàn)S中CHO和H原子分別吸附在T2位和T1鄰位。該反應(yīng)過(guò)程的活化能為0.72 eV，同時(shí)反應(yīng)熱為-0.13 eV.而CH2O在Cu(100)表面和Cu(111)表面C—H鍵斷裂的活化能分別為0.85 eV[30]和1.15 eV[31].使CH2O通過(guò)C—O鍵斷裂生成CH2和O需要克服2.15 eV的活化能，同時(shí)反應(yīng)熱為1.51 eV.由此可知，與C—H鍵斷裂相比，C—O鍵斷裂不是CH2O裂解的可能路徑。此外反應(yīng)速率結(jié)果顯示C—H鍵斷裂的速率(3.47×102s-1)是C—O鍵斷裂速率(4.54×10-24s-1)的1026倍。

2.2.4CHO的裂解

CHO在T2位點(diǎn)生成后，可進(jìn)一步通過(guò)C—H斷裂生成CO和H，或通過(guò)C—O鍵斷裂生成CH和O，勢(shì)能圖如圖6所示。對(duì)于C—H鍵的斷裂，該過(guò)程從吸附在T2位點(diǎn)的CHO中C—H鍵的伸縮振動(dòng)開始。C—H鍵長(zhǎng)從IS中的0.112 9 nm伸長(zhǎng)到TS中的0.130 0 nm，O—C—H鍵角從IS中的122.55°增大到TS中232.55°.在FS中，H原子吸附于T2位點(diǎn)，而CO移動(dòng)到T3位點(diǎn)。該過(guò)程的活化能為0.69 eV，反應(yīng)熱為0.14 eV.而CHO在Cu(100)表面C—H鍵斷裂的活化能為0.84 eV[30]。此外，C—O鍵斷裂的活化能為2.64 eV.與CHO脫氫相比，由于C—O鍵斷裂的活化能遠(yuǎn)高于C—H鍵斷裂的活化能，CHO更易發(fā)生脫氫反應(yīng)。CHO中C—H鍵和C—O鍵斷裂的速率分別為5.95×103s-1和1.30×10-32s-1.

圖6 Cu@CNT表面CH3O后續(xù)裂解的勢(shì)能圖Fig.6 Potential energy profile of CH3O subsequent decomposition on Cu@CNT

綜上所述，甲醇在Cu@CNT上裂解的最優(yōu)路徑是CH3OH→CH3O→CH2O→CHO→CO，限速步驟是CH3OH→CH3O.

3 結(jié)論

本文運(yùn)用密度泛函理論探索了甲醇在Cu@CNT表面的裂解反應(yīng)機(jī)理，得到以下結(jié)論：甲醇穩(wěn)定吸附于T1位點(diǎn)，其他物種則穩(wěn)定吸附于T2位點(diǎn)。對(duì)于甲醇的初始裂解，其更易通過(guò)O—H鍵斷裂生成CH3O，其后CH3O逐步脫氫最終生成CO，其中O—H鍵的斷裂是限速步驟。且在Cu@CNT表面，甲醇裂解主反應(yīng)路徑中各基元反應(yīng)的活化能分別為1.01 eV、0.73 eV、0.72 eV和0.69 eV.甲醇在Cu(100)表面裂解的主路徑與甲醇在Cu@CNT表面裂解的主路徑相同，各基元反應(yīng)的活化能分別為1.60 eV、3.08 eV、0.85 eV和0.84 eV，故可知Cu@CNT催化劑的催化活性高于Cu(100)催化劑的催化活性。此外，甲醇裂解的最終產(chǎn)物CO在Cu@CNT和Cu(100)表面的吸附能分別為-0.47 eV和-1.04 eV，故Cu@CNT具有更高的CO耐受性。不同于活性位點(diǎn)裸露于反應(yīng)介質(zhì)中的Cu(100)催化劑，Cu@CNT催化劑中碳納米管的封裝作用使活性中心和反應(yīng)介質(zhì)隔離，反應(yīng)過(guò)程中無(wú)積碳的生成。由此可知，碳納米管封裝金屬粒子有效地增強(qiáng)了金屬催化劑的活性和穩(wěn)定性，是甲醇裂解的高效穩(wěn)定催化劑。