3. 杂质分布基础：高斯分布模型及其在离子注入中的适用性与局限性
离子注入后，杂质在衬底中的深度分布是决定器件电学特性的核心因素。在注入能量与剂量固定的前提下，杂质浓度随深度的变化通常呈现为一种对称的“钟形”曲线。在早期工艺建模中，高斯分布模型因其数学形式简洁、物理意义明确，被广泛用作描述注入杂质分布的一阶近似工具。
核心概念：高斯分布模型假设注入离子在靶材中经历的核阻止与电子阻止过程具有统计随机性，且能量损失事件相互独立，从而使得杂质在投影射程（Rp）附近呈正态分布。
3.1 高斯分布模型的数学表达
对于垂直入射的单能离子束，杂质浓度 N(x) 随深度 x 的变化由下式给出：
N(x) = (Φ / (√(2π) · ΔRp)) · exp[ - (x - Rp)² / (2 · ΔRp²) ]其中：
- Φ —— 注入剂量（atoms/cm²），决定曲线下的总面积。
- Rp —— 投影射程（Projected Range），即杂质浓度达到峰值时的平均深度。
- ΔRp —— 投影射程的标准偏差（Straggle），表征分布展宽程度。
该模型的关键参数 Rp 和 ΔRp 通常通过蒙特卡洛模拟（如 SRIM/TRIM）或实验测量获得，并随注入能量与离子-靶材组合而变化。
工艺提示：在低能量（< 10 keV）注入轻离子（如 B⁺）时，高斯分布仍能给出可接受的近似；但在高能量或重离子注入时，分布会明显偏离对称高斯形状。
3.2 高斯分布模型的结构化图示
下图展示了高斯分布模型中杂质浓度与深度的关系，并标注了关键参数 Rp、ΔRp 以及峰值浓度 Nmax 的位置。
3.3 高斯分布模型的适用条件
高斯分布模型在以下工艺场景中具有较好的适用性：
| 适用场景 | 具体条件 | 典型示例 |
|---|---|---|
| 低能量轻离子注入 | 能量 < 30 keV，离子质量 < 31 amu | B⁺ → Si，能量 10–20 keV |
| 非晶靶材 | 靶材无长程有序结构，沟道效应可忽略 | 预非晶化后的硅衬底 |
| 低剂量注入 | 剂量 < 10¹⁴ cm⁻²，浓度远低于固溶度 | 阈值电压调整注入 |
| 快速工艺估算 | 无需精确尾部分布，仅需峰值与结深 | 早期工艺设计或教学演示 |
工程经验：在 0.18 μm 及以上技术节点，高斯分布常用于快速估算源/漏延伸区（SDE）的结深；但在先进节点（≤ 28 nm）中，必须使用更精确的分布模型。
3.4 高斯分布模型的局限性
实际离子注入过程远比高斯分布描述的复杂。以下列出该模型的主要局限性：
- 分布不对称性（Skewness）：高能或重离子注入时，杂质分布呈现明显的“拖尾”现象——前缘陡峭、后缘平缓，高斯分布无法捕捉这种偏态。
- 沟道效应（Channeling）：在单晶靶材中，部分离子沿晶格通道深入穿透，形成深尾分布，高斯模型完全忽略这一机制。
- 浓度依赖的扩散效应：高温退火过程中，高浓度梯度会导致扩散系数变化，高斯分布无法预测退火后的再分布。
- 多能量/多角度注入：实际工艺常采用多次注入以形成平坦分布，高斯模型仅适用于单能单次注入。
3.5 高斯分布与真实分布的对比
下表总结了高斯分布模型与真实注入分布之间的关键差异：
| 特征 | 高斯分布模型 | 真实注入分布 |
|---|---|---|
| 分布形状 | 完全对称 | 通常呈正偏态（右拖尾） |
| 峰值位置 | 严格位于 Rp | 可能因沟道效应略浅 |
| 尾部行为 | 指数衰减，无长尾 | 存在指数型或幂律型长尾 |
| 对能量的依赖性 | 仅通过 Rp 和 ΔRp 体现 | 高阶矩（偏度、峰度）也随能量变化 |
| 退火后预测能力 | 无法直接使用 | 需耦合扩散模型 |
3.6 超越高斯：高阶分布模型简介
为了克服高斯分布的局限性，工艺仿真中引入了更精确的分布模型：
- Pearson IV 分布：通过引入偏度（Skewness）和峰度（Kurtosis）两个高阶矩，能够准确描述非对称分布和拖尾特征，是目前 TCAD 仿真中的标准模型。
- 双 Pearson 分布：将注入分布拆分为“主峰”和“沟道尾”两部分，分别用两个 Pearson 函数拟合，适用于单晶衬底中的沟道效应校正。
- 蒙特卡洛（Monte Carlo）模拟：直接模拟每个离子的碰撞历史，精度最高，但计算成本大，常用于模型校准和特殊工艺开发。
注意：高斯分布虽然简单，但在先进工艺（如 7 nm 及以下节点）中可能导致结深预测偏差达 30%–50%，从而影响器件阈值电压和短沟道效应的评估。务必根据工艺窗口选择合适的分布模型。
3.7 高斯分布与高阶模型的结构对比
下图对比了高斯分布与真实分布（以 Pearson IV 为例）在形状上的差异，突出显示非对称拖尾效应。
3.8 小结
高斯分布模型为理解离子注入杂质分布提供了直观的数学框架，在低能量、轻离子、非晶靶材等受限条件下仍具实用价值。然而，随着工艺节点不断微缩，沟道效应、高能拖尾和多能量注入等复杂因素使得高斯模型的局限性日益突出。现代 TCAD 仿真已普遍采用 Pearson IV 或蒙特卡洛方法，以更准确地预测杂质分布，从而支撑精确的器件设计与工艺优化。