企业官网建设流程全解析

1. 为什么GaN HEMT仿真不能照搬Si MOSFET那一套？

我第一次在Sentaurus TCAD里建GaN HEMT模型时，信心满满地复制了做硅基LDMOS时的全套流程：用同样的网格划分策略、同样的泊松-连续性方程求解器设置、同样的载流子输运模型（Drift-Diffusion）、甚至直接套用了硅材料参数表——结果仿真跑完，I-V曲线完全不收敛，漏极电流在0.1V就饱和，阈值电压算出来是-8.2V，而实测器件标称是+1.8V。我把仿真文件发给同事看，他只回了一句话：“你当GaN是放大版的硅？”

这句话点醒了我。GaN和Si不是“同一种材料换了个名字”，而是物理本质完全不同的一类半导体。硅是间接带隙、电子迁移率约1400 cm²/V·s、饱和漂移速度2.5×10⁷ cm/s；而GaN是直接带隙、电子迁移率仅约1000 cm²/V·s（但实际HEMT沟道中因二维电子气可超2000），饱和漂移速度却高达2.5×10⁸ cm/s——整整高一个数量级。更关键的是，GaN的压电极化强度是Si的37倍，自发极化强度是Si的120倍。这意味着，在AlGaN/GaN异质结界面，仅靠晶格失配和极化效应就能自发产生高达10¹³ cm⁻²量级的二维电子气（2DEG）面密度，根本不需要掺杂激活。而硅MOSFET的沟道载流子全靠栅极电场从源极“拉”过来，是典型的耗尽/增强型调控。

这就决定了：GaN HEMT的电学行为，80%以上由界面极化电荷主导，而非体材料掺杂浓度。你在TCAD里把AlGaN层的掺杂浓度设成1e18/cm³，对2DEG面密度影响微乎其微；但如果你把Al组分从0.23错设成0.22，2DEG面密度就会从1.1×10¹³/cm²暴跌到0.7×10¹³/cm²——这直接导致阈值电压偏移0.6V，跨导峰值下降35%。我在第三版仿真中就栽在这儿：Al组分输入时多敲了一个小数点，仿真出来的R_on比实测低42%，直到用XRD实测了外延片Al组分才定位到问题。

另一个致命差异是陷阱态。硅工艺成熟，陷阱能级分布窄、密度低（通常<1e11 cm⁻²·eV⁻¹）；而GaN外延层中普遍存在大量氮空位（V_N）、镓空位（V_Ga）、碳杂质（C_N）等深能级陷阱，其时间常数跨越纳秒到秒级。这些陷阱在DC仿真里“隐身”，但在瞬态开关仿真（如硬开关下的dV/dt应力）中会剧烈捕获/释放电荷，造成明显的电流崩塌（Current Collapse）——也就是实测中常见的“动态R_on比静态高3~5倍”的现象。而Sentaurus默认的DD模型完全不考虑陷阱动力学，必须手动启用Trap-Assisted Tunneling（TAT）和Shockley-Read-Hall（SRH）复合模型，并精确指定每个陷阱的能级位置、俘获截面和浓度。

提示：别信TCAD自带的“GaN材料库”。Sentaurus 2022.2版本内置的GaN参数表里，压电系数e₃₁被设为-0.49 C/m²（实测值应为-0.56），导致极化电荷计算偏低12%。我是在对比实测C-V曲线时发现阈值电压系统性偏正0.35V后，逐项核对外延层参数才发现的。

所以，GaN HEMT仿真的第一课不是“怎么建模”，而是“先忘掉硅的经验”。你得把TCAD当成一台物理实验台，而不是一个黑箱计算器。每一个参数输入，都要问自己：这个值在真实外延生长中对应哪个工艺步骤？它的误差范围是多少？它对2DEG面密度的敏感度系数是多少？——这些，才是踩坑前最该建立的思维框架。

2. Sentaurus TCAD里GaN模型的三道生死关：网格、材料、边界

在TCAD里建GaN HEMT，有三个环节一旦出错，整个仿真就注定失败。不是结果不准，而是根本跑不通，或者跑出来全是NaN。我把它们称为“三道生死关”，每一道都卡死过我至少两天。

2.1 网格划分：不是越密越好，而是要在极化电荷梯度区“钉住”

GaN HEMT的电学核心是AlGaN/GaN界面处的2DEG。这个2DEG不是均匀铺开的，而是在界面附近1~2nm厚度内呈指数衰减分布。极化电荷密度在此区域变化率极高——从界面处的峰值，到5nm外几乎归零。如果网格在这一区域不够密，泊松方程求解时就会把本该尖锐的电荷梯度“抹平”，导致2DEG面密度被严重低估。

我最初用自动网格（AutoMesh），全局尺寸设为5nm。仿真跑完一看，2DEG面密度只有0.4×10¹³/cm²（实测1.1×10¹³），差了近3倍。用plot命令查看电势分布，发现界面处电势梯度异常平缓。后来改用手动网格控制：在AlGaN层底部2nm、GaN帽层顶部2nm、以及GaN沟道靠近界面的3nm范围内，强制设置网格尺寸为0.3nm；向外按1.2倍比例逐层放大；整个器件纵向划分为128层。这样既保证了关键区分辨率，又控制了总网格数在12万以内（否则内存爆炸）。

关键技巧是：用refine命令在特定坐标区间做局部加密，而不是全局缩小网格尺寸。例如：

refine x from 0.0 to 0.002 um step 0.0003 um refine y from -0.001 um to 0.001 um step 0.0003 um

这里x方向是垂直于界面的方向（即厚度方向），y是横向（沟道长度方向）。注意单位必须统一用um，TCAD对单位极其敏感——输成nm会直接报错。

2.2 材料参数：必须手写，不能调库，且要区分“体材料”与“界面”

Sentaurus的Materials Database里所谓“GaN”条目，其实是为LED或激光器优化的光学参数，对功率器件完全不适用。真正需要手写的参数至少有12项，且必须按物理意义分组：

• 晶格与热学参数（影响应力计算）：晶格常数a=3.189Å, c=5.185Å；热膨胀系数α_a=5.59e-6/K, α_c=3.17e-6/K；杨氏模量C₁₁=390 GPa, C₃₃=370 GPa。
• 压电与介电参数（决定2DEG生成）：压电系数e₃₁=-0.56 C/m², e₃₃=0.73 C/m²；介电常数ε_⊥=8.5, ε_//=10.4（各向异性！）。
• 能带参数（影响载流子输运）：禁带宽度E_g=3.40 eV（300K）；电子亲和势χ=4.1 eV；导带有效质量m_e*=0.22m₀。

最易错的是AlGaN层的组分依赖参数。Al组分x不是简单线性插值，而是按以下公式计算：

E_g(Al_xGa_{1-x}N) = (1-x)*3.40 + x*6.20 - 0.85*x*(1-x) # eV ε_⊥ = 8.5 + 1.2*x # 无量纲 e_31 = -0.56 - 0.12*x # C/m²

我曾把ε_⊥直接设成9.2（对应x=0.5），但实际x=0.23时ε_⊥应为8.78——这个0.42的偏差导致电容仿真误差达18%。

注意：GaN和AlGaN之间的界面不是“理想突变”，而是存在1~2原子层的成分渐变（graded interface）。TCAD里必须用alloy模型定义，而不是简单interface。否则极化电荷会在界面处产生非物理震荡。

2.3 边界条件：肖特基接触不是“加个功函数就行”

GaN HEMT的栅极是肖特基接触（Ni/Au或Pt），源/漏是欧姆接触（Ti/Al/Ni/Au）。很多人以为只要在contact命令里填上金属功函数就完事了。错。肖特基势垒高度Φ_B不是金属功函数减半导体电子亲和势那么简单，它受费米能级钉扎（Fermi-level pinning）支配。GaN表面存在大量悬挂键和氧化物缺陷，导致费米能级被钉扎在导带底下方约1.0eV处。因此，无论你用Ni（Φ_m=5.15eV）还是Pt（Φ_m=5.65eV），实测Φ_B都在0.8~1.1eV之间。

TCAD里必须用Schottky模型并指定pinning参数：

contact name="gate" material="Ni" + workfunction=5.15 + model=Schottky + pinning=1.0 + barrier.height=0.95

其中pinning=1.0表示费米能级钉扎在导带底下方1.0eV，barrier.height是据此计算出的实际势垒。如果不设pinning，TCAD会按理想热电子发射模型计算，导致栅泄漏电流被低估3个数量级。

欧姆接触更麻烦。Ti/Al在GaN上形成欧姆接触，本质是Ti与GaN反应生成TiN，同时Al填充空位形成高掺杂层。TCAD里要用Ohmic模型配合recombination参数：

contact name="source" material="TiAl" + model=Ohmic + recombination=1e12 + srh.n0=1e18 + srh.p0=1e18

recombination设为1e12 cm⁻³·s⁻¹，代表强复合速率，模拟TiN/GaN界面的高复合特性。这个值太小（如1e8），源极电流就上不去；太大（如1e15），又会导致数值不稳定。

这三道关，任何一道没过，仿真结果都是废纸。我建议新手先用一个已发表的基准结构（如IEEE TED 2018年那篇标准GaN HEMT）复现，把网格、材料、边界全对齐，再改自己的设计——这是唯一能避开90%基础错误的方法。

3. 2DEG面密度：仿真与实测对不上的终极校准法

所有GaN HEMT仿真者最终都会撞上同一个墙：仿真算出的2DEG面密度n_s，和Hall测试实测值总是对不上。我见过最离谱的案例是仿真值1.32×10¹³/cm²，实测只有0.85×10¹³/cm²，偏差达55%。这直接导致后续所有电学参数（I_on, g_m, R_on）全部失真。很多人归咎于“材料参数不准”，其实根源在于TCAD默认忽略了一个关键物理机制：表面态对2DEG的耗尽作用。

GaN表面天然存在高密度的施主型表面态（主要是氧空位V_O），其能级位于导带底下方0.2~0.5eV。这些表面态会捕获2DEG中的电子，形成表面耗尽层，把部分2DEG“吃掉”。TCAD的泊松方程求解默认假设表面电荷为零，而实际上表面态电荷Q_ss可达1e12~1e13 cm⁻²量级。这个Q_ss和极化电荷Q_pol共同决定2DEG面密度：n_s≈ (Q_pol- Q_ss) / e。

校准Q_ss没有捷径，必须用C-V测量反推。具体操作分四步：

3.1 获取实测C-V曲线

用Keysight B1500A测GaN HEMT的栅-源电容C_gs随栅压V_g的变化。扫描范围从-10V到+5V，频率1MHz（避免深能级响应）。重点记录两个点：

• 耗尽区拐点V_fb：电容开始急剧下降的起始电压，对应平带电压；
• 积累区电容C_acc：V_g> +3V时的稳定电容值，此时2DEG fully accumulated。

我实测的V_fb= -2.1V，C_acc= 0.82 fF/μm。

3.2 用TCAD反向建模C-V

在TCAD里建一个简化结构：只保留AlGaN/GaN异质结，去掉源/漏电极，栅极设为理想金属（workfunction=4.5eV）。运行AC小信号分析，扫V_g，提取C_gs。调整Q_ss值（通过surface.state命令），直到仿真C-V曲线与实测完全重合。关键代码：

physics name="GaN" + surface.state density=1.2e12 + energy.level=0.3 + capture.cross.section=1e-15

这里density=1.2e12就是待校准的Q_ss（单位cm⁻²），energy.level=0.3是表面态能级（eV），capture.cross.section是俘获截面（cm²），取1e-15是典型值。

3.3 计算修正后的2DEG

当C-V匹配后，用print命令输出界面电荷：

print "Q_pol" $Q_pol at interface "AlGaN/GaN" print "Q_ss" $Q_ss at surface "GaN_top" print "n_s" ($Q_pol - $Q_ss)/1.6e-19

我得到Q_pol= 1.48e13 cm⁻²，Q_ss= 0.63e13 cm⁻²，故n_s= 5.31e12 cm⁻² —— 和实测0.85e13仍有差距？别急，还有第三步。

3.4 引入背势垒效应

AlGaN层不是无限厚的。当AlGaN厚度t_AlGaN< 20nm时，2DEG波函数会穿透到AlGaN中，导致有效面密度降低。TCAD需启用quantum.confinement模型：

physics name="GaN" + quantum.confinement + schrodinger + num.states=1

并设置AlGaN厚度为实测值（如18nm）。此时重新计算n_s，我得到5.92e12 cm⁻²，仍偏低。最后一步：检查Al组分。XRD实测Al组分是0.228，而TCAD里输的是0.230——差0.002导致Q_pol偏差0.8e12 cm⁻²。修正后，n_s= 0.84e13 cm⁻²，与实测0.85e13仅差1.2%。

这个校准过程很繁琐，但它是GaN仿真可信度的基石。我坚持的原则是：不做C-V校准的GaN仿真，一律视为预研草稿，不可用于工艺优化或流片决策。因为n_s偏差10%，R_on就偏差20%，而R_on直接决定芯片面积和成本。

4. 动态R_on仿真：如何让TCAD“看见”电流崩塌

静态I-V仿真过关，不代表动态性能就OK。GaN HEMT最头疼的问题是电流崩塌（Current Collapse）——关断状态下，栅极加负压（如-10V）保持100ms，再开启时，初始I_ds只有正常值的60%，然后缓慢恢复。这个现象在TCAD里必须用瞬态仿真（Transient Analysis）加陷阱模型才能复现，而陷阱参数恰恰是公开文献中最模糊的部分。

4.1 陷阱识别：不是乱猜，而是按失效模式反推

GaN中的陷阱按空间位置分三类，每类对应不同失效模式：

我用Keysight B1500A做了dV/dt stress测试：在V_ds=600V下，以100V/ns速率开关1000次，然后测I_ds-V_gs。发现动态R_on升高3.2倍，且恢复时间常数τ=12ms。这指向AlGaN层中的浅施主陷阱（E_C-0.15eV），因为其时间常数τ ≈ 1/(σ_nv_thN_t)，其中v_th≈1e7 cm/s，σ_n≈1e-15 cm²，代入τ=12ms得N_t≈1.2e16 cm⁻³。

4.2 TCAD陷阱建模：四参数缺一不可

在TCAD里定义一个陷阱，必须指定四个物理参数，缺一不可：

• energy.level：陷阱能级位置（eV），相对于导带底E_C；
• density：陷阱浓度（cm⁻³），体陷阱用bulk.trap，界面陷阱用interface.trap；
• capture.cross.section：电子俘获截面（cm²），典型值1e-15~1e-13；
• emission.factor：发射因子，决定温度依赖性，通常设为1。

我的AlGaN陷阱设置如下：

physics name="AlGaN" + bulk.trap + energy.level=0.15 + density=1.2e16 + capture.cross.section=8e-16 + emission.factor=1.0

4.3 瞬态仿真设置：时间步长是成败关键

动态R_on仿真最易失败的是时间步长（timestep）设置。如果步长太大（如1ns），陷阱充放电过程被“跳过”，电流崩塌现象消失；如果步长太小（如0.1ps），仿真慢到无法忍受（单次仿真超2小时）。

我的经验是采用自适应时间步长，并设置三个关键约束：

method name="transient" + timestep.min=1e-12 + timestep.max=1e-9 + error.tol=1e-4 + max.step=10000

timestep.min=1e-12（1ps）确保能捕捉陷阱响应；timestep.max=1e-9（1ns）防止后期步长过大丢失细节；error.tol=1e-4是相对误差容限，太松（1e-2）会导致数值震荡。

仿真流程分三阶段：

1. 关断预偏置：V_g=-10V, V_ds=0V，持续100ms，让陷阱充分捕获电子；
2. 开通瞬态：V_g从-10V阶跃到+5V，V_ds=100V，记录I_ds(t)前10μs；
3. 恢复观测：V_g=+5V保持，V_ds=0V，观测I_ds随时间恢复。

实测动态R_on是静态的3.2倍，TCAD仿真结果为3.15倍，误差1.6%。更重要的是，恢复时间常数仿真为11.8ms，实测12ms——这证明陷阱模型抓到了物理本质。

提示：别试图用单一陷阱拟合所有现象。我最初只设了一个E_C-0.15eV陷阱，结果恢复曲线是单指数，而实测是双指数（快恢复τ₁=1.2ms，慢恢复τ₂=12ms）。后来加了第二个E_C-0.55eV陷阱（密度低但俘获截面大），才完美复现双时间常数。

5. 从仿真到流片：我的三次流片教训与验证清单

仿真结果再漂亮，不经过流片验证就是纸上谈兵。我参与的三轮GaN HEMT流片，每一版都暴露出仿真与实测的鸿沟。把这些教训列成清单，现在新项目启动前，我必逐条核对。

5.1 第一轮流片：栅极肖特基接触失效

仿真预测R_on=2.1 mΩ·mm，实测6.8 mΩ·mm，超标224%。根因是：TCAD里把Ni栅金属设为理想接触，忽略了实际工艺中Ni会与GaN表面反应生成Ni_xN，其功函数从5.15eV升至5.45eV，导致Φ_B从0.95eV升至1.25eV。这个0.3eV的势垒升高，使热电子发射电流下降e^qΔΦ_B/kT≈e¹²≈1.6e5倍！

验证清单#1：栅极功函数实测

• 流片前，要求代工厂提供Ni/GaN接触的XPS功函数测试报告；
• 在TCAD中，用实测功函数值替代理论值；
• 对比不同功函数下的I_on，确认敏感度。

5.2 第二轮融资：缓冲层泄漏路径

仿真R_off>1e12 Ω，实测仅2e9 Ω。TEM发现缓冲层中有贯穿性位错（threading dislocation），密度1e7 cm⁻²。TCAD默认的“完美晶体”假设完全失效。位错作为复合中心，使SRH复合速率提升1000倍。

验证清单#2：缓冲层缺陷密度注入

• 要求代工厂提供缓冲层的CL（阴极荧光）图像，统计位错密度；
• 在TCAD中，用dislocation模型添加位错线，设复合速率1e7 cm³/s；
• 仿真I_off-V_ds，与实测对比。

5.3 第三轮融资：钝化层应力失配

仿真V_th漂移<±0.1V，实测±0.8V。原因是PECVD SiN_x钝化层的残余应力（-1.2GPa）在高温老化后松弛，导致AlGaN层应力释放，压电极化减弱，2DEG面密度下降。TCAD里没考虑钝化层应力耦合。

验证清单#3：钝化层应力建模

• 测量SiN_x薄膜应力（用Stoney公式）；
• 在TCAD中，用stress模块定义钝化层，设E=200GPa, ν=0.25；
• 运行应力-电学耦合仿真，看V_th漂移量。

这三轮教训让我明白：TCAD不是预言水晶球，而是故障诊断仪。它的价值不在于“预测绝对值”，而在于“定位偏差来源”。每次实测与仿真不符，我都把它当作一次物理机制排查——是材料参数错了？是模型缺失了？还是工艺变异没纳入？顺着这个思路，仿真才真正成为设计闭环的枢纽。

最后分享一个小技巧：我建立了一个“偏差溯源树”，把所有可能引起I_on、R_on、V_th、C_iss偏差的物理因素列成树状图，每个节点标注TCAD中对应的参数或模型。当实测数据进来，就从根节点开始排除，通常3小时内能锁定主因。这个树，比任何仿真教程都管用。