第3章雪崩击穿机理

3、雪崩击穿机理：碰撞电离、倍增因子M、电离率积分

各位同学，今天我们来聊聊雪崩击穿的核心机理。说实话，这个知识点我当年学的时候也觉得挺抽象的，直到后来在项目中亲手测到一次雪崩击穿，才真正理解了它的物理本质。

雪崩击穿，说白了就是载流子在强电场下像滚雪球一样倍增的过程。你想想看，一个电子在电场中加速，撞到晶格上产生新的电子-空穴对，新产生的载流子继续加速、继续碰撞……最后电流就像雪崩一样暴涨。嗯，这就是名字的由来。

3.1 碰撞电离：雪崩的种子

碰撞电离是雪崩击穿的起点。一个载流子在电场中获得足够动能后，与晶格原子发生碰撞，把价带电子撞到导带上去，产生一对新的电子-空穴对。

这里有个关键参数——电离率，用 α 和 β 表示。α 是电子的电离率，β 是空穴的电离率。单位是 cm⁻¹，意思是一个载流子漂移 1 厘米能产生多少对电子-空穴对。

电离率的经验公式（Chynoweth 模型）：

α = A_n · exp(-B_n / E)
β = A_p · exp(-B_p / E)

其中 E 是电场强度，A_n、B_n、A_p、B_p 是材料相关的常数。

我在项目中遇到过一个问题：仿真出来的击穿电压总是比实测低 20% 左右。查了半天，发现是电离率参数用的是默认值，没针对我们的工艺修正。后来换了实测提取的参数，结果就对上了。所以啊，电离率参数一定要校准，不然仿真就是自欺欺人。

3.2 倍增因子 M：雪崩的放大倍数

倍增因子 M 定义很简单：

M = I / I₀

I 是发生雪崩后的电流，I₀ 是初始电流（没有雪崩倍增时的电流）。M 越大，说明雪崩效应越强。

当 M → ∞ 时，就发生了击穿。这个临界条件就是：

∫₀^W α · dx = 1

这个积分式就是电离率积分，W 是耗尽层宽度。当电离率积分等于 1 时，一个载流子进入耗尽层，最终会产生无穷多个载流子——击穿发生了。

我个人习惯用这个公式快速估算击穿电压：

M = 1 / (1 - ∫α·dx)

当分母趋近于 0 时，M 趋近于无穷大。这个形式在数值仿真中特别好用。

3.3 电离率积分：击穿判据

电离率积分是判断雪崩击穿的核心判据。对于单边突变结，电场分布是三角形的，电离率积分可以写成：

∫₀^W α(E(x)) · dx = 1

其中电场分布 E(x) 是：

E(x) = E_m · (1 - x/W)

E_m 是最大电场强度，在冶金结界面处。

我曾经犯过一个错误：用均匀电场近似去算电离率积分，结果击穿电压算出来差了 30%。后来才意识到，电场分布的形状对电离率积分影响极大。三角形分布和矩形分布，积分结果能差好几倍。

避坑指南：

电离率对电场极其敏感，E 变化 10%，α 能变化好几倍
温度升高，电离率下降，击穿电压会升高（正温度系数）
不同材料的电离率参数差异很大，Si、GaAs、SiC 的参数不能混用

3.4 知识体系总览

下面这张图是我自己整理的雪崩击穿机理的知识框架，你看一眼就能把整个逻辑串起来：

这张图把雪崩击穿的逻辑链条理得很清楚：从碰撞电离出发，经过倍增因子和电离率积分，最终落到击穿电压上。中间掺杂浓度、温度、材料、结结构这些因素都会影响最终结果。

3.5 实用小工具：电离率积分计算

下面我给出一段 Python 代码，用来计算单边突变结的电离率积分。我自己做仿真时经常用这个做快速验证：

import numpy as np

def ionization_integral(Nd, V, material='Si'):
    """
    计算单边突变结的电离率积分
    
    参数:
        Nd: 掺杂浓度 (cm^-3)
        V: 反向偏压 (V)
        material: 材料类型
    
    返回:
        integral: 电离率积分值
    """
    # 材料参数 (Chynoweth 模型)
    params = {
        'Si':  {'An': 7.03e5, 'Bn': 1.23e6, 'Ap': 1.58e6, 'Bp': 2.04e6},
        'GaAs': {'An': 2.04e5, 'Bn': 5.56e5, 'Ap': 2.04e5, 'Bp': 5.56e5},
    }
    
    p = params[material]
    q = 1.6e-19
    eps = 11.7 * 8.85e-14  # Si 介电常数
    
    # 耗尽层宽度
    W = np.sqrt(2 * eps * V / (q * Nd))
    # 最大电场
    Em = 2 * V / W
    
    # 数值积分
    x = np.linspace(0, W, 1000)
    E = Em * (1 - x / W)
    alpha = p['An'] * np.exp(-p['Bn'] / E)
    beta = p['Ap'] * np.exp(-Bp / E)
    
    integral = np.trapz(alpha + beta, x)
    return integral

# 示例：计算 5V 偏压下的电离率积分
result = ionization_integral(1e16, 5)
print(f"电离率积分 = {result:.4f}")

使用建议：

先用这个脚本快速扫一遍击穿电压范围，再用 TCAD 做精细仿真
电离率积分接近 1 时，步长要加密，不然数值误差会很大
不同材料的参数差异很大，别用 Si 的参数去算 GaAs

好了，雪崩击穿机理就讲到这里。核心就三件事：碰撞电离产生种子载流子，倍增因子 M 衡量放大倍数，电离率积分等于 1 就是击穿判据。下次你在仿真中看到电流突然暴涨，就知道是电离率积分过 1 了。

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