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 HEMTAtlas的优势在于它内置了极化模型可以处理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仿真流程我通常分成五步走结构定义确定衬底材料Si、SiC或蓝宝石、缓冲层厚度、GaN沟道厚度、AlGaN势垒层厚度和Al组分、栅极位置和尺寸。掺杂与极化设置GaN HEMT通常不需要刻意掺杂沟道但缓冲层可能需要弱p型或深能级补偿。极化参数必须和材料组分匹配。电极定义与接触设置源漏欧姆接触、栅极肖特基接触接触功函数和界面态需要合理设置。物理模型选择根据仿真目标选择漂移-扩散、能量平衡或者蒙特卡洛模型。求解与后处理先跑平衡态再跑I-V曲线最后提取转移特性、输出特性、跨导、击穿电压等参数。这个流程看起来线性实际操作中经常需要回头调整。比如跑I-V时发现电流太小可能是极化电荷算错了发现跨导异常可能是栅极漏电没控制好。每一次回头都是一次对器件物理的重新理解。2. Atlas中GaN HEMT结构搭建的核心细节2.1 材料参数与极化模型设置GaN HEMT的核心结构是AlGaN/GaN异质结。在Atlas中定义材料时不能只用默认的GaN参数需要手动设置关键材料属性。以下是我常用的一套参数供参考参数GaNAlNAl0.25Ga0.75N禁带宽度 (eV)3.426.24.1电子亲和势 (eV)4.11.93.5相对介电常数9.58.59.2电子迁移率 (cm²/V·s)1200300800饱和速度 (cm/s)2.5e71.5e72.0e7极化模型的设置是重中之重。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 LOC0.0 SPACING0.05 X.MESH LOC0.5 SPACING0.02 X.MESH LOC1.0 SPACING0.05 Y.MESH LOC0.0 SPACING0.01 Y.MESH LOC0.02 SPACING0.001 Y.MESH LOC0.025 SPACING0.0005 Y.MESH LOC0.05 SPACING0.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 NAMEgate WORKFUNCTION4.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 SELBGaN的碰撞电离系数比硅低击穿电场更高。Atlas中的SELB模型需要设置A和B参数对于GaN常用的值大约是A1.0e8 cm⁻¹B1.5e7 V/cm。这些参数不是固定的不同文献给出的值有差异需要根据实验数据标定。3.2 陷阱模型与动态特性模拟GaN HEMT最让人头疼的问题之一是电流崩塌current collapse也就是动态导通电阻增大。这个现象的根源是缓冲层陷阱和表面态。Atlas中可以用TRAPS语句定义陷阱TRAPS REGION1 NUM2 TRAP TYPEACCEPTOR LEVEL0.6 DENSITY1e17 CROSS1e-15 TRAP TYPEDONOR LEVEL0.4 DENSITY5e16 CROSS1e-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 MATERIALGaN TC.CONST1.3 MATERIAL MATERIALGaN HEAT.CAP3.0这里TC.CONST是热导率单位W/(cm·K)GaN大约1.3HEAT.CAP是热容单位J/(cm³·K)GaN大约3.0。自热效应会显著降低饱和电流因为迁移率随温度升高而下降。我做过对比同样的结构不考虑自热时饱和电流约1.2 A/mm考虑自热后降到0.9 A/mm左右降幅约25%。这个差异在功率应用中不能忽略。如果仿真目标只是直流特性自热模块可以不开用等温仿真加快速度。但如果要评估功率器件的热稳定性自热模块必须开。4. 仿真结果分析与常见问题排查4.1 转移特性与阈值电压提取转移特性是GaN HEMT最基本的特性曲线。Atlas中通过扫描栅压固定漏压记录漏极电流SOLVE VGATE0.0 VSTEP-0.1 VFINAL-5.0 NAMEgate这里从0V扫到-5V步长-0.1V。GaN HEMT通常是耗尽型阈值电压为负。提取阈值电压的方法有多种我常用的是线性外推法在转移特性线性区做切线切线与栅压轴的交点即为阈值电压。实测中常见的问题是阈值电压偏移。如果仿真得到的阈值电压比预期偏负很多可能的原因包括极化电荷计算偏大导致2DEG浓度偏高需要更负的栅压才能耗尽。栅极功函数设置不当肖特基势垒偏低。AlGaN势垒层厚度或Al组分设置错误。排查时可以先检查2DEG浓度在平衡态下提取异质结界面处的电子浓度如果超过1.5e13 cm⁻²说明极化电荷偏大需要检查Al组分和应变参数。4.2 输出特性与导通电阻分析输出特性是漏压扫描固定栅压记录漏极电流SOLVE VDRAIN0.0 VSTEP0.5 VFINAL10.0 NAMEdrainGaN HEMT的输出特性通常分为线性区和饱和区。线性区的斜率对应导通电阻饱和区的平坦度反映输出电阻。导通电阻的组成包括源漏接触电阻、沟道电阻、漂移区电阻。仿真中如果导通电阻偏大可能的原因有源漏接触电阻设置过大。沟道迁移率偏低检查迁移率模型参数。栅极长度过长沟道电阻增加。我遇到过一次导通电阻异常偏大的情况排查后发现是源漏电极的RESISTANCE参数设成了1e-3 Ω·cm²而实际应该是1e-5 Ω·cm²量级。改过来后导通电阻恢复正常。4.3 常见收敛问题与排查速查表Atlas仿真GaN HEMT时收敛问题是最常见的障碍。以下是我整理的一份速查表问题现象可能原因解决方法初始平衡态不收敛极化电荷过大网格太粗减小极化电荷加密界面网格栅压扫描时发散栅极漏电模型太激进关闭隧穿模型先跑通再开漏压扫描时电流突跳碰撞电离模型参数不当调整SELB参数减小电压步长瞬态仿真不收敛陷阱时间常数太小增大捕获截面减小时间步长自热仿真温度异常热导率或热容设置错误检查材料热参数确认单位实操心得Atlas的收敛问题很多时候不是模型本身的问题而是网格和步长的问题。我习惯在METHOD语句里设置ITLIMIT50和DELTA0.1给求解器更多迭代次数和更小的初始步长。如果还是不行就用AUTO语句让Atlas自动调整。4.4 仿真与实验偏差的标定思路仿真和实验有偏差是正常的关键是如何系统性地缩小偏差。我的标定思路是先标转移特性调整栅极功函数和AlGaN参数让阈值电压和跨导峰值对上。再标输出特性调整迁移率模型参数和接触电阻让导通电阻和饱和电流对上。最后标动态特性调整陷阱参数让瞬态响应和电流崩塌对上。每一步只调一个参数调完记录避免多个参数同时改导致无法追溯。这个过程可能需要反复迭代但每次迭代都会让你对器件物理的理解更深一层。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语句里设置THREADS4可以显著加快速度。合理设置电压步长步长太大容易发散步长太小浪费时间。我通常先用0.5V步长跑一遍确认没问题后再用0.1V步长跑精细曲线。5.3 后续扩展方向GaN HEMT仿真还有很多可以深入的方向。比如射频特性仿真用小信号AC分析提取截止频率和最大振荡频率。击穿特性仿真用碰撞电离模型和热模型耦合评估击穿电压和雪崩耐量。可靠性仿真用陷阱模型和热模型耦合评估长期工作下的参数退化。电路级仿真把器件仿真结果提取成紧凑模型导入SPICE做电路级仿真。这些方向我还在陆续探索后续有新的心得再分享。仿真这件事说到底是用物理模型去逼近真实器件每一次和实验的对齐都是一次对物理理解的修正。踩过的坑越多对器件的直觉就越准。
GaN HEMT仿真探索：基于Silvaco Atlas的极化模型与陷阱效应实操
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 HEMTAtlas的优势在于它内置了极化模型可以处理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仿真流程我通常分成五步走结构定义确定衬底材料Si、SiC或蓝宝石、缓冲层厚度、GaN沟道厚度、AlGaN势垒层厚度和Al组分、栅极位置和尺寸。掺杂与极化设置GaN HEMT通常不需要刻意掺杂沟道但缓冲层可能需要弱p型或深能级补偿。极化参数必须和材料组分匹配。电极定义与接触设置源漏欧姆接触、栅极肖特基接触接触功函数和界面态需要合理设置。物理模型选择根据仿真目标选择漂移-扩散、能量平衡或者蒙特卡洛模型。求解与后处理先跑平衡态再跑I-V曲线最后提取转移特性、输出特性、跨导、击穿电压等参数。这个流程看起来线性实际操作中经常需要回头调整。比如跑I-V时发现电流太小可能是极化电荷算错了发现跨导异常可能是栅极漏电没控制好。每一次回头都是一次对器件物理的重新理解。2. Atlas中GaN HEMT结构搭建的核心细节2.1 材料参数与极化模型设置GaN HEMT的核心结构是AlGaN/GaN异质结。在Atlas中定义材料时不能只用默认的GaN参数需要手动设置关键材料属性。以下是我常用的一套参数供参考参数GaNAlNAl0.25Ga0.75N禁带宽度 (eV)3.426.24.1电子亲和势 (eV)4.11.93.5相对介电常数9.58.59.2电子迁移率 (cm²/V·s)1200300800饱和速度 (cm/s)2.5e71.5e72.0e7极化模型的设置是重中之重。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 LOC0.0 SPACING0.05 X.MESH LOC0.5 SPACING0.02 X.MESH LOC1.0 SPACING0.05 Y.MESH LOC0.0 SPACING0.01 Y.MESH LOC0.02 SPACING0.001 Y.MESH LOC0.025 SPACING0.0005 Y.MESH LOC0.05 SPACING0.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 NAMEgate WORKFUNCTION4.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 SELBGaN的碰撞电离系数比硅低击穿电场更高。Atlas中的SELB模型需要设置A和B参数对于GaN常用的值大约是A1.0e8 cm⁻¹B1.5e7 V/cm。这些参数不是固定的不同文献给出的值有差异需要根据实验数据标定。3.2 陷阱模型与动态特性模拟GaN HEMT最让人头疼的问题之一是电流崩塌current collapse也就是动态导通电阻增大。这个现象的根源是缓冲层陷阱和表面态。Atlas中可以用TRAPS语句定义陷阱TRAPS REGION1 NUM2 TRAP TYPEACCEPTOR LEVEL0.6 DENSITY1e17 CROSS1e-15 TRAP TYPEDONOR LEVEL0.4 DENSITY5e16 CROSS1e-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 MATERIALGaN TC.CONST1.3 MATERIAL MATERIALGaN HEAT.CAP3.0这里TC.CONST是热导率单位W/(cm·K)GaN大约1.3HEAT.CAP是热容单位J/(cm³·K)GaN大约3.0。自热效应会显著降低饱和电流因为迁移率随温度升高而下降。我做过对比同样的结构不考虑自热时饱和电流约1.2 A/mm考虑自热后降到0.9 A/mm左右降幅约25%。这个差异在功率应用中不能忽略。如果仿真目标只是直流特性自热模块可以不开用等温仿真加快速度。但如果要评估功率器件的热稳定性自热模块必须开。4. 仿真结果分析与常见问题排查4.1 转移特性与阈值电压提取转移特性是GaN HEMT最基本的特性曲线。Atlas中通过扫描栅压固定漏压记录漏极电流SOLVE VGATE0.0 VSTEP-0.1 VFINAL-5.0 NAMEgate这里从0V扫到-5V步长-0.1V。GaN HEMT通常是耗尽型阈值电压为负。提取阈值电压的方法有多种我常用的是线性外推法在转移特性线性区做切线切线与栅压轴的交点即为阈值电压。实测中常见的问题是阈值电压偏移。如果仿真得到的阈值电压比预期偏负很多可能的原因包括极化电荷计算偏大导致2DEG浓度偏高需要更负的栅压才能耗尽。栅极功函数设置不当肖特基势垒偏低。AlGaN势垒层厚度或Al组分设置错误。排查时可以先检查2DEG浓度在平衡态下提取异质结界面处的电子浓度如果超过1.5e13 cm⁻²说明极化电荷偏大需要检查Al组分和应变参数。4.2 输出特性与导通电阻分析输出特性是漏压扫描固定栅压记录漏极电流SOLVE VDRAIN0.0 VSTEP0.5 VFINAL10.0 NAMEdrainGaN HEMT的输出特性通常分为线性区和饱和区。线性区的斜率对应导通电阻饱和区的平坦度反映输出电阻。导通电阻的组成包括源漏接触电阻、沟道电阻、漂移区电阻。仿真中如果导通电阻偏大可能的原因有源漏接触电阻设置过大。沟道迁移率偏低检查迁移率模型参数。栅极长度过长沟道电阻增加。我遇到过一次导通电阻异常偏大的情况排查后发现是源漏电极的RESISTANCE参数设成了1e-3 Ω·cm²而实际应该是1e-5 Ω·cm²量级。改过来后导通电阻恢复正常。4.3 常见收敛问题与排查速查表Atlas仿真GaN HEMT时收敛问题是最常见的障碍。以下是我整理的一份速查表问题现象可能原因解决方法初始平衡态不收敛极化电荷过大网格太粗减小极化电荷加密界面网格栅压扫描时发散栅极漏电模型太激进关闭隧穿模型先跑通再开漏压扫描时电流突跳碰撞电离模型参数不当调整SELB参数减小电压步长瞬态仿真不收敛陷阱时间常数太小增大捕获截面减小时间步长自热仿真温度异常热导率或热容设置错误检查材料热参数确认单位实操心得Atlas的收敛问题很多时候不是模型本身的问题而是网格和步长的问题。我习惯在METHOD语句里设置ITLIMIT50和DELTA0.1给求解器更多迭代次数和更小的初始步长。如果还是不行就用AUTO语句让Atlas自动调整。4.4 仿真与实验偏差的标定思路仿真和实验有偏差是正常的关键是如何系统性地缩小偏差。我的标定思路是先标转移特性调整栅极功函数和AlGaN参数让阈值电压和跨导峰值对上。再标输出特性调整迁移率模型参数和接触电阻让导通电阻和饱和电流对上。最后标动态特性调整陷阱参数让瞬态响应和电流崩塌对上。每一步只调一个参数调完记录避免多个参数同时改导致无法追溯。这个过程可能需要反复迭代但每次迭代都会让你对器件物理的理解更深一层。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语句里设置THREADS4可以显著加快速度。合理设置电压步长步长太大容易发散步长太小浪费时间。我通常先用0.5V步长跑一遍确认没问题后再用0.1V步长跑精细曲线。5.3 后续扩展方向GaN HEMT仿真还有很多可以深入的方向。比如射频特性仿真用小信号AC分析提取截止频率和最大振荡频率。击穿特性仿真用碰撞电离模型和热模型耦合评估击穿电压和雪崩耐量。可靠性仿真用陷阱模型和热模型耦合评估长期工作下的参数退化。电路级仿真把器件仿真结果提取成紧凑模型导入SPICE做电路级仿真。这些方向我还在陆续探索后续有新的心得再分享。仿真这件事说到底是用物理模型去逼近真实器件每一次和实验的对齐都是一次对物理理解的修正。踩过的坑越多对器件的直觉就越准。