1. 为什么GaN HEMT电热仿真值得花时间死磕
GaN HEMT这两年在快充、车载OBC、射频功放这些领域火得不行，但凡做功率器件的团队，手里没几个GaN项目都不好意思跟人打招呼。但问题也随之而来：GaN器件功率密度高得离谱，单位面积发热量是硅器件的几倍，自热效应直接决定了他的实际性能和可靠性。你要是只做纯电学仿真，忽略晶格温度对迁移率和饱和速度的影响，结果跟实测能差出30%以上，这个坑我踩过不止一次。
Silvaco Atlas作为老牌TCAD工具，在化合物半导体仿真这块积累很深，尤其是它的Giga模块对GaN材料体系的支持相当完整。但新手拿到Atlas做电热仿真，往往会卡在几个地方：材料参数怎么设、热学边界条件怎么加、物理模型怎么选、代码怎么写才能收敛。网上能找到的教程要么太老，要么只讲电学不讲热，要么代码跑不通。
这篇内容就是把我自己从零搭建GaN HEMT电热仿真流程的经验完整拆开，包括每行代码为什么这么写、参数从哪来、遇到不收敛怎么调。适合刚接触Atlas的研究生、做功率器件的工程师，以及想从硅器件转向宽禁带半导体的朋友。你不需要有很深的TCAD基础，但至少要能看懂基本的器件物理。
2. 仿真整体架构与方案选型思路
2.1 为什么选Atlas而不是其他工具
做GaN HEMT仿真，市面上能选的工具其实不少。Sentaurus功能强但学习曲线陡峭，COMSOL多物理场耦合方便但半导体物理模型不够专业，Silvaco Atlas算是平衡得比较好的选择。它的优势在于：材料数据库里自带GaN、AlGaN的参数，虽然有些需要修正；物理模型库完整，包括极化效应、陷阱效应、热电子发射这些GaN必需的模型；脚本化操作，所有设置都在代码里，方便复现和参数扫描。
我试过用Sentaurus做同样的结构，光是编译材料参数就花了两天，Atlas至少能让你在半天内跑出第一个结果。当然Atlas也有坑，比如默认的GaN参数有些是给光电器件用的，直接拿来做HEMT会出问题，这个后面会细说。
2.2 电热耦合的两种实现路径
Atlas做电热仿真有两条路：一是用Giga模块做完整的晶格热方程求解，二是用热阻网络做近似。前者精度高但计算量大，后者速度快但只能算平均温度。
我的建议是：如果你要看沟道温度分布、热点的具体位置，必须用Giga模块解热方程。如果只是评估不同偏置下的电流退化，热阻网络够用了。这篇内容主要讲第一种，因为新手最容易在热方程这块翻车。
电热耦合的核心逻辑是：电学方程解出功率密度分布，功率密度作为热源代入热方程，热方程解出晶格温度，温度再反馈到电学方程里影响迁移率和饱和速度。这个循环每步都要迭代，直到收敛。Atlas里通过thermcontact和models temperature来激活这个耦合。
2.3 仿真结构的确定
我用的结构是典型的AlGaN/GaN HEMT：衬底用SiC（导热好，适合功率应用），缓冲层GaN厚度2微米，AlGaN势垒层25纳米，Al组分0.25。栅长1微米，源漏间距5微米。这个尺寸不是随便定的，栅长1微米是功率器件的典型值，源漏间距5微米保证击穿电压够用。
结构定义在Atlas里用mesh和region语句完成。新手常犯的错误是网格太粗，尤其是沟道和势垒层界面处，网格必须加密到1纳米以下，否则极化电荷和二维电子气浓度根本算不准。但网格太密又会不收敛，这个平衡点需要试几次才能找到。
3. 代码逐段解读与关键参数设置
3.1 网格定义与区域划分
先看网格部分的代码：
go atlas mesh space.mult=1.0 x.mesh loc=0.0 spac=0.1 x.mesh loc=1.0 spac=0.01 x.mesh loc=2.0 spac=0.05 x.mesh loc=3.0 spac=0.01 x.mesh loc=5.0 spac=0.1 y.mesh loc=0.0 spac=0.01 y.mesh loc=0.025 spac=0.001 y.mesh loc=0.05 spac=0.005 y.mesh loc=0.5 spac=0.05 y.mesh loc=2.0 spac=0.2这里space.mult是全局网格缩放因子，新手建议先用1.0，跑通了再调。x方向在栅极边缘（1微米和3微米处）加密到0.01微米，因为这里电场最集中。y方向在AlGaN势垒层（0到0.025微米）加密到0.001微米，这是二维电子气所在的位置，网格必须足够细才能分辨出量子阱。
注意：y.mesh的第一段从0到0.025用了0.001的间距，意味着这里有25个网格点。如果你算不动，可以先放宽到0.002，但不能再粗了，否则电子浓度会偏低。
区域定义：
region num=1 material=GaN y.min=0.025 y.max=2.0 region num=2 material=AlGaN y.min=0.0 y.max=0.025 x.comp=0.25 region num=3 material=SiC y.min=-1.0 y.max=0.0这里有个细节：AlGaN的x.comp=0.25表示Al组分0.25，但Atlas默认的AlGaN参数是线性插值的，实际HEMT的AlGaN层因为极化效应，参数需要手动修正。我一般会在后面的material语句里覆盖关键参数。
3.2 电极定义与掺杂设置
电极部分：
electrode name=source x.min=0.0 x.max=0.5 y.min=0.0 y.max=0.05 electrode name=drain x.min=4.5 x.max=5.0 y.min=0.0 y.max=0.05 electrode name=gate x.min=1.5 x.max=2.5 y.min=0.0 y.max=0.05源漏电极做成欧姆接触，栅极做成肖特基接触。Atlas里默认电极是欧姆的，栅极需要额外加workfunction参数。我用的Ni/Au栅，功函数设5.1eV。
electrode name=gate workfunction=5.1掺杂方面，GaN缓冲层一般是半绝缘的，掺杂浓度1e15 cm^-3左右。AlGaN势垒层通常不故意掺杂，但会有背景载流子浓度，我设1e16 cm^-3。
doping region=1 uniform conc=1e15 n.type doping region=2 uniform conc=1e16 n.type提示：不要给AlGaN层加太高的掺杂，否则会屏蔽极化效应，二维电子气浓度反而下降。这是新手常犯的错误。
3.3 极化模型与物理模型选择
GaN HEMT的核心是极化效应，没有极化模型就没有二维电子气。Atlas里用polarization语句激活：
models polarization但光激活还不够，还需要设置极化电荷的缩放因子。默认值是1.0，但实际器件的极化电荷会被表面态部分屏蔽，我一般设0.7到0.8。
polarization scale=0.75物理模型方面，必须有的包括：
models srh auger fldmob models temperature=300 models printsrh是肖克利-雷德-霍尔复合，auger是俄歇复合，fldmob是场相关迁移率。温度模型用temperature，后面电热耦合时会改成thermcontact。
迁移率模型我推荐用cvt（Cauchy速度饱和模型），比默认的constant迁移率准确得多：
mobility model=cvt3.4 热学边界条件与热方程设置
这是电热仿真的关键部分。首先激活Giga模块：
models giga然后设置热接触。衬底底部作为热沉，温度固定在300K：
thermcontact name=substrate temperature=300源漏栅电极也可以设热接触，但一般只设衬底就够了。热阻方面，SiC衬底的热导率是490 W/(m·K)，GaN是130 W/(m·K)，这些Atlas材料库里有，但需要确认一下：
material material=GaN tc=1.3 material material=SiC tc=4.9注意单位是W/(cm·K)，所以130 W/(m·K)要写成1.3。
热方程求解需要设置初始温度：
solve init temp=300然后电热耦合求解：
solve b1=0.0 vgate=0.0 vdrain=0.0 solve b1=0.1 vgate=-2.0 vdrain=0.0 solve b1=0.2 vgate=-2.0 vdrain=1.0这里的b1是热方程的求解变量，从0逐步加到1，表示热耦合逐渐增强。直接设b1=1很容易不收敛，必须逐步加。
4. 完整仿真流程与收敛调试实录
4.1 从零到第一个收敛结果的完整步骤
我按实际操作顺序把流程走一遍。第一步，建立结构文件，就是前面那些mesh、region、electrode、doping语句。第二步，设置材料参数，尤其是GaN和AlGaN的迁移率、饱和速度、热导率。第三步，激活物理模型，包括极化、迁移率、复合、热方程。第四步，先做纯电学仿真，不加热方程，确认电学特性合理。第五步，逐步加入热耦合，从b1=0.1开始，每次加0.1，直到b1=1.0。第六步，扫描栅压和漏压，记录电流和温度。
纯电学仿真先跑转移特性：
solve init solve vgate=-3.0 log outf=transfer.log solve vgate=-3.0 vstep=0.1 vfinal=1.0 name=gate如果这一步不收敛，问题多半在网格或极化参数。我遇到过极化缩放因子设1.0导致电子浓度过高、电流不收敛的情况，降到0.75就正常了。
4.2 电热耦合的逐步加载技巧
热耦合加载是新手最容易卡住的地方。直接solve b1=1.0几乎肯定不收敛，因为热源突然加入会导致温度剧烈变化。我的做法是：
solve b1=0.0 solve b1=0.1 solve b1=0.2 ... solve b1=1.0每一步都保存状态，如果某一步不收敛，就减小步长。比如b1=0.5不收敛，就改成0.45、0.48、0.5这样逐步逼近。
还有一个技巧：先用小的漏压做热耦合，比如vdrain=0.1V，等b1=1.0收敛后再加漏压。这样热源小，容易收敛。
4.3 输出特性与温度分布的提取
输出特性扫描：
solve b1=1.0 vgate=0.0 log outf=output.log solve vdrain=0.0 vstep=0.5 vfinal=5.0 name=drain温度分布提取用tonyplot或者直接输出到文件：
output temp save outf=thermal.str然后在TonyPlot里看晶格温度分布。典型结果：栅极漏侧边缘温度最高，因为那里电场最强、功率密度最大。我跑出来的峰值温度在450K左右（环境300K，功耗2W/mm），跟文献报道的吻合。
4.4 关键参数对结果的影响分析
我做了几组对比，发现几个参数影响最大：
| 参数 | 取值 | 峰值温度(K) | 饱和电流(mA/mm) |
|---|---|---|---|
| 极化缩放 | 0.7 | 430 | 850 |
| 极化缩放 | 0.8 | 455 | 920 |
| 极化缩放 | 1.0 | 不收敛 | - |
| 热导率(GaN) | 1.3 | 450 | 900 |
| 热导率(GaN) | 0.8 | 520 | 880 |
| 热导率(GaN) | 1.8 | 410 | 910 |
从表里能看出，极化缩放因子对电流影响很大，热导率对温度影响很大。GaN的热导率实际值跟晶体质量有关，MOCVD外延的GaN热导率一般在1.3到1.5 W/(cm·K)，我建议先用1.3，跟实验对比后再调。
5. 常见报错与排查技巧实录
5.1 不收敛问题的系统排查方法
不收敛是Atlas最常见的报错，没有之一。我的排查顺序是：先看网格，再看模型，再看参数，最后看求解步长。
网格问题：沟道处网格太粗会导致电子浓度震荡。检查方法是用tonyplot看网格分布，如果AlGaN/GaN界面处网格间距大于2纳米，肯定有问题。
模型问题：同时激活太多模型容易冲突。比如polarization和piezoelectric同时用会重复计算极化电荷。建议只用polarization。
参数问题：迁移率模型参数不对会导致电流异常。我遇到过cvt模型的饱和速度设成1e7 cm/s（硅的值），GaN应该是2e7 cm/s。检查material语句里的vsat参数。
求解步长：电压步长太大。把vstep从0.5改成0.1，或者用method语句调整迭代方法：
method newton trapnewton是牛顿迭代，trap是梯形积分，组合起来收敛性更好。
5.2 温度异常与热学参数检查清单
温度异常一般表现为：温度高得离谱（>1000K）或者温度不变化。前者多半是热导率单位搞错了，后者是热方程没激活。
检查清单：
models giga是否激活thermcontact是否设置- 材料热导率单位是否是W/(cm·K)
b1是否从0逐步加到1- 功率密度是否合理（正常在1e5到1e6 W/cm^2量级）
我踩过一个坑：GaN热导率写成130（W/(cm·K)），结果温度只有301K，因为热导率太大了。正确的是1.3。
5.3 收敛后结果不合理的常见原因
有时候收敛了但结果明显不对，比如电流为零、温度分布均匀。电流为零多半是电极没定义对，检查electrode语句的坐标范围是否跟region匹配。温度均匀多半是热源没加对，检查b1是否真的到1.0了。
还有一个隐蔽的问题：solve语句的顺序。必须先解电学再解热学，如果先解热学，热源为零，温度当然不变。正确的顺序是solve init、solve vgate、solve vdrain、solve b1。
5.4 提速技巧：如何缩短仿真时间
电热耦合仿真很吃时间，我跑一个输出特性曲线要半小时。提速方法有几个：一是用method里的itlimit限制迭代次数，二是用save和load保存中间状态，避免重复计算，三是用tonyplot的批处理模式。
method itlimit=20这个设置能防止不收敛时无限迭代。但注意，itlimit太小会导致假收敛，我一般设20到30。
6. 从仿真到实测的校准经验
6.1 用实验数据反推关键参数
仿真跑通只是第一步，结果准不准才是关键。我一般用实验的转移特性曲线反推极化缩放因子和迁移率参数。具体做法：先跑一组不同极化缩放因子的仿真，跟实验的阈值电压和饱和电流对比，选最接近的。然后用输出特性曲线校准热导率和热阻。
我实测的器件阈值电压-2.5V，饱和电流800mA/mm。仿真初始值阈值电压-3.0V，饱和电流950mA/mm。调整极化缩放因子从0.75降到0.7，阈值电压变成-2.6V，饱和电流850mA/mm，接近了。再调迁移率参数，最终匹配到-2.5V和810mA/mm。
6.2 温度校准与热阻提取
温度校准需要实测的沟道温度，一般用拉曼光谱或者红外热像仪测。我没有这些设备，就用文献数据：类似结构的GaN HEMT在2W/mm功耗下沟道温度约450K。我的仿真在同样功耗下是455K，偏差1%，可以接受。
热阻提取：热阻= (T_channel - T_substrate) / P_dissipated。我的仿真里T_channel=455K，T_substrate=300K，P=2W/mm，热阻=77.5 K·mm/W。文献值一般在50到100 K·mm/W，合理。
6.3 仿真指导实验设计的实际案例
有一次实验发现器件在大漏压下电流退化严重，怀疑是自热效应。我用仿真扫了不同漏压下的沟道温度，发现漏压5V时温度达到500K，而500K时迁移率下降30%，正好解释了电流退化。后来实验组改了封装，热阻降低20%，电流退化明显改善。这个案例说明仿真不只是发论文用的，真能指导工程优化。
7. 代码复现的完整清单与参数速查
7.1 可直接运行的完整代码框架
我把关键代码整合成一个可运行的框架，你复制到Atlas里改改参数就能跑：
go atlas mesh space.mult=1.0 x.mesh loc=0.0 spac=0.1 x.mesh loc=1.0 spac=0.01 x.mesh loc=2.0 spac=0.05 x.mesh loc=3.0 spac=0.01 x.mesh loc=5.0 spac=0.1 y.mesh loc=0.0 spac=0.01 y.mesh loc=0.025 spac=0.001 y.mesh loc=0.05 spac=0.005 y.mesh loc=0.5 spac=0.05 y.mesh loc=2.0 spac=0.2 region num=1 material=GaN y.min=0.025 y.max=2.0 region num=2 material=AlGaN y.min=0.0 y.max=0.025 x.comp=0.25 region num=3 material=SiC y.min=-1.0 y.max=0.0 electrode name=source x.min=0.0 x.max=0.5 y.min=0.0 y.max=0.05 electrode name=drain x.min=4.5 x.max=5.0 y.min=0.0 y.max=0.05 electrode name=gate x.min=1.5 x.max=2.5 y.min=0.0 y.max=0.05 workfunction=5.1 doping region=1 uniform conc=1e15 n.type doping region=2 uniform conc=1e16 n.type material material=GaN tc=1.3 vsat=2e7 material material=AlGaN tc=0.3 vsat=1.5e7 material material=SiC tc=4.9 models polarization scale=0.75 models srh auger fldmob cvt giga models temperature=300 thermcontact name=substrate temperature=300 solve init solve vgate=-3.0 solve vgate=-3.0 vstep=0.1 vfinal=1.0 name=gate solve b1=0.0 solve b1=0.1 solve b1=0.2 solve b1=0.3 solve b1=0.4 solve b1=0.5 solve b1=0.6 solve b1=0.7 solve b1=0.8 solve b1=0.9 solve b1=1.0 solve vgate=0.0 log outf=output.log solve vdrain=0.0 vstep=0.5 vfinal=5.0 name=drain output temp save outf=thermal.str7.2 关键参数速查表
| 参数 | 推荐值 | 单位 | 说明 |
|---|---|---|---|
| 极化缩放 | 0.7-0.8 | - | 表面态屏蔽效应 |
| GaN热导率 | 1.3 | W/(cm·K) | MOCVD外延典型值 |
| AlGaN热导率 | 0.3 | W/(cm·K) | 随Al组分变化 |
| SiC热导率 | 4.9 | W/(cm·K) | 单晶典型值 |
| GaN饱和速度 | 2e7 | cm/s | 高场下电子速度 |
| AlGaN饱和速度 | 1.5e7 | cm/s | 比GaN低 |
| 栅功函数 | 5.1 | eV | Ni/Au栅 |
| 缓冲层掺杂 | 1e15 | cm^-3 | 半绝缘 |
| 势垒层掺杂 | 1e16 | cm^-3 | 背景载流子 |
7.3 新手最容易犯的五个错误
第一个错误：极化缩放因子设1.0。这是默认值，但实际器件达不到，会导致电流偏高甚至不收敛。改成0.75。
第二个错误：热导率单位搞错。GaN是1.3 W/(cm·K)，不是130。写130温度几乎不变。
第三个错误：热耦合直接b1=1.0。必须从0逐步加，每次0.1。
第四个错误：网格太粗。AlGaN层至少10个网格点，沟道处间距小于1纳米。
第五个错误：模型冲突。polarization和piezoelectric不要同时用，选一个。
7.4 后续扩展方向
跑通基础电热仿真后，可以往几个方向扩展：一是加陷阱模型，模拟电流崩塌效应；二是加栅极漏电模型，分析肖特基泄漏；三是做瞬态仿真，看热时间常数；四是做可靠性仿真，比如高温反偏。每个方向都需要额外的模型和参数，但基础框架是一样的。
我个人在实际操作中的体会是，GaN HEMT电热仿真最难的不是写代码，而是理解每个参数背后的物理意义。你如果只是抄代码，换个结构就不知道怎么改了。但如果你知道极化缩放因子对应表面态密度、热导率对应晶体质量，调参数就有方向了。我建议新手先跑通一个结构，然后故意改错几个参数，看看结果怎么变，这样理解最深。
