2. BSIM6 模型物理基础:表面势模型原理、沟道电荷模型、阈值电压与体效应

BSIM6 是面向先进 CMOS 工艺的紧凑型 MOSFET 模型,其核心物理基础建立在 表面势 (Surface Potential, ψs) 的精确求解之上。与早期基于阈值电压 (Vth) 的模型不同,BSIM6 通过连续、解析的表面势表达式,实现了从亚阈值区到强反型区的平滑过渡,从而更准确地描述沟道电荷、电流以及各类物理效应。本章将深入剖析 BSIM6 的三大物理支柱:表面势模型、沟道电荷模型以及阈值电压与体效应。

2.1 表面势模型原理

表面势 ψs 是半导体表面相对于体硅的静电势,它是所有 MOSFET 物理的核心。BSIM6 采用基于 泊松方程 的解析解,避免了传统模型中分段函数的非物理不连续性。

2.1.1 基本方程与迭代求解

在长沟道近似下,表面势由以下隐式方程决定:

Vgs - Vfb - ψs = γ * √(ψs + φt * exp((ψs - 2*φf - Vch) / φt))

其中:

  • Vgs:栅源电压
  • Vfb:平带电压
  • γ:体效应系数 (γ = √(2*q*εsi*Nsub) / Cox)
  • φf:费米势
  • φt:热电压 (kT/q)
  • Vch:沟道准费米势

BSIM6 通过 牛顿-拉夫森 (Newton-Raphson) 迭代法求解此方程,并在源端 (Vch=0) 和漏端 (Vch=Vds) 分别计算表面势 ψs0 和 ψsL。为了加速仿真,模型提供了高效的初始猜测算法,确保迭代在 3-5 步内收敛。

关键点: BSIM6 的表面势求解是 连续且解析 的,这保证了跨工作区域(亚阈值、线性、饱和)的导数连续性,对于模拟电路(如放大器、混频器)的失真分析至关重要。

2.1.2 表面势与沟道电势的关系

下图展示了沿沟道方向的表面势分布,这是理解沟道电荷和电流的基础。

BSIM6 表面势分布 · 沟道电荷 · 阈值电压与体效应 沿沟道方向的表面势分布(饱和区) 表面势 ψs (V) x=0 x=L 夹断点 ψs0 ψs(x) ψsL 源端 ψs0 由 Vgs 决定,基本不变 Vds 增大,ψsL 变化很小 体偏压对 ψs–Vgs 的影响 ψs (V) Vgs (V) Vbs = 0 V Vbs = -0.5 V Vbs = -1.0 V 负 Vbs 使曲线右移,阈值电压 Vth 增大 ① 表面势求解 Vgs - Vfb - ψs = γ·√(ψs + φt·e^n) 牛顿-拉夫森迭代求解 端点 ψs0 / ψsL,3~5 步收敛 ② 沟道电荷模型 Qinv = -Cox·(Vgs - Vfb - ψs) 多晶硅栅耗尽 + 量子修正 沿沟道积分得漏源电流 Ids ③ 阈值电压与体效应 Vth = Vfb + 2φf + γ·√(2φf - Vbs) K1 / K2 / K3、VBSC 参数化 自动包含 SCE 与 DIBL 效应
物理洞察: 在饱和区,靠近漏端的表面势 ψsL 几乎不再随 Vds 增加而增加(夹断效应),而源端 ψs0 保持恒定。BSIM6 通过精确计算这两个端点的表面势,自然地描述了沟道长度调制效应 (CLM)。

2.2 沟道电荷模型

沟道电荷是连接表面势与电流的核心桥梁。BSIM6 采用 基于电荷的模型 (Charge-Based Model),直接对沟道中的反型层电荷密度 Qinv 进行积分,从而得到漏源电流。

2.2.1 反型层电荷密度

在表面势 ψs 已知的前提下,反型层电荷密度由下式给出:

Qinv(x) = -Cox * (Vgs - Vfb - ψs(x))

其中 Cox 是单位面积栅氧化层电容。BSIM6 进一步考虑了 多晶硅栅耗尽 (Poly Depletion) 和 量子力学效应 (QME) 对 Qinv 的修正,使得模型在超薄栅氧 (EOT < 2nm) 下依然准确。

2.2.2 电荷薄层近似与电流积分

基于漂移-扩散输运理论,漏源电流 Ids 通过对沟道电荷从源端到漏端进行积分得到:

Ids = μ * W / L * ∫ψs0ψsL Qinv(ψs) * dψs

BSIM6 利用解析的 Qinv-ψs 关系,直接计算出闭合形式的电流表达式,避免了数值积分带来的仿真开销。

BSIM6 沟道电荷模型关键参数
参数名 描述 典型值范围 物理效应
TOXE 等效栅氧化层物理厚度 1.0 – 5.0 nm 决定 Cox 和栅控能力
VFB 平带电压 -1.0 – 0.0 V 影响阈值电压和表面势起点
MOB 迁移率模型选择标志 0 / 1 / 2 控制迁移率退化效应
QMFACTOR 量子力学效应修正因子 0.5 – 1.0 修正薄栅氧下的电荷分布
模型配置注意: 在 BSIM6 中,参数 MOB 和 QMFACTOR 的取值会显著影响沟道电荷的精度。对于模拟设计,建议将 MOB 设置为 2(高阶迁移率模型),并开启 QMFACTOR 的自动校准功能。

2.3 阈值电压与体效应

虽然 BSIM6 是表面势模型,但阈值电压 Vth 作为一个重要的宏观参数,仍然用于模型验证和电路设计中的快速估算。BSIM6 中的 Vth 不再是模型的核心变量,而是由表面势推导出的 衍生量。

2.3.1 阈值电压的定义与提取

在 BSIM6 中,阈值电压被定义为 当表面势 ψs 等于 2*φf 时的栅极电压。通过求解表面势方程,可以得到:

Vth = Vfb + 2*φf + γ * √(2*φf - Vbs)

其中 Vbs 是体偏压。这个表达式与经典的长沟道 Vth 公式一致,但 BSIM6 通过表面势的精确求解,自动包含了短沟道效应 (SCE) 和漏致势垒降低 (DIBL) 对 Vth 的影响。

2.3.2 体效应:衬底偏置的影响

体效应描述了衬底偏压 Vbs 对阈值电压和沟道电流的调制作用。BSIM6 通过 体效应对表面势的直接影响 来建模,而不是使用传统的体效应系数 γ。

下图展示了不同 Vbs 下表面势随 Vgs 的变化曲线,清晰地体现了体效应。

不同体偏压下的表面势 vs. 栅压 Vgs (V) ψs (V) 2*φf 判据电平 Vbs = 0 V Vbs = -0.5 V Vbs = -1.0 V 负 Vbs → 曲线右移 阈值电压 Vth 增大

2.3.3 体效应的参数化表示

BSIM6 通过以下参数对体效应进行精细建模:

  • K1, K2:一阶和二阶体效应系数,控制 Vth 随 Vbs 的变化率。
  • K3:短沟道体效应修正系数,用于窄沟道器件。
  • VBSC:体效应从强反型到弱反型的过渡系数。
SPICE 调用示例: 在 BSIM6 模型中,体效应参数通常通过 .MODEL 语句定义。以下是一个典型的 NMOS 模型卡片段:
.MODEL NMOS_BSIM6 NMOS

+ LEVEL = 14

+ TOXE = 1.8E-9

+ VFB = -0.55

+ K1 = 0.45

+ K2 = -0.08

+ K3 = 1.5

+ VBSC = 0.02

2.4 本章小结

BSIM6 的物理基础建立在 表面势的精确求解 之上,通过连续的 ψs 表达式,统一了亚阈值、线性与饱和区的模型描述。沟道电荷模型直接由表面势推导,确保了电流计算的物理一致性。阈值电压作为衍生量,自然包含了体效应、短沟道效应等复杂物理机制。掌握这些物理基础,是正确配置 BSIM6 参数和进行高效 SPICE 仿真的前提。