二维数值仿真软件MEDICI概述.pdf

《二维数值仿真软件MEDICI概述.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、二维数值仿真软件MEDICI概述1MEDICI功能简介[28]Medici是一款二维数值仿真软件，是由先驱公司研发的。它对势能场和载子的二维分布建模，通过解泊松方程和电子、空穴的电流连续性等方程来获取特定偏置下的电学特性。用该软件可以对双极型、MOS型等半导体器件进行模拟，这个程序通过解二极管和双极型三极管以及和双载流子有关的电流效应（诸如闩锁效应）的电流连续性方程和泊松方程来分析器件。Medici也能分析单载流子起主要作用的器件，例如：MOSFET，JFET，MESFET。另外，MEDICI还可以被用来分析器件在瞬态情况下的变化。在亚微米器件模

2、拟中，MEDICI通过联解电子和空穴的能量平衡和其他的器件方程，可以对深亚微米的器件进行模拟。像热载流子和速度过冲等效应在MEDICI中都已经考虑了，并能够对它们的影响进行分析。2MEDICI的基本物理描述2.1基本方程MEDICI的主要功能就是分别对静电势Ψ、电子浓度n和空穴浓度p求解三大类自连续的微分方程，包括泊松方程、连续性方程和波尔兹曼输运理论（即电流密度方程）。1、泊松方程：半导体器件的电学行为由泊松方程控制。2+−ερ∇Ψ=−q(pnNN−+−)−（1-1）DAs2、连续性方程：电子和空穴的连续性方程也控制着电学行为。∂n1→→=∇⋅

3、J−(UG−)=F(,,)Ψnpnnnn∂tq（1-2）∂−p1→→=∇⋅J−(UG−)=F(,,)Ψnppppp∂tq（1-3）→3、波尔兹曼输运理论：在波尔兹曼输运理论中，公式（1-2）中的J和公式n→（1-3）中的J可以被描述成关于载流子浓度和电子及空穴的准费米势的两个方p程。J=−∇qnµφ（1-4）nnnJ=−∇qnµφ（1-5）ppp同样，J和J可以被写成关于Ψ、n和p的方程式，由漂移电流和扩散电np流两部分组成。J=qµEnqD+∇n（1-6）nnnnJ=qµEpqD−

4、∇p（1-7）pppp其中，µ和µ分别是电子和空穴的迁移率，D和D是电子和空穴的扩散npnp系数，假设服从波尔兹曼载流子统计规律并忽略能带变窄效应带来的影响。2.2三种基本的复合方式现在MEDICI支持求解以下三种体内复合，分别是间接复合（也叫肖特基-里德-霍尔复合，即SRH复合），俄歇复合和直接复合（也称为带到带复合或者光学复合），公式（1-2）和（1-3）中的Un和Up代表电子和空穴的复合率，其为正值时，代表载流子复合，为负值时则代表载流子产生。则：UUUU===++UU（1-8）npSRHAugerdir其中：2pnn−iU=SRHETRA

5、P−ETRAPττ[nn+exp()]++[pnexp()]pinikTkT（1-9）2U=CDIRECTnpn.()−（1-10）diri2222U=AUGNpn()()−+nnAUGPnp−pn（1-11）Augerii上述公式中，ni代表本征载流子浓度，τn和τp代表电子和空穴寿命，并可能和浓度相关，参数ETRAP代表陷阱能级Et和本征费米能级Ei之差，AUGN和AUGP为与材料相关的特定常数。MATERIAL语句可以修改ETRAP，AUGN，AUGP和C.DIRECT的值。2.3其它物理描述除以上的基本方程及三种基本复合方式之外，MEDI

6、CI还支持以下物理描述，包括：与浓度相关的寿命、与晶格温度相关的寿命、复合导致的隧穿、肖特基势垒隧穿、波尔兹曼统计、禁带与有效态密度、重掺杂导致的禁带变窄、应力导致的禁带变化、费米-狄拉克统计、杂质不完全离化、低温模拟、MOSFET反型层中的量子效应、修正的局部密度近似量子模型等，相关内容可参考MEDICI手册。2.4MEDICI编程中的关键MEDICI编程（以跑反向耐压为例）主要包括了定义变量、网格描述、结构及材料描述、杂质分布描述、求解所用的模型、方法及特定技巧、求解条件的设定等部分组成。其中，结构及材料描述、杂质分布描述、求解条件的设定等部

7、分是由器件本身决定的，变量定义是个人喜好的问题，这几个部分的编程是相对比较简单的。网格描述是编程的重点，这决定了程序的运行效率。在求解的模型、方法、技巧乃至求解条件都一样的前提下，一套好的网格可以在不影响结果精度的前提下大大缩短程序运行的时间，并且提高程序的收敛性。求解所用的模型、方法及技巧的选择是编程的难点，这直接影响程序的可靠性，在不同的情况下，对模型的选择不仅必然影响结果的可靠性，而且会严重影响到程序的收敛性。以下我们将对网格描述及模型取舍作进一步的讨论。2.4.1网格描述三个基本注意点αN首先，格点数目N对仿真时间有直接影响，求解一次的算

8、法步骤正比于p，pα的取值通常在1.5~2之间变化。因此，格点增加，运算时间大幅度增加。其次，器件的各个区域的电学行为不一致，所以需要在