半导体专业试验补充silvaco器件仿真精.docx
《半导体专业试验补充silvaco器件仿真精.docx》由会员分享，可在线阅读，更多相关《半导体专业试验补充silvaco器件仿真精.docx（25页珍藏版）》请在冰豆网上搜索。
半导体专业试验补充silvaco器件仿真精
实验2PN结二极管特性仿真
1、实验内容
（1PN结穿通二极管正向I-V特性、反向击穿特性、反向恢复特性等仿真。
（2结构和参数:
PN结穿通二极管的结构如图1所示，两端高掺杂，n-为耐压层低掺杂，具体参数：
器件宽度4卩m器件长度20卩m耐压层厚度16卩m,p区厚度2卩m,n区厚度2^m。
掺杂浓度：
p+区浓度为1xi019cm-3,n+区浓度为1K1019cm-3,耐压层参考浓度为5X1015cm-3。
0W
p'
n
■
+n
图1普通耐压层功率二极管结构
2、实验要求
（1掌握器件工艺仿真和电气性能仿真程序的设计
（2掌握普通耐压层击穿电压与耐压层厚度、浓度的关系
3、实验过程
#启动Athena
goathena
#器件结构网格划分
linexloc=0.0spac=0.4
linexloc=4.0spac=0.4
lineyloc=0.0spac=0.5
lineyloc=2.0spac=0.1
lineyloc=10spac=0.5
lineyloc=18spac=0.1
lineyloc=20spac=0.5
#初始化Si衬底；
initsiliconc.phos=5e15orientation=100two.d
#沉积铝；
depositalumthick=1.1div=10
#电极设置
electrodename=anodex=1
electrodename=cathodebackside
#输出结构图
structureoutf=cb0.str
tonyplotcb0.str
#启动Atlas
goatlas
#结构描述
dopingp.typeconc=1e20x.min=0.0x.max=4.0y.min=0y.max=2.0uniform
dopingn.typeconc=1e20x.min=0.0x.max=4.0y.min=18y.max=20.0uniform
#选择模型和参数
modelscvtsrhprint
methodcarriers=2
impactselb
#选择求解数值方法
methodnewton
#求解
solveinit
logoutf=cb02.log
solvevanode=0.03
solvevanode=0.1vstep=0.1vfinal=5name=anode
#画出IV特性曲线
tonyplotcb02」og
#退出
quit
图2为普通耐压层功率二极管的仿真结构。
正向I-V特性曲线如图3所示，导通
电压接近0.8V
图2普通耐压层功率二极管的仿真结构
B
图3普通耐压层功率二极管的正向I-V特性曲线
运用雪崩击穿的碰撞电离模型，加反向偏压，刚开始步长小一点，然后逐渐加大步长。
solvevanode=-0.1vstep=-0.1vfinal=-5name=anode
solvevanode=-5.5vstep=-0.5vfinal=-20name=anode
solvevanode=-22vstep=-2vfinal=-40name=anode
solvevanode=-45vstep=-5vfinal=-240name=anode
求解二极管反向IV特性，图4为该二极管的反向I-V特性曲线。
击穿时的纵向电场分布如图5
所示,最大电场在结界面处，约为2.5X105V?
CM,在耐压层中线性减小到80000
V?
cm-1o
ULAS
图4普通耐压层功率二极管的反向I-V特性曲线
ATLA^
9-
9-
SrSop1dr
图5普通耐压层功率二极管击穿时的电场分布
导通的二极管突加反向电压，需要经过一段时间才能恢复反向阻断能力。
电路图如图6所示。
设t=0前电路已处于稳态，ld=IfO。
t=0时，开关K闭合，二极管从导通向截止过渡。
在一段时间内，电流Id以diO/dt=-Ur/L的速率下降。
在一段时间内电流Id会变成负值再逐渐恢复到零。
仿真时先对器件施加一个1V的正向
偏压,然后迅速改变电压给它施加一个反向电压增大到2V。
solvevanode=1
logoutf=cj2_1」og
solvevcathode=2.0ramptime=2.0e-8tstop=5.0e-7tstep=1.Oe-1O
反向恢复特性仿真时，也可以采用如图7的基本电路，其基本原理为：
在初始时刻,电阻R1的值很小，电阻R2的值很大，例如可设R1为1X10-3Q,R2为1X106Q电感
L1可设为3nH;电压源及电流源也分别给定一个初始定值v1,i1;那么由于R2远大于
R1,则根据KCL可知，电流i1主要经过R1支路,即i1的绝大部分电流稳定的流过二极管，二极管正向导通,而R2支路几乎断路，没有电路流过。
然后，在短暂的时间内，使电阻R2的阻值骤降。
此时，电阻器R2作为一个阻源，其阻值在极短的时间间隔内以指数形式从1X106◎下降到1X10-3Q。
这一过程本质上是使与其并联的连在二极管阳极的电流源i1短路，这样电流i1几乎全部从R2支路流过,而二极管支路就没有i1的分流，此刻电压源v1开始起作用，二极管两端就被施加了反偏电压，由于这些过程都在很短的时间内完成，因而能够很好的实现二极管反向恢复特性的
模拟。
反向恢复特性仿真图如图8所示,PN结功率二极管的反向恢复时间约为
50ns。
图6反向恢复特性测试原理电路图
11疋电懺淤1
二极管
独立电压源V1
图7二极管反向恢复特性模拟电路图
Zk
a>
AKAS
AnodeCkJirerrt（A）
图8器件反向恢复特性曲线
实验3PN结终端技术仿真
1、实验内容
由于PN结在表面的曲率效应，使表面的最大电场常大于体内的最大电场，器件的表面易击穿，采用终端技术可使表面最大电场减小，提高表面击穿电压。
场限环和场板是功率器件中常用的两种终端技术。
场限环技术是目前功率器件中被大量使用的一种终端技术。
其基本原理是在主结表面和衬底之间加反偏电压后，主结的PN结在反向偏压下形成耗尽层，并随着反向偏置电压的增加而增加。
当偏置电压增加到一定值是，主结的耗尽层达到环上，如图1所示，这样就会使得有一部分电压有场环分担，将主结的电场的值限制在临界击穿电压以内，这将显著的减小主结耗尽区的曲率，从而增加击穿电压。
阳极
图1场限环
场板结构在功率器件中被广泛应用。
场板结构与普通PN结的区别在于场板结
构中PN区引线电极横向延伸到PN区外适当的距离。
而普通PN结的P区引线电极的横向宽度一般不超过P扩散区的横向尺寸。
PN结反向工作时,P区相对于N型衬底加负电位。
如果场板下边的二氧化硅层足够厚，则这个电场将半导体表面的载流子排斥到体内，使之表面呈现出载流子的耗尽状态，如图2所示，就使得在同样电压作用下，表面耗尽层展宽，电场减小,击穿电压得到提高。
2、实验要求
（1场限环特性仿真
场限环:
击穿电压200V,设计3个环，环的宽度依次为6、5、5、5卩m间距为4、
5、6卩m外延层浓度为1X1015cm-3,观察表面电场。
（2场板特性仿真
场板:
氧化层厚度111m结深111m场板长度分别为0叩、2叩、4叩、6叩、8叩、101m外延层浓度为1X1015cm-3观察表面电场。
005“E
1*cm1C.W\im
图2场板
3、场板的应用实例：
场板对大功率GaNHEMT击穿电压的影响
（1内容
（aGaNHEMT的工作机理、击穿特性刻画以及对场板结构的GaNHEMT击穿特性的进行仿真分析。
（b结构和参数：
场板结构的GaNHEMT的结构尺寸及掺杂浓度如图3所示
'H
Insulator.对
nAIGaN1x10ie/
1x101B/cm30,04|im
0.05um
n-GaN1x1015/cm30.18ym
图3场板结构的大功率GaNHEMT
（2要求
（a掌握定义一个完整半导体器件结构的步骤，并能对其电性能进行仿真研究。
（b理解场板技术对器件击穿电压提高的作用原理并能结合仿真结果给出初步
分析。
（3实验过程
#启动internal,定义结构参数
#场板长度从1um增大到2.25um步长为0.25um通过改变I取值来改变场板长度setl=1.0
#drain-gatedistanee
setLdg=5.1
#fieldplatethickness
sett=1.77355
#AIGaNcompositionfraction
setxc=0.295
#settraplifetime
setlt=1e-7
setlight=1e-5
#meshlocationsbasedonfieldplategeometry
setxl=0.9+$l
setxd=0.9+$Ldg
sety1=0.3+$t
sety2=$y1+0.02
sety3=$y2+0.04
sety4=$y2+0.18
#启动二维器件仿真器
goatlas
meshwidth=1000
#网格结构
x.ml=0.0s=0.1
x.m1=0.05s=0.05
x.ml=0.5s=0.05
x.m1=0.9s=0.025
x.ml=（0.9+$xl/2s=0.05
x.ml=$xls=0.025
x.ml=（$xl+$xd/2s=0.25
x.ml=$xd-0.05s=0.05
x.ml=$xds=0.05
#
y.ml=0.0s=0.1000
y.ml=0.3s=0.1000
y.ml=$y1s=0.0020
y.ml=$y2s=0.0020
y.ml=$y3s=0.0100y.ml=$y4s=0.0500
#devicestructure
#POLAR.SCALEischosentomatchcalibratedvalues
#of2DEGchargeconcentration
regionnum=1mat=SiNy.min=0y.max=$y1
regionnum=2mat=AIGaNy.min=$y1y.max=$y2donors=1e16p=$xcpolarcalc.strainpolar.scale=-0.5
regionnum=3mat=GaNy.min=$y2y.max=$y4donors=1e15polarcalc.strainpolar.scale=-0.5
#
electname=sourcex.max=0y.min=$y1y.max=$y3
electname=drainx.min=6.0y.min=$y1y.max=$y3
electname=gatex.min=0.5x.max=0.9y.min=0.3y.max=$y1
electname=gatex.min=0.5x.max=$xly.min=0.3y.max=0.3
#
dopinggaussiancharacteristic=0.01conc=1e18n.typex.left=0.0\
x.right=0.05y.top=$y1y.bottom=$y3ratio.lateral=0.01direction=ydopinggaussiancharacteristic=0.01conc=1e18n.typex.left=$xd-0.05\
x.right=$xdy.top=$y1y.bottom=$y3ratio.lateral=0.01direction=y
ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff
7fTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfTf
#KMparameterset
ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff7fTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfTf
materialmaterial=GaNeg300=3.4align=0.8permitt=9.5\
mun=900mup=10vsatn=2e7nc300=1.07e18nv300=1.16e19\
real.index=2.67imag.index=0.001\
taun0=$lttaup0=$lt
materialmaterial=AIGaNaffinity=3.82eg300=3.96align=0.8permitt=9.5\
mun=600mup=10nc300=2.07e18nv300=1.16e19\
real.index=2.5imag.index=0.001\
taun0=$lttaup0=$lt
ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff
7fTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfT7ffTfT7ffTfTfT7ffTfTf
modelprintfermifldmobsrh
impactmaterial=GaNselban仁2.9e8an2=2.9e8bn仁3.4e7bn2=3.4e7\
即仁2.9e8即2=2.9e8bp仁3.4e7bp2=3.4e7
#
contactname=gatework=5.23
#人为引进光照以利于实现阻断状态下仿真收敛，这是仿真研究击穿的常用手段
beamnumber=1x.o=0y.o=$y4+0.1angle=270wavelength=0.3
#
outputcon.bandval.bandband.paramchargee.mobh.mobflowlinesqss
#IdVg特性求解
solve
logoutf=ganfetex02_O.log
solvevdrain=0.05
solvevstep=-0.2vfinal=-2name=gate
solvevstep=-0.1vfinal=-4name=gate
logoff
saveoutfile=ganfetex02_0.str
extractinitinfile="ganfetex02_0.log"
extractname="Vpinchoff"xintercept（maxslope（curve（v."gate",i."drain"
#IdVd击穿曲线
methodautonrgcarr.itlimit=10clim.dd=1e3clim.eb=1e3nblockit=25solveinit
#turnonopticalsourcetohelpinitiatebreakdown
##人为引进光照以利于实现阻断状态下仿真收敛
solveb1=$lightindex.check
#
solvensteps=10vfinal=$Vpinchoffname=gateb1=$light
logoutf=ganfetex02_$'index'.log
solvevstep=0.1vfinal=1name=drainb1=$light
solvevstep=1vfinal=10name=drainb1=$light
solvevstep=2vfinal=20name=drainb1=$light
solvevstep=5vfinal=1200name=drainb1=$lightcname=draincompl=0.5
#changetocurrentcontacttoresolvebreakdown
contactname=draincurrent
solve
solveimultistep=1.1ifinal=1name=drain
#
saveoutfile=ganfetex02_$'index'.str
#
extractinitinfile="ganfetex02_1$'index'.log"
extractname="a"slope（maxslope（curve（i."drain",v."drain"
extractname="b"xintercept（maxslope（curve（i."drain",v."drain"
extractname="Vdmax"max（curve（i."drain",v."drain"
extractname="ldmax"x.valfromcurve（i."drain”,v."drain"wherey.val=$Vdmax
extractname="Vd1"$Vdmax-20
extractname="Id1"y.valfromcurve（v."drain",i."drain"wherex.val=$Vd1
extractname="c"gradfromcurve（v."drain",i."drain"wherex.val=$Vdmax
extractname="d"$ldmax-$c*$Vdmax
extractname="Vbr"（$b-$d/（$c-（1/$a
extractname="ls"$b+$Vbr/$a
tonyplotganfetex02_1.strganfetex02_2.strganfetex02_3.strganfetex02_4.strganfetex02_5.strganfetex02_6.str-setganfetex02_1.set
tonyplot-overlayganfetex02_1.logganfetex02_2.logganfetex02_3」og
ganfetex02_4.logganfetex02_5.logganfetex02_6.log-setganfetex02_0.set
quit
图4-9为不同场板长度下半导体层中碰撞离化率的分布图。
正向I-V特性曲线
如图5所示,导通电压接近0.8Vo
图4场板长度L=1um的沟道中电子碰撞产生率模拟分布
0
图5场板长度L=1.25um的沟道中电子碰撞产生率模拟分布
気
J
图6场板长度L=1.5um的沟道中电子碰撞产生率模拟分布
图7场板长度L=1.75um的沟道中电子碰撞产生率模拟分布
—
图8场板长度L=2um的沟道中电子碰撞产生率模拟分布
FieWLengUi-20mcrw
图9场板长度L=2.25um的沟道中电子碰撞产生率模拟分布
图10-12是半导体中电场强度分布随场板长度的变化。
1
2
3
4
5
6
1000000
2000000
3000000
4000000
5000000
EIectric
FieId（V/m
L（um
图10不同场板长度的沟道中总电场分布
4
5
6
-5000000
-4000000
-3000000
-2000000
-1000000
EFieldX（V/m
L（um
图11不同场板长度的沟道中X电场分布
0123456
-4000000
-3000000
-2000000
-100000001000000
2000000
EFieldY（V/mL（um
图12不同场板长度的沟道中丫电场分布
图13是Id-Vd击穿曲线，可以清楚看到击穿电压从l=1um时的300V左右增大了l=2.25um时的800V以上。
图13不同长度的场板在关断情况下的输出I-V特性
通过对电场分布和碰撞离化率分布的分析知道，场板变长一方面会减弱漏端电场峰值，但另一方面也使发生碰撞离化的区域增大，所以这种构型的场板不是越长越好。
实验4短沟道MOS晶体管特性仿真
升级会员 