当前位置：首页 > news >正文

从论文到仿真：手把手复现GaN二极管声子辅助隧穿效应的Silvaco实现

news 2026/5/12 7:48:50

从理论到实践：GaN二极管声子辅助隧穿效应的Silvaco仿真全解析

在半导体器件研究中，氮化镓(GaN)肖特基二极管的反向漏电流机制一直是学术界和工业界关注的焦点。传统热电子发射理论往往无法完全解释实验观察到的温度依赖性现象，而声子辅助隧穿模型为这一难题提供了新的理论视角。本文将带您深入理解Pipinys理论模型的核心思想，并逐步演示如何通过Silvaco Atlas工具实现从理论公式到仿真验证的完整闭环。

1. 理论基础与模型构建

声子辅助隧穿理论由Pipinys和Lapeika在2006年首次提出，其核心观点是：在反向偏置条件下，电子并非单纯通过热激发跨越势垒，而是借助声子相互作用实现隧穿过程。这一机制特别适用于解释宽禁带半导体如GaN中的异常漏电流现象。

1.1 关键物理参数解析

模型中包含四个核心参数，它们共同决定了隧穿概率：

参数符号	物理意义	典型值范围	单位
PIP.NT	界面陷阱密度	1e13-1e14	cm^-3
PIP.ET	陷阱能级深度	0.7-0.9	eV
PIP.ACC	电子-声子耦合系数	1.5-2.5	无量纲
PIP.OMEGA	特征声子能量	0.06-0.08	eV

在Silvaco实现中，这些参数通过MATERIAL语句进行定义：

material material=GaN align=0.68 eab=0.16 \ taun0=1e-9 taup0=2e-8 nsrhn=4e18 nsrhp=4e18 \ d.tunnel=4e-6 me.sbt=0.20 mh.sbt=0.20 \ augn=1e-31 augp=1e-31 arichn=26 \ pip.nt=1.5e13 pip.et=0.8 pip.acc=2.0 pip.omega=0.07

1.2 模型激活与边界条件设置

要使声子辅助隧穿模型生效，必须在接触定义中显式启用PIPINYS参数：

contact name=anode workfun=5.18 pipinys surf.rec me.tunnel=0.222

其中：

workfun=5.18设定金属功函数为5.18eV
me.tunnel=0.222表示有效隧穿质量比为0.222
surf.rec启用表面复合速度模型

2. Silvaco仿真环境搭建

2.1 器件结构与网格划分

构建GaN肖特基二极管需要准确定义几何结构和掺杂分布。以下是一个典型的网格定义示例：

mesh width=1e6 x.mesh loc=0.0 spac=0.5 x.mesh loc=1.0 spac=0.5 y.mesh loc=0.0 spac=0.0025 y.mesh loc=1.0 spac=0.1 region num=1 material=GaN electrode num=1 name=anode top material=alu electrode num=2 name=cathode bottom doping region=1 uniform n.type conc=5e16

关键参数说明：

y方向近界面处采用0.0025μm的精细网格，以准确求解势垒区电场
衬底掺杂设为5×10^16 cm^-3的n型掺杂
顶部电极为铝(alu)，底部为欧姆接触

2.2 物理模型与求解器配置

完整的物理模型组合对仿真精度至关重要：

mobility material=GaN gansat.n models srh auger albrct fermi temp=$temp ust print method newton carriers=1 maxtrap=10

此处启用了：

氮化物饱和迁移率模型(gansat.n)
SRH复合、Auger复合、带间隧穿(albrct)模型
费米-狄拉克统计(fermi)
通用肖特基隧穿(ust)模型

3. 温度依赖性分析实现

3.1 多温度点扫描设置

研究温度效应需要通过sweep命令实现自动参数扫描：

go internal load infile=diodeex10_aux.in sweep parameter=tempera type=list data="125, 150, 200, 250, 300, 350, 400"

温度范围覆盖125K至400K，模拟典型实验测量条件。

3.2 反向偏置扫描策略

在每个温度点下，采用分段扫描策略获取完整的I-V特性：

solve init solve vanode=0.1 solve vanode=0.2 solve name=anode vstep=-0.1 vfinal=-4.0 solve vanode=-4.5 solve vanode=-5.0 name=anode vstep=-1.0 vfinal=-10 solve vanode=-12.0 name=anode vstep=-2.0 vfinal=-20

这种分段方法在低偏置区采用0.1V的精细步长，在高偏置区逐渐增大步长，兼顾精度与效率。

4. 结果提取与可视化分析

4.1 关键参数提取技术

通过extract命令获取最大电场和阴极电流：

extract name="maxEx_VR04" 1e2*max(probe."maxfieldY") extract name="I_VR04" max(i."cathode")

其中：

maxfieldY探头监测y方向电场分布
电流提取自阴极接触点

4.2 多温度结果对比展示

使用TonyPlot进行数据可视化：

tonyplot -overlay diodeex10_temp_125.log \ diodeex10_temp_150.log \ diodeex10_temp_200.log \ diodeex10_temp_250.log \ diodeex10_temp_300.log \ diodeex10_temp_350.log \ diodeex10_temp_400.log -set diodeex10_1.set

典型输出应包括：