📚 课程目录第2章 SPICE模型基础

2. SPICE 模型基础：仿真原理、MOSFET 模型层级与参数分类

SPICE（Simulation Program with Integrated Circuit Emphasis）是集成电路设计中最核心的电路仿真工具。理解 SPICE 的仿真原理、MOSFET 模型层级结构以及模型参数的分类，是进行 BSIMProPlus 统计建模与参数提取的前提。本章将系统梳理这些基础知识，为后续的统计建模实操打下坚实的理论基础。

2.1 SPICE 仿真原理

SPICE 的核心任务是对由半导体器件（如 MOSFET、BJT、二极管、电阻、电容等）组成的电路网表进行数值求解，计算节点电压和支路电流。其仿真流程可概括为以下步骤：

网表解析：读取用户编写的电路网表文件（.cir / .sp），识别器件类型、连接关系、模型名称及参数值。
模型评估：根据每个器件的模型名称，调用对应的模型方程（如 BSIM4 的 I-V 和 C-V 方程），代入模型参数和当前节点电压，计算器件的电流、电荷及其对电压的导数（即电导、跨导、电容等）。
矩阵建立：利用改进节点分析法（MNA），将每个器件的贡献（导纳矩阵元素和电流源向量）组装到全局的稀疏矩阵方程 G·v = i 中。
数值求解：对非线性方程组采用 Newton-Raphson 迭代法进行线性化求解；对瞬态分析则结合数值积分（如梯形法、向后欧拉法）进行时域离散。
收敛判断与输出：检查迭代是否收敛，若收敛则输出指定节点的电压或支路电流波形；若不收敛则调整步长或报错。

  关键理解：SPICE 仿真的精度和速度，直接取决于器件模型的物理准确性（模型层级）和数值算法的鲁棒性。BSIMProPlus 的核心工作就是提取高精度的模型参数，使 SPICE 仿真结果与硅片测试数据高度吻合。

2.2 MOSFET 模型层级

MOSFET 模型从早期简单的解析公式发展到现代高度复杂的物理模型，其层级（Level）直接决定了仿真的精度和适用范围。SPICE 中通过 .MODEL 语句的 LEVEL 参数来指定模型类型。

2.2.1 早期模型：Level 1 / 2 / 3

这些模型诞生于 20 世纪 80 年代，适用于长沟道器件（沟道长度 L > 2 μm），计算速度快但物理机制简化，无法准确描述短沟道效应。

模型	核心方程	主要特点	局限性
Level 1	平方律模型	最简单的 MOS 模型，仅考虑一阶沟道长度调制	忽略亚阈值、迁移率退化、速度饱和
Level 2	基于几何的解析模型	引入体电荷效应、迁移率退化、亚阈值电流	方程复杂，收敛性差，短沟道精度不足
Level 3	半经验模型	引入经验参数拟合短沟道效应，收敛性优于 Level 2	物理意义弱，外推能力差

工程提示：在现代先进工艺（如 28nm 及以下）中，Level 1/2/3 模型已完全被淘汰。但在教学或早期工艺（如 0.35μm 以上）的简单电路仿真中，Level 3 模型仍偶有使用。

2.2.2 现代工业标准模型：BSIM3 / BSIM4

BSIM（Berkeley Short-channel IGFET Model）系列模型由加州大学伯克利分校开发，是当前业界最广泛使用的 MOSFET 紧凑模型。BSIMProPlus 的核心任务就是针对 BSIM3/4 模型进行参数提取与统计建模。

BSIM3v3 (Level 8)： 发布于 1996 年，首次系统性地解决了短沟道效应（如 DIBL、Vth 滚降、亚摆幅退化）的建模问题。采用单方程统一模型，覆盖从亚阈值到强反型的所有工作区。参数数量约 100-150 个。
BSIM4 (Level 54)： 发布于 2000 年，是 BSIM3 的全面升级，专为深亚微米（0.13μm 及以下）和纳米级工艺设计。主要改进包括：
- 引入 源/漏寄生电阻（Rds） 的物理模型。
- 改进 噪声模型（1/f 噪声和热噪声）。
- 支持 NQS（非准静态） 效应模型，适用于高频应用。
- 更精确的 栅极隧穿电流（Igate） 模型。
- 参数数量约 200-300 个，灵活性极高。

2.3 模型参数分类

BSIM 模型参数数量庞大（BSIM4 超过 300 个），为了高效地进行参数提取和统计建模，必须对参数进行科学分类。BSIMProPlus 中通常将参数分为以下几类：

分类	典型参数	物理含义	提取优先级
工艺/几何参数	TOX, XJ, LINT, WINT	栅氧厚度、结深、沟道长度/宽度偏移量	高（由工艺文件或测试结构直接给出）
阈值电压参数	VTH0, K1, K2, DVT0, DVT1	长沟道阈值电压、体效应系数、短沟道 Vth 退化	高（直接影响开关特性）
迁移率参数	U0, UA, UB, UC	低场迁移率、垂直电场退化系数	高（决定驱动电流）
输出电导/饱和参数	VSAT, PCLM, PDIBL1, PDIBL2	饱和速度、沟道长度调制、DIBL 效应	高（决定输出阻抗和增益）
寄生参数	RDSW, RDSWMIN, CJ, CJSW	源漏寄生电阻、结电容	中（影响高频和 I/O 特性）
噪声参数	NOIA, NOIB, NOIC, AF, KF	闪烁噪声和热噪声系数	低（仅 RF/模拟电路关注）

统计建模注意事项：在 BSIMProPlus 中进行统计建模时，并非所有参数都需要被设为统计变量。通常只选择对电路性能影响最大、且工艺波动敏感的“核心参数”（如 VTH0, U0, VSAT, RDSW 等）进行统计建模。过多的统计变量会导致模型复杂度过高，且难以收敛。

2.4 SPICE 网表中的模型调用示例

以下是一个典型的 SPICE 网表片段，展示了如何调用 BSIM4 模型（Level 54）并定义 MOSFET 器件：

* BSIM4 模型调用示例


.lib 'bsim4_models.lib' nom





* 定义 NMOS 模型


.MODEL nmos_tt nmos (


+ LEVEL   = 54


+ VERSION = 4.8


+ TOX     = 2.2e-9


+ VTH0    = 0.35


+ U0      = 0.045


+ VSAT    = 1.1e5


+ RDSW    = 200


+ ...     * 其他参数


)





* 定义 PMOS 模型


.MODEL pmos_tt pmos (


+ LEVEL   = 54


+ VERSION = 4.8


+ TOX     = 2.2e-9


+ VTH0    = -0.35


+ U0      = 0.018


+ VSAT    = 0.9e5


+ RDSW    = 250


+ ...     * 其他参数


)





* 实例化 MOSFET


M1 drain gate source bulk nmos_tt W=1u L=0.18u AD=0.5p AS=0.5p





* 仿真控制


VDD vdd 0 DC 1.8V


VG gate 0 DC 0.9V


VD drain 0 DC 0.1V


VS source 0 DC 0V


VB bulk 0 DC 0V





.DC VD 0 1.8 0.01


.PROBE


.END

  BSIMProPlus 角色：上述网表中的模型参数（如 VTH0, U0, VSAT 等）正是 BSIMProPlus 通过优化算法从测试数据中提取出来的。统计建模则是在这些标称值（Typical）基础上，引入统计分布（如高斯分布），生成 Corner 模型或 Monte Carlo 模型，用于评估工艺波动对电路良率的影响。