1、SilvacoTCADSilvacoTCAD器件器件仿真仿真3 3这一讲主要内容材料特性设置物理模型设置特性获取结果分析从例子hemtex01.in看整个流程*Silvaco学习2材料参数状态Material，设置材料参数材料参数和物理模型的选取有关，常用的参数及说明如下：*Silvaco学习3材料参数*Silvaco学习4材料参数*Silvaco学习5材料参数*Silvaco学习6能带参数能带结构示意图*Silvaco学习71221()bmsEcEcEgEgalign能带不连续参数align的说明要在结构文件中查看能带，需添加语句Output con.band val.band band.p
2、ara材料参数例子*Silvaco学习8material material=InGaAs align=0.36 eg300=0.75 nc300=2.1e17 nv300=7.7e18 copt=9.6e-11 material material=InP affinity=4.4 align=0.36 eg300=1.35 nc300=5.7e17 nv300=1.1e19 copt=1.2e-10 material region=1 taun0=5.0e-10 taup0=1.0e-9 vsatn=2.5e7 mun0=4000 mup0=200impact selb material=In
3、GaAs an2=5.15e7 ap2=9.69e7 bn2=1.95e6 bp2=2.27e6 impact selb material=InP an2=1e7 ap2=9.36e6 bn2=3.45e6 bp2=2.78e6Material taun0=1.e-9 taup0=1.e-9 f.conmun=hemtex01_interp.libmaterial align=0.6光学参数在光电特性仿真中材料的光学参数（折射率实部和虚部）尤为重要。两种方法可设置材料的光学参数1、C解释器编写参数文件2、添加、更改材料的折射率文件内的信息文件路径：X:sedatoolslibAtlas5.14
4、.0.Rcommon*.n为折射率实部，*.k为折射率虚部*Silvaco学习9C解释器可通过C解释器编辑函数来描述材料参数C解释器模板路径X:sedatoolslibAtlas5.14.0.Rcommontemplate例子hemtex01.in 中 “f.conmun=hemtex01_interp.lib”*Silvaco学习10物理模型models，计算时所采用的物理模型物理模型是指有物理意义的方程物理描述在手册的3.6节Physics Models部分*Silvaco学习1100300300tmunlntmuplpTmunTmup低场迁移率模型的物理方程及参数低场迁移率模型中可用户定
5、义的参数低场迁移率模型中可用户定义的参数状态参数默认值单位MobilityMun1000cm2/VsMobilityMup500cm2/VsMobilityTmun1.5MobilityTmup1.5物理模型推荐的模型MOSFETs类型：srh，cvt，bgnBJT，thyristors等：Klasrh，klaaug，kla，bgn击穿仿真：Impact，selb*Silvaco学习12Models conmob fldmob srh auger temp=300 print例句：Impact selbModel bgn fldmob srh*Silvaco学习13界面特性interface
6、y.max=0.1 qf=1e11interface x.min=4 x.max=4 y.min=0.25 y.max=0.1 qf=1e11 s.n=1e4 s.p=1e4Interface定义界面态电荷（密度cm-2），s.n和s.p分别为电子和空穴的表面复合速率计算方法Method，仿真时的数值计算方法主要的方法有newton，gummel，block*Silvaco学习14method gummel newton trap itlimit=20 maxtrap=6 报错信息：“trap times more than 4 times”指计算不收敛。需要注意的几个地方：Itlimit是迭
7、代的次数上限；maxtrap设置步长折回的次数；电流边界情形不能用Gummel特性获取加偏执是用solve状态先需要设置数据保存在日志文件，之后才可以用Tonyplot显示出来*Silvaco学习15Log outfile=*.logSolve init Solve vdrain=0.1Solve vgate=0.05 vstep=0.05 vfinal=1.0 name=gateSolve ibase=1e-6简单例句：特性获取*Silvaco学习16solve vdrain=0.1solve vdrain=0.2 solve vdrain=2.0solve vdrain=0.1solve
8、vgate=0.1 vstep=0.1 vfinal=3.0 name=gatesolve vcollector=2solve vbase=0.0 vstep=0.1 vfinal=2 name=basecontact name=base common=collectorsolve vbase=0.0 vstep=0.1 vfinal=2 name=baseI-V 特性：转移特性：GP 特性：特性获取*Silvaco学习17impact selbmethod trap climit=1e - 4 maxtrap=10#solve initsolve vbase=0.025solve vbase
9、=0.05solve vbase=0.2#contact name=base currentsolve ibase=3.e - 15#log outfile=gp.log mastersolve vcollector=0.2solve vstep=0.2 vfinal=5 name=collectorsolve vstep=0.5 vfinal=10 name=collector compl=5.e - 11 p=3CE击穿特性：特性获取*Silvaco学习18solve initsolve vbase=0.05 vstep=0.05 vfinal=0.8 name=basecontact n
10、ame=base current#solve ibase=2.e - 6save outf=bjt_ib_1.str mastersolve ibase=4.e - 6save outf=bjt_ib_2.str master#load inf=bjt_ib_1.str masterlog outf=bjt_ib_1.logsolve vcollector=0.0 vstep=0.25 vfinal=5.0 name=collector#load inf=bjt_ib_2.str masterlog outf=bjt_ib_2.logsolve vcollector=0.0 vstep=0.2
11、5 vfinal=5.0 name=collector电流控制性器件的输出特性：特性获取*Silvaco学习19solve initsolve vgate=1.1 outf=solve_tmp1solve vgate=2.2 outf=solve_tmp2#load infile=solve_tmp1log outf=mos_ids_1.logsolve name=drain vdrain=0 vstep=0.3 vfinal=3.3#load infile=solve_tmp2log outf=mos_ids_2.logsolve name=drain vdrain=0 vstep=0.3
12、vfinal=3.3电压控制型器件的输出特性：特性获取*Silvaco学习20solve vbase=0 vstep=0.05 vfinal=2.0 name=base ac freq=1e6solve vbase=0.7 ac freq=1e9 fstep=1e9 nfstep=10solve vbase=0.7 ac freq=1e6 fstep=2 mult.f nfstep=10solve prev ac freq=1e6 fstep=2 mult.f nfstep=10solve vgate=0 vstep=0.05 vfinal=1 name=gate ac freq=1e6 fs
13、tep=2 mult.f nfsteps=10某频率下的特性（CV）：电容随频率的变化（1G11GHz）：（1M1.024GHz）直流偏置和交流频率都改变：特性获取*Silvaco学习21solve vgate=1.0 ramptime=1e8 tstop=1e7 tstep=1e10瞬态特性：瞬态仿真参数的说明*Silvaco学习22特性获取log outf=ac.log s.param inport=gate outport=drain width=100solve ac.analysis direct frequency=1.e9 fstep=2.e9 nfsteps=20log out
14、f=ac.log s.param inport=gate outport=drain in2port=source out2port=source width=100 rin=100solve ac.analysis direct frequency=1.e9 fstep=2.e9 nfsteps=20S参数仿真（电流增益随频率变化的特性）：*Silvaco学习23特性获取models srh fldmob bz=1.0model temperature=300solve init#probe name=hall1 x=0.0 y=0.5 potentialprobe name=hall2 x
15、=1.0 y=0.5 potentialprobe name=reference x=0.5 y=0.5 potential#log outf=hall_effect.logsolve name=anode vanode=0.0 vstep=0.05 vfinal=1.0霍尔效应仿真：*Silvaco学习24特性获取beam num=1 x.orign=5 y.orign=2 angle=90 wavelenght=.8beam num=1 x.orign=5 y.orign=2 angle=90 wavel.start=.5 wavel.end=1.7 wave.num=13beam num
16、=2 x.orign=1 y.orign=1 angle=90 wavelength=1.5 back.refl front.refl reflect=5 min.power=0.01beam num=3 x.orign=2 y.orign=0.5 angle=90 wavelength=0.9 rays=101 gaussian mean=0 xsigma=0.25光电效应仿真时光线的设置：波长、入射点位置和角度：一定范围的波长：前后背面的反射：光强为高斯分布：*Silvaco学习25特性获取beam num=1 x.origin=2.5 y.origin=1.0 angle=90.0 wa
17、velength=.4#solve initsolve vcathode=0.05 vstep=0.05 vfinal=1.0 name=cathode#log outf=spectral_Response.log mastersolve prev b1=2 lambda=0.5solve prev b1=2 lambda=0.6solve prev b1=2 lambda=0.7solve prev b1=2 lambda=0.8solve prev b1=2 lambda=0.9solve prev b1=2 lambda=1.0solve prev b1=2 lambda=1.2光电效应
18、仿真：*Silvaco学习26特性获取models bgn srh conmob fldmob method newton# set up initial DC biassolve vcathode=0.1 vstep=0.1 vfinal=3.0 name=cathode # set up the track of the radiation for SEUsingleeventupset entry=2.0,1.5,0. exit=1.,1.,4 radius=0.05 density=1e18method halfimpl dt.min=3.e-12log outf=seu_0.logs
19、olve tfinal=1.e-9 tstep=3.e-12SEU仿真：三维器件仿真三维器件仿真流程和语法可参照二维仿真位置坐标是三维数据物理模型少数会有差别*Silvaco学习27*Silvaco学习28仿真结果分析仿真结果的表现形式1、包含有电学、光学、热学信息等的结构文件（可由Tonyplot 导出数据）2、仿真得到的日志文件（可由Tonyplot 导出数据）3、抽取的结果（*.dat 或是在实时输出窗口显示）*Silvaco学习29从结构文件中导出数据Format:Tonyplot User Data*Silvaco学习30从日志文件导出数据Format:Comma Separated Values 谢谢！*31Silvaco学习器件仿真，Silvaco仿真至此全部结束 交流时间*32Silvaco学习
