4. 宽禁带材料物理模型选择

在 Sentaurus 中仿真 SiC、GaN 等宽禁带器件时,物理模型的选择直接决定结果是否可信。宽禁带材料具有禁带宽度大(SiC 约 3.26 eV,GaN 约 3.4 eV)、临界击穿场强高(SiC 约 2~3 MV/cm)、热导率差异大、掺杂激活率低等特征,这些特征使许多在硅中“默认可用”的模型在宽禁带中必须重新审视。本节围绕四类核心模型——漂移-扩散(DD)、热力学(Thermodynamic)、迁移率模型与复合模型——说明其适用边界、关键参数以及校准方法。

一、漂移-扩散与热力学模型的取舍

漂移-扩散模型假设载流子与晶格处于热平衡,仅用准费米能级梯度驱动输运,计算量小、收敛性好,适合沟道输运、低场区、常规 I-V 特性提取。但在宽禁带功率器件中,高场区(如 SiC MOSFET 的 JFET 区、GaN HEMI 的栅漏漂移区)载流子温度显著高于晶格温度,DD 模型会低估碰撞电离与速度饱和效应,导致击穿电压偏高、导通电阻偏低。

热力学模型在 DD 基础上引入晶格热流方程,求解载流子与晶格的能量平衡,能反映自热效应与温度分布。对于 SiC 器件,由于功率密度高、热导率虽高但芯片面积小,自热不可忽略;对于 GaN HEMT,二维电子气(2DEG)沟道局部温升可达数百开尔文,必须启用热力学模型并耦合热边界条件。选择原则是:若器件工作电流密度超过 10 A/cm² 或存在明显高场区,优先使用热力学模型;若仅做低场参数提取或快速趋势对比,DD 模型足够。

宽禁带器件物理模型选择与耦合关系 漂移-扩散模型 (DD) 准费米能级驱动 · 热平衡假设 适用:低场区、快速趋势对比 热力学模型 (Thermodynamic) 耦合晶格热流方程 · 自热效应 适用:高功率密度、高场区 迁移率模型 Caughey-Thomas · 高场饱和 掺杂/温度/场强依赖 复合模型 SRH · Auger · 辐射复合 决定少子寿命与漏电流 模型耦合链 电场分布 E(x) 载流子温度 T_c(x) 迁移率 μ(E,T) 复合率 R(n,p,T) 关键校准参数 Caughey-Thomas vsat, β, μ0, α SRH 寿命 τn, τp, Et, Nt Auger 系数 Cn, Cp (10⁻³¹ cm⁶/s) 热边界条件 Rth, T_amb, 热沉 碰撞电离 α(E), 临界场 E_c

二、迁移率模型的适用性与参数校准

宽禁带材料的迁移率模型不能直接沿用硅参数。SiC 的体迁移率较低(4H-SiC 电子约 800~1000 cm²/V·s,空穴约 100~120 cm²/V·s),且各向异性明显;GaN 的 2DEG 迁移率可达 1500~2000 cm²/V·s,但受合金无序散射和界面粗糙度影响大。Sentaurus 中常用的 Caughey-Thomas 模型形式为:

μ(E) = μ0 / [1 + (μ0·E / vsat)^β]^(1/β)

其中 μ0 为低场迁移率,vsat 为饱和速度,β 为拟合指数。对 4H-SiC 电子,典型值 μ0 ≈ 900 cm²/V·s,vsat ≈ 2.0×10⁷ cm/s,β ≈ 1.2;对 GaN 2DEG,μ0 ≈ 1800 cm²/V·s,vsat ≈ 2.5×10⁷ cm/s,β ≈ 1.0。校准方法:先用霍尔测试或低场 I-V 提取 μ0,再用高场饱和区(如 HEMT 输出特性饱和段)拟合 vsat 与 β。注意温度依赖:μ0(T) = μ0(300K)·(T/300)^(-γ),SiC 的 γ 约 1.5~2.0,GaN 约 1.0~1.5。

三、复合模型的适用性与参数校准

宽禁带材料的少子寿命通常较短(SiC 约 0.1~1 μs,GaN 约 1~10 ns),SRH 复合主导低注入区,Auger 复合在高注入区(如大电流导通)不可忽略。Sentaurus 中 SRH 模型需要指定 τn、τp、陷阱能级 Et 和浓度 Nt;Auger 模型需要 Cn、Cp。对 4H-SiC,典型 Cn ≈ 5×10⁻³¹ cm⁶/s,Cp ≈ 2×10⁻³¹ cm⁶/s;GaN 的 Auger 系数更小,约 10⁻³²~10⁻³¹ cm⁶/s 量级。辐射复合在 GaN 基 LED 中必须启用,系数 B 约 10⁻¹⁰~10⁻⁹ cm³/s。

校准策略:用二极管反向恢复或开路电压衰减(OCVD)提取 SRH 寿命;用大注入条件下的 I-V 曲线拟合 Auger 系数;用光致发光(PL)或电致发光(EL)效率曲线提取辐射复合系数。注意:宽禁带材料中 SRH 寿命对温度敏感,τ(T) = τ(300K)·(T/300)^(-1.5) 可作为一阶近似。

参数校准顺序建议:先校准迁移率(低场 μ0 → 高场 vsat、β),再校准 SRH 寿命,最后校准 Auger 与辐射复合系数。每步用独立实验数据验证,避免多参数同时拟合导致不唯一。

四、模型组合与仿真设置示例

以 4H-SiC MOSFET 导通特性仿真为例,推荐模型组合:

物理模型Sentaurus 关键字关键参数示例
漂移-扩散eQuantumPotential 关闭默认
热力学Thermodynamic热沉 300 K,Rth = 1 K/W
迁移率Philips Unified / Caughey-Thomasμ0=900, vsat=2e7, β=1.2
SRHSRHτn=τp=0.5 μs, Et=0.5 eV
AugerAugerCn=5e-31, Cp=2e-31 cm⁶/s
碰撞电离Okuto / van Overstraetenα 参数按 SiC 拟合

仿真后应检查:导通电阻随温度的变化趋势(正温度系数?)、饱和电流是否受 vsat 限制、关断漏电流是否由 SRH 主导。若导通电阻随温度下降,说明迁移率温度模型参数 γ 偏小;若漏电流远大于预期,检查 SRH 寿命是否设置过短或陷阱能级是否接近禁带中央。

常见错误:直接使用 Sentaurus 默认硅参数仿真 SiC/GaN,会导致击穿电压偏低、导通电阻偏小、漏电流偏大。务必逐项替换为宽禁带材料参数,并用实验数据校准。

综上,宽禁带器件仿真中模型选择的核心逻辑是:低场用 DD,高场/自热用热力学;迁移率必须用 Caughey-Thomas 或 Philips 并校准 vsat 与 β;复合模型至少启用 SRH,大注入时加 Auger,光电器件加辐射复合。参数校准应分步、独立验证,避免过拟合。