半导体外延层晶格失配度的计算精.docx
- 文档编号:26706455
- 上传时间:2023-06-21
- 格式:DOCX
- 页数:18
- 大小:275KB
半导体外延层晶格失配度的计算精.docx
《半导体外延层晶格失配度的计算精.docx》由会员分享,可在线阅读,更多相关《半导体外延层晶格失配度的计算精.docx(18页珍藏版)》请在冰豆网上搜索。
半导体外延层晶格失配度的计算精
第30卷第1期
2006年2月南昌大学学报(理科版
JournalofNanchangUniversity(NaturalScienceVol.30No.1
Feb.2006
文章编号:
1006-0464(200601-0063-05
半导体外延层晶格失配度的计算
何菊生,张 萌,肖祁陵
(南昌大学材料科学与工程学院,江西南昌 330047
3
摘 要:
比较了晶格失配度的各种定义,建议统一使用同一定义外延生长层和衬底的二维晶格失配度的计算,最后讨论了结合,正三角形晶格和长方形晶格匹配,,宽比接近匹配比时,整体匹配较好。
关键词:
晶格失配度;体结构;;:
12 (,即晶格错配度,是描述衬底和外延膜晶格匹配的参量。
晶格失配和热膨胀系数失配不同程度地影响晶体的外延生长,在外延层中产生大量缺陷,甚至无法生长单晶,影响器件的性能和寿命。
对晶体材料的失配度研究也就成为人们最为关心的热点之一。
然而,各种文献并没有采用统一的失配度定义,导致对于同样的衬底和外延层却有不同的失配度。
例如,GaN在蓝宝石衬底上的失配度有16%15%
[3]
[1]
就有三个不同的失配度。
对六方GaN在六方蓝宝石衬底上的失配度分别用这三个定义计算,便得到16%,13.8%,15%的结果。
用定义1和定义2计算的结果相差较大时,有的文献便采用较小的值。
此外,材料的晶格常数随材料本身的缺陷和掺杂、衬底的选择、生长工艺不同而不同,不同的文献常有不同的实验值,但这些对失配度的影响不大,主要还是采用不同的定义带来的失配度大小的差别。
一般认为,|δ|<5%为完全共格,|δ|=5%~25%为
[7]
半共格界面,|δ|>25%完全失去匹配能力。
14%
[2]
以及
等不同说法。
因此,了解失配度的定义,掌
握其计算方法,在研究工作中显得非常重要。
2 晶格失配度的计算
211 简化模型
1 晶格失配度的定义
半导体外延膜与衬底间的失配度有三种定义。
设无应力时衬底的晶格常数为as,外延薄膜的晶格
常数为ae,晶格失配度为δ,
δ定义1 1δ定义2 2a-aaeas-ae
as
[6][5][4]
绝大部分半导体材料具有金刚石结构、闪锌矿
结构(β相或纤锌矿(α相结构。
金刚石和闪锌矿结构是由两套面心立方结构沿对角线方向平移1/4对角线长套构而成的,纤锌矿结构则是由两套六方密堆积结构沿C轴平移5/8C套构而成的。
结构虽然复杂,但是失配度只涉及生长面上衬底和外延膜的两个二维晶格面,因此为简化,暂不考虑它们的体结构。
金刚石与闪锌矿结构的(111面和纤锌矿结构的(0001面或(0002面原子排列完全一样,最小的二维晶格都是正三角形,我们不妨称其为正三角形晶格,而金刚石与闪锌矿结构的(100面最小二维晶格都是正方形,我们称其为正方形晶格,如图1所示。
此外,还有正交晶系和正方晶系的长方形晶格等。
设晶格常数为a,正三角形晶格的原子列原子排
2|as-ae|δ定义3 3
as+ae
定义1的分母是外延膜的晶格常数,这是许多文献采用的定义。
也有不少文献采用定义2,用衬底的晶格常数做分母。
定义3的分母既不是衬底,也不是外延膜,而是两者的平均晶格常数,这是求半导体异质结的失配度定义,有些文献也将它用于一般的外延膜和衬底间。
这样,同样的外延膜和衬底
收稿日期:
2005-09-20
作者简介:
何菊生(1966-,男,硕士生;3通讯作者:
张 萌(1961-,女,教授,博士生导师,TieGang_zm@126.com1
・64・南昌大学学报(理科版2006年
二维晶格结构相同时,某一晶向的失配度和相应的垂直方向的失配度相同,即晶格匹配相同,只要这一方向晶格匹配,相应的垂直方向就自然匹配。
一般地,各文献上的失配度基本上是指衬底和外延膜二维晶格相同时的失配度。
计算失配度的方法是:
①晶面匹配。
最常见的是:
(0001hex∥(0001hex(两正三角形晶格匹配,(100fcc∥(100fcc(两正方形晶格匹配;其次是:
(0001hex(111fcc(;最后
图1
半导体外延层生长面上的两种最重要的二维晶格
列周期如图2,周期有a,a,a2a/22a,如图3。
是,,
。
晶格25%时,直接确定匹配晶向,〈1120〉∥〈1120〉,〈100〉∥〈100〉,〈1120〉∥〈110〉。
否则,晶向匹配使晶格30°旋转(正三角形晶格间匹配或45°旋转(正方形晶格间匹配,即〈1120〉∥〈1100〉,〈100〉∥〈110〉,这样有助于减小失配度,降低界面能。
③计算失配度。
一般情况下,直接根据晶格常数计算;出现原子列旋转时,其中一晶格常数须改成能和另一原子排列周期相匹配的原子排列周期。
对于性质不同的晶面匹配,两垂直方向的失配度不同。
例如,GaN属六方结构时,原子排列周期可为
a,a,2a,即3.185,5.517,6.370;而为立方结构
图2 正三角形晶格的原子排列周期
时,原子排列周期可为a/2,a,2a,即3.21,4.54,
6.42。
它和部分衬底的失配度可列表计算如表1:
(其中未注明文献来源的为作者的计算值
需要指出的是,失配度必须在已知外延层和衬
图3 正三角形晶格的原子排列周期
212 衬底和外延膜的二维晶格相同时的失配度
绝大多数情况下,外延生长六方(立方晶体时常采用体结构相同的六方(立方晶体,这是因为生
长面上的二维晶格都是正三角形(正方形时晶格容易匹配。
外延生长六方晶体还可选择立方晶体的(111面,也是这个原因。
衬底和外延层接触面的
外延层/衬底
GaN/Al2O3GaN/ZnOGaN/6HSiCGaN/SiGaN/SiGaN/GaAs
底的晶体结构下才能确定。
GaN等氮族化合物半导体材料既可有闪锌矿结构又可有纤锌矿结构。
研究表明,外延层的初始稳定晶向主要取决于外延材料与衬底材料的晶体结构及跨越界面的化学键的特
[12]
性,而不依赖于具体的生长过程。
实验表明,生长在(111Si面或(111GaAs面上的GaN基本上是纤
(晶格常数单位:
10-10m
表1 GaN在常见衬底材料上的失配度
晶面匹配
(0001GaN∥(0001Al2O3(0001GaN∥(0001ZnO(0001GaN∥(0001SiC(0001GaN∥(111Si(100GaN∥(100Si(100GaN∥(100GaAs
原子排列周期
5.517/4.7583.185/3.2523.185/3.0803.185/3.8394.54/5.434.54/5.653
晶向匹配
〈1100〉〈1120〉GaN∥Al2O3〈112〉〈112〉GaN∥ZnO〈1120〉〈1120〉GaN∥ZnO〈1120〉〈110〉GaN∥Si〈100〉〈100〉GaN∥Si〈100〉〈100〉GaN∥GaAs
δ1
16%
[1][8][9]
δ2
14%
[2]
2.1%3.4%17%
2.1%3.3%21%20%25%
[10]
17%20%
[11]
第1期 何菊生等:
半导体外延层晶格失配度的计算・65・
锌矿结构,而生长在(001Si面或(001GaAs面上的GaN基本上是闪锌矿结构,即c-GaN。
也就是说,衬底的结构对生长出的外延层的结构有极其重要的影响。
从上面表中容易发现,随失配度的增大,δ1和
δ-2的差也增大。
这可以从它们的关系图线δδ21
δ=1直观看出。
δ1、2可同时为正(差表示为δ2-δ1或同时为负(差表示为|δ2|-|δ1|
。
图中失配δ度之差δ2-δ1和|2|-|δ1|都随δ1的增大而增大。
δ在δ1<5%时,δ1、2近似相等;δ1<10%时,相差1%以内;δ1<15%时,相差3%以内,5%。
因此,的δ1。
原子列方向旋转了30°,晶向匹配为〈1010〉ZnO∥
〈100〉=Si,该方向上失配度δ1=
5.632
-3.6%。
生长面内垂直方向上,〈1210〉ZnO∥〈010〉Si,该方向上ZnO的原子排列周期分别是是原方向的3倍,即5.632(3.2523,9.756(5.632。
原子排列周期匹配为:
9.756/10.860,失配度δ′=11%。
1=
9.ZnOSi(按(0001〈101110ZnO18%,本文ZnO在Si(100面Si(111面上生长的原因。
2.3.2 长方形晶格与正三角形晶格匹配
β-FeSi2属正交晶系,a=0.9879nm,b=0.7799nm,c=0.7939nm,在Si衬底上β-FeSi2单晶的XRD图谱上同时出现了Si(111和β-Fe2
[14]
Si2(220或(202的衍射峰,说明了它们的晶面匹配是(220β-FeSi2∥(111Si,或(202β-FeSi2∥
(111Si。
①(220β-FeSi2∥(111Si时,(220β-FeSi2长方形晶格的原子排列周期(边长为0.7839nm,
22
.9879+0.7799=1.259(nm,而能与之相匹
图4
配的Si原子排列周期为0.5432=0.7678(nm,垂直方向为0.7678(nm,即1.330nm,周期匹配为0.7839/0.7678,1.259/1.330,所以(220〈001〉〈110〉β-FeSi2∥(111Si,
δ=-2.1%1a0.7839
在垂直方向上,(220〈110〉β-FeSi2∥(111〈112〉Si,
2
2
2
2
213 二维晶格不同时的失配度
尽管材料生长优先考虑体结构相同或生长面二维晶格相同的衬底,但是,我们仍经常看到材料在体
结构不同或生长面二维晶格不同的衬底上生长的报
[13]
道,例如,孙一军等人在立方GaAs(100衬底上制备出了单相六方GaN薄膜,李春延等人在硅衬底上制备出了正交晶系的β-FeSi2单晶,贺洪波小组
[15]
[14]
和叶志镇小组
[16]
分别在Si(100上用磁控
溅射法制备出了质量优良的ZnO薄膜,等等。
下面我们分别分析ZnO/Si(100和β-FeSi2/Si的失配度。
2.3.1 正方形晶格与正三角形晶格匹配
由于ZnO是六方晶体,ZnO/Si(100生长面上ZnO是正三角形二维晶格,所以这是正三角形晶格和正方形晶格匹配。
晶面匹配为(0001ZnO∥(100Si。
ZnO的原子排列周期分别是3.252〈1120〉,5.632(即3.252〈1010〉,能和Si的晶格常数5.(ZnO的原子排列周期匹配的是5.632,430〈100〉
δ′1a.9879+0.7799
=5.7%
②同理,(202〈010〉β-FeSi2∥Si(111〈101〉Si,δ.6%,1b=-1
在垂直方向上,(202〈101〉β-FeSi2∥Si(111〈112〉.5%Si,δ1b=5
显然,β-FeSi2在Si(111面生长,晶格匹配很好,有利于生长高质量的材料。
下面是GaN在低晶格失配材料上的失配度。
其中LiAlO2是正方晶体,晶格常数a=5.170,c=6.260;LiGaO2是正交晶体,晶格常数a=5.402,b=6.
・66・
-10
南昌大学学报(理科版2006年
372,c=5.007,单位都是10m。
由正三角形晶格的原子排列规律可知,最适宜与之匹配的长方形晶格的长宽比为1.732和1.155(即2/
3/,我们不妨称之为匹配长宽比。
上面β
-FeSi2晶胞的对角面的长宽比为1.274,表中LiA2lO2和LiGaO2则分别为1.211和1.180。
从失配度
外延层/衬底
GaN/LiAlO2
角度看,低失配长方形晶格衬底的选择应同时遵循
两个原则:
①长方形长宽比应接近匹配长宽比1.732或1.155;②宽应接近正三角形晶格的最近邻或次近邻原子间距。
2.3.3 外延层与蓝宝石R面匹配
R面(0112是个超长矩形组成的二维晶格,是
匹配
失配度δ1
.26.1.7%6.8%
表2 GaN在重要低衬底材料上的失配度
体结构六方/正方六方/生长面二维晶格正三角形/长方形(0001〈1120〈001〉GaN∥(LilO2(0001〈〈010GaN(LiAlO2〈1120〉GaNLiGaO2(1〉〈100〉GaN∥(001LiGaO2
GaN/LiGaO2
6.370/6.3720.03%5.516/5.402
2.1%
较大,,也可能生长
[17]
出比C面(0001质量更好的纤锌矿GaN。
GaN的优先生长面是(0001,晶面匹配是(0001GaN∥(0112Al2O3,R面宽方向为〈2110〉,原子排列周期为4.785,相应的GaN为3.1853,方向为〈1010〉,因此,〈1010〉〈2110〉GaN∥Al2O3,δ1=13.3%。
垂直方向上,〈1210〉〈0111〉GaN∥Al2O3,周期匹配为(3.185×5/15.384,由此得失配度δ′1=3.4%。
从晶格失配度角度来说,比和C面匹配更
图5
(200、θ=41.5°Si(400的两个衍射峰分别位于2
好。
θ=69.2°和2处,2
3 晶格失配度的测定
θ311 利用2扫描
当衬底和外延层都为立方结构时,失配度的大
θ扫描图谱来确定。
方法小可以通过X射线衍射2
是从图谱中找出两匹配晶面所对应的衍射峰(或多
θ级衍射峰,根据这两个衍射峰对应的2即可确定。
设衬底和外延层的晶格常数分别为a1和a2,衍射峰对应的晶面间距分别为d1和d2,由布拉格方程
θλ,2d2sinθλ,对于一级衍射峰,失配有2d1sin1=2=度的大小为
θθa1-a2d1-d2sin2-sin1
δ===
θa2d2sin1
θθ特别地,当θ,则δ=1、2相差不大时,设θ1-θ2=Δ
θθsin2-sin1
θθ≈-Δ・ctg1。
θsin1对于多级衍射峰,如图5中3C-SiC在Si衬底
上外延,得到3C-SiC的二级衍射峰和Si的四级衍射峰,也很容易确定它们的失配度。
图中3C-SiC
d1d2
θλ,2sinθλ,其中d1sin1=2=24
和d2是3C-SiC(100、Si(100对应的面间距。
δ
a1-a2d1-d220%。
与理论值基本一致。
a2
d2
失配度超过25%时,要考虑晶格旋转,根据匹配的原子排列周期重新确定。
θΦ扫描312 利用2扫描、
一般情况下,外延层和衬底的晶面和晶向匹配,都需通过XRD四圆衍射仪测绘出极射赤平投影来
[19]
确定。
如图所示,是纤锌矿GaN在未经氮化或氮化不充分的的蓝宝石衬底上外延得到的部分极图。
图(a(0001GaN∥(0001Al2O3,显然晶面旋转了30°,〈1010〉〈2110〉GaN∥Al2O3,这是两正三角形晶格匹配,δ1=
3.1853
-10
=13.7%。
图(b外
-10延层(1012GaN是长方形晶格,长7.571×10m,
宽3.185×10
m。
晶面匹配为(1012GaN∥
(0001Al2O3,晶向匹配为〈1210〉〈0110〉GaN∥Al2O3,
第1期 何菊生等:
半导体外延层晶格失配度的计算・67・
周期匹配为(3.185×3/(4.785
δ1=
=13.3%;垂直方向上晶向匹
3.185×3
热膨胀失配将严重地影响材料质量,外延层较厚时
甚至产生龟裂,热失配度的研究也成为了极重要的课题。
尽管如此,从理论和实验上研究材料的晶格失配度仍然是材料生长技术研究的中心环节之一。
气相外延生长要求失配度小于10%,液相外延生长要求失配度小于1%,而用于光电器件的异质结的失配度更应小于0.1%。
超薄应变层、缓变组分层一定程度上解决了晶格失配问题,但彻底解决GaN、ZnO。
配为〈1011〉GaN∥〈2110〉Al2O3,周期匹配为((1012GaN宽/(4.785×2,所以,δ′1=22
(3.1853+5.185
2
2
=-26.4%,该
晶向匹配很差,但没有影响密排列(宽方向和衬底
:
[1]SW,LeeSS,etal.TheGrowthofβ-Li2
GaO2FilmsUsingNovelSinglePrecursors.JournalofCrystalGrowth,2001(226:
481-487.
[2] LeeJJ,ParkYS,YangCS,etal.MBEGrowthofWurtz2
iteGaNonLaAlO3(100Substrate.JournalofCrystalGrowth,2000(213:
33-39.
[3] 沈晓明,付 羿,冯 淦,等.用侧向外延法降低立方
相GaN中的层错密度[J].半导体学报,2002,23(10:
1093.
[4] 吴自勤,王兵主.薄膜生长[M].北京:
科学出版社,
2003.第六章,262,122,341.
[5] 冯 端.固体物理学大辞典[M].北京:
高等教育出版
图6 GaN/Al2O3(0001匹配关系的极图
社,1995:
2,566.
[6] 余金中,王杏华.异质结构和量子结构.物理[J],
2001,30(3:
171.
[7] 潘金生,仝健民,田民波.材料科学基础[M].北京:
清
匹配。
大量事实说明,正三角形晶格和长方形晶格匹配,长方形晶格的宽列原子的匹配具有优先性,不受长列原子匹配的影响。
这是因为随着化学键的增长,原子间的化学作用迅速减弱,长列原子对宽列原子影响受到很大限制。
华大学出版社.2004:
415.
[8] 李庚伟,吴正龙,邵素珍,等.氧离子束辅助激光束淀
积生长ZnO/Si的研究[J].材料导报,2005,19(2:
109.
[9] LuJ,HaworthL,WestwoodDI,etal.InitialStagesofMo2
lecular-beamEpitaxyGrowthofGaNon6H-SiC(0001.ApplPhysLett,2001,78(8:
1080.
[10] Seong-Hwanjiang,Seung-JaeLee,Iu-SeokSeo,et
al.CharacteristicsofGaN/Si(111EpitaxyGrownUsingAl0.1Ga0.9N/AlNCompositeNucleationLayersHavingDifferentThicknessesofAlN.JournalofCrystalGrowth241(2002:
289-296.
[11] 韩培德,杨海峰,杨 辉,等.分子束外延GaN/GaAs
4 结束语
本文系统地分析了失配度的基本计算方法,对
于三元或四元固溶体半导体同样适用。
计算的关键之一在于晶向匹配,特别是有否晶格旋转。
实验表明,晶格旋转与晶格常数、外延层厚度、生长温度等都有关。
这个问题有时很复杂,例如,ZnO/Al2O3中,除出现取向关系为〈1120〉〈1010〉ZnO∥Al2O3的主畴外,还会出现〈1120〉〈1120〉旋转ZnO∥Al2O3的30°畴;再如,cGaN/GaAs,如出现45°的晶格旋转,则失配度从20%降为13%,或从25%降为11%,有利于降低界面自由能,但至今未见有关报道。
半导体材料生长是一项系统工程,失配度研究在其中占有一席之地。
由于环保、经济、技术等方面的原因,晶格失配始终未能彻底解决。
材料的生长通常在高温下进行,降至常温时外延层和衬底间的
[20]
立方相异质结构的电子显微分析.半导体学报[J],
1999,120(1:
58.
[12] 冯 端,师昌绪,刘治国.材料科学导论.北京:
化学
工业出版社,2002:
588.
[13] 孙一军,李爱珍,齐 鸣.立方GaAs(100衬底上制
备单相GaN薄膜[J].半导体学报,2001,3,22(3:
315.
(下转第102页
・102・南昌大学学报(理科版2006年[10] 唐启义,冯明光.实用统计分析及其DPS数据处理系统[M].北京:
科学出版社,2002:
626-627.ApplicationGreySystemTheorytoForecastingandAnalysisonAffectecFactorsofAtmosphericNO2Concentration——TakeChangchenDistrictFoshanCityasExample—(DepartmentofResourceandEnvironment,FoshanUniversity,Foshan528000,China参考文献:
ingtheatosphericNO2concentrationofCheanchendistrictareformthedensityofmotorvehicleandtheamountofmmotorvehicle.AndwiththeGrayModelG(1,1,theatmosphericNO
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 半导体 外延 晶格 失配 计算