📚 课程目录第2章热输运物理基础

2. 热输运物理基础

宽禁带半导体（GaN、SiC、β-Ga₂O₃、AlN、金刚石等）的热输运与硅器件有本质差异：它们往往同时具备高德拜温度、强极性键合、宽禁带、可高掺杂或可形成二维电子气等特征，导致热量的载体不再只是单一机制。要理解器件自热、结温分布与热失效，必须先厘清“谁在导热、沿哪个方向导、在什么尺度上还服从傅里叶定律”。

一、三种热载流子：声子、电子、光子

在固体中，热能的输运由三类准粒子承担。第一类是声子，即晶格振动的量子化，是绝缘体和半导体中占绝对主导的热载体。宽禁带材料键合强、原子质量轻（如 AlN、SiC、金刚石），声子群速度高、德拜温度高，因此本征热导率可以很高；但一旦引入同位素无序、点缺陷、位错、界面或合金散射，声子平均自由程迅速下降，热导率随之塌陷。

第二类是电子。在重掺杂 n 型 GaN 或 SiC 中，电子既导电也导热，遵循维德曼–弗兰兹关系 κ_e = L·σ·T（L 为洛伦兹数，约 2.44×10⁻⁸ W·Ω/K²）。但宽禁带材料掺杂效率低、载流子浓度通常低于 10¹⁹ cm⁻³，电子热导率一般只占总热导率的百分之几，只有在极高掺杂或金属电极、互连层中才不可忽略。

第三类是光子热载流子，即辐射传热。在高温（>800 K）或大功率光电器件（LED、激光器）中，电子–空穴复合产生的光子、以及热辐射本身会携带能量。对 GaN 基 LED，光子逃逸与再吸收、以及“热光子”在芯片内部的多次反射，会改变局部热流路径。在常规功率器件工作温度下，光子贡献通常远小于声子，但在光泵浦、激光退火、瞬态热反射实验中必须单独建模。

二、热导率张量：各向异性的数学表达

热导率在一般晶体中不是标量，而是二阶张量 κ_ij。傅里叶定律的矢量形式为 q_i = −κ_ij ∂T/∂x_j，即热流方向不一定与温度梯度方向平行。对宽禁带材料，这一点尤其重要：

4H-SiC：六方晶系，κ_a（垂直 c 轴）≈ 490 W/(m·K)，κ_c（沿 c 轴）≈ 390 W/(m·K)，各向异性约 25%。
GaN：纤锌矿结构，κ_a ≈ 200 W/(m·K)，κ_c ≈ 180 W/(m·K)，差异较小但外延层与衬底界面处仍会造成热流偏折。
β-Ga₂O₃：单斜晶系，κ 在 [010]、[100]、[001] 方向分别约为 27、13、11 W/(m·K)，各向异性超过 2 倍，是热管理设计的核心难点。

工程上处理张量有两种做法：一是坐标变换，把晶体主轴坐标系下的对角张量旋转到器件坐标系；二是直接测量或仿真得到 κ_xx、κ_yy、κ_zz 及交叉项。对 β-Ga₂O₃ 功率器件，若忽略各向异性，结温预测误差可达 30% 以上。

三、傅里叶框架：适用条件与边界

傅里叶定律 q = −κ∇T 是宏观经验定律，其成立前提是：特征长度 L 远大于热载流子平均自由程 Λ，且时间尺度远大于声子弛豫时间 τ。对 SiC、GaN 在室温下，Λ 约 10–100 nm，因此当器件特征尺寸大于 1 μm、热源调制频率低于 1 GHz 时，傅里叶框架足够准确。稳态热仿真、常规热阻提取、封装级热分析都建立在此框架上。

但宽禁带器件正不断逼近极限：GaN HEMI 的二维电子气沟道厚度仅几纳米，SiC MOSFET 的沟道长度可小于 0.5 μm，β-Ga₂O₃ 的纳米膜厚度可到 100 nm 以下。此时 L 与 Λ 可比，傅里叶定律失效，必须引入非傅里叶框架。

四、非傅里叶导热：从 BTE 到热波

非傅里叶导热主要有三条技术路线：

玻尔兹曼输运方程（BTE）：直接求解声子分布函数 f(r, k, t) 的演化，包含漂移、扩散与散射项。对稳态问题，常用弛豫时间近似（RTA）或离散坐标法（DOM）数值求解。BTE 能自然给出弹道–扩散转变、界面热阻与尺寸效应，是纳米尺度热输运的“金标准”。
双曲热传导模型（Cattaneo–Vernotte）：在傅里叶定律中引入弛豫时间 τ：τ ∂q/∂t + q = −κ∇T。当热源调制频率 ω 满足 ωτ ~ 1 时，热流不再瞬时响应温度梯度，出现“热波”传播，温度场以有限速度传播而非无限扩散。
声子流体动力学：当声子–声子正常散射远强于倒逆散射时，声子集体表现为黏性流体，出现 Poiseuille 型热流与第二声。在石墨、金刚石等超高 κ 材料中已被实验观测到。

判断流程（工程实用）：先估算 Λ（可由 κ、C_v、v_g 估算：Λ ≈ 3κ/(C_v·v_g)），再与器件最小特征尺寸 L_min 比较。若 L_min/Λ > 10，用傅里叶；1 < L_min/Λ < 10，用 BTE 或双曲模型修正；L_min/Λ < 1，必须用弹道输运或蒙特卡洛方法。

五、教学示例：GaN HEMT 沟道热输运估算

假设一个 GaN HEMT 的二维电子气沟道，工作温度 T = 400 K，GaN 的 κ_c ≈ 160 W/(m·K)，比热 C_v ≈ 2.5×10⁶ J/(m³·K)，声子群速度 v_g ≈ 6000 m/s。估算声子平均自由程：

Λ ≈ 3κ / (C_v · v_g)
  = 3 × 160 / (2.5e6 × 6000)
  ≈ 3.2 × 10⁻⁸ m = 32 nm

若沟道热源长度 L_h ≈ 0.25 μm = 250 nm，则 L_h/Λ ≈ 7.8，处于过渡区。此时用纯傅里叶仿真会低估峰值温度约 10–15%，因为弹道效应使热源附近热阻增大。修正方法：在傅里叶模型中引入“有效热导率” κ_eff = κ / (1 + Λ/L_h)，得到 κ_eff ≈ 160 / (1 + 32/250) ≈ 142 W/(m·K)，再代入仿真，结温预测更接近 BTE 结果。

注意：上述数值为教学示例中的合理假设，并非针对某一具体器件的实测结果。实际器件的 κ、Λ 会随温度、掺杂浓度、缺陷密度和应力状态变化，必须通过时域热反射（TDTR）、3ω 法或拉曼热测量等手段标定。

六、从物理到工程：三条实用结论

先判机制，再选模型。不要默认所有宽禁带器件都用傅里叶。β-Ga₂O₃ 和纳米膜 GaN 必须考虑各向异性与弹道效应。
张量方向要与热流方向对齐。在版图设计时，让高 κ 晶向对准主热流路径，可显著降低峰值结温。
界面热阻常是瓶颈。GaN/SiC、GaN/Si、Ga₂O₃/金刚石等异质界面的热边界电阻（TBR）可达 10⁻⁸–10⁻⁷ m²·K/W，等效于数百纳米厚的低 κ 层，必须在热阻网络中显式建模。

掌握声子、电子、光子三类载流子的相对贡献，理解热导率张量的方向性，并能根据特征尺寸与平均自由程之比选择傅里叶或非傅里叶框架，是后续分析宽禁带器件自热、热失效与热管理方案的基础。