李凡生，余小英，房 慧, ，黃燦勝，王如志

(1. 廣西民族師范學(xué)院 物理與電子工程學(xué)院，崇左 532200； 2.北京工業(yè)大學(xué) 材料科學(xué)與工程學(xué)院，北京 100124)

1 引 言

作為一種n型半導(dǎo)體材料，氧化鋅的激子束縛能和帶隙較大，可見光透過率高，有優(yōu)異的常溫發(fā)光性能，在液晶顯示板、發(fā)光二極管和薄膜晶體管等器件中均有應(yīng)用. 氧化鋅具有巖鹽、類石墨和六方三種結(jié)構(gòu)，其中六方纖鋅礦結(jié)構(gòu)氧化鋅最為穩(wěn)定，在光電、催化、壓電和熱電等領(lǐng)域具有大量的應(yīng)用[1-5]. 本征氧化鋅體材料的載流子濃度約為1.7×1017，但是其載流子遷移率較高[3]. 文獻(xiàn)報(bào)道中通過摻雜和結(jié)構(gòu)優(yōu)化來改善其體材料的載流子性質(zhì)，優(yōu)化其光電性能[6-8]. 此外，低維化的氧化鋅也可作為一種功能材料在相關(guān)領(lǐng)域發(fā)揮作用. 李等人制備了氧化鋅基薄膜并研究了其光電性能. 結(jié)果表明主族元素?fù)诫s氧化鋅薄膜的可見光透過率保持在 80%以上，其電性能也受制備條件影響，比如升高樣品退火溫度，載流子濃度和遷移率均增大，導(dǎo)電性能也增強(qiáng)[8]. 氧化鋅塊體材料和低維材料在制備過程中，經(jīng)常會(huì)出現(xiàn)缺陷，如Zn和O的空位缺陷，Zn和O的位置互換，以及Zn和O的間隙缺陷等，然而在低維氧化鋅材料尤其是二維膜材料方面的報(bào)道相對(duì)較少，對(duì)其高原子密度晶面的電子學(xué)性質(zhì)方面未見報(bào)道，在氧空位缺陷對(duì)低維材料電子結(jié)構(gòu)影響方面也未見報(bào)道. 本文在密度泛函理論計(jì)算分析方法的基礎(chǔ)上，系統(tǒng)研究了本征氧化鋅(110)二維膜材料的電子結(jié)構(gòu)性質(zhì)，并研究了氧空位對(duì)氧化鋅(110)二維膜材料電子結(jié)構(gòu)性質(zhì)的影響.

2 計(jì)算模型與過程

纖鋅礦ZnO塊體材料呈六角對(duì)稱，空間群為P63mc，晶格參數(shù)a=b=3.249 ?，c=5.229 ?，α=β=90°，γ=120°[3]. 為了得到ZnO基(110)二維膜材料的結(jié)構(gòu)模型，對(duì)其體材料超晶胞進(jìn)行(110)晶面的切割，進(jìn)而得到ZnO基(110)二維膜材料結(jié)構(gòu)，分子式為Zn16O16，并在(110)二維膜兩側(cè)設(shè)置厚度均1nm的真空層，以隔斷庫侖相互作用，在此基礎(chǔ)上進(jìn)行相關(guān)計(jì)算分析. 圖1所示為ZnO基(110)二維膜結(jié)構(gòu)示意圖，圖中標(biāo)出了Zn原子和O原子的位置. 自從Kohn等人建立密度泛函理論之后，固體材料電子性質(zhì)的量子化學(xué)計(jì)算變得容易. 本論文工作建立在密度泛函理論框架之內(nèi)，周期性勢(shì)場(chǎng)作用下的電子結(jié)構(gòu)計(jì)算分析采用贗勢(shì)法進(jìn)行，并采用Vanderbilt 的超軟勢(shì)，將Zn和O的價(jià)電子取為Zn(3d104s2)和O(2s22p4)，所有電子的波函數(shù)采用平面波表示. 過程中將內(nèi)層d態(tài)電子的強(qiáng)關(guān)聯(lián)相互作用考慮進(jìn)去，并將此能量設(shè)置為1.5 eV. 電子之間的交換關(guān)聯(lián)作用采用廣義梯度法近似(PBE，Perdew Burke Ernzerhof)分別展開[9-11]. 計(jì)算過程中收斂精度設(shè)置為，平面波矢組基矢截?cái)嗄?00 eV，能量收斂精度0.002 meV/atom，平面波電子波函數(shù)的傅里葉變換網(wǎng)格密度設(shè)置為48×90×54，采用Monkhorst-pack函數(shù)進(jìn)行布里淵區(qū)k點(diǎn)的平均采樣，其k點(diǎn)網(wǎng)格為1×1×1.

圖 1 ZnO基(110)二維膜結(jié)構(gòu)模型和晶面示意圖Fig. 1 Schematic pictures of 3-dimension and the plane for two dimensional (110) thin film of ZnO

3 結(jié)果與討論

化合物的生成焓是指由元素單質(zhì)反應(yīng)生成該化合物所釋放的能量，其值的大小表明這種化合物的熱力學(xué)穩(wěn)定性，也表明這種化合物的制備難易程度. 一般來說形成焓越小，該化合物就越容易形成，其熱力學(xué)穩(wěn)定性就越高. 以化學(xué)式為AmBn的化合物生成過程為例，其生成焓可以通過如下公式計(jì)算：

(1)

上式中，ΔH為生成炩，EAB是反應(yīng)生成的化合物晶胞系統(tǒng)的總能，EA和EB分別為A元素單質(zhì)和B元素單質(zhì)的平均能量，mA和nB分別是化合物晶胞中A元素和B元素的數(shù)量. 由表1可知，本征ZnO基(110)二維膜的形成焓為-4.18 eV，氧空位的ZnO基(110)二維膜的形成焓為-4.02 eV，結(jié)果表明本征二維膜比氧空位的ZnO基(110) 二維膜穩(wěn)定性高.

空位形成能是指體系形成該氧空位所需要的能量，可以表征形成空位的難易程度，其值越小，形成該空位缺陷相對(duì)越容易[10-11]. 以ZnO中形成氧空位為例，其空位形成能可以通過如下公式計(jì)算：

Ef=EV-EI+μO-μV

(2)

上式中，Ef為氧空位形成能，EV為氧空位膜體系的總能，EI為本征ZnO膜體系的總能，μO和μV分別為氧和空位的化學(xué)勢(shì). 由表1可知，ZnO基(110)二維膜中氧空位的形成能為-1.79 eV，此處的負(fù)值表明ZnO基(110)二維膜有失去氧的傾向.

表1 ZnO基(110)二維膜的生成焓、空位形成能、帶隙和電荷分布

Table 1 Formation enthalpy, vacancy formation energy, band gap and charge population of two dimensional (110) thin film of ZnO

Formationenthalpy/eVFormationenergy/eVBandgap/eVCharge/eIntrinsicZnO(110)film-4.18/2.30Zn0.93O-0.93OvacantZnO(110)film-4.02-1.791.87Zn0.84/0.93/0.95/0.96O-0.93/-0.95/-0.96

圖2給出ZnO本征(110)與O空位(110)膜材料的能帶結(jié)構(gòu). 從圖中可以看出，ZnO本征(110)膜的能帶中導(dǎo)帶態(tài)和淺價(jià)帶態(tài)能級(jí)分布變化較為劇烈，而深價(jià)帶態(tài)能級(jí)分布變化較為平緩，表明兩個(gè)區(qū)域能級(jí)具有不同的有效質(zhì)量. 氧空位(110)膜的能帶中導(dǎo)帶態(tài)和價(jià)帶態(tài)能級(jí)分布變化均較為平緩，表明這些能級(jí)有效質(zhì)量均較大[5]. 從圖中還能看出，氧空位的(110)膜導(dǎo)帶向下移動(dòng)，并且導(dǎo)帶中的能級(jí)產(chǎn)生了簡(jiǎn)并，表明這些電子定域性的增強(qiáng).

圖 2 ZnO本征(110)與O空位(110)膜材料的能帶結(jié)構(gòu)Fig. 2 Band structures for intrinsic (110) and O vacant (110) thin film of ZnO

圖3給出ZnO本征(110)與O空位(110)膜材料的微細(xì)能帶結(jié)構(gòu). 從圖中可以看出，ZnO本征(110)膜的能帶中導(dǎo)帶底能級(jí)有3個(gè)能谷，其能量分別為2.3 eV，2.31 eV和2.842 eV，其能級(jí)具有正的有效質(zhì)量，其為多能谷材料. 價(jià)帶頂能級(jí)中有2個(gè)能峰，分別位于費(fèi)米能和-0.005 eV，其能級(jí)具有負(fù)的有效質(zhì)量. 其導(dǎo)帶底和價(jià)帶頂均位于G點(diǎn)，因此ZnO本征(110)膜為直接帶隙型材料，帶隙寬度為2.3 eV. 氧空位(110)膜材料的能帶結(jié)構(gòu)中，導(dǎo)帶底能級(jí)有2個(gè)能谷，分別位于1.877 eV和1.88 eV，可以看出這些位置的能級(jí)有效質(zhì)量比本征膜大幅度增大，這些位置的電子速度普遍較低. 其價(jià)帶頂能級(jí)中也有2個(gè)能峰，分別位于-0.027 eV和費(fèi)米能處，可以看出，G點(diǎn)處的價(jià)帶空穴有效質(zhì)量大于本征態(tài)膜. 其導(dǎo)帶底和價(jià)帶頂均位于不同點(diǎn)，因此氧空位(110)膜為間接帶隙型材料，帶隙寬度為1.877 eV. 對(duì)比它們的帶隙寬度可以看出，氧空位使(110)膜的帶間躍遷勢(shì)壘降低，導(dǎo)電性增強(qiáng).

圖 3 ZnO本征(110)與O空位(110)膜材料的微細(xì)能帶結(jié)構(gòu)Fig. 3 Detailed band structures for intrinsic (110) and O vacant (110) thin film of ZnO

圖4給出了ZnO本征(110)與O空位(110)膜材料的態(tài)密度. 從圖中可以看出，ZnO本征(110)膜和O空位(110)膜材料態(tài)密度曲線在價(jià)帶中均有3個(gè)峰值，表明這些位置電子定域性較強(qiáng)，與能帶結(jié)構(gòu)結(jié)果一致.圖5給出了ZnO本征(110)與O空位(110)膜材料的微細(xì)態(tài)密度曲線. 從圖中看出s態(tài)電子對(duì)ZnO本征(110)膜費(fèi)米能處態(tài)密度貢獻(xiàn)小于p態(tài)和d態(tài)電子，而在導(dǎo)帶態(tài)密度中s態(tài)電子對(duì)態(tài)密度貢獻(xiàn)大于p態(tài)和d態(tài)電子. 氧空位(110)膜費(fèi)米能處p態(tài)電子對(duì)態(tài)密度貢獻(xiàn)大于s態(tài)和d態(tài)電子，而在導(dǎo)帶態(tài)密度中s態(tài)電子對(duì)態(tài)密度貢獻(xiàn)大于p態(tài)和d態(tài)電子. ZnO的(110)膜經(jīng)過氧空位之后，s電子參與形成價(jià)帶頂，其能級(jí)可形成空穴. 圖6分別給出了ZnO本征(110)與O空位(110)膜材料中Zn和O的分波態(tài)密度曲線. 可以看出(110)膜Zn中的d電子在費(fèi)米能上貢獻(xiàn)最大，s電子貢獻(xiàn)為零. 而氧空位(110)膜的Zn中的s電子在費(fèi)米能上貢獻(xiàn)最大，d電子貢獻(xiàn)大幅度降低，聯(lián)系圖5的分析可知，(110)膜產(chǎn)生氧空位之后，Zn的s電子參與形成其價(jià)帶頂能級(jí). 氧原子中的s和d電子對(duì)費(fèi)米能上態(tài)密度的貢獻(xiàn)趨勢(shì)不變.

圖 4 ZnO本征(110)與O空位(110)膜材料的態(tài)密度Fig. 4 Densities of states for intrinsic (110) and O vacant (110) thin film of ZnO

圖 5 ZnO本征(110)與O空位(110)膜材料的微細(xì)態(tài)密度Fig. 5 Detailed densities of states for intrinsic (110) and O vacant (110) thin film of ZnO

圖7給出了ZnO本征(110)膜材料中s、p和d電子的分波態(tài)密度曲線. 可以看出所有電子均在-1 eV附近具有一個(gè)峰位，表明電子在-1 eV附近定域化. ZnO本征(110)膜材料中s、p和d電子在費(fèi)米能上分波態(tài)密度取值分別為0，2.34和0.14，其中p電子貢獻(xiàn)最大，其參與形成價(jià)帶頂能級(jí). 圖8給出了氧空位ZnO(110)膜材料中s、p和d電子的分波態(tài)密度曲線. 可以看出所有電子均在-1.5 eV和費(fèi)米能附近具有峰位，表明電子在這些能量附近定域化. 氧空位(110)膜材料中s、p和d電子在費(fèi)米能上分波態(tài)密度取值分別為1.59，1.11和0.19，其中s電子貢獻(xiàn)最大，其參與形成價(jià)帶頂能級(jí)，這與上述圖5和圖6的分析結(jié)果一致.

圖9分別給出了ZnO本征(110)與O空位(110)膜材料的電子密度圖. 從圖中可以看出，本征(110)膜中Zn-O鍵長(zhǎng)有兩種，距離較近的Zn-O鍵電子云共用明顯，距離較遠(yuǎn)的Zn-O鍵庫侖作用明顯. O空位(110)膜中的Zn-O鍵的電子 云情況有模糊化傾向. 結(jié)合Zn-O鍵的分布數(shù)據(jù)可以看出，ZnO本征(110)膜材料中其長(zhǎng)度分別為0.17641 nm和0.18258 nm，其數(shù)量分別為0.58和0.4，距離較近的Zn-O鍵傾向于共價(jià)型結(jié)合，距離較遠(yuǎn)的Zn-O鍵傾向于離子型結(jié)合. O空位(110)膜材料中Zn-O鍵數(shù)量分別為0.49，0.51和0.52以及0.45，0.46，0.48，0.49和0.5，表現(xiàn)為混合型結(jié)合傾向的增大.

圖 6 ZnO本征(110)與O空位(110)膜材料中Zn和O的分波態(tài)密度Fig. 6 Partial densities of states of Zn and O for intrinsic (110) and O vacant (110) thin film of ZnO

圖 7 ZnO本征(110)膜材料中s、p和d電子的分波態(tài)密度

圖 8 O空位(110)膜材料中s、p和d電子的分波態(tài)密度

4 結(jié) 論

圖 9 ZnO本征(110)與O空位(110)膜材料的電子密度Fig. 9 Electron densities for intrinsic (110) and O vacant (110) thin film of ZnO

本文利用密度泛函理論方法系統(tǒng)研究了本征氧化鋅(110)二維膜材料與氧空位氧化鋅(110)二維膜材料的形成和電子結(jié)構(gòu)性質(zhì). ZnO的本征(110)二維膜的形成焓為-4.18 eV，氧空位的(110)二維膜的形成焓為-4.02 eV，本征二維膜穩(wěn)定性較高，ZnO的(110)二維膜有失去氧的傾向. 氧空位的(110)膜導(dǎo)帶向下移動(dòng)，并且導(dǎo)帶中的能級(jí)產(chǎn)生了簡(jiǎn)并. 本征ZnO的(110)膜為直接帶隙型多能谷材料，帶隙寬度為2.3 eV，氧空位(110)膜為間接帶隙型多能谷材料，帶隙寬度為1.877 eV. 氧空位(110)膜材料的能帶結(jié)構(gòu)中，導(dǎo)帶底能級(jí)有2個(gè)能谷，分別位于1.877 eV和1.88 eV，這些位置的能級(jí)有效質(zhì)量比本征膜大幅度增大，電子速度普遍較低. ZnO的(110)膜產(chǎn)生氧空位之后，Zn的s電子參與形成其價(jià)帶頂能級(jí). 氧空位(110)膜材料中Zn-O鍵的混合型結(jié)合傾向增大.