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)
突变结掺杂分布与空间电荷区 x 浓度 (log) ND (施主) NA (受主) -xp xn 0 (冶金结) 耗尽层宽度 W = xp + xn +qND -qNA 空间电荷区

4.2 线性缓变结(Linearly Graded Junction)

线性缓变结的掺杂浓度在冶金结附近随位置线性变化。这种分布常见于高温长时间扩散工艺,杂质分布接近余误差函数或高斯分布,在结附近可近似为线性梯度。

线性缓变结的数学描述:

掺杂净浓度:N(x) = ND(x) - NA(x) = a·x

其中 a 为掺杂浓度梯度(cm-4),x=0 处为冶金结位置。

线性缓变结掺杂分布 x 净掺杂 N(x) P型区 (NA > ND) N型区 (ND > NA) -W/2 +W/2 0 (冶金结) 斜率 a = dN/dx

耗尽层宽度公式推导

4.3 突变结耗尽层宽度

基于耗尽层近似,结合泊松方程和边界条件,可推导出突变结的耗尽层宽度公式。

推导步骤:

  1. 泊松方程: d²ψ/dx² = -ρ(x)/εs
  2. 电荷分布:
    • P区 (-xp ≤ x < 0):ρ(x) = -qNA
    • N区 (0 < x ≤ xn):ρ(x) = +qND
  3. 电中性条件: NA·xp = ND·xn
  4. 边界条件: 在 x = -xp 和 x = xn 处,电场 E = 0
  5. 积分得到电场分布,再积分得到电势分布
  6. 利用内建电势 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 突变结电荷分布

在耗尽层近似下,突变结的空间电荷区电荷分布呈现矩形形状:

突变结空间电荷区电荷分布 x ρ(x) +qND -qNA -xp xn 0 ρ=0 ρ=0

电荷分布的关键特征:

  • 矩形分布: 电荷密度在耗尽区内为常数,在边界处突变为零
  • 电中性约束: 总正电荷量 = 总负电荷量,即 qNDxn = qNAxp
  • 不对称扩展: 重掺杂侧耗尽层窄,轻掺杂侧耗尽层宽

4.6 线性缓变结电荷分布

线性缓变结的电荷分布呈三角形:

线性缓变结空间电荷区电荷分布 x ρ(x) -W/2 +W/2 0 -q·a·x +q·a·x ρ=0

线性缓变结电荷分布特征:

  • 三角形分布: 电荷密度随位置线性变化,在冶金结处为零
  • 对称性: 耗尽区关于冶金结(x=0)对称扩展
  • 梯度决定斜率: 电荷密度斜率等于 q·a,a 越大电荷分布越陡峭
耗尽层近似的局限性:
  • 忽略了耗尽区边界处的载流子拖尾效应(Debye 拖尾)
  • 在正向偏压或低掺杂条件下,自由载流子浓度不可忽略
  • 对于超薄耗尽层(纳米尺度),量子效应开始显现
  • TCAD 仿真中通常使用更精确的费米-狄拉克统计和泊松方程自洽求解

4.7 电荷分布与电场、电势的关系

空间电荷区的电荷分布直接决定了电场和电势的分布形态:

物理量 突变结 线性缓变结
电荷分布 ρ(x) 矩形(常数) 三角形(线性)
电场分布 E(x) 三角形(线性) 抛物线(二次)
电势分布 ψ(x) 抛物线(二次) 三次曲线
最大电场位置 冶金结 x=0 处 冶金结 x=0 处
TCAD 仿真验证:

在 Sentaurus TCAD 或 Silvaco 中,通过绘制 DopingConcentration、SpaceCharge、ElectricField 和 Potential 的一维切面图,可以直观验证耗尽层近似下电荷-电场-电势的积分关系。实际仿真结果会显示耗尽区边界处的平滑过渡,与理想近似模型存在细微差异。