γ-Fe2O3抗As2O3中毒能力的分子模擬

周文波,牛勝利*,劉思彤,王 棟,韓奎華,王永征

γ-Fe2O3抗As2O3中毒能力的分子模擬

周文波1,牛勝利1*,劉思彤1,王 棟2,韓奎華1,王永征1

(1.山東大學能源與動力工程學院,高效節(jié)能及儲能技術與裝備山東省工程實驗室,山東 濟南 250061；2.佐治亞理工學院布魯克貝爾可持續(xù)發(fā)展學院,佐治亞 亞特蘭大 30332)

采用密度泛函理論研究了γ-Fe2O3表面As2O3的吸附以及摻雜改性提高抗As2O3中毒性能的作用機理.計算了As2O3在完整以及O缺陷γ-Fe2O3(001)表面的吸附性能,包括吸附位點、吸附結構、吸附能、PDOS等.同時建立了Mo、Ti、Mg摻雜的γ-Fe2O3模型,探討了助劑摻雜對抗砷中毒能力的作用機制,并考慮了摻雜量的影響.結果表明:As2O3傾向于以O端化學吸附在γ-Fe2O3(001)表面Feoct位,吸附過程發(fā)生強烈的相互作用和電荷轉移.當表面存在O缺陷時,As2O3的吸附能得到提高.Mo、Ti、Mg傾向于摻雜在Feoct位,增強了對As2O3的吸附能力,并且增大Mo的摻雜量可以強化As2O3的吸附.As2O3傾向于與活性較強的Mo、Ti、Mg發(fā)生反應,從而保護活性Fe位不受砷中毒,Ti和Mg的摻雜還抑制了相鄰Fe位對As2O3的吸附.Mo、Ti、Mg的摻雜還促進了催化劑表面對NH3的吸附,增強了表面酸性強度,有利于SCR反應.Mo、Ti、Mg原子的摻雜有利于提高γ-Fe2O3催化劑的抗砷中毒性能.

As2O3；吸附；γ-Fe2O3；O缺陷；Mo、Ti、Mg摻雜；密度泛函理論

煤燃燒產生的氮氧化物會引起一系列環(huán)境問題,如酸雨和臭氧破壞等[1-3].面對巨大的環(huán)境壓力和嚴格的排放標準,NH3選擇性催化還原(NH3-SCR)仍然是主流的NO控制技術[4-5].商業(yè)化的V2O5- WO3/TiO2催化劑存在反應溫度窗口窄、抗毒性差以及釩的生物毒性等諸多缺點[6-8].因此,設計更合理、更環(huán)保的中低溫SCR催化劑勢在必行.

由于來源廣泛,環(huán)境友好,催化性能優(yōu)異等優(yōu)點,鐵氧化物催化劑已迅速引起研究人員的關注.Liu等[9]比較了α-Fe2O3和γ-Fe2O3的催化活性和選擇性,發(fā)現NO和NH3分子都更容易吸附在γ-Fe2O3表面,在150～300℃時γ-Fe2O3表現出了比α-Fe2O3更高的催化活性.張信莉等[10]制備了Mn摻雜的γ-Fe2O3催化劑,發(fā)現Mn的摻雜拓寬了脫硝活性窗口,并實現了在125～200℃效率保持100%.Xu等[11]在γ-Fe2O3催化劑中引入Mg,提高了250～350℃溫度范圍內的催化性能,NO轉化率超過90%.上述研究充分表明,γ-Fe2O3可以成為良好的中低溫SCR催化劑,但實際煙氣中的一些有害成分會導致催化劑失活.

煙氣中,重金屬砷主要以As2O3氣體形式存在,通常被認為是對SCR催化劑產生不利影響的最重要的形態(tài)[12-13].As2O3與催化劑結合阻止氨的吸附或阻塞催化劑內部通道,從而抑制NH3和NO在SCR催化劑上的活化[14].同時,As2O3能夠破壞催化劑表面Lewis酸位點和Br?nsted酸位點,減弱酸性位點對NH3的吸附.Li等[15]發(fā)現As2O3含量達到2.9%時,催化劑明顯失活,在250～500℃時NO轉化率下降到40%以下.Peng等[16]研究了As2O3對釩基催化劑失活的影響,認為主要是由于Lewis酸位點和不穩(wěn)定的Br?nsted酸位點數量減少所致.Li等[17]研究了砷對γ-Fe2O3催化劑的毒害行為,發(fā)現在325℃時催化劑上的NO轉換率損失了44%.因此,迫切需要提高SCR催化劑的抗砷中毒能力,以提高其使用壽命和性能.

另一方面,助劑摻雜可以有效改善SCR催化劑的抗砷中毒性能.Peng等[18]比較了MoO3和WO3摻雜的V2O5/TiO2催化劑的抗砷中毒機理.研究發(fā)現,氧化砷更容易與活性較高的MoO3反應,從而保護催化劑上的活性物種.Li等[19]在γ-Fe2O3催化劑中引入Mo后發(fā)現Mo在200～400℃溫度范圍內提高了脫硝活性,并阻止了活性鐵和氣相砷的結合,對催化劑的抗砷中毒起很大的作用.通過Mg和Ti摻雜改性的γ-Fe2O3催化劑在100～350℃的催化活性和抗砷性能都明顯提高[20].

由于密度泛函理論(DFT)具有在分子尺度上探測反應機理的強大功能,越來越多的研究者將其應用于催化劑砷吸附與中毒研究.Wu等[21]利用DFT探究了WO3和MoO3摻雜對釩基催化劑與As2O3相互作用的影響,發(fā)現As2O3易與MoO3發(fā)生反應,從而保護活性成分V2O5.Xing等[22]用DFT模擬發(fā)現NH3和As2O3分子在CuO(111)表面的吸附主要以物理吸附為主,As2O3比NH3更容易吸附在CuO表面.Hu等[23]利用分子模擬發(fā)現As2O3的吸附取決于所處的位置,物理吸附是由于阻礙了物理吸附向化學吸附轉變的高能壘所致.然而,采用DFT研究γ-Fe2O3表面As2O3的吸附及抗中毒性能的研究很少,特別是在吸附過程中相互作用的具體微觀解釋還鮮有報道.

本研究采用密度泛函理論(DFT)方法研究了As2O3在γ-Fe2O3(001)表面以及缺陷表面的吸附性能和Mo、Ti、Mg摻雜對γ-Fe2O3抗砷中毒的影響.研究結果可為γ-Fe2O3催化劑的開發(fā)和耐砷中毒機理的探索提供理論參考.

1 計算模型和方法

1.1 計算模型

本研究建立了γ-Fe2O3的晶胞模型(圖1).γ- Fe2O3為反尖晶石型結構,晶胞參數為===8.405?,===90?[24].每個晶胞含有32個O陰離子和21.33個Fe陽離子,其中八面體Fe表示為Feoct,四面體Fe表示為Fetet[25].由于γ-Fe2O3(001)面具有最強的反應活性,因此被選作反應面[26-27].如圖2所示,模型選取7層原子,底部4層原子固定,頂部3層原子松弛[28].設置了P(2×2)的平面模型,以及為了防止兩個周期性平板之間的相互作用,設置了15?的真空層.

圖1 γ-Fe2O3的晶胞模型

圖2 γ-Fe2O3(001)表面結構

圖3 幾何優(yōu)化的As2O3

1.2 計算方法

所有DFT計算均基于CASTEP程序包[29].擴展平面波采用具有Perdew-Burke-Ernzerhof交換相關函數的廣義梯度近似(GGA-PBE)[30-31].截止能量為400eV,k點網格為2×2×1.離子核用超軟贗勢描述,并考慮自旋極化.結構優(yōu)化收斂精度:能量為2.0×10-5eV/atom,最大應力為0.05eV/?和最大位移為0.002?,自洽場SCF收斂標準為2.0×10-6eV/atom.在γ-Fe2O3中,陽離子占據在四面體和八面體位置,由于四面體晶格和八面體晶格中陽離子的自旋取向為反平行排列,導致亞鐵磁行為[32-33].計算允許自旋極化,以反映磁赤鐵礦γ-Fe2O3的磁性,磁矩設置遵循亞鐵磁性排列,四面體和八面體位點上Fe的磁矩方向相反,并采用高自旋態(tài).

優(yōu)化后γ-Fe2O3的晶格常數為8.357?,與實驗數據相符[24].在尺寸為10?′10?′10?晶格中進行As2O3的幾何優(yōu)化,優(yōu)化結果如圖3所示.根據計算結果,As-O鍵長為1.862?,As-O-As角為79.365°,分別與參考值1.836?和80.6°非常接近[34].

ads定義為表面吸附能,計算見式(1):

ads=tot-ad-catal(1)

式中:ads、tot、ad和catal分別代表吸附能、氣體吸附在γ-Fe2O3(001)表面的系統(tǒng)總能、氣體分子的能量和催化劑的能量,單位eV.

Δ定義為As2O3與γ-Fe2O3(001)表面的電荷轉移,計算見式(2)[35]:

Δ=ad-gas(2)

式中: Δ、ad和gas分別表示轉移的電荷量、吸附態(tài)和氣態(tài)小分子的電荷,單位e.

氧空位形成能計算公式如式(3)[36]:

v=OV+ 0.5O2-surf(3)

式中:v、OV、O2和surf分別代表氧空位形成能、氧缺陷表面的能量、氧分子的能量以及完美表面的總能量,單位eV.

dope定義為表面結構的摻雜能,比較不同摻雜表面的穩(wěn)定性,計算見式(4)[37]:

dope=(M-Fe2O3)-(Fe2O3)+(Fe)-(M)(4)

式中:dope、(M-Fe2O3)、(Fe2O3)、(Fe)和(M)分別代表摻雜形成能、摻雜單M(Mo、Ti、Mg)原子的γ-Fe2O3表面的能量、γ-Fe2O3表面的能量、Fe和M單質晶胞內的單個原子的能量,單位eV.

2 結果與討論

2.1 γ-Fe2O3(001)表面As2O3的吸附

2.1.1 吸附構型與吸附能 根據γ-Fe2O3(001)表面的空間位置和對稱性,本文研究了As2O3分子在γ-Fe2O3(001)表面Feoct位、Fetet位、O位、hollow位和橋位5個初始吸附位點上的吸附[38].Hu等[23]通過對As2O3小分子進行分子軌道分析,發(fā)現As2O3分子的As原子端和O原子端,均可作為吸附中心.通過幾何優(yōu)化,根據結構特征得到了6種As2O3分子在γ-Fe2O3(001)表面的穩(wěn)定構型(圖4).

從表1可知,As2O3分子擇優(yōu)吸附的結合能大小順序為:ads(a)>ads(c)>ads(f)>ads(d)>ads(b)>ads(e).圖4中的構型(a)最穩(wěn)定,As2O3分子以其O端吸附在八面體Feoct位,形成了鍵長為2.202?的O-Fe鍵.構型(a)的吸附能為-0.75eV,電荷轉移0.19e,表明其吸附過程為放熱的化學吸附.此外,在構型(c)和(f)中,初始吸附位點為橋位,但優(yōu)化后As2O3分子發(fā)生移動,最終分別以O端和N端吸附在Feoct位.構型(c)、(f)和(d)的吸附能分別為-0.69,-0.51和-0.50eV,屬于化學吸附,其中As2O3分子的O原子端比As原子表現出更強的化學活性.在構型(b)中,As2O3分子弱化學吸附在hollow位,吸附過程為放熱反應、釋放能量0.41eV.在構型(e)中,As2O3分子在Fetet位的吸附,并未形成新的化學鍵,吸附能為-0.23eV,屬于物理吸附.從吸附能分析結果表明,As2O3分子在γ-Fe2O3表面的吸附既有物理吸附又有化學吸附,與不同的吸附位點有關,Feoct原子是最佳吸附活性位.As2O3更傾向于以其O端與催化劑表面的Feoct原子進行鍵合.課題組前期的研究表明,γ-Fe2O3催化劑表面的八面體鐵原子(Feoct)具有最好的氣體吸附能力,其中NH3主要吸附到Feoct位形成吸附態(tài)NH3,其中NO的最佳吸附方式是以N端彎曲的形式吸附在Feoct頂位[39].催化劑吸附As2O3后,活性位會被部分覆蓋,阻止了反應物的吸附,抑制了SCR反應.As2O3與活性鐵中心結合,占據活性位點,造成失活,與Li等[17]的實驗結果一致.

圖4 As2O3在γ-Fe2O3表面的吸附構型

表1 As2O3分子吸附在γ-Fe2O3表面的優(yōu)化構型

注:O-Feoct代表As2O3的O原子與Feoct原子之間的距離?As-Fe代表As2O3的As原子與Fe(Feoct或Fetet)原子之間的距離.

2.1.2 Mulliken電荷轉移 通過計算Mulliken電荷布局分析,可以反映吸附過程中電荷分布與轉移情況,為闡明其結構變化和吸附能提供了支持[40].從表2可知,As2O3分子在γ-Fe2O3表面的吸附過程中,除物理吸附的構型(e)外,其他構型的Mulliken電荷轉移均為正數,表明吸附過程電子主要從As2O3分子轉移到γ-Fe2O3表面.正如構型(a)～(c)中,當As2O3分子以O端吸附在催化劑表面時,轉移的電荷主要由O原子提供,As原子的幾乎沒有電荷轉移.當As2O3分子以As端吸附在催化劑表面時,與催化劑鍵合的As1原子得到電子,其他原子失去電子.As2O3分子最穩(wěn)定的吸附構型(a),是以其O端吸附在Feoct原子頂位.在構型(a)的吸附過程中,電荷轉移為0.19e,O1和、O2和O3分別向基底表面提供了0.05和0.08和0.05e的電子.As2O3分子中的所有O原子都失去了電子,尤其是O2原子,這是O-Fe鍵形成的主要因素.

表2 As2O3吸附在催化劑表面的Mulliken電荷布局

2.1.3 偏態(tài)密度(PDOS)分析 如圖5所示,在大約-5.2、-2.8、-1和2.2eV時,O原子的p軌道和Feoct的d軌道顯示出雜化產生的共振峰.這表明O原子和Feoct陽離子雜化形成穩(wěn)定的化學鍵以及穩(wěn)定的吸附構型.

圖5 構型a中O和Fe原子的PDOS圖

O?Fe分別為As2O3分子吸附在γ-Fe2O3表面的O原子以及表面的Feoct原子,實線?劃線和點劃線分別為s軌道?p軌道和d軌道

綜合以上分析可知,As2O3分子在γ-Fe2O3(001)表面的吸附方式,既有物理吸附也有化學吸附. As2O3分子更傾向于以其O端吸附在八面體Feoct位.As2O3分子吸附過程,會與γ-Fe2O3(001)表面發(fā)生強烈的相互作用和電荷轉移.氣態(tài)As2O3會與γ-Fe2O3上的活性鐵中心結合,占據活性位點,導致催化劑活性下降,引起催化劑失活.

2.2 O缺陷表面As2O3的吸附

2.2.1 O缺陷的γ-Fe2O3(001)表面模型 固體材料表面通常含有空位缺陷,其中氧空位缺陷是作為催化或者化學吸附的重要活性位點,進而影響催化活性[41-42].此外,氧空位缺陷是最普遍和廣泛研究的陰離子缺陷,在金屬氧化物表面氧空位形成能相對較低[25].首先,在圖2的γ-Fe2O3(001)表面結構的基礎上,去除原結構表面Feoct與Fetet原子之間的一個氧原子,建立催化劑的氧空位模型,如圖6中的O defect位置所示.氧空位的引入會使初始模型發(fā)生扭曲,導致空位附近的原子(Fe和O)出現一些顯著移動.計算得出的v表明形成氧空位缺陷的可行性,γ-Fe2O3(001)(ov)的氧空位形成能v為3.56eV,表示形成氧空位缺陷時,需要外界提供能量.后續(xù)將該模型作為氧缺陷的γ-Fe2O3(001)表面結構并作為計算對象.

圖6 氧缺陷的γ-Fe2O3表面結構

圖7 As2O3在氧缺陷的γ-Fe2O3表面的吸附構型

2.2.2 As2O3在O缺陷表面的吸附 在圖7中,構型(a)的結合形式幾乎與圖4的構型(a)相同,As2O3以O端吸附在O缺陷表面的Feoct位.構型(a)的吸附能為-0.88eV(表3),相較于在完美催化劑表面的吸附,吸附釋放的能量增加了0.13eV.在這種構型下,Mulliken電荷轉移為0.20e,表明在吸附過程中,一些電子從As2O3分子轉移到基底.在構型(b)中,As2O3以As端吸附在O缺陷表面的Feoct位時,相較于完美表面,其吸附能的絕對值降低.在構型(c)中,As2O3以O端吸附在O缺陷的空位處,出現了新的吸附形式,O原子最終與四面體Fetet原子鍵合,形成了鍵長為1.983?的O-Fe鍵.該吸附過程顯著提高了催化劑表面As2O3分子的吸附能力,其吸附能-1.83eV,電荷轉移為0.12e,屬于化學吸附過程.綜上可知,相對于完整催化劑表面,當催化劑γ-Fe2O3(001)表面存在O空位缺陷時,As2O3分子在γ-Fe2O3(001)表面的吸附能力增強.這與已有的研究結果一致:催化劑表面O空位缺陷的存在,有利于As2O3的吸附[43].

表3 As2O3在氧缺陷的γ-Fe2O3表面的優(yōu)化參數

2.3 M(Mo、Ti、Mg)摻雜對γ-Fe2O3(001)表面As2O3和NH3的吸附的影響

2.3.1 M(Mo、Ti、Mg)摻雜γ-Fe2O3(001)的表面結構 已有研究表明,助劑改性可以有效改善γ-Fe2O3抗砷中毒性能[19-20].由圖2可知,γ-Fe2O3(001)表面存在2種不同配位類型(八面體Feoct和四面體Fetet)的Fe原子.因此,為了構建摻雜的γ-Fe2O3(001)表面,分別使用不同的摻雜劑M(代表Mo、Ti、Mg)替換2種Fe原子(Feoct和Fetet).

如圖8所示,通過結構優(yōu)化得到了M-γFe2O3(001)表面模型,并計算了所有類型的M-γFe2O3(001)表面的摻雜形成能.分別對比Mo原子、Ti原子和Mg原子的不同位置的摻雜結構,發(fā)現在八面體位置的摻雜結構的比四面體位置更穩(wěn)定,Mo、Ti、Mg更傾向的摻雜位置均是Feoct位.在后續(xù)研究中,將以Mo、Ti和Mg原子摻雜在Feoct位的M-γFe2O3(001)進行分析.

圖8 摻雜改性后催化劑表面結構俯視圖

2.3.2 M(Mo、Ti、Mg)摻雜對As2O3吸附的影響 通過Mo、Ti、Mg原子摻雜,γ-Fe2O3(001)表面增加了新的活性位點,同時會影響相鄰Feoct位點對As2O3的吸附性能.基于上文的研究結果,As2O3更傾向于以其O端與γ-Fe2O3催化劑表面的Feoct原子鍵合.如圖9a和b所示,As2O3吸附在Mo位點和相鄰Feoct位點上,吸附能分別為-1.30和-0.75eV(表4).吸附過程的電荷轉移分別為0.22和0.19e.通過與圖4(a)結構γ-Fe2O3(001)表面的吸附比較,發(fā)現Mo的鄰近Feoct位對As2O3吸附能沒有變化,Mo位點對As2O3氣體吸附能大大提高,說明摻雜的Mo原子更容易與As2O3氣體分子結合,Mo的活躍性比Fe高.通過比較As2O3分子電荷轉移的變化,可以發(fā)現Mo點上獲得更多的電子,相互作用更加強烈.由于As2O3與Mo相互作用更強,As2O3更傾向于與Mo結合,從而達到保護活性組分Fe位點的作用.這與Li等[19]實驗結果一致:在γ-Fe2O3中引入Mo,發(fā)現Mo提高了脫硝活性,并阻止了活性鐵和砷的結合,As2O3對活性組分破壞減弱,對催化劑的抗砷中毒起很大的作用.

如圖9c和d所示,As2O3吸附在Ti位點和相鄰Feoct位點上,吸附能分別為-1.13和-0.71eV.在吸附過程的電荷轉移量分別為0.22和0.18e.相比之下,As2O3在摻Ti的γ-Fe2O3表面的吸附能明顯高于未摻Ti的表面,增加了0.38eV,并且電荷轉移量更多,As2O3的活化程度增強.在Ti摻雜體系中,Ti原子的相鄰Feoct位點對As2O3氣體的吸附能降低.這表明摻雜的Ti不僅可以優(yōu)先吸附As2O3,保護活性組分Fe位點,還抑制了活性Fe位點對As2O3的吸附.Ti原子摻雜改性可以提高γ-Fe2O3催化劑的抗砷性能.

如圖9e和f所示,As2O3吸附在Mg位點和相鄰Feoct位點上,吸附能分別為-1.09和-0.60eV.吸附過程的電荷轉移分別為0.09和0.18e.與Ti摻雜類似,通過Mg的摻雜,As2O3不僅更傾向于吸附在摻雜的Mg位點,從而保護活性Feoct位點,還降低了As2O3在Feoct位點的吸附能.與未摻雜的γ-Fe2O3相比,Mg摻雜對相鄰Feoct位As2O3的吸附能降低了0.15eV,降低效果更加明顯.通過Mg摻雜可以提升γ-Fe2O3抗砷能力.

綜合比較Mo、Ti、Mg對As2O3吸附的影響,發(fā)現As2O3吸附能順序:Mo-γFe2O3(Mo)>Ti- γFe2O3(Ti)>Mg-γFe2O3(Mg)>Mo-γFe2O3(Fe)>Ti- γFe2O3(Fe)>Mg-γFe2O3(Fe),顯然Mo的摻雜對As2O3吸附能力更強.綜上所述,Mo、Ti、Mg原子摻雜改性均可提高對As2O3的吸附能力,保護活性Fe位點,此外通過Ti和Mg摻雜還可以抑制活性Fe位點對As2O3的吸附.Mo、Ti和Mg均是良好的抗砷中毒摻雜助劑.

圖9 As2O3在M-γFe2O3(001)表面的優(yōu)化吸附模型

表4 As2O3在M-γFe2O3(001)表面吸附的優(yōu)化參數

2.3.3 M(Mo、Ti、Mg)摻雜對NH3吸附的影響 如圖10所示,首先,可以看出,NH3在γ-Fe2O3(001)表面的吸附形式是化學吸附,其吸附能為-1.01eV.還值得注意的是,在Mo、Ti、Mg原子摻雜以后,NH3在γ-Fe2O3(001)表面路易斯酸位點的吸附能值顯著增加,分別增加了0.5,0.25和0.13eV.其中,通過Mo的摻雜提高催化劑酸強度最強.當NH3在Mo、Ti、Mg摻雜的Lewis酸位點上的吸附能值增加,表明NH3的吸附受到促進,吸附強度增強,這將有利于催化劑上NH3的活化.因此可以得出結論,Mo、Ti、Mg的摻雜促進了催化劑表面酸性強度的增加,有利于NH3-SCR反應.

圖10 NH3在γ-Fe2O3及M摻雜表面的吸附構型

圖11 As2O3在2M-γFe2O3(001)表面的吸附模型

2.3.4 摻雜量對As2O3吸附的影響 如圖11a～d所示,As2O3以O端吸附在雙Mo原子摻雜位置不同的2Mo-γFe2O3表面的Mo位或Fe位,摻雜位置依次是兩個Mo原子分別摻雜在Fe的兩個八面體位、兩個四面體位以及一個四面體位和一個八面體位.對于結構(a),雙Mo原子摻雜在表面雙Feoct位,電荷轉移到催化劑表面0.21e(表5),As2O3分子吸附能為-1.45eV,相較于單Mo原子摻雜在Feoct位的吸附能增大了0.15eV,這表明通過增大Mo的摻雜量可以強化As2O3的吸附.對于結構(c和d),雙Mo原子摻雜在Feoct和Fetet位,As2O3在Mo位和相鄰Feoct位,吸附能分別為-1.20和-0.75eV,電荷轉移分別為0.24和0.17e,并與雙Mo原子摻雜在雙Feoct相比較,Mo位點吸附能明顯下降.

表5 As2O3在2M-γFe2O3(001)表面吸附的優(yōu)化參數

如圖11e～h所示,As2O3在不同結構的2Ti-γFe2O3表面的吸附能依次為-1.14,-0.69,-1.24和-0.71eV,電荷轉移分別為0.24,0.20,0.24,0.18e.對于結構(e),雙Ti原子摻雜在雙Feoct位,相較于未摻雜的γ-Fe2O3表面,提高了As2O3在催化劑表面的吸附能.但是,相較于單Ti原子摻雜,吸附能變化很少,增加了0.01eV.對于結構(f),雙Ti原子摻雜在雙Fetet位,降低了As2O3在催化劑表面活性Fe位的吸附能,有利于抗砷中毒.對于結構(g和h),雙Ti原子摻雜在Feoct和Fetet位,與單Ti原子摻比,砷在Ti位點吸附能增加,Fe位點吸附能下降.Ti原子摻雜量的增加,可以更好的保護活性Fe位點,提高抗砷性能.

如圖11i～l所示,As2O3在不同結構的2Mg- γFe2O3表面的吸附能依次為-0.94,-0.71,-1.08和-0.57eV,電荷轉移分別為0.08,0.15,0.09,0.19e.對于結構(i),雙Mg原子摻雜在催化劑表面的雙Feoct位,相較于單Mg原子摻雜在Feoct位的吸附能降低了0.15eV.對于結構(k和l),雙Mg原子摻雜在Feoct和Fetet位,As2O3會優(yōu)先吸附在活性Mg位并在相鄰Feoct位的吸附能出現明顯降低,保護活性組分.

綜上可知,增加Mo原子在Feoct位點的摻雜量,可以強化As2O3的吸附.通過增加Ti原子的摻雜量,可以增加As2O3在Ti位點吸附能,并降低在Fe位點吸附性能,可以更好的保護活性Fe位點,從而提高催化劑的抗As2O3性能.雙Mg原子摻雜在雙Feoct位,相較于單Mg原子摻雜在Feoct位的吸附能下降.

3 結論

3.1 As2O3分子在γ-Fe2O3表面的吸附既有物理吸附又有化學吸附.吸附能的不同與催化劑表面的吸附位點有關,Feoct原子是最佳吸附活性位.As2O3更傾向于以其O端吸附在γ-Fe2O3表面的Feoct位.

3.2 當γ-Fe2O3(001)表面存在O空位缺陷時,催化劑表面與As2O3分子之間存在明顯的電荷轉移和強相互作用,對As2O3分子的吸附能力得到了顯著增強.

3.3 Mo、Ti、Mg三種改性原子傾向的摻雜位置均是γ-Fe2O3(001)表面的Feoct位.通過Mo、Ti、Mg的摻雜均可提高對As2O3的吸附能力,從而保護活性Fe位點,并且Ti和Mg的摻雜還抑制了相鄰活性Fe位點對As2O3的吸附,此外增大Mo的摻雜量可以強化As2O3的吸附.通過Mo、Ti、Mg的摻雜可以增加NH3在γ-Fe2O3(001)表面的吸附,提高了催化劑表面酸強度,可以促進NH3-SCR反應.通過引入Mo、Ti、Mg對γ-Fe2O3催化劑提高抗砷中毒性能有很好的效果.

[1] 卞若愚,安忠義,李啟超,等.O3-NH3協(xié)同活性焦脫硫脫硝的均相預反應特性研究 [J]. 中國環(huán)境科學,2021,41(10):4476-4483.

Bian R U,An Z Y,Li Q C,et al. Characteristics of simultaneous removal of NOand SO2by O3-NH3synergy. [J]. China Environmental Science,2021,41(10):4476-4483.

[2] Shi J,Zhang Z H,Chen M X,et al. Promotion effect of tungsten and iron co-addition on the catalytic performance of MnO/TiO2for NH3-SCR of NO[J]. Fuel,2017,210:783-789.

[3] 蔣春來,宋曉暉,鐘悅之,等.2010～2015年中國燃煤電廠NO排放特征[J]. 中國環(huán)境科學,2018,38(8):2903-2910.

Jiang C L,Song X H,Zhong Y Z,et al. Characteristics of NOemissions from coal-fired power plants in China from 2010 to 2015 [J]. China Environmental Science,2018,38(8):2903-2910.

[4] Li Y H,Deng J L,Song W Y,et al. Nature of Cu species in Cu-SAPO-18catalyst for NH3-SCR: Combination of experiments and DFT calculations [J]. The Journal of Physical Chemistry C,2016,120(27):14669-14680.

[5] 付金艷,王振峰,白心蕊,等.γ-Al2O3酸性修飾稀土尾礦NH3-SCR脫硝性能[J]. 中國環(huán)境科學,2020,40(9):3741-3747.

Fu J Y,Wang Z F,Bai X R,et al. Denitration performance of NH3-SCR from γ-Al2O3acid modified rare earth tailings [J]. China Environmental Science,2020,40(9):3741-3747.

[6] 劉 晶,熊志波,周 飛,等.新型鈰鎢鈦復合氧化物催化還原脫硝機理 [J]. 中國環(huán)境科學,2018,38(5):1670-1676.

Liu J,Xiong Z B,Zhou F,et al. The NH3-SCR mechanism of a novel cerium-tungsten-titanium mixed oxide catalyst prepared through the hydrothermal co-precipitation method modified by H2O2complex [J]. China Environmental Science,2018,38(5):1670-1676.

[7] Chen H F,Xia Y,Huang H,et al. Highly dispersed surface active species of Mn/Ce/TiW catalysts for high performance at low temperature NH3-SCR [J]. Chemical Engineering Journal,2017,330: 1195-1202.

[8] 卿夢霞,劉 亮,尹子駿,等.商用V/W/Ti系脫硝催化劑表面SO3生成的反應機理[J]. 中國環(huán)境科學,2021,41(7):3161-3168.

Qing M X,Liu L,Yin Z J,et al. Generation mechanism of SO3on the surface of commercial V/W/Ti DeNOcatalysts [J]. China Environmental Science,2021,41(7):3161-3168.

[9] Liu C X,Yang S J,Ma L,et al. Comparison on the performance of α-Fe2O3and γ-Fe2O3for selective catalytic reduction of nitrogen oxides with ammonia [J]. Catalysis Letters,2013,143(7):697-704.

[10] 張信莉,王 棟,陳蓮芳,等.Mn摻雜對磁性γ-Fe2O3低溫SCR脫硝活性的影響[J]. 工程熱物理學報,2014,35:995-998.

Zhang X L,Wang D,Chen L F et al. Influence of Mn doping on magnetic γ-Fe2O3catalysts for selective catalytic reduction at low temperature [J]. Journal of Engineering Thermophysics,2014,35:995-998.

[11] Xu L T,Niu S L,Lu C M,et al. NH3-SCR performance and characterization over magnetic iron-magnesium mixed oxide catalysts [J]. Korean Journal of Chemical Engineering,2017,34:1576-1583.

[12] 張 月,王春波,劉慧敏,等.金屬氧化物吸附劑干法脫除氣相As2O3實驗研究[J]. 燃料化學學報,2015,43(4):476-482.

Zhang Y,Wang C B,Liu H M et al. Removal of gas-phase As2O3in dry process by metal oxide adsorbents [J]. Journal of Fuel Chemistry and Technology,2015,43(4):476-482.

[13] 龔泓宇,胡紅云,劉慧敏,等.燃煤過程中砷的遷移轉化及控制技術綜述[J]. 中國電機工程學報,2020,40(22):15.

Gong H Y,Hu H Y,Liu H M et al. Review of arsenic transformation and emission control during coal combustion [J]. Proceedings of the CSEE,2020,40(22):15.

[14] Ren S,Li S L,Su Z H,et al. Poisoning effects of KCl and As2O3on selective catalytic reduction of NO with NH3over Mn-Ce/AC catalysts at low temperature [J]. Chemical Engineering Journal,2018,351:540-547.

[15] Li X,Li J H,Peng Y,et al. Regeneration of commercial SCR Catalysts: Probing the existing forms of arsenic oxide [J]. Environmental Science & Technology,2015,49(16):9971-9978.

[16] Peng Y,Li J H,Si W Z,et al. Insight into deactivation of commercial SCR catalyst by arsenic: An experiment and DFT study [J]. Environmental Science & Technology,2014,48(23):13895-13900.

[17] Li X Y,Chen J,Xiao Y,et al. Insight into the homogenous and heterogeneous transformation behavior of arsenic on commercial V2O5-WO3-TiO2and novel γ-Fe2O3catalysts during selective catalytic reduction of NO[J]. Fuel,2021,301:121051.

[18] Peng Y,Si W Z,Li X,et al. Comparison of MoO3and WO3on arsenic poisoning V2O5/TiO2catalyst: DRIFTS and DFT study [J]. Applied Catalysis B: Environmental,2016,181:692-698.

[19] Li X Y,Chen J,Lu C M,et al. Performance of Mo modified γ-Fe2O3catalyst for selective catalytic reduction of NOwith ammonia: Presence of arsenic in flue gas [J]. Fuel,2021,294:120552.

[20] Li X Y,Chen J,Chen S Y,et al. Performance of Mg-Ti modified iron-based catalyst in NH3-SCR of NO at the presence of arsenic: Influence of oxygen and temperature [J]. Journal of Industrial and Engineering Chemistry,2021,101:387-396.

[21] Wu Y W,Zhou X Y,Mi T G,et al. Effect of WO3and MoO3doping on the interaction mechanism between arsenic oxide and V2O5-based SCR catalyst: A theoretical account [J]. Molecular Catalysis,2021,499:111317.

[22] Xing J Y,Wang C B,Si T,et al. Adsorption mechanism and competitive adsorption of As2O3and NH3molecules on CuO(111) surface: a DFT study [J]. Journal of Molecular Modeling,2021,27(6):178.

[23] Hu P B,Weng Q Y,Li D L,et al. Research on the removal of As2O3by γ-Al2O3adsorption based on density functional theory [J]. Chemosphere,2020,257:127243.

[24] J?rgensen J E,Mosegaard L,Thomsen L E,et al. Formation of γ-Fe2O3nanoparticles and vacancy ordering: An in situ X-ray powder diffraction study [J]. Journal of Solid State Chemistry,2007,180(1): 180-185.

[25] Jian W,Wang S P,Zhang H X,et al. Disentangling the role of oxygen vacancies on the surface of Fe3O4and γ-Fe2O3[J]. Inorganic Chemistry Frontiers,2019,6(10):2660-2666.

[26] Ren D D,Gui K T,Gu S C. Comparison of sulfur poisoning resistance of Ce/Mn doped γ-Fe2O3(0 0 1) surface in NH3-SCR reaction with DFT method [J]. Applied Surface Science,2021,561:149847.

[27] Xie C Y,Sun Y L,Zhu B Z. The promoting mechanism of doping Mn,Co,and Ce on gas adsorption property and anti-SO2oxidation over γ-Fe2O3(001) surface: A density functional theory study [J]. Colloids and Surfaces A: Physicochemical and Engineering Aspects,2021,628:127218.

[28] Guo P,Guo X,Zheng C G. Roles of γ-Fe2O3in fly ash for mercury removal: Results of density functional theory study [J]. Applied Surface Science,2010,256(23):6991-6996.

[29] Segall M D,Lindan P J D,Probert M J,et al. First-principles simulation: ideas,illustrations and the CASTEP code [J]. Journal of Physics: Condensed Matter,2002,14(11):2717-2744.

[30] Song Z J,Wang B,Yu J,et al. Effect of Ti doping on heterogeneous oxidation of NO over Fe3O4(1 1 1) surface by H2O2: A density functional study [J]. Chemical Engineering Journal,2018,354:517-524.

[31] Zhang Y,Liu J. Density functional theory study of arsenic adsorption on the Fe2O3(001) surface [J]. Energy & Fuels,2019,33(2):1414-1421.

[32] Bentarcurt Y L,Calatayud M,Klapp J,et al. Periodic density functional theory study of maghemite (001) surface. Structure and electronic properties [J]. Surface Science,2018,677:239-253.

[33] Grau-Crespo R,Al-Baitai A Y,Saadoune I,et al. Vacancy ordering and electronic structure of γ-Fe2O3(maghemite): a theoretical investigation [J]. Journal of Physics: Condensed Matter,2010,22 (25):255401.

[34] Da Hora G C A,Longo R L,Da Silva J B P. Calculations of structures and reaction energy profiles of As2O3and As4O6species by quantum chemical methods [J]. International journal of quantum chemistry,2012,112(20):3320-3324.

[35] Lyu Z K,Niu S L,Lu C M,et al. A density functional theory study on the selective catalytic reduction of NO by NH3reactivity of α-Fe2O3(0 0 1) catalyst doped by Mn,Ti,Cr and Ni [J]. Fuel,2020,267:117147.

[36] Maitarad P,Han J,Namuangruk S,et al. Theoretical guidance and experimental confirmation on catalytic tendency of M-CeO2(M = Zr,Mn,Ru or Cu) for NH3-SCR of NO [J]. Molecular simulation,2017,43(13-16):1240-1246.

[37] 楊 濤,曹 蕃,劉利軍,等.摻雜Ce/Zr對γ-Al2O3(110)表面的影響[J]. 燃燒科學與技術,2017,23(6):542-546.

Yang T,Cao F,Liu L J et al. Impact of Ce/Zr Doping on γ-Al2O3(110) surface [J]. Journal of Combustion Science and Technology,2017,23(6):542-546.

[38] Ren D D,Gui K T. Study of the adsorption of NH3and NOon the nano-γFe2O3(001) surface with density functional theory [J]. Applied Surface Science,2019,487:171-179.

[39] 周文波,牛勝利,王 棟,等.鈦改性對γ-Fe2O3選擇催化還原脫硝性能強化機制的分子模擬研究[J]. 燃料化學學報,2020,48(10): 1224-1235.

Zhou W B,Niu S L,Wang D,et al. Promoting effect of Ti in the Ti-modified γ-Fe2O3catalyst on its performance in the selective catalytic reduction of NO with ammonia,a DFT calculation study [J]. Journal of Fuel Chemistry and Technology,2020,48(10):1224-1235.

[40] Li Z P,Niu S L,Han K H,et al. Investigation into influences of methanol pre-adsorption on CaO(100) surface in transesterification for biodiesel production with molecular simulation [J]. Applied Catalysis A: General,2021,609:117908.

[41] Li F F,Shi C M,Wang X F,et al. The important role of oxygen defect for NO gas-sensing behavior of α-Fe2O3(0 0 1) surface: Predicted by density functional theory [J]. Computational Materials Science,2018,146:1-8.

[42] Feng Y C,Wang N N,Guo X. Density functional theory study on improved reactivity of alkali-doped Fe2O3oxygen carriers for chemical looping hydrogen production [J]. Fuel,2019,236:1057-1064.

[43] Zhang K H,Hu L T,Wang C F,et al. Middle-low-temperature oxidation and adsorption of arsenic from flue gas by Fe–Ce-based composite catalyst [J]. Chemosphere,2022,288:132425.

[44] Peng Y,Yu W W,Su W K,et al. An experimental and DFT study of the adsorption and oxidation of NH3on a CeO2catalyst modified by Fe,Mn,La and Y [J]. Catalysis Today,2015,242:300-307.

[45] Yang S J,Li J H,Wang C Z,et al. Fe-Ti spinel for the selective catalytic reduction of NO with NH3: Mechanism and structure–activity relationship [J]. Applied Catalysis B: Environmental,2012,117- 118:73-80.

[46] Man I C,Soriga S G,Parvulescu V. Effect of Ca and Sr in MgO(100) on the activation of methanol and methyl acetate [J]. Catalysis Today,2018,306:207-214.

致謝：本論文的科學計算得到了山東大學的高性能計算云平臺的計算支持和幫助.

Molecular simulation study on the anti-As2O3poisoning ability of γ-Fe2O3.

ZHOU Wen-bo1,NIU Sheng-li1*,LIU Si-tong1,WANG Dong2,HAN Kui-hua1,WANG Yong-zheng1

(1.Shandong Engineering Laboratory for High-efficiency Energy Conservation and Energy Storage Technology & Equipment,School of Energy and Power Engineering,Shandong University,Jinan 250061,China；2.Brook Byers Institute for Sustainable Systems and School of Civil and Environmental Engineering,Georgia Institute of Technology,Atlanta 30332,United States).,2022,42(8)：3600～3609

The γ-Fe2O3catalyst,due to its advantages in low cost and high catalytic performance,is thought to be a promising medium-low temperature SCR catalyst,but the As2O3in the flue gas likely becomes seriously deactivated. In this study,the density functional theory was used to characterize the adsorption of As2O3on the γ-Fe2O3surface as well as the mechanism of doping modification to improve the anti-As2O3poisoning performance. The adsorption properties of As2O3on the intact and O-deficient γ-Fe2O3(001) surfaces were examined,including adsorption site,adsorption structure,adsorption energy,PDOS,etc. At the same time,the catalyst model of γ-Fe2O3doped with Mo,Ti,and Mg was established to understand the mechanism of doping additives on improving the resistance to arsenic poisoning. The results show that the As2O3tends to be chemically adsorbed on Feoctsites on the γ-Fe2O3(001) surface with the O-terminus,and strong interaction and electron transfer occur during the adsorption process. When there are O defects on the surface,the adsorption energy of As2O3molecules increases. Mo,Ti,and Mg tend to be doped in Feoctsites,which thus enhances the adsorption capacity of As2O3. Increasing the doping amount of Mo can promote the adsorption of As2O3. As2O3tends to react with the more active Mo,Ti,and Mg,thereby protecting the active Fe sites from arsenic poisoning. The doping of Ti and Mg also inhibits the adsorption of As2O3on adjacent Fe sites. The doping of Mo,Ti,and Mg also promotes the adsorption of NH3on the catalyst surface and elevates the acidity of the surface,which is beneficial to the SCR reaction and to improving the anti-arsenic poisoning performance of the γ-Fe2O3catalyst.

As2O3；adsorption；γ-Fe2O3；O defect；Mo,Ti,and Mg doping；density functional theory (DFT)

X511

A

1000-6923(2022)08-3600-10

2022-01-19

山東省重大科技創(chuàng)新工程項目(2019JZZY020305)

* 責任作者,教授,nsl@sdu.edu.cn

周文波(1997-),男,山東青島人,山東大學碩士研究生,主要從事燃燒與污染物控制等領域研究.發(fā)表論文2篇.