超晶格半导体材料的光磁电效应I精.docx
《超晶格半导体材料的光磁电效应I精.docx》由会员分享,可在线阅读,更多相关《超晶格半导体材料的光磁电效应I精.docx(16页珍藏版)》请在冰点文库上搜索。
超晶格半导体材料的光磁电效应I精
第26卷 第9期2005年9月
半 导 体 学 报
CHINESEJOURNALOFSEMICONDUCTORS
Vol.26 No.9
Sep.,2005
罗诗裕 男,1940年出生,教授,主要从事凝聚态物理与半导体超晶格方面的研究.Email:
Luoshy@
2004212209收到,2005205209定稿
ν2005中国电子学会
超晶格半导体材料的光磁电效应(Ⅰ
罗诗裕 邵明珠
(东莞理工学院,东莞 523106
摘要:
从Shockley2read统计出发,引入载流子寿命与浓度的相关性,描述了超晶格半导体载流子的输运特征,将载流子的输运方程化为二阶非线性方程,并用双参数摄动法找到了方程的一般解.在二阶近似下,计算了半导体材料的短路电流和光导电流,进一步揭示了大信号情况下光磁电效应的非线性特征.关键词:
超晶格;半导体;光磁电效应;非线性
PACC:
7240;7280
中图分类号:
O472 文献标识码:
A 文章编号:
025324177(20050921744205
1 引言
超晶格微结构是半导体物理与器件研究的前沿
领域之一.超晶格量子阱的主要特征是载流子(电子或空穴的运动在生长方向上受到限制,因此,量子阱的光电性质不同于体材料而呈现出许多新特点.超晶格量子阱的能带结构与体材料不同,而且有效质量的各向异性可以差几个量级,子能带的带隙还可以任意调节.超晶格量子阱低维结构是设计、制造下一代超大规模集成电路的物理基础.超晶格物理代表着半导体物理的一个全新层次.一方面由于组分、掺杂和结构可以在原子尺度上人为控制,这就为设计、制造新一代固体器件提供了技术基础;而且由于超晶格微结构中呈现出的新现象和新效应,又为它的应用提供了广阔前景.例如,利用量子阱中维度限制引起的激子吸收饱和现象制成通道快、高频性能好、能耗低、可保温工作的激子型光学双稳器件,以及开关速度可达ps量级的新一代量子器件等.
值得注意的是,广泛使用的半导体材料硅对微电子技术的发展作出过重大贡献,但是由于它的载流子迁移率低、能带结构为间接带隙,使得它的应用受到很大限制.前者使它的运行速度无法和GaAs相比;后者使它不能发出可见光.于是,如何改进硅基材料的性能以实现大规模的光电集成,一直是许
多科学家追求的目标.改进硅特性的方法有两种:
(1杂质工程,通过掺杂引入新的复合中心;(2能带工程,通过调整硅的能带结构(比如超晶格,使它的间接带隙变为准直接带隙,并使它的能带展宽到可见光的范围,而带边缘的弯曲则可以改变电子、空穴有效质量,从而使载流子迁移率上升.
超晶格的最大几何特点是它的晶格形变,而超晶格的光磁电特性又直接与它的形变有关,因此,通过控制或调节形变超晶格的形变参数就可以得到不同光磁电效应的半导体材料.我们曾对超晶格的沟道效应及其形变参数作过具体分析[1~6].注意到,要研究超晶格的光磁电效应就必须研究载流子输运.基于该点,本工作将分成两部分讨论,第一部分(本文首先讨论一般半导体材料中的光磁电效应;第二部分(下一篇文章将结果移植到超晶格,并考虑到超晶格材料的宏观周期性边界条件,进一步讨论超晶格光磁电效应的非线性特征.
在研究光磁电效应(下面简称PME效应的早期工作中,通常都假定载流子的寿命与它的浓度无关,因而描写载流子输运的连续性方程是一个线性二阶微分方程,结果对于光强比较弱(即低注入或低浓度情况符合得比较好.随着测量技术的不断提高,人们发现,即使在低注入情况下,PME效应也可观察到非线性特征;而在大注入情况下,则只有非线性微分方程才能描写.于是有人从经验出发,假定了不同形式的寿命2浓度相关性,从而改善了理论同实
第9期罗诗裕等:
超晶格半导体材料的光磁电效应(Ⅰ
验的拟合程度.然而,这类唯象分析缺乏必要的理论
依据.好在人们早就从Shockley—read统计出发导
出了载流子寿命τ与浓度Δξ/ξ0之间的依赖关系.
但因函数τ(Δξ/ξ0比较复杂,文献[7]仅将寿命τ按
浓度展开,保留一次项,从而将连续性方程化为含
(Δξ/ξ0
2
项的二阶非线性微分方程,成功地揭示了PME效应的非线性特征.同线性相比,允许的注入信号已有了明显的增加.但是,因它只保留了一次项,说明它仍然只能在比较小的信号下才成立.文献[8]作了进一步处理.但没有完全考虑载流子寿命的参数相关性,给结果带来了一定的局限性.基于上述考虑,本文直接从Shockley2read统计出发,引入载流子寿命与浓度的依赖性,把载流子的连续性方程化为二阶非线性微分方程.并用双参数摄动法找到了方程的一般解[9].在二级近似下,具体计算了半导体材料的短路电流和光导电流.讨论了这两个量与注入强度和吸收系数之间的关系,揭示了大注入情况下半导体的光磁电效应的非线性特征.
2 运动方程
在稳恒状态下,半导体中载流子输运满足如下
连续性方程
1q
ΔJ-Δξτ=-g(x,y,z(1其中
Δξ是载流子浓度(当Δξ=Δn时,表示过剩电子浓度;当Δ
ξ=Δp时,表示过剩空穴浓度;J是载流子密度矢量(J=Jn表示电子电流密度矢量;J
=Jp表示空穴电流密度矢量;τ是少子寿命(τ=τn是电子寿命;τ=τp是空穴寿命;q是载流子电荷(绝对值;g(x,y,z是载流子的产生率.
方程(1是一个三维问题,引入下列假设可将问题简化.为不失一般性,我们假定:
(1局部电中性,即半导体内部处处满足条件Δn=Δp;(2Hall角比较小;(3产生率只是一个坐标(比如y的函数.适当选择坐标系,使入射光沿y方向,磁场沿z方向,则短路电流沿x方向,于是,连续性方程(1化为
1q×dJdy-Δξτ
=-g(y
(2其中 J=Jy,且由公式
J=qD
d(Δξdy
(3给出,D是扩散系数.由Shockley—read统计,可导
出少子寿命τ与相对浓度Δξ/ξ0之间的关系
τ=
τ0(1+αΔξ
ξ0
1+(
ε+αΔξξ0
(4
其中
ε=τ0-τ∞
τ0(1+c(5α=
τ∞τ0(1+c
(6其中 τ0和τ∞分别表示无限小和无限大注入强度下的载流子寿命;c=p0/n0,而n0和p0分别是样品放在暗处的电子浓度和空穴浓度.
产生率可用如下指数规律表达:
g(y=Iκexp(-κy
(7
其中 I是入射光强度;κ是材料的吸收系数.将
(3,(4和(7式代入(2式,可得
D
d2
ξdy2-1
τ0
×
1+(
ε+αΔξ
ξ0
1+α
Δξξ0
Δξ=-Iκexp(-κy
(8
相应的边界条件为
D
d
Δξdy
y=0
=s
Δξ(0(9D
d
Δξdy
y=w
=0
(10
其中 s是表面复合速度;w是样品厚度.引入无量纲的量
W=w/L,Y=y/L,S=sw/D,u=Δξ/ξ0
(11
其中
L
2
=τ0D(12
由方程(8和(9可得无量纲的边值关系
d2udY2-1+(
ε+du1+αuu=Nexp(-μY(13
和
du
dY
Y=0
=au(0(14dudY
Y=W
=0
(15
其中
a=S/W, N=-IκL2
/(Dξ0,
μ=κL=K/W, K=κw
(16
且Y∈(0,W.
5
471
半 导 体 学 报第26卷
3 双参数摄动解
方程(13是一个复杂的二阶非线性微分方程,
不存在严格的解析解,且近似解也很难得到.我们曾用了几种近似方法作过尝试,发现只有摄动解的近似程度比较高,收敛速度快.下面我们用双参数摄动法进行求解.注意到少子寿命与参数ε和α有关,我们把它选为方程(13的独立变数,其目的是试图找到如下形式的双参数摄动解:
u(Y,α,
ε=∑
∞
i,j=0
αiεj
uij
(17
将上式代入方程(13,并注意到边界条件(14和
(15式,比较同次幂,可得零级解u00满足方程
u″
00-u00=Nexp(-μY
(18其中 u′=du/dY.一级解u01和u10满足方程
u″
01-u01=u2
00
(19
u″10-u10=u200+u00[Nexp(-μY-u″
00](20
二级解u02,u20和u11满足方程
u″
02-u02=2u01u00
(21u″20-u20=2u10u00+u10[Nexp(-μY-2u″
00]
(22
u″
11-u11=u01[Nexp(-μY+2u00]+
2u00u10-u00u″01-u01u″
00
(23
相应的边界条件由
u′ij|
Y=0=auij(0u′
ij|Y=w
=0
(24
给出,其中(i,j=0,1以及i=0,j=2;j=0,i=2.而通解uij满足方程
u″
ij-uij
=Hij(unm,u″
nm
(i,j≥0,而m,n小于i,j(25
边界条件由(24式给出.
一般解可表示为
uij=αijexp(Y+βijexp(-Y+Iij(Y
(26其中
αij=
-(1+aGij
(1+aexp(W-(1-aexp(-W
(27
βij=
-(1-aGij
(1+aexp(W-(1-aexp(-W
(28
Iij=
∫
Y
0Hij
(tsh(Y-tdt(29
Gij=
∫W
H
ij
(tch(W-tdt(30
由(26~(30式可求得零级解u00为
u00=α00exp(Y+β00exp(-Y+I00(Y(31其中
α00=
-(1+aG00
(1+aexp(W-(1-aexp(-W
(32
β00=
-(1-aG00
(1+aexp(W-(1-aexp(-W
(33
G00=
Nμμ2
-1ch(W-N
μ2-1
sh(W- Nμ
μ2
-1
exp(-μW
(34
I00=N
∫
Y
exp(-μtsh(Y-tdt
=
Nμ
μ2
-1sh(Y+Nμ2-1
exp(-μY- N
μ2-1
ch(Y
(35一级解u01为
u01=α01exp(Y+β01exp(-Y+I01(Y(36其中 一级解u10为
u10=0
(37 由于一级解还不能反映参数α的影响.必须进
一步找出系统的二级解,它的一般表达式为
u02=α02exp(Y+β02exp(-Y+
2
∫
Y
u01(tu
00
(tsh(Y-tdt
u20=0
(38
u11=α11exp(Y+β11exp(-Y+
∫
Y
[-u3
00(t]sh(Y-tdt
(39
令(27和(28式中的i=0,j=2,即可得(38式中的α02和β02;令(27和(28式中的i=j=1,即可得到(39式中的α11和β11(因太冗长,具体表达式不再一一列出.由(31,(36~(39式,可将二阶的双参数摄动解表示为
u=u00+εu01+ε2
u02+αεu11(40可以证明,当0≤Y≤W时,级数(17式绝对一致收敛,只要条件
|ε|,|α|<
1
4u
max
(Y
(41
满足.下面我们将在二级近似下讨论PME效应的短路电流和光导电流.
6
471
第9期罗诗裕等:
超晶格半导体材料的光磁电效应(Ⅰ
4 结果和讨论
在利用PME效应测定载流子寿命和表面复合
速度时,常常需要测定材料的光导电流IPc.相应的光电导由公式
ΔG=q(μn+μp
∫
W
0Δξ(ydy
(42
给出,注意到u=Δξ/ξ0,上式可改写为
ΔGR=ΔG
q(μn+μpξ0L
=∫
W
0u(YdY
(43
其中 u(Y由(17式给出,在二级近似下由(40式
给出.将(40式代入(43式,完成积分可得材料的
“光电导”
ΔGR.为了比较,我们选择了与文献[7,8]相同的一组参数.图1给出了吸收系数K=10时,“光电导”
ΔGR与入射光强度N0之间的关系.选择了三组不同参数:
ε=α=0;ε=-10,α=0和ε=-12,α=5,其中ε=α=0描述的是线性情况.由图1可以看出,当浓度N0<10-3时(对应于光强比较弱,线性近似比较好.但是,随着光强的增加,比如
N0>5×10-3时,“光电导”
ΔGR表现出了明显的非线性特征.计算表明,当材料的吸收系数增加时,非线性更加明显.可见,对于强吸收样品,载流子寿命与入射光强度的相关性变得十分重要
.
图1 “光电导”
ΔGR与入射光强度N0之间的关系Fig.1 RelationshipbetweenΔGRandN0
根据过剩载流子沿y方向的分布函数u(y,可
将PME短路电流Isc表示为
Isc=-qD(μn+μpB[
Δξ(w-Δξ(0](44其中 μn和μp分别为电子和空穴迁移率.注意到u=Δξ/ξ0,上式可化为
IR=
Isc
qD(μn+μpB
=ξ0[u(0-u(W](45
对于电子,u=Δn/n0,ξ0=n0,u(Y由(17式给出,在二级近似下由(40式给出.计算表明,在强注入情
况下,短路电流IR同样表现出与“光电导”ΔGR类似的非线性特征.
结果表明,引入载流子寿命与浓度的相关性,可
以很好地描述半导体载流子的输运特征,并用双参数摄动法找到了二阶非线性输运方程的一般解,在二阶近似下,计算了半导体材料的光导电流和短路电流,进一步揭示了大信号情况下半导体光磁电效应的非线性特征[7,8].参考文献
[1] LuoSY,ShaoMZ.Dislocationmodelforstrainedsuperlat2
ticeanddechannellingeffectsofaparticles.ChineseJournalofSemiconductors,2003,24(5:
485(inChinese[罗诗裕,邵明
珠.形变超晶格的位错模型与粒子的退道效应.半导体学报,
2003,24(5:
485]
[2] LuoSY,ShaoMZ.Channellingeffectsforstrainedsuperlat2
ticeandphaseplanarcharacteristicsofsystem.NuclearPysicsReview,2002,19(4:
407(inChinese[罗诗裕,邵明珠.形变
超晶格的沟道效应与系统的相平面特征.原子核物理评论,
2002,19(4:
407]
[3] ShaoMZ,LuoSY.Sine2squaredpotentialandmotionbehav2
iourofchannelingparticlesforstrainedsuperlattice.ChineseJournalofSemiconductors,1993,14(4:
353(inChinese[邵
明珠,罗诗裕.正弦平方势与形变超晶格的沟道特征.半导体学报,1993,14(4:
353]
[4] LuoSY,ShaoMZ,HuXD.Averagefieldideaandsingle
particlemodelfor22dimensioncrystallizationbeams(Ⅰ.HighEnergyPhysicsandNuclearPhysics,2004,28(1:
96(inChinese[罗诗裕,邵明珠,胡西多.二维晶化束的平均场概念与单粒子模型(Ⅰ
.高能物理与核物理,2004,28(1:
96][5] LuoSY,TanYM,ShaoMZ,etal.Motiondampinginchan2
nelingeffecsandchaoticbehaviourofsestem.ActaPhysicaSinica,2004,53(4:
1157(inChinese[罗诗裕,谭永明,邵明
珠,等.沟道效应的运动阻尼与系统走向混沌的临界特征.物
理学报,2004,53(4:
1157]
[6] LuoSY,ShaoMZ,DengCL.Sine2squaredpotentialandcha2
oticbehaviouforstrainedsuperlattice.ChniseJournalofSem2
7
471
半 导 体 学 报第26卷
iconductors,2005,26(2:
299(inChinese[罗诗裕,邵明珠,
邓成良.正弦平方势与形变超晶格的混沌行为.半导体学报,
2005,26(2:
299]
[7] AugelliV,VasanelliL,LeoM,etal.Atheoryofthephoto2
magnetoelectriceffectwithinjection2level2dependentlifetime.
JApplPhys,1980,51:
2784
[8] LeiseS.Photo2conductanceofGaAs.GSIReport,2001:
1
[9] NayfehAH.Introductiontoperturbationtechniques.John
Wiley&Sons,1981
Photo2Magneto2ElectricEffectsforStrainedSuperlattice(Ⅰ
LuoShiyuandShaoMingzhu
(DongguanUniversityofTechnology,Dongguan 523106,China
Abstract:
Theequationforcontinuityofthechargecarriesthroughasemiconductorsuperlatticeisreducedtoasecondordernonlineardifferntialequationbyusingtheinjection2leveldependentlifetimederivedfromShockley2readstatistics.Thegeneralsolutionisfoundbyusingthetwo2parametersperturbationmethod.ThePMEshort2circuitcurrentIscandthephoto2conduct2anceΔGinasemiconductorarecalculatedinthesecondoderapproximation.Thenonlinearpropertiesofthephoto2magneto2elec2triceffectinthecaseofalargesignalarefurtherreversed.
Keywords:
superlattice;semiconductor;photo2magneto2electriceffect;nonlinearity
PACC:
7240;7280
ArticleID:
025324177(20050921744205
LuoShiyu male,wasbornin1940,professor.Heengagedintheresearchoncondensationphysics.Email:
Luoshy@
Received9December2004,revisedmanuscriptreceived9May2005ν2005ChineseInstituteofElectronics8471