2. 半导体物理基础回顾：能带理论、载流子统计、漂移扩散方程、PN结基础
器件仿真与模型验证的底层逻辑建立在半导体物理的三大支柱之上：能带理论解释材料特性，载流子统计描述电荷分布，漂移扩散方程刻画输运过程。PN结作为最基础的结构单元，是验证上述理论正确性的试金石。本节将系统回顾这些核心概念，为后续的仿真建模提供理论锚点。
2.1 能带理论：从原子轨道到能带结构
孤立原子的电子处于分立的能级上。当大量原子聚集形成晶体时，原子轨道发生交叠，分立能级展宽为能带。相邻能带之间的区域称为禁带（Band Gap），其宽度 \(E_g\) 决定了材料的基本电学性质。
关键概念： 导带底 \(E_C\)、价带顶 \(E_V\)、禁带宽度 \(E_g = E_C - E_V\)。费米能级 \(E_F\) 是电子占据概率为 1/2 的能级，是统计分布的参考基准。
仿真提示： 在 Sentaurus 或 Silvaco 中，能带参数（如 \(E_g\)、电子亲和能 \(\chi\)）直接输入材料模型文件。温度依赖关系通常采用 Varshni 公式：\(E_g(T) = E_g(0) - \frac{\alpha T^2}{T + \beta}\)。
2.2 载流子统计：费米-狄拉克分布与玻尔兹曼近似
电子在能级上的分布遵循费米-狄拉克统计。在非简并半导体（\(E_C - E_F \gg kT\)）中，可简化为玻尔兹曼近似，得到电子和空穴浓度的解析表达式：
| 载流子类型 | 浓度表达式 | 有效状态密度 |
|---|---|---|
| 电子 \(n_0\) | \(n_0 = N_C \exp\left(-\frac{E_C - E_F}{kT}\right)\) | \(N_C = 2\left(\frac{2\pi m_n^* kT}{h^2}\right)^{3/2}\) |
| 空穴 \(p_0\) | \(p_0 = N_V \exp\left(-\frac{E_F - E_V}{kT}\right)\) | \(N_V = 2\left(\frac{2\pi m_p^* kT}{h^2}\right)^{3/2}\) |
仿真陷阱： 当掺杂浓度超过 \(10^{19} \text{cm}^{-3}\) 时，玻尔兹曼近似失效，必须使用完整的费米-狄拉克积分。此时需在仿真器中开启“Fermi-Dirac Statistics”模型，否则会严重高估载流子浓度。
质量作用定律：在热平衡下，\(n_0 p_0 = n_i^2\)，其中 \(n_i\) 为本征载流子浓度。对于硅，\(n_i \approx 1.5 \times 10^{10} \text{cm}^{-3}\)（300K）。
2.3 漂移扩散方程：电流的两种驱动力
半导体中的电流由两种物理机制驱动：漂移（电场驱动）和扩散（浓度梯度驱动）。总电流密度方程为：
\[ J_n = q \mu_n n E + q D_n \nabla n \] \[ J_p = q \mu_p p E - q D_p \nabla p \]其中，\(\mu\) 为迁移率，\(D\) 为扩散系数，通过爱因斯坦关系关联：\(D = \frac{kT}{q} \mu\)。
| 输运机制 | 驱动力 | 方向 | 典型应用场景 |
|---|---|---|---|
| 漂移电流 | 电场 \(E\) | 电子逆电场方向，空穴顺电场方向 | MOSFET 沟道导通、电阻区 |
| 扩散电流 | 浓度梯度 \(\nabla n\) 或 \(\nabla p\) | 从高浓度区指向低浓度区 | PN 结耗尽区边缘、双极晶体管基区 |
仿真实现： 在 TCAD 仿真中，漂移扩散方程与泊松方程耦合求解。泊松方程 \(\nabla^2 \psi = -\frac{q}{\epsilon}(p - n + N_D^+ - N_A^-)\) 提供电场分布，漂移扩散方程更新载流子浓度。两者迭代至收敛。
2.4 PN结基础：从冶金结到耗尽层
PN结是半导体器件中最基本的结结构。当P型与N型半导体接触时，由于载流子浓度差，发生扩散运动，在界面附近形成空间电荷区（耗尽层），并建立内建电势 \(V_{bi}\)。
2.4.1 平衡态PN结
- 内建电势：\(V_{bi} = \frac{kT}{q} \ln\left(\frac{N_A N_D}{n_i^2}\right)\)
- 耗尽层宽度：\(W = \sqrt{\frac{2\epsilon}{q} \left(\frac{1}{N_A} + \frac{1}{N_D}\right) V_{bi}}\)
- 最大电场：\(E_{max} = \frac{2V_{bi}}{W}\)
2.4.2 正偏与反偏特性
| 偏置状态 | 耗尽层变化 | 电流机制 | I-V 关系 |
|---|---|---|---|
| 零偏 | 平衡宽度 \(W_0\) | 漂移 = 扩散，净电流为零 | \(I = 0\) |
| 正偏 (\(V_F\)) | 耗尽层变窄 | 扩散电流主导 | \(I = I_S \left(e^{qV_F/kT} - 1\right)\) |
| 反偏 (\(V_R\)) | 耗尽层展宽 | 漂移电流主导（漏电流） | \(I \approx -I_S\)（饱和） |
模型验证要点： 在仿真中验证PN结时，需检查：(1) 零偏时内建电势是否与理论值一致；(2) 正偏I-V曲线斜率是否为 \(q/kT\)（约 60 mV/decade）；(3) 反偏击穿电压是否与掺杂浓度和结深匹配。
2.5 关键参数与仿真输入对照表
以下表格总结了器件仿真中需要输入的半导体物理参数，以及对应的典型值（以硅为例，300K）：
| 参数 | 符号 | 典型值 | 仿真输入位置 |
|---|---|---|---|
| 禁带宽度 | \(E_g\) | 1.12 eV | 材料定义文件 |
| 本征载流子浓度 | \(n_i\) | \(1.5 \times 10^{10} \text{cm}^{-3}\) | 自动计算（基于 \(E_g\) 和有效质量） |
| 电子迁移率 | \(\mu_n\) | 1350 \(\text{cm}^2/\text{V·s}\) | 迁移率模型参数 |
| 空穴迁移率 | \(\mu_p\) | 480 \(\text{cm}^2/\text{V·s}\) | 迁移率模型参数 |
| 相对介电常数 | \(\epsilon_r\) | 11.7 | 材料定义文件 |
| 电子亲和能 | \(\chi\) | 4.05 eV | 材料定义文件（用于异质结） |