企业官网建设流程全解析

1. GaN HEMT仿真探索的背景与整体思路

1.1 为什么选择GaN HEMT作为研究对象

功率半导体器件这几年变化很快，硅基器件在高压高频场景下逐渐逼近材料极限，宽禁带半导体开始走到台前。GaN（氮化镓）和SiC（碳化硅）是两条主流路线，其中GaN HEMT（高电子迁移率晶体管）凭借高电子迁移率、高击穿电场和高饱和漂移速度，在射频功率放大和高效电源转换领域占据了一个很特殊的位置。

我最初接触GaN HEMT是因为一个电源项目，需要在几百千赫兹到兆赫兹的开关频率下保持高效率，同时耐压要覆盖几百伏。硅基MOSFET在这个频段开关损耗已经很难压下去，SiC MOSFET虽然耐压够，但高频特性不如GaN。于是转向GaN HEMT，而器件选型和结构优化离不开TCAD仿真，Silvaco Atlas就是这条路上的主力工具。

这里说的“探索”，不是漫无目的地试参数，而是带着明确目标去理解器件内部的物理过程：二维电子气（2DEG）是怎么形成的、极化效应贡献了多少、栅极漏电从哪里来、缓冲层陷阱如何影响动态导通电阻。这些问题如果只靠看论文，理解始终隔一层，只有自己在Atlas里把结构搭出来、把物理模型打开、把曲线跑出来，才能形成真正的直觉。

1.2 Silvaco Atlas在GaN器件仿真中的定位

Silvaco Atlas是一个基于物理的二维/三维器件仿真器，核心是求解泊松方程和载流子连续性方程，配合漂移-扩散模型或者更高级的能带模型。对于GaN HEMT，Atlas的优势在于它内置了极化模型，可以处理AlGaN/GaN异质结界面处的自发极化和压电极化，这是GaN HEMT区别于普通MOSFET的关键。

很多人第一次用Atlas做GaN仿真会踩一个坑：直接照搬硅器件的结构定义方式，结果发现根本看不到2DEG。原因在于GaN HEMT的导电沟道不是靠掺杂形成的，而是靠极化效应在异质结界面处感应出来的。如果不开启极化模型，或者极化参数设置不对，仿真出来的器件就是一块“死”的半导体，电流几乎为零。

Atlas里和GaN相关的关键模块包括：

• 极化模型：通过POLARIZATION语句开启，需要指定CALC.STRAIN和相关的压电常数。
• 陷阱模型：GaN缓冲层和表面态对器件特性影响极大，需要用TRAPS语句定义。
• 迁移率模型：GaN的电子迁移率在高场下会饱和，需要开启FLDMOB。
• 热模型：大功率工作时自热效应显著，需要耦合GIGA或者HEAT模块。

这些模块不是全部打开就好，每开一个都会增加收敛难度。我的经验是分阶段来：先跑通基本结构，确认2DEG浓度合理，再逐步加陷阱、加自热、加击穿模型。

1.3 整体仿真流程设计

一个完整的GaN HEMT仿真流程，我通常分成五步走：

1. 结构定义：确定衬底材料（Si、SiC或蓝宝石）、缓冲层厚度、GaN沟道厚度、AlGaN势垒层厚度和Al组分、栅极位置和尺寸。
2. 掺杂与极化设置：GaN HEMT通常不需要刻意掺杂沟道，但缓冲层可能需要弱p型或深能级补偿。极化参数必须和材料组分匹配。
3. 电极定义与接触设置：源漏欧姆接触、栅极肖特基接触，接触功函数和界面态需要合理设置。
4. 物理模型选择：根据仿真目标选择漂移-扩散、能量平衡或者蒙特卡洛模型。
5. 求解与后处理：先跑平衡态，再跑I-V曲线，最后提取转移特性、输出特性、跨导、击穿电压等参数。

这个流程看起来线性，实际操作中经常需要回头调整。比如跑I-V时发现电流太小，可能是极化电荷算错了；发现跨导异常，可能是栅极漏电没控制好。每一次回头都是一次对器件物理的重新理解。

2. Atlas中GaN HEMT结构搭建的核心细节

2.1 材料参数与极化模型设置

GaN HEMT的核心结构是AlGaN/GaN异质结。在Atlas中定义材料时，不能只用默认的GaN参数，需要手动设置关键材料属性。以下是我常用的一套参数，供参考：

极化模型的设置是重中之重。Atlas中通过以下语句开启：

POLARIZATION CALC.STRAIN

然后需要在材料定义中指定压电常数和弹性常数。对于AlGaN，其自发极化系数和压电系数随Al组分变化，Atlas支持用ALCOMP参数来插值。如果Al组分设成0.25，但极化参数用的是AlN的值，算出来的2DEG浓度会偏高很多。

我实测过一组对比：Al组分0.25、势垒层25nm时，极化模型正确设置的情况下，2DEG面密度大约在1e13 cm⁻²量级；如果忘记开极化模型，2DEG浓度直接掉到1e10以下，基本等于没有沟道。

注意：Atlas的极化模型对网格精度很敏感。异质结界面处的网格必须足够密，通常建议在界面上下各2nm范围内至少有5个网格点，否则极化电荷的积分会出错。

2.2 网格划分策略与收敛性权衡

网格划分是Atlas仿真中最考验经验的一环。网格太粗，物理量变化剧烈的区域算不准；网格太密，矩阵规模爆炸，收敛时间成倍增加。GaN HEMT的网格需要重点关注三个区域：

• 异质结界面：2DEG就分布在这里，网格必须密。
• 栅极边缘：电场集中，尤其是栅漏之间的峰值电场。
• 沟道下方缓冲层：如果考虑缓冲层陷阱和漏电，这里也需要适当加密。

我的做法是先用较粗的网格跑通结构，确认能收敛后，再在关键区域逐步加密。Atlas支持X.MESH和Y.MESH语句手动定义网格线，也可以用QGRID自动加密。对于GaN HEMT，我通常手动定义界面处的网格，因为自动加密有时候会在不需要的地方浪费节点。

一个典型的网格设置片段：

X.MESH LOC=0.0 SPACING=0.05 X.MESH LOC=0.5 SPACING=0.02 X.MESH LOC=1.0 SPACING=0.05 Y.MESH LOC=0.0 SPACING=0.01 Y.MESH LOC=0.02 SPACING=0.001 Y.MESH LOC=0.025 SPACING=0.0005 Y.MESH LOC=0.05 SPACING=0.002

这里Y方向在0.02到0.025微米之间是异质结界面，网格间距设到0.5nm，保证极化电荷的精确积分。

2.3 电极定义与接触特性处理

GaN HEMT的源漏接触通常是欧姆接触，栅极是肖特基接触。Atlas中定义电极时，需要指定接触类型和功函数。

源漏欧姆接触相对简单，用CONTACT语句指定NAME和WORKFUNCTION即可。但实际GaN器件的欧姆接触电阻并不为零，如果忽略接触电阻，仿真出来的导通电阻会偏小。我通常会在源漏区域加一个薄的重掺杂层来模拟接触电阻，或者直接在电极上设置RESISTANCE参数。

栅极肖特基接触是难点。肖特基势垒高度直接影响栅极漏电和阈值电压。对于Ni/Au栅极，势垒高度大约在0.8到1.0 eV之间。Atlas中设置：

CONTACT NAME=gate WORKFUNCTION=4.9

这里4.9 eV对应Ni的功函数。但实际器件的栅极特性还受界面态影响，如果仿真中栅极漏电和实验对不上，可以在栅极下方加一层薄绝缘层，或者用INTERFACE语句定义界面态。

实操心得：栅极漏电对网格和模型非常敏感。如果发现栅极电流随栅压指数上升的斜率不对，先检查栅极下方的网格是否足够密，再检查是否开启了隧穿模型。Atlas中的TUNNEL模型可以模拟Fowler-Nordheim隧穿和直接隧穿，但开启后收敛会变难，建议先用较简单的模型跑通，再逐步加复杂度。

3. 物理模型选择与参数标定实操

3.1 迁移率模型与速度饱和效应

GaN的电子迁移率在低场下可以到1200 cm²/V·s以上，但在高场下会饱和，饱和速度约2.5e7 cm/s。Atlas中默认的迁移率模型是常数迁移率，对于GaN HEMT远远不够。必须开启场依赖迁移率模型：

MODEL FLDMOB

这个模型会根据电场强度调整迁移率，模拟速度饱和效应。但FLDMOB模型本身也有多种参数，比如VSAT（饱和速度）、BETAN（拟合参数），需要根据材料调整。

我对比过不同迁移率模型对输出特性的影响：只开常数迁移率时，饱和电流偏高，输出曲线在高漏压下的平坦度不够；开启FLDMOB后，饱和电流下降约15%，输出曲线的平坦度明显改善，更接近实验数据。

如果仿真目标包含击穿特性，还需要开启碰撞电离模型：

IMPACT SELB

GaN的碰撞电离系数比硅低，击穿电场更高。Atlas中的SELB模型需要设置A和B参数，对于GaN，常用的值大约是A=1.0e8 cm⁻¹，B=1.5e7 V/cm。这些参数不是固定的，不同文献给出的值有差异，需要根据实验数据标定。

3.2 陷阱模型与动态特性模拟

GaN HEMT最让人头疼的问题之一是电流崩塌（current collapse），也就是动态导通电阻增大。这个现象的根源是缓冲层陷阱和表面态。Atlas中可以用TRAPS语句定义陷阱：

TRAPS REGION=1 NUM=2 TRAP TYPE=ACCEPTOR LEVEL=0.6 DENSITY=1e17 CROSS=1e-15 TRAP TYPE=DONOR LEVEL=0.4 DENSITY=5e16 CROSS=1e-14

这里定义了两个陷阱：一个受主陷阱，能级距导带0.6 eV，浓度1e17 cm⁻³；一个施主陷阱，能级0.4 eV，浓度5e16 cm⁻³。这些参数需要根据DLTS测试或者瞬态电流测试来标定。

陷阱模型对收敛性影响很大。我的经验是：先不加陷阱跑通直流特性，确认基本结构没问题后，再逐步加陷阱。加陷阱时先加一个，跑通后再加第二个。如果一次性加太多，收敛失败的概率很高。

注意：陷阱的捕获截面（CROSS）参数也很关键。如果截面设得太大，陷阱充放电时间常数会很小，瞬态仿真时步长需要设得很小，否则会错过关键变化。我通常先用1e-15 cm²量级的截面跑，再根据瞬态响应调整。

3.3 自热效应与热网络耦合

GaN HEMT在大功率工作时自热效应显著，沟道温度可以比环境温度高几十甚至上百摄氏度。Atlas中可以用GIGA模块或者HEAT模块来耦合热效应。

GIGA模块是Atlas自带的非等温仿真模块，可以同时求解电学和热学方程。使用时需要定义热导率、热容等参数：

MATERIAL MATERIAL=GaN TC.CONST=1.3 MATERIAL MATERIAL=GaN HEAT.CAP=3.0

这里TC.CONST是热导率，单位W/(cm·K)，GaN大约1.3；HEAT.CAP是热容，单位J/(cm³·K)，GaN大约3.0。

自热效应会显著降低饱和电流，因为迁移率随温度升高而下降。我做过对比：同样的结构，不考虑自热时饱和电流约1.2 A/mm，考虑自热后降到0.9 A/mm左右，降幅约25%。这个差异在功率应用中不能忽略。

如果仿真目标只是直流特性，自热模块可以不开，用等温仿真加快速度。但如果要评估功率器件的热稳定性，自热模块必须开。

4. 仿真结果分析与常见问题排查

4.1 转移特性与阈值电压提取

转移特性是GaN HEMT最基本的特性曲线。Atlas中通过扫描栅压，固定漏压，记录漏极电流：

SOLVE VGATE=0.0 VSTEP=-0.1 VFINAL=-5.0 NAME=gate

这里从0V扫到-5V，步长-0.1V。GaN HEMT通常是耗尽型，阈值电压为负。提取阈值电压的方法有多种，我常用的是线性外推法：在转移特性线性区做切线，切线与栅压轴的交点即为阈值电压。

实测中常见的问题是阈值电压偏移。如果仿真得到的阈值电压比预期偏负很多，可能的原因包括：

• 极化电荷计算偏大，导致2DEG浓度偏高，需要更负的栅压才能耗尽。
• 栅极功函数设置不当，肖特基势垒偏低。
• AlGaN势垒层厚度或Al组分设置错误。

排查时可以先检查2DEG浓度：在平衡态下提取异质结界面处的电子浓度，如果超过1.5e13 cm⁻²，说明极化电荷偏大，需要检查Al组分和应变参数。

4.2 输出特性与导通电阻分析

输出特性是漏压扫描，固定栅压，记录漏极电流：

SOLVE VDRAIN=0.0 VSTEP=0.5 VFINAL=10.0 NAME=drain

GaN HEMT的输出特性通常分为线性区和饱和区。线性区的斜率对应导通电阻，饱和区的平坦度反映输出电阻。

导通电阻的组成包括：源漏接触电阻、沟道电阻、漂移区电阻。仿真中如果导通电阻偏大，可能的原因有：

• 源漏接触电阻设置过大。
• 沟道迁移率偏低，检查迁移率模型参数。
• 栅极长度过长，沟道电阻增加。

我遇到过一次导通电阻异常偏大的情况，排查后发现是源漏电极的RESISTANCE参数设成了1e-3 Ω·cm²，而实际应该是1e-5 Ω·cm²量级。改过来后导通电阻恢复正常。

4.3 常见收敛问题与排查速查表

Atlas仿真GaN HEMT时，收敛问题是最常见的障碍。以下是我整理的一份速查表：

实操心得：Atlas的收敛问题很多时候不是模型本身的问题，而是网格和步长的问题。我习惯在METHOD语句里设置ITLIMIT=50和DELTA=0.1，给求解器更多迭代次数和更小的初始步长。如果还是不行，就用AUTO语句让Atlas自动调整。

4.4 仿真与实验偏差的标定思路

仿真和实验有偏差是正常的，关键是如何系统性地缩小偏差。我的标定思路是：

1. 先标转移特性：调整栅极功函数和AlGaN参数，让阈值电压和跨导峰值对上。
2. 再标输出特性：调整迁移率模型参数和接触电阻，让导通电阻和饱和电流对上。
3. 最后标动态特性：调整陷阱参数，让瞬态响应和电流崩塌对上。

每一步只调一个参数，调完记录，避免多个参数同时改导致无法追溯。这个过程可能需要反复迭代，但每次迭代都会让你对器件物理的理解更深一层。

5. 从仿真到设计迭代的经验总结

5.1 结构参数对性能的影响规律

经过多轮仿真，我总结了一些GaN HEMT结构参数对性能的影响规律：

• AlGaN势垒层厚度：增厚势垒层可以增加2DEG浓度，但也会增加栅极到沟道的距离，降低跨导。存在一个最优厚度，通常在20到30nm之间。
• Al组分：提高Al组分可以增加极化电荷，提高2DEG浓度，但也会增加合金散射，降低迁移率。Al组分0.2到0.3是常见范围。
• 栅极长度：缩短栅极长度可以提高跨导和截止频率，但也会加剧短沟道效应。对于射频应用，栅极长度通常在0.1到0.5微米。
• 缓冲层厚度：增厚缓冲层可以降低衬底漏电，但也会增加热阻。对于Si衬底，缓冲层通常在1到3微米。

这些规律不是绝对的，具体器件需要具体优化。但有了这些规律，至少知道往哪个方向调。

5.2 仿真效率提升的实用技巧

Atlas仿真GaN HEMT通常比较耗时，一个完整的I-V扫描可能需要几十分钟甚至几个小时。以下是我常用的效率提升技巧：

• 分阶段仿真：先跑粗网格，确认趋势后再跑细网格。粗网格的结果虽然精度不够，但趋势通常是对的。
• 保存中间状态：用SAVE语句保存平衡态和关键偏置点的状态，后续仿真可以从这些状态继续，不用从头跑。
• 并行计算：Atlas支持多线程，在METHOD语句里设置THREADS=4可以显著加快速度。
• 合理设置电压步长：步长太大容易发散，步长太小浪费时间。我通常先用0.5V步长跑一遍，确认没问题后再用0.1V步长跑精细曲线。

5.3 后续扩展方向

GaN HEMT仿真还有很多可以深入的方向。比如：

• 射频特性仿真：用小信号AC分析提取截止频率和最大振荡频率。
• 击穿特性仿真：用碰撞电离模型和热模型耦合，评估击穿电压和雪崩耐量。
• 可靠性仿真：用陷阱模型和热模型耦合，评估长期工作下的参数退化。
• 电路级仿真：把器件仿真结果提取成紧凑模型，导入SPICE做电路级仿真。

这些方向我还在陆续探索，后续有新的心得再分享。仿真这件事，说到底是用物理模型去逼近真实器件，每一次和实验的对齐都是一次对物理理解的修正。踩过的坑越多，对器件的直觉就越准。