2. 理想二极管方程：肖克利方程推导、热电压、反向饱和电流
理想二极管方程是 PN 结直流特性的数学基础，由肖克利（Shockley）于 1949 年基于扩散理论推导得出。该方程描述了在稳态、小注入、无产生-复合、无串联电阻等理想条件下，PN 结电流与外加电压之间的解析关系。理解其推导过程、热电压的物理意义以及反向饱和电流的构成，是后续进行 SPICE 参数提取与拟合的前提。
2.1 肖克利方程推导
推导基于以下核心假设：
- 耗尽层近似：空间电荷区边界陡峭，耗尽区外为电中性区。
- 玻尔兹曼统计：载流子浓度服从玻尔兹曼分布。
- 小注入条件：注入的少数载流子浓度远小于平衡多数载流子浓度。
- 无产生-复合：耗尽区内无载流子产生与复合，电流完全由扩散电流主导。
推导的关键步骤：
- 边界条件：在外加偏压 \(V\) 下，耗尽层边界处的少数载流子浓度由玻尔兹曼关系给出： \[ n_p(0) = n_{p0} e^{qV/kT}, \quad p_n(0) = p_{n0} e^{qV/kT} \] 其中 \(n_{p0}\) 和 \(p_{n0}\) 分别为平衡时 P 区的电子浓度和 N 区的空穴浓度。
- 扩散方程：在中性区，少数载流子满足稳态扩散方程： \[ \frac{d^2 \Delta n_p}{dx^2} = \frac{\Delta n_p}{L_n^2} \] 其中 \(\Delta n_p = n_p - n_{p0}\)，\(L_n\) 为电子扩散长度。
- 求解电流密度：求解扩散方程并代入边界条件，得到 P 区电子扩散电流密度和 N 区空穴扩散电流密度，求和后得到总电流密度： \[ J = J_s \left( e^{qV/kT} - 1 \right) \]
最终得到肖克利方程（理想二极管方程）：
肖克利方程（理想二极管方程）
\[ I = I_S \left( e^{\frac{V}{n V_T}} - 1 \right) \] 其中：
\(I\) —— 通过 PN 结的电流 (A)
\(I_S\) —— 反向饱和电流 (A)
\(V\) —— 外加偏压 (V)
\(n\) —— 理想因子（理想情况下 \(n=1\)）
\(V_T\) —— 热电压 (V)
\[ I = I_S \left( e^{\frac{V}{n V_T}} - 1 \right) \] 其中：
\(I\) —— 通过 PN 结的电流 (A)
\(I_S\) —— 反向饱和电流 (A)
\(V\) —— 外加偏压 (V)
\(n\) —— 理想因子（理想情况下 \(n=1\)）
\(V_T\) —— 热电压 (V)
💡 推导要点：肖克利方程的本质是扩散电流与边界载流子浓度的指数关系。在正向偏压下，\(e^{V/(nV_T)} \gg 1\)，电流随电压指数增长；在反向偏压下，\(e^{V/(nV_T)} \to 0\)，电流趋近于 \(-I_S\)。
2.2 热电压 \(V_T\)
热电压是半导体物理中的一个基本参数，其表达式为：
\[ V_T = \frac{kT}{q} \]其中：
- \(k = 1.380649 \times 10^{-23} \, \text{J/K}\) —— 玻尔兹曼常数
- \(T\) —— 温度 (K)
- \(q = 1.602176634 \times 10^{-19} \, \text{C}\) —— 元电荷
在室温（\(T = 300 \, \text{K}\)）下，热电压的典型值为：
\[ V_T = \frac{1.380649 \times 10^{-23} \times 300}{1.602176634 \times 10^{-19}} \approx 0.02585 \, \text{V} \approx 25.85 \, \text{mV} \]热电压的温度依赖性对 SPICE 仿真至关重要。下表列出了不同温度下的 \(V_T\) 值：
| 温度 \(T\) (K) | 温度 \(T\) (°C) | 热电压 \(V_T\) (mV) |
|---|---|---|
| 250 | -23.15 | 21.54 |
| 300 | 26.85 | 25.85 |
| 350 | 76.85 | 30.16 |
| 400 | 126.85 | 34.47 |
⚠️ 温度敏感性：在 SPICE 仿真中，热电压随温度线性变化。对于高精度仿真，必须指定仿真温度，否则默认 300 K 下的 \(V_T\) 可能导致显著的电流计算误差，尤其在高温或低温应用场景中。
2.3 反向饱和电流 \(I_S\)
反向饱和电流 \(I_S\) 是肖克利方程中的关键参数，它决定了 PN 结在反向偏置下的漏电流水平，并直接影响正向导通电压。其物理表达式为：
\[ I_S = A \left( \frac{q D_n n_{p0}}{L_n} + \frac{q D_p p_{n0}}{L_p} \right) \]其中：
- \(A\) —— PN 结横截面积
- \(D_n, D_p\) —— 电子和空穴的扩散系数
- \(L_n, L_p\) —— 电子和空穴的扩散长度
- \(n_{p0}, p_{n0}\) —— 平衡少数载流子浓度
由于 \(n_{p0} \propto n_i^2 / N_A\)，\(p_{n0} \propto n_i^2 / N_D\)，且本征载流子浓度 \(n_i^2 \propto T^3 e^{-E_g/kT}\)，因此 \(I_S\) 对温度和禁带宽度高度敏感：
\[ I_S(T) = I_S(T_0) \left( \frac{T}{T_0} \right)^3 e^{\frac{E_g}{k} \left( \frac{1}{T_0} - \frac{1}{T} \right)} \]下表总结了影响 \(I_S\) 的主要因素：
| 影响因素 | 变化趋势 | 对 \(I_S\) 的影响 |
|---|---|---|
| 温度 \(T\) ↑ | 升高 | \(I_S\) 指数级增大（每 10°C 约翻倍） |
| 禁带宽度 \(E_g\) ↑ | 增大 | \(I_S\) 指数级减小 |
| 掺杂浓度 \(N_A, N_D\) ↑ | 升高 | \(I_S\) 减小（少数载流子浓度降低） |
| 结面积 \(A\) ↑ | 增大 | \(I_S\) 线性增大 |
2.4 理想二极管方程结构图
下图展示了肖克利方程中各个物理量之间的关系与信号流向：
2.5 SPICE 中的理想二极管模型
在 SPICE 中，理想二极管模型直接基于肖克利方程实现。其核心参数包括：
| SPICE 参数名 | 符号 | 描述 | 默认值 |
|---|---|---|---|
| IS | \(I_S\) | 反向饱和电流 | 1e-14 A |
| N | \(n\) | 理想因子 | 1.0 |
| EG | \(E_g\) | 禁带宽度 | 1.11 eV (Si) |
| XTI | \(X_{TI}\) | \(I_S\) 温度指数系数 | 3.0 |
SPICE 中二极管直流特性的核心代码实现如下：
** SPICE 理想二极管模型定义
.MODEL Dideal D (
+ IS = 1e-14 ; 反向饱和电流 (A)
+ N = 1.0 ; 理想因子
+ EG = 1.11 ; 禁带宽度 (eV)
+ XTI = 3.0 ; 温度指数系数
)
** 调用示例：D1 anode cathode Dideal
D1 1 0 Dideal
** 直流扫描分析
.DC V1 0 1 0.01
.PROBE
.END
🔬 工程意义：在 SPICE 仿真中，理想二极管模型是复杂二极管模型（如 SPICE Gummel-Poon 模型）的基础。通过调整 IS 和 N 参数，可以拟合实际二极管的直流特性曲线。对于硅二极管，典型 IS 值在 \(10^{-14}\) 到 \(10^{-16}\) A 之间，而肖特基二极管的 IS 通常更大（\(10^{-6}\) 到 \(10^{-8}\) A）。
2.6 理想因子 \(n\) 的物理意义
理想因子 \(n\) 反映了实际二极管与理想肖克利方程的偏离程度：
- \(n = 1\)：纯扩散电流主导，完全符合肖克利理论。
- \(1 < n < 2\)：存在耗尽层复合电流（复合-产生电流）贡献。
- \(n > 2\)：可能存在大注入效应、串联电阻效应或表面漏电流。
在 SPICE 参数提取中，\(n\) 通常通过拟合正向 I-V 特性的中间区域（中电流区）获得。下图展示了不同 \(n\) 值对 I-V 特性的影响：
💡 参数提取技巧：在 SPICE 参数拟合中，通常先在中电流区（0.3V~0.7V 对于硅二极管）通过半对数坐标下的斜率提取 \(n\)，然后通过截距提取 \(I_S\)。这种方法可以有效分离两个参数的影响，提高拟合精度。