4、MOSFET老化模型：阈值电压漂移、迁移率退化、亚阈值摆幅变化，以及老化模型参数提取
MOSFET的老化（Aging）是集成电路长期可靠性面临的核心挑战。在先进工艺节点下，偏置温度不稳定性（BTI，包括NBTI和PBTI）以及热载流子注入（HCI）效应会导致器件参数随时间发生漂移，进而影响电路时序、功耗和噪声容限。本章将深入剖析MOSFET老化模型的三大核心参数——阈值电压漂移（ΔVth）、迁移率退化（Δμ）和亚阈值摆幅变化（ΔSS），并系统讲解如何从实测数据中提取老化模型参数，为后续的电路级可靠性仿真奠定基础。
上图展示了MOSFET老化模型的整体框架：左侧为三大退化机制（NBTI、PBTI、HCI），中间为对应的核心老化参数（ΔVth、Δμ、ΔSS），右侧为模型构建方法，底部则为从应力实验到模型验证的完整参数提取流程。这张图贯穿了本章的全部核心内容，后续各节将逐一展开。
4.1 核心原理：三大老化参数的物理本质与退化机制
MOSFET老化参数的漂移并非孤立事件，而是源于器件内部微观缺陷的积累。理解这些参数的物理本质，是建立准确老化模型的前提。
（1）阈值电压漂移（ΔVth）：阈值电压是MOSFET开启的关键电压。老化导致的ΔVth主要来源于两个方面：一是界面陷阱（Nit）——Si/SiO2界面处因应力产生的悬挂键和缺陷，它们会俘获载流子，改变有效界面电荷密度；二是氧化层陷阱（Not）——栅氧化层内部（尤其是高k介质如HfO2中）预存的缺陷位点，在电场作用下俘获载流子。对于PMOS器件，NBTI是主导机制，表现为负向栅压应力下ΔVth正向漂移（值增大）；对于NMOS器件，PBTI在高k介质中更为显著。HCI效应则主要在沟道近漏端产生界面损伤，导致局部阈值电压升高。
（2）迁移率退化（Δμ）：迁移率反映了载流子在沟道中的输运能力。老化过程中，界面陷阱和氧化层电荷会增加载流子的库仑散射和表面粗糙度散射，从而降低有效迁移率。迁移率退化通常与ΔVth呈幂律关系，即Δμ/μ0 = -B·(ΔVth)^α，其中α通常在0.5~1.0之间。迁移率退化直接导致饱和电流Idsat下降和跨导gm降低，影响电路的速度性能。
（3）亚阈值摆幅变化（ΔSS）：亚阈值摆幅SS定义为亚阈值区漏电流变化一个数量级所需的栅压变化量，其理想值为60mV/dec（室温）。SS与界面陷阱密度Nit直接相关：SS = (kT/q)·ln10·(1 + Cd/Cox + q·Nit/Cox)。当Nit增加时，SS会退化（增大），意味着器件关断特性变差，静态功耗增加。ΔSS是衡量界面退化程度的重要指标，也是区分BTI与HCI贡献的辅助手段。
知识结构表：三大老化参数对比
| 参数 | 主导机制 | 物理来源 | 对电路影响 | 典型退化量级 |
|---|---|---|---|---|
| ΔVth | NBTI/PBTI/HCI | 界面陷阱+氧化层陷阱 | 时序漂移、噪声容限下降 | 10~100 mV（10年） |
| Δμ | 库仑散射增强 | 界面电荷散射 | 驱动电流下降、速度变慢 | 2%~10% |
| ΔSS | 界面陷阱增加 | Nit增大 | 关断漏电增加、静态功耗上升 | 5~30 mV/dec |
4.2 判断方法与操作流程：老化模型参数提取的标准步骤
老化模型参数提取是一个严谨的工程过程，需要结合应力实验、电学测试和数值拟合。以下是标准化的五步流程：
- 应力实验设计：确定应力条件（Vgs、Vds、温度）和时间点。BTI应力通常采用恒定栅压（如Vgs = Vdd×1.1），HCI应力则需同时施加较高的Vds。温度一般选取85°C、105°C、125°C等加速条件。时间点需覆盖从1秒到10^4秒的多个数量级，以便拟合时间幂律指数。
- 电学特性测试：在每个应力时间点暂停应力，快速测量I-V特性。测量需在极短时间内完成（通常<1ms），以减少恢复效应（BTI的退火现象）。关键测量包括：转移特性曲线（Id-Vgs）、输出特性曲线（Id-Vds）和亚阈值区特性。
- 参数提取：从I-V数据中提取ΔVth、Δμ、ΔSS。ΔVth常用恒定电流法（如Id = 100nA×W/L）或线性外推法提取；Δμ通过比较应力前后的最大跨导gm_max变化来估算；ΔSS则从亚阈值区Id-Vgs曲线的斜率计算。
- 模型拟合：将提取的退化量随时间、电压、温度的变化关系拟合到老化模型中。常用模型为幂律时间依赖：ΔVth(t) = A·(Vgs - Vth0)^n·exp(-Ea/kT)·t^m，其中A为工艺相关常数，n为电压加速指数（通常1~3），Ea为激活能（BTI约0.1~0.2eV），m为时间幂律指数（约0.15~0.25）。
- 验证与校准：将拟合得到的模型参数代入电路仿真器（如HSPICE、Spectre），与实测的环形振荡器频率退化或SRAM静态噪声容限退化数据进行对比，验证模型精度。若偏差超过5%，需重新检查提取流程或调整模型形式。
经验提醒：BTI恢复效应是参数提取的最大干扰源。建议采用“快速测量”技术（如On-the-fly方法），即在应力不中断的情况下直接测量线性区漏电流，从而避免恢复效应。若必须中断应力，应保证测量时间一致，并在数据后处理时进行恢复修正。
4.3 案例拆解：28nm工艺NMOS器件PBTI老化参数提取实战
本案例基于某28nm高k金属栅工艺平台，对NMOS器件进行PBTI老化实验，提取其老化模型参数，并验证模型精度。
背景与约束：测试器件尺寸W/L = 1μm/28nm，栅氧化层为HfO2/SiO2叠层。应力条件：Vgs = 1.2V（工作电压1.0V），Vds = 0V（BTI条件），温度T = 125°C。应力时间点：1s、10s、100s、1000s、10000s。每个时间点暂停应力，在Vds = 50mV线性区测量Id-Vgs曲线，提取ΔVth。
初始现象与数据：在10000s应力后，ΔVth达到约45mV，且随时间呈幂律增长。提取的ΔVth(t)数据如下表所示：
| 应力时间 (s) | ΔVth (mV) | Δgm/gm0 (%) | ΔSS (mV/dec) |
|---|---|---|---|
| 1 | 8.2 | 1.1 | 2.3 |
| 10 | 15.6 | 2.0 | 4.1 |
| 100 | 24.8 | 3.2 | 6.5 |
| 1000 | 34.1 | 4.4 | 8.8 |
| 10000 | 45.3 | 5.8 | 11.2 |
分析判断过程：首先，对ΔVth与时间取双对数坐标，拟合得到斜率m ≈ 0.18，符合PBTI的典型时间指数。其次，ΔSS的变化表明界面陷阱也在增加，但ΔSS/ΔVth比值较小（约0.25 mV/dec per mV），说明氧化层陷阱（Not）主导了ΔVth，而界面陷阱贡献较小。最后，Δgm/gm0与ΔVth的关系近似为Δgm/gm0 ≈ 0.13×ΔVth，表明迁移率退化与阈值电压漂移呈线性相关。
实际操作步骤：
- 使用半导体参数分析仪（如Keysight B1500A）进行应力与测量自动切换，编写测试程序实现“应力-测量-恢复”循环。
- 在每个应力时间点，先快速测量Id-Vgs（Vds=50mV），扫描范围0~1.2V，步长10mV，测量时间<5ms。
- 提取ΔVth：采用恒定电流法（Id = 100nA × W/L），记录应力前后Vgs的差值。
- 提取Δgm：从Id-Vgs曲线计算跨导gm = dId/dVgs，取峰值gm_max，计算退化百分比。
- 提取ΔSS：在亚阈值区（Id从1nA到100nA）拟合log(Id)-Vgs直线，计算斜率倒数即为SS。
- 将数据导入MATLAB或Python，使用最小二乘法拟合ΔVth(t) = A·t^m，得到A = 8.2 mV/s^0.18，m = 0.18。
- 改变应力电压（0.9V、1.0V、1.1V、1.2V）重复实验，拟合电压加速指数n ≈ 2.1。
- 改变温度（85°C、105°C、125°C）重复实验，拟合激活能Ea ≈ 0.12 eV。
结果对比与验证：将提取的模型参数代入老化模型，预测10年（3.15×10^8 s）后的ΔVth ≈ 82mV。将该预测值输入到SRAM单元仿真中，发现静态噪声容限（SNM）退化约12%，与文献报道的28nm工艺数据吻合。同时，将模型外推至不同工作电压（0.8V~1.2V），预测结果与实测的加速老化数据偏差小于5%。
风险与代价：本实验的代价是测试时间较长（10000s应力约3小时），且高温应力可能对器件造成不可逆损伤。此外，快速测量对探针台和测试系统的时序要求极高，需确保测量延迟一致，否则会引入恢复效应误差。
复盘经验：PBTI在NMOS高k器件中主导ΔVth，且氧化层陷阱贡献大于界面陷阱。时间幂律指数m约为0.18，电压加速指数n约为2.1，激活能Ea约为0.12eV。这些参数可作为该工艺节点老化模型的初始值，但需注意不同批次晶圆可能存在工艺波动，建议在多个批次上重复提取。
4.4 常见误区与总结复盘
在MOSFET老化模型参数提取过程中，工程师常陷入以下误区：
- 误区一：忽略BTI恢复效应。若测量延迟过长，ΔVth会被低估30%~50%。解决方法是采用快速测量或On-the-fly方法，并严格记录测量时序。
- 误区二：将ΔVth与Δμ混为一谈。两者虽然相关，但物理机制不同。ΔVth主要影响阈值点，Δμ影响线性区和饱和区的电流斜率。在模型拟合时，应分别提取，避免交叉耦合。
- 误区三：忽视温度激活能Ea的工艺依赖性。不同工艺的Ea差异较大（NBTI约0.1~0.2eV，HCI约-0.1~0.1eV），不能直接套用文献值，必须通过变温实验提取。
- 误区四：模型外推过度。幂律模型在短时间（<10^4s）内拟合良好，但外推到10年寿命时可能偏差较大。建议结合物理模型（如R-D模型）进行交叉验证，并考虑饱和效应（即ΔVth不会无限增长）。
- 误区五：忽略器件尺寸效应。小尺寸器件（如W/L < 100nm）的随机电报噪声（RTN）会干扰ΔVth提取，建议使用较大面积器件进行参数提取，再通过统计模型修正到目标尺寸。
风险提示：老化模型参数提取的最终目的是用于电路级可靠性仿真。若参数提取不准确，可能导致芯片在寿命末期出现时序违规或功能失效。建议在每次流片后，使用测试芯片的环形振荡器或SRAM阵列进行老化校准，确保模型预测与实际硅片退化一致。
总结复盘：本章系统阐述了MOSFET老化模型的三大核心参数——ΔVth、Δμ、ΔSS的物理本质、退化机制和提取方法。通过28nm工艺的PBTI案例，展示了从应力实验设计到模型验证的完整流程。关键要点包括：理解BTI与HCI的微观差异、掌握快速测量技术抑制恢复效应、正确拟合时间/电压/温度依赖关系、以及通过电路级验证确保模型精度。在实际工程中，老化模型参数提取需要结合工艺特性、测试条件和电路需求进行迭代优化，才能为可靠性仿真提供可信的输入。