一种肖特基二极管的仿真模型的计算方法及系统
阅读说明：本技术 一种肖特基二极管的仿真模型的计算方法及系统 (Calculation method and system of simulation model of Schottky diode ) 是由 祁路伟 张德海 孟进 于 2019-10-30 设计创作，主要内容包括：本发明公开了一种肖特基二极管的仿真模型的生成方法及系统,所述方法包括：将肖特基二极管的轻掺杂的外延层上方设置SiO<Sub>2</Sub>介质层,介质层上方定义为肖特基接触；根据MIS结构的能带图计算介质层的厚度和介质层与外延层界面处的陷阱电荷密度；利用Sentaurus TCAD工具绘制出肖特基二极管的结构图形；为每一个区域指定材料以及掺杂浓度；添加接触作为器件的电极,添加肖特基接触作为器件的阳极,欧姆接触为器件的阴极；将肖特基二极管的结构图形进行剖分和网格化。利用该生成方法得到的模型的仿真结果与实际结果吻合,为后续建立准确的高频模型打下坚实基础。(The invention discloses a method and a system for generating a simulation model of a Schottky diode, wherein the method comprises the following steps: arranging SiO above the lightly doped epitaxial layer of the Schottky diode 2 A Schottky contact is defined above the dielectric layer; calculating the thickness of the dielectric layer and the trap charge density at the interface of the dielectric layer and the epitaxial layer according to the energy band diagram of the MIS structure; drawing a structural graph of the Schottky diode by using a Sentaurus TCAD tool; specifying a material and a doping concentration for each region; adding a contact as an electrode of the device, adding a Schottky contact as an anode of the device, and using an ohmic contact as a cathode of the device; and (4) subdividing and meshing the structural graph of the Schottky diode. The simulation result of the model obtained by the generation method is consistent with the actual result, and a solid foundation is laid for the subsequent establishment of an accurate high-frequency model.)
技术领域
本发明涉及太赫兹领域，具体涉及一种肖特基二极管的仿真模型的计算方法及系统。
背景技术
太赫兹(THz)波是频率介于0.3THz～3THz范围内的电磁波，其波长短，是很好的宽带信息载体，在高分辨率光谱学、成像系统、通信领域有很大的应用前景；太赫兹波能够穿透生物组织并不会对生物细胞产生有害电离，特别适合生物医学成像，并有取代X射线成像的趋势；太赫兹波的光子能量能够与物质基本单元之间的能量产生共振，在物质科学研究领域具有很广阔的应用前景。随着太赫兹波的广泛应用，对高功率、高频率倍频电路的需求越来越大。
砷化镓(GaAs)肖特基二极管在毫米波和太赫兹波频段的应用中占据重要的地位。由于目前肖特基二极管的器件模型与倍频电路相关设计的仿真模型不同步，在高频段的设计中，倍频电路的实际测试结果与仿真结果差距较大，为高频段、高功率太赫兹倍频单片电路的研究造成了困难。因此，需要建立准确的砷化镓(GaAs)肖特基二极管器件模型，从而为其高频段、高功率的设计与应用提供参考和依据。
由于实际的肖特基二极管器件在肖特基接触区域不存在介质，即肖特基接触为半导体和金属Ti的直接接触。但是由于半导体的边缘存在晶格损伤以及在半导体和金属的接触处存在一定浓度的界面态，该损伤和界面态会引起一定的泄漏电流，使实际的器件理想因子大于1。采用实际的肖特基二极管结构(肖特基金属/肖特基接触/外延层/缓冲层)仿真器件的电学特性时，仿真的结果是理想的状态：理想因子n＝1，不能预测出实际器件电流的构成。由于半导体的边缘存在化学键的断裂、金属和半导体接触时存在晶格损伤导致一定密度的界面电荷，该损伤和缺陷电荷会造成一定的泄漏电流即隧穿电流，使实际的器件理想因子大于1。
发明内容
本发明的目的在于克服上述技术缺陷，将损伤和缺陷电荷因素等效为在器件的肖特基金属与半导体接触之间存在一定厚度的介质层，并为该区域赋予一定密度的界面电荷，从而准确预测实际的电流分布即隧穿电流所占热电子发射电流的比值，并在SentaurusTCAD中得以验证。
为了实现上述目的，本发明提出了一种肖特基二极管的仿真模型的生成方法，所述方法包括：
将肖特基二极管的轻掺杂的外延层上方设置SiO2介质层，介质层上方定义为肖特基接触；
根据MIS结构的能带图计算介质层的厚度，以及介质层与外延层界面处的陷阱电荷密度；
介质层与外延层界面处的陷阱电荷密度；
利用Sentaurus TCAD工具绘制出肖特基二极管的结构图形；
为每一个区域指定材料以及掺杂浓度；
添加接触作为器件的电极，添加肖特基接触作为器件的阳极，欧姆接触为器件的阴极；
将肖特基二极管的结构图形进行剖分和网格化。
作为上述方法的一种改进，所述介质层的厚度、介质层与外延层界面处的陷阱电荷密度的计算步骤包括：
根据MIS结构的能带图得到如下关系：
qΦm＝qΔ+qλ+qΦB0 (1)
qΦB0＝qψS(V)+qVn-qV (2)
其中，Φm为金属功函数，ΦB0为势垒高度，Δ为界面层上的电势，λ为半导体的电子亲合能，ψS(V)为内建电势，Vn为300K时热电压，k为玻尔兹曼常数，T为温度；q为单位电荷，V为偏压；
由(1)和(2)得到：
Δ＝Φm-λ-ψS(V)-Vn+V (3)
则介质层与外延层界面处的陷阱电荷密度Qit为：
Qit为界面陷阱电荷密度，Dit(0)为界面陷阱密度，Eg为禁带宽度，
结合(3)式得到：
热平衡状态下，半导体耗尽层中形成的空间电荷密度Qsc(V)为：
其中，ND为轻掺杂层掺杂浓度，WD为耗尽层厚度，εs为半导体的介电常数；
固定电荷Qf为：
Qf＝qNf
其中，Qf为固定电荷密度，单位为C/cm-2，Nf为固定电荷密度，单位为/cm-2；
界面层上的电势Δ通过金属和半导体表面电荷运用高斯定理得到：
其中，Qm为金属表面电荷密度，δ为界面层厚度，εi为界面层的介电常数；
代入上式公式得到：
将其代入(3)中得到：
则理想因子n(V)为：
结合(*1)和(*2)，求得介质层与外延层界面处的陷阱电荷密度Qit＝Qit0和介质层的厚度δ＝δ0。
作为上述方法的一种改进，所述肖特基二极管的材料为GaAs、氮化镓或碳化硅。
本发明还提出了一种肖特基二极管的仿真模型的生成系统，所述系统包括：
计算模块，用于将肖特基二极管的轻掺杂的外延层上方设置SiO2介质层，介质层上方定义为肖特基接触；根据MIS结构的能带图计算介质层的厚度和介质层与外延层界面处的陷阱电荷密度；
结构图形绘制模块，用于利用Sentaurus TCAD工具绘制出肖特基二极管的结构图形；
材料设置模块，用于为每一个区域指定材料以及掺杂浓度；
电极设置模块，用于添加接触作为器件的电极，添加肖特基接触作为器件的阳极，欧姆接触为器件的阴极；
结构剖分模块，用于将肖特基二极管的结构图形进行剖分和网格化。
作为上述系统的一种改进，所述介质层的厚度、介质层与外延层界面处的陷阱电荷密度的计算步骤包括：
根据MIS结构的能带图得到如下关系：
qΦm＝qΔ+qλ+qΦB0 (1)
qΦB0＝qψS(V)+qVn-qV (2)
其中，Φm为金属功函数，ΦB0为势垒高度，Δ为界面层上的电势，λ为半导体的电子亲合能，ψS(V)为内建电势，Vn为300K时热电压，
由(1)和(2)得到：
Δ＝Φm-λ-ψS(V)-Vn+V (3)
则介质层与外延层界面处的陷阱电荷密度Qit为：
Qit为界面陷阱电荷密度，Dit(0)为界面陷阱密度，Eg为禁带宽度，
结合(3)式得到：
热平衡状态下，半导体耗尽层中形成的空间电荷密度Qsc(V)为：
其中，ND为轻掺杂层掺杂浓度，WD为耗尽层厚度，εs为半导体的介电常数；
固定电荷Qf为：
Qf＝qNf
其中，Qf为固定电荷密度，单位为C/cm-2，Nf为固定电荷密度，单位为/cm-2；
界面层上的电势Δ通过金属和半导体表面电荷运用高斯定理得到：
其中，Qm为金属表面电荷密度，δ为界面层厚度，εi为界面层的介电常数；
代入上式公式得到：
将其代入(3)中得到：
则理想因子n(V)为：
结合(*1)和(*2)，求得介质层与外延层界面处的陷阱电荷密度Qit＝Qit0和介质层的厚度δ＝δ0。
作为上述系统的一种改进，所述肖特基二极管的材料为GaAs、氮化镓或碳化硅。
本发明的优势在于：
1、本发明的方法将损伤和缺陷电荷因素等效为在器件的肖特基金属与半导体接触之间存在一定厚度的介质层，并为该区域赋予一定密度的界面电荷，从而准确预测实际的电流分布即隧穿电流所占热电子发射电流的比值；
2、本发明的方法基于Sentaurus TCAD器件仿真软件，考虑到肖特基接触界面处的陷阱电荷和界面态，提出一种肖特基二极管器件仿真模型的生成方法，利用该生成方法得到的模型的仿真结果与实际结果吻合，为后续建立准确的高频模型打下坚实基础。
附图说明
图1为用SDE工具中的三维结构图；
图2为本发明的肖特基二极管MIS结构的纵向解剖图；
图3为具有一定距离的界面层的金属和n型半导体接触的详细能带图；
图4为正向IV特性对比图；
图5为反偏特性对比图。
具体实施方式
下面结合附图对本发明的技术方案进行详细的说明。
本发明基于Sentaurus TCAD器件仿真软件，考虑到肖特基接触界面处的陷阱电荷和界面态，提出一种基于MIS结构的肖特基二极管器件模型的计算方法，利用该计算方法得到的模型的仿真结果与实际结果吻合，为后续建立准确的高频模型打下坚实基础。
在Sentaurus TCAD中完成器件的电学特性仿真，首先利用SDE工具画出器件的结构图(如图1所示)：
1、绘制出该器件的结构图形；
2、为每一个区域指定材料以及掺杂浓度；
3、添加接触作为器件的电极，本实施例中添加肖特基接触作为器件的阳极，欧姆接触(电阻接触)作为器件的阴极；
4、将该结构剖分网格以便于利用模型求解。
本实施例中，GaAs衬底厚度为50um；重掺杂的缓冲层厚度为1.5um，掺杂浓度为5E18/cm3，掺杂类型为n型；轻掺杂的外延层厚度为300nm，掺杂浓度为2E17/cm3，掺杂类型为n型；外延层上方是介质层，本发明中采用SiO2，厚度为δ0，介质层与外延层界面处界面处定义密度为Qit的陷阱电荷，类型为受主；介质层上方定义肖特基接触。如图2所示。
综上所述，为了分析器件的物理机制包括理想因子的决定因素以及漏电流的来源并仿真出与实际器件完全一致的电学特性，将肖特基接触处的界面态和晶格损伤等效成为一定厚度的介质层，该介质层的厚度和介质层与外延层界面处的陷阱电荷密度可以计算得到。
图3为具有一定距离的界面层的金属和n型半导体接触的详细能带图，在该结构处施加一定的偏压V根据MIS结构的能带图得到如下关系：
根据MIS结构的能带图得到如下关系：
qΦm＝qΔ+qλ+qΦB0 (1)
qΦB0＝qψS(V)+qVn-qV (2)
其中，Φm为金属功函数，ΦB0为势垒高度，Δ为界面层上的电势，λ为半导体的电子亲合能，ψS(V)为内建电势，Vn为300K时热电压，
由(1)和(2)得到：
Δ＝Φm-λ-ψS(V)-Vn+V (3)
则介质层与外延层界面处的陷阱电荷密度Qit为：
Qit为界面陷阱电荷密度，Dit(0)为界面陷阱密度，Eg为禁带宽度，
结合(3)式得到：
热平衡状态下，半导体耗尽层中形成的空间电荷密度Qsc(V)为：
其中，ND为轻掺杂层掺杂浓度，WD为耗尽层厚度，εs为半导体的介电常数；
固定电荷Qf为：
Qf＝qNf
其中，Qf为固定电荷密度，单位为C/cm-2，Nf为固定电荷密度，单位为/cm-2；
界面层上的电势Δ通过金属和半导体表面电荷运用高斯定理得到：
其中，Qm为金属表面电荷密度，δ为界面层厚度，εi为界面层的介电常数；
代入上式公式得到：
将其代入(3)中得到：
则理想因子n(V)为：
结合(*1)和(*2)，求得介质层与外延层界面处的陷阱电荷密度Qit＝Qit0和介质层的厚度δ＝δ0。
将计算得到得结果代入仿真计算中，得到的仿真曲线。与实际测试曲线对比，如图4和图5所示。
本发明结合实际的肖特基与半导体接触界面处的陷阱电荷和界面态，利用Sentaurus TCAD工具提出一种基于MIS结构的GaAs肖特基二极管器件仿真模型的生成方法，利用该方法得到的仿真结果与实际结果吻合，为后续建立准确的高频模型提供理论依据。该方法同样适用于其他材料的肖特基二极管，包括氮化镓、碳化硅。
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
- 上一篇：一种医用注射器针头装配设备
- 下一篇：指令处理方法、装置及相关产品