底琳佳， 戴顯英， 苗東銘， 吳淑靜， 郝 躍

(1. 西安電子科技大學(xué) 微電子學(xué)院，陜西 西安 710071；2. 西安電子科技大學(xué) 寬帶隙半導(dǎo)體國家重點(diǎn)學(xué)科實(shí)驗(yàn)室，陜西 西安 710071)

單軸應(yīng)變鍺(Ge)具有載流子遷移率提升高、與硅(Si)集成工藝兼容等優(yōu)勢，在Si基高遷移率器件與集成電路應(yīng)用方面引起了人們的廣泛關(guān)注[1-3]．能帶結(jié)構(gòu)是進(jìn)一步研究應(yīng)變Ge電學(xué)、光學(xué)性能的理論基礎(chǔ)，而k·p方法不僅具有高效的計(jì)算效率，且能夠準(zhǔn)確地描繪出在布里淵區(qū)給定點(diǎn)附近的導(dǎo)帶和價帶結(jié)構(gòu)，通常用于計(jì)算載流子的有效質(zhì)量和運(yùn)輸特性．

目前，針對應(yīng)變Ge能帶結(jié)構(gòu)和運(yùn)輸特性的研究主要集中在價帶和空穴遷移率增強(qiáng)方面[4-5]，關(guān)于導(dǎo)帶結(jié)構(gòu)，特別是與電子散射和運(yùn)輸密切相關(guān)的有效質(zhì)量的理論研究相對較少．此外，已有研究表明雙軸應(yīng)變Ge可以從間接帶隙轉(zhuǎn)變?yōu)橹苯訋栋雽?dǎo)體，并具有與Ⅲ-Ⅴ化合物半導(dǎo)體相比擬的光增益，可用于光電子器件的設(shè)計(jì)[6-7]．與雙軸應(yīng)變相比，單軸應(yīng)變具有一些顯著的優(yōu)點(diǎn)，例如，在高垂直電場下維持空穴遷移率增益，為納米級短溝道器件提供大的驅(qū)動電流提升等[8]．因此，有關(guān)單軸應(yīng)變Ge是否存在帶隙轉(zhuǎn)換以及單軸應(yīng)力如何影響電子輸運(yùn)的研究具有重要的意義．

筆者采用30k·p方法，計(jì)算了不同單軸應(yīng)力(方向、類型和量級)作用下Ge的能帶結(jié)構(gòu)，定量分析了Ge導(dǎo)帶能谷分裂與偏移、縱向和橫向電子有效質(zhì)量及電子態(tài)密度有效質(zhì)量隨應(yīng)力的變化情況．

1 30 k·p方法

應(yīng)變系統(tǒng)的薛定諤方程可由下式表征[9]:

(1)

其中，p=-i(h/(2π))，為動量算符，m為自由電子質(zhì)量，V0(r)為未形變晶格的周期勢，i、j代表x、y、z，ε為應(yīng)變張量．將布洛赫波函數(shù)φnk(r)= exp(ik·r)unk(r) 代入式(1)，可得

式(2)中的哈密頓量即為包含應(yīng)變效應(yīng)的k·p哈密頓量．

如果未形變系統(tǒng)哈密頓量的特征值En0及特征函數(shù)un0已知，如下式所示:

(5)

可將式(3)中的k·p乘積項(xiàng)及式(4)看作微擾量，并以已知的特征函數(shù)為基展開unk，根據(jù)微擾理論可以得到一系列k·p久期方程:

式(6)中的前兩項(xiàng)用于建立未形變系統(tǒng)的k·p哈密頓矩陣，第3項(xiàng)用于建立應(yīng)變引起的微擾矩陣．

對于Ge，選取來自[000]、(2π/a)[111]和(2π/a)[200]平面波的15個Γ態(tài)[10]，并根據(jù)群論選擇規(guī)則可將由式(6)中第2項(xiàng)決定的動量矩陣參數(shù)減小到10個:

圖1 30 k·p模型原理圖

此外，由于自旋軌道相互作用的存在，根據(jù)15個Γ態(tài)所屬原子態(tài)的不同，可將其分為3類[11]:Γ1l、Γ2′l、Γ1u和Γ2′u屬于s原子態(tài)，自旋軌道相互作用使它們變?yōu)槎群啿B(tài);Γ12′屬于d原子態(tài)，自旋軌道相互作用，使其變?yōu)樗亩群啿B(tài);Γ25′l、Γ25′u和Γ15屬于p原子態(tài)，自旋軌道相互作用使它們分裂為一組四度簡并態(tài)和一組二度簡并態(tài)，相應(yīng)的分裂能分別為Δ25′l、Δ25′u和Δ15，且在Γ25′l和Γ25′u之間存在非零的耦合系數(shù)Δ25u，25l．圖1給出了30k·p模型的原理圖，并列出了相應(yīng)的Γ態(tài)及其對應(yīng)的原子態(tài)波函數(shù)、動量矩陣參數(shù)、自旋軌道相互作用系數(shù)．

因此，根據(jù)以上Γ態(tài)可以建立包含應(yīng)變效應(yīng)的30k·p哈密頓矩陣，因?yàn)榫仃囆问较鄬?fù)雜，文中并未列出，相應(yīng)的矩陣可在文獻(xiàn)[12]中找到．根據(jù)已有的單軸應(yīng)變張量模型[13]，通過對角化 30k·p哈密頓矩陣可得到單軸應(yīng)變Ge的能帶結(jié)構(gòu)．計(jì)算過程中各參數(shù)的具體數(shù)值列在表1和表2中，由式(7)決定的應(yīng)變微擾矩陣系數(shù)和文獻(xiàn)[12]中的一致．

表1 鍺在Γ點(diǎn)處的能量本征值及自旋軌道相互作用系數(shù)

表2 鍺的動量矩陣參數(shù)

圖2 弛豫Ge(實(shí)線)和[111]方向2 GPa單軸應(yīng)變Ge(虛線)能帶結(jié)構(gòu)

2 單軸應(yīng)變Ge帶隙特性

圖2給出了弛豫Ge(實(shí)線)和沿[111]方向的2 GPa單軸應(yīng)變Ge(虛線)的能帶圖．對于弛豫Ge，計(jì)算得到的L和Δ能谷能量分別為 0.662 eV 和 0.851 eV．從圖2可以看出，單軸應(yīng)變Ge的各能級均發(fā)生了明顯的偏移和翹曲，說明帶隙寬度和有效質(zhì)量會隨著應(yīng)力而變化．此外，由于應(yīng)力的引入，晶格對稱性遭到破壞，原有的L和Δ能谷簡并度將發(fā)生改變．

圖3 單軸應(yīng)變Ge導(dǎo)帶各能谷隨應(yīng)力的變化情況

值得注意的是，在單軸應(yīng)力的作用下，Ge導(dǎo)帶帶邊能級會發(fā)生改變．圖4給出了不同方向單軸應(yīng)力作用下Ge導(dǎo)帶各能谷帶隙寬度隨應(yīng)力的變化情況，可以更直觀反應(yīng)帶邊能級的轉(zhuǎn)換．由圖4(a)可見，當(dāng)沿[001]方向的單軸壓應(yīng)力大于 2.1 GPa 后，導(dǎo)帶帶邊能級變?yōu)棣つ芄?，但Ge仍為間接帶隙半導(dǎo)體;而在張應(yīng)力作用下，由于Γ能谷帶隙寬度隨應(yīng)力增大而減小的趨勢遠(yuǎn)大于L能谷，使得在應(yīng)力增大到一定程度后，Γ能谷將成為導(dǎo)帶帶邊能級，根據(jù)計(jì)算結(jié)果，當(dāng)張應(yīng)力大于 4.8 GPa 后，Ge從原來的間接帶隙半導(dǎo)體變?yōu)橹苯訋栋雽?dǎo)體，禁帶寬度為 0.402 eV．類似地，由圖4(c)可見，沿[111]方向的單軸張應(yīng)力大于 6.2 GPa 后，Ge變?yōu)橹苯訋栋雽?dǎo)體，禁帶寬度為 0.244 eV．

圖4 單軸應(yīng)變Ge導(dǎo)帶各能谷帶隙寬度隨應(yīng)力的變化情況

文中結(jié)果與文獻(xiàn)[13]結(jié)果不同，這是由于文獻(xiàn)[13]未考慮應(yīng)力作用下高能量能級對導(dǎo)帶能谷的影響，且其形變勢參數(shù)并不充分[6]．此外，文獻(xiàn)[14]采用第一性原理對張應(yīng)變Ge能帶結(jié)構(gòu)進(jìn)行了計(jì)算，發(fā)現(xiàn)[001]和[111]方向單軸張應(yīng)力均可使Ge由間接帶隙轉(zhuǎn)變成直接帶隙．文獻(xiàn)[15]報道的[001]方向懸浮Ge納米線結(jié)構(gòu)在產(chǎn)生的5.7% (5.87 GPa) 的單軸張應(yīng)變作用下，已轉(zhuǎn)變?yōu)橹苯訋栋雽?dǎo)體．這些相關(guān)結(jié)論可以間接驗(yàn)證文中所得結(jié)果的準(zhǔn)確性．

3 單軸應(yīng)變Ge電子有效質(zhì)量

電子有效質(zhì)量取決于其所處狀態(tài)下的能帶結(jié)構(gòu)，可由下式表征:

(8)

根據(jù)單軸應(yīng)變Ge能帶結(jié)構(gòu)，采用二次函數(shù)擬合方法可以得到Ge導(dǎo)帶各能谷電子有效質(zhì)量．為了驗(yàn)證計(jì)算的有效性，表3列出了實(shí)驗(yàn)測量[16]、贗勢法[17]、緊束縛法[18]以及文中 30k·p方法計(jì)算得到的弛豫Ge導(dǎo)帶能谷電子有效質(zhì)量．表中，ml表示縱向有效質(zhì)量，mt表示橫向有效質(zhì)量，可以看出，所得結(jié)果基本一致．

表3 弛豫Ge導(dǎo)帶各能谷電子有效質(zhì)量

圖5分別給出了不同方向單軸應(yīng)力作用下Ge導(dǎo)帶L和Δ能谷縱向、橫向有效質(zhì)量隨應(yīng)力的變化情況．單軸應(yīng)力破壞了晶體對稱性，使導(dǎo)帶能谷產(chǎn)生分裂和偏移，有效質(zhì)量也會隨之改變．總體而言，導(dǎo)帶L能谷在各方向單軸應(yīng)力作用下，縱向、橫向有效質(zhì)量變化并不明顯，而Δ能谷在[110]、[111]單軸應(yīng)力作用下有大幅改變，這意味著剪切應(yīng)力對于Δ能谷的作用比L能谷更為顯著．

圖5 單軸應(yīng)變Ge導(dǎo)帶L、Δ能谷縱向和橫向有效質(zhì)量隨應(yīng)力的變化情況

態(tài)密度有效質(zhì)量與有效狀態(tài)密度、載流子散射幾率密切相關(guān)，也是應(yīng)變Ge的重要參數(shù)．應(yīng)變Ge電子態(tài)密度有效質(zhì)量可由下式給出[19]:

(9)

其中，mc=(mlmt2)1/3為導(dǎo)帶某一能谷的態(tài)密度有效質(zhì)量，ΔEc，split為導(dǎo)帶能谷分裂能，kB為玻爾茲曼常數(shù)，a和b分別表示能級分裂后產(chǎn)生的能量較低能谷的簡并度和能量較高能谷的簡并度，例如[111]單軸張應(yīng)力作用下，對于Δ能谷，a、b分別為6、0，對于L能谷，a、b分別為3、1．

圖6 單軸應(yīng)變Ge導(dǎo)帶電子態(tài)密度有效質(zhì)量

圖6為不同方向單軸應(yīng)力作用下Ge導(dǎo)帶L和Δ能谷的電子態(tài)密度有效質(zhì)量隨應(yīng)力的變化情況．由圖6(a)可見，沿[110]和[111]方向的單軸應(yīng)力可以有效降低L能谷電子態(tài)密度有效質(zhì)量，并在應(yīng)力大于 1 GPa 后，趨近為常數(shù)，其中沿[111]方向的單軸壓應(yīng)力可使L能谷電子態(tài)密度有效質(zhì)量取得最小值 0.22m0．類似地，由圖6(b)可見，Δ能谷的電子態(tài)密度有效質(zhì)量在沿[001]方向單軸壓應(yīng)力大于 1 GPa 后取得最小值 0.52m0．

態(tài)密度有效質(zhì)量的顯著降低，是由于單軸應(yīng)力使導(dǎo)帶能谷發(fā)生分裂，能級簡并度降低，隨著分裂能的增大，電子越來越集中于低能量能谷．態(tài)密度有效質(zhì)量的減小，可以有效降低載流子的散射幾率，提高遷移率，從這一角度考慮，為了得到更高的電子遷移率，沿[111]方向的單軸壓應(yīng)變Ge可作為N溝道金屬氧化物半導(dǎo)體(N-channel Metal Oxide Semiconductor，NMOS)器件導(dǎo)電溝道材料的最佳選擇．

4 結(jié) 束 語

文中采用30k·p方法，對沿[001]、[110]、[111]方向的單軸應(yīng)變Ge能帶結(jié)構(gòu)進(jìn)行了計(jì)算．單軸應(yīng)力的引入使得各能級發(fā)生不同程度的翹曲與偏移，進(jìn)而改變了單軸應(yīng)變Ge的帶隙特性和電子有效質(zhì)量．結(jié)果表明:在[001]、[111]方向單軸張應(yīng)力作用下，Ge由間接帶隙轉(zhuǎn)變?yōu)橹苯訋?，帶隙性質(zhì)的轉(zhuǎn)變使得單軸應(yīng)變Ge可用于光電集成器件; 導(dǎo)帶L和Δ能谷電子縱向、橫向有效質(zhì)量并不明顯地依賴于單軸應(yīng)力，但沿[111]和[001]方向的單軸壓應(yīng)力可分別使L和Δ能谷態(tài)密度有效質(zhì)量取得最小值，這有利于減小電子散射幾率，提升電子遷移率．所得相關(guān)結(jié)論可為設(shè)計(jì)高性能單軸應(yīng)變鍺器件應(yīng)力和溝道方向選擇提供理論參考．

[1] LIU Y, NIU J B, WANG H J, et al. Strained Germanium Quantum Well PMOSFETs on SOI with Mobility Enhancement by External Uniaxial Stress[J]. Nanoscale Research Letters, 2017, 12(1): 120-124.

[2] LEE C F, HE R Y, CHEN K T, et al. Strain Engineering for Electron Mobility Enhancement of Strained Ge NMOSFET with SiGe Alloy Source/Drain Stressors[J]. Microelectronic Engineering, 2015, 138(C): 12-16.

[3] LIU J S, CLAVEL M B, HUDAIT M K. An Energy-efficient Tensile-strained Ge/InGaAs TFET 7T SRAM Cell Architecture for Ultralow-voltage Applications[J]. IEEE Transactions on Electron Devices, 2017, 64(5): 2193-2200.

[4] MA J L, FU Z F, LIU P, et al. Hole Mobility Enhancement in Uniaxial Stressed Ge Dependence on Stress and Transport Direction[J]. Science China: Physics, Mechanics and Astronomy, 2014, 57(10): 1860-1865.

[5] 白敏, 宣榮喜, 宋建軍, 等. 壓應(yīng)變Ge/(001)Si1-xGex空穴散射與遷移率模型[J]. 物理學(xué)報, 2015, 64(3): 038501-038506.

BAI Min, XUAN Rongxi, SONG Jianjun, et al. Hole Scattering and Mobility in Compressively Strained Ge/(001)Si1-xGex[J]. Acta Physica Sinica, 2015, 64(3): 038501-038506.

[6] KURDI M E, FISHMAN G, SAUVAGE S, et al. Band Structure and Optical Gain of Tensile-strained Germanium Based on a 30 Bandk·pFormalism[J]. Journal of Applied Physics, 2010, 107(1): 013710-013717.

[7] TANI K, SAITO S, ODA K, et al. Room-temperature Direct Band-gap Electroluminescence from Germanium(111)-Fin Light-emitting Diodes[J]. Japanese Journal of Applied Physics, 2017, 56(3): 032102.

[8] DAI X Y, SHAO C F, HAO Y. Uniaxially Strained Silicon on Insulator with Wafer Level Prepared by Mechanical Bending and Annealing[J]. Applied Physics Express, 2013, 6(8): 081302.

[9] SUN Y, THOMPSON S E, NISHIDA T. Strain Effect in Semiconductors:Theory and Device Applications[M]. New York: Springer, 2010: 118.

[10] CARDONA M, PQLLAK F H. Energy-band Structure of Germanium and Silicon: thek·pMethod[J]. Physical Review, 1966, 142(2): 530-543.

[11] RICHARD S, ANIEL F, FISHMAN G. Energy-band Structure of Ge, Si, and GaAs: a Thirty-bandk·pMethod[J]. Physical Review B, 2004, 70(23): 18976-18984.

[12] RIDEAU D, FERAILLE M, CIAMPOLINI L, et al. Strained Si, Ge, and Si1-xGexAlloys Modeled with a First-principles-optimized Full-zonek·pMethod[J]. Physical Review B, 2006, 74(19): 195208-195227.

[13] 戴顯英, 李金龍, 郝躍. 應(yīng)變鍺的導(dǎo)帶結(jié)構(gòu)計(jì)算與分析[J]. 西安電子科技大學(xué)學(xué)報, 2014, 41(2): 120-124.

DAI Xianying, LI Jinlong, HAO Yue. Calculation and Analysis of the Conduction Band-structure of Strained Germanium[J]. Journal of Xidian University, 2014, 41(2): 120-124.

[14] LIU L, ZHANG M, HU L, et al. Effect of Tensile Strain on the Electronic Structure of Ge: a First-principles Calculation[J]. Journal of Applied Physics, 2014, 116(11): 113105.

[15] SUKHDEO D S, NAM D, KANG J H, et al. Direct Bandgap Germanium-on-Silicon Inferred from 5.7%100Uniaxial Tensile Strain[J]. Photonics Research, 2014, 2(3): A8-A13.

[16] MADELUNG O, R?SSLER U, SCHULZ M. Semiconductors·Group Ⅳ Elements, Ⅳ-Ⅳ and Ⅲ-Ⅴ Compounds. Part b - Electronic, Transport, Optical and Other Properties[M]. Berlin: Springer, 2002: 2854.

[17] KIM J, FISCHETTI M V. Electronic Band Structure Calculations for Biaxially Strained Si, Ge, and Ⅲ-Ⅴ Semiconductors[J]. Journal of Applied Physics, 2010, 108(1): 013710.

[18] NIQUET Y M, RIDEAU D, TAVERNIER C, et al. Onsite Matrix Elements of the Tight-binding Hamiltonian of a Strained Crystal: Application to Silicon, Germanium, and Their Alloys[J]. Physical Review B, 2009, 79(24): 245201-245213.

[19] SONG J J, CHAO Y, ZHANG H M, et al. Longitudinal, Transverse, Density-of-states, and Conductivity Masses of Electrons in (001), (101) and (111) Biaxially-strained-Si and Strained-Si1-xGex[J]. Science China: Physics, Mechanics and Astronomy, 2012, 55(11): 2033-2037.