2. 半导体物理基础:PN结原理、MOSFET结构、载流子输运、阈值电压概念

SPICE 模型的核心是对半导体器件物理行为的数学描述。要准确提取模型参数并验证模型精度,必须深入理解器件内部的工作原理。本章将系统回顾构成现代集成电路的两大基石——PN 结与 MOSFET 的物理基础,重点阐述载流子输运机制与阈值电压的物理本质,为后续的模型参数提取奠定坚实的理论基础。

2.1 PN结原理

PN 结是所有半导体器件中最基本的结构,由 P 型半导体与 N 型半导体通过冶金结接触形成。其核心物理特性包括内建电势、耗尽层与电流-电压特性。

2.1.1 内建电势与耗尽层

当 P 型与 N 型半导体接触时,由于载流子浓度梯度,空穴从 P 区向 N 区扩散,电子从 N 区向 P 区扩散。扩散运动在结区留下固定的电离杂质离子(P 区的负离子与 N 区的正离子),形成空间电荷区,即耗尽层。耗尽层内建电场阻止载流子的进一步扩散,最终达到热平衡状态,形成内建电势 \( V_{bi} \)。

内建电势的表达式为:

V_bi = (k * T / q) * ln(N_A * N_D / n_i^2)

其中,\( k \) 为玻尔兹曼常数,\( T \) 为温度,\( q \) 为电子电荷,\( N_A \) 与 \( N_D \) 分别为受主与施主杂质浓度,\( n_i \) 为本征载流子浓度。

关键物理量: 耗尽层宽度 \( W_d \) 随外加反向偏压增大而展宽,随掺杂浓度升高而变窄。这一特性直接决定了 PN 结的电容行为,是 SPICE 模型中结电容模型(如二极管模型中的 CJ 参数)的物理基础。

2.1.2 电流-电压特性

PN 结的电流-电压关系由 Shockley 理想二极管方程描述:

I = I_S * [exp(q * V / (n * k * T)) - 1]

其中,\( I_S \) 为反向饱和电流,\( n \) 为理想因子(通常在 1 到 2 之间),\( V \) 为外加电压。正向偏置时电流呈指数增长,反向偏置时电流趋近于 \( -I_S \)。

SPICE 建模要点: 实际 PN 结存在串联电阻、大注入效应、产生-复合电流等非理想因素。SPICE 模型通过参数 RS(串联电阻)、ISR(复合电流系数)和 N(发射系数)来修正理想方程,使模型更贴合实测数据。
PN 结能带图与耗尽层结构(热平衡态) P 型区 空穴多子 受主离子 N_A⁻ N 型区 电子多子 施主离子 N_D⁺ 耗尽层 内建电场 空间电荷区 费米能级 E_F(连续) 导带底 E_C(弯曲) 价带顶 E_V(弯曲) V_bi 内建电势

2.2 MOSFET结构与工作原理

金属-氧化物-半导体场效应晶体管(MOSFET)是数字与模拟集成电路的核心器件。其基本结构包含源极、漏极、栅极和衬底四个端子。

2.2.1 基本结构与类型

MOSFET 按导电沟道类型分为 NMOS(N 沟道)和 PMOS(P 沟道)。以 NMOS 为例,在 P 型衬底上形成两个重掺杂 N⁺ 区作为源极和漏极,栅极位于源漏之间的沟道上方,由金属(或多晶硅)与薄栅氧化层构成。

MOSFET 类型对比
参数NMOSPMOS
衬底类型P 型N 型
源/漏掺杂N⁺P⁺
载流子类型电子(迁移率高)空穴(迁移率低)
阈值电压极性正(增强型)负(增强型)
典型应用数字逻辑、高速开关负载管、互补结构

2.2.2 沟道形成与工作区域

当栅极施加正电压 \( V_{GS} \) 时,栅氧化层下方的 P 型衬底表面发生反型,形成 N 型导电沟道,连接源极与漏极。根据偏置条件,MOSFET 工作于三个主要区域:

  • 截止区: \( V_{GS} < V_{TH} \),沟道未形成,漏电流极小(亚阈值泄漏电流除外)。
  • 线性区(三极管区): \( V_{GS} > V_{TH} \) 且 \( V_{DS} < V_{GS} - V_{TH} \),沟道连续,漏电流随 \( V_{DS} \) 线性增加。
  • 饱和区: \( V_{GS} > V_{TH} \) 且 \( V_{DS} \ge V_{GS} - V_{TH} \),沟道在漏端夹断,漏电流基本饱和,由 \( V_{GS} \) 控制。
模型验证注意: 不同工作区域的电流方程形式差异显著。SPICE 模型(如 BSIM4)通过平滑函数实现区域间的连续过渡,避免导数不连续导致的收敛问题。参数提取时必须覆盖所有工作区域,确保模型全局精度。
NMOS 结构剖面图(增强型) P 型衬底(Substrate) 源极 N⁺ 漏极 N⁺ 栅氧化层 SiO₂ 栅极(Gate) N 型沟道(反型层) S D G B

2.3 载流子输运机制

半导体器件中的电流源于载流子(电子与空穴)的定向运动。主要的输运机制包括漂移运动与扩散运动。

2.3.1 漂移运动

在外加电场 \( \vec{E} \) 作用下,载流子获得定向速度,形成漂移电流。漂移电流密度为:

J_drift = q * (n * μ_n + p * μ_p) * E

其中,\( n \) 和 \( p \) 分别为电子与空穴浓度,\( \mu_n \) 和 \( \mu_p \) 分别为电子与空穴迁移率。迁移率受掺杂浓度、温度、晶格散射和表面粗糙度等因素影响,是 SPICE 模型中重要的物理参数。

2.3.2 扩散运动

当载流子浓度存在空间梯度时,载流子从高浓度区向低浓度区扩散,形成扩散电流。扩散电流密度为:

J_diff = q * D_n * (dn/dx) - q * D_p * (dp/dx)

其中,\( D_n \) 和 \( D_p \) 为扩散系数,与迁移率通过爱因斯坦关系关联:\( D/\mu = kT/q \)。

SPICE 模型关联: 在 MOSFET 沟道中,漂移与扩散同时存在。线性区以漂移为主,亚阈值区以扩散为主。BSIM4 模型通过统一的电流方程涵盖两种机制,参数如 U0(低场迁移率)和 VSAT(饱和速度)直接影响模型精度。

2.3.3 迁移率退化效应

在高垂直电场下,沟道载流子受到栅氧化层界面散射,迁移率显著下降。常见的迁移率退化模型为:

μ_eff = μ_0 / [1 + θ * (V_GS - V_TH)]

其中,\( \theta \) 为迁移率退化系数,是 SPICE 模型中的关键参数之一。

载流子输运机制:漂移与扩散 漂移运动 驱动力:外加电场 E 载流子沿电场方向运动 电流密度:J = q·n·μ·E 迁移率 μ 受散射影响 对比 扩散运动 驱动力:浓度梯度 dn/dx 载流子从高浓度→低浓度 电流密度:J = q·D·dn/dx 扩散系数 D = μ·kT/q

2.4 阈值电压概念

阈值电压 \( V_{TH} \) 是 MOSFET 最重要的电学参数之一,定义为形成强反型沟道所需的最小栅极电压。它决定了器件的开启条件,直接影响电路的功耗、速度和噪声容限。

2.4.1 阈值电压的物理组成

阈值电压由多个物理分量构成:

V_TH = V_FB + 2 * φ_F + (sqrt(2 * q * ε_si * N_A * (2 * φ_F + V_SB))) / C_ox

其中:

  • \( V_{FB} \): 平带电压,由功函数差和氧化层固定电荷决定。
  • \( 2\varphi_F \): 强反型所需表面势,\( \varphi_F = (kT/q) \cdot \ln(N_A/n_i) \)。
  • 体效应项: 与衬底偏压 \( V_{SB} \) 相关,体现衬偏调制效应。
  • \( C_{ox} \): 单位面积栅氧化层电容,\( C_{ox} = \varepsilon_{ox}/t_{ox} \)。

2.4.2 影响阈值电压的关键因素

阈值电压影响因素与 SPICE 参数对应
物理因素影响趋势SPICE 模型参数
栅氧化层厚度 \( t_{ox} \) 增大\( V_{TH} \) 增大TOX
衬底掺杂浓度 \( N_A \) 升高\( V_{TH} \) 增大NDEP / NSUB
衬底偏压 \( V_{SB} \) 增大\( V_{TH} \) 增大(体效应)VTH0 / K1 / K2
温度升高\( V_{TH} \) 减小KT1 / KT2
短沟道效应(SCE)\( V_{TH} \) 减小DVT0 / DVT1
窄沟道效应(NCE)\( V_{TH} \) 增大W0 / DVT0W
参数提取实践: 阈值电压的提取通常采用恒定电流法(\( V_{GS} \) 在 \( V_{DS} \) 较小时对应 \( I_D = 100 \, \text{nA} \times W/L \) 时的值)或线性外推法(\( I_D-V_{GS} \) 曲线最大斜率点切线外推至 \( I_D=0 \))。提取结果需结合不同沟道长度和宽度的测试数据,以分离短沟道和窄沟道效应。
阈值电压提取方法:线性外推法 V_GS (V) I_D 切线(最大斜率) V_TH I_D = 0 外推点对应的 V_GS 即为 V_TH

2.5 本章知识结构总结

以下表格归纳了本章核心知识点及其在 SPICE 模型中的对应关系,帮助建立物理原理与模型参数之间的桥梁。

半导体物理基础与 SPICE 模型参数映射
物理概念核心公式/原理关键 SPICE 参数模型验证关注点
PN 结内建电势\( V_{bi} = \frac{kT}{q} \ln\left(\frac{N_A N_D}{n_i^2}\right) \)VJ (结电势)正向 I-V 拟合精度
PN 结耗尽层电容\( C_j = \frac{C_{j0}}{(1 - V/V_J)^M} \)CJO, M, FCC-V 曲线匹配度
MOSFET 沟道电流(线性区)\( I_D = \mu C_{ox} \frac{W}{L} (V_{GS} - V_{TH}) V_{DS} \)U0, VTH0低 V_DS 下 I-V 线性度
MOSFET 沟道电流(饱和区)\( I_D = \frac{1}{2} \mu C_{ox} \frac{W}{L} (V_{GS} - V_{TH})^2 (1 + \lambda V_{DS}) \)U0, VTH0, LAMBDA饱和区输出电导
迁移率退化\( \mu_{eff} = \frac{\mu_0}{1 + \theta (V_{GS} - V_{TH})} \)UA, UB, UC高 V_GS 下电流下降趋势
体效应\( V_{TH} = V_{TH0} + \gamma (\sqrt{2\varphi_F + V_{SB}} - \sqrt{2\varphi_F}) \)K1, K2不同 V_SB 下 V_TH 变化
亚阈值电流\( I_D \propto \exp\left(\frac{V_{GS} - V_{TH}}{n kT/q}\right) \)NFACTOR, CIT亚阈值摆幅 SS 精度
学习路径提示: 本章建立的物理概念是后续章节的基础。第 3 章将深入 PN 结二极管模型与参数提取,第 4 章至第 6 章将系统讲解 MOSFET 模型(Level 1/2/3 到 BSIM4)的方程体系与参数提取流程。建议在学习过程中反复对照本章的物理原理,理解每个模型参数的物理来源。