3. Black方程参数提取与电迁移寿命外推
Black方程与寿命外推的物理含义
电迁移失效的加速寿命通常用 Black 方程描述：
其中 MTTF 为中位失效时间（h），J 为电流密度（A/cm²），T 为绝对温度（K），k = 8.617×10⁻⁵ eV/K，n 为电流密度指数，Ea 为激活能（eV），A 为与材料、结构、工艺相关的常数。本案例的任务是：用一组加速老化试验数据，同时拟合出 n 与 Ea，再把结果外推到实际工作条件（Juse、Tuse）下的寿命。
关键点在于：n 与 Ea 必须由两组独立变量分别驱动。若只改变温度而电流密度恒定，只能得到 Ea；若只改变电流密度而温度恒定，只能得到 n。因此试验矩阵必须同时覆盖两个维度。
加速试验矩阵与示例数据
下表为教学示例数据（假设值，用于演示拟合流程，非实测结果）。材料体系为 Al-Cu 互连 / SiO₂ 介质，测试结构为等宽等长的金属条，失效判据为电阻上升 20%。
| 组号 | J (MA/cm²) | T (°C) | T (K) | MTTF (h) |
|---|---|---|---|---|
| 1 | 1.0 | 200 | 473.15 | 1200 |
| 2 | 2.0 | 200 | 473.15 | 310 |
| 3 | 4.0 | 200 | 473.15 | 78 |
| 4 | 2.0 | 250 | 523.15 | 62 |
| 5 | 2.0 | 300 | 573.15 | 14 |
组 1–3 固定温度、变电流密度，用于提取 n；组 2、4、5 固定电流密度、变温度，用于提取 Ea。组 2 是两个维度的交点，起锚定作用。
参数提取：两步对数线性回归
对 Black 方程两边取自然对数：
第一步：固定 T = 473.15 K，用组 1–3 拟合 n。
J = [1.0, 2.0, 4.0] # MA/cm^2
MTTF = [1200, 310, 78] # h
x = ln(J) = [0.0000, 0.6931, 1.3863]
y = ln(MTTF)= [7.0901, 5.7366, 4.3567]
斜率 = (y3 - y1) / (x3 - x1)
= (4.3567 - 7.0901) / (1.3863 - 0.0000)
= -2.7334 / 1.3863
= -1.9717
n = 1.97 （电流密度指数，接近经典值 2）用组 1–3 的中间点校验：x₂ = 0.6931，预测 y₂ = 7.0901 − 1.9717×0.6931 = 5.7235，实测 5.7366，偏差约 0.013（≈1.3%），说明该温度下幂律关系成立。
第二步：固定 J = 2.0 MA/cm²，用组 2、4、5 拟合 Ea。
T = [473.15, 523.15, 573.15] # K
MTTF = [310, 62, 14] # h
x = 1/T = [0.0021135, 0.0019115, 0.0017447]
y = ln(MTTF) = [5.7366, 4.1271, 2.6391]
斜率 = (y3 - y1) / (x3 - x1)
= (2.6391 - 5.7366) / (0.0017447 - 0.0021135)
= -3.0975 / (-0.0003688)
= 8398.6
Ea = 斜率 × k = 8398.6 × 8.617e-5 = 0.7237 eV用组 4 校验：x₂ = 0.0019115，预测 y₂ = 5.7366 + 8398.6×(0.0019115 − 0.0021135) = 5.7366 − 1.6965 = 4.0401，实测 4.1271，偏差约 0.087（≈2.1%），Arrhenius 关系基本成立。
常数 A 的求解与外推
把 n = 1.97、Ea = 0.7237 eV 代入任一组数据求 A。以组 2 为例：
ln(A) = ln(MTTF) + n·ln(J) - Ea/(k·T)
= 5.7366 + 1.97×0.6931 - 0.7237/(8.617e-5 × 473.15)
= 5.7366 + 1.3654 - 17.752
= -10.650
A = exp(-10.650) = 2.37e-5 （h·(A/cm^2)^n 量纲）外推到实际工作条件，例如 Juse = 0.5 MA/cm²、Tuse = 100 °C = 373.15 K：
ln(MTTF_use) = ln(A) - n·ln(J_use) + Ea/(k·T_use)
= -10.650 - 1.97×ln(0.5) + 0.7237/(8.617e-5 × 373.15)
= -10.650 + 1.3654 + 22.508
= 13.223
MTTF_use = exp(13.223) ≈ 5.5e5 h ≈ 63 年该结果属于由示例数据直接推出的外推值，用于演示方法，不代表真实器件的实测寿命。
拟合结果与判读逻辑
结果判读与工程判断
- n ≈ 2：与经典电迁移理论（n = 2）一致，说明该温度区间内失效由电子风力主导，幂律模型适用。若拟合得到 n 明显偏离 2（如 n < 1.5 或 n > 3），需检查是否存在焦耳热自加热、电流拥挤或测试结构边缘效应。
- Ea ≈ 0.72 eV：处于 Al 基互连晶格扩散与晶界扩散之间的典型区间。若 Ea 偏低（< 0.5 eV），提示表面/界面扩散通道占主导，可能对应介质界面质量差或钝化层不完整。
- 外推跨度：从 200 °C 外推到 100 °C，温度应力比约 1.27，寿命放大倍数约 1.8×10³。外推跨度越大，n 与 Ea 的微小误差被指数放大，因此建议外推温度不低于实际工作温度 + 20 °C，且至少保留一个中间温度点做校验。
- 置信区间：上述为点估计。工程上应给出 n 与 Ea 的 95% 置信区间，并用最坏边界值计算保守寿命。若每组样本数 ≥ 15，可用对数正态分布拟合中位寿命并给出 σ。