做功率器件这行几乎避不开英飞凌CoolGaN。数据手册翻来覆去就那些数导通电阻、BV击穿电压、Qg栅极电荷、开关损耗。但真正动手去优化器件结构、评估下一代方案时光有数据手册完全不够——场板为什么要这么设计p-GaN掺多少才能把阈值电压顶到1V以上缓冲层里那点碳掺杂浓度对动态损耗影响有多大这些东西用万用表和示波器永远看不到。Silvaco TCAD能干的事就是把AlGaN/GaN HEMT一层层建出来把极化电荷、电场分布、载流子输运这些底层物理摊开看。这篇文章我从Silvaco TCAD的安装写起完整走一遍CoolGaN这类增强型GaN功率器件的建模、BV评估、Qg提取和开关损耗仿真流程顺便把仿真里最容易翻车的几个地方也一并说了。1. 为什么用Silvaco TCAD拆解CoolGaN这类器件1.1 数据手册给的是结果TCAD给的是过程任何一颗商用GaN功率器件到工程师手里时已经被封装好、测试完我们能接触到的只是最终的外特性。数据手册告诉你BV是650V但不会告诉你这个耐压是靠哪一层电场尖峰撑住的手册告诉你Qg是5.8 nC但不会告诉你Miller平台附近的Cgd非线性有多严重手册告诉你开关损耗是X微焦耳但不会告诉你在哪个电流密度下、哪个温度下开始恶化。TCAD仿真恰恰补上这一层。器件结构可以通过TCAD重建载流子浓度、能带图、电场分布这些都是仿真里的增量输出。这意味着你可以做几件实际测试很难做的事单独看某一层的作用。比如把缓冲层掺杂浓度改一个数量级观察对BV和动态Ron的影响。把失效过程“放慢”。击穿前那一刻的电场尖峰在哪、碰撞电离率哪里最高能看到导致失效的物理根因。做结构优化。在流片之前用TCAD筛选场板尺寸、AlGaN厚度、p-GaN掺杂浓度等参数。1.2 CoolGaN的通用层结构与TCAD抽象思路英飞凌的CoolGaN有很多具体型号但从器件物理角度讲它的核心是增强型E-modep-GaN栅极HEMT。TCAD建模不需要100%复刻厂家的工艺细节只需要抓住对电学特性影响最大的几何和掺杂参数。一个典型的CoolGaN类器件从下往上长这样层材料典型厚度在仿真中的作用衬底高阻硅/SiC数百微米机械支撑参与衬底偏置效应成核层AlN100~200 nm释放应力异质外延过渡缓冲层GaN常掺C2~4 μm纵向耐压、泄露和动态Ron的关键势垒层AlGaN15~25 nm提供极化电荷形成2DEG栅极层p-GaNMg掺杂70~150 nm耗尽栅下2DEG实现E-mode钝化SiN100~300 nm表面钝化、场板介质TCAD仿真里我一般保留衬底、GaN缓冲层、AlGaN势垒层、p-GaN这几个核心区域把成核层简化掉或者合并到缓冲层底部。原因很简单成核层的细节对器件表面附近的电场分布影响很小但对网格数量和收敛时间影响很大第一版仿真没必要为它买单。材料参数上GaN的禁带宽度取3.4 eV介电常数9.5左右AlGaN的禁带宽度按铝组分线性插值再加个弓形修正。Al组分通常取0.2~0.25这个值直接决定极化电荷量和2DEG密度。2. Silvaco TCAD安装与结构建模从环境准备到2DEG出现2.1 环境准备和License常见问题Silvaco TCAD到手后的第一道槛通常是环境。不是软件本身难装而是License和版本兼容容易把人卡住。先说下载安装的基本流程从Silvaco官网申请试用License或者用公司正版授权拿到license文件后把SFLMSilvaco Floating License Manager起起来设置好环境变量然后启动DeckBuild。Linux环境是主力Windows虽然有版本但很多仿真脚本、数据处理还是在Linux下跑得最稳。最常遇到的几个问题License起不来。先ping一下License服务器再确认端口有没有被防火墙挡。环境变量LM_LICENSE_FILE或SFLM_LICENSE要指向实际license路径。命令行启动后找不到可执行文件。检查PATH里有没有包含Atlas、TonyPlot这些工具的bin目录。DeckBuild能打开但仿真报错“no such feature”。多半是license feature不全找管理员确认授权了哪个模块。GaN HEMT仿真比较依赖ATLAS的Blaze模块。安装完成后用一个最简单的二极管结构跑通一遍确认环境没问题再往上搭GaN HEMT。千万不要一上来就建全结构不然出问题时分不清是自己脚本错了还是环境有问题。一个实用习惯把deckbuild -run -outfile或者命令行提交跑成批量脚本让仿真结果和日志自动归档。跑参数扫描时这样效率高得多。2.2 用ATLAS语句搭建GaN HEMT结构Silvaco TCAD做器件仿真核心工具是ATLAS。我用它搭GaN HEMT结构时基本思路是先定义网格再定义区域和电极然后掺杂、设置材料和物理模型最后加偏置扫描。下面是一个简化的ATLAS输入结构展示搭建增强型p-GaN HEMT骨架的方式go atlas # 网格定义横向x方向纵向y方向 mesh x.mesh0.05 space.multi1.0 mesh y.mesh0.002 space.multi20 mesh x.mesh0.10 space.multi20 # 区域从上到下分别是p-GaN、AlGaN、GaN缓冲层 region number1 materialGaN y.min0.44 y.max3.00 region number2 materialAlGaN y.min0.27 y.max0.44 region number3 materialGaN y.min0.15 y.max0.27 # 电极 electrode namesource x.min0.0 x.max1.0 y.min0.10 y.max0.15 electrode namegate x.min4.0 x.max6.0 y.min0.00 y.max0.15 electrode namedrain x.min18.0 x.max20.0 y.min0.10 y.max0.15这里的一个关键决策是网格密度。AlGaN/GaN异质结界面附近的网格间距我通常压到0.002 μm量级因为2DEG就聚集在这个界面上纵向网格太粗会直接把载流子浓度分布抹平。横向方向栅极边缘和场板边缘必须加密电场尖峰就在那里。掺杂方面缓冲层p型掺杂浓度1e16 cm⁻³量级有时要加C补偿陷阱。AlGaN势垒层不故意掺杂极化电荷会在界面自然诱导n型导电。p-GaN栅极Mg掺杂浓度在1e17~1e19 cm⁻³这个范围调。2.3 极化电荷与p-GaN耗尽让E-mode成立的建模关键GaN HEMT和普通MOSFET最大的区别在于导电沟道不是靠栅压“感应”出来的而是AlGaN/GaN界面的极化效应天然形成了一层高浓度的二维电子气。TCAD建模必须把这个物理机制表达出来否则器件静态特性就对不上。极化来源有两个自发极化spontaneous polarization和压电极化piezoelectric polarization。AlGaN势垒层受到GaN缓冲层的晶格失配应力两种极化在AlGaN/GaN界面形成正的净固定电荷吸引电子聚集形成2DEG。在Silvaco TCAD里通常用一个界面电荷来描述这个效应典型值在 1e13 cm⁻² 量级。这个值直接决定导通电流和阈值电压是校准的第一步。对于E-mode器件p-GaN栅极层的建模更讲究。p-GaN和栅极金属形成的结在零偏压下会耗尽栅极下方的2DEG让器件默认关断。仿真里需要给p-GaN区域设置合适的p型掺杂浓度并且用金属功函数work function控制栅极接触。p-GaN掺杂浓度越高阈值电压越高但栅极漏电也越难模拟这是需要反复迭代的一组参数。建完结构,第一件事不是急着扫BV而是先跑输出特性Id-Vg和Id-Vd确认阈值电压在合理范围增强型器件通常0.5~2V最大饱和电流量级大致对得上目标关断态漏电没有异常偏高。如果这些基础特性都离谱后面BV和Qg的结果没有参考价值。3. 击穿电压BV仿真网格、碰撞电离和缓冲层陷阱的联合作用3.1 漏极扫描与击穿判据BV仿真的电路接法很简单源极接地栅极加0V或负压让器件处于关断态漏极电压从0V逐步往上扫描观察漏极漏电流。操作上有两个细节容易忽视第一扫描步进。粗扫可以每10V一步但接近预计击穿电压时一定要加密到1V甚至0.5V一步不然击穿点捕捉不准。我习惯用动态步进声明让Silvaco根据电流变化率自动调整步长。第二击穿判据。所谓“击穿电压”不是某个物理量一跳变就认定的通常看漏极电流是否超过某个标准比如1 μA/μm 或 1 mA/mm。不同文献和产品标准不一样写文章或做报告时一定要注明用的什么判据否则数据没有可比性。3.2 碰撞电离模型与网格收敛性BV仿真必须打开碰撞电离模型。Silvaco TCAD里的常用选择是Selberherr模型impact selb这个模型用载流子的碰撞电离系数来描述雪崩击穿。GaN材料和硅的碰撞电离参数差异很大需要单独设置GaN的电子和空穴碰撞电离系数。很多初学仿真一上来直接默认参数结果BV比实际值偏高一倍以上大概率就是这里错了。碰撞电离是强非线性过程对网格非常敏感。栅极靠近漏极一侧的电场尖峰区域横向网格至少要加密到0.02 μm纵向在AlGaN/GaN界面附近0.005 μm以内。判断网格是否收敛有个笨但有效的办法把关键区域网格密度做2倍加密如果BV变化在5%以内基本可以认为网格足够。变化还很大的话继续加密。3.3 缓冲层陷阱和衬底偏置对BV的显著影响CoolGaN这类GaN on Si器件的BV并不是单纯看AlGaN/GaN异质结的纵向耐压。缓冲层的质量很大程度上决定器件能扛多高的电压。工艺上缓冲层会掺碳C或铁Fe来补偿残留施主但这些深能级陷阱在高压下会电离影响纵向电场分布和动态导通电阻。TCAD仿真的处理方式是在缓冲层中间加一个深受主陷阱能级能级位置通常在价带顶上方0.5~0.9 eV的位置浓度根据材料质量在 1e16~1e18 cm⁻³ 之间调节。这个参数对BV和动态Ron影响非常大但它在数据手册上不存在只能靠TCAD加量测反馈来校准。衬底偏置效应也值得专门提一句。GaN on Si器件的衬底是导电的高阻硅衬底电位实际会影响背面电场分布。很多TCAD仿真图省事把衬底悬空或直接忽略结果低压时问题不大高压时漏电路径和实际不符。正确做法是给衬底电极加明确的偏置电位比如接地或接源极和实际应用一致。3.4 场板在仿真里的建模作用CoolGaN这类650V器件的耐压设计很大程度靠场板来缓解栅极边缘的电场集中。TCAD里可以看到没有场板时电场峰值像一把尖刀加上场板后电场峰值被摊开成两个或者多个峰整体面积不变但峰值高度下降BV因此提高。仿真时场板就是把电极图形延长覆盖到钝化层上再接对应的电位。栅场板接栅电位源场板接源电位。效果上源场板通常对动态性能改善更明显栅场板对静态BV更直接。模拟时可以分别扫描栅场板长度和源场板长度看BV和Cgd反馈电容的变化趋势。往往BV提高1%对应的场板长度增加会拖累开关速度好几个百分点这就是器件级优化的取舍。4. 栅极电荷Qg提取瞬态电流法和AC电容积分法4.1 为什么Qg不能“直接读”出来Qg是电路级的参数不是器件某个稳态下的天然输出。它代表把栅极从关断电压充到开通电压、完成整个开通过程所需要的总电荷量。TCAD仿真的原始输出里没有“Qg”这个变量需要自己构造测试条件再积分。数据手册的Qg通常在指定测试条件比如从Vgs0到10V、Vds400V下测得的。你不能随便选一个偏置点说这就是Qg。仿真想做得好就要高度还原这个测试条件。4.2 恒定栅极电流法提取Qg这是最直观的一种方法。思路是给栅极一个恒定的充电电流持续充电同时监测Vgs随时间的变化。因为电流恒定栅极电压平台区Miller平台持续的时间越长说明反馈电荷越多。TCAD里的做法是设计一个瞬态仿真先把器件偏置到目标漏极电压比如Vds400V栅极处于关断电位。从外部接一个恒流源到栅极电流大小选在1mA量级太大太小都会影响平台分辨力。瞬态仿真跑起来记录Vgs(t)。把Vgs从初始值充到目标值比如10V所需的时间乘上恒定电流就是总Qg。Vgs(t)波形在Miller平台那一段会变得非常平对应Cgd剧烈变化、栅极电流几乎都用来给反馈电容充电的阶段。平台的长短直接对应Qg的重要组成。把Ig固定这个条件想清楚QgIg×t整个提取过程就是一次单位换算。用这个方法好处是物理含义清楚和数据手册的测试原理一致适合做产品对标。缺点是对数值噪声比较敏感尤其是平台区域稍微有点振荡就会让时间读数不准所以时间步长要控制得足够小。4.3 AC小信号法提取Ciss/Crss再积分第二条路线是先把电容曲线求出来再积分。这个方法更接近电路工程师的直觉Qg其实是Ciss和Crss在充放电过程中的积分结果。在Silvaco ATLAS里做AC小信号分析在栅极加一个1MHz左右的小信号激励在不同Vgs和Vds偏置下提取Y参数再换算成Cgs、Cgd、Cds。CissCgsCgdCrssCgd。把电容从关断态偏置扫描到开通态偏置对Vgs积分就得到Qg。这个方法能额外看到Cgd的非线性变化——Miller平台区间的本质就是Cgd从很高的值急剧下降到很低的值。有了电容曲线你不仅能算Qg还能诊断它主要是由哪一段电压区间贡献的。这对仿真驱动优化很有用如果发现Cgd在Vds接近0时特别大说明场板设计和钝化层厚度还有调整空间。4.4 Qg读数的实操注意事项仿真Qg往往比数据手册偏低一些这是正常的主要原因没有算封装寄生电感的贡献没有完整模拟栅极Pad和内部引线的寄生电容p-GaN栅极层中的受主态、施主态会贡献额外的电荷。所以对标数据手册时重点看趋势不要盯着绝对值差10个点就怀疑自己模型。重点观察几个趋势对不对Vds越高Qg越大还是越小、Vgs平台电压随温度怎么走、场板加长后Qg涨多少。这些趋势对了模型就是可信的。5. 开关损耗仿真混合模式下的双脉冲测试复现5.1 单器件仿真到混合模式仿真的跨越前面的BV和Qg都可以在纯器件仿真里完成因为偏置条件相对简单。开关损耗就不同了它天然是一个电路行为栅极驱动器、负载电感、续流二极管都在共同作用单器件瞬态仿真很难真实还原。Silvaco TCAD提供了混合模式MixedMode仿真能力把ATLAS器件模型和SPICE电路元件搭在同一个网表里。GaN HEMT的器件部分用TCAD描述外部的栅极电阻、负载电感、续流二极管、电压源用电路元件描述。这样既能保留器件的物理细节又能模拟真实电路工况。5.2 双脉冲测试的电路搭建开关损耗的行业标准测量方法是双脉冲测试Double Pulse TestDPT。TCAD复现DPT的思路一样直流母线电压VDD设为典型值比如400V负载电感串在母线和器件漏极之间电感值选到能在第一脉冲期间把电流升到目标值续流二极管并联在电感两端给感性电流提供续流路径栅极驱动器由一个脉冲电压源串联栅极电阻构成输出两段式栅极脉冲。第一段脉冲把电感电流充到目标电流然后关断在死区时间内电流通过续流二极管“续流”。第二段短脉冲打开器件让器件在目标电流和母线电压同时存在的情况下完成开关动作。真正的导通损耗和关断损耗就是在第二段脉冲开通时刻和第一段脉冲关断时刻测的。5.3 从波形中计算Eon和Eoff仿真跑完直接从波形处理软件里导出Vds(t)、Id(t)然后做乘积积分开通损耗Eon ∫Vds×Id dt积分区间从栅源电压上升到目标值的10%开始到Vds下降到目标值的10%结束关断损耗Eoff ∫Vds×Id dt积分区间从Vds上升到目标值的10%开始到Id下降到目标值的10%结束这个积分定义不是唯一的IEEE标准里给出的定义和各家产品数据手册的定义有细微差异。仿真报告里一定要写清楚用了哪个定义否则别人复现不出同样的数。TCAD混合模式仿真的一个现实问题2D仿真的电流是按单位宽度计算的。你建的是一个截面电流结果默认是单位宽度1 μm的电流。如果仿真目标是一个10 mm宽的大功率器件所有电流、电荷、损耗都要乘上这个宽度因子。Silvaco里可以用area参数来设置有效面积或宽度让外部电路的电流和实际器件匹配。这一步忘记处理开关损耗数据会差好几个数量级属于最常见的初级错误。另外续流二极管最好用真实的TCAD器件模型而不是理想二极管模型这看目标。如果是纯做GaN器件本身的开关损耗参数优化用理想续流二极管可以排除外部变量问题更清楚。如果想复现整个模块的实际损耗那续流管也要建模尤其要看它的反向恢复行为。但GaN肖特基二极管基本没有反向恢复对于很多应用理想模型已经能给出足够参考。6. 参数校准与常见坑让仿真贴近实测数据6.1 最影响结果的三类参数TCAD建模做到最后考验的不是会不会操作性而是会不会校准。一套未经校准的GaN HEMT仿真即使定性趋势对定量指标也可能差得很远。根据我自己的经验下面三类参数对结果影响最大参数类型具体参数主要影响的指标极化相关界面电荷密度、AlGaN Al组分、AlGaN厚度2DEG密度、Id、Vth输运相关迁移率模型参数、饱和速度Id-Vg、Ron、开关电流陷阱相关缓冲层陷阱能级与浓度、表面陷阱BV、动态Ron、关断漏电极化相关参数决定“有没有沟道、有多少电流”是最底层的东西。它的标定依据通常是霍尔测试或者C-V测试得到的2DEG面密度。界面电荷密度微调一个5%Id就可能差10%。输运参数决定“电流能跑多快”。GaN的峰值迁移率和饱和速度和硅差很多低温、高温下的变化趋势也要注意。迁移率模型里那个峰值迁移率参数建议用同批次晶圆的TLM传输线法数据来标。陷阱参数是最难搞的因为它对静态特性影响不明显但对动态特性影响巨大。缓冲层陷阱主要影响高压下的漏电和动态Ron退化。表面陷阱对电流崩塌影响很大。TCAD里陷阱的参数有能级位置、俘获截面、浓度三个维度排列组合起来非常让人头疼。实际操作中我只挑一到两个主要陷阱能级来拟合再加一个表面态密度作为调节项不然参数自由度太多收敛反而没有意义。6.2 我踩过的几个坑网格依赖问题。第一次做BV仿真出来的结果比量测高20%而且怎么调陷阱参数都不对最后发现是漏极边缘网格不够密电场峰值被网格光滑效应抹掉了。BV仿真的网格绝对不是“差不多就行”一定要做网格收敛性检查。直接照抄论文材料参数。GaN材料参数在不同文献里差异很大尤其是碰撞电离系数和饱和速度。同一套结构用A论文的参数和B论文的参数BV能差50V。最好的办法是找到最接近自己工艺平台的那篇文献先从它的参数开始再用自己手里的实测数据修正。忽略了温度边界条件。功率器件仿真自热效应self-heating不能一直关着。TCAD默认的等温条件反映的是短脉冲测量结果连续工作或大电流下必须打开热边界模型让晶格温度随功率升高。否则输出特性在高压大电流段会偏乐观开关损耗趋势也会判断错方向。用2D仿真直接对标3D封装器件的绝对数值。这个问题在前面的宽度因子那里提过实际使用中太容易忽视了。不光是电流还包括寄生电容、电感效应这些在2D TCAD里根本没有。如果要对标封装后的产品TCAD结果只适合校准芯片物理不适合直接和外部的系统级波形做数值对比。6.3 我建议的参数校准流程最后分享一套我自己在用的校准流程思路是先从稳态到瞬态从低压到高压逐级逼近用转移特性Id-Vg调界面电荷密度和p-GaN掺杂浓度先把阈值电压和峰值电流对上。用输出特性Id-Vd调迁移率模型和源漏接触电阻把线性区和饱和区形态校准好。用电容-电压C-V数据调AlGaN厚度和钝化层厚度确保Cgs/Cgd的量级正确。用关断态漏电调缓冲层陷阱浓度用BV调碰撞电离系数和缓冲层厚度。用QG实测值做最终验证必要时微调栅极区域的网格和陷阱参数。每一步只动一到两个参数不要试图一次拟合所有数据。TCAD仿真搭结构可能只用半天真正费时间的全在校准这一步。但校准做好了后面再做结构优化、电场分析、温度外推结果才有参考价值否则仿真做得再花哨也只是画了一堆漂亮的错误曲线。我用这套流程跑过不止一个增强型GaN方案从结构建模到BV、Qg、开关损耗三个维度的评估整体下来的经验是TCAD仿真最值钱的部分从来不是那张结构截图而是你能用它快速回答“如果把这里改一下器件关键指标会怎么变”。流片很贵仿真便宜这正是器件工程师值得在这套工具上投入时间的原因。