2. SPICE 模型基础回顾：MOSFET / BJT / Diode 等基本器件模型方程，模型参数物理意义，模型层级（Level 1/2/3/BSIM 等）
SPICE 模型是连接半导体物理与电路仿真的数学桥梁。本章系统回顾三大核心器件（MOSFET、BJT、Diode）的基础模型方程、关键参数的物理含义以及模型层级演进，为后续实测数据校准与自动化提取流程奠定理论基础。
学习目标：掌握 Level 1/2/3 与 BSIM 系列的核心差异；理解模型参数与工艺/物理量的对应关系；能够根据应用场景选择合适的模型层级。
2.1 二极管（Diode）基础模型
二极管是最简单的半导体器件，其 SPICE 模型基于 Shockley 方程，并加入寄生效应与高阶修正。
2.1.1 核心方程
理想二极管电流方程为：
I_D = IS * (exp(V_D / (N * V_T)) - 1)其中：
- IS：饱和电流（A），与结面积、掺杂浓度相关
- N：发射系数（理想因子），通常在 1.0 ~ 2.0 之间
- V_T：热电压，V_T = kT/q ≈ 0.02585 V @ 300K
2.1.2 关键模型参数
| 参数名 | SPICE 关键字 | 物理意义 | 典型值范围 |
|---|---|---|---|
| 饱和电流 | IS | 反向偏置漏电流 | 1e-14 ~ 1e-12 A |
| 发射系数 | N | 非理想性因子 | 1.0 ~ 2.0 |
| 串联电阻 | RS | 接触与体电阻 | 0.1 ~ 10 Ω |
| 结电容 | CJO | 零偏结电容 | 0.1 ~ 10 pF |
| 击穿电压 | BV | 反向击穿电压 | 5 ~ 100 V |
校准要点：实测 I-V 曲线中，低电流区（< 1 μA）主要校准 IS 和 N；高电流区（> 1 mA）主要校准 RS；C-V 曲线用于提取 CJO 和 VJ（内建电势）。
2.2 BJT（双极型晶体管）基础模型
BJT 模型以 Gummel-Poon 模型为基准，SPICE 中常用的是 GP 模型及其简化版本。
2.2.1 Gummel-Poon 核心方程
集电极电流（正向有源区）：
I_C = IS * exp(V_BE / (NF * V_T)) * (1 + V_CE / VAF)基极电流：
I_B = IS / BF * exp(V_BE / (NF * V_T)) + ISE * exp(V_BE / (NE * V_T))2.2.2 关键模型参数
| 参数名 | SPICE 关键字 | 物理意义 | 典型值范围 |
|---|---|---|---|
| 饱和电流 | IS | 传输饱和电流 | 1e-16 ~ 1e-14 A |
| 正向电流增益 | BF | 最大 β_F | 50 ~ 500 |
| Early 电压 | VAF | 基区宽度调制效应 | 20 ~ 200 V |
| 正向渡越时间 | TF | 基区少子渡越时间 | 1 ~ 100 ps |
| 发射极串联电阻 | RE | 发射极接触电阻 | 0.1 ~ 5 Ω |
常见误区：BF 并非恒定值，实测 β 随 I_C 和 V_CE 变化。Gummel-Poon 模型通过 ISE、NE、IKF 等参数描述 β 随电流的变化，校准时应覆盖低、中、高三个电流区。
2.3 MOSFET 模型层级演进
MOSFET 模型从简单的物理解析模型发展到高度精确的表面势基模型，层级划分反映了精度与计算复杂度的权衡。
2.3.1 模型层级对比
| 模型层级 | 核心方程基础 | 关键物理效应 | 参数数量 | 适用工艺节点 |
|---|---|---|---|---|
| Level 1 | Shichman-Hodges 平方律 | 仅沟道长度调制 | ~15 | > 5 μm |
| Level 2 | 几何解析模型 | 体效应、亚阈值导电 | ~30 | 2 ~ 5 μm |
| Level 3 | 半经验半解析 | 短沟道效应、窄沟道效应 | ~40 | 1 ~ 2 μm |
| BSIM3v3 | 阈值电压基经验模型 | 迁移率退化、速度饱和、DIBL | ~100 | 0.35 ~ 0.13 μm |
| BSIM4 | 阈值电压基 + 源/漏寄生 | STI 应力、栅极漏电、NQS | ~200 | 90 nm ~ 28 nm |
| BSIM-CMG | 表面势基 | FinFET 多栅、量子效应 | ~300 | 14 nm 及以下 |
2.3.2 Level 1 模型方程（基础参考）
线性区（V_DS < V_GS - V_TH）：
I_DS = KP * (W/L) * [(V_GS - V_TH) * V_DS - 0.5 * V_DS^2] * (1 + LAMBDA * V_DS)饱和区（V_DS ≥ V_GS - V_TH）：
I_DS = 0.5 * KP * (W/L) * (V_GS - V_TH)^2 * (1 + LAMBDA * V_DS)2.3.3 关键参数物理意义
| 参数 | SPICE 关键字 | 物理意义 | 提取来源 |
|---|---|---|---|
| 阈值电压 | VTO | 零偏阈值电压 | I_D-V_GS 线性外推 |
| 跨导参数 | KP | μ₀ * C_ox | I_D-V_GS 斜率 |
| 沟道长度调制 | LAMBDA | 饱和区输出电导 | I_D-V_DS 饱和区斜率 |
| 体效应系数 | GAMMA | 衬底偏置对 V_TH 的影响 | V_TH vs V_BS 曲线 |
| 表面势 | PHI | 强反型表面势 (2φ_F) | 工艺参数计算 |
2.4 BSIM 系列模型核心特性
BSIM（Berkeley Short-channel IGFET Model）系列是工业界最广泛使用的 MOSFET 模型族，其核心优势在于对短沟道效应和寄生效应的精确建模。
2.4.1 BSIM3v3 与 BSIM4 关键差异
| 特性 | BSIM3v3 | BSIM4 |
|---|---|---|
| 栅极漏电模型 | 无 | 直接隧穿、Fowler-Nordheim |
| 源/漏寄生电阻 | 简单 RDS | 偏置依赖 RDSW |
| STI 应力效应 | 无 | SA, SB, SD 参数 |
| NQS 非准静态 | ELM 选项 | 内置 NQS 模型 |
| 噪声模型 | 简单 1/f | 统一噪声模型 |
| 自热效应 | 无 | SHMODEL 选项 |
2.4.2 BSIM4 核心方程框架
阈值电压模型：
V_TH = VTH0 + K1 * (sqrt(PHI - V_BS) - sqrt(PHI)) - K2 * V_BS
+ DVT0 * (exp(-DVT1 * L_eff / (2 * L_t)) + 2 * exp(-DVT1 * L_eff / L_t)) * (PHI - V_BS)
- DVT2 * V_DS其中 DIBL（漏致势垒降低）效应通过 DVT2 项显式建模。
自动化提取优势：BSIM4 参数数量约 200 个，手动提取需要 2-4 周。通过自动化流程，可将提取周期压缩至 2-3 天，且参数相关性分析更系统化。
2.5 模型选择决策流程
在实际项目中，模型层级的选择直接影响仿真精度与收敛性。以下决策流程可作为参考：
- 工艺节点评估：若 ≤ 28 nm，直接选择 BSIM-CMG（FinFET）或 BSIM4（Planar）；若 > 0.5 μm，Level 3 即可满足需求。
- 应用场景判断：模拟/RF 电路需要高精度 BSIM4/BSIM-CMG；数字电路标准单元库可使用 BSIM3v3 以提升仿真速度。
- 参数提取资源：若缺乏完整测试数据，可先用 Level 3 快速建模，后续再升级至 BSIM 系列。
- EDA 工具兼容性：确认目标仿真器支持的模型版本（如 HSPICE 支持 BSIM4.8.0，Spectre 支持 BSIM-CMG 110.0）。