3. BSIM 模型核心方程
BSIM (Berkeley Short-channel IGFET Model) 系列模型是业界标准的紧凑型 MOSFET 模型。其核心方程精确描述了晶体管从亚阈值区到强反型区、从线性区到饱和区的电学行为。本章深入剖析 BSIM 模型中的四个核心方程：阈值电压模型、迁移率模型、饱和速度模型以及输出电导模型。理解这些方程是进行参数提取与性能优化的基础。
3.1 阈值电压模型 (Vth Model)
阈值电压是 MOSFET 从弱反型过渡到强反型的临界栅压。BSIM 模型中的阈值电压不仅包含长沟道器件的本征值，还引入了短沟道效应 (SCE)、窄沟道效应 (NWE)、衬偏效应 (Body Effect) 以及漏致势垒降低 (DIBL) 的修正。
BSIM4 中阈值电压的基本表达式为：
Vth = VTH0 + K1 * (sqrt(PHI - Vbs) - sqrt(PHI))
- K2 * Vbs
+ K3 * (Tox / (W0 + W)) * PHI
- DVT0 * exp(-DVT1 * Leff / (2 * Ltw))
* (Vbi - PHI)
- (DVT2 * exp(-DVT1 * Leff / (2 * Ltw)) + (ETA0 + ETAB * Vbs)) * Vds其中关键参数及其物理含义如下表所示：
| 参数 | 物理含义 | 典型影响 |
|---|---|---|
| VTH0 | 长沟道、零偏置下的本征阈值电压 | 决定器件开关的基本电压点 |
| K1, K2 | 体效应系数 (Body Effect Coefficients) | 描述衬底偏压对 Vth 的调制 |
| DVT0, DVT1 | 短沟道效应系数 (SCE Coefficients) | 沟道长度减小时，Vth 下降的程度 |
| ETA0, ETAB | DIBL 效应系数 | 漏极电压升高导致 Vth 降低 |
| K3 | 窄沟道效应系数 (NWE Coefficient) | 沟道宽度减小时，Vth 的变化 |
Vds = 50mV (线性区) 和 Vds = Vdd (饱和区) 两组数据。通过对比不同沟道长度 (L) 和不同衬底偏压 (Vbs) 下的 Vth 变化，可以分别提取 SCE 参数和体效应参数。
下图展示了短沟道效应和 DIBL 效应对阈值电压的影响机制：
3.2 迁移率模型 (Mobility Model)
载流子迁移率 (μ) 直接影响沟道电流的大小。在 BSIM 模型中，迁移率受到垂直电场、横向电场、掺杂浓度以及温度等多种因素的退化。BSIM4 使用一个统一的迁移率退化模型，其核心是 有效迁移率 (μeff) 的概念。
迁移率退化的主要机制包括：
- 声子散射 (Phonon Scattering)： 由晶格热振动引起，随温度升高而加剧。
- 库仑散射 (Coulomb Scattering)： 由电离杂质和界面电荷引起，在低垂直电场下占主导。
- 表面粗糙度散射 (Surface Roughness Scattering)： 由 Si-SiO₂ 界面的不平整引起，在高垂直电场下占主导。
BSIM4 中迁移率模型的基本形式为：
μeff = μ0 / (1 + (UA + UC * Vbs) * Eeff + UB * Eeff²)其中，Eeff 是有效垂直电场，μ0 是低场迁移率，UA、UB、UC 是迁移率退化系数。
| 参数 | 物理含义 | 主导散射机制 |
|---|---|---|
| μ0 | 低场迁移率 (Low-field Mobility) | 基础值，受掺杂和温度影响 |
| UA | 一阶迁移率退化系数 | 声子散射 (与 Eeff 线性相关) |
| UB | 二阶迁移率退化系数 | 表面粗糙度散射 (与 Eeff² 相关) |
| UC | 体偏相关的迁移率退化系数 | 库仑散射 (受耗尽层电荷影响) |
UA 和 UB 同时拟合所有栅压下的电流。通常，UA 主导中等栅压区 (Vgs 略高于 Vth)，而 UB 主导高栅压区 (Vgs 接近 Vdd)。建议分区域进行拟合。
下图展示了不同散射机制对迁移率的综合影响：
3.3 饱和速度模型 (Velocity Saturation Model)
当沟道中的横向电场 (Ey) 达到临界值时，载流子的漂移速度不再随电场线性增加，而是趋于饱和。这一现象称为 速度饱和 (Velocity Saturation)，它是短沟道器件中电流饱和的主要原因，也是区分长沟道和短沟道行为的关键。
BSIM 模型使用一个分段或平滑函数来描述载流子速度与电场的关系：
v = (μeff * Ey) / (1 + (Ey / Esat)^BEXP)^(1/BEXP)其中：
v是载流子漂移速度。μeff是有效迁移率。Ey是沿沟道方向的横向电场。Esat是饱和电场，定义为Esat = 2 * Vsat / μeff，其中Vsat是饱和速度。BEXP是速度饱和的过渡指数，通常取 1 或 2。
Ey << Esat 时，v ≈ μeff * Ey (线性区)；当 Ey >> Esat 时，v ≈ Vsat (饱和区)。速度饱和导致短沟道器件的饱和电流 (Idsat) 与沟道长度 (L) 成反比，而非长沟道器件的平方关系。
速度饱和对 I-V 特性的影响可以通过以下对比表理解：
| 特性 | 长沟道器件 (无速度饱和) | 短沟道器件 (速度饱和主导) |
|---|---|---|
| 饱和电流 Idsat | ∝ (Vgs - Vth)² / L | ∝ Cox * W * Vsat * (Vgs - Vth) |
| 跨导 gm | ∝ (Vgs - Vth) / L | ∝ Cox * W * Vsat (趋于常数) |
| 饱和电压 Vdsat | ≈ Vgs - Vth | ≈ L * Esat (更小) |
下图对比了有无速度饱和效应下的 I-V 特性曲线：
3.4 输出电导模型 (Output Conductance Model)
输出电导 (gds) 定义为 gds = ∂Ids / ∂Vds，它反映了漏源电压对沟道电流的调制能力。在理想情况下，饱和区的输出电导应为零，但实际器件中存在多种效应导致 gds 非零，主要包括：
- 沟道长度调制 (CLM)： 随着 Vds 增加，漏极耗尽区向源极扩展，有效沟道长度缩短，导致电流增加。
- DIBL： 高 Vds 降低阈值电压，间接增加电流。
- 衬底电流诱导体效应 (SCBE)： 高 Vds 下，碰撞电离产生的衬底电流改变体电势，调制 Vth。
BSIM4 的输出电导模型通过一个统一的表达式来捕捉这些效应：
gds = gds0 + Δgds_CLM + Δgds_DIBL + Δgds_SCBE其中各分量的物理机制和关键参数如下：
| 分量 | 物理机制 | 关键 BSIM 参数 | 对 gds 的贡献 |
|---|---|---|---|
| gds0 | 本征输出电导 (理想情况) | — | 基础值，通常很小 |
| Δgds_CLM | 沟道长度调制 | PCLM | 随 Vds 增加而增加，在长沟道中显著 |
| Δgds_DIBL | 漏致势垒降低 | PDIBL1, PDIBL2 | 在短沟道中占主导，与 L 成反比 |
| Δgds_SCBE | 衬底电流诱导体效应 | PSCBE1, PSCBE2 | 在高 Vds 下突然增大，导致 gds 回升 |
gds vs. Vds 曲线。在中等 Vds 区域，gds 的下降趋势主要由 CLM 和 DIBL 决定；而在高 Vds 区域，gds 的回升则主要由 SCBE 引起。通过观察 gds 曲线的“U”形或“V”形特征，可以判断各效应的相对强弱。
下图展示了输出电导随 Vds 的变化及各效应的贡献区域：