PN结耗尽层电场分布与击穿机制:TCAD仿真分析 - 第4节
4. 耗尽层近似:突变结与线性缓变结
耗尽层近似(Depletion Approximation)是分析PN结空间电荷区最核心的简化模型。其基本假设为:在耗尽区内,自由载流子(电子和空穴)完全耗尽,仅留下固定的电离施主和受主电荷;在耗尽区外,半导体保持电中性,电场为零。根据掺杂分布的不同,耗尽层近似可分为突变结和线性缓变结两种基本类型。
4.1 突变结(Abrupt Junction)
突变结是指P区和N区的掺杂浓度在冶金结界面处发生阶跃式变化。实际工艺中,低能量离子注入或浅扩散形成的PN结可近似为突变结。
- 掺杂分布:ND(x) 和 NA(x) 在 x=0 处突变
- 空间电荷区电荷密度:ρ(x) = q(ND - NA)
- 耗尽层宽度与偏压呈平方根关系:W ∝ √(Vbi + VR)
4.2 线性缓变结(Linearly Graded Junction)
线性缓变结的掺杂浓度在冶金结附近随位置线性变化。这种分布常见于高温长时间扩散工艺,杂质分布接近余误差函数或高斯分布,在结附近可近似为线性梯度。
掺杂净浓度:N(x) = ND(x) - NA(x) = a·x
其中 a 为掺杂浓度梯度(cm-4),x=0 处为冶金结位置。
耗尽层宽度公式推导
4.3 突变结耗尽层宽度
基于耗尽层近似,结合泊松方程和边界条件,可推导出突变结的耗尽层宽度公式。
推导步骤:
- 泊松方程: d²ψ/dx² = -ρ(x)/εs
- 电荷分布:
- P区 (-xp ≤ x < 0):ρ(x) = -qNA
- N区 (0 < x ≤ xn):ρ(x) = +qND
- 电中性条件: NA·xp = ND·xn
- 边界条件: 在 x = -xp 和 x = xn 处,电场 E = 0
- 积分得到电场分布,再积分得到电势分布
- 利用内建电势 Vbi 求解耗尽层宽度
W = xp + xn = √[ (2εs/q) · (1/NA + 1/ND) · (Vbi + VR) ]
其中:
- εs:半导体介电常数(硅为 11.7ε0)
- q:电子电荷量(1.602×10-19 C)
- Vbi:内建电势
- VR:外加反向偏压
单边突变结(One-sided Abrupt Junction): 当一侧掺杂浓度远高于另一侧时(如 N+P 结),耗尽层主要向轻掺杂侧扩展,公式可简化为:
W ≈ √[ (2εs·(Vbi + VR)) / (q·NB) ]
其中 NB 为轻掺杂侧的掺杂浓度(衬底浓度)。
4.4 线性缓变结耗尽层宽度
对于线性缓变结,净掺杂分布为 N(x) = a·x,泊松方程变为:
d²ψ/dx² = -q·a·x / εs
通过两次积分并利用边界条件,得到耗尽层宽度:
W = [ (12εs·(Vbi + VR)) / (q·a) ]1/3
| 参数 | 突变结 | 线性缓变结 |
|---|---|---|
| 耗尽层宽度 W | ∝ (Vbi + VR)1/2 | ∝ (Vbi + VR)1/3 |
| 最大电场 Emax | ∝ (Vbi + VR)1/2 | ∝ (Vbi + VR)2/3 |
| 结电容 Cj | ∝ (Vbi + VR)-1/2 | ∝ (Vbi + VR)-1/3 |
| 典型应用 | 浅结、低能量注入 | 深结、高温扩散 |
空间电荷区电荷分布
4.5 突变结电荷分布
在耗尽层近似下,突变结的空间电荷区电荷分布呈现矩形形状:
电荷分布的关键特征:
- 矩形分布: 电荷密度在耗尽区内为常数,在边界处突变为零
- 电中性约束: 总正电荷量 = 总负电荷量,即 qNDxn = qNAxp
- 不对称扩展: 重掺杂侧耗尽层窄,轻掺杂侧耗尽层宽
4.6 线性缓变结电荷分布
线性缓变结的电荷分布呈三角形:
线性缓变结电荷分布特征:
- 三角形分布: 电荷密度随位置线性变化,在冶金结处为零
- 对称性: 耗尽区关于冶金结(x=0)对称扩展
- 梯度决定斜率: 电荷密度斜率等于 q·a,a 越大电荷分布越陡峭
- 忽略了耗尽区边界处的载流子拖尾效应(Debye 拖尾)
- 在正向偏压或低掺杂条件下,自由载流子浓度不可忽略
- 对于超薄耗尽层(纳米尺度),量子效应开始显现
- TCAD 仿真中通常使用更精确的费米-狄拉克统计和泊松方程自洽求解
4.7 电荷分布与电场、电势的关系
空间电荷区的电荷分布直接决定了电场和电势的分布形态:
| 物理量 | 突变结 | 线性缓变结 |
|---|---|---|
| 电荷分布 ρ(x) | 矩形(常数) | 三角形(线性) |
| 电场分布 E(x) | 三角形(线性) | 抛物线(二次) |
| 电势分布 ψ(x) | 抛物线(二次) | 三次曲线 |
| 最大电场位置 | 冶金结 x=0 处 | 冶金结 x=0 处 |
在 Sentaurus TCAD 或 Silvaco 中,通过绘制 DopingConcentration、SpaceCharge、ElectricField 和 Potential 的一维切面图,可以直观验证耗尽层近似下电荷-电场-电势的积分关系。实际仿真结果会显示耗尽区边界处的平滑过渡,与理想近似模型存在细微差异。