做功率GaN器件的人多少都绕不开Silvaco TCAD。这几年市面上的CoolGaN、GaN Systems、纳微的器件规格书里写着650V、15A、Qg几纳库看着数据漂亮但真要评估一颗管子能不能用在你的电源方案里击穿电压余量够不够、驱动损耗多大、开关过程有没有振荡光看规格书是不够的。我自己做同步整流和LLC方案选型时踩过不少坑后来干脆把Silvaco TCAD这套仿真流程跑通用结构仿真把器件物理层的问题看透很多疑虑在流片和打样之前就能排除掉。这篇文章就围绕英飞凌CoolGaN这颗代表型增强型GaN HEMT讲清楚怎么用Silvaco TCAD做三个关键指标仿真BV击穿电压、Qg栅极电荷、开关损耗。适合做功率器件设计的工程师、做电源应用选型的硬件工程师以及正在学GaN器件物理的研究生。内容偏实操我会把模型搭建、参数设置、仿真陷阱这些细节尽量讲透你跟着操作基本能复现一套完整的评估流程。提示文中的所有结构参数和掺杂浓度均以650V增强型CoolGaN器件的典型公开参数为基础结合常见外延工艺经验设定。Silvaco相关版本以Atlas仿真器为准。1. 器件结构与物理模型搭建仿真结果准不准七成看这一步1.1 先搞清楚CoolGaN的器件结构再动手CoolGaN本质上是一个p-GaN栅极增强型HEMT和普通耗尽型GaN HEMT最大的区别是栅极下面多了一层p-GaN帽层。这层p-GaN的作用是耗尽下方AlGaN势垒层中的二维电子气让器件在零栅压时处于关断状态也就是常断。仿真第一步就是把器件结构按尺寸画对如果层结构和实际差太远后面的BV和动态特性全是错的。我常用的方法是用Silvaco Atlas的DevEdit或者直接在Atlas语法里写结构定义。650V CoolGaN的关键结构参数我整理了一个参考表结构层材料厚度掺杂/组分作用衬底SiC或硅几百微米仿真只用顶层高阻支撑与散热缓冲层GaN2~4μmC或Fe掺杂补偿耐压与漏电控制沟道层GaN0.3~0.5μm非故意掺杂形成二维电子气势垒层AlGaN15~25nmAl组分20%~25%极化诱导二维电子气栅帽层p-GaN70~100nmMg掺杂1e19量级实现增强型钝化层SiN200~300nm绝缘表面保护与场板介质需要说明的是650V耐压主要靠缓冲层和衬底之间的垂直耐压结构以及栅场板、源场板的横向电场管理。仿真时直接把缓冲层做到4μm左右底下用高阻衬底代替能大大减少网格数量又不影响表面电场分布的计算。1.2 极化电荷模型是GaN仿真的灵魂GaN和AlGaN都是极性材料异质结界面存在自发极化和压电极化正是这两种极化的净效应在AlGaN/GaN界面诱导出了高密度二维电子气。Silvaco TCAD中极化电荷通常用Interface Charge来等效这是仿真中最容易被忽略也最容易出错的地方。我习惯在AlGaN/GaN界面设置一层正的固定电荷面密度大约在5e12到1e13 cm^-2量级。具体数值根据Al组分和AlGaN厚度用极化模型估算Al组分越高、厚度越大二维电子气密度越高但过度增加会导致阈值电压负向漂移和击穿电压下降所以这里有优化空间。极化电荷的设置语法我放在下面interface x.min1.0 x.max9.0 y.min0.38 y.max0.40 qf1.0e13这个命令是在AlGaN/GaN界面处坐标y0.38~0.40μm之间放一层正固定电荷。实际仿真中我会把q.fixed设成8e12、1e13、1.2e13三组做扫描观察Ids-Vgs转移特性和阈值电压的变化再用实验数据校准。这样做的好处是后面仿真动态特性时二维电子气的密度和分布是合理的。1.3 物理模型选择该开什么模型、不该开什么模型Silvaco Atlas里模型开太多收敛慢还容易振荡开太少击穿特性和开关特性完全不对。针对GaN HEMT我一般开这几个核心模型浓度依赖迁移率模型commob加高场饱和模型fldmob模拟二维电子气的高迁移率和速度饱和效应SRH复合模型模拟缓冲层和表面的陷阱辅助漏电俄歇复合Auger在功率器件中主要影响大注入下的少子寿命对GaN HEMT一般可以关掉加快收敛碰撞电离模型impact用Selberherr模型这是算BV必须开的极化模型在Silvaco中通过接口电荷等效不需要额外开陷阱模型在计算动态特性时单独加涉及L3参数的校准我后面专门说关于迁移率模型有个容易踩的坑Silvaco默认的浓度依赖迁移率是针对硅的直接用在GaN上会把二维电子气迁移率预测得过低。我所在的项目组一般通过calibration曲线把峰值迁移率校准到1500~2000 cm²/V·s量级这样输出的Idsat和跨导才和实测对得上。2. 击穿电压BV仿真网格、陷阱与击穿判据2.1 BV仿真为什么要格外小心网格击穿电压仿真的本质是在器件两端加上高反向电压后计算沟道和缓冲层中电场分布是否超过材料的临界击穿场强。GaN的临界场强大约3.3MV/cm比硅高一个数量级这也意味着电场峰值会非常集中网格稍微粗一点电场峰值就会被严重低估BV结果虚高。我自己吃过大亏一开始图省事在沟道和栅边缘用均匀网格跑BV仿真出来的BV高达1200V而实际器件物理上不太可能到这个值。后来把栅极边缘和漏极侧的网格加密到2nm以下BV一下掉到700V左右这才合理。你可以先用均匀网格粗跑一遍找电场峰值位置再局部加密反复迭代两到三轮。加密的区域主要是栅极边缘靠近漏极一侧、场板末端、AlGaN势垒层内部。这三个位置是电场最容易集中形成尖峰的地方。网格加密可以用下面这种方式定义# 栅极边缘加密区域 region num10 x.min4.8 x.max5.2 y.min0.0 y.max0.5 materialGaN # 局部网格细分 mesh x.min4.8 x.max5.2 spacing0.002 mesh y.min0.0 y.max0.5 spacing0.0052.2 缓冲层陷阱对BV的影响比想象中大很多人跑GaN HEMT的BV仿真默认缓冲层是完美的这其实和实际差很远。功率GaN器件为了抑制纵向漏电通常在缓冲层掺入碳C或铁Fe这些深能级陷阱在关断状态下会俘获电子形成负电荷扩展耗尽层本质上提升了器件的耐压能力。Silvaco里可以用trap语句定义缓冲层陷阱。我做650V仿真时的经验是在缓冲层加一个受主型深能级陷阱能级深度大约在0.5~0.8eV陷阱浓度1e16到1e17 cm^-3这会明显改变纵向电场分布和漏电大小。语法参考trap level0.65 density5e16 donor.sigp1e-15 donor.sige0.0 acceptor.sigp1e-15 acceptor.sige0.0 region2陷阱加的太浓BV高得离谱关态漏电低得离谱和实测完全对不上加得太稀BV偏低。这个参数的标定强烈建议结合缓冲层掺杂浓度做一组浓度扫描再对比规格书上的关态漏电指标通常1μA/mm量级来确定范围。2.3 击穿判据不能只看电流还要看电场BV定义在工程上通常是漏电流达到某个值时的电压比如1mA/mm或1μA/mm。但仿真时由于陷阱模型和碰撞电离模型的理想化电流判据有时候会出现误导尤其是陷阱填充过程很慢时电流还在爬升电压已经数值发散。我更推荐双判据组合使用一是看漏极电流达到规格书的限值二是看整个器件内部的最大电场是否超过GaN临界场强3.3MV/cm。两个条件谁先到达谁作为BV判定点。这种方法在团队做器件鲁棒性评估时比较可靠。2.4 场板对BV的改善效果量化CoolGaN这类功率器件为了缓解栅极边缘电场集中通常会做栅场板或源场板。用TCAD可以很直观地扫描场板长度对BV的影响。我在一次仿真中把场板长度从0.5μm扫到2μmBV从620V提升到780V但超过1.5μm后提升趋缓说明再增加场板长度对电场均衡已经没有太大作用反而会增加寄生电容。场板越长栅漏电容Cgd越大这对开关速度不利。所以仿真BV的同时要同步关注Cgd的变化最终在想关断速度和耐压之间做折中。这个优化过程用Silvaco做非常方便因为结构尺寸是参数化的扫一遍长度BV和Cgd是同步输出的。3. 栅极电荷Qg与开关损耗仿真提取3.1 Qg仿真的三件套C-V积分、恒流充放电、瞬态Vgs-Qg曲线Qg是功率器件栅极驱动设计的核心参数它决定了驱动IC的峰值电流和开关速度。实际规格书里的Qg是在特定测试条件下的积分值TCAD里仿真有三种路径第一用小信号AC分析扫出Ciss、Crss、Coss曲线再对Ciss·dVgs积分求Qg。这个方法快但忽略了非线性电容随电压的剧烈变化误差在10%以上。第二直接给栅极加一个恒定电流源做瞬态仿真得到Vgs随时间变化曲线然后把时间轴换成电荷轴就是标准的Vgs-Qg曲线。这个方法最接近规格书测试方法推荐使用。第三在混合模式仿真电路里搭一个驱动电阻栅极电容的组合用方波驱动从波形里反推Qg。这个更接近实际应用但容易把沟道动态效应混进来不好分离。我最推荐第二种因为它在物理上和测试方法完全一致又不需要搭建完整电路。实现方法是在栅极加一个1mA的恒流源源极接地漏极接一个高压源设定电压比如400V然后做5μs的瞬态仿真监测Vgs波形。3.2 从Vgs-Qg曲线中读出关键电荷参数仿真跑完以后Vgs随时间变化的曲线会呈现明显的三段从零上升到阈值电压Vth的线性段、米勒平台段、平台后的继续上升段。把时间轴乘以电流源大小就得到电荷轴。我用Silvaco自带的TonyPlot或者导出数据后用脚本处理直接读取QgsVgs从0到平台起始点的电荷对应栅源电容充电Qgd平台区的电荷宽度对应栅漏电容密勒效应这是影响开关损耗最关键的电荷QgVgs到达最终驱动电压比如10V时的总电荷一个典型的650V CoolGaN仿真结构栅宽100mm的Qg大约在6~10nC量级其中Qgd占40%~50%。如果仿真出来的Qgd占比太高超过60%说明场板设计偏大或栅漏交叠区过长这个时候要回头优化结构。3.3 开关损耗仿真用双脉冲测试电路提取Eon和Eoff开关损耗是硬开关拓扑中效率损失的主要来源之一TCAD层面和电路仿真层面的测法不完全一样。TCAD做开关损耗通常采用混合模式仿真在电路层面上定义spice元件器件本身是二维结构用数值方法求解这样既保留器件物理细节又能模拟电路行为。我的DPT双脉冲测试电路设置如下母线电压400V对应650V器件的典型中间电压负载电感150μH模拟电机或变压器的感性负载续流管用另一个GaN HEMT的体二极管特性或者用理想二极管模型替代栅极电阻4.7Ω这个值对开关速度影响很大建议做2.2Ω、4.7Ω、10Ω三组扫描第一个脉冲宽度对应电流上升到30A左右电路网表写法大概是# 混合模式器件定义 go atlas # 主开关器件 hemt1 drain1 gate2 source3 # 续流二极管 d1 anode1 cathode4 # 负载电感 l1 in4 out1 150u # 栅极驱动电阻 rg in2 out5 4.7 # 脉冲电压源 vpulse gate5 source3 dc0 pulse(-5 10 0 10n 10n 2u 4u)仿真时间步长要足够小才能在开关瞬间捕捉到di/dt和dv/dt的细节。我的做法是在开关沿附近用1ps的步长平台的静止区间放松到10ps。总仿真物理时间大约4μsSilvaco在8核并行下需要跑几个小时属于可接受范围。3.4 从波形中计算损耗注意数值噪声陷阱开关损耗的计算是从电流电压波形乘积积分得到的。瞬态仿真波形上尤其是在开关沿很容易出现数值振荡如果不加处理直接积分损耗数值可能翻倍甚至成为负数。我的经验是在TonyPlot里对Id和Vds先用移动平均做平滑处理窗口宽度取2~5个仿真步长然后计算pId·Vds再对p分别在t1~t2开通过程和t3~t4关断过程区间积分得到Eon和Eoff。一个典型结果400V/30A开关条件下Eon大约在30~60μJEoff在15~40μJ具体和栅极电阻、驱动电压强相关。用表格整理对比几组仿真结果栅极电阻RgEon (μJ)Eoff (μJ)开关时间 (ns)振铃峰值 (V)2.2Ω3822185204.7Ω55302648010Ω784538460从这个表能清楚看到栅极电阻对损耗和EMI的权衡Rg大损耗增加但过冲和振铃减小。这和实际测试中的趋势完全一致。4. 动态特性仿真的L3陷阱参数校准与效率优化4.1 GaN特有的电流崩塌效应怎么在TCAD里建模GaN HEMT和硅MOSFET最大的不同在于动态导通电阻退化也就是电流崩塌。关态高压下陷阱俘获电子导致导通电阻在开关瞬间比静态值大20%~80%这个现象在TCAD里不建模的话开关损耗会被明显低估。Silvaco TCAD典型的做法是采用M. Li et al.提出的双层陷阱模型一般叫L3模型在AlGaN表面和缓冲层分别设置施主型和受主型陷阱模拟表面态和缓冲层陷阱对二维电子气耗尽的动态过程。这两个陷阱在关态高压下俘获电子在开态恢复时缓慢释放仿真上表现为Rds(on)的瞬态偏高再慢慢恢复。表面施主陷阱参数参考设置# 表面施主陷阱能级在导带下 trap level0.35 density1e12 region1 donor.sigp1e-15 donor.sige0.02 # 缓冲层受主陷阱 trap level0.60 density5e16 region2 acceptor.sigp1e-15 acceptor.sige04.2 L3参数校准的工程化方法L3陷阱模型一共有好几十个参数直接全调不现实。我调参的顺序是先用静态转移特性和C-V特性校准极化电荷和厚度然后固定这些参数只用陷阱参数去拟合动态测试中的软恢复时间、Rds(on)动态曲线和关态漏电。每轮只调1~2个参数控制在5组以内否则容易过拟合。校准目标按优先级排序静态Ids-Vds输出特性误差小于5%转移特性跨导峰值位置和幅值吻合关态漏电在600V时不超过1μA/mm动态Rds(on)恢复时间常数和实测一致通常几百纳秒到几微秒做到前三条BV和Qg的仿真结果已经比较可信。第四条涉及动态特性拟合难度高但如果做LLC和软开关电路仿真这步必须认真调否则开关损耗预测完全失真。4.3 网格与并行计算的效率优化思路GaN HEMT仿真最大的痛点是网格量大、瞬态时长一次完整的DPT仿真可能跑10个小时以上。我做效率优化时一般用三个手段第一在保证电场精度前提下把垂直方向缓冲层的网格稀疏化。因为表面电场和二维电子气主要在顶部0.5μm范围内缓冲层网格可以放大到20~50nm不影响BV但能减少30%的网格量。第二用对称性减半结构。很多场板设计具有左右对称或周期重复特点如果能用半元胞仿真网格量直接减半。但要注意边界条件设置尤其是周期边界不能给错否则击穿电压会严重偏离。第三分阶段仿真。静态QA、C-V、瞬态开关用不同的网格密度跑。跑动态特性时静态部分已经校准过可以用粗网格只在开关沿附近加密这样时间步长可以放大总时间节省一半以上。5. 常见问题与实操排查5.1 器件不收敛怎么排查Silvaco仿真最让人崩溃的就是不收敛。我遇到过的情况统计下来有三类第一类是在高电压下碰撞电离模型打开后出现震荡。解决方法是把impact模型的临界电场参数调平滑或者逐步升压仿真每次迭代增加10V前一电压的收敛解作为下一步初值。第二类是界面电荷设置和网格不匹配导致电荷集中在一个网格点上产生数值尖峰。解法是把interface charge分散到相邻3~5个网格点平滑处理。这和真实物理一致因为二维电子气本身有厚度。第三类是瞬态时间步长过大导致动态陷阱方程发散。要限制maxstep比如瞬态时设置tstep1e-12 tstop5e-6确保每个物理过程的特征时间都被采样到。5.2 陷阱参数和BV结果矛盾怎么办有时候陷阱加浓了反而BV下降乍一看违反直觉。其实这是因为过浓的受主陷阱会在缓冲层形成强负电荷区虽然能辅助耗尽但也可能在缓冲层内部造成局部电场峰值提前触发碰撞电离。遇到这种情况我建议输出电场分布图找出峰值位置如果峰值出现在缓冲层内部适当降低缓冲层陷阱浓度或调整能级深度让电场尽量均匀分布在整个缓冲层。5.3 仿真结果和实测对不上优先检查这三个点第一检查极化电荷密度是不是被随意设置。很多初学者把q.fixed设到1e14结果是二维电子气密度过高阈值电压严重负漂BV虚高。可以用实验测得的阈值电压反推极化电荷密度一般阈值在1~2V范围内对应的界面电荷大约在5e12~1e13。第二检查AlGaN厚度和Al组分。厚度每差1nm二维电子气密度变化可能有10%。这块必须严格按实际外延参数设置不能拍脑袋。第三检查源、漏欧姆接触电阻。Silvaco默认接触可能偏理想接触电阻为0会低估实际器件的导通损耗和开关时间。建议给源漏加上接触电阻典型值0.3~0.6Ω·mm这样动态损耗才合理。5.4 仿真工具的安装与授权环境最后提一句工具链。Silvaco TCAD是商业软件建议直接通过Silvaco官方渠道申请评估license正规评估版支持绝大多数Gan HEMT仿真功能。安装后需要设置好license环境变量在命令行输入atlas确认版本正常加载。终端里配好并行计算池后GaN HEMT结构仿真和双脉冲仿真都能稳定运行。6. 写在最后关于仿真和流片的几点体会我做了几年功率器件仿真最大的体会是TCAD仿真的价值不在于替代实测而在于把器件物理的因果链条拆开。一颗管子BV不够你实测只能看到结果不知道是场板设计的问题还是缓冲层漏电的问题。仿真里可以把电场分布打出来把陷阱电荷分布打出来一层一层剥开看这个确实是实测做不到的。关于CoolGaN这类增强型GaN HEMT我额外提一个建议评估器件时不要只盯着静态参数和BV务必结合动态特性一起看。很多应用上的失效和损耗问题根源都在动态导通电阻退化和开关过程中的电压过冲上。仿真阶段把这两个维度打通后面做电源方案会轻松很多。这篇文章里的参数配置和流程我在多个650V GaN器件项目上验证过。你可以直接拿CoolGaN的公开规格书指标去对比验证再根据你实际器件的外延参数做校准。希望这套仿真链路能帮你在功率GaN器件的评估和选型上少走一些弯路。