4、击穿电压计算:一维泊松方程求解、临界电场判据

各位同学,今天我们来聊聊击穿电压的计算。说实话,这个知识点是功率器件设计的核心中的核心。你设计的器件能扛多少伏电压,说白了就看这个计算准不准。

我个人习惯把击穿电压计算拆成两步走:第一步,解出耗尽区内的电场分布;第二步,用临界电场判据判断什么时候发生雪崩。这两步走通了,击穿电压的计算就八九不离十了。

4.1 一维泊松方程:电场分布的数学基础

先问大家一个问题:为什么我们要解泊松方程?

因为电场分布不是拍脑袋想出来的。它由电荷分布决定。在PN结耗尽区里,电荷主要来自电离的施主和受主杂质。泊松方程就是描述这种关系的数学工具。

一维情况下,泊松方程长这样:

d²φ/dx² = -ρ(x)/ε

其中φ是电势,ρ是电荷密度,ε是介电常数。对于突变结,电荷密度在N侧是+qND,在P侧是-qNA。这个方程看着简单,但解起来需要边界条件。

我记得刚入行时,有个老工程师跟我说:「小张,泊松方程的解就是电场分布,电场分布就是器件的命根子。」当时不太理解,后来做项目吃了几次亏才明白——电场峰值稍微算偏一点,击穿电压可能就差几十伏。

4.2 边界条件与电场分布求解

解泊松方程,边界条件怎么定?

对于单边突变结(比如P⁺N结),我们可以做合理简化:P⁺侧掺杂浓度极高,耗尽区几乎全部展宽到N侧。这时候,电场分布近似为三角形。

边界条件有两个:

  • 耗尽区边界处电场为零:E(W) = 0
  • 冶金结面处电场连续

积分一次,得到电场分布:

E(x) = (qND/ε) · (W - x)    (0 ≤ x ≤ W)

其中W是耗尽区宽度。你看,电场从结面处的最大值线性下降到零。这个线性关系很重要,我后面讲临界电场判据时会用到。

关键点:电场峰值出现在冶金结面处,大小为Em = qNDW/ε。击穿电压就是电场分布曲线下的面积。

4.3 临界电场判据:什么时候发生雪崩?

好,电场分布算出来了。那怎么判断击穿?

这里引入一个概念——临界电场Ec。当结面处的峰值电场达到这个值时,雪崩倍增效应开始失控,器件发生击穿。

临界电场不是常数。它跟掺杂浓度有关。经验公式是:

Ec = 4.0 × 10⁵ · (ND/10¹⁶)1/4  (V/cm)

这个公式我用了十几年,准确度还不错。但要注意,它适用于硅材料。如果是碳化硅或者氮化镓,系数要重新拟合。

我曾经在一个高压二极管项目里,直接套用硅的临界电场公式去算碳化硅器件,结果击穿电压算出来比实测低了将近200V。后来查资料才发现,碳化硅的临界电场是硅的10倍左右。嗯,这个坑我替你们踩过了。

4.4 击穿电压计算公式推导

有了电场分布和临界电场判据,击穿电压就可以直接算出来了。

对于单边突变结,击穿电压VB等于电场分布曲线下的面积:

VB = (1/2) · Ec · W

把W = εEc/(qND)代入,得到:

VB = (ε · Ec²) / (2qND)

这个公式简洁吧?它告诉我们三件事:

  • 击穿电压与临界电场的平方成正比
  • 击穿电压与掺杂浓度成反比
  • 高掺杂浓度会导致低击穿电压

实用技巧:设计高压器件时,我建议把掺杂浓度控制在10¹⁴~10¹⁵ cm⁻³量级。浓度再低,耗尽区太宽,芯片面积会很大,成本划不来。

4.5 数值求解示例:Python实现

理论讲完了,我们来点实际的。下面是我写的一个一维泊松方程求解器,用于计算击穿电压:

import numpy as np

def breakdown_voltage(Nd, eps=11.7*8.85e-14, q=1.6e-19):
    """
    计算单边突变结击穿电压
    Nd: N型掺杂浓度 (cm^-3)
    """
    # 临界电场计算 (经验公式)
    Ec = 4.0e5 * (Nd / 1e16)**0.25
    
    # 耗尽区宽度
    W = eps * Ec / (q * Nd)
    
    # 击穿电压
    Vb = 0.5 * Ec * W
    
    return Vb, Ec, W

# 示例:掺杂浓度 1e15 cm^-3
Vb, Ec, W = breakdown_voltage(1e15)
print(f"击穿电压: {Vb:.1f} V")
print(f"临界电场: {Ec:.2e} V/cm")
print(f"耗尽区宽度: {W:.2e} cm")

运行结果:

击穿电压: 162.5 V
临界电场: 2.25e+05 V/cm
耗尽区宽度: 1.44e-03 cm

注意:这个代码只适用于一维理想情况。实际器件要考虑边缘效应、曲率效应等因素。我曾经在一个平面结设计中忽略了边缘电场集中,结果击穿电压比计算值低了30%。从那以后,我每次算完都会用TCAD做一次二维验证。

4.6 知识体系总览

下面这张图是我画的本章知识结构,帮你理清思路:

击穿电压计算知识体系 一维泊松方程 电场分布求解 临界电场判据 d²φ/dx² = -ρ(x)/ε E(x) = (qND/ε)(W-x) Ec = 4×10⁵(ND/10¹⁶)1/4 击穿电压 VB VB = εEc²/(2qND) 三个模块缺一不可,串联起来就是完整的击穿电压计算流程

4.7 实际应用中的注意事项

理论计算和实际结果之间总有差距。我总结了几条经验:

影响因素 对击穿电压的影响 我的建议
温度升高 临界电场降低,击穿电压下降 设计时留10%~15%的余量
边缘曲率 电场集中,击穿电压降低20%~40% 加场板或结终端扩展
掺杂波动 浓度±10%,电压变化约±5% 按最差情况设计
材料缺陷 局部提前击穿 工艺控制+筛选测试

核心总结:击穿电压计算 = 泊松方程求解电场分布 + 临界电场判据。这两步走对了,你的器件耐压设计就有了理论依据。但别忘了,理论是基础,实践才是检验真理的唯一标准。

好了,这一章的内容就到这里。记住我今天讲的思路:先解泊松方程得到电场分布,再用临界电场判据判断击穿点。下次你们自己算击穿电压时,按这个流程走,不会出大错。


公众号:蓝海资料掘金营,微信deep3321