尧图网页设计 
1. GaN HEMT仿真为什么不能照搬硅基器件的那套流程1.1 从一次失败的仿真复现说起几年前我第一次接手GaN HEMT的TCAD仿真任务心态很轻松——之前做硅基LDMOS和FinFET的仿真流程已经很熟了换个材料参数、改改结构尺寸应该跑得通。结果第一版跑出来的转移特性曲线跟实测差了将近一个数量级阈值电压偏了2V以上跨导峰值位置完全对不上。排查了两天才意识到问题出在哪儿我压根没处理极化效应把AlGaN/GaN异质结当成了普通的半导体异质结来处理。这个教训很典型。硅基器件里PN结或者MOS结构的沟道载流子靠掺杂和栅压控制物理图像清晰。但GaN HEMT的核心导电机制完全不同——二维电子气2DEG的来源不是掺杂而是自发极化和压电极化在AlGaN/GaN界面处产生的极化电荷。你不显式地把这层物理加进去仿真器就完全看不到沟道在哪儿。所以这篇内容主要面向两类人一是刚接触GaN器件仿真、有半导体物理基础但没踩过坑的工程师二是做过硅基TCAD仿真、想迁移到宽禁带半导体领域的技术人员。我会围绕Silvaco TCAD这套工具链把极化效应的建模逻辑、物理模型的选择依据、高频特性提取的完整流程拆开来讲尽量做到你照着能复现。1.2 极化效应到底在仿真里扮演什么角色要理解为什么极化效应是GaN HEMT仿真的命门得先搞清楚2DEG是怎么来的。AlGaN和GaN都是纤锌矿结构本身没有中心对称性所以沿c轴方向存在自发极化。当AlGaN势垒层外延生长在GaN缓冲层上时晶格失配还会产生额外的压电极化。这两部分极化强度在异质结界面处不连续形成一个面极化电荷密度σ_pol。这个电荷密度可以到10^13 cm^-2量级对应的面电子密度足以在没有任何掺杂的情况下形成高浓度2DEG。在Silvaco TCAD里处理这个问题的核心思路是在AlGaN/GaN界面处显式定义极化电荷。具体做法有两种路径——一种是用INTERFACE语句直接指定固定电荷密度另一种是用Atlas的极化模型让仿真器根据材料组分自动计算。前者简单粗暴但不够灵活后者更接近物理本质但需要正确设置参数。我个人的经验是如果你的Al组分是固定的、结构也不复杂用固定电荷近似完全够用计算速度快、收敛性好。但如果你要做Al组分渐变或者多沟道结构就必须用极化模型自动计算否则每个界面都要手动调参数工作量巨大且容易出错。1.3 高频特性仿真和直流仿真的本质区别很多人做GaN HEMT仿真只关注直流特性——转移特性、输出特性、跨导。这些当然重要但GaN HEMT真正的战场在射频功率领域高频特性才是核心指标。高频仿真和直流仿真的本质区别在于你必须考虑载流子的输运延迟和电容效应。直流仿真里沟道里的电子被认为能瞬间响应栅压变化但在高频下电子从源端运动到漏端需要时间栅极电容的充放电也需要时间。这些延迟效应直接决定了器件的截止频率fT和最大振荡频率fmax。在Silvaco里做高频仿真核心是小信号AC分析。你需要在直流工作点的基础上叠加一个小幅度的正弦信号然后提取Y参数再转换成h参数或者S参数最后算出fT和fmax。这个过程听起来简单但实际操作中有很多细节——比如AC分析的频率范围怎么设、栅极电阻怎么提取、寄生电容怎么剥离——每一个都会显著影响最终结果。2. Silvaco Atlas中极化模型的参数设置与物理模型选型2.1 材料参数文件里最容易漏掉的几个关键量Silvaco Atlas自带了一套GaN和AlGaN的材料参数但默认值不一定适合你的仿真场景。我建议在material语句里至少显式覆盖以下几个参数电子迁移率模型GaN的体迁移率大概在1000-1200 cm²/V·s但2DEG的迁移率受合金散射、界面粗糙度散射影响实际值可能在1500-2000 cm²/V·s。如果你用的是默认的常数迁移率模型仿真出来的跨导会偏低。饱和速度GaN的电子饱和速度约2.5×10^7 cm/s比硅高出一倍多。这个参数直接影响高频下的电流驱动能力。介电常数GaN的静态介电常数约9.5AlGaN随Al组分变化大概在9.0-9.5之间。这个参数影响电容计算进而影响fT。极化系数包括自发极化系数和压电系数。这些参数在Silvaco的默认材料文件里可能不完整需要手动补充。我一般会在仿真文件开头单独定义一个material区域把所有关键参数集中写在一起方便后续调试和修改。这样做的另一个好处是当你需要切换不同Al组分时只需要改一个地方。2.2 极化电荷的两种建模方式及其适用场景前面提到了两种极化建模方式这里展开说一下具体怎么操作。方式一固定界面电荷这是最简单的做法。在AlGaN/GaN界面处用INTERFACE语句指定电荷密度INTERFACE QF-1.0e13这里的负号表示电子积累。电荷密度的大小需要根据Al组分和应变状态来估算。对于Al组分0.25的AlGaN/GaN异质结极化电荷密度大概在1.0-1.2×10^13 cm^-2。这种方式的优点是收敛快、计算资源消耗少。缺点是如果你改了Al组分或者势垒层厚度需要重新估算电荷密度不够自动化。方式二Atlas内置极化模型Atlas提供了POLARIZATION模型可以自动计算极化电荷。你需要设置MODELS POLARIZATION然后在材料参数里指定压电系数和自发极化系数。这种方式更物理但需要注意几个坑极化模型对网格精度很敏感界面处的网格必须足够密否则电荷分布会不准确。如果你的结构里有钝化层或者栅极绝缘层极化模型可能会在多个界面处都产生电荷需要仔细检查。极化模型的收敛性比固定电荷差可能需要调整数值方法参数。我的建议是初期用固定电荷快速验证结构设计和掺杂分布确认基本电学特性合理后再切换到极化模型做精细仿真。这样能节省大量调试时间。2.3 哪些物理模型是必须开的哪些可以关掉GaN HEMT仿真涉及的物理机制比硅器件复杂但也不是所有模型都要开。下面这张表是我总结的模型开关建议物理模型是否必须原因注意事项极化模型必须决定2DEG来源网格要密收敛性差场依赖迁移率必须高场下速度饱和影响输出特性饱和区SRH复合建议开影响关态泄漏参数需要校准陷阱模型视情况影响电流崩塌和动态特性参数难以准确获取热力学模型视情况大功率下自热效应显著计算量大隧穿模型视情况影响栅极泄漏薄势垒时重要对于初次仿真我建议先开极化模型场依赖迁移率SRH复合跑通直流特性后再逐步加入其他模型。一次性全开容易导致不收敛而且出了问题很难定位是哪个模型引起的。3. 从结构建模到直流特性收敛的完整调试链路3.1 网格划分哪里该密哪里可以粗网格划分是TCAD仿真的基本功但GaN HEMT对网格的要求比硅器件更苛刻。原因是极化电荷只存在于原子级厚度的界面处如果网格太粗电荷会被“抹平”到几个纳米的范围内导致2DEG浓度计算错误。我的经验规则是AlGaN/GaN界面垂直方向网格步长不超过0.5 nm最好在0.2-0.3 nm。势垒层内部垂直步长1-2 nm即可因为势垒层通常只有20-30 nm厚。沟道区域垂直步长1-2 nm水平方向可以根据栅极长度适当放宽。栅极边缘水平方向需要加密因为电场集中在这里。缓冲层可以逐步放宽到10-50 nm因为这里的物理量变化缓慢。一个常见的错误是全局用均匀网格。这样要么界面处不够密要么整体网格量爆炸。正确的做法是用MESH语句配合X.MESH和Y.MESH做非均匀划分。另外提醒一点网格太密也会导致不收敛。因为网格越密方程组的条件数越差数值误差越大。所以不要盲目加密够用就行。3.2 掺杂和极化电荷的配合逻辑GaN HEMT通常是 unintentionally dopedUID结构沟道和缓冲层不需要刻意掺杂。但有几个地方需要注意AlGaN势垒层如果是UID载流子完全靠极化产生。如果做了调制掺杂比如Si delta doping需要显式定义掺杂分布。GaN缓冲层实际外延中会有背景载流子浓度通常在10^15-10^16 cm^-3量级。这个背景浓度会影响缓冲层的泄漏和击穿特性建议在仿真中体现。Cap层很多GaN HEMT表面有一层几纳米厚的GaN cap用于改善欧姆接触和抑制表面态。这层cap会影响表面极化电荷的计算不能忽略。极化电荷和掺杂的关系是极化电荷提供2DEG的来源掺杂影响2DEG的分布和缓冲层行为。两者需要配合设置不能只调一个。3.3 收敛性调试从“仿真发散”到“稳定收敛”的排查顺序“仿真发散”是TCAD用户最常遇到的问题之一。GaN HEMT因为极化模型和非线性迁移率模型的存在收敛难度比硅器件更高。我总结了一套排查顺序第一步检查网格质量。界面处网格是否足够密有没有长宽比过大的网格这些都会导致数值不稳定。第二步检查物理模型组合。是不是一次性开了太多模型试着先关掉陷阱模型和热力学模型只保留极化迁移率SRH看能否收敛。第三步调整数值方法。Atlas提供了多种数值求解方法可以尝试减小初始电压步长METHOD语句里的CLIMIT参数切换到Newton-Raphson方法的阻尼版本增加最大迭代次数第四步检查边界条件。电极的功函数设置是否合理欧姆接触的掺杂是否足够高这些都会影响收敛。第五步分步仿真。先只跑零偏压确认能收敛然后逐步加电压每一步都保存结果。这样即使后面发散了也能知道是从哪一步开始出问题的。我个人的习惯是在正式仿真之前先跑一个极简版本——只有基本结构、只开极化模型、只加一个栅压扫描。确认这个能收敛后再逐步加复杂度。这样能把问题定位到最小的范围。4. 高频小信号仿真从Y参数到fT/fmax的提取细节4.1 小信号AC分析的正确设置方式高频仿真的核心是AC小信号分析。在Silvaco Atlas里你需要先跑一个直流仿真确定工作点通常是最大跨导对应的栅压。在直流工作点基础上用SOLVE语句加上AC参数指定频率范围和步数。仿真器会自动计算Y参数。关键设置包括SOLVE AC FREQ1e9 NFREQ100这表示从1 GHz开始取100个频率点。频率范围的选择很重要——太低比如1 MHz以下看不到电容效应太高比如100 GHz以上可能超出模型适用范围。我一般从1 GHz扫到50 GHz覆盖大多数GaN HEMT的工作频段。还有一个容易忽略的点AC分析的幅度必须足够小通常默认是1V但对于小信号分析来说太大了。建议在SOLVE语句里显式指定ACMAG0.01或者更小确保工作在线性区。4.2 寄生参数剥离什么时候需要怎么做仿真直接提取的Y参数包含了器件的本征部分和寄生部分。如果你关心的是本征fT需要把寄生电阻和寄生电容剥离掉。但如果你只是想看整体性能直接用仿真结果算fT/fmax也可以。剥离的方法通常有两种方法一Cold FET法。在器件完全关断的状态下Vgs远低于阈值Vds0测量S参数提取寄生电容和电阻。这种方法在实测中很常用在仿真中也可以复现。方法二解析法。根据器件的几何尺寸和材料参数直接计算寄生电容栅-源电容、栅-漏电容和寄生电阻源电阻、漏电阻。这种方法精度取决于参数准确性。在仿真中我一般先用解析法估算寄生参数然后在AC仿真结果中手动剥离看看本征fT和整体fT差多少。如果差距不大小于20%说明寄生效应不严重如果差距很大就需要优化器件结构了。4.3 fT和fmax的计算公式及仿真验证截止频率fT的定义是电流增益降到1时的频率。计算公式是fT gm / (2π × (Cgs Cgd))其中gm是跨导Cgs和Cgd分别是栅-源和栅-漏电容。最大振荡频率fmax的定义是功率增益降到1时的频率。计算公式更复杂fmax fT / (2 × sqrt((Rg Rs Ri) / Rds 2π × fT × Rg × Cgd))其中Rg是栅极电阻Rs是源极电阻Ri是沟道充电电阻Rds是输出电阻。在仿真中验证这两个公式的方法是从AC仿真结果中提取Y参数然后转换成h参数h21就是电流增益找到|h21|1对应的频率就是fT。同样从Y参数转换成Mason增益U找到|U|1对应的频率就是fmax。我建议在仿真报告中同时给出两种方法的结果——公式计算值和仿真提取值。如果两者吻合说明参数提取正确如果不吻合需要检查是哪个参数出了问题。5. 极化效应仿真中最容易踩的五个坑5.1 坑一把极化电荷当成掺杂来处理这是新手最常犯的错误。极化电荷是界面电荷不是体掺杂。如果你在AlGaN层里加了一层高浓度掺杂来“模拟”2DEG虽然电学特性看起来差不多但物理图像完全错了。这样做会导致阈值电压对势垒层厚度的依赖关系不对高频电容特性不对无法正确反映极化效应随温度的变化正确的做法是在界面处定义电荷而不是在体材料里掺杂。5.2 坑二忽略表面态对2DEG的耗尽实际GaN HEMT的表面存在大量陷阱态这些表面态会俘获电子耗尽部分2DEG导致实际阈值电压比理论值偏正。如果你在仿真中不包含表面态算出来的阈值电压会偏负跨导也会偏大。在Silvaco里可以通过在表面定义陷阱来实现。陷阱的能级和密度需要根据实验数据校准没有通用值。我一般先用文献里的典型值比如Ec-0.5eV密度10^12-10^13 cm^-2然后根据实测的阈值电压微调。5.3 坑三高频仿真时忘记更新工作点AC小信号分析必须在一个确定的直流工作点上进行。如果你改了栅压或者漏压但没有重新跑直流仿真AC结果就是错的。我见过有人直接在一个旧的直流结果上跑AC结果fT算出来差了30%。正确的流程是每次改变偏置条件后先跑直流仿真确认工作点再跑AC仿真。在Silvaco的输入文件里这意味着SOLVE语句的顺序不能乱。5.4 坑四迁移率模型参数直接照搬硅器件GaN的迁移率特性跟硅完全不同。硅的迁移率随电场增加而下降的规律比如Caughey-Thomas模型在GaN里需要重新校准参数。如果你直接用硅的参数高场下的速度饱和行为会算错导致输出特性的饱和区电流偏差很大。我建议从文献里找GaN专用的迁移率模型参数或者用实测的I-V曲线反推。如果实在找不到至少要把饱和速度改成GaN的值2.5×10^7 cm/s这个参数对高频特性影响最大。5.5 坑五不检查电荷守恒极化模型计算完成后一定要检查整个器件的电荷守恒。具体做法是在平衡状态下把所有电荷包括极化电荷、掺杂离子、自由载流子积分看总和是否接近零。如果偏差很大说明极化电荷设置有问题。这个检查在Silvaco里可以通过输出电荷分布来实现。虽然麻烦一点但能避免很多后续的诡异问题。6. 从仿真到实测的校准思路让数据对得上6.1 先校准直流再校准射频仿真和实测对不上是常态关键是要有系统的校准流程。我的建议是第一步校准阈值电压。阈值电压主要受势垒层厚度、Al组分、表面态影响。先调这几个参数让仿真阈值电压跟实测一致。第二步校准最大跨导。跨导主要受迁移率和栅极电容影响。如果阈值电压已经对了但跨导不对优先调迁移率参数。第三步校准输出电阻。输出电阻反映沟道长度调制效应和缓冲层泄漏。如果输出特性曲线太平或者太斜需要调缓冲层的背景掺杂或者陷阱参数。第四步校准fT和fmax。直流特性都对了之后再跑AC仿真。如果fT对不上检查栅极电容和跨导是否准确如果fmax对不上检查寄生电阻是否合理。这个流程的核心逻辑是直流特性是基础射频特性是上层建筑。直流不对射频肯定不对直流对了射频不一定对但至少有了调试的基础。6.2 参数提取中的不确定度处理TCAD仿真里有很多参数是没法直接测量的比如陷阱密度、界面态密度、极化系数。这些参数的不确定度会传递到最终结果里。我的处理方式是对关键参数做敏感性分析。比如把陷阱密度在合理范围内变化±50%看阈值电压变化多少。如果变化很小说明这个参数不关键可以用典型值如果变化很大就需要想办法通过实验数据来约束。另一个技巧是用多个实验数据点同时约束。比如同时用转移特性和电容-电压特性来校准比只用一条曲线更可靠。6.3 什么时候该放弃仿真直接做实验仿真不是万能的。有些问题仿真很难回答比如陷阱的具体能级和密度分布表面态的具体来源和化学性质器件退化机制这些问题需要实验手段比如DLTS、XPS、TEM来回答。仿真的价值在于在实验之前筛选方案在实验之后解释现象。如果你发现仿真调了很久都对不上实测可能不是仿真参数的问题而是物理模型本身不适用。这时候该做的是回到实验获取更多信息而不是继续在仿真里死磕。7. 写在最后一些个人体会做GaN HEMT仿真这些年最大的感受是极化效应不是“一个选项”而是“整个仿真是否成立的前提”。你可以不开热力学模型、不开陷阱模型但极化模型必须开而且必须开对。我见过太多人在这上面栽跟头包括我自己。另一个体会是高频仿真比直流仿真更需要“手感”。直流仿真里参数和结果的关系比较直观——迁移率大了跨导就大阈值电压偏了调势垒层厚度。但高频仿真里fT和fmax受多个参数耦合影响有时候调一个参数fT升了但fmax降了。这时候需要理解背后的物理机制而不是盲目试参数。最后分享一个小技巧在仿真文件里加注释。把每个参数的来源、每个模型的开关理由、每次调试的改动都写清楚。GaN HEMT仿真文件动辄几百行过一个月再回来看没有注释根本记不住当时为什么这么设。这个习惯看起来不起眼但能省下大量重复调试的时间。 GaN HEMT仿真避坑指南：极化效应建模与高频特性提取
1. GaN HEMT仿真为什么不能照搬硅基器件的那套流程1.1 从一次失败的仿真复现说起几年前我第一次接手GaN HEMT的TCAD仿真任务心态很轻松——之前做硅基LDMOS和FinFET的仿真流程已经很熟了换个材料参数、改改结构尺寸应该跑得通。结果第一版跑出来的转移特性曲线跟实测差了将近一个数量级阈值电压偏了2V以上跨导峰值位置完全对不上。排查了两天才意识到问题出在哪儿我压根没处理极化效应把AlGaN/GaN异质结当成了普通的半导体异质结来处理。这个教训很典型。硅基器件里PN结或者MOS结构的沟道载流子靠掺杂和栅压控制物理图像清晰。但GaN HEMT的核心导电机制完全不同——二维电子气2DEG的来源不是掺杂而是自发极化和压电极化在AlGaN/GaN界面处产生的极化电荷。你不显式地把这层物理加进去仿真器就完全看不到沟道在哪儿。所以这篇内容主要面向两类人一是刚接触GaN器件仿真、有半导体物理基础但没踩过坑的工程师二是做过硅基TCAD仿真、想迁移到宽禁带半导体领域的技术人员。我会围绕Silvaco TCAD这套工具链把极化效应的建模逻辑、物理模型的选择依据、高频特性提取的完整流程拆开来讲尽量做到你照着能复现。1.2 极化效应到底在仿真里扮演什么角色要理解为什么极化效应是GaN HEMT仿真的命门得先搞清楚2DEG是怎么来的。AlGaN和GaN都是纤锌矿结构本身没有中心对称性所以沿c轴方向存在自发极化。当AlGaN势垒层外延生长在GaN缓冲层上时晶格失配还会产生额外的压电极化。这两部分极化强度在异质结界面处不连续形成一个面极化电荷密度σ_pol。这个电荷密度可以到10^13 cm^-2量级对应的面电子密度足以在没有任何掺杂的情况下形成高浓度2DEG。在Silvaco TCAD里处理这个问题的核心思路是在AlGaN/GaN界面处显式定义极化电荷。具体做法有两种路径——一种是用INTERFACE语句直接指定固定电荷密度另一种是用Atlas的极化模型让仿真器根据材料组分自动计算。前者简单粗暴但不够灵活后者更接近物理本质但需要正确设置参数。我个人的经验是如果你的Al组分是固定的、结构也不复杂用固定电荷近似完全够用计算速度快、收敛性好。但如果你要做Al组分渐变或者多沟道结构就必须用极化模型自动计算否则每个界面都要手动调参数工作量巨大且容易出错。1.3 高频特性仿真和直流仿真的本质区别很多人做GaN HEMT仿真只关注直流特性——转移特性、输出特性、跨导。这些当然重要但GaN HEMT真正的战场在射频功率领域高频特性才是核心指标。高频仿真和直流仿真的本质区别在于你必须考虑载流子的输运延迟和电容效应。直流仿真里沟道里的电子被认为能瞬间响应栅压变化但在高频下电子从源端运动到漏端需要时间栅极电容的充放电也需要时间。这些延迟效应直接决定了器件的截止频率fT和最大振荡频率fmax。在Silvaco里做高频仿真核心是小信号AC分析。你需要在直流工作点的基础上叠加一个小幅度的正弦信号然后提取Y参数再转换成h参数或者S参数最后算出fT和fmax。这个过程听起来简单但实际操作中有很多细节——比如AC分析的频率范围怎么设、栅极电阻怎么提取、寄生电容怎么剥离——每一个都会显著影响最终结果。2. Silvaco Atlas中极化模型的参数设置与物理模型选型2.1 材料参数文件里最容易漏掉的几个关键量Silvaco Atlas自带了一套GaN和AlGaN的材料参数但默认值不一定适合你的仿真场景。我建议在material语句里至少显式覆盖以下几个参数电子迁移率模型GaN的体迁移率大概在1000-1200 cm²/V·s但2DEG的迁移率受合金散射、界面粗糙度散射影响实际值可能在1500-2000 cm²/V·s。如果你用的是默认的常数迁移率模型仿真出来的跨导会偏低。饱和速度GaN的电子饱和速度约2.5×10^7 cm/s比硅高出一倍多。这个参数直接影响高频下的电流驱动能力。介电常数GaN的静态介电常数约9.5AlGaN随Al组分变化大概在9.0-9.5之间。这个参数影响电容计算进而影响fT。极化系数包括自发极化系数和压电系数。这些参数在Silvaco的默认材料文件里可能不完整需要手动补充。我一般会在仿真文件开头单独定义一个material区域把所有关键参数集中写在一起方便后续调试和修改。这样做的另一个好处是当你需要切换不同Al组分时只需要改一个地方。2.2 极化电荷的两种建模方式及其适用场景前面提到了两种极化建模方式这里展开说一下具体怎么操作。方式一固定界面电荷这是最简单的做法。在AlGaN/GaN界面处用INTERFACE语句指定电荷密度INTERFACE QF-1.0e13这里的负号表示电子积累。电荷密度的大小需要根据Al组分和应变状态来估算。对于Al组分0.25的AlGaN/GaN异质结极化电荷密度大概在1.0-1.2×10^13 cm^-2。这种方式的优点是收敛快、计算资源消耗少。缺点是如果你改了Al组分或者势垒层厚度需要重新估算电荷密度不够自动化。方式二Atlas内置极化模型Atlas提供了POLARIZATION模型可以自动计算极化电荷。你需要设置MODELS POLARIZATION然后在材料参数里指定压电系数和自发极化系数。这种方式更物理但需要注意几个坑极化模型对网格精度很敏感界面处的网格必须足够密否则电荷分布会不准确。如果你的结构里有钝化层或者栅极绝缘层极化模型可能会在多个界面处都产生电荷需要仔细检查。极化模型的收敛性比固定电荷差可能需要调整数值方法参数。我的建议是初期用固定电荷快速验证结构设计和掺杂分布确认基本电学特性合理后再切换到极化模型做精细仿真。这样能节省大量调试时间。2.3 哪些物理模型是必须开的哪些可以关掉GaN HEMT仿真涉及的物理机制比硅器件复杂但也不是所有模型都要开。下面这张表是我总结的模型开关建议物理模型是否必须原因注意事项极化模型必须决定2DEG来源网格要密收敛性差场依赖迁移率必须高场下速度饱和影响输出特性饱和区SRH复合建议开影响关态泄漏参数需要校准陷阱模型视情况影响电流崩塌和动态特性参数难以准确获取热力学模型视情况大功率下自热效应显著计算量大隧穿模型视情况影响栅极泄漏薄势垒时重要对于初次仿真我建议先开极化模型场依赖迁移率SRH复合跑通直流特性后再逐步加入其他模型。一次性全开容易导致不收敛而且出了问题很难定位是哪个模型引起的。3. 从结构建模到直流特性收敛的完整调试链路3.1 网格划分哪里该密哪里可以粗网格划分是TCAD仿真的基本功但GaN HEMT对网格的要求比硅器件更苛刻。原因是极化电荷只存在于原子级厚度的界面处如果网格太粗电荷会被“抹平”到几个纳米的范围内导致2DEG浓度计算错误。我的经验规则是AlGaN/GaN界面垂直方向网格步长不超过0.5 nm最好在0.2-0.3 nm。势垒层内部垂直步长1-2 nm即可因为势垒层通常只有20-30 nm厚。沟道区域垂直步长1-2 nm水平方向可以根据栅极长度适当放宽。栅极边缘水平方向需要加密因为电场集中在这里。缓冲层可以逐步放宽到10-50 nm因为这里的物理量变化缓慢。一个常见的错误是全局用均匀网格。这样要么界面处不够密要么整体网格量爆炸。正确的做法是用MESH语句配合X.MESH和Y.MESH做非均匀划分。另外提醒一点网格太密也会导致不收敛。因为网格越密方程组的条件数越差数值误差越大。所以不要盲目加密够用就行。3.2 掺杂和极化电荷的配合逻辑GaN HEMT通常是 unintentionally dopedUID结构沟道和缓冲层不需要刻意掺杂。但有几个地方需要注意AlGaN势垒层如果是UID载流子完全靠极化产生。如果做了调制掺杂比如Si delta doping需要显式定义掺杂分布。GaN缓冲层实际外延中会有背景载流子浓度通常在10^15-10^16 cm^-3量级。这个背景浓度会影响缓冲层的泄漏和击穿特性建议在仿真中体现。Cap层很多GaN HEMT表面有一层几纳米厚的GaN cap用于改善欧姆接触和抑制表面态。这层cap会影响表面极化电荷的计算不能忽略。极化电荷和掺杂的关系是极化电荷提供2DEG的来源掺杂影响2DEG的分布和缓冲层行为。两者需要配合设置不能只调一个。3.3 收敛性调试从“仿真发散”到“稳定收敛”的排查顺序“仿真发散”是TCAD用户最常遇到的问题之一。GaN HEMT因为极化模型和非线性迁移率模型的存在收敛难度比硅器件更高。我总结了一套排查顺序第一步检查网格质量。界面处网格是否足够密有没有长宽比过大的网格这些都会导致数值不稳定。第二步检查物理模型组合。是不是一次性开了太多模型试着先关掉陷阱模型和热力学模型只保留极化迁移率SRH看能否收敛。第三步调整数值方法。Atlas提供了多种数值求解方法可以尝试减小初始电压步长METHOD语句里的CLIMIT参数切换到Newton-Raphson方法的阻尼版本增加最大迭代次数第四步检查边界条件。电极的功函数设置是否合理欧姆接触的掺杂是否足够高这些都会影响收敛。第五步分步仿真。先只跑零偏压确认能收敛然后逐步加电压每一步都保存结果。这样即使后面发散了也能知道是从哪一步开始出问题的。我个人的习惯是在正式仿真之前先跑一个极简版本——只有基本结构、只开极化模型、只加一个栅压扫描。确认这个能收敛后再逐步加复杂度。这样能把问题定位到最小的范围。4. 高频小信号仿真从Y参数到fT/fmax的提取细节4.1 小信号AC分析的正确设置方式高频仿真的核心是AC小信号分析。在Silvaco Atlas里你需要先跑一个直流仿真确定工作点通常是最大跨导对应的栅压。在直流工作点基础上用SOLVE语句加上AC参数指定频率范围和步数。仿真器会自动计算Y参数。关键设置包括SOLVE AC FREQ1e9 NFREQ100这表示从1 GHz开始取100个频率点。频率范围的选择很重要——太低比如1 MHz以下看不到电容效应太高比如100 GHz以上可能超出模型适用范围。我一般从1 GHz扫到50 GHz覆盖大多数GaN HEMT的工作频段。还有一个容易忽略的点AC分析的幅度必须足够小通常默认是1V但对于小信号分析来说太大了。建议在SOLVE语句里显式指定ACMAG0.01或者更小确保工作在线性区。4.2 寄生参数剥离什么时候需要怎么做仿真直接提取的Y参数包含了器件的本征部分和寄生部分。如果你关心的是本征fT需要把寄生电阻和寄生电容剥离掉。但如果你只是想看整体性能直接用仿真结果算fT/fmax也可以。剥离的方法通常有两种方法一Cold FET法。在器件完全关断的状态下Vgs远低于阈值Vds0测量S参数提取寄生电容和电阻。这种方法在实测中很常用在仿真中也可以复现。方法二解析法。根据器件的几何尺寸和材料参数直接计算寄生电容栅-源电容、栅-漏电容和寄生电阻源电阻、漏电阻。这种方法精度取决于参数准确性。在仿真中我一般先用解析法估算寄生参数然后在AC仿真结果中手动剥离看看本征fT和整体fT差多少。如果差距不大小于20%说明寄生效应不严重如果差距很大就需要优化器件结构了。4.3 fT和fmax的计算公式及仿真验证截止频率fT的定义是电流增益降到1时的频率。计算公式是fT gm / (2π × (Cgs Cgd))其中gm是跨导Cgs和Cgd分别是栅-源和栅-漏电容。最大振荡频率fmax的定义是功率增益降到1时的频率。计算公式更复杂fmax fT / (2 × sqrt((Rg Rs Ri) / Rds 2π × fT × Rg × Cgd))其中Rg是栅极电阻Rs是源极电阻Ri是沟道充电电阻Rds是输出电阻。在仿真中验证这两个公式的方法是从AC仿真结果中提取Y参数然后转换成h参数h21就是电流增益找到|h21|1对应的频率就是fT。同样从Y参数转换成Mason增益U找到|U|1对应的频率就是fmax。我建议在仿真报告中同时给出两种方法的结果——公式计算值和仿真提取值。如果两者吻合说明参数提取正确如果不吻合需要检查是哪个参数出了问题。5. 极化效应仿真中最容易踩的五个坑5.1 坑一把极化电荷当成掺杂来处理这是新手最常犯的错误。极化电荷是界面电荷不是体掺杂。如果你在AlGaN层里加了一层高浓度掺杂来“模拟”2DEG虽然电学特性看起来差不多但物理图像完全错了。这样做会导致阈值电压对势垒层厚度的依赖关系不对高频电容特性不对无法正确反映极化效应随温度的变化正确的做法是在界面处定义电荷而不是在体材料里掺杂。5.2 坑二忽略表面态对2DEG的耗尽实际GaN HEMT的表面存在大量陷阱态这些表面态会俘获电子耗尽部分2DEG导致实际阈值电压比理论值偏正。如果你在仿真中不包含表面态算出来的阈值电压会偏负跨导也会偏大。在Silvaco里可以通过在表面定义陷阱来实现。陷阱的能级和密度需要根据实验数据校准没有通用值。我一般先用文献里的典型值比如Ec-0.5eV密度10^12-10^13 cm^-2然后根据实测的阈值电压微调。5.3 坑三高频仿真时忘记更新工作点AC小信号分析必须在一个确定的直流工作点上进行。如果你改了栅压或者漏压但没有重新跑直流仿真AC结果就是错的。我见过有人直接在一个旧的直流结果上跑AC结果fT算出来差了30%。正确的流程是每次改变偏置条件后先跑直流仿真确认工作点再跑AC仿真。在Silvaco的输入文件里这意味着SOLVE语句的顺序不能乱。5.4 坑四迁移率模型参数直接照搬硅器件GaN的迁移率特性跟硅完全不同。硅的迁移率随电场增加而下降的规律比如Caughey-Thomas模型在GaN里需要重新校准参数。如果你直接用硅的参数高场下的速度饱和行为会算错导致输出特性的饱和区电流偏差很大。我建议从文献里找GaN专用的迁移率模型参数或者用实测的I-V曲线反推。如果实在找不到至少要把饱和速度改成GaN的值2.5×10^7 cm/s这个参数对高频特性影响最大。5.5 坑五不检查电荷守恒极化模型计算完成后一定要检查整个器件的电荷守恒。具体做法是在平衡状态下把所有电荷包括极化电荷、掺杂离子、自由载流子积分看总和是否接近零。如果偏差很大说明极化电荷设置有问题。这个检查在Silvaco里可以通过输出电荷分布来实现。虽然麻烦一点但能避免很多后续的诡异问题。6. 从仿真到实测的校准思路让数据对得上6.1 先校准直流再校准射频仿真和实测对不上是常态关键是要有系统的校准流程。我的建议是第一步校准阈值电压。阈值电压主要受势垒层厚度、Al组分、表面态影响。先调这几个参数让仿真阈值电压跟实测一致。第二步校准最大跨导。跨导主要受迁移率和栅极电容影响。如果阈值电压已经对了但跨导不对优先调迁移率参数。第三步校准输出电阻。输出电阻反映沟道长度调制效应和缓冲层泄漏。如果输出特性曲线太平或者太斜需要调缓冲层的背景掺杂或者陷阱参数。第四步校准fT和fmax。直流特性都对了之后再跑AC仿真。如果fT对不上检查栅极电容和跨导是否准确如果fmax对不上检查寄生电阻是否合理。这个流程的核心逻辑是直流特性是基础射频特性是上层建筑。直流不对射频肯定不对直流对了射频不一定对但至少有了调试的基础。6.2 参数提取中的不确定度处理TCAD仿真里有很多参数是没法直接测量的比如陷阱密度、界面态密度、极化系数。这些参数的不确定度会传递到最终结果里。我的处理方式是对关键参数做敏感性分析。比如把陷阱密度在合理范围内变化±50%看阈值电压变化多少。如果变化很小说明这个参数不关键可以用典型值如果变化很大就需要想办法通过实验数据来约束。另一个技巧是用多个实验数据点同时约束。比如同时用转移特性和电容-电压特性来校准比只用一条曲线更可靠。6.3 什么时候该放弃仿真直接做实验仿真不是万能的。有些问题仿真很难回答比如陷阱的具体能级和密度分布表面态的具体来源和化学性质器件退化机制这些问题需要实验手段比如DLTS、XPS、TEM来回答。仿真的价值在于在实验之前筛选方案在实验之后解释现象。如果你发现仿真调了很久都对不上实测可能不是仿真参数的问题而是物理模型本身不适用。这时候该做的是回到实验获取更多信息而不是继续在仿真里死磕。7. 写在最后一些个人体会做GaN HEMT仿真这些年最大的感受是极化效应不是“一个选项”而是“整个仿真是否成立的前提”。你可以不开热力学模型、不开陷阱模型但极化模型必须开而且必须开对。我见过太多人在这上面栽跟头包括我自己。另一个体会是高频仿真比直流仿真更需要“手感”。直流仿真里参数和结果的关系比较直观——迁移率大了跨导就大阈值电压偏了调势垒层厚度。但高频仿真里fT和fmax受多个参数耦合影响有时候调一个参数fT升了但fmax降了。这时候需要理解背后的物理机制而不是盲目试参数。最后分享一个小技巧在仿真文件里加注释。把每个参数的来源、每个模型的开关理由、每次调试的改动都写清楚。GaN HEMT仿真文件动辄几百行过一个月再回来看没有注释根本记不住当时为什么这么设。这个习惯看起来不起眼但能省下大量重复调试的时间。