半导体专业实验补充silvaco器件仿真Word文件下载.docx

半导体专业实验补充silvaco器件仿真Word文件下载.docx

《半导体专业实验补充silvaco器件仿真Word文件下载.docx》由会员分享，可在线阅读，更多相关《半导体专业实验补充silvaco器件仿真Word文件下载.docx（30页珍藏版）》请在冰豆网上搜索。

＃启动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.log

＃退出

quit

图2为普通耐压层功率二极管的仿真结构.正向I-V特性曲线如图3所示，导通电压接近0.8V。

图2普通耐压层功率二极管的仿真结构

图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.5×

105V•cm—1,在耐压层中线性减小到80000V•cm—1。

图4普通耐压层功率二极管的反向I-V特性曲线

图5普通耐压层功率二极管击穿时的电场分布

导通的二极管突加反向电压,需要经过一段时间才能恢复反向阻断能力。

电路图如图6所示.设t=0前电路已处于稳态,Id=If0。

t=0时，开关K闭合，二极管从导通向截止过渡。

在一段时间内，电流Id以di0/dt=-Ur/L的速率下降.在一段时间内电流Id会变成负值再逐渐恢复到零。

仿真时先对器件施加一个1V的正向偏压，然后迅速改变电压给它施加一个反向电压增大到2V.

solvevanode=1

logoutf=cj2_1。

log

solvevcathode=2。

0ramptime=2.0e-8tstop=5。

0e—7tstep=1.0e-10

反向恢复特性仿真时，也可以采用如图7的基本电路,其基本原理为：

在初始时刻，电阻R1的值很小，电阻R2的值很大，例如可设R1为1×

10—3

，R2为1×

106

；

电感L1可设为3nH;

电压源及电流源也分别给定一个初始定值v1，i1；

那么由于R2远大于R1，则根据KCL可知，电流i1主要经过R1支路，即i1的绝大部分电流稳定的流过二极管，二极管正向导通，而R2支路几乎断路，没有电路流过。

然后,在短暂的时间内，使电阻R2的阻值骤降。

此时,电阻器R2作为一个阻源，其阻值在极短的时间间隔内以指数形式从1×

下降到1×

10-3

。

这一过程本质上是使与其并联的连在二极管阳极的电流源i1短路，这样电流i1几乎全部从R2支路流过，而二极管支路就没有i1的分流，此刻电压源v1开始起作用,二极管两端就被施加了反偏电压，由于这些过程都在很短的时间内完成，因而能够很好的实现二极管反向恢复特性的模拟。

反向恢复特性仿真图如图8所示，PN结功率二极管的反向恢复时间约为50ns。

图6反向恢复特性测试原理电路图

图7二极管反向恢复特性模拟电路图

图8器件反向恢复特性曲线

实验3PN结终端技术仿真

由于PN结在表面的曲率效应,使表面的最大电场常大于体内的最大电场，器件的表面易击穿，采用终端技术可使表面最大电场减小,提高表面击穿电压。

场限环和场板是功率器件中常用的两种终端技术。

场限环技术是目前功率器件中被大量使用的一种终端技术。

其基本原理是在主结表面和衬底之间加反偏电压后，主结的PN结在反向偏压下形成耗尽层，并随着反向偏置电压的增加而增加.当偏置电压增加到一定值是，主结的耗尽层达到环上,如图1所示，这样就会使得有一部分电压有场环分担，将主结的电场的值限制在临界击穿电压以内,这将显著的减小主结耗尽区的曲率，从而增加击穿电压。

图1场限环

场板结构在功率器件中被广泛应用.场板结构与普通PN结的区别在于场板结构中PN区引线电极横向延伸到PN区外适当的距离。

而普通PN结的P区引线电极的横向宽度一般不超过P扩散区的横向尺寸。

PN结反向工作时,P区相对于N型衬底加负电位。

如果场板下边的二氧化硅层足够厚，则这个电场将半导体表面的载流子排斥到体内，使之表面呈现出载流子的耗尽状态，如图2所示，就使得在同样电压作用下，表面耗尽层展宽，电场减小，击穿电压得到提高。

（1）场限环特性仿真

场限环：

击穿电压200V，设计3个环，环的宽度依次为6、5、5、5μm，间距为4、5、6μm，外延层浓度为1×

1015cm-3，观察表面电场。

（2）场板特性仿真

场板：

氧化层厚度1μm，结深1μm，场板长度分别为0μm、2μm、4μm、6μm、8μm、10μm，外延层浓度为1×

图2场板

3、场板的应用实例：

场板对大功率GaNHEMT击穿电压的影响

（1）内容

（a）GaNHEMT的工作机理、击穿特性刻画以及对场板结构的GaNHEMT击穿特性的进行仿真分析。

（b）结构和参数:

场板结构的GaNHEMT的结构尺寸及掺杂浓度如图3所示。

图3场板结构的大功率GaNHEMT

（2）要求

（a）掌握定义一个完整半导体器件结构的步骤，并能对其电性能进行仿真研究。

（b）理解场板技术对器件击穿电压提高的作用原理并能结合仿真结果给出初步分析。

（3）实验过程

＃启动internal，定义结构参数

#场板长度从1um增大到2。

25um，步长为0.25um,通过改变l取值来改变场板长度

setl=1.0

＃drain—gatedistance

setLdg=5.1

＃fieldplatethickness

sett=1.77355

#AlGaNcompositionfraction

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

#启动二维器件仿真器

meshwidth=1000

#网格结构

x。

ml=0。

0s=0.1

05s=0.05

x.ml=0。

5s=0.05

9s=0。

025

ml=（0。

9+$xl）/2s=0.05

x.ml=$xls=0.025

ml=（＄xl+$xd）/2s=0。

25

ml=＄xd-0.05s=0.05

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.0100

y.ml=$y4s=0。

0500

＃devicestructure

#POLAR.SCALEischosentomatchcalibratedvalues

#of2DEGchargeconcentration

regionnum=1mat=SiNy。

min=0y。

max=$y1

regionnum=2mat=AlGaNy.min=$y1y。

max=$y2donors=1e16x。

comp=＄xcpolarcalc。

strainpolar.scale=-0.5

regionnum=3mat=GaNy.min=＄y2y.max=＄y4donors=1e15polarcalc.strainpolar.scale=-0。

#

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=y

left=＄xd—0。

05\

x。

right=＄xdy。

top=＄y1y.bottom=$y3ratio.lateral=0。

＃####＃#＃#＃＃＃##＃##＃#＃＃###＃##＃＃#＃＃＃#####＃##＃#＃#＃#＃#＃####＃######＃##＃＃#

#KMparameterset

#＃#＃#＃###＃＃＃＃＃##＃#＃＃＃＃＃##＃#＃##＃#＃#＃＃＃＃＃##＃＃＃＃###＃#＃#＃＃＃#＃######＃＃#＃

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=AlGaNaffinity=3。

82eg300=3。

96align=0.8permitt=9.5\

mun=600mup=10nc300=2。

real。

index=2.5imag.index=0.001\

＃＃＃##＃#＃＃##＃＃＃＃#＃＃＃＃#＃＃####＃＃＃＃＃#＃##＃##＃#＃＃##＃＃＃#＃＃＃#＃＃##＃#####＃＃＃＃

modelprintfermifldmobsrh

impactmaterial=GaNselban1=2。

9e8an2=2。

9e8bn1=3.4e7bn2=3.4e7\

ap1=2。

9e8ap2=2。

9e8bp1=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_0。

solvevdrain=0。

05

solvevstep=-0。

2vfinal=—2name=gate

solvevstep=-0.1vfinal=-4name=gate

logoff

saveoutfile=ganfetex02_0。

extractinitinfile="

ganfetex02_0.log"

extractname=”Vpinchoff"

xintercept（maxslope（curve（v.”gate”，i。

”drain”）））

#IdVd击穿曲线

methodautonrgcarr.itlimit=10clim。

dd=1e3clim。

eb=1e3nblockit=25

solveinit

＃turnonopticalsourcetohelpinitiatebreakdown

＃#人为引进光照以利于实现阻断状态下仿真收敛

solveb1=＄lightindex.check

solvensteps=10vfinal=$Vpinchoffname=gateb1=$light

logoutf=ganfetex02_＄’index'

solvevstep=0.1vfinal=1name=drainb1=＄light

solvevstep=1vfinal=10name=drainb1=$light

solvevstep=2vfinal=20name=drainb1=$light

solvevstep=5vfinal=1200name=drainb1=＄lightcname=draincompl=0。

＃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"

）））

b"

xintercept（maxslope（curve（i."

drain”,v。

drain”）））

Vdmax”max（curve（i。

drain”））

extractname=”Idmax”x。

valfromcurve（i."

，v。

）wherey.val=$Vdmax

extractname=”Vd1"

$Vdmax—20

extractname=”Id1"

y。

valfromcurve（v。

，i.”drain”）wherex。

val=$Vd1

c"

gradfromcurve（v.”drain"

，i.”drain”）wherex.val=$Vdmax

d”$Idmax—$c＊$Vdmax

Vbr”（$b-＄d）/（$c-（1/＄a））

Is"

＄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.logganfetex02_4。

logganfetex02_5.logganfetex02_6.log—setganfetex02_0。

set

图4—9为不同场板长度下半导体层中碰撞离化率的分布图。

正向I—V特性曲线如图5所示，导通电压接近0。

8V。

图4场板长度L=1um的沟道中电子碰撞产生率模拟分布

图5场板长度L=1.25um的沟道中电子碰撞产生率模拟分布

图6场板长度L=1.5um的沟道中电子碰撞产生率模拟分布

图7场板长度L=1.75um的沟道中电子碰撞产生率模拟分布

图8场板长度L=2um的沟道中电子碰撞产生率模拟分布

图9场板长度L=2.25um的沟道中电子碰撞产生率模拟分布

图10—12是半导体中电场强度分布随场板长度的变化.

图10不同场板长度的沟道中总电场分布

图11不同场板长度的沟道中X电场分布

图12不同场板长度的沟道中Y电场分布

图13是Id—Vd击穿曲线，可以清楚看到击穿电压从l=1um时的300V左右增大了l=2。

25um时的800V以上.

图13不同长度的场板在关断情况下的输出I-V特性

通过对电场分布和碰撞离化率分布的分析知道，场板变长一方面会减弱漏端电场峰值，但另一方面也使发生碰撞离化的区域增大，所以这种构型的场板不是越长越好。

实验4短沟道MOS晶体管特性仿真

（1）短沟道LDD—MOS晶体管结构定义.

（2）转移特性、输出特性。

（3）结构和参数：

器件结构下图所示，宽度1。

2μm，衬底为P型、厚度0.8μm、浓度1×

1014cm-3、晶向<

100>

，栅氧化层厚度13nm，栅为n+掺杂多晶硅.

（2）改变表面浓度，改变栅氧化层厚度,观察阈值电压变化。

#启动Athena

＃器件结构网格划分;

linexloc=0.2spac=0.006

linexloc=0.4spac=0。

006

linexloc=0.6spac=0.01

0spac=0.002

2spac=0.005

5spac=0.05

lineyloc=0.8spac=0。

15

（建议定义左边一半）

#初始化；

＃栅氧化，干氧11分钟，温度950。

diffustime=11temp=950dryo2press=1。

00hcl.pc=3

提取栅氧化层厚度，extractname=”Gateoxide"

thicknessmaterial=”Sio-2"

mat.occno=1x.val=0。

#阈值电压调整；

implantborondose=9.5e11energy=10crystal

提取表面浓度

#淀积多晶硅；

depopolythick=0。

2divi=10

#定义多晶硅栅

etchpolyleftp1.x=0.35

＃多晶硅氧化，湿氧,900度，3分钟；

methodfermicompress

diffusetime=3temp=900weto2press=1。

＃多晶硅掺杂

implantphosphordose=3.0e13energy=20crystal

＃侧墙的形成

淀积氧化层:

depooxidethick=0。

12divisions=10

干法刻蚀：

etchoxidedrythick=0.12

#源漏砷注入，快速退火

implantarsenicdose=5.0e15energy=50crystal

methodfermi

diffusetime=1temp=900nitropress=1.0

#金属化

etchoxideleftp1.x=0。

2

depositaluminthick=0。

03divi=2

etchaluminrightp1。

x=0.18

＃提取器件参数：

结深，源漏方块电阻,侧墙下的方块电阻，阈值电压

＃extractfinalS/DXj

extractname=”nxj"

xjsiliconmat。

occno=1x。

val=0。

1junc.occno=1

＃extracttheN++regionssheetresistance

extractname=”n++sheetrho”sheet。

resmaterial=”Silicon"

mat。

occno=1x.val=0.05region。

occno=1

＃extractthesheetrhounderthespacer，oftheLDDregion

extractname=”lddsheetrho"

sheet。

resmaterial=”Silicon”\

mat.occno=1x.val=0.3region.occ