📚 课程目录第4章焦耳热与热迁移耦合

4. 焦耳热与热迁移耦合

在微互连结构中，电迁移与热迁移往往不是两个独立过程，而是通过焦耳热紧密耦合在一起。电流通过焊点、铜柱或金属线时，电阻产生焦耳热，使局部温度升高；温度升高又改变电阻率、扩散系数和应力分布，反过来影响电流密度分布与原子输运路径。理解这一耦合链条，是判断失效位置和寿命的关键。

一、焦耳热如何建立温度梯度

焦耳热功率密度可写为 q = J²ρ，其中 J 为电流密度，ρ 为电阻率。由于互连结构中电流密度分布不均匀，焦耳热分布也不均匀。以典型倒装焊点为例，电流从铜柱进入焊料时，在入口附近发生电流集聚，局部电流密度可达平均值的 2～5 倍，焦耳热相应放大 4～25 倍。这种非均匀加热在焊点内部形成温度梯度，典型值在 10～100 K/mm 量级。

温度梯度的方向由热沉位置和电流路径共同决定。若芯片侧为热沉、基板侧散热较差，则焊点顶部温度较低、底部温度较高，温度梯度指向基板方向。这一方向直接决定热迁移的原子输运方向：原子倾向于向低温区迁移，即从热端向冷端扩散。

二、温度梯度与热应力分布

温度梯度不仅驱动原子迁移，还在材料界面产生热应力。由于焊料、铜柱和基板的热膨胀系数不同（SAC305 约 21 ppm/K，Cu 约 17 ppm/K，FR-4 约 15～18 ppm/K），温度变化时各层变形不协调，在界面处产生剪切和拉伸应力。典型焊点中，热应力可达 20～80 MPa，在角部区域甚至超过 100 MPa。

热应力通过两种方式影响质量输运：一是改变原子迁移的驱动力，应力梯度本身可作为原子扩散的驱动力；二是通过应力诱导的空位浓度变化，改变扩散系数。在高温端，压应力使空位浓度降低，原子迁移受阻；在低温端，拉应力使空位浓度升高，原子更容易沉积形成小丘。

三、质量输运的耦合方程

在电迁移与热迁移耦合作用下，原子净通量可表示为：

J_atom = (C·D/kT) · [ Z*·e·ρ·J + Q*·(∇T/T) + Ω·∇σ ]

其中：
  C   — 原子浓度
  D   — 扩散系数，D = D₀·exp(-Ea/kT)
  Z*  — 有效电荷数（电子风驱动力）
  e   — 元电荷
  ρ   — 电阻率
  J   — 电流密度
  Q*  — 热迁移驱动力（传输热）
  ∇T  — 温度梯度
  Ω   — 原子体积
  ∇σ  — 应力梯度

该方程表明，原子通量由三部分叠加：电子风力（第一项）、热迁移力（第二项）和应力梯度力（第三项）。在焊点中，电子风力通常指向阳极，热迁移力指向冷端，两者方向可能相同或相反，取决于温度梯度方向与电流方向的关系。

四、耦合效应的实际表现与判断

以典型 Cu/SAC305/Cu 焊点为例，假设电流密度为 1×10⁴ A/cm²，焊点直径 100 μm，环境温度 125°C，焦耳热使焊点中心温度升高约 15～30 K。若芯片侧为热沉，温度梯度约为 50 K/mm，方向从基板指向芯片。此时电子风力指向芯片侧（电流从芯片流向基板时），热迁移力也指向芯片侧，两者叠加加速原子向芯片侧迁移，导致芯片侧界面出现空洞，基板侧出现小丘。

判断耦合效应是否主导失效，可观察以下特征：

空洞位置：若空洞集中在电流入口附近且靠近冷端，说明电子风力与热迁移力叠加；若空洞出现在热端，则热迁移可能起主导作用。
温度分布：红外热像或有限元仿真显示焊点内部温差超过 10 K 时，热迁移贡献不可忽略。
电流密度依赖性：若失效时间与 1/J² 成正比，说明焦耳热主导；若与 1/J 成正比，说明电子风力主导。
应力分布：通过微区 X 射线衍射或拉曼光谱测量界面应力，若应力梯度与温度梯度方向一致，说明热应力耦合显著。

工程建议：在互连设计中，降低焦耳热耦合效应的优先措施包括：增大焊点截面积以降低电流密度；优化散热路径使温度梯度最小化；选择热膨胀系数匹配的材料组合；在焊点周围增加铜柱或散热通孔，将热量快速导出。

五、耦合效应的量化评估示例

以下给出一个简化的耦合通量估算示例，用于判断电子风力与热迁移力的相对大小。假设焊点材料为 SAC305，工作温度 T = 400 K，电流密度 J = 1×10⁴ A/cm²，电阻率 ρ = 12 μΩ·cm，有效电荷数 Z* = -5，传输热 Q* = -0.1 eV，温度梯度 ∇T = 50 K/mm。

电子风力项：F_em = Z*·e·ρ·J
  = (-5)×(1.6×10⁻¹⁹)×(12×10⁻⁸)×(1×10⁴)
  ≈ -9.6×10⁻²² N

热迁移力项：F_tm = Q*·(∇T/T)
  = (-0.1)×(1.6×10⁻¹⁹)×(50×10³/400)
  ≈ -2.0×10⁻¹⁸ N

比值：F_tm / F_em ≈ 2.1×10³

该示例表明，在给定参数下热迁移力远大于电子风力，说明温度梯度是主导驱动力。但需注意，实际焊点中温度梯度受散热条件影响很大，若散热良好使 ∇T 降至 5 K/mm，则热迁移力下降一个数量级，电子风力可能重新成为主导。因此，耦合效应的判断必须结合具体热边界条件。

注意：上述计算为简化估算，未考虑应力梯度项和温度对扩散系数的非线性影响。实际工程分析应使用有限元耦合仿真，将电-热-力-扩散多物理场同时求解。

综上，焦耳热与热迁移的耦合是微互连失效分析中不可回避的核心问题。电流密度分布决定焦耳热分布，焦耳热建立温度梯度，温度梯度驱动热迁移并产生热应力，热应力又反过来改变原子输运路径。只有将这四个环节作为一个整体来理解，才能准确预测失效位置、判断主导机制并制定有效的可靠性提升方案。