2. 工艺仿真基础：氧化、扩散、离子注入、退火等基本工艺步骤的物理模型与 Sentaurus Process 命令入门

集成电路制造的本质是在硅片上通过一系列精确控制的物理化学过程，构建出具有特定电学特性的微纳结构。TCAD 工艺仿真的核心任务，就是在虚拟环境中复现这些过程。本章将深入剖析氧化、扩散、离子注入与退火四大基础工艺步骤的物理模型，并同步介绍其在 Sentaurus Process 中的对应命令与参数设置，为后续的器件结构搭建与校准打下坚实基础。

2.1 氧化工艺：从 Deal-Grove 模型到 Sentaurus 实现

热氧化是硅平面工艺中最基础、最关键的步骤之一，用于生长高质量的 SiO₂ 薄膜，作为栅介质、场氧或注入掩蔽层。其物理机制是氧化剂（O₂ 或 H₂O）通过扩散穿过已生长的氧化层，到达 Si-SiO₂ 界面与硅反应。

2.1.1 Deal-Grove 模型

Deal-Grove 模型是描述热氧化动力学的经典解析模型，它将氧化过程分为三个阶段：

• 气体输运： 氧化剂从气相输运到氧化层表面。
• 固相扩散： 氧化剂穿过已生长的 SiO₂ 层扩散到 Si-SiO₂ 界面。
• 界面反应： 氧化剂在界面与硅发生化学反应，生成新的 SiO₂。

该模型给出的氧化层厚度 \(x_{ox}\) 与时间 \(t\) 的关系为：

\(x_{ox}^2 + A \cdot x_{ox} = B \cdot (t + \tau)\)

其中，\(B\) 为抛物线速率常数（扩散控制），\(B/A\) 为线性速率常数（反应控制），\(\tau\) 为初始氧化层厚度对应的偏移时间。下表总结了干氧与湿氧的关键差异：

2.1.2 Sentaurus Process 氧化命令

在 Sentaurus Process 中，氧化过程通过 Diffuse 命令结合 Temperature 和 Oxidizing 参数实现。以下是一个典型的干氧氧化脚本片段：

# 定义初始衬底
line x location=0.0 spacing=0.01 tag=left
line x location=1.0 spacing=0.01 tag=right
line y location=0.0 spacing=0.002 tag=top
line y location=1.0 spacing=0.1 tag=bottom
region silicon xlo=left xhi=right ylo=top yhi=bottom
init concentration=1e15 phosphorus

# 干氧氧化：900°C，30分钟
diffuse temperature=900 time=0.5 \
        oxidizing ambient=O2 pressure=1.0

# 湿氧氧化：950°C，60分钟
diffuse temperature=950 time=1.0 \
        oxidizing ambient=H2O pressure=1.0

2.2 扩散工艺：Fick 定律与掺杂分布控制

扩散是杂质原子在高温驱动下从高浓度区域向低浓度区域迁移的过程。在 IC 制造中，扩散主要用于形成深结、阱区或作为注入后退火过程中的再分布。

2.2.1 Fick 扩散模型

扩散过程由 Fick 第一定律和第二定律描述：

• Fick 第一定律： 扩散通量 \(J = -D \cdot \nabla C\)，其中 \(D\) 为扩散系数，\(C\) 为杂质浓度。
• Fick 第二定律： \(\frac{\partial C}{\partial t} = \nabla \cdot (D \nabla C)\)，描述浓度随时间的变化。

扩散系数 \(D\) 强烈依赖于温度，遵循 Arrhenius 关系：\(D = D_0 \cdot \exp(-E_a / kT)\)。不同杂质在硅中的扩散行为差异显著：

2.2.2 Sentaurus Process 扩散命令

Sentaurus Process 中的扩散通过 Diffuse 命令执行，可以激活多种物理模型：

# 预淀积扩散：恒定表面浓度
diffuse temperature=950 time=30 \
        phosphorus concentration=1e20

# 推进扩散：恒定总掺杂量
diffuse temperature=1100 time=120

# 激活高级物理模型（浓度依赖扩散）
pdbSet Diffuse Phosphorus cDiffModel 1
diffuse temperature=1000 time=60

2.3 离子注入：精确掺杂的利器

离子注入通过高能离子束将杂质原子直接轰入硅衬底，具有低温工艺、精确剂量控制、各向异性等优势，是现代超大规模集成电路中主要的掺杂手段。

2.3.1 注入物理模型

离子注入的最终掺杂分布由以下因素决定：

• 射程与投影射程： 离子在靶材中的平均穿透深度，由注入能量和离子质量决定。
• 纵向离散 (Straggle)： 由于能量损失随机性导致的分布展宽，用 \(\Delta R_p\) 表示。
• 横向离散： 离子在垂直于入射方向上的散射，影响掩蔽层边缘的横向掺杂分布。
• 沟道效应： 离子沿硅晶格通道方向穿透更深，导致分布尾部拖长。通常通过倾斜注入或预非晶化抑制。

经典的高斯分布模型：\(N(x) = N_{peak} \cdot \exp\left[-\frac{(x - R_p)^2}{2\Delta R_p^2}\right]\)

2.3.2 Sentaurus Process 注入命令

# 硼离子注入：能量 30 keV，剂量 1e13 cm⁻²
implant boron dose=1e13 energy=30 tilt=7 rotation=30

# 磷离子注入：使用双 Pearson 分布模型
implant phosphorus dose=5e14 energy=50 \
        tilt=0 rotation=0 \
        pearson=1

# 砷离子注入：考虑非晶化效应
implant arsenic dose=1e15 energy=80 \
        tilt=7 rotation=30 \
        amorph=1

2.4 退火工艺：激活与修复

离子注入会在硅晶格中造成大量损伤，并使杂质原子处于间隙位置，不具备电活性。退火工艺通过高温热处理实现两个目标：晶格损伤修复和杂质电激活。

2.4.1 退火物理机制

• 固相外延再生长 (SPER)： 在 500-600°C 范围内，非晶层以单晶衬底为籽晶进行外延再生长，修复注入损伤。
• 杂质激活： 杂质原子从间隙位置移动到替位位置，成为电活性掺杂中心。激活率与退火温度和杂质浓度相关。
• 扩散与再分布： 退火过程中杂质同时发生扩散，导致掺杂分布展宽。快速热退火 (RTA) 可在实现高激活率的同时抑制扩散。

2.4.2 Sentaurus Process 退火命令

# 常规炉管退火：修复损伤 + 激活
diffuse temperature=950 time=30

# 快速热退火 (RTA)：高温短时间
diffuse temperature=1050 time=0.1 \
        rampRate=100

# 尖峰退火 (Spike Anneal)：极高温极短时间
diffuse temperature=1100 time=0.01 \
        rampRate=200

# 两步退火：低温预退火 + 高温激活
diffuse temperature=600 time=60
diffuse temperature=1000 time=10

2.5 综合案例：CMOS 阱区工艺仿真脚本

以下是一个完整的 CMOS 阱区工艺仿真脚本，整合了本章所学的所有基础工艺步骤：

# ============================================
# CMOS 双阱工艺仿真脚本 (Sentaurus Process)
# ============================================

# 1. 定义网格与衬底
line x location=0.0 spacing=0.05 tag=left
line x location=2.0 spacing=0.05 tag=right
line y location=0.0 spacing=0.005 tag=top
line y location=5.0 spacing=0.5 tag=bottom
region silicon xlo=left xhi=right ylo=top yhi=bottom
init concentration=1e15 boron

# 2. 场氧化层生长 (湿氧)
diffuse temperature=950 time=60 \
        oxidizing ambient=H2O pressure=1.0

# 3. N-well 注入 (磷)
mask name=nwell_mask left=0.0 right=1.0
implant phosphorus dose=2e13 energy=80 \
        tilt=7 rotation=30 mask=nwell_mask

# 4. P-well 注入 (硼)
mask name=pwell_mask left=1.0 right=2.0
implant boron dose=1e13 energy=50 \
        tilt=7 rotation=30 mask=pwell_mask

# 5. 阱区推进退火
diffuse temperature=1100 time=180

# 6. 去除氧化层
struct tdr=well_formation.tdr

# 输出掺杂分布
select z=Phosphorus
writeSpectral profile_nwell.dat
select z=Boron
writeSpectral profile_pwell.dat

2.6 本章小结