1. 半导体物理基础回顾:能带理论、载流子浓度、迁移率与电阻率的关系

在深入 SPICE 模型参数提取之前,必须牢固掌握半导体物理的三大基石:能带理论解释固体中电子的允许与禁止状态;载流子浓度决定导电粒子的数量;迁移率衡量载流子在电场下的运动能力。三者通过电导率公式 σ = q (n μn + p μp) 紧密关联,是 SPICE 模型中所有电学特性的物理根源。

SPICE 建模视角: 能带结构决定了阈值电压与漏电流的物理上限;载流子浓度直接对应 MOS 反型层电荷密度;迁移率则是 BSIM 模型中最重要的退化效应参数之一。本节内容将直接服务于后续的模型参数提取与仿真验证。

1.1 能带理论:固体中电子的量子态

在孤立原子中,电子占据分立的能级。当大量原子周期性排列形成晶体时,原子轨道发生交叠,分立能级展宽为能带。相邻能带之间的能量区间称为禁带(Band Gap, Eg),电子无法在其中稳定存在。

  • 价带(Valence Band, Ev): 零度下被电子完全填满的能带,电子被共价键束缚,不参与导电。
  • 导带(Conduction Band, Ec): 价带之上的空带或部分填充带,电子进入导带后可在晶格中自由移动,形成电流。
  • 禁带宽度 Eg = Ec - Ev: 决定本征载流子浓度与器件工作温度范围的关键参数。
关键参数: 硅(Si)在 300K 时 Eg = 1.12 eV;锗(Ge)Eg = 0.67 eV;砷化镓(GaAs)Eg = 1.42 eV。禁带宽度越大,本征载流子浓度越低,器件漏电越小,但开启电压也越高。
能带结构示意图(本征半导体) 导带 (Ec) 禁带 (Eg = Ec - Ev) 价带 (Ev) EF 本征半导体:费米能级 EF 位于禁带中央附近

1.2 载流子浓度:电子与空穴的统计分布

在热平衡状态下,导带电子浓度 n0 与价带空穴浓度 p0 由费米-狄拉克分布与态密度共同决定。对于非简并半导体,可用玻尔兹曼近似简化:

n₀ = N_c · exp[-(E_c - E_F) / (k_B T)]
p₀ = N_v · exp[-(E_F - E_v) / (k_B T)]

其中 Nc 与 Nv 分别为导带与价带的有效态密度,kB 为玻尔兹曼常数,T 为温度。

本征载流子浓度是半导体最重要的本征参数之一:

n_i = √(N_c N_v) · exp[-E_g / (2 k_B T)]

在 300K 下,硅的 ni ≈ 1.5 × 1010 cm-3。掺杂后,多数载流子浓度近似等于掺杂浓度,少数载流子浓度由 ni2 = n0 p0 决定。

常见半导体材料在 300K 下的本征参数
材料 Eg (eV) Nc (cm-3) Nv (cm-3) ni (cm-3)
Si 1.12 2.8 × 1019 1.04 × 1019 1.5 × 1010
Ge 0.67 1.04 × 1019 6.0 × 1018 2.4 × 1013
GaAs 1.42 4.7 × 1017 7.0 × 1018 2.1 × 106
SPICE 参数提取注意: 本征载流子浓度 ni 是温度敏感参数,在 SPICE 模型中通常通过 NI 或 NSUB 参数设定。温度每升高 10°C,ni 约增加一倍,直接导致漏电流指数上升。在提取结漏电模型参数时,必须精确控制温度。

1.3 迁移率:载流子在电场下的运动能力

迁移率(μ)定义为载流子在单位电场强度下的平均漂移速度:vd = μ E。迁移率受多种散射机制限制:

  • 晶格散射(声子散射): 随温度升高而增强,μ ∝ T-3/2。
  • 电离杂质散射: 随掺杂浓度升高而增强,μ ∝ Nd-1/2 或 Na-1/2。
  • 表面粗糙度散射: 在 MOS 反型层中占主导,随垂直电场增强而加剧。
  • 载流子-载流子散射: 在高注入条件下(如双极型器件)不可忽略。

在 SPICE 模型中,迁移率退化效应通过 U0(低场迁移率)、UA、UB、UC 等参数描述。BSIM4 模型使用以下经验公式:

μ_eff = μ₀ / [1 + (UA + UC · Vbseff) · (Vgsteff + 2Vth) / TOXE + UB · (Vgsteff + 2Vth)² / TOXE²]

其中 μ0 为低场迁移率,TOXE 为栅氧化层电学厚度,Vgsteff 为有效栅过驱动电压。

迁移率随掺杂浓度与温度的变化趋势 掺杂浓度 N (cm-3) 迁移率 μ (cm²/V·s) μn (300K) μp (300K) μn (400K) 高 μ 低 μ 迁移率随掺杂浓度升高而降低,随温度升高而降低(晶格散射增强)

1.4 电阻率与电导率:宏观电学特性的物理根源

电阻率 ρ 与电导率 σ 互为倒数,是半导体宏观导电能力的直接度量。对于均匀掺杂的半导体:

σ = q (n μn + p μp)
ρ = 1 / σ = 1 / [q (n μn + p μp)]

其中 q = 1.602 × 10-19 C 为基本电荷量。对于 N 型半导体(n ≫ p),可简化为 ρ ≈ 1 / (q n μn);对于 P 型半导体(p ≫ n),ρ ≈ 1 / (q p μp)。

SPICE 模型中的电阻率应用: 在 SPICE 中,电阻率通过 RHO 参数或薄层电阻 RSH 定义。扩散电阻、阱电阻、多晶硅电阻的建模均基于此公式。参数提取时,通常通过 TLM(传输线模型)测试结构测量薄层电阻,再反推掺杂浓度与迁移率。
不同掺杂浓度下 N 型硅的电阻率(300K)
掺杂浓度 Nd (cm-3) 电子迁移率 μn (cm²/V·s) 电阻率 ρ (Ω·cm) 典型应用
1 × 1014 1350 ≈ 46.3 高阻衬底
1 × 1016 1200 ≈ 0.52 轻掺杂阱
1 × 1018 350 ≈ 0.018 源/漏重掺杂
1 × 1020 100 ≈ 0.0006 多晶硅栅极
电阻率 ρ 与掺杂浓度 N 的关系(N 型 Si, 300K) 掺杂浓度 Nd (cm-3) 电阻率 ρ (Ω·cm) 1014 1015 1016 1017 1018 ρ ∝ 1/(N·μ) 高阻区 低阻区 掺杂浓度越高,电阻率越低;但迁移率退化使曲线斜率在高掺杂区变缓

1.5 电导率公式在 SPICE 模型中的统一表达

在 SPICE 模型中,无论是 MOS 沟道电流、双极型集电极电流还是电阻模型,最终都回归到电导率公式的积分形式。以 MOS 晶体管为例,沟道电流的漂移分量可写为:

IDS = W · μ · Qinv · dV/dy

其中 Qinv 为反型层电荷面密度,由栅压与能带弯曲决定。将 Qinv = Cox (VGS - Vth) 代入并沿沟道积分,即得到经典的萨支唐(Sah)方程:

IDS = (W/L) · μ · Cox · [(VGS - Vth)VDS - VDS²/2]

该方程清晰地展示了能带理论(通过 Vth)、载流子浓度(通过 Cox 与 VGS - Vth)以及迁移率(μ)如何共同决定器件的宏观电学特性。

参数提取实践: 在 SPICE 参数提取流程中,通常先通过 CV 测量提取氧化层厚度与掺杂分布(能带与载流子浓度相关),再通过 IV 测量提取迁移率退化参数与串联电阻。电阻率公式 ρ = 1/(q n μ) 是连接 CV 与 IV 数据的桥梁。

1.6 本章小结:物理量与 SPICE 模型参数的对应关系

半导体物理量与 SPICE 模型参数的映射
物理量 符号 SPICE 模型参数 典型值(Si, 300K)
禁带宽度 Eg EG 1.12 eV
本征载流子浓度 ni NI 1.5 × 1010 cm-3
低场电子迁移率 μn0 U0 ≈ 670 cm²/V·s (NMOS)
低场空穴迁移率 μp0 U0 ≈ 250 cm²/V·s (PMOS)
衬底掺杂浓度 Nsub NSUB 1 × 1015 ~ 1 × 1017 cm-3
薄层电阻 Rsh RSH 1 ~ 1000 Ω/□