1. 项目概述从极化效应到高频开关的仿真主线很多同学拿到Silvaco TCAD第一反应是找个现成的MOSFET例子跑通流程。真开始做GaN HEMT时才发现问题远不是“跑通”那么简单——同样的结构、同样的电压算出来的漏电流可能差两个数量级阈值电压偏得离谱动态导通电阻怎么调都压不下去。我这些年用Silvaco TCAD做GaN功率器件仿真最大的体会是GaN HEMT的设计和仿真核心压根不在网格和求解器而在物理机制。极化效应决定了二维电子气的形成二维电子气又决定了一切静态和动态特性从阈值电压、跨导、导通电阻到高频小信号参数、开关损耗全都挂在同一条物理链路上。这篇文章我会从AlGaN/GaN异质结的极化电荷讲起把Silvaco TCAD里做GaN HEMT的完整流程拆开包括结构搭建、网格设计、极化模型设置、直流特性提取、高频小信号AC分析、动态Ron仿真以及常见的陷阱参数调试方法。适合三类人看刚接触TCAD仿真、准备把GaN功率器件作为研究课题的学生已经在跑Silvaco但结果总跟实测对不齐的工程师还有做高频开关电源、射频前端需要理解器件内部机理的电路设计者。我不打算写成一本文档翻译更多是把那些手册里不写、但实际调试时一定会踩的坑一个个讲清楚。1.1 仿真对象不是“器件”而是极化电荷传统硅MOSFET的导电沟道靠栅压感性积累电子仿真时只要把掺杂、氧化层厚度、功函数设对结果八九不离十。GaN HEMT完全是另一回事。AlGaN势垒层和GaN沟道层之间没有掺杂导电的二维电子气来自自发极化与压电极化的净效应。AlGaN层的自发极化强度比GaN大加上AlGaN受到GaN衬底拉伸应变产生的压电极化在异质结界面形成固定正极化电荷。为了维持电中性界面附近的导带向下弯曲电子被限制在极薄的三角形量子阱里形成密度高达1e13 cm^-2量级的二维电子气。所以在TCAD里第一件事不是画网格而是确认极化电荷算得对不对。如果极化模型没开、参数没设对后面的所有仿真都没有实际意义。反而可以这样理解GaN HEMT的仿真本质上就是对极化电荷、陷阱电荷和能带弯曲三者平衡关系的求解。二维电子气只是这只方程组的“解”。1.2 静态指标和动态指标必须放在同一条链路上看GaN HEMT的卖点有两个一个是高频率一个是高开关速度。但这两条路对器件结构的要求并不完全一致。做射频功放时希望栅长短、跨导高、电容小尽量把fT和fmax顶上去做功率开关时更关心击穿电压、导通电阻、动态Ron和电流崩塌。用Silvaco做设计验证时不能只跑一条转移特性曲线就算完成而是要把直流、小信号AC、瞬态开关放在同一个器件模型下连起来看。我自己的习惯是建立一张“指标链路表”先定Al组分、势垒厚度、栅长、缓冲层掺杂等基本结构参数再用TCAD输出二维电子气面密度、阈值电压、跨导峰值、特征频率、动态Ron增量。每调一个参数就同时看静态和动态两方面的变化。很多初学者会忽略陷阱对高频参数的影响实际上缓冲层陷阱和表面态不仅拖累动态Ron也会通过改变有效输出电导和跨导频率响应的方式把fmax拉下来。这些只能靠TCAD联合求解单纯依赖解析公式很难描述。2. 仿真环境与器件结构搭建把GaN HEMT在Silvaco里“长”出来2.1 软件安装和License要点先说环境。Silvaco TCAD是商业软件License请务必走正规渠道学生可以找导师或学院申请教育授权公司用户走官方商务渠道。不要图省事去用来路不明的破解包一方面软件稳定性没法保证出现网格报错你根本分不清是模型问题还是安装问题另一方面TCAD仿真结果要写进论文或产品报告来源不明的工具链天然就是学术和工程风险。安装本身不复杂Linux和Windows都有对应版本。我更推荐在Linux服务器上跑大规模参数扫描顺手用脚本批量改结构参数确实比Windows界面下逐个点方便得多。安装完第一件事不是直接跑到GaN例子而是先跑一遍自带验证算例确许可执行文件和License服务通信正常。Silvaco自带的deckbuild和tonyplot是日常最常用的两个入口前者管命令行和仿真流程后者看结构、能带、载流子分布和各类曲线。初次上手时我习惯在DeckBuild里把命令逐行敲一遍遇到报错可以直接定位到具体语句。2.2 结构尺寸和网格设计沟道附近的网格决定成败GaN HEMT的典型结构从上往下大致是钝化层Si3N4等、AlGaN势垒层、AlN插入层可选、GaN沟道/缓冲层。源漏是欧姆接触栅极一般做肖特基接触。在Silvaco ATLAS里建结构时我用的是二维截面横向x方向从源到漏纵向y方向从衬底底到器件表面。以一颗典型功率GaN HEMT为例我会把GaN缓冲层厚度设为3 μmAlGaN势垒厚度设为20~25 nmAl组分设在0.25左右栅长取0.5~1 μm源漏间距根据耐压需求放到3~5 μm。ATLAS里最基础的命令是这样组织的go atlas mesh x.mesh loc0.0 spacing0.20 x.mesh loc0.5 spacing0.05 x.mesh loc1.0 spacing0.02 x.mesh loc1.5 spacing0.005 x.mesh loc2.0 spacing0.01 x.mesh loc3.0 spacing0.05 y.mesh loc0.0 spacing0.30 y.mesh loc2.8 spacing0.10 y.mesh loc3.0 spacing0.002 y.mesh loc3.03 spacing0.0005 y.mesh loc3.05 spacing0.002注意看y方向的网格间距。在GaN缓冲层主体里网格可以放松到0.3 μm甚至更粗但在AlGaN/GaN异质结界面附近网格必须密到0.5~2 nm量级。原因很简单二维电子气被限制在量子阱里厚度只有几纳米载流子浓度峰值和能带弯曲都发生在这个区域。网格太疏算出来的二维电子气峰值浓度会被严重低估阈值电压也会漂移。横向x方向则要在栅极边缘和源漏接触边界加密因为这些地方电场集中击穿特性和动态Ron仿真对栅边缘的网格密度非常敏感。结构区域定义可以用类似下面的方式region num1 materialGaN y.min0.0 y.max3.0 region num2 materialAlGaN x.comp0.25 y.min3.0 y.max3.025 region num3 materialSi3N4 y.min3.025 y.max3.05这样定义出来的结构和实际外延片基本对应底部是GaN缓冲/沟道中间是AlGaN势垒表面是钝化层。源漏和栅极通过电极语句落在AlGaN表面electrode namesource x.min0.0 x.max0.5 y.min3.0 y.max3.02 electrode namegate x.min1.2 x.max1.7 y.min3.0 y.max3.02 electrode namedrain x.min2.5 x.max3.0 y.min3.0 y.max3.02源漏接触可能还需要在AlGaN表面下方加n型重掺杂区模拟欧姆接触金属烧结后对势垒层的破坏和合金渗透这里可以先用高浓度n型掺杂加金属欧姆接触来简化。栅极材料用金属功函数来表征肖特基接触的高度典型值在4.6~5.2 eV之间具体按工艺校准。2.3 材料参数与极化模型的落地Al组分、应变与界面电荷材料参数我习惯分成两个层次处理。第一层是基础参数包括带隙、电子亲和能、介电常数、迁移率模型参数这些Silvaco材料库里大多有默认值但对AlGaN这个三元合金软件通常需要根据Al组分进行插值所以x.comp0.25这个参数很关键它会同时影响带隙、极化强度和能带偏移。第二层是极化参数这是GaN HEMT仿真的灵魂。在Silvaco里极化效应有两种落地方式。第一种是用MODELS语句里的极化模型让软件根据Al组分和应变状态自动计算自发极化与压电极化再在异质结界面生成等效固定电荷。第二种是更直接的等效界面电荷法在INTERFACE语句里给AlGaN/GaN界面指定一个固定电荷密度比如QF1e13单位是cm^-2。我个人的建议是前期快速验证结构时用第二种简单、收敛快做精细设计或研究应变工程时用第一种因为它能反映出Al组分变化、AlN插入层、应力释放等工艺细节对极化电荷的影响这是固定电荷法做不到的。极化电荷的量级和方向务必要算清楚。AlGaN与GaN的自发极化强度差大约在0.01~0.03 C/m^2量级压电极化会随Al组分和应变速增加。按典型Al0.25Ga0.75N/GaN结构估算界面净极化正电荷面密度大约在1e13 cm^-2左右这和实测霍尔得到的二维电子气面密度基本同一量级。如果仿真算出来的二维电子气面密度要么接近零、要么超过2e13 cm^-2优先检查极化参数有没有设错方向、Al组分有没有在区域定义里正确传递。注意在比较不同文献的极化参数时一定要确认符号约定。有些文献用负值表示自发极化方向指向衬底有些把压电极化常数定义成压缩应变为正直接照抄很容易得到反向的界面电荷仿真结果自然全乱。3. 直流特性仿真转移、输出、阈值电压与校准3.1 转移特性与阈值电压提取结构建好、极化模型打开之后第一步先跑转移特性。常见的做法是固定漏极电压比如Vds10 V把栅极电压从-6 V扫到2 V步长取0.1 V输出Id随Vgs的变化曲线。在Silvaco里用SOLVE语句逐点扫描或者用LOG语句把扫出来的电流存成文件最后在TonyPlot里查看。阈值电压的提取方式有几种可以用线性区外推法也可以取跨导峰值处切线与横轴的交点。GaN HEMT通常是耗尽型器件阈值电压在-3~-1 V附近。仿真时如果发现阈值电压严重偏正比如变成1 V以上我一般先怀疑极化界面电荷偏小或者栅金属功函数设得太高。反过来如果阈值电压偏负太多比如-6 V以下大概率是极化电荷偏大或者势垒层厚度设得太厚。另外要注意跨导曲线。GaN HEMT的峰值跨导通常出现在Vgs接近0 V附近典型值在100~300 mS/mm量级。如果扫描出来的峰值跨导太低先查网格在沟道界面够不够密再查迁移率模型参数。很多时候不是物理模型错了而是网格平均化了二维电子气的真实分布导致有效沟道迁移率被稀释。3.2 输出特性、导通电阻与饱和行为输出特性是固定几个Vgs值扫描Vds从0到20 V或更高。GaN HEMT的输出特性在线性区斜率对应导通电阻Ron饱和区则反映沟道夹断和速度饱和效应。对功率器件来说我最关心的是低Vds区的线性段这里算出来的Ron可以直接和实测的导通电阻做对比。仿真的导通电阻通常由三部分组成源接触电阻、沟道电阻、漏接触电阻。TCAD默认会把接触电阻算得很理想实测里欧姆接触电阻往往是总电阻的大头。所以对比仿真和实测Ron时我会先在TCAD里用TLM结构单独提取接触电阻率再把固定串联电阻加到源漏两端而不是直接拿理想结构硬碰硬。输出特性还有一个常见的坑Vds较大时GaN高场迁移率模型选择不当饱和电流可能比实测高很多。GaN沟道里的峰值电场很容易超过200 kV/cm必须使用带速度饱和的迁移率模型。Silvaco里常用的有场依赖迁移率模型和能带输运模型前者适合快验证后者精度更高但计算慢。建议至少保证高场饱和速度在1e7 cm/s量级否则开关瞬态仿真的电流斜率会失真。4. 高频性能仿真fT和fmax的提取与优化思路4.1 小信号AC分析的偏置点选择高频特性不是靠加大Vds硬扫出来的而是要在特定直流偏置点上做小信号AC分析。我会先把器件偏置到饱和区典型条件如Vgs0 V、Vds10~15 V让它工作在跨导较高、输出电导相对稳定的区域。然后加一个小信号正弦激励频率从1 MHz扫到100 GHz得到不同频率下的Y参数或S参数再从中提取电流增益和功率增益。Silvaco里的AC分析可以用类似下面的方式solve vdrain15 vgate0 log outfileac.log solve ac freq1.0e6 fstep1.25 nstep40 an1这里的核心逻辑是小信号参数由直流工作点处的偏导数决定因此AC仿真结果能不能信前提是直流通导点偏得准。如果Vgs设得太接近阈值跨导很小算出来的fT没参考意义如果Vds太小器件还在线性区电流增益曲线会怪异得像一个放大器但根本没放大能力。4.2 从Y参数中拆解寄生参数提取fT和fmax传统做法是画电流增益H21和功率增益MAG/MSG的曲线分别找它们降到0 dB对应的频率。经验公式也可以做估算fT约等于gm除以2π乘以(Cgs加Cgd)fmax还额外受栅电阻Rg和Cgd影响。TCAD里的优势是可以直接看内部电荷和电容分布把寄生拆开。有一次我把一个GaN HEMT的fmax调得很低一开始以为是陷阱问题后来在TCAD里输出局部电荷分布发现栅漏之间的电容被钝化层厚度和场板结构拉得很高。这个在实测里很难一眼定位但在TCAD里通过扫描场板长度能很清楚看到Cgd随场板长度变化的趋势。所以高频参数仿真不光是“算一个数”更重要的是告诉设计者“瓶颈在哪”。4.3 高频优化栅长、势垒厚度和寄生电阻的权衡高频仿真最常见的结构优化是缩栅长。栅长从1 μm缩到0.25 μmfT通常会涨3到5倍但栅长缩短也会带来短沟道效应加剧、阈值电压漂移、输出电导上升。在Silvaco里做参数扫描时我会同时看阈值电压漂移量和fT的变化找一个“保持阈值漂移在可接受范围、fT尽量高”的栅长区间。另一个容易被忽略的点是势垒厚度。AlGaN势垒减薄栅控能力增强对短沟道效应有抑制但极化电荷也会随之减小导致二维电子气密度下降、导通电阻变大。所以对射频器件势垒厚度通常在15~20 nm之间做折中对功率开关可以适当加厚势垒换更低的Ron但要接受阈值电压更负、更难驱动。仿真中把这个权衡关系用一组曲线摆出来比在实测中反复流片试错要高效得多。5. 开关性能仿真动态Ron、电流崩塌与瞬态开关波形5.1 动态Ron的核心测量思路GaN功率器件的开关损耗跟静态Ron不是一回事。器件从高电压关断状态切换到导通状态后导通电阻可能在微秒到毫秒量级内明显上升这就是动态Ron退化业内也叫电流崩塌。物理根源是关态高压下陷阱俘获电子开通瞬间这些被俘获的电子还没有全部释放挤占了二维电子气的通道空间导致导通电阻变高。在Silvaco里做动态Ron仿真关键是复制出“先高压关断再低压开通”的时序。典型操作流程是先把器件偏置在Vgs-6 V、Vds100 V的关态让陷阱有足够时间俘获电子然后快速把栅压切到0 V、把Vds降到1 V左右测量导通电流用这个电流反算出的Ron和静态Ron做对比。两个Ron的比值就是动态Ron退化倍数这是功率GaN器件最核心的可靠性指标之一。仿真时需要在陷阱模型里显式加入俘获和发射动力学不能只用固定界面电荷。Silvaco的陷阱语句可以定义陷阱能级、俘获截面、密度和空间分布关键是能级位置要选对。我做缓冲层陷阱校准时通常把能级放在导带底以下0.8 eV附近密度在1e16到1e18 cm^-3范围扫描看哪个参数组合能复现实测动态Ron的变化规律。5.2 陷阱参数和电场分布如何联动动态Ron的仿真结果对陷阱空间分布极其敏感。表面态和钝化层界面的陷阱主要影响栅漏区域靠近表面的电场会在表面形成“虚栅”导致开态电流被局部耗尽。缓冲层陷阱则会在高Vds下被沟道热电子注入影响整个沟道下方的电场分布。两种陷阱在TCAD里可以用不同的空间区域定义来区分。我当时调一个电流崩塌严重恶化的器件实测动态Ron涨了40%静态Ron完全正常。初始仿真里只加了缓冲层陷阱结果动态Ron几乎没变化后来把表面陷阱也加进去动态Ron退化趋势才和实测对上。这里有个经验表面陷阱对电流崩塌的作用往往比缓冲层陷阱更直接因为表面陷阱靠近栅漏高电场区俘获电荷对沟道耗尽作用更强。所以做动态Ron仿真时不要只盯着缓冲层表面态一定要一起校。5.3 用MixedMode做开关波形与损耗估算要评估实际开关损耗需要把器件放进一个电路里跑瞬态。Silvaco的MixedMode允许把一个TCAD器件模型和外部电路元件联立外围接栅极驱动电阻、电感负载、续流二极管等。这样就能看到真实的开关波形计算开通损耗和关断损耗。我自己常用的是一个简单的双电阻测试电路高压母线串电感GaN HEMT做下管续流二极管做上管。先让器件关断电感储存能量再给栅极一个脉冲观察Vds下降和Id上升的交叠区域。开关损耗就是Vds和Id乘积在切换时间内的时间积分。TCAD仿真里的好处是能同时输出器件内部的电场和电荷分布判断损耗是从沟道电阻来的还是从电容充放电来的。这对优化栅极驱动电阻和外部栅极回路很有帮助。如果做驱动电路设计栅极驱动信号的质量会直接影响开关波形。我在搭GaN测试板时经常在驱动前级加二阶低通有源滤波器做信号调理滤掉PWM信号里的高频噪声避免误导通。这部分虽然属于电路级设计不归TCAD管但仿真时把驱动信号理想化、实测却带着毛刺两者对不上也是常有的事。建议把驱动信号上升沿的过冲和振铃都纳入瞬态仿真的输入条件。6. 常见问题与排查技巧实录6.1 不收敛的排查顺序TCAD仿真不收敛是常态尤其是GaN HEMT这种强极化、强电场、多陷阱耦合的结构。第一次遇到不收敛别急着怀疑软件坏了按下面顺序查网格检查。沟道界面附近网格尺寸和相邻网格的比例是否突变过大ATLAS对网格的长宽比和相邻间距有隐含要求突变太大直接导致矩阵奇异性。把沟道附近y方向间距收紧到0.5~2 nm往往能解决一大半问题。初始条件。高电压扫描时从零偏置直接跳到100 V很容易不收敛。加一个足够长的ramp或者先用小电压建立初始解再逐级加电压。物理模型。打开的模型越多非线性越强越难收敛。对不同的问题分开建模直流用基础模型动态Ron才加陷阱击穿再加碰撞电离。陷阱的初始化。陷阱俘获和发射本身是强非线性过程尤其在高电场下时间常数跨越量级很大。做瞬态仿真时时间步长要从极小值开始再逐步放大。6.2 二维电子气面密度偏少或偏多先查这几项如果仿真出来的最大电子浓度剖面在AlGaN/GaN界面处没有形成尖锐峰值或者面密度只有1e12 cm^-2量级按优先级排查极化模型是否打开这是最大概率原因。Al组分x.comp是否按预期传给了电荷计算很多结构定义错误不是发生在区域语句里而是发生在模型参数传递时。界面网格是否足够密网格太疏会把二维电子气剖面“抹平”峰值自然低。表面受主态和缓冲层陷阱是否过多陷阱太多会消耗二维电子气导致面密度下降。如果面密度超过2e13 cm^-2就看是不是极化电荷算重复了既开了极化模型又在界面额外加了固定电荷双份电荷叠加结果当然过高。6.3 仿真和实测对不齐怎么办校准才是TCAD仿真真正花时间的环节。我的做法是把校准分成三步第一步校直流先调极化电荷和栅功函数让阈值电压和实测对齐再调接触电阻和迁移率让Ron和饱和电流对齐。第二步校小信号调Cgs、Cgd和寄生电阻让fT和fmax对到同一个量级。第三步校动态调陷阱能级、密度和俘获截面让动态Ron退化倍数和实测尽量一致。三步之间有耦合可能需要反复迭代但顺序不能乱跳过直流直接调动态参数最后一定是伪拟合。注意TCAD仿真参数不是越多越好。模型里的每个自由参数都应该有物理依据或实验支撑。参数一旦超过实测数据的约束能力结论就容易变成“调参调出来的”而不是“仿真设计出来的”。这一点在写论文和做工程评审时尤其要注意。7. 一点实用体会算起来我拿Silvaco TCAD做GaN HEMT也有不少年头了。回头看最值钱的不是学会了多少条命令而是建立了一种“先用简单模型理解物理再用复杂模型逼近现实”的做事方式。初学者最容易掉进去的坑是一上来就把所有模型全打开、所有参数全调一遍结果算出来的东西连自己都解释不了。我现在每做一个GaN HEMT仿真都会先问三个问题极化电荷的源在哪里陷阱对静态和动态特性的影响分别是什么仿真结论对应到工艺或电路设计下一步该改什么想清楚这三个问题仿真才不会变成自娱自乐的数字游戏。Silvaco TCAD本身只是个工具真正推动设计迭代的是对GaN器件物理的理解以及把仿真和实测反复对标的那股笨功夫。