2. 理想二极管方程：Shockley方程推导与假设条件

Shockley 方程是描述 PN 结直流电流-电压特性的核心解析模型，也是 SPICE 中二极管模型（D 模型）的数学基础。本节从物理机制出发，逐步推导该方程，并明确其成立所需的理想化假设条件。

2.1 推导起点：PN 结的能带与载流子分布

在热平衡状态下，PN 结的扩散电流与漂移电流相互抵消，净电流为零。当施加正向偏压 \( V_D \) 时，势垒高度降低，非平衡载流子注入中性区，形成扩散电流。推导 Shockley 方程的关键物理量包括：

• 内建电势 \( V_{bi} \)：由掺杂浓度和温度决定，\( V_{bi} = \frac{kT}{q} \ln\left( \frac{N_A N_D}{n_i^2} \right) \)
• 耗尽层宽度 \( W \)：随偏压变化，\( W = \sqrt{ \frac{2\varepsilon_s}{q} \left( \frac{1}{N_A} + \frac{1}{N_D} \right) (V_{bi} - V_D) } \)
• 少子浓度边界条件：在耗尽层边缘，少子浓度与偏压呈指数关系

2.2 少子扩散方程与边界条件

以 N 型中性区（P 区注入电子）为例，稳态少子扩散方程为：

d²Δn_p(x) / dx² = Δn_p(x) / L_n²

其中 \( \Delta n_p(x) = n_p(x) - n_{p0} \) 为过剩电子浓度，\( L_n = \sqrt{D_n \tau_n} \) 为电子扩散长度。通解形式为：

Δn_p(x) = A * exp(-x / L_n) + B * exp(x / L_n)

边界条件由耗尽层边缘（x = 0）和欧姆接触（x → ∞）给出：

代入边界条件，得到 N 区中的过剩电子分布：

Δn_p(x) = n_{p0} [exp(qV_D / kT) - 1] * exp(-x / L_n)

2.3 扩散电流密度与总电流

在耗尽层边缘（x = 0）处，电子扩散电流密度为：

J_n(0) = q * D_n * (dΔn_p / dx)|_{x=0}
       = (q * D_n * n_{p0} / L_n) * [exp(qV_D / kT) - 1]

同理，P 区空穴扩散电流密度为：

J_p(0) = (q * D_p * p_{n0} / L_p) * [exp(qV_D / kT) - 1]

总电流密度为两者之和，乘以结面积 \( A \) 即得 Shockley 方程：

I_D = I_S * [exp(qV_D / (n * kT)) - 1]

其中：
I_S = A * q * [ (D_n * n_{p0} / L_n) + (D_p * p_{n0} / L_p) ]
    = A * q * n_i² * [ (D_n / (L_n * N_A)) + (D_p / (L_p * N_D)) ]

2.4 理想化假设条件

Shockley 方程的推导基于以下 7 条理想化假设，理解这些假设对于后续工程调试至关重要：

1. 突变耗尽层近似：耗尽层边界清晰，中性区无电场，载流子仅通过扩散运动。
2. 小注入条件：注入的少子浓度远小于多子浓度，即 \( \Delta n_p \ll p_{p0} \)，保证中性区电中性。
3. 耗尽层内无复合：通过耗尽层的电流全部由扩散电流贡献，忽略耗尽层内的复合电流（SRH 复合）。
4. 玻尔兹曼统计：载流子服从非简并统计，准费米能级在耗尽层内保持恒定。
5. 均匀掺杂：P 区和 N 区的掺杂浓度均匀，无梯度效应。
6. 一维模型：电流仅沿垂直于结平面的方向流动，忽略边缘效应。
7. 恒温条件：温度固定，不考虑自热效应。

2.5 理想与实际的偏差：SVG 对比图

下图展示了理想 Shockley 方程预测的 I-V 曲线与实际 PN 结测量曲线之间的典型偏差区域：

I_S(T) = I_S(T₀) * (T / T₀)³ * exp[ (E_g / (n * k)) * (1/T₀ - 1/T) ]