1. 半导体物理基础:能带理论、载流子浓度、迁移率与电阻率、PN结原理

深亚微米CMOS工艺的器件建模,其根基在于半导体物理。BSIM模型中的阈值电压、迁移率退化、漏致势垒降低(DIBL)等核心效应,均源于能带结构、载流子输运与PN结行为。本章从四个基石概念出发,建立从物理原理到SPICE模型参数的映射关系。

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

在孤立原子中,电子占据分立的能级。当大量原子聚集形成晶体时,原子轨道发生交叠,分立能级展宽为允带,允带之间为禁带(Band Gap, \(E_g\))。对于半导体,价带(Valence Band, \(E_v\))被电子填满,导带(Conduction Band, \(E_c\))为空,两者之间为禁带。硅在300K下的禁带宽度 \(E_g = 1.12\ \text{eV}\)。

关键概念:费米能级 \(E_F\) 是电子占据概率为1/2的能级。本征半导体中 \(E_F\) 位于禁带中央;N型掺杂使 \(E_F\) 向导带靠近;P型掺杂使 \(E_F\) 向价带靠近。
能带结构示意图(本征半导体) 导带 (E_c) 禁带 (E_g = 1.12 eV) 价带 (E_v) E_F E_c E_v

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

半导体中的导电载流子包括电子(导带中)和空穴(价带中)。在热平衡下,载流子浓度由费米-狄拉克分布与态密度决定。对于非简并半导体,可用玻尔兹曼近似简化:

参数符号公式(300K)典型值(Si)
导带有效态密度\(N_c\)\(2(2\pi m_n^* kT/h^2)^{3/2}\)\(2.8 \times 10^{19}\ \text{cm}^{-3}\)
价带有效态密度\(N_v\)\(2(2\pi m_p^* kT/h^2)^{3/2}\)\(1.04 \times 10^{19}\ \text{cm}^{-3}\)
本征载流子浓度\(n_i\)\(\sqrt{N_c N_v} \exp(-E_g/2kT)\)\(1.5 \times 10^{10}\ \text{cm}^{-3}\)
电子浓度(N型)\(n_0\)\(N_D\)(完全电离)\(10^{15} \sim 10^{18}\ \text{cm}^{-3}\)
空穴浓度(P型)\(p_0\)\(N_A\)(完全电离)\(10^{15} \sim 10^{18}\ \text{cm}^{-3}\)
BSIM建模关联:阈值电压 \(V_{th}\) 与沟道掺杂浓度 \(N_{ch}\) 直接相关。\(N_{ch}\) 越高,\(V_{th}\) 越大,但载流子迁移率也会因杂质散射而下降。BSIM4中通过 VTH0 和 K1、K2 参数体现这一关系。
载流子浓度与掺杂浓度关系(300K) N型半导体 n₀ ≈ N_D p₀ = n_i² / N_D 多数载流子:电子 少数载流子:空穴 P型半导体 p₀ ≈ N_A n₀ = n_i² / N_A 多数载流子:空穴 少数载流子:电子 n·p = n_i²

1.3 迁移率与电阻率:载流子输运特性

迁移率(\(\mu\))定义为单位电场下载流子的平均漂移速度:\(\mu = v_d / E\),单位为 \(\text{cm}^2/(\text{V·s})\)。电阻率(\(\rho\))与迁移率的关系为:

\[ \rho = \frac{1}{q(n\mu_n + p\mu_p)} \]

其中 \(q\) 为电子电荷(\(1.6 \times 10^{-19}\ \text{C}\))。对于N型硅,\(n \gg p\),故 \(\rho \approx 1/(q n \mu_n)\)。

散射机制对迁移率的影响温度依赖性BSIM建模参数
晶格散射(声子)高温下主导,\(\mu \propto T^{-3/2}\)负温度系数MU0、KT1
电离杂质散射高掺杂下显著,\(\mu \propto N^{-1/2}\)弱正温度系数NDEP、NSD
表面粗糙度散射高垂直电场下主导(MOSFET沟道)弱UA、UB、UC
深亚微米注意:在短沟道MOSFET中,迁移率受垂直电场(栅压)和横向电场(漏压)共同调制。BSIM4使用 mobMod 选择迁移率模型,默认采用 mobMod = 1(基于垂直电场的通用迁移率退化模型)。
迁移率与电阻率关系(N型硅,300K) 迁移率 \(\mu_n\) vs 掺杂 \(N_D\) 低掺杂:\(\mu_n \approx 1350\ \text{cm}^2/\text{V·s}\) 中掺杂:\(\mu_n \approx 500\ \text{cm}^2/\text{V·s}\) 高掺杂:\(\mu_n \approx 100\ \text{cm}^2/\text{V·s}\) 掺杂↑ → 杂质散射↑ → 迁移率↓ 电阻率 \(\rho\) vs 掺杂 \(N_D\) 轻掺杂:\(\rho \approx 10\ \Omega\cdot\text{cm}\) 中等掺杂:\(\rho \approx 0.1\ \Omega\cdot\text{cm}\) 重掺杂:\(\rho \approx 0.001\ \Omega\cdot\text{cm}\) 掺杂↑ → 电阻率↓(导电性↑) \(\rho = 1/(q n \mu_n)\)

1.4 PN结原理:耗尽区与内建电势

PN结是CMOS器件的基本结构单元(源/漏与衬底、阱与衬底)。当P型与N型半导体接触时,多数载流子相互扩散,在界面附近形成耗尽区(空间电荷区),并建立内建电势 \(V_{bi}\):

\[ V_{bi} = \frac{kT}{q} \ln\left(\frac{N_A N_D}{n_i^2}\right) \]

对于硅PN结,典型 \(V_{bi}\) 在0.6V~0.9V之间。耗尽区宽度 \(W_{dep}\) 随反向偏压增大而展宽:

\[ W_{dep} = \sqrt{\frac{2\varepsilon_s}{q} \left(\frac{1}{N_A} + \frac{1}{N_D}\right) (V_{bi} + V_R)} \]
BSIM建模关联:源/漏与衬底之间的PN结电容(CJ、CJSW)直接影响电路速度。BSIM4通过 XJ(结深)、VJ(内建电势)、MJ(梯度因子)等参数建模结电容随偏压的变化。
PN结结构:耗尽区与能带弯曲 P型区 空穴:多数载流子 N_A = 10¹⁶ cm⁻³ N型区 电子:多数载流子 N_D = 10¹⁶ cm⁻³ 耗尽区 空间电荷 无自由载流子 内建电势 V_bi ≈ 0.7V E_c (P侧) E_c (N侧) 能带弯曲:势垒阻止多数载流子继续扩散

1.5 从物理到SPICE:BSIM模型参数映射

上述半导体物理量直接对应BSIM4模型中的核心参数。下表列出关键映射关系:

物理量BSIM4参数物理含义典型值(45nm工艺)
禁带宽度 \(E_g\)EG禁带宽度,影响本征载流子浓度1.12 eV
沟道掺杂浓度 \(N_{ch}\)NDEP沟道耗尽区掺杂浓度\(2 \times 10^{18}\ \text{cm}^{-3}\)
迁移率 \(\mu_0\)MU0低场迁移率300 \(\text{cm}^2/\text{V·s}\) (NMOS)
内建电势 \(V_{bi}\)VJ源/漏结内建电势0.7 V
结深 \(X_j\)XJ源/漏扩散结深30 nm
耗尽区宽度因子MJ结电容梯度因子0.5
建模实践提示:在SPICE中提取BSIM参数时,首先通过工艺文档获取掺杂浓度、结深、栅氧化层厚度等物理参数,作为初始值(init 参数)。然后通过优化算法拟合I-V和C-V测试数据,调整 VTH0、U0、VSAT 等关键参数。物理参数与模型参数的一致性越高,模型的外推预测能力越强。
从半导体物理到BSIM模型参数映射流程 半导体物理参数 \(N_A, N_D, \mu, E_g, X_j\) BSIM4模型参数 NDEP, MU0, VJ, XJ SPICE网表 .MODEL NMOS BSIM4 物理参数 → 模型初始值 → 优化拟合 → 最终BSIM模型 物理一致性越高,模型的外推预测能力越强