2. SiC材料特性与物理模型选择

在Sentaurus中做SiC功率器件仿真,第一步不是画结构,而是把材料参数和物理模型选对。4H-SiC与6H-SiC虽然都是SiC,但晶型堆垛顺序不同,导致带隙、迁移率各向异性、碰撞电离系数等关键参数存在明显差异。工程上绝大多数功率器件(MOSFET、JBS、SBD)使用4H-SiC,6H-SiC更多出现在早期研究或特定高压场景。本节围绕晶型结构、各向异性迁移率、不完全电离、带隙变窄、碰撞电离系数、热导率与温度模型,说明它们在Sentaurus材料文件中的配置方法。

一、晶型结构与各向异性坐标约定

4H-SiC为六方晶系,晶格常数a≈3.073 Å、c≈10.053 Å,堆垛顺序为ABAC,带隙约3.26 eV;6H-SiC堆垛顺序为ABCACB,带隙约3.02 eV。两者都是间接带隙半导体,但4H-SiC电子迁移率更高、各向异性更弱,因此成为主流。在Sentaurus中,晶型本身不直接作为关键字,而是通过各向异性参数体现。器件坐标约定为:c轴沿y方向(或z方向,取决于结构定义),载流子沿c轴与垂直c轴方向迁移率不同。4H-SiC电子迁移率垂直c轴约1000 cm²/(V·s)、沿c轴约800 cm²/(V·s);空穴迁移率垂直c轴约120 cm²/(V·s)、沿c轴约60 cm²/(V·s)。这些数值是示例量级,实际应以工艺与文献标定为准。

二、各向异性迁移率模型

Sentaurus中迁移率各向异性通过AnisotropicMobility或材料参数中的方向分量实现。常用做法是在材料文件里定义muMin、muMax、Nref、alpha等参数,并分别给出平行与垂直c轴的分量。例如:

Material = "4H-SiC" {
  # 各向异性迁移率(示例参数,需按实际标定)
  AnisotropicMobility {
    eMuMax = (1000, 800)   # (垂直c轴, 平行c轴) cm^2/Vs
    hMuMax = (120, 60)
    eMuMin = (40, 30)
    hMuMin = (15, 10)
    Nref   = 1e17
    alpha  = 0.7
  }
}

若器件沟道沿c轴,则迁移率取平行分量;若沟道垂直c轴,则取垂直分量。Sentaurus在求解时会根据电流方向自动插值。需要注意:各向异性迁移率必须与eMobility模型配合,不能只写各向异性而忽略温度与掺杂依赖。

三、不完全电离模型

SiC中Al、B、N等掺杂剂离化能较大(Al受主约0.2 eV,N施主约0.05 eV),室温下并非全部电离。Sentaurus提供IncompleteIonization模型,需在物理模型段启用:

Physics {
  IncompleteIonization
  eDensity = IncompleteIonization
  hDensity = IncompleteIonization
}

材料文件中需给出离化能E_d、E_a和简并因子。例如4H-SiC中N施主E_d=0.05 eV,Al受主E_a=0.2 eV。若忽略该模型,仿真在低温和高掺杂区会高估载流子浓度,导致阈值电压和导通电阻偏差。建议在沟道掺杂>1e17 cm⁻³或温度<300 K时务必启用。

四、带隙变窄与碰撞电离系数

重掺杂下带隙变窄(BGN)会改变本征载流子浓度。Sentaurus中通过BandgapNarrowing模型启用,材料文件需提供Bgn参数。4H-SiC的BGN系数通常取Bgn = 1.5e-2 eV·cm量级,具体随掺杂浓度变化。碰撞电离系数决定击穿电压,Sentaurus使用ImpactIonization模型,常用Okuto或vanOverstraeten形式。4H-SiC电子与空穴电离系数差异大,需分别给出a、b、Ecrit等参数。例如:

ImpactIonization {
  Model = Okuto
  eA = 1.0e7
  eB = 1.5e7
  hA = 1.0e7
  hB = 3.0e7
}

这些参数为示例,实际应参考器件击穿电压标定。若击穿电压仿真偏高,可适当降低eB或hB。

五、热导率与温度模型

SiC热导率约3.7 W/(cm·K)(4H-SiC,室温),随温度升高而下降,近似k(T)=k300*(T/300)^(-1.2)。Sentaurus中通过ThermalConductivity模型配置,材料文件写:

ThermalConductivity {
  k300 = 3.7
  alpha = -1.2
}

温度模型还需启用Temperature依赖的迁移率、饱和速度等。Sentaurus中常用Temperature关键字激活,并在材料文件中给出beta、Tnom等。若做电热耦合仿真,还需在Physics中加Thermodynamic或Hydrodynamic。

4H-SiC / 6H-SiC 材料参数与Sentaurus模型映射 晶型结构 4H-SiC: ABAC Eg ≈ 3.26 eV a=3.073 Å, c=10.053 Å 6H-SiC: ABCACB Eg ≈ 3.02 eV 主流功率器件选4H Sentaurus 物理模型 AnisotropicMobility IncompleteIonization BandgapNarrowing ImpactIonization ThermalConductivity Temperature 材料参数文件 Material = "4H-SiC" eMuMax = (1000,800) hMuMax = (120,60) E_d = 0.05 eV E_a = 0.2 eV k300 = 3.7 W/cmK alpha = -1.2 参数影响链与仿真结果判读 各向异性迁移率 沟道方向选择 不完全电离 载流子浓度修正 带隙变窄 本征浓度变化 碰撞电离 击穿电压 热导率 温度分布 结果判读:阈值电压偏移 → 检查不完全电离;击穿电压偏差 → 检查碰撞电离系数 导通电阻异常 → 检查各向异性迁移率方向;温升过高 → 检查热导率温度系数

六、配置流程与检查要点

实际配置时,建议按以下顺序操作:先在Material段定义4H-SiC或6H-SiC,写入各向异性迁移率、离化能、BGN系数、碰撞电离参数、热导率;再在Physics段启用IncompleteIonization、BandgapNarrowing、ImpactIonization、ThermalConductivity和Temperature;最后在Solve段设置温度扫描或电热耦合。检查时重点看:迁移率方向是否与沟道一致、离化能是否与掺杂剂匹配、碰撞电离系数是否导致击穿电压合理、热导率温度系数是否使温升收敛。若仿真不收敛,可先关闭不完全电离和BGN,逐步加入模型定位问题。

建议在项目目录下建立4H-SiC.par和6H-SiC.par两个材料文件,通过Material关键字切换,避免每次修改主命令文件。参数标定可参考Sentaurus自带MaterialDB中的SiC条目,再根据实际器件测试数据微调。