應(yīng)變與ZnSe/ZnCdSe雙勢壘對電子隧穿的壓力影響研究

錢岙軻

摘 要： 運用求解任意勢中波函數(shù)與轉(zhuǎn)移矩陣相結(jié)合的方法，計算得到流體靜壓下應(yīng)變ZnCdSe/ZnSe雙勢壘電子隧穿的共振能級、波函數(shù)和透射系數(shù)?？紤]晶格常數(shù)、有效質(zhì)量及體彈性模量等參量的壓力效應(yīng)對電子隧穿的影響。數(shù)值結(jié)果顯示，與無靜壓比較，勢壘的高度增加，則電子有效質(zhì)量減小，導(dǎo)致共振峰向高能區(qū)移動，但峰值高度不變。此外，給出勢阱有無靜壓的兩種情形下的波形圖，通過對比，可以進(jìn)一步看出靜壓對ZnCdSe/ZnSe電子隧穿的影響。

關(guān)鍵詞： 電子隧穿;靜壓;雙勢壘;透射系數(shù);電子波函數(shù)

中圖分類號： TB????? 文獻(xiàn)標(biāo)識碼： A????? doi：10.19311/j.cnki.1672-3198.2019.12.093

1 引言

近年來，由Ⅱ-Ⅵ族寬禁帶半導(dǎo)體材料組成的量子異質(zhì)結(jié)構(gòu)具有很大的潛在應(yīng)用價值，如高亮度藍(lán)/綠光發(fā)光二極管和激光二極管，而且利用共振隧穿可大幅度提高光子器件和電子器件的工作速度和效率，這使得異質(zhì)結(jié)構(gòu)中量子隧穿問題成為物理學(xué)界關(guān)注的熱點之一.研究共振隧穿現(xiàn)象的基本問題之一，是計算電子穿越勢壘的透射系數(shù)。

Messiah曾于1961年給出了一種計算電子穿越任意勢壘之透射系數(shù)的WKB方法，但此方法用于突變勢情形，結(jié)果不夠理想，且不能解釋Chang，Esaki和Tsu關(guān)于半導(dǎo)體異質(zhì)結(jié)構(gòu)中電子的共振隧穿。隨后，Chandra等人對于平面摻雜三角勢，通過解薛定諤方程，得到了透射系數(shù)的數(shù)值結(jié)果。之后，Ando等人采用多階式近似和改進(jìn)的轉(zhuǎn)移矩陣法討論了電子穿越任意勢的共振隧穿性質(zhì)，但此方法不易求得波函數(shù)。班士良等人曾將求解任意勢中波函數(shù)與轉(zhuǎn)移矩陣相結(jié)合，給出計算電子穿越任意勢的透射系數(shù)、反射系數(shù)及波函數(shù)之較為簡單的數(shù)值計算方法。

基于上述的研究，考慮靜壓和晶格失配，運用求解任意勢中波函數(shù)與轉(zhuǎn)移矩陣相結(jié)合的方法，計算出了流體靜壓下應(yīng)變ZnCdSe/ZnSe雙勢壘電子隧穿的波函數(shù)和透射系數(shù)，并通過仿真模擬，得到了靜壓對電子隧穿的影響。

2 理論模型與計算

考慮電子在雙勢壘一維勢的x方向的運動，將雙壘異質(zhì)結(jié)分成三個區(qū)域，在1區(qū)（xb）中，勢場為常數(shù)，記為V01和V03; 在2區(qū)（a≤x≤b）中，勢場V（x）是任意有限雙壘勢。電子波函數(shù)SymbolyA@i（x）（i=1，2，3）可以寫為：

ψ1 x =A1eik1x+B1eik1x （1）

ψ2 x =A2f x +B2g x? （2）

ψ3 x =A3eik3x+B3e-ik3x （3）

考慮到波函數(shù)ψi（x）（i=1，2，3）及其一階導(dǎo)數(shù)在邊界處的連續(xù)性，根據(jù)反射系數(shù)R和透射系數(shù)T的定義

R=（t21/t11）2 T=（1/t11）2（k3/k1）

又因為R+T=1，波函數(shù)的系數(shù)滿足以下關(guān)系：

A1B1? =

eik1a e-ik1a

ik1eik1a -ik1e-ik2a? -1

f a? g a

m1 a? m2 a? f′ a?? m1 a? m2 a? g′ a???? A2B2? =

a11 a12a21 a22

A2B2??? （4）

因為區(qū)域3中沒有反射波的存在，所以B3=0，則各區(qū)域中電子的非歸一化波函數(shù)可寫為：

Ψ1 A1 =eik1x+ t21 t11 e-ik1x （5）

Ψ2 A1 =?? 1 a11a22-a21a12?? a22-a12 t21 t11? f x +? t21 t11 a11-a21 g x????? （6）

Ψ3 A1 = 1 t11 eik3x （7）

將區(qū)域2分成j個子區(qū)域，則會有j+1個節(jié)點，每個子區(qū)域的寬度h=（b-a）/j，公式2中的f（x）和g（x）均記為y，則滿足：

- h-2 2? d dx? 1 m2 x?? dy dx +V x y=- h-2 2? dz dx +V x y=Ey （8）

由上式可以得到：

z′= dz dx =- 2 h-2? E-V x? y=F x，y x??? （9）

dy dx =m x z x =G x，z x??? （10）

利用Rung-Kutts（RK）數(shù)值解法求得波函數(shù)f（x）和g（x）在x=a b 點的平面波初始解為

f1eik2 a a，f'1=ik2 a eik2 a a，g1=e-ik2 a a，g'1=-ik2 a eik2 a a （11）

k2 a = 2m2（a）（E-V（a）） /h-

由此可得到波函數(shù)f（x）和g（x）。

在壓力下，材料的晶格常數(shù)用默納漢公式來表示，而三元混晶ZnCdSe在ZnSe界面還有晶格失配的影響，即雙軸應(yīng)變影響。計算80k下ZnSe和ZnCdSe的禁帶寬度Eg：

Eg（T1）-Eg（T2） T1-T2 = dE dT? （12）

Eg=Eg0+ acvΔV V -δEs（P） （13）

Eg（ZnCdSe）=（1-x）Eg（ZnSe）+xEg（CdSe）-bx（1-x） （14）

式中，δEs（P）=b 1+ 2C12 C11? ε（b）xx.acv是禁帶寬度的形變勢，b是切形變勢常量。

根據(jù)公式，可得電子的勢壘高度Ve：

Ve=0.70 Eg ZnSe -Eg ZnCdSe?? （15）

根據(jù)文獻(xiàn)[8]，得到電子有效質(zhì)量：

m0 m*e Ee? =1+ Ep 3?? 2 Eg+Ee + 1 Eg+Ee+△0?? （16）

式中，Ep是Kane勢關(guān)聯(lián)矩陣元所表示的質(zhì)量，Δ0是自旋軌道分裂能量，依據(jù)上式，通過數(shù)值計算可以得到在溫度80K下ZnSe和ZnCdSe的電子有效質(zhì)量。

3 數(shù)值計算結(jié)果

通過數(shù)值模擬計算，得到了透射系數(shù)隨電子能量變化的關(guān)系圖，如圖1所示，圖中波峰表明共振隧穿的位置。

由圖1可知，沒有靜壓時，在能量從0～0.12eV的變化范圍內(nèi)，共有3個共振峰，數(shù)值結(jié)果表明，共振峰所對應(yīng)的能量值分別是E01=0.01355800eV，E02 =0.0529720eV，E03=0.1154920eV，分別對應(yīng)于電子的基態(tài)、第一激發(fā)態(tài)、第二激發(fā)態(tài)。從圖2可以看出，在加了靜壓的時候，能量從0～0.12eV的變化范圍內(nèi)，共有2個共振峰，分別對應(yīng)于能量E01=0.0154000eV，E02=0.0592000eV，由其對應(yīng)的波形圖可以看出它們分別表示電子的基態(tài)和第一激發(fā)態(tài)的能量。從圖3可以看出，我們假設(shè)靜壓下電子的有效質(zhì)量是個固定值，在能量從0～0.12eV的變化范圍內(nèi)，共有3個共振峰，數(shù)值結(jié)果表明，共振峰所對應(yīng)的能量值分別是E01=0.01355800eV，E02=0.0529720eV，E03=0.1154920eV，分別對應(yīng)于電子的基態(tài)、第一激發(fā)態(tài)和第二激發(fā)態(tài)的能量。

通過對比圖1和圖2可知，對于有靜壓存在的雙勢壘結(jié)構(gòu)，使得共振峰向高能區(qū)移動，且共振峰的尖銳程度也有所變化。在圖2中，對應(yīng)于能量 E03的共振峰已經(jīng)移到了高能區(qū)，在這里由于阱寬度的選擇而沒有顯示;第一激發(fā)態(tài)與第二激發(fā)態(tài)的能量區(qū)間范圍內(nèi)透射系數(shù)增大，說明電子在此能量范圍內(nèi)的透射幾率增加。由圖1和圖3比較可知，共振峰的個數(shù)沒有變化，且兩圖中共振峰所對應(yīng)的能量值是相等的，即共振峰沒有移動。由此可得，共振峰向高能區(qū)移動是由電子有效質(zhì)量減小引起的。我們將波函數(shù)寫為：

Ψn/A1=Ψ1/A1+Ψ2/A1+Ψ3/A1=Ψnr+iΨni

（下角標(biāo)r表示實部，i表示虛部）電子在空間范圍內(nèi)出現(xiàn)的概率寫為 Ψn/A1 2（n=0，1分別代表基態(tài)、第一激發(fā)態(tài)），分別給出波形和幾率圖：

圖4和圖5分別給出無或有靜壓時電子的基態(tài)波形圖。圖6為假設(shè)靜壓下電子有效質(zhì)量是固定值的電子基態(tài)波形圖。圖4～6中，實線1代表基態(tài)電子出現(xiàn)的幾率，虛線2、3分別對應(yīng)于電子基態(tài)波函數(shù)的實、虛部波形圖，圖7和圖8分別為無、有靜壓時電子第一激發(fā)態(tài)波形圖。圖9為假設(shè)靜壓下電子有效質(zhì)量是固定值的電子第一激發(fā)態(tài)波形圖。圖7～9中，實線1代表第一激發(fā)態(tài)電子出現(xiàn)的幾率，虛線2、3分別對應(yīng)于電子第一激發(fā)態(tài)波函數(shù)的實、虛波形圖。

由圖4～9可知，有無靜壓，波函數(shù)在-30-130A之間都呈現(xiàn)良好的對稱性，它也表明對稱性以及邊界匹配性影響波函數(shù)，波函數(shù)隨著能量而變化。 通過圖5與圖6、圖8與圖9的對比可知，波函數(shù)虛部波形圖有明顯的變化，說明電子有效質(zhì)量變化引起波函數(shù)虛部的變化。目前結(jié)果表明，對于波函數(shù)的幾率不僅實部對其有貢獻(xiàn)，而且虛部也起著非常重要的作用。這種有趣的量子力學(xué)效應(yīng)有待進(jìn)一步的討論 。應(yīng)該指出的是，本文采用的一維勢可極易推廣到三維情形。

4 結(jié)論

通過數(shù)值模擬得出以下結(jié)論：靜壓使得雙勢壘結(jié)構(gòu)的共振峰向高能區(qū)移動，共振峰的尖銳程度變小。共振峰向高能區(qū)移動是由電子有效質(zhì)量減小引起的。對稱性以及邊界匹配性影響波函數(shù)，波函數(shù)隨著能量而變化，電子有效質(zhì)量變化引起波函數(shù)虛部的變化。

參考文獻(xiàn)

[1] Messiah? A.quantum Mechanics[M].North-Holland：Amsterdam，1961.

[2]Chang LL.EsakiL，TsuR.Resonant tunneling in semiconductor double barrier[J].App l.Phys.Lett，1974，24（12）：593-594.

[3]Chanda A，Eastman LF.Quantum mechanical reflect ion at triangular”planar2doped”potential barriers for transmission[J].J.Appl.Phys，1982，53：9165.

[4]AndoY，Itoh T.Calculation of transmission tunneling current across arbitrary potential barriers[J].J.Appl.Phys，1987，61（4）：1497-1502.

[5]Ban SL，Hasbun JE，Liang XX.A novelmethod for quantum transmission across arbitray potential barriers[J].J L um in，2000，87-89：369-371.

[6]Tuchman Judah A，Kim Sang sig，Sui Zhifeng，et al.Exciton photoluminescence in strained and unstrained ZnSe under hydrostatic pressure[J].Phys.Rev.B，1999，46（20）：13371-13378.