1. 为什么GaN HEMT的仿真不能只靠“理想模型”糊弄过去做功率器件仿真的朋友大概都有这种体会拿到一个GaN HEMT的仿真任务第一反应是打开TCAD工具画结构、定掺杂、加偏压跑出来的IV曲线看着挺像那么回事但一到高频小信号或者开关瞬态结果就跟实测对不上。问题出在哪绝大多数情况下是极化效应没有被正确处理。GaN HEMT和传统的硅基MOSFET在物理机制上有本质区别。硅器件靠掺杂浓度和栅极电压控制沟道载流子而GaN HEMT的核心导电沟道——也就是二维电子气2DEG——主要来源于自发极化和压电极化在AlGaN/GaN异质结界面处产生的极化电荷。这些极化电荷不是掺杂进去的而是由晶体结构本身的对称性破缺和晶格失配应变决定的。换句话说你如果不把极化物理模型配置对2DEG的浓度就是错的后面所有电学特性都是空中楼阁。Silvaco TCAD在化合物半导体仿真领域积累了很多年它的极化模型和压电模型是内置的但“内置”不等于“自动正确”。很多新手直接选默认参数跑结果2DEG面密度算出来跟文献差了百分之三四十然后抱怨工具不准。实际上极化缩放因子、应变弛豫度、表面态模型这几个关键参数需要根据你的外延结构逐项确认。这篇文章面向的是已经有一定TCAD操作基础、正在做GaN HEMT电学仿真或者准备入手高频特性分析的工程师和研究生。我会从极化效应的物理配置讲起一路拆到高频小信号仿真的设置细节把Silvaco TCAD里那些“文档没写清楚但实际很要命”的地方逐个点出来。你不需要有很深的半导体物理背景但至少要能看懂TCAD的基本命令结构。2. 极化效应在Silvaco TCAD中到底怎么配才靠谱2.1 自发极化与压电极化的物理来源拆解先把这个事情说清楚。GaN是纤锌矿结构沿c轴方向没有 inversion symmetry所以即使没有外加应变晶体内部也会存在一个非零的电偶极矩这就是自发极化。AlGaN的自发极化系数和GaN不同两者叠在一起在异质结界面就会形成一个极化电荷面密度。压电极化则是另一回事。AlGaN和GaN的晶格常数不一样AlGaN层在生长时被GaN衬底“拽着”产生面内应变。这个应变通过压电张量转化为额外的极化贡献。Al组分越高晶格失配越大压电极化越强。两者加起来总极化电荷在AlGaN/GaN界面处吸引自由电子形成2DEG。这就是为什么GaN HEMT不需要故意掺杂就能有很高的沟道电子浓度——通常能做到1e13 cm^-2量级。在Silvaco TCAD中你需要确保以下物理模型被激活models polariz piezo这两个关键词分别对应自发极化和压电极化。但光加这两个词还不够还需要在材料定义区域指定极化相关的参数。2.2 材料参数配置中的关键陷阱Silvaco的材料数据库里自带GaN和AlGaN的参数但默认的AlGaN参数通常是针对某个特定Al组分的。如果你的结构用的是Al0.25Ga0.75N而数据库默认是Al0.3Ga0.7N那极化电荷算出来就会有偏差。我的做法是手动定义AlGaN的组分相关参数。在material语句中用x.comp指定Al组分然后让工具根据线性插值或者bowing公式计算极化系数。具体来说material materialAlGaN x.comp0.25 \ polar.scale1.0 \ piezo.scale1.0这里的polar.scale和piezo.scale是两个非常关键的缩放因子。在理想情况下它们都等于1但实际仿真中由于表面态、钝化层、以及外延质量的影响有效极化电荷会比理论值低。很多文献报道的“极化衰减因子”在0.5到0.8之间。你需要根据你的器件结构有没有SiN钝化、栅极是不是绝缘栅结构来调整这两个因子。注意不要为了凑2DEG浓度就随便改缩放因子。正确的做法是先算理论值再跟Hall测试或者C-V测试的实测数据对比反推出合理的缩放因子。如果偏差超过20%先检查结构参数是不是写错了。2.3 应变弛豫度的处理方式AlGaN层如果比较厚超过临界厚度就会发生应变弛豫。弛豫之后压电极化会减弱。Silvaco TCAD里可以通过relax参数来指定弛豫度从0完全应变到1完全弛豫。对于典型的AlGaN/GaN HEMTAlGaN厚度在20-30 nm范围通常假设完全应变即弛豫度为0。但如果你的结构里有厚AlGaN层或者多量子阱结构就必须考虑弛豫。实际操作中我一般会先跑一个完全应变的版本看看2DEG浓度是否合理。如果算出来偏高再逐步增加弛豫度。这个过程需要跟XRD摇摆曲线或者Raman光谱的应变测试结果对照。2.4 表面态与费米能级钉扎的影响这是最容易被忽略的一块。GaN HEMT的2DEG浓度不仅取决于极化电荷还取决于表面态对费米能级的钉扎位置。表面态会捕获部分极化电荷释放的电子降低有效2DEG浓度。在Silvaco中你可以在表面处定义interface语句指定表面态密度和中性能级位置interface qf1e13 y.min0 y.max0 \ neutral0.5这里的qf是表面态密度neutral是中性能级相对于禁带的位置。对于AlGaN表面中性能级通常在导带底下1.5-2.0 eV左右。这个参数直接影响到阈值电压和2DEG浓度调不对的话器件开启电压会差好几伏。我个人的经验是先用默认值跑一遍然后把qf从1e12逐步加到1e13观察阈值电压的移动。通常当qf在5e12到2e13之间时仿真结果跟实测的阈值电压比较吻合。当然具体值取决于你的钝化工艺。3. 从结构定义到网格划分的实操细节3.1 外延结构定义的顺序与坐标规范Silvaco TCAD的结构定义用的是mesh和region语句。对于GaN HEMT典型结构从下往上依次是GaN缓冲层、GaN沟道层、AlGaN势垒层、可选的GaN帽层、钝化层。坐标定义上我习惯把y0放在AlGaN/GaN界面处这样2DEG的位置就在y0附近查看结果比较方便。具体写法mesh space.mult1.0 x.mesh loc0.0 spac0.05 x.mesh loc1.0 spac0.05 y.mesh loc-2.0 spac0.1 y.mesh loc-0.5 spac0.01 y.mesh loc0.0 spac0.002 y.mesh loc0.025 spac0.002 y.mesh loc0.03 spac0.01 y.mesh loc0.1 spac0.05注意在y0附近网格间距要足够细建议在2 nm以下。因为2DEG的厚度只有几个纳米网格太粗的话电子浓度峰值会被“抹平”算出来的浓度偏低。3.2 掺杂定义中的常见错误GaN HEMT的GaN缓冲层通常有背景掺杂一般是碳掺杂或者铁掺杂来获得高阻态。在Silvaco中你可以用doping语句定义均匀掺杂或者高斯分布。一个常见的错误是忘记定义AlGaN层的掺杂。AlGaN势垒层通常是非故意掺杂的但如果你不显式定义Silvaco可能会继承上一层材料的掺杂设置导致势垒层里出现不该有的自由载流子。我的做法是在每个region语句之后紧跟一个doping语句明确指定该区域的掺杂类型和浓度。即使是本征材料也写一个doping conc1e15避免继承错误。3.3 电极定义与接触类型选择GaN HEMT的源漏接触通常是欧姆接触栅极根据器件类型可以是肖特基栅或者MIS栅。在Silvaco中用electrode语句定义electrode namesource x.min0 x.max0.1 y.min0.025 y.max0.03 electrode namedrain x.min0.9 x.max1.0 y.min0.025 y.max0.03 electrode namegate x.min0.4 x.max0.6 y.min0.03 y.max0.03对于欧姆接触需要设置workfunction或者直接用ohmic关键字。对于肖特基栅需要指定金属功函数。GaN HEMT中常用的栅金属是Ni/Au功函数大约在5.1 eV左右。提示如果你的仿真中出现“接触处载流子浓度异常”或者“电流不收敛”先检查电极的workfunction设置。肖特基接触的功函数差0.2 eV阈值电压就会移动将近0.2 V。4. 直流特性仿真从2DEG浓度到转移特性4.1 零偏压下的2DEG浓度提取结构建好之后第一步是跑零偏压下的解。在solve语句中先不做任何偏压扫描直接求解泊松方程和连续性方程solve init solve previous save outfhemt_zero.str然后可以用TonyPlot打开结果文件查看导带底和电子浓度的空间分布。在AlGaN/GaN界面处你应该能看到一个尖锐的电子浓度峰值典型值在1e19到1e20 cm^-3量级积分得到的面密度在8e12到1.5e13 cm^-2之间。如果这个值明显偏低比如只有1e12大概率是极化模型没生效或者缩放因子设得太小。如果明显偏高超过2e13检查是不是表面态密度设得太低或者Al组分写高了。4.2 转移特性扫描的设置要点转移特性是GaN HEMT最核心的直流特性。在Silvaco中用solve语句做栅压扫描solve vgate0.0 namegate solve vgate-1.0 namegate solve vgate-2.0 namegate ...扫描步长建议先粗后细。先用0.5 V步长扫一遍找到阈值电压的大致位置然后在阈值附近用0.1 V步长细扫。漏极电压一般设在3-5 V保证器件工作在饱和区。但要注意GaN HEMT在高漏压下会有自热效应和陷阱效应如果只做直流仿真这些效应是体现不出来的。所以直流仿真得到的跨导和阈值电压跟脉冲IV测试的结果更接近。4.3 跨导与阈值电压的提取方法跨导gm的定义是漏极电流对栅极电压的导数。在Silvaco中你可以用extract语句自动提取extract namegm_max max(curve(vgate, -deriv(idrain, vgate)))阈值电压的提取方法有好几种最常用的是线性外推法在gm最大的点做切线与vgate轴的交点就是阈值电压。Silvaco的extract语句也支持这种提取方式。我个人的经验是仿真得到的阈值电压通常比实测偏负0.3-0.5 V这主要是因为仿真中没有考虑栅极金属与AlGaN之间的界面态。如果你在栅极下方加一层薄介质比如2 nm的SiN阈值电压会正移更接近实测值。5. 高频小信号仿真从S参数到截止频率5.1 小信号AC仿真的基本设置高频特性仿真是GaN HEMT TCAD中最考验设置功底的部分。Silvaco中做小信号仿真需要在直流工作点的基础上叠加一个小幅度的AC信号。首先你需要把器件偏置在合适的直流工作点。通常选在gm最大的栅压和Vds10 V左右。然后solve ac freq1e9这会在1 GHz频率下求解小信号响应。但单频点不够你需要扫频solve ac freq1e9 fstep1e9 nfreq10这会从1 GHz扫到10 GHz步长1 GHz。5.2 S参数提取与去嵌入Silvaco可以直接输出Y参数或者S参数。但要注意仿真得到的S参数是包含寄生参数的。GaN HEMT的寄生电容和寄生电阻对高频特性影响很大特别是截止频率fT和最大振荡频率fmax。去嵌入的方法通常是先仿真一个“开路”结构去掉本征沟道提取寄生电容再仿真一个“短路”结构提取寄生电感。然后用Y参数减法得到本征参数。在Silvaco中你可以通过定义多个电极和不同的偏置条件来模拟这些结构。这个过程比较繁琐但如果不做去嵌入提取的fT会偏低20%-30%。5.3 截止频率fT和最大振荡频率fmax的计算fT的定义是电流增益降到1时的频率。在Silvaco中你可以从Y参数计算h21h21 |Y21 / Y11|然后找到|h21|1对应的频率。fmax则是最大稳定增益MSG或者最大可用增益MAG降到1时的频率。对于典型的GaN HEMT栅长0.25 um时fT大约在30-50 GHzfmax在60-100 GHz。如果你的仿真结果明显低于这个范围先检查寄生参数是不是没去嵌入再检查栅极电阻是不是设得太大。注意高频仿真中网格尺寸对结果影响很大。栅极下方的网格如果太粗本征延迟时间会被高估fT偏低。建议栅极下方网格间距不超过5 nm。6. 常见问题与排查技巧实录6.1 仿真不收敛的典型原因与对策GaN HEMT仿真不收敛是家常便饭。最常见的原因有三个一是极化电荷太大导致初始猜测的解偏离太远二是网格太粗特别是异质结界面处三是偏压步长太大。对策方面我的习惯是先用method语句设置较小的阻尼因子和较多的迭代次数method newton trap itlimit50然后在solve语句中用ramp方式逐步加偏压而不是一步到位。比如从0 V到-5 V分成20步每步0.25 V。如果还是不收敛可以试试先解泊松方程solve poisson再解漂移扩散方程solve drift最后耦合求解。6.2 2DEG浓度异常偏低的排查路径如果你发现2DEG浓度比预期低很多按以下顺序排查排查项可能问题检查方法极化模型未激活polariz/piezo查看models语句Al组分设置错误或未设置查看material语句缩放因子polar.scale太小尝试设为1.0表面态qf太大降低qf值网格界面处太粗加密y方向网格应变弛豫relax参数不为0检查relax设置6.3 高频仿真中S参数异常的调试方法S参数异常通常表现为S21在低频就下降很快或者S11出现不合理的谐振。这时候先检查电极的寄生参数设置。Silvaco中可以通过resist和capacit语句给电极加寄生电阻和电容。另外AC仿真的频率步长也会影响结果。如果步长太大谐振峰可能被跳过。建议在谐振频率附近用更细的步长。6.4 收敛性与精度的平衡技巧TCAD仿真永远在收敛性和精度之间做权衡。我的经验是直流仿真可以适当放宽收敛容差tol设为1e-5换取更快的收敛速度但高频仿真必须收紧容差1e-8否则S参数的相位误差会很大。还有一个技巧先用粗网格跑通整个流程确认物理模型没问题再逐步加密网格做最终仿真。不要一上来就用最细的网格那样调试起来非常痛苦。7. 一些踩过坑之后才明白的经验极化缩放因子不是拟合参数。很多人把它当成万能旋钮2DEG浓度不对就调它。但缩放因子应该反映的是真实的物理效应——表面态、钝化质量、外延缺陷。如果你发现需要把缩放因子调到0.3才能匹配实测那大概率是结构参数或者表面态设置有问题而不是极化本身的问题。高频仿真的去嵌入比直流仿真重要得多。直流仿真中接触电阻的影响可以通过校准工作函数来补偿。但高频下寄生电容和电感会直接改变S参数的形状。如果你不做去嵌入提取的fT和fmax只能作为相对参考不能跟实测对比。网格密度对高频特性的影响远大于对直流特性的影响。直流特性主要取决于总电荷量网格粗一点总电荷量变化不大。但高频特性取决于载流子的输运延迟网格粗了延迟时间被高估fT就偏低。我一般会在栅极下方和漏极边缘用0.5-1 nm的网格其他地方可以放宽到5-10 nm。表面态模型的选择会影响阈值电压的长期稳定性。如果你只做单次仿真可能感觉不到。但如果你做可靠性仿真比如BTI或者热载流子表面态的动态捕获和释放就会导致阈值电压漂移。Silvaco中有专门的陷阱动力学模型但配置起来比较复杂建议先把直流和高频的基础仿真做扎实再考虑这些进阶内容。最后分享一个实用的小技巧在Silvaco中你可以用save语句在扫描过程中保存中间结果文件然后用TonyPlot批量对比不同偏压下的能带图和浓度分布。这对于理解器件物理非常有帮助比只看最终的IV曲线直观得多。