3. 杂质扩散基础:扩散方程、扩散系数、扩散机制

杂质扩散是半导体工艺中决定结深和杂质分布形态的核心物理过程。在离子注入后的退火步骤中,注入产生的晶格损伤与杂质原子的热运动共同决定了最终的掺杂分布。理解扩散的数学描述、物理参数及微观机制,是优化注入工艺窗口、实现精确杂质分布控制的理论基础。

3.1 扩散方程:菲克定律与工艺建模

杂质扩散行为在宏观上由菲克定律描述。在半导体工艺中,通常考虑一维情况,并假设扩散系数为常数或与浓度相关。

  • 菲克第一定律:稳态扩散下,扩散通量 J 与浓度梯度成正比。
    J = -D · ∂C/∂x
    其中 D 为扩散系数 (cm²/s),C 为杂质浓度 (cm⁻³),x 为深度。
  • 菲克第二定律:非稳态扩散,描述浓度随时间的变化。
    ∂C/∂t = ∂/∂x (D · ∂C/∂x)
    当 D 为常数时,简化为:∂C/∂t = D · ∂²C/∂x²
工艺建模要点: 在 TCAD 仿真中,扩散方程通常与点缺陷(空位、自间隙)的耦合方程联立求解,以准确模拟瞬态增强扩散(TED)效应。

3.2 扩散系数:温度与激活能

扩散系数 D 强烈依赖于温度,遵循阿伦尼乌斯关系:

D = D₀ · exp(-Eₐ / kT)

  • D₀:指前因子(频率因子),与晶格振动频率和跃迁尝试次数相关。
  • Eₐ:扩散激活能 (eV),代表杂质原子跃迁所需克服的能量势垒。
  • k:玻尔兹曼常数 (8.617 × 10⁻⁵ eV/K)。
  • T:温度 (K)。

下表列出了硅中几种常见杂质在无损伤条件下的本征扩散系数参数:

杂质元素 D₀ (cm²/s) Eₐ (eV) 典型温度范围 (°C)
硼 (B) 0.76 3.46 900 – 1200
磷 (P) 3.85 3.66 900 – 1200
砷 (As) 0.066 3.44 1000 – 1200
锑 (Sb) 0.214 3.65 1100 – 1300
工艺提示: 注入损伤会引入大量点缺陷,使表观扩散系数比本征值高出数个数量级。TCAD 仿真中必须使用与点缺陷耦合的扩散模型,而非简单的本征扩散系数。

3.3 扩散机制:间隙、空位与 Kick-out

杂质在硅晶格中的扩散主要通过点缺陷(空位和自间隙原子)介导。根据杂质与点缺陷的相互作用方式,可分为以下几种机制:

  1. 空位机制 (Vacancy Mechanism):杂质原子占据晶格格点,当邻近出现空位时,杂质跃迁进入空位。适用于磷、砷、锑等替代式杂质。
  2. 间隙机制 (Interstitial Mechanism):杂质原子在晶格间隙位置之间直接迁移。适用于锂、钠等小尺寸杂质。
  3. 间隙-替代交换机制 (Interstitial-Substitutional Exchange):杂质在间隙位和替代位之间转换,包括:
    • Kick-out 机制:间隙杂质(如硼间隙)将替代位硅原子踢出,自身占据格点,产生自间隙原子。
    • Frank-Turnbull 机制:间隙杂质与空位结合,进入替代位。
关键概念: 硼在硅中主要通过 Kick-out 机制扩散,其扩散速率受自间隙原子浓度控制。注入损伤产生的过量自间隙原子是导致硼瞬态增强扩散(TED)的根本原因。

3.4 扩散机制与点缺陷耦合关系图

下图展示了硅中杂质扩散的主要机制及其与点缺陷的耦合关系:

硅中杂质扩散机制与点缺陷耦合关系 硅中杂质扩散机制与点缺陷耦合关系 硅晶格格点 空位 (V) 自间隙 (I) 替代杂质 (Aₛ) 间隙杂质 (Aᵢ) 空位扩散 间隙扩散 Kick-out: Aᵢ + Siₛ → Aₛ + I Frank-Turnbull: Aᵢ + V → Aₛ 注入损伤产生大量自间隙 (I) 和空位 (V), 显著增强杂质扩散,尤其是通过 Kick-out 机制的硼扩散。

3.5 扩散机制对比与工艺影响

扩散机制 主导杂质 依赖的点缺陷 对注入损伤的敏感性 工艺影响
空位机制 P, As, Sb 空位 (V) 中等 损伤退火后空位浓度升高,加速扩散
间隙机制 Li, Na (少量) 自间隙 (I) 低 在硅工艺中不常见
Kick-out 机制 B, In 自间隙 (I) 高 注入损伤导致严重的 TED,结深难以控制
Frank-Turnbull 机制 Au, Pt (深能级) 空位 (V) 中等 用于寿命控制,受空位浓度影响
联合优化思路: 在注入工艺窗口拓展中,需根据杂质种类选择退火条件。对于硼注入,采用快速热退火 (RTA) 或激光退火,缩短高温时间,抑制 TED;对于磷/砷注入,可通过预非晶化注入 (PAI) 控制空位分布,实现陡峭结。

3.6 扩散方程在 TCAD 中的实现

在 TCAD 仿真工具(如 Sentaurus Process)中,杂质扩散的求解基于以下耦合方程组:

∂C_A/∂t = ∇·(D_A ∇C_A) + 生成项 - 复合项
∂C_I/∂t = ∇·(D_I ∇C_I) + G_I - R_I
∂C_V/∂t = ∇·(D_V ∇C_V) + G_V - R_V
  • C_A:杂质浓度
  • C_I, C_V:自间隙和空位浓度
  • D_A, D_I, D_V:对应的扩散系数
  • G, R:点缺陷的生成率和复合率,与注入损伤相关
仿真实践: 设置硼注入后的退火工艺时,必须激活“点缺陷耦合扩散”模型,并输入注入损伤分布(+1 模型或 Monte Carlo 损伤分布),否则仿真结果将严重偏离实际。