1. GaN HEMT器件仿真探索的背景与整体思路1.1 为什么选择GaN HEMT作为研究对象功率半导体这几年变化很快硅基器件在中低压领域已经逼近材料极限想再往上突破宽禁带材料就成了绕不开的方向。GaN氮化镓和SiC碳化硅是两条主流路线SiC适合高压大功率GaN则在高频、高效率场景里优势明显。GaN HEMT高电子迁移率晶体管凭借其高电子迁移率、高击穿电场和高工作温度在快充、射频功放、数据中心电源等场合越来越常见。我这次探索的出发点很直接想搞清楚一个GaN HEMT从结构参数到电学特性之间的定量关系。实验流片成本太高一轮下来周期长、费用大所以用TCAD工具做前期仿真验证是性价比最高的路径。Silvaco的Atlas模块是我比较熟悉的器件仿真工具它的物理模型库比较全对III-V族化合物半导体的支持也成熟尤其是极化效应这块处理AlGaN/GaN异质结的时候有现成的模型可以用。这个探索适合谁看如果你正在做功率器件仿真、宽禁带半导体研究或者刚接触Silvaco Atlas想找一个完整的实操案例那这篇内容应该能帮到你。我会从结构搭建、材料参数设置、物理模型选择一直讲到特性曲线分析和常见报错处理尽量把每个环节的“为什么”讲清楚。1.2 整体仿真流程的设计考量做器件仿真不能上来就点运行得先想清楚整个流程怎么走。我的思路是分四步结构定义、材料与模型配置、求解设置、结果分析。这四步的顺序不能乱因为后面的每一步都依赖前面的正确性。结构定义阶段要确定器件的几何尺寸包括栅长、栅漏间距、势垒层厚度、缓冲层厚度这些。这些参数不是随便填的它们直接决定了器件的电场分布和电流能力。比如栅漏间距增大击穿电压会提高但导通电阻也会增大这是一个典型的trade-off。材料与模型配置阶段是GaN HEMT仿真的核心难点。和硅器件不同GaN HEMT的二维电子气2DEG不是靠掺杂形成的而是靠AlGaN/GaN异质结的自发极化和压电极化效应。所以必须正确设置极化模型和表面态参数否则根本得不到合理的2DEG浓度。求解设置阶段涉及偏压扫描策略。我一般先用零偏压做初始求解然后逐步加漏压和栅压。直接上高偏压很容易不收敛这是新手最常踩的坑。结果分析阶段主要看输出特性曲线Id-Vd、转移特性曲线Id-Vg、跨导、阈值电压这些关键指标再和文献数据做对比验证。1.3 工具选型为什么是Silvaco Atlas市面上做器件仿真的工具不少Sentaurus、COMSOL、Silvaco各有各的强项。我选Silvaco Atlas主要基于几个考虑。第一Atlas对III-V族材料的极化模型支持比较完善。它内置了polarization模型可以自动计算自发极化和压电极化产生的电荷密度这对GaN HEMT来说是刚需。如果用不支持的模型就得手动算极化电荷再作为固定电荷加进去麻烦且容易出错。第二Atlas的语法相对直观。它的输入文件是结构化的文本格式器件结构、材料参数、物理模型、求解方式分块写逻辑清晰调试起来方便。相比之下有些工具的脚本化程度太高改一个参数要翻半天文档。第三学习资料多。Silvaco在国内高校和研究所的普及率比较高遇到问题容易找到参考案例。Atlas自带的示例库里有HEMT的模板可以作为起点修改。当然Atlas也有不足的地方比如自热效应的耦合求解比较慢二维仿真做久了容易卡。但对于入门级的结构探索和趋势分析它完全够用。2. GaN HEMT核心结构与物理模型拆解2.1 AlGaN/GaN异质结的极化效应原理GaN HEMT最核心的物理机制就是极化效应。这个东西不理解清楚仿真参数就没法设对。简单说AlGaN和GaN这两种材料都是纤锌矿结构本身不具备中心对称性所以存在自发极化。当AlGaN层外延生长在GaN层上时晶格常数不匹配还会产生压电极化。这两种极化加在一起在AlGaN/GaN界面处形成一个很强的极化电场。这个电场的作用是把GaN一侧的电子吸引到界面附近形成一个很薄但浓度很高的电子层也就是二维电子气2DEG。2DEG的浓度可以到1e13 cm^-2量级迁移率也很高因为电子和施主杂质在空间上是分离的减少了电离杂质散射。在Atlas里极化效应的处理方式是在models语句中加上polarization。同时需要在材料参数里设置好AlGaN的Al组分、晶格常数、压电常数和自发极化系数。这些参数Atlas的材料库里有默认值但Al组分不同的时候合金材料的参数需要按线性插值或者bowing公式计算。注意极化电荷的符号很容易搞反。如果仿真出来的2DEG浓度是零或者负的大概率是极化方向设反了。检查polarization模型的方向设置和材料层的生长顺序。2.2 器件结构参数的确定与依据我这次仿真的基准结构是这样的衬底用GaN缓冲层GaN厚度2微米势垒层AlGaN厚度25纳米Al组分0.25栅长1微米栅漏间距5微米栅源间距1微米。这些参数的选择有讲究。势垒层厚度25纳米是一个折中值。太薄了极化效应弱2DEG浓度不够太厚了栅控能力下降而且应力弛豫的风险增大。文献里常见的范围是20到30纳米我取中间值。Al组分0.25也是常见选择。Al组分越高极化效应越强2DEG浓度越高但同时晶格失配也越大材料质量越难保证。0.25到0.3是比较务实的区间。栅漏间距5微米是为了兼顾击穿电压和导通电阻。如果要追求高耐压可以加大到10微米以上但导通电阻会明显上升。我这次主要看基本特性5微米够用了。栅长1微米属于长栅器件短沟道效应不明显仿真相对容易收敛。如果要做射频应用栅长会缩到0.1到0.25微米那时候短沟道效应和速度饱和就得重点考虑。2.3 关键物理模型的选择逻辑Atlas里物理模型很多全开的话仿真速度慢而且容易不收敛。我的原则是按需开启只选和GaN HEMT直接相关的。必开的模型包括polarization极化效应、srh肖克利-雷德-霍尔复合、fermi费米-狄拉克统计、impact碰撞电离。polarization是基础没有它就没有2DEG。srh用于描述耗尽区的复合电流。fermi在载流子浓度很高的时候更准确2DEG区域浓度高需要用费米统计。impact用于击穿特性分析如果只看输出特性可以暂时不开。可选模型包括mobility相关的模型。GaN材料的迁移率模型比较复杂Atlas提供了cvt、analytical等几种。我一般用cvt模型它考虑了横向电场和纵向电场对迁移率的影响比较适合HEMT结构。不推荐的模型bgn带隙窄化在GaN里效果不明显可以不开。auger复合在GaN里也不是主导机制除非做高电流密度仿真否则没必要开。实操心得模型不是越多越好。每多开一个模型求解器就多一组方程要解收敛难度和仿真时间都会增加。先用最小模型集跑通再根据结果逐步添加。3. Silvaco Atlas仿真实操全流程3.1 网格划分策略与代码实现网格划分是仿真里最容易被忽视但影响很大的环节。网格太粗结果不准网格太细仿真跑不动。我的策略是分区加密在栅极边缘、异质结界面和漏极边缘这些电场变化剧烈的地方用细网格在衬底和电极之间的体区域用粗网格。具体来说异质结界面的纵向网格间距控制在0.5纳米以内因为2DEG的厚度只有几纳米网格粗了根本分辨不出来。栅极边缘的横向网格间距控制在10纳米以内因为那里是电场峰值区域。体GaN区域的网格可以放宽到50到100纳米。下面是我用的网格定义代码片段mesh space.mult1.0 x.mesh loc0.0 spac0.05 x.mesh loc1.0 spac0.01 x.mesh loc2.0 spac0.05 x.mesh loc6.0 spac0.1 x.mesh loc7.0 spac0.05 y.mesh loc0.0 spac0.001 y.mesh loc0.025 spac0.0005 y.mesh loc0.03 spac0.001 y.mesh loc2.0 spac0.05这段代码里x.mesh定义横向网格y.mesh定义纵向网格。spac参数控制网格间距单位是微米。可以看到在y0.025微米也就是异质结界面附近网格间距设到了0.0005微米也就是0.5纳米。注意网格数量直接决定仿真时间。二维仿真里网格数超过5万之后求解时间会明显上升。如果只是做趋势分析没必要追求极致精度适当放宽网格可以大幅缩短迭代周期。3.2 材料参数与极化模型配置材料参数配置是GaN HEMT仿真的关键步骤。Atlas的材料库里有GaN和AlN的默认参数但AlGaN合金的参数需要自己定义。我用的材料定义代码如下material materialGaN user.groupsemiconductor \ eg3003.42 affinity4.1 permittivity9.5 \ mun1000 mup30 nc3002.3e18 nv3004.6e19 material materialAlGaN user.groupsemiconductor \ eg3003.421.0*x \ affinity4.1-1.0*x \ permittivity9.5-0.5*x \ mun300 mup10这里x代表Al组分。带隙、亲和能、介电常数都随Al组分线性变化。迁移率我设得比较保守AlGaN的迁移率远低于GaN因为合金散射很强。极化模型的配置在models语句里models polarization srh fermi cvtpolarization开启极化效应srh开启复合模型fermi开启费米统计cvt开启迁移率模型。还需要设置极化电荷的缩放因子。Atlas里有个polarization.scale参数默认是1.0。如果仿真出来的2DEG浓度和实验差太多可以微调这个因子。但一般情况下不建议乱改先检查材料和结构参数是否正确。3.3 电极定义与偏压扫描设置电极定义决定了器件的测试方式。GaN HEMT一般有源极、漏极、栅极三个电极有时候还有衬底电极。electrode namesource x.min0.5 x.max1.5 y.min0.0 y.max0.0 electrode namegate x.min2.5 x.max3.5 y.min0.0 y.max0.0 electrode namedrain x.min6.5 x.max7.5 y.min0.0 y.max0.0源极和漏极是欧姆接触栅极是肖特基接触。在Atlas里欧姆接触默认处理肖特基接触需要设置功函数。contact namegate workfunction5.1这里功函数设5.1电子伏特对应镍金的肖特基势垒。实际器件的栅极金属功函数会影响阈值电压仿真里要设对。偏压扫描我分两步走。先做转移特性扫描漏压固定在某个值栅压从负扫到正solve init solve vdrain0.1 solve vgate-3.0 vstep0.1 vfinal1.0 namegate save outftransfer.plt再做输出特性扫描栅压固定漏压从零扫到正solve init solve vgate0.0 solve vdrain0.0 vstep0.5 vfinal10.0 namedrain save outfoutput.plt实操心得扫描步长不要设太大。栅压步长0.1伏、漏压步长0.5伏是比较稳妥的。步长太大容易跳过关键拐点曲线不光滑步长太小仿真时间成倍增加。如果遇到不收敛先把步长减半试试。3.4 仿真结果提取与曲线分析仿真跑完之后用TonyPlot打开plt文件看曲线。我主要关注几个指标。转移特性曲线看阈值电压和跨导。阈值电压是漏电流达到某个判据比如1毫安每毫米时的栅压。跨导是转移曲线的斜率反映栅压对漏电流的控制能力。我这次仿真的阈值电压大概在-2.5伏左右跨导峰值在100毫西门子每毫米量级和文献里的数据基本吻合。输出特性曲线看饱和电流和导通电阻。在栅压0伏时饱和电流大概在0.5安培每毫米导通电阻在5欧姆毫米左右。这些数值和同结构器件的实验数据在一个量级说明模型配置是合理的。2DEG浓度可以通过提取异质结界面处的电子浓度来验证。在TonyPlot里切一条纵向截面看界面处的电子浓度峰值。我这次仿真得到的峰值浓度在1e19 cm^-3量级积分到厚度方向大概对应1e13 cm^-2的面密度符合预期。如果结果不对排查顺序是这样的先看2DEG有没有形成没有的话检查极化模型再看阈值电压是否合理不对的话检查栅极功函数和势垒层厚度最后看饱和电流偏小的话检查迁移率模型和源漏电阻。4. 常见问题与排查技巧实录4.1 收敛失败的原因与解决路径收敛问题是Atlas仿真里最让人头疼的。我踩过的坑包括初始求解不收敛、偏压扫描中途报错、网格太密导致矩阵求解失败。初始求解不收敛最常见的原因是材料参数或模型配置有矛盾。比如极化电荷设得太大导致初始电势分布极端求解器找不到解。解决办法是先用最简单的模型集只开polarization和srh跑零偏压跑通了再逐步加模型和偏压。偏压扫描中途报错通常是步长太大。Atlas的求解器是牛顿迭代法步长太大的话初值离解太远迭代发散。解决办法是减小步长或者用method语句调整迭代参数method newton trap maxiters100maxiters控制最大迭代次数默认是50遇到难收敛的情况可以加到100甚至200。网格太密导致的收敛问题比较隐蔽。网格数太多的话矩阵维度大数值精度问题会导致求解失败。这时候需要适当放宽网格尤其是体区域的网格。注意收敛失败的时候不要盲目改参数。先看报错信息Atlas会提示是哪一步、哪个变量不收敛。根据提示定位问题比瞎试效率高得多。4.2 仿真结果与实验偏差的分析方法仿真和实验有偏差是正常的关键是要知道偏差从哪来。如果阈值电压偏负可能是栅极功函数设低了或者势垒层厚度设薄了。阈值电压的近似公式是Vth φb - ΔEc - qNssd / ε其中φb是肖特基势垒ΔEc是导带偏移Nss是表面态密度d是势垒层厚度。任何一个参数变化都会影响阈值电压。如果饱和电流偏小可能是迁移率设低了或者源漏接触电阻设大了。GaN HEMT的迁移率受很多因素影响仿真里用的常数迁移率模型和实际有差距。可以尝试用场依赖的迁移率模型或者根据实验数据反推一个等效迁移率。如果击穿电压偏小可能是碰撞电离系数设得不对或者电场分布不合理。GaN的碰撞电离系数比硅低很多Atlas默认的参数可能需要调整。另外栅极边缘的电场集中效应会降低击穿电压实际器件里会用场板来缓解仿真里也可以加场板结构。4.3 提升仿真效率的实用技巧仿真效率直接影响迭代速度。我总结了几条实用的。第一先用粗网格跑趋势再用细网格跑精确值。趋势分析不需要很高的精度粗网格跑得快能快速判断参数变化的方向。确定方向后再用细网格做精确计算。第二保存中间结果。Atlas支持在扫描过程中保存多个plt文件这样如果后面某一步不收敛不用从头再跑。我一般在每个偏压点都保存一次。第三合理使用method语句。不同的求解方法收敛速度不一样。newton适合大多数情况gummel适合初值很差的时候。可以先用gummel跑几步再切换到newton。第四并行计算。Atlas支持多线程求解在method语句里加threads4可以调用四个核心。对于网格数大的仿真并行能明显缩短时间。常见问题可能原因排查方法解决措施初始求解不收敛模型太多或参数矛盾逐步减少模型先用最小模型集扫描中途报错步长太大查看报错位置减小步长或加maxiters2DEG浓度为零极化模型未开或方向反检查models语句开启polarization并检查方向阈值电压偏负功函数低或势垒层薄对比文献值调整功函数或厚度饱和电流偏小迁移率低或接触电阻大检查材料参数调整迁移率模型击穿电压偏小碰撞电离参数不对检查impact模型调整电离系数4.4 从仿真到实验的衔接建议仿真做得再好最终还是要和实验对照。我建议在仿真阶段就把实验能测的量都提取出来这样流片之后可以直接对比。实验上能测的包括转移特性、输出特性、C-V特性、击穿电压、导通电阻。仿真里对应的提取方法要提前想好。比如C-V特性仿真里可以通过交流小信号分析来提取需要在solve语句里加ac分析。另外仿真里用的结构参数要和版图一致。栅长、栅漏间距这些在版图里是多少仿真里就设多少。不要仿真用一个值流片用另一个值那样对比就没有意义了。表面态密度是一个很难确定的参数。实验上可以通过栅极漏电和电流崩塌效应来间接推断仿真里可以扫描不同的Nss值找到和实验最吻合的那个。我一般从1e12到1e13 cm^-2扫一遍看哪个值对应的阈值电压和实验最接近。5. 结构参数扫描与器件性能优化5.1 势垒层厚度对2DEG浓度的影响势垒层厚度是影响2DEG浓度最直接的参数之一。我做了从15纳米到35纳米的扫描步长5纳米。结果是厚度从15纳米增加到25纳米时2DEG浓度从大约6e12 cm^-2增加到1e13 cm^-2增长很明显。但从25纳米继续增加到35纳米时浓度增长变缓只增加到1.1e13 cm^-2左右。这说明极化效应在薄势垒层时更敏感厚到一定程度后趋于饱和。这个趋势和理论分析一致。2DEG浓度和势垒层厚度的关系不是线性的而是先快后慢。原因是极化电荷产生的电场随厚度增加而增强但同时势垒层内的应变弛豫也会随厚度增加而加剧部分抵消了极化效应。从器件设计角度25纳米是一个比较合适的值。再薄的话2DEG浓度不够导通电阻大再厚的话栅控能力下降而且材料生长难度增加。5.2 Al组分对阈值电压的调控规律Al组分的变化会同时影响极化效应和导带偏移对阈值电压的影响比较复杂。我扫描了Al组分从0.15到0.35的范围。结果是Al组分增加时阈值电压向负方向移动。Al组分0.15时阈值电压约-1.8伏Al组分0.35时约-3.2伏。这个趋势的原因是Al组分增加极化效应增强2DEG浓度增加需要更负的栅压才能耗尽沟道。同时Al组分增加也会增大导带偏移进一步降低阈值电压。但Al组分不能无限增加。超过0.3之后AlGaN和GaN的晶格失配增大材料容易出现裂纹和缺陷实验上很难做出高质量的外延片。所以仿真里可以扫到0.35看趋势但实际设计一般不超过0.3。5.3 栅漏间距与击穿电压的权衡栅漏间距是决定击穿电压的关键参数。我扫描了3微米到10微米的范围。结果是击穿电压随栅漏间距增加而提高但不是线性的。3微米时击穿电压约80伏5微米时约150伏10微米时约280伏。增长速率逐渐放缓因为击穿电压主要受栅极边缘电场集中效应控制间距增大到一定程度后电场分布趋于均匀继续增大间距的收益递减。同时导通电阻随栅漏间距增加而线性增大。5微米时的导通电阻比3微米时大约60%。所以这是一个典型的trade-off要高耐压就得牺牲导通性能。实际设计中会根据应用需求来选。快充应用追求低导通电阻栅漏间距会取小一些功率电源追求高耐压栅漏间距会取大一些。仿真可以快速给出不同间距下的击穿电压和导通电阻帮助设计者做决策。5.4 钝化层与表面态的影响初探表面态是GaN HEMT里一个很特殊的问题。AlGaN表面存在大量的悬挂键会俘获电子形成虚拟栅效应导致电流崩塌。仿真里可以通过设置表面态密度来模拟这个效应。我在AlGaN表面加了一层固定电荷密度从1e12到1e13 cm^-2扫描。结果是表面态密度增加时阈值电压向正方向移动饱和电流下降。这是因为表面态俘获了电子减少了2DEG浓度。当表面态密度达到1e13 cm^-2时饱和电流下降了约30%。实际器件里会用氮化硅钝化层来降低表面态密度。仿真里可以通过在AlGaN表面加一层钝化层材料并设置较低的界面态密度来模拟钝化效果。我对比了有钝化和无钝化的情况有钝化时饱和电流恢复了约20%说明钝化确实有效。这个探索还比较初步表面态的物理机制很复杂包括俘获时间常数、发射时间常数等Atlas里可以用trap语句来更精确地描述。后续可以在这方面深入。