3. 迁移率退化原理

在SOI器件中,自热效应引发的温升会显著加剧晶格振动(声子散射),从而降低载流子迁移率。理解这一退化机制是进行TCAD热仿真配置的理论基础。本节将深入剖析迁移率退化的物理本质、工程近似方法(Mathiessen定则)以及TCAD中常用的温度依赖模型。

3.1 晶格振动加剧与散射机制

载流子在半导体中的输运受到多种散射机制的共同作用。当器件温度因自热效应升高时,晶格原子的热振动幅度增大,导致声子散射(尤其是光学声子散射)的散射概率显著增加。这种散射机制直接破坏了载流子的准动量,从而降低其迁移率。

核心物理图像: 温度升高 → 晶格振动加剧 → 声子密度增加 → 载流子-声子散射增强 → 平均自由程缩短 → 迁移率下降。

在SOI器件中,主要的散射机制包括:

  • 晶格散射(声子散射): 主导机制,随温度升高而增强,迁移率与T-n(n≈1.5~2.5)成正比。
  • 电离杂质散射: 在重掺杂区域显著,随温度升高而减弱(载流子热运动加快,屏蔽效应增强)。
  • 表面粗糙度散射: 在SOI的埋氧层界面和栅氧化层界面处显著,对温度不敏感,但在薄体SOI中影响巨大。
  • 载流子-载流子散射: 在高注入条件下(如饱和区)不可忽略。
散射机制与温度依赖关系 温度 (K) 散射率 (1/s) 晶格散射 (声子) 电离杂质散射 表面粗糙度散射 总散射率 (Mathiessen) 注:总散射率由Mathiessen定则近似,取各散射率之和。温度升高时,晶格散射主导迁移率退化。

3.2 Mathiessen定则:迁移率退化的工程近似

Mathiessen定则提供了一个简洁的工程近似方法,用于计算多种散射机制同时存在时的总迁移率。该定则假设各散射机制相互独立,总散射率等于各散射率之和,从而总迁移率由各迁移率分量的倒数之和的倒数给出。

Mathiessen定则公式:
1 / μ_total = 1 / μ_lattice + 1 / μ_impurity + 1 / μ_surface + 1 / μ_carrier + ...
其中 μ_lattice 为晶格散射限制迁移率,μ_impurity 为电离杂质散射限制迁移率,依此类推。

在TCAD仿真中,Mathiessen定则通常被嵌入到迁移率模型中。对于SOI器件的自热效应分析,最关键的项是晶格散射限制迁移率 μ_lattice,因为它对温度最为敏感。

Mathiessen定则中各迁移率分量的温度依赖特性
散射机制 迁移率分量符号 温度依赖关系 SOI中的重要性
晶格散射(声子) μlattice ∝ T-1.5 ~ T-2.5 极高(自热效应核心)
电离杂质散射 μimpurity ∝ T1.5 中等(重掺杂源漏区)
表面粗糙度散射 μsurface 弱温度依赖 极高(薄体SOI)
载流子-载流子散射 μcarrier ∝ T0.5 ~ T1.0 中等(高注入条件)
Mathiessen定则的局限性: 该定则假设各散射机制完全独立,这在物理上并不严格成立。例如,在强电场下,载流子能量分布偏离平衡态,各散射机制之间存在耦合。TCAD中更精确的模型(如蒙特卡洛方法)可以绕过这一近似,但计算成本较高。在工程仿真中,Mathiessen定则仍是最常用的折中方案。

3.3 温度依赖模型:TCAD中的实现

在TCAD工具(如Sentaurus、Silvaco、COMSOL)中,迁移率的温度依赖通常通过经验公式或物理模型实现。以下是最常用的几种模型:

3.3.1 幂律模型(Power-Law Model)

这是最简单的温度依赖模型,适用于晶格散射主导的区域:

μ(T) = μ(T0) * (T / T0)^(-α)

其中 T0 为参考温度(通常300K),α 为温度指数。对于电子,α ≈ 1.5~2.0;对于空穴,α ≈ 2.0~2.5。该模型在TCAD中常用于快速校准。

3.3.2 Caughey-Thomas模型(C-T模型)

Caughey-Thomas模型是TCAD中最广泛使用的迁移率模型之一,它同时考虑了掺杂浓度和温度的影响:

μ(N, T) = μ_min + (μ_max - μ_min) / (1 + (N / N_ref)^γ)

其中 μ_min、μ_max、N_ref 和 γ 均为温度依赖的参数。温度依赖通过以下方式引入:

μ_max(T) = μ_max(T0) * (T / T0)^(-α)
μ_min(T) = μ_min(T0) * (T / T0)^(-β)
N_ref(T) = N_ref(T0) * (T / T0)^(δ)
γ(T) = γ(T0) * (T / T0)^(-ε)

3.3.3 高场饱和模型(Canali模型)

当横向电场很强时,载流子速度会饱和。Canali模型将低场迁移率与高场饱和速度结合:

v(E) = μ_eff * E / (1 + (μ_eff * E / v_sat)^β)^(1/β)

其中饱和速度 v_sat 也随温度变化:

v_sat(T) = v_sat(T0) * (1 - A * (T - T0))

对于硅,A ≈ 0.87×10-3 K-1。

TCAD迁移率模型选择流程 开始:选择迁移率模型 考虑自热效应? 否 使用常温模型 (幂律模型或C-T模型) 是 选择温度依赖模型 (C-T模型+温度参数) 设置温度依赖参数 α, β, δ, ε, v_sat(T) 等 执行热-电耦合仿真 提示:对于SOI自热效应仿真,必须启用温度依赖的迁移率模型,并正确设置温度指数参数。

3.4 迁移率退化对SOI器件性能的影响

迁移率退化直接导致SOI器件的驱动电流下降和跨导退化。在自热效应严重的薄体SOI器件中,迁移率退化可能使饱和电流降低10%~30%。以下表格总结了不同工作条件下迁移率退化的典型影响:

自热效应导致的迁移率退化对SOI器件性能的影响
工作条件 典型温升 (ΔT) 迁移率退化比例 主要影响
低VDS、低VGS 5~15 K 2%~5% 可忽略
高VDS、中等VGS 30~60 K 10%~20% 驱动电流下降明显
高VDS、高VGS(饱和区) 60~120 K 20%~35% 跨导退化、负微分电阻
脉冲工作(短脉冲) 10~30 K 3%~10% 瞬态热效应,需考虑热时间常数
TCAD仿真实践要点:
  • 在Sentaurus中,使用 Physics(Material="Silicon") { Mobility(PhuMob) } 启用温度依赖的迁移率模型。
  • 对于SOI器件,建议同时启用 Mobility(HighFieldSaturation) 和 Mobility(Enormal) 以考虑高场和垂直电场效应。
  • 校准温度指数 α 时,建议参考实验数据或文献值(电子α≈1.5,空穴α≈2.0)。
  • 在热仿真中,确保 Thermodynamic 模型已启用,以耦合电学和热学方程。

3.5 小结