Ni，Pd，Pt負載α-MoC催化水煤氣變換反應(yīng)理論

史 肖，鄒雪燕，黃 偉，左志軍

(1.太原理工大學(xué) 省部共建煤基能源清潔高效利用國家重點實驗室,山西 太原 030024; 2.太原理工大學(xué) 煤科學(xué)與技術(shù)教育部重點實驗室,山西 太原 030024)

水煤氣變換反應(yīng)(WGS，CO(g)+H2O(g)→CO2(g)+H2(g))是CO去除和H2制備的途徑[1]，也是合成氣轉(zhuǎn)化過程中的重要反應(yīng)，如甲醇合成[2]和費托合成[3]。常見的WGS反應(yīng)催化劑有Au，Cu，Pd，Ni和Pt基等催化劑[4-8]。為了提高催化劑的活性，通常通過添加助劑[9-12]或改變載體[13-15]的方法。

過渡金屬碳化物是一種金屬間填充型化合物，是由碳原子填隙式融進過渡性金屬的晶格中形成。由于其表面性質(zhì)和催化活性類似于Pt等貴金屬，目前廣泛應(yīng)用于催化加氫、脫氫、WGS和異構(gòu)化等反應(yīng)[16-17]。如，Au/α-MoC用于WGS反應(yīng)時，反應(yīng)速率和CO轉(zhuǎn)化率分別是Cu/Zn/Al2O3催化劑的354倍和17.8倍[18]，同時Au/α-MoC催化劑的活性明顯優(yōu)于Au/β-MoC。α-MoC對WGS反應(yīng)起到了重要的作用，因此，筆者選取了α-MoC作為載體。

理論計算結(jié)果顯示，α-MoC(111)載體對WGS反應(yīng)沒有催化效果，但是能夠促進H2O的解離。當其負載Au和Cu后能夠提高WGS反應(yīng)在2種催化劑表面的活性[18-19]。Pd，Ni，Pt基催化劑也是WGS反應(yīng)的常用催化劑，為了進一步了解α-MoC載體在WGS反應(yīng)中的作用和WGS反應(yīng)機理，筆者使用密度泛函理論(DFT)和動力學(xué)蒙特卡羅(KMC)方法系統(tǒng)的研究這3種金屬負載α-MoC催化劑上的WGS反應(yīng)過程，研究不同反應(yīng)溫度時的反應(yīng)速率和催化轉(zhuǎn)換頻率(TOF)，并對WGS反應(yīng)催化劑的優(yōu)化提供指導(dǎo)。

1 計算方法和模型

計算使用VASP軟件[20-22]，采用了GGA-PBE泛函求解電子交換相關(guān)能[23]。平面波截斷能取值為 415 eV，布里淵區(qū)的k點選擇 3×3×1[24]。在結(jié)構(gòu)優(yōu)化的過程中當力的變化小于0.1 eV/nm和能量的變化小于1.0×105eV/atom時達到收斂標準。采用CI-NEB方法進行過渡態(tài)搜索，對達到收斂標準的過渡態(tài)進行頻率分析，通過唯一虛頻確定過渡態(tài)[25]。

首先對金屬晶胞M(M=Ni，Pd，Pt)和α-MoC晶胞進行了優(yōu)化，優(yōu)化后對應(yīng)的的晶格常數(shù)分別為aNi=0.351 3 nm，aPd=0.393 7 nm，aPt=0.396 7 nm，aα-MoC=0.437 0 nm，分別與實驗值aNi=0.352 4 nm，aPd=0.389 1 nm，aPt=0.392 4 nm和aMoC=0.427 8 nm接近[26-27]。因此，選擇的計算方法和參數(shù)對本文的計算體系是合理的。α-MoC(111)采用7層3×3的平板模型，1.5 nm的真空層。由于H2O的快速解離，α-MoC表面的Mo位點被O占據(jù)[19]，如圖1所示，圖1中top,hol,bri為吸附位點。計算過程中固定底面兩層，其余原子弛豫。金屬M4簇負載存在2種構(gòu)型:平面簇構(gòu)型(2d rhombicM4)和四面體簇(3d-tetrahedralM4)，對應(yīng)的吸附能分別為Ni4(2d:-8.41 eV;3d:-7.63 eV)，Pd4(2d:-8.56 eV;3d:-6.77 eV)，Pt4(2d:-10.20 eV;3d:-9.32 eV)。計算結(jié)果表明，平面簇構(gòu)型更穩(wěn)定，故僅考慮了2d-M4負載在α-MoC(111)表面。

圖1 M4/α-MoC(111)的側(cè)視圖和俯視圖Fig.1 Top and side view of M4/α-MoC(111)

2 結(jié)果與討論

對于水煤氣變換反應(yīng)，反應(yīng)機理主要有2種:① 氧化還原機理:CO直接與H2O分解產(chǎn)生的氧結(jié)合，生成CO2;② 甲酸鹽途徑或者羧酸鹽途徑:CO與H2O分解產(chǎn)生的羥基結(jié)合，生成甲酸鹽或羧酸鹽中間體，最后分解成CO2[28-30]。反應(yīng)過程中所涉及到的中間體，反應(yīng)物和產(chǎn)物的最穩(wěn)定吸附構(gòu)型如圖2所示。

圖2和表1顯示，反應(yīng)物、中間體和產(chǎn)物在3類催化劑表面的最優(yōu)吸附位相近，只有*H在Pt4/α-MoC(111)和*CO2在Ni4/α-MoC(111)的吸附例外。對中間體來說(*O，*OH，*CO2，*H2O，*CHOO，*COOH)，它們在Pd4/α-MoC(111)和Pt4/α-MoC(111)的吸附能相近，吸附穩(wěn)定性低于對應(yīng)物種吸附在Ni4/α-MoC(111)。

圖2 WGS反應(yīng)所涉及的各中間體在Ni4/α-MoC(111)的穩(wěn)定吸附構(gòu)型Fig.2 Adsorption configurations of intermediates involved in the WGS on the surfaces of Ni4/α-MoC(111)

表1 M4/α-MoC(111)面WGS反應(yīng)過程中涉及到的中間體，反應(yīng)物和產(chǎn)物的最優(yōu)吸附位和吸附能Table 1 Adsorption energy and adsorption site of the possible intermediates,reactants and productions during the WGS reaction on M4/α-MoC(111)

反應(yīng)物*CO吸附在Ni4/α-MoC(111)和Pd4/α-MoC(111)的吸附能相近(-2.23和-2.21 eV)，而在Pt4/α-MoC(111)面的吸附穩(wěn)定性低于Ni4/α-MoC(111)和Pd4/α-MoC(111)面。相比CO在Ni(111)，Pd(111)和Pt(111)的吸附能，載體α-MoC提高了CO在Ni和Pd的吸附穩(wěn)定性，但是降低了在Pt上的吸附穩(wěn)定性[31-32]。對*H2O來說，載體α-MoC一定程度上提高了其在Ni，Pd和Pt表面的吸附穩(wěn)定性[31-32]。

優(yōu)化后，產(chǎn)物*CO2遠離Pd4/α-MoC(111)和Pt4/α-MoC(111)面，但是其吸附在Ni4/α-MoC(111)的bri位，這是由于*CO2在Ni4/α-MoC(111)面的吸附穩(wěn)定性強于Pd4/α-MoC(111)和Pt4/α-MoC(111)面。總體來說，相較Ni(111)，Pd(111)和Pt(111)，α-MoC(111)提高了CO2在3種催化劑表面的吸附穩(wěn)定性[31-32]。對于H2來說，其在3個表面的吸附能力很弱，載體對其吸附基本沒有影響。

*H2O的解離(*H2O+*→*OH+*H)是WGS反應(yīng)的第1步，在Ni4/α-MoC(111)上*H2O解離的能壘為1.16 eV。對于氧化還原機理，*O生成存在2種方式:①*OH直接裂解(*OH+*→*O+*H)能壘Ea為0.95 eV;②*OH歧化反應(yīng)(*OH+*OH→*H2O+*O)能壘為0.37 eV。最后，*CO氧化生成*CO2需克服1.62 eV的能壘。對于甲酸鹽和羧酸鹽路徑，*CO與*OH反應(yīng)，生成*CHOO(Ea=2.84 eV)或*COOH(Ea=1.61 eV)。結(jié)果表明，在Ni4/α-MoC(111)上氧化還原機理容易發(fā)生，速控步驟為*CO氧化。

在Pd4/α-MoC(111)上，*H2O克服1.23 eV能壘解離生成*OH 和*H。當發(fā)生氧化還原反應(yīng)時，*OH歧化反應(yīng)的能壘為0.35 eV，而*OH直接解離的能壘為1.02 eV。最后，*CO氧化生成*CO2的能壘為0.64 eV。在發(fā)生甲酸鹽或羧酸鹽路徑中，*CO與*OH反應(yīng)生成*CHOO和*COOH的能壘分別為1.64和1.28 eV。因此，在Pd4/α-MoC(111)上，WGS反應(yīng)的氧化路徑為主要路徑，H2O的解離為速控步驟。

在Pt4/α-MoC(111)上，H2O的解離能壘為0.86 eV。當發(fā)生氧化還原反應(yīng)時，*OH+*→*O+*H和*OH+*OH→*H2O+*O的反應(yīng)能壘分別為1.28和0.35 eV，*CO氧化生成*CO2的能壘為0.81 eV。在發(fā)生甲酸鹽或羧酸鹽路徑中，*CO與*OH反應(yīng)生成*CHOO和*COOH的能壘分別為1.72和0.66 eV。由于生成*COOH的能壘遠遠低于生成*HCOO的能壘，因此*CO和*OH反應(yīng)優(yōu)先生成*COOH。最后，經(jīng)*COOH與*OH反應(yīng)生成*CO2和*H2O(Ea=0.42 eV)。

基于此，根據(jù)WGS反應(yīng)的正逆反應(yīng)能壘和不同溫度時的指前因子，利用KMC研究了WGS反應(yīng)在Ni4/α-MoC(111)，Pd4/α-MoC(111)和Pt4/α-MoC(111)上的反應(yīng)路徑和催化轉(zhuǎn)換頻率。對于反應(yīng)物和產(chǎn)物的吸附和脫附，考慮了熵的貢獻[33]。KMC模型和計算參數(shù)等詳細信息見文獻[34]。

KMC結(jié)果顯示，標準大氣壓(n(CO)∶n(H2O)=1∶1，物質(zhì)的量之比)和反應(yīng)溫度為500 K時，在Ni4/α-MoC(111)和Pd4/α-MoC(111)催化劑上，水煤氣變換的反應(yīng)機理為氧化還原機理如圖3所示，圖3中,TS為過渡態(tài)。即CO(g)+H2O(g)→*CO+*H2O→*CO+*OH+*H→*CO+*O+*2H→*CO2+*H2→CO2(g)+H2(g)。需要指出的是，盡管*OH+*OH→*H2O+*O的反應(yīng)能壘遠遠小于*OH+*→*O+*H的反應(yīng)能壘，然而由于高的*H2O→*OH+*H反應(yīng)能壘導(dǎo)致表面*OH覆蓋度過低。因此，KMC結(jié)果顯示*OH的直接解離是*O產(chǎn)生的反應(yīng)途徑。在Pt4/α-MoC(111)催化劑上，水煤氣變換的反應(yīng)路徑為羧酸鹽路徑如圖4所示，即CO(g)+2H2O(g)→*CO+2*H2O→*CO+2*OH+2*H→*COOH+*OH+2*H→*CO2+*H2O+*H2→CO2(g)+H2O(g)+H2(g)。由于H2和CO2的TOF相等，因此圖5僅列出了H2的TOF(反應(yīng)條件為一個大氣壓下，CO與H2O物質(zhì)的量之比為1∶1)?？傮w來說，H2的TOF隨反應(yīng)溫度的升高而升高。H2在Ni4/α-MoC(111)和Pd4/α-MoC(111)的TOF偏低。在500 K時催化劑單位活性位點上轉(zhuǎn)換次數(shù)約為0和0.056 s-1，這是由于WGS反應(yīng)過程中反應(yīng)能壘較高造成的。對于Pt4/α-MoC(111)，值得注意的是，400 K時，H2在催化劑單位活性位點上轉(zhuǎn)換次數(shù)約為0，這是因為Pt4/α-MoC(111)上CO的脫附能較高(1.07 eV)，活性位點被CO所覆蓋。隨著反應(yīng)溫度的升高，CO的脫附能(500 K下脫附能為0.83 eV)降低，H2在催化劑單位活性位點上的TOF約為6.3 s-1，WGS反應(yīng)發(fā)生。結(jié)果表明，*CO的強穩(wěn)定性不利于WGS反應(yīng)的進行。當反應(yīng)溫度為500 K時，H2的TOF的排序是:Pt4/α-MoC(111)?Pd4/α-MoC(111)>Ni4/α-MoC(111)。因此，當3種金屬負載在α-MoC載體上時，Pt4/α-MoC(111)催化劑活性最好。同時，其反應(yīng)活性高于Pt/Al2O3(TOF≈0.7 s-1)和Pt/TiO2(TOF≈0.2 s-1)[35]，這表明Pt/α-MoC是WGS反應(yīng)的高活性催化劑。

圖3 WGS反應(yīng)在Ni4/α-MoC(111)，Pd4/α-MoC(111)上反應(yīng)勢能Fig.3 Potential energy of WGS on Ni4/α-MoC(111) and Pd4/α-MoC(111)

圖4 WGS反應(yīng)在Pt4/α-MoC(111)上反應(yīng)勢能Fig.4 Potential energy of WGS on Pt4/α-MoC(111)

圖5 WGS反應(yīng)在M4/α-MoC(111)上的H2的TOFFig.5 H2 TOF for WGS reaction on M4/α-MoC(111)

3 結(jié) 論

(1)通過密度泛函理論和動力學(xué)蒙特卡洛方法探究了過渡金屬Ni,Pd,Pt負載α-MoC催化劑的WGS反應(yīng)的反應(yīng)機理和活性。

(2)DFT計算結(jié)果發(fā)現(xiàn)，在Ni4/α-MoC(111)和Pd4/α-MoC(111)上，WGS反應(yīng)的反應(yīng)路徑為氧化還原路徑為CO(g)+H2O(g)→*CO+*H2O→*CO+*OH+*H→*CO+*O+*2H→*CO2+*H2→CO2(g)+H2(g);在Pt4/α-MoC(111)催化劑上，WGS反應(yīng)通過羧酸鹽路徑發(fā)生，即CO(g)+2H2O(g)→*CO+2*H2O→*CO+2*OH+2*H→*COOH+*OH+2*H→*CO2+*H2O+*H2→CO2(g)+H2O(g)+H2(g)。

(3)動力學(xué)模擬結(jié)果表明，Pt4/α-MoC(111)具有較高的催化活性。標準大氣壓下，反應(yīng)溫度在400～500 K，H2O與CO的物質(zhì)的量之比為1時，H2的TOF排序:Pt4/α-MoC(111)?Pd4/α-MoC(111)>Ni4/α-MoC(111)，Pt/α-MoC是WGS反應(yīng)的高活性催化劑。