PN结耗尽层电场分布与击穿机制：TCAD仿真分析 - 第1节
1.1 PN结基本概念
PN结是半导体器件物理的基石，由P型半导体与N型半导体通过冶金学方法结合而成。在界面处，由于载流子浓度梯度，空穴从P区向N区扩散，电子从N区向P区扩散，形成空间电荷区（耗尽层），并建立内建电场。
核心要点： PN结的I-V特性、电容特性及击穿行为决定了二极管、双极型晶体管（BJT）、金属-氧化物-半导体场效应晶体管（MOSFET）等几乎所有半导体器件的性能边界。
1.1.1 热平衡状态下的能带图
在热平衡下，PN结的费米能级（EF）处处相等。P区价带顶靠近费米能级，N区导带底靠近费米能级，界面处能带弯曲形成势垒，势垒高度为内建电势差 Vbi。
1.2 耗尽层定义与空间电荷区
耗尽层（Depletion Region）又称空间电荷区，是PN结界面附近载流子被耗尽、仅剩固定电离杂质的区域。其宽度 Wdep 随外加偏压变化，是决定结电容和击穿电压的关键参数。
| 参数 | 符号 | 典型值（硅PN结） | 物理意义 |
|---|---|---|---|
| 内建电势 | Vbi | 0.7 ~ 0.9 V | 零偏时势垒高度 |
| 耗尽层宽度 | Wdep | 0.1 ~ 10 μm | 空间电荷区总宽度 |
| 最大电场 | Emax | 105 ~ 106 V/cm | 冶金结界面处峰值电场 |
| 击穿电压 | VBR | 10 ~ 1000 V | 发生雪崩击穿时的反向偏压 |
1.2.1 耗尽层近似
在耗尽层近似下，空间电荷区内的自由载流子浓度近似为零，电荷密度仅由电离施主（ND+）和电离受主（NA-）贡献。泊松方程简化为：
d²ψ/dx² = -ρ(x)/ε_s
其中：
- ψ：电势
- ρ(x)：空间电荷密度
- ε_s：半导体介电常数
对于单边突变结（P⁺N）：
- P区：ρ = -qN_A (x_p < x < 0)
- N区：ρ = +qN_D (0 < x < x_n)1.3 击穿机制概述
当PN结承受反向偏压时，耗尽层内电场增强。当反向电压超过临界值，反向电流急剧增大，发生击穿。主要击穿机制包括：
- 雪崩击穿（Avalanche Breakdown）： 载流子在强电场中获得足够动能，通过碰撞电离产生电子-空穴对，形成倍增效应。这是大多数功率器件的主要击穿机制。
- 齐纳击穿（Zener Breakdown）： 在重掺杂PN结中，耗尽层极窄，强电场直接破坏共价键，产生隧道电流。通常发生在击穿电压 < 5V 的场合。
- 热击穿（Thermal Breakdown）： 焦耳热导致结温升高，本征载流子浓度增加，进一步增大漏电流，形成正反馈。属于二次效应，通常需要结合热仿真分析。
工程提示： 在实际功率器件设计中，雪崩击穿电压通常由耗尽层内电场积分决定。TCAD仿真可以精确提取 Emax 位置和数值，指导场板、场限环等终端结构设计。
1.4 TCAD仿真在器件设计中的重要性
TCAD（Technology Computer-Aided Design）通过数值求解半导体基本方程（泊松方程、载流子连续性方程、输运方程），实现对器件电学特性、物理场分布和工艺影响的精确预测。
1.4.1 TCAD仿真的核心价值
| 应用场景 | 传统方法 | TCAD仿真优势 |
|---|---|---|
| 耗尽层电场分布 | 解析近似（突变结/线性缓变结） | 任意掺杂分布下的精确电场剖面 |
| 击穿电压预测 | 经验公式（误差 > 20%） | 电离积分法，误差 < 5% |
| 温度效应 | 手工修正系数 | 自洽热耦合仿真 |
| 工艺波动分析 | 流片试错（成本高） | DOE虚拟实验，快速优化 |
1.4.2 TCAD仿真流程
- 结构定义： 设置器件几何尺寸、掺杂分布（高斯/阶梯/扩散）、电极位置。
- 物理模型选择： 激活迁移率模型、复合模型（SRH、Auger）、碰撞电离模型（van Overstraeten-de Man、Okuto-Crowell）。
- 网格划分： 在耗尽层、结界面等关键区域加密网格，保证电场梯度分辨率。
- 求解与收敛： 采用Newton或Gummel迭代法，设置偏压扫描步长。
- 后处理分析： 提取电场分布曲线、电离积分、击穿点、电流密度等。
常见误区： 初学者常忽略网格独立性验证。耗尽层内电场梯度极大，若网格过疏，Emax 可能被低估 30% 以上，导致击穿电压预测偏高。建议在结界面处设置最小网格尺寸 < 10 nm。
1.4.3 本课程使用的仿真工具
本课程基于 Sentaurus TCAD（Synopsys）平台，涵盖 SDE（结构编辑器）、SDevice（器件仿真器）、Inspect（可视化工具）等模块。所有案例均提供完整的命令脚本和参数设置，便于学员复现。
; SDE 脚本示例：创建PN结结构
(sdegeo:create-rectangle (position 0.0 0.0 0.0) (position 1.0 1.0 0.0) "Silicon" "region_1")
(sdedr:define-constant-profile "Const.P_5e17" "BoronActiveConcentration" 5e17)
(sdedr:define-constant-profile-region "Const.P_5e17" "region_1")
(sdedr:define-constant-profile "Const.N_1e15" "PhosphorusActiveConcentration" 1e15)
(sdedr:define-constant-profile-region "Const.N_1e15" "region_1")通过本课程的学习，学员将掌握从解析模型到TCAD仿真的完整方法论，能够独立分析PN结耗尽层电场分布，并准确预测雪崩击穿电压。