3、杂质分布基础:杂质浓度分布函数、高斯分布与Pearson分布模型
在离子注入工艺中,杂质在衬底中的最终分布决定了器件的电学特性。本章系统阐述杂质浓度分布的核心数学描述,重点对比高斯分布与Pearson分布的适用场景与建模方法,为后续注入能量阶梯优化提供理论基础。
3.1 杂质浓度分布函数概述
杂质浓度分布函数 N(x) 描述了注入离子在衬底深度方向上的浓度变化。其基本形式由注入剂量 Φ(单位:cm⁻²)和投影射程 Rp 决定。实际分布受核阻止与电子阻止共同作用,呈现非对称特征。
核心参数定义:
- Φ — 注入剂量(atoms/cm²),控制总杂质数量。
- Rp — 投影射程(nm),离子在衬底中的平均穿透深度。
- ΔRp — 射程离散(straggle),分布宽度的标准差。
- γ — 偏度(skewness),衡量分布的非对称程度。
- β — 峰度(kurtosis),描述分布尾部的拖尾特性。
| 模型名称 | 适用场景 | 参数数量 | 对称性 | 尾部精度 |
|---|---|---|---|---|
| 高斯分布 | 轻离子、低能量、快速估算 | 2(Rp, ΔRp) | 完全对称 | 低 |
| Pearson IV 分布 | 重离子、高能量、精确工艺仿真 | 4(Rp, ΔRp, γ, β) | 非对称 | 高 |
| 双Pearson 分布 | 多峰注入、通道效应修正 | 8(两组参数) | 可描述双峰 | 极高 |
3.2 高斯分布模型
高斯分布是最简单的杂质分布模型,假设离子在衬底中的碰撞过程完全随机,且能量损失均匀。其数学表达式为:
N(x) = (Φ / (√(2π) · ΔRp)) · exp[ - (x - Rp)² / (2 · ΔRp²) ]
高斯分布仅需两个参数:Rp 和 ΔRp。该模型在以下条件下近似有效:
- 注入离子为轻元素(如硼 B⁺、磷 P⁺ 在低能量下)。
- 注入能量较低(< 50 keV),核阻止占主导且散射各向同性。
- 衬底为非晶态或预非晶化处理,通道效应可忽略。
工程限制: 高斯分布无法描述注入分布的非对称拖尾(如硼注入的“通道尾巴”),在深亚微米器件中会导致结深预测偏差达 30% 以上。实际工艺中,高斯模型仅用于快速手算或初步参数扫描。
3.3 Pearson 分布模型
Pearson 分布族是一组由四个矩(均值、方差、偏度、峰度)定义的连续概率分布。在离子注入仿真中,Pearson IV 型是最常用的变体,其概率密度函数满足微分方程:
dN(x)/dx = (x - a) · N(x) / (b₀ + b₁·x + b₂·x²)
其中系数 a、b₀、b₁、b₂ 由 Rp、ΔRp、γ、β 四个参数唯一确定。Pearson IV 的解析解形式复杂,通常通过数值积分或查表法计算。
3.4 高斯分布与 Pearson 分布的选择准则
| 判断条件 | 推荐模型 | 原因 |
|---|---|---|
| |γ| < 0.1 且 β ≈ 3.0 | 高斯分布 | 分布近似对称,高斯误差可接受 |
| 0.1 ≤ |γ| ≤ 0.8 | Pearson IV | 非对称显著,必须修正 |
| |γ| > 0.8 或 β 偏离 3.0 较大 | 双Pearson 或 联合仿真 | 单峰模型失效,需多峰叠加 |
| 注入能量 > 200 keV | Pearson IV | 高能注入尾部拖尾严重 |
工艺经验: 对于硼(B⁺)注入,即使能量低至 30 keV,偏度 γ 也常达到 -0.3 ~ -0.5,此时必须使用 Pearson IV 模型。对于砷(As⁺)或锑(Sb⁺)等重离子,Pearson IV 是标准选择。
3.5 分布参数提取与工程实践
实际工程中,分布参数通常通过以下途径获取:
- SRIM/TRIM 蒙特卡洛仿真:直接模拟离子在衬底中的输运过程,统计输出 Rp、ΔRp、γ、β。
- 实验 SIMS 数据拟合:对二次离子质谱(SIMS)实测浓度曲线进行矩分析,提取四矩参数。
- 工艺仿真器内置数据库:如 TCAD 工具(Sentaurus、ATHENA)预置了常见离子/能量组合的 Pearson 参数表。
参数提取精度要求:
- Rp 和 ΔRp 的误差应控制在 ±5% 以内。
- 偏度 γ 的误差应小于 ±0.1,峰度 β 的误差应小于 ±0.3。
- 对于高能注入(> 200 keV),建议使用双 Pearson 模型或引入通道效应修正项。
3.6 分布模型对注入能量阶梯设计的影响
在注入能量阶梯优化中,分布模型的准确性直接影响阶梯层数的选择和能量间隔的设定:
- 高斯模型:适用于快速估算,但会低估尾部浓度,导致阶梯层数偏少、结深偏浅。
- Pearson IV 模型:准确描述非对称拖尾,可精确计算各能量层的浓度叠加效果,避免过补偿或欠补偿。
- 多峰分布:当注入能量跨度大(如从 10 keV 到 500 keV)时,需考虑不同能量层的分布交叠,使用 Pearson 模型可精确计算总杂质剖面。