Fe摻雜α-Bi2O3光電性質(zhì)的第一性原理研究

熊智慧，李志西，尹亞慶，蒲朝波，曾體賢，安辛友

(1.成都師范學(xué)院物理與工程技術(shù)學(xué)院，成都 611130; 2.西華師范大學(xué)物理與空間科學(xué)學(xué)院，南充 637002)

0 引 言

近年來，光催化半導(dǎo)體材料因其在環(huán)境污染治理、太陽能轉(zhuǎn)換等領(lǐng)域的巨大應(yīng)用前景而備受關(guān)注[1]。TiO2性能穩(wěn)定、催化活性高，是備受歡迎的光催化材料之一。但其禁帶寬度較大(約3.2 eV)[2]，只能吸收能量較高的紫外光才能響應(yīng)，對(duì)太陽光的利用率較低(僅為4%)，因此TiO2的實(shí)際應(yīng)用受到了很大程度的限制[3-5]。

Bi2O3是一種新型半導(dǎo)體材料，具有α、β、γ 和δ四種主要的晶型，其中α-Bi2O3在常溫常壓下最為穩(wěn)定，具有優(yōu)異的物理和化學(xué)性質(zhì)，是一種應(yīng)用前景良好的新型可見光催化劑[6-8]。但其禁帶寬度較寬，可見光吸收利用率低，電子-空穴對(duì)復(fù)合率高，可見光催化活性低[9]。研究發(fā)現(xiàn)，金屬離子摻雜可有效改善上述問題，進(jìn)而達(dá)到提高半導(dǎo)體光催化活性的目的。華中師范大學(xué)李玉潔課題組利用第一性原理研究了Cu摻雜α-Bi2O3[10]，在禁帶中引入了雜質(zhì)能級(jí)并減小了禁帶寬度，理論上證明了Cu摻雜可以提高α-Bi2O3光催化活性。哈爾濱理工大學(xué)Ding等[11]利用檸檬酸法成功地制備了Ag摻雜α-Bi2O3納米片，通過亞甲基藍(lán)氧化脫色對(duì)α-Bi2O3的光催化活性進(jìn)行觀察發(fā)現(xiàn)：與純?chǔ)?Bi2O3相比，Ag摻雜α-Bi2O3出現(xiàn)了紅移的現(xiàn)象，帶隙減小，從而說明了Ag摻雜增強(qiáng)了α-Bi2O3的光催化活性。蘭州理工大學(xué)戴劍鋒課題組以理論計(jì)算結(jié)果為出發(fā)點(diǎn)，制備得到Pr 摻雜α-Bi2O3的納米催化劑顆粒并對(duì)其可見光催化活性進(jìn)行了評(píng)價(jià)，Pr 摻雜α-Bi2O3表現(xiàn)出了更高的可見光催化活性[12]。內(nèi)蒙古工業(yè)大學(xué)張慶宏課題組[13]以鈦酸丁酯為鈦源，F(xiàn)e(NO3)3為鐵源，采用溶膠-凝膠、單軸靜電紡絲聯(lián)合技術(shù)成功制備了Fe/TiO2納米粉末，通過UV-Vis 漫反射光譜分析發(fā)現(xiàn)Fe摻雜TiO2后吸收閾值增大，吸收光譜發(fā)生明顯的紅移，帶隙變窄， TiO2的催化活性及其在可見光區(qū)的吸收得到了有效提高。北京師范大學(xué)Yin課題組采用溶劑熱合成法成功制備了低含量Fe摻雜Bi2O3[14]，在不同退火溫度下，F(xiàn)e摻雜Bi2O3的晶相由Bi2O2CO3轉(zhuǎn)變?yōu)閱涡宝?Bi2O3和立方Bi25FeO40，其α-Bi2O3與Bi25FeO40的混合物在可見光照射下表現(xiàn)出高的光催化活性。

由于α-Bi2O3光催化活性的影響因素較為復(fù)雜, 因此金屬離子摻雜α-Bi2O3的改性機(jī)理還有待進(jìn)一步明確。為了研發(fā)出新型高效的可見光催化劑，進(jìn)一步提高α-Bi2O3的可見光催化活性，本文對(duì)Fe摻雜α-Bi2O3體系進(jìn)行了第一性原理研究。通過計(jì)算摻雜體系的晶體結(jié)構(gòu)、電子結(jié)構(gòu)和光學(xué)性質(zhì)，從量子化學(xué)的角度對(duì)其光催化活性的變化機(jī)理進(jìn)行探討，為新型高效可見光催化材料的制備及其改性設(shè)計(jì)提供理論依據(jù)。

1 結(jié)構(gòu)模型和計(jì)算方法

1.1 結(jié)構(gòu)模型

α-Bi2O3單胞結(jié)構(gòu)中含有12個(gè)氧原子和8個(gè)鉍原子，如圖1所示。單斜α-Bi2O3空間群為P21/C，空間體積為0.329 35 nm3。α-Bi2O3晶格常數(shù)為a=0.584 4 nm，b=0.815 7 nm，c=0.750 3 nm，α=γ=90°，β=112.97°。為了減小邊緣效應(yīng)對(duì)體系能量的影響，本文采用如圖2所示的包含40個(gè)原子的2×1×1超胞模型計(jì)算研究Fe摻雜α-Bi2O3的電子結(jié)構(gòu)性質(zhì)，用一個(gè)Fe原子替代α-Bi2O3超胞中的一個(gè)Bi原子實(shí)現(xiàn)摻雜，摻雜的原子濃度為2.5%。

圖1 α-Bi2O3單胞結(jié)構(gòu)模型Fig.1 Model of α-Bi2O3 unit cell

圖2 Fe摻雜α-Bi2O3超胞模型Fig.2 Model of Fe doped α-Bi2O3 supercell

1.2 計(jì)算方法

本文基于密度泛函理論(DFT)中的HSE06雜化泛函方法計(jì)算了未摻雜α-Bi2O3和Fe摻雜α-Bi2O3的電子結(jié)構(gòu)和光學(xué)性質(zhì)，并通過能帶結(jié)構(gòu)研究分析其在可見光波段的催化活性及機(jī)理。計(jì)算過程中采用VASP(Viennaabinitiosimulation package)程序包對(duì)未摻雜和Fe摻雜α-Bi2O3的晶胞進(jìn)行幾何優(yōu)化。計(jì)算過程中的電子組態(tài)分別為Bi(6s26p3)、O(2s22p4)和Fe(3d64s2)，價(jià)電子與離子實(shí)之間的相互作用選用投影綴加平面波(PAW)勢(shì)來描述[15]。未摻雜α-Bi2O3體系Monkhorst-Pack K點(diǎn)網(wǎng)格[16]取為6×4×6，截?cái)嗄茉O(shè)置為500 eV；Fe摻雜α-Bi2O3體系，K點(diǎn)網(wǎng)格采用3×2×3，截?cái)嗄鼙３?00 eV不變。為了確保體系能量得到精確收斂，自洽計(jì)算體系總能收斂小于5×10-5eV/atom, 原子的公差偏移小于5×10-5nm，晶體內(nèi)應(yīng)力偏差小于0.02 GPa，殘余應(yīng)力低于0.1 eV/nm。

2 結(jié)果與討論

2.1 晶體結(jié)構(gòu)參數(shù)

表1為幾何優(yōu)化后本征α-Bi2O3和Fe摻雜α-Bi2O3的晶格常數(shù)。本文計(jì)算得到本征α-Bi2O3晶格常數(shù)與實(shí)驗(yàn)數(shù)據(jù)接近[17]，誤差較小，說明本文所建α-Bi2O3模型有效，參數(shù)設(shè)置較為合理。與本征α-Bi2O3相比，F(xiàn)e原子摻雜使得體系的體積減小，但其晶格常數(shù)c/a變化較小，說明Fe摻雜后超晶胞體系結(jié)構(gòu)穩(wěn)定。由于Fe的離子半徑(0.055 nm)小于Bi離子半徑(0.103 nm)，F(xiàn)e-O鍵的鍵長(0.187 4 nm)小于Bi-O鍵的鍵長(0.218 0 nm)，根據(jù)量子化學(xué)理論可知晶胞體積應(yīng)當(dāng)減小，結(jié)果表明Fe摻雜對(duì)α-Bi2O3的結(jié)構(gòu)影響較小，為進(jìn)一步實(shí)驗(yàn)研究提供了理論依據(jù)。

表1 幾何優(yōu)化后本征α-Bi2O3和Fe摻雜后的晶格常數(shù)Table 1 Lattice constants of instrinsic and Fe doped α-Bi2O3 after structure optimized

2.2 電子結(jié)構(gòu)性質(zhì)

2.2.1 本征α-Bi2O3態(tài)密度

圖3為本征α-Bi2O3的總態(tài)密度和分態(tài)密度圖，費(fèi)米能級(jí)位于價(jià)帶頂0 eV處，由圖3可知，2.69～7.1 eV為導(dǎo)帶，主要由Bi6p態(tài)電子軌道組成，O2p態(tài)電子軌道有少量貢獻(xiàn)；-5.7～0 eV為價(jià)帶，主要由O2p和Bi6p態(tài)電子軌道組成，Bi6s態(tài)電子軌道有少量貢獻(xiàn)。本征α-Bi2O3存在一個(gè)較大的禁帶寬度Eg，其值約為2.69 eV，與實(shí)驗(yàn)測(cè)量值2.75 eV[18]非常接近，表明采用HSE06雜化泛函方法可以精確計(jì)算體系的能量分布。

2.2.2 Fe摻雜α-Bi2O3態(tài)密度

圖4為Fe摻雜α-Bi2O3總態(tài)密度和分態(tài)密度圖。由圖4可知，F(xiàn)e摻雜后禁帶寬度為2.34 eV，該值比本征α-Bi2O3的禁帶寬度減小了0.35 eV，禁帶寬度的變化源于體系摻雜后電子結(jié)構(gòu)的變化。Fe摻雜α-Bi2O3后在價(jià)帶頂附近引入了Fe3d軌道雜質(zhì)能級(jí)，F(xiàn)e3d軌道與O2p軌道在價(jià)帶頂附近雜化， 使得價(jià)帶頂能級(jí)向高能量方向發(fā)生極微弱移動(dòng)，與未摻雜α-Bi2O3相比，其移動(dòng)量為0.02 eV，且Fe4s態(tài)和3d態(tài)對(duì)導(dǎo)帶底幾乎沒有貢獻(xiàn)。導(dǎo)帶底由O2p態(tài)和Bi6p態(tài)共同決定，但Fe摻雜后它們的能帶分布向低能區(qū)移動(dòng)，即導(dǎo)帶向低能區(qū)移動(dòng)，與未摻雜α-Bi2O3相比，其移動(dòng)量為0.33 eV。因此，禁帶寬度變窄，有助于光電子激發(fā)躍遷。

圖3 本征α-Bi2O3態(tài)密度Fig.3 Densities of states of instrinsic α-Bi2O3

圖4 Fe摻雜α-Bi2O3態(tài)密度Fig.4 Densities of states of Fe doped α-Bi2O3

2.3 光學(xué)性質(zhì)

凝聚態(tài)物理中通常用復(fù)介電函數(shù)ε(ω)=ε1(ω)+iε2(ω)來系統(tǒng)描述固體宏觀光學(xué)響應(yīng)。體系光子-電子間相互作用可以用基態(tài)電子時(shí)間相關(guān)微繞來計(jì)算。通過占據(jù)態(tài)和非占據(jù)態(tài)波函數(shù)的動(dòng)量矩陣元和Kramers-Kronig關(guān)系可以分別得到介電虛部ε2(ω)和介電實(shí)部ε1(ω)[19-21]：

(1)

(2)

其中，e為電子電荷，m為電子質(zhì)量，ω為光子入射光子頻率，ω′為高頻介電常數(shù)，M為動(dòng)量躍遷矩陣，i和j分別為躍遷初態(tài)和終態(tài)，fi為第i態(tài)的費(fèi)米分布函數(shù)，Ei為第i態(tài)電子能量，k為波矢，P代表主值積分。光學(xué)吸收系數(shù)I(ω)和電子能量損失譜L(ω)可以通過介電函數(shù)的虛部和實(shí)部計(jì)算得到，關(guān)系式如下所示[19-21]：

(3)

(4)

圖5為本征及Fe摻雜α-Bi2O3的復(fù)介電函數(shù)虛部曲線。由圖可知，本征和Fe摻雜α-Bi2O3的介電響應(yīng)位于紫外區(qū)域，但Fe的摻入有效增強(qiáng)了可見光波段的介電響應(yīng)。圖6為本征及Fe摻雜α-Bi2O3的光學(xué)吸收譜，選取的波長范圍為300～800 nm。本征α-Bi2O3的吸收帶邊在550 nm附近，與實(shí)驗(yàn)結(jié)果一致。Fe摻雜α-Bi2O3后，吸收帶邊位于650 nm左右，且有更大的光吸收系數(shù)，說明Fe摻雜α-Bi2O3能有效提高其可見光吸收效率。在650 nm以后本征α-Bi2O3幾乎沒有吸收，而Fe摻雜體系仍然具有一定的光吸收能力，且吸收邊緣明顯延伸到紅外波段，即發(fā)生了紅移。上述結(jié)果與本文Fe摻雜后體系的禁帶寬度減小相對(duì)應(yīng)，說明Fe摻雜α-Bi2O3后光吸收范圍更大，有利于提高 α-Bi2O3的光催化性能。

圖5 本征及Fe摻雜α-Bi2O3的介電函數(shù)譜Fig.5 Dielectric function curves of instrinsic and Fe doped α-Bi2O3

圖6 本征及Fe摻雜α-Bi2O3的光吸收曲線Fig.6 Optical absorption spectra of instrinsic and Fe doped α-Bi2O3

3 結(jié) 論

本文基于密度泛函理論對(duì)本征及Fe摻雜α-Bi2O3體系進(jìn)行結(jié)構(gòu)優(yōu)化，并采用雜化泛函 HSE06方法對(duì)它們的電子結(jié)構(gòu)和光學(xué)性質(zhì)進(jìn)行了研究。結(jié)果表明：Fe摻雜α-Bi2O3對(duì)其結(jié)構(gòu)影響較小；本征α-Bi2O3的禁帶寬度為2.69 eV，F(xiàn)e摻雜后禁帶寬度減小了0.35 eV，其值為2.34 eV；對(duì)其電子結(jié)構(gòu)研究表明，F(xiàn)e摻雜α-Bi2O3后，F(xiàn)e的3d軌道能量低于Bi的6p軌道的能量使得導(dǎo)帶下移、禁帶寬度減小；對(duì)其光學(xué)性質(zhì)研究得出Fe摻雜α-Bi2O3的光吸收帶發(fā)生紅移，擴(kuò)大了光吸收范圍，提高了光吸收效率，進(jìn)而使可見光催化活性得到了增強(qiáng)。