本征缺陷對β-Ga2O3光催化性質(zhì)影響的第一性原理研究

張靜林， 李家印， 張 龍， 楊蓮紅， 朱應(yīng)濤

(1.昌吉學(xué)院 計算機工程系， 昌吉 831100； 2. 昌吉學(xué)院 物理系， 昌吉 831100)

1 引 言

氫氣作為解決能源危機和環(huán)境污染的終極能源之一，受到科研人員的廣泛關(guān)注. 其中，借助太陽能的光催化技術(shù)是制備氫氣的一種重要手段. 傳統(tǒng)半導(dǎo)體如TiO2，SrTiO3，因其低成本，熱、化學(xué)、光穩(wěn)定性高和無毒等性質(zhì)，成為光催化材料研究的熱點. 但是其寬帶隙(大約3.0 eV)，只能吸收占太陽光總能量4%的紫外光，光吸收效率低，另外光生電子空穴對具有較高的復(fù)合率，因此這些因素限制了該類半導(dǎo)體在光催化領(lǐng)域的實際應(yīng)用. 為了提高光催化的效率，人們主要從兩方面的工作展開：1)對傳統(tǒng)半導(dǎo)體進行電子結(jié)構(gòu)的調(diào)控，以提高光吸收率和光生電子空穴對分離率(比如，摻雜金屬或者非金屬，構(gòu)建異質(zhì)結(jié)等)；2)尋找一種新型具有可見光吸收和電子空穴復(fù)合率低的半導(dǎo)體.

Ga2O3，其穩(wěn)定相是β相，其他相均為亞穩(wěn)相，在溫度高時都會轉(zhuǎn)化成β相[1]. 空間群是C2/m，其帶隙在4.2-4.9 eV，本征激發(fā)只能吸收紫外光，并且是深紫外，因此該材料的研究主要集中在紫外探測器、氣體傳感器等. Ga2O3作為光催化材料也引起了人們極大的興趣，但因為其寬的帶隙，理想的Ga2O3光催化材料只對UV-C光有響應(yīng). 而Binet等人[2]報道了β-Ga2O3熒光譜中出現(xiàn)藍光的峰，并將此歸因于本征缺陷(O空位)，這也為拓展光吸收的光催化研究提供了可能. 由于Ga的d10的電子組態(tài)的特點[3, 4]，導(dǎo)帶中應(yīng)該是Ga原子的s或者p軌道，因此應(yīng)該具有很好的光催化活性分解水產(chǎn)生H2和O2，Takashi等人[5]報道了Ga2O3和NiO負載的Ga2O3在深紫外區(qū)域分解水，并且解釋NiO的負載是有必要的. 除了光催化分解水之外，有毒的有機廢物處理也受到廣泛的關(guān)注，Hou等人[6]發(fā)現(xiàn)多孔β-Ga2O3在紫外光區(qū)光催化分解苯的活性優(yōu)于商業(yè)化的TiO2P25，可在空氣氣氛下將苯及其衍生物分解成為CO2，并且具有較好的穩(wěn)定性. 隨后他們又研究了α、β、γ相Ga2O3具有分解有機易揮發(fā)有機物的能力，得出β相Ga2O3的活性優(yōu)于其他兩相，并且從實驗上解釋了不同相光催化活性的差異源于不同的幾何和電子結(jié)構(gòu). 由于β相Ga2O3導(dǎo)帶中電子高的還原性質(zhì)，Zhao等人[7]發(fā)現(xiàn)其可以有效降解全氟辛酸并且氟離子同時進入溶液中. 另外，β相Ga2O3可以UV-C選擇性分解甲苯成為CO2，這是TiO2做不到的，作者將光催化活性高的原因歸納為高的比表面積、小的晶粒、結(jié)晶度高和定域的內(nèi)電場等. 作為一種寬帶隙的光催化劑，Ga2O3在紫外區(qū)表現(xiàn)出較高的光催化活性，太陽能的利用率不高，因此如何拓展Ga2O3的光吸收范圍也引起了人們的關(guān)注. Ga2O3與中等帶隙的In2O3固溶體[8]表現(xiàn)出高效光催化分解水的能力，可以歸因于帶隙的減小和帶邊位置的調(diào)節(jié). 摻雜也是提高光催化活性的一個重要手段，Ryou等人[9]利用水熱法制備了Sn-摻雜的Ga2O3納米結(jié)構(gòu)，表現(xiàn)出高的光催化活性，并將其解釋為Sn作為在禁帶中捕獲點提高了電子和空穴對的分離. 對于非金屬摻雜，N[10]和Si[11]的摻雜導(dǎo)致Ga2O3的吸收邊紅移，從而提高了其光催化的活性. Zhang等人[12]制備了Ga2O3納米帶，并指出在熒光光譜中有藍光的發(fā)射峰，并將原因解釋為在納米帶的制備過程中產(chǎn)生的Ga缺陷、O缺陷導(dǎo)致的電子結(jié)構(gòu)的改變. 納米結(jié)構(gòu)材料出現(xiàn)的缺陷態(tài)和帶邊態(tài)可能使得光吸收譜發(fā)生紅移，因此人們做了大量工作去研究納米結(jié)構(gòu)的Ga2O3(比如，納米線、納米帶、納米片)的相關(guān)性質(zhì)[13-17]. 其中，Tien等人[18]制備了β相Ga2O3納米帶，光催化活性提高歸因于缺陷出現(xiàn)導(dǎo)致的光吸收效率的改善. 理論上，人們對Ga2O3的幾何結(jié)構(gòu)、電子結(jié)構(gòu)、力學(xué)性質(zhì)、振動光譜及熱力學(xué)性質(zhì)做了大量的研究[19-22]，對于β-Ga2O3及其本征缺陷引起光催化性質(zhì)改變機理的解釋較少. 因此，本文有必要通過第一性原理去計算含本征缺陷的Ga2O3的幾何和電子性質(zhì)，解釋缺陷對β-Ga2O3光催化性質(zhì)的影響，闡釋光催化機理.

2 計算方法和細節(jié)

本文采用密度泛函的第一性原理方法計算了Ga2O3的相關(guān)幾何結(jié)構(gòu)和電子性質(zhì)，所有的計算是在VASP(Vienna Ab initial Simulation Package)和CRYSTAL14程序中完成的. 鑒于傳統(tǒng)廣義梯度算法(GGA)低估禁帶寬度的缺點，模擬中采用雜化泛函(比如B3LYP等)可以得到符合實驗的禁帶寬度. 而基于平面波的程序采用傳統(tǒng)的GGA方法在幾何優(yōu)化中存在優(yōu)勢，可以得到接近實驗值的幾何結(jié)構(gòu)參數(shù)，但是采用雜化泛函計算將非常的耗時，而CRYSTAL14采用高斯基組，在計算雜化泛函中的Fock部分相對較為容易，并且在考慮對稱性的情況計算耗時較少，因此可以迅速的得到體系的電子結(jié)構(gòu)，但是其對于缺少對稱性的體系優(yōu)化幾何結(jié)構(gòu)較為困難. 因此，該工作中采用傳統(tǒng)GGA利用VASP軟件進行含缺陷體系的(Ga空位和O空位)Ga2O3超胞的幾何結(jié)構(gòu)優(yōu)化，在幾何優(yōu)化的基礎(chǔ)上，采用雜化泛函基于CRYSTAL14軟件進行電子結(jié)構(gòu)等的計算. 下面將對兩種方法采用的參數(shù)進行說明：

1)VASP，采用綴加平面波贗勢方法，采用PBE(Perdew-Burke-Ernzerhof)交換關(guān)聯(lián)函數(shù)和超軟贗勢. 平面波作為波函數(shù)展開的基組，平面波的數(shù)量是通過動量截斷值(Ecutoff)來控制的，此Ecutoff值為500 eV. 對于幾何優(yōu)化，將所有原子上的力作為收斂的判據(jù)，設(shè)置為0.01 eV/?，對于電子結(jié)構(gòu)計算結(jié)束的標志是當能量的變化值ΔE≤1.0×10-5eV. 采用Monkhorst-Pack方案在布里淵區(qū)取點實現(xiàn)關(guān)于倒格矢空間的積分，對于β-Ga2O3塊體網(wǎng)格取點大小為4×16×8，不包含Γ點. 對于含有Ga和O空位體系采用1×2×1的超胞(含有40個原子)來計算，結(jié)構(gòu)如圖1所示. 而對于空位缺陷，在平面波的計算就是去除原子Ga或O.

圖1 β-Ga2O3幾何結(jié)構(gòu)圖，藍色大球和紅色小球分別代表的是Ga原子和O原子. 圖中標出了Ga的兩種配位方式，與O形成四面體和八面體結(jié)構(gòu). 結(jié)構(gòu)圖由VESTA軟件繪出.Fig. 1 Geometric structure of β-Ga2O3, the large blue and small red spheres represent Ga and O atoms respectively. Two coordinated types of Ga are shown in the figure, forming tetrahedral and octahedral structures with O. The structure diagram is drawn by VESTA software.

2)CRYSTAL14，采用雜化泛函交換關(guān)聯(lián)函數(shù)(B3LPY)，用全電子高斯函數(shù)線性組合單電子的晶體軌道，基組分別為：Ga_pob_DZVP_rev2[511s-31p-1d]和O_pob_DZVP_rev2[63311s-5311p-41d][23]. 實空間的庫侖和交換積分計算使用五個參數(shù)控制(TOL1-TOL5)，本文的參數(shù)值分別為10-8，10-8，10-8，10-8和10-16. 倒格矢空間取樣對于塊體采用8×8×8網(wǎng)格，對于1×2×1超胞采用4×4×4網(wǎng)格塊體的晶格幾何優(yōu)化采用BFGS(Broyden-Fletcher-Goldfarb-Shanno)方案，收斂標準設(shè)置力的最大值和均方根分別小于4.5×10-4和3×10-4a.u.，而位移的最大值和均方根小于1.8×10-3和1.2×10-3，為CRYSTAL14幾何優(yōu)化的默認值，電子結(jié)構(gòu)計算能量收斂參數(shù)設(shè)置為10-9a.u..

3 結(jié)果與討論

3.1 理想的Ga2O3

首先，為了驗證方法的可靠性并與含缺陷體系對比，本文計算了理想的β-Ga2O3幾何結(jié)構(gòu)和相應(yīng)的性質(zhì)，部分結(jié)果列于表1中. 理想的β-Ga2O3具有空間群C2/m，是單斜晶體. 晶格參數(shù)是a=12.42 ?，b=3.07 ?，c=5.86 ?，和β=103.74°(CRYSTAL14計算結(jié)果)，與實驗值[24, 25]和其他理論值[19-21]非常接近. 利用Brich-Murnaghan方程對E-V關(guān)系進行了擬合(如圖2所示)，得到了理想Ga2O3的體模量為171.75 GPa，與其他的理論計算相符合[22]，因此其具有優(yōu)良的力學(xué)性能. 在β-Ga2O3中，四配位的Ga與O形成四面體結(jié)構(gòu)和六配位的Ga與O形成八面體結(jié)構(gòu)，優(yōu)化的Ga4c-O鍵長分別為1.85和1.89 ?，優(yōu)化的Ga6c-O鍵長分別為1.96、2.04和2.10 ?. 因此，Ga和O并沒有形成正四面體和正八面體結(jié)構(gòu)，其中八面體結(jié)構(gòu)中的鍵長稍長于四面體的Ga-O，可以推測八面體中Ga和O之間的相互作用略弱于四面體，這也在Mulliken布居分析中得到驗證，Ga4c靜電荷為1.650|e|，相比于Ga6c(1.615|e|)，四配位的Ga失去更多的電子給O原子，Ga和O的離子鍵應(yīng)該更強，還有Ga4c與O之間的重疊布居也大于Ga6c和O，說明了四面體中的Ga-O結(jié)合強于八面體. 圖3給出了理想的Ga2O3的能帶結(jié)構(gòu)，導(dǎo)帶底(conduction band minimum(CBM))位于Γ點，其價帶頂(valence band maximum(VBM))出現(xiàn)在D和E之間的某個點(不過該點和Γ點的能量相差約0.02 eV)，因此屬于間接帶隙結(jié)構(gòu)，禁帶寬度為4.73 eV，和其他的理論值相符合，也驗證了其在深紫外區(qū)光作為探測器的可能性. 從能帶圖可以看出，VBM附近能帶較為平坦，CBM附近的能帶較為色散，表明其具有較低導(dǎo)帶電子有效質(zhì)量和較高的價帶空穴的有效質(zhì)量，電子和空穴質(zhì)量較大的差異可能導(dǎo)致載流子的復(fù)合率降低，在不考慮光吸收的情況下，β-Ga2O3或許在光催化領(lǐng)域也有潛在的應(yīng)用. 從總態(tài)密度(TDOS)和投影態(tài)密度(PDOS)(如圖3所示)可以看出，價帶頂附近主要由O的2p軌道組成，有少量Ga的s、p和d出現(xiàn)在價帶比較深的位置，而CBM附近主要由Ga的s軌道組成，由于s軌道的離域性，因此在能帶上表現(xiàn)出較為色散的特征，Ga的d軌道由于d10的特點，沒有電子得失，出現(xiàn)遠離VBM的位置，并且表現(xiàn)定域的特點，而鑒于變形的四面體和八面體結(jié)構(gòu)，d軌道并沒有出現(xiàn)明顯t2g和eg的分裂，唯一例外是在八面體結(jié)構(gòu)中有dx2-y2和dz2有混合，但是沒有明顯三下二下的分裂.

圖2 β-Ga2O3的能量Energy(eV)和體積V(?3)的關(guān)系.Fig.2 Relationship between energy (eV) and volume V (?3) of β-Ga2O3

圖3 β-Ga2O3能帶結(jié)構(gòu)(a)和總態(tài)密度(TDOS)及投影態(tài)密度(PDOS)(b)，紅色點線代表費米能級，并設(shè)置為能量的零點.Fig. 3 The band structure (a) and total density of states (TDOS) and projected density of states (PDOS) (b) of β- Ga2O3. The red dot line represents Fermi level which is set as the zero point of energy.

3.2 含Ga空位缺陷的Ga2O3

為了確定Ga和O空位缺陷的穩(wěn)定性，本文計算了缺陷形成能Ef，計算公式如下：

Ef=Edef-Ebulk+μGa(μO)

其中Edef是含Ga或O空位的超胞和理想Ga2O3的總能，μGa，μO是Ga和O元素的化學(xué)勢，其中μO在富氧的條件下是通過計算O2的總能得到，即，μO= 0.5EO2，對于貧氧的條件，即富Ga情況，uO=(EGa2O3-2EGa)/3，其中EGa2O3是理想的Ga2O3的總能量，EGa是bulk的金屬Ga的總能除以原胞中Ga原子的數(shù)目，對于μGa反之亦然. 雖然四配位的Ga與O的相互作用比較強，但是由于配位數(shù)低，構(gòu)造空位缺陷需要斷鍵少，此處選擇移除四配位Ga原子來模擬，計算的Ga空位缺陷形成能在富氧下是4.32 eV，而貧氧下為9.00 eV，O空位在富氧下為4.04 eV，而貧氧下為0.92 eV，可以看出無論在富氧或者貧氧的情況，O空位的缺陷形成能要低于Ga空位的缺陷形成能，缺陷形成能Ef越低缺陷體系越穩(wěn)定，因此，推測O空位形成的可能性更大一些，由于O空位的形成，Ga2O3可能會表現(xiàn)出N型半導(dǎo)體的特征.

為了研究Ga空位缺陷對于Ga2O3體系的影響，本文計算了Ga空位體系幾何結(jié)構(gòu)和電子結(jié)構(gòu).Ga空位的出現(xiàn)雖然使得晶格參數(shù)有所減小，但是并沒有引起很大的改變，最大的改變沿c軸，僅為0.09 ?. 對于Ga-O鍵長，只考慮了空位附近的原子的弛豫，鍵長有相應(yīng)的縮小，并沒有引起很大的變形. 因此，Ga空位的出現(xiàn)并未導(dǎo)致晶格參數(shù)和坐標明顯的改變. 電子結(jié)構(gòu)方面，本文計算了含Ga缺陷Ga2O3的Mulliken電荷、能帶結(jié)構(gòu)和態(tài)密度. 從Mulliken電荷布居分析看出，含Ga空位的Ga2O3中Ga原子的Mulliken電荷數(shù)目相較于理想情況增大，這個很容易理解是由于Ga原子的失去，多余的電子是由其他Ga原子提供. 重疊布居分析可以得出，某些Ga-O之間重疊布居增加，說明Ga-O之間的作用也有所增強，這可以解釋某些鍵長的縮短. 含Ga空位的Ga2O3表現(xiàn)出順磁性的特征，非成對電子表現(xiàn)出O 2p的特征. 圖4給出了含Ga空位缺陷的Ga2O3能帶結(jié)構(gòu)，對于α自旋，CBM仍然在Γ點，VBM在Y點，保持了原來的直接帶隙的特點. CBM附近能帶比較色散，VBM附近能帶比較定域，禁帶寬度為5.08 eV，其較理想的體系有所藍移，在禁帶中沒有出現(xiàn)雜質(zhì)能級，而VBM相對于β自旋有所下降，這也是禁帶寬度增加的原因. 對于β自旋，在禁帶中出現(xiàn)了三條缺陷能帶，并且均為空帶，都未被電子占據(jù)，因此電子可能從Ga2O3的VBM躍遷到缺陷能級，電子躍遷能為1.93 eV，相對于理想體系的4.73 eV的禁帶寬度，有顯著的減小，因此空位缺陷的出現(xiàn)可能會導(dǎo)致材料的光吸收向可見光拓展，從而提高了光吸收效率. 另外，本文計算了含Ga的Ga2O3的總的能態(tài)密度(TDOS)和投影態(tài)密度(PDOS)，如圖4(b)所示. 從DOS和PDOS圖上可以看出，和理想的Ga2O3一樣，含Ga缺陷的Ga2O3的VBM只要是由O的2p軌道組成，而CBM是由Ga的4s軌道組成，因此α自旋能帶保持了理想能帶結(jié)構(gòu)的色散特征. 只是α自旋的能帶O 2p較β向低能區(qū)下移. 對于β自旋能帶，Ga空位導(dǎo)致了禁帶中缺陷能態(tài)的產(chǎn)生，并且這些缺陷能態(tài)主要是O 2p軌道組成. 因此，電子躍遷就是從O 2p到O 2p，由于宇稱的相似性，躍遷幾率可能不會太大，另外，O 2p之間的電子躍遷可能會導(dǎo)致電子空穴容易復(fù)合，不利于光催化反應(yīng)的進行.

表1 計算的理想、Ga缺陷、O缺陷的Ga2O3幾何結(jié)構(gòu)(晶格參數(shù)、鍵長)和帶隙，理想Ga2O3的體模量和Ga缺陷、O缺陷的Ga2O3的缺陷形成能.

3.3 含O空位缺陷的Ga2O3

為構(gòu)建O空位缺陷體系，考慮斷鍵的數(shù)目，本文選擇移除低配位O. 優(yōu)化的晶格參數(shù)a=12.43 ?，b=6.14 ?，c=5.83 ?，β=103.96°，和理想的Ga2O3相比基本沒有改變. 為了觀測O空位缺陷對局部幾何結(jié)構(gòu)影響，首先考慮缺陷附近的Ga-O鍵長的情況，Ga4c-O鍵長為1.94和2.18 ?，相對于理想的情況有所伸長，最大的伸長量為0.29 ?，而Ga6c-O鍵長在1.94-2.04?之間，和理想的體系的鍵長改變不大，而其他位置基本上沒明顯的改變. 因此，和Ga空位一樣，缺陷并沒有引起幾何結(jié)構(gòu)明顯的改變. 由于O原子的失去，會產(chǎn)生兩個多余的電子，兩個電子可能就會遷移到缺陷附近的Ga上，從Mulliken電荷分析可以看出，O空位附近的四配位和六配位Ga的凈電荷分別為0.840|e|和0.138|e|，也驗證了關(guān)于電子轉(zhuǎn)移的猜測. 其他Ga的凈電荷接近于理想Ga2O3. 為了研究電子躍遷和光吸收情況，本文計算了含O空位的Ga2O3的能帶結(jié)構(gòu)和態(tài)密度及投影態(tài)密度(如圖5所示)，從能帶結(jié)構(gòu)圖上可以看出，因為O空位的引入，在Ga2O3出現(xiàn)了缺陷能級，并且導(dǎo)帶邊下移，此處將費米能級視為VBM，即VBM在E和C點之間，而CBM在Γ點，因此保持了間接帶隙的特點，光躍遷的能量為2.92 eV，這也能解釋熒光實驗觀察到藍光峰的原因. 對于光吸收，Γ→Γ直接躍遷的帶隙是4.30 eV，相對于理想的4.73 eV，有所減小，光吸收有紅移. 和理想的Ga2O3一樣，VBM附近的能帶保持了定域的特點，而CBM附近的能帶仍然離域，載流子的有效質(zhì)量有較大的差異，也降低了載流子在遷移過程中的復(fù)合. 另外從DOS圖(如圖5所示)可以看出，CBM主要是Ga的4s和4p軌道占據(jù)，而價帶主要是O的2p，少量的Ga 4s和4p軌道組成. 而缺陷能帶(VBM附近的能帶)主要是Ga 3d，4s，4p和O 2p軌道組成，表明有部分電荷轉(zhuǎn)移到Ga的4s和4p軌道，這也驗證了O空位出現(xiàn)多余的電子轉(zhuǎn)移到了Ga離子上. 和含Ga空位的體系不同，電子躍遷是從缺陷態(tài)(d，s，p)到導(dǎo)帶(4s，4p)在宇稱上是被允許的. 因此，由于間接帶隙、電子和空穴有效質(zhì)量的差異和減小的躍遷能的特點，含O空位的Ga2O3光吸收的效率提高和光生載流子復(fù)合率降低，從而有可能提高其光催化的效率.

圖4 含Ga空位Ga2O3的能帶結(jié)構(gòu)(a)和態(tài)密度及投影態(tài)密度(b)，(a)左圖對應(yīng)α自旋，(a)右圖對應(yīng)β自旋. 紅色點線代表費米能級，并設(shè)置為能量的零點.Fig.4 The band structure (a) and thedensity of states and projected density of states (b) of Ga2O3 with Ga vacancy, the left figure in (a) corresponds to α spin, the right figure in (a) corresponds to β spin. The red dot line represents the Fermi level which is set as the zero point of the energy.

圖5 含O缺陷的Ga2O3的能帶結(jié)構(gòu)(a)和態(tài)密度及投影態(tài)密度(b). 紅色點線代表費米能級，并設(shè)置為能量的零點.Fig. 5 The band structure (a), density of states (DOS) and projected density of states (b) of Ga2O3 with O-defects. The red dots line represents Fermi level which is set as the zero point of energy.

4 結(jié) 論

基于密度泛函理論，采用標準(PBE)和雜化泛函(B3LYP)，利用VASP和CRYSTAL14軟件包計算了理想、含Ga空位和含O空位的幾何結(jié)構(gòu)和電子結(jié)構(gòu). 計算的理想的Ga2O3的晶格參數(shù)和實驗相近，驗證了計算的可靠性. 計算的Ga2O3的體模量為171.75 GPa，表明其具有很好的力學(xué)性質(zhì). 理想的Ga2O3屬于間接帶隙結(jié)構(gòu)，其禁帶寬度為4.73 eV，能響應(yīng)紫外光，價帶頂平坦的能帶和導(dǎo)帶底色散的能帶表明有較大的價帶空穴和較小的導(dǎo)帶電子. 含O空位的缺陷無論在貧氧還是在富氧的條件下，形成能都小于含Ga缺陷的Ga2O3，意味著Ga2O3中空位缺陷以O(shè)空位為主，Ga2O3表現(xiàn)出n型半導(dǎo)體特點. 對于含Ga缺陷的Ga2O3，空位的出現(xiàn)并沒有對幾何結(jié)構(gòu)有顯著的影響，O空位的缺陷在禁帶中引入了缺陷能級，電子躍遷能也降低為1.93 eV. 對于含O空位的Ga2O3體系，在空位周圍的原子結(jié)構(gòu)有一定的變形，鍵長的最大改變值為0.29 ?，晶格參數(shù)的改變很小，在禁帶中出現(xiàn)雜質(zhì)能級，電子的躍遷能為2.92 eV，相對于理想的Ga2O3有顯著減小，光響應(yīng)拓展到可見光區(qū)，因此Ga2O3由于其價帶和導(dǎo)帶的顯著不同，在光催化材料領(lǐng)域有潛在的研究價值. 由于其只能吸收紫外光的特點，未來關(guān)于Ga2O3光催化的研究應(yīng)該是在保持其高效的紫外光催化活性的同時拓展光吸收范圍.