薄膜工艺中的数学｜CVD的反应与传输数学

CVD工艺的核心竞争力在于精确预测和控制薄膜的生长行为。Arrhenius方程把温度和速率的关系压缩成一条直线，Fick扩散定律量化了前驱体从气相到表面的传输阻力，一维耗尽模型解释了LPCVD炉管中膜厚沿轴向递减的物理根源，Knudsen扩散-反应模型则预测了高深宽比沟槽中台阶覆盖率的极限。本章从这四个模型出发，给出最小够用的推导与典型数值，便于在配方调整和工艺评估时做定量判断。

一、CVD反应动力学

CVD的沉积速率本质上由化学反应速率决定。理解反应动力学，是工艺工程师预测和调控薄膜生长的基础。

Arrhenius关系与活化能

CVD沉积速率与温度的关系遵循Arrhenius方程。在表面反应控制区域内，沉积速率 G 可以表示为：

G = G₀ · exp[−Eₐ / (R·T)]

其中 G₀ 为指前因子（pre-exponential factor），Eₐ 为表观活化能（apparent activation energy，单位 J/mol），R = 8.314 J/(mol·K) 为气体常数，T 为绝对温度（K）。对上式两边取自然对数：

ln(G) = ln(G₀) − Eₐ / (R·T)

以 ln(G) 对 1/T 作图即为Arrhenius图。在表面反应控制区域内，数据将排列成一条斜率为 −Eₐ/R 的直线。通过线性回归即可从实验数据中提取活化能。在半对数坐标上，横轴通常标为 1000/T （K⁻¹），以便直接读取温度。斜率越陡，说明活化能越高，沉积速率对温度越敏感。

典型CVD工艺的活化能参考值(精简)

从表中可以看到，LPCVD工艺的活化能通常在1.5 eV以上，沉积速率对温度非常敏感。PECVD由于等离子体提供了额外的能量来激活反应，表观活化能显著降低，沉积速率对温度的依赖大幅减弱。这一差异直接影响了两类工艺对温度控制精度的要求：LPCVD通常需要 ±1 °C 以内的温度均匀性，而PECVD对温度波动的容忍度宽得多。

两种速率控制区域

CVD过程中，前驱体分子必须先从气相传输到晶圆表面，然后才能在表面发生化学反应。整个过程由两个串联步骤组成。第一步是气相传输：前驱体分子通过扩散穿越边界层到达表面，传输速率正比于气相浓度与表面浓度的差值。第二步是表面反应：前驱体分子在表面吸附、分解并形成薄膜，反应速率取决于表面浓度和温度。

设气相主体中前驱体浓度为 C_g，晶圆表面处的浓度为 C_s。传输通量 J_T 和反应通量 J_R 分别为：

J_T = h_g · (C_g − C_s)

J_R = k_s · C_s

其中 h_g 为传质系数（mass transfer coefficient），k_s 为表面反应速率常数。在稳态条件下，传输通量等于反应通量（J_T = J_R），由此可以消去表面浓度 C_s，得到总沉积通量：

J = [k_s · h_g / (k_s + h_g)] · C_g

这个表达式具有清晰的物理含义。定义Damköhler数 Da = k_s / h_g，代入可以改写为：

J = h_g · C_g · Da / (1 + Da)

当 Da ≪ 1 时（即 k_s ≪ h_g），表面反应速度远慢于传输速度，整个过程由表面反应控制。此时 J ≈ k_s · C_g，沉积速率对温度高度敏感（通过Arrhenius关系），但对流场和边界层厚度不敏感。这对应LPCVD在反应控制区间的典型工作状态。当 Da ≫ 1 时（即 k_s ≫ h_g），表面反应极快，前驱体一到达表面就立刻被消耗，整个过程由气相传输控制。此时 J ≈ h_g · C_g，沉积速率主要取决于气体流速和扩散系数，对温度的依赖很弱。这对应APCVD在传质控制区间的典型工作状态。

工程判据 

对于LPCVD工艺，通常选择在表面反应控制区域（Da ≪ 1）下操作。原因在于温度可以精确控制（偏差小于 ±1 °C），因此沉积速率的均匀性容易保证。在传质控制区域操作时，均匀性取决于气体流场，控制难度较高。 

图1 CVD沉积速率的Arrhenius图与两种控制区域

二、气相传输与边界层模型

当沉积过程受传质控制，或处于反应与传质的过渡区域时，就必须定量描述前驱体分子从气相主体到晶圆表面的传输过程。本节建立基于Fick扩散定律和边界层理论的传输模型。

Fick扩散定律

Fick第一定律描述了浓度梯度驱动的稳态扩散过程：

J = −D · (dC/dx)

其中 J 为扩散通量（mol/（m²·s）），D 为扩散系数（m²/s），dC/dx 为浓度梯度。负号表示扩散方向从高浓度指向低浓度。

气体的扩散系数与温度和压力密切相关。对于理想气体，扩散系数可以表示为：

D ∝ T^(3/2) / P

即扩散系数正比于温度的3/2次方，反比于压力。这意味着降低压力可以显著提高气体扩散系数。在室温（300 K）、1 atm下，硅烷在氢气中的扩散系数约为 D ≈ 0.2 cm²/s。仅考虑降压效果（温度不变），将压力从760 Torr降至0.5 Torr，扩散系数增大约1520倍。在LPCVD的典型工作条件下（600 °C, 0.5 Torr），同时叠加温度效应后，扩散系数可达 D ≈ 1500 cm²/s，比常压室温条件增大约7500倍。正是这种扩散系数的巨幅提升，使得LPCVD能够实现优良的片内均匀性。

关键结论  

降低总压是CVD工艺从传质控制转向表面反应控制的最有效手段。LPCVD正是利用低压条件下扩散系数的大幅增加，使传质阻力可以忽略，从而让沉积速率完全由温度精确控制。

沉积速率的完整表达

综合反应动力学和传输模型，对于一阶表面反应，沉积速率可以写成：

G = [k_s · D / (δ · k_s + D)] · C_g / N

其中 N 为薄膜材料的原子数密度（atoms/cm³），δ 为边界层厚度，传质系数 h_g ≈ D/δ。这个表达式清楚地体现了两种限制机制的竞争： 传输阻力 1/h_g 和反应阻力 1/k_s 串联相加，总阻力由较大的那项主导，即较慢的步骤成为速率控制步骤。

三、LPCVD立式炉的气体动力学

LPCVD立式炉（或卧式炉管）可以同时处理100片以上的晶圆，是产量极高的批量工艺设备。然而，当前驱体气体沿炉管流动时，会在每一片晶圆表面发生反应而被消耗，导致气相浓度沿流动方向逐渐降低。这种气体耗尽（depletion）效应直接影响片间（wafer-to-wafer）厚度均匀性。

一维耗尽模型

建立如下简化模型： 假设炉管沿轴向为一维流动，晶圆均匀分布在长度为 L 的沉积区域内。设沿轴向坐标 z 处的前驱体浓度为 C（z），入口处（z = 0）的浓度为 C₀。对长度为 dz 的微元进行质量守恒分析。流入的前驱体通量为 Q·C（z），流出的通量为 Q·C（z+dz）。在该微元内，前驱体被晶圆消耗，消耗速率正比于反应速率常数 k_s、浓度 C(z) 和单位长度内的有效反应面积 S（cm²/cm）。对于一阶反应，稳态质量守恒方程为：

Q · dC/dz = −k_s · S · C(z)

这是一个标准的一阶线性常微分方程。以初始条件 C(0) = C₀ 求解，得到浓度沿轴向的分布为指数衰减：

C(z) = C₀ · exp(−k_s·S·z / Q)

由于沉积速率正比于当地浓度（一阶反应），第 j 片晶圆处的沉积速率为：

G(j) = k_s · C₀ · exp(−k_s·S·j·d / Q)

其中 d 为相邻晶圆的间距，j 为从进气口计算的晶圆序号。可以看出，沉积速率从进气端到排气端呈指数递减。该模型假设温度沿炉管恒定。在实际操作中，如果使用温度梯度补偿（即排气端温度略高于进气端），需要同时考虑反应速率常数 k_s 的空间变化。

图2 LPCVD立式炉气体耗尽模型示意图 

片间均匀性与补偿策略

定义片间厚度均匀性为：

U = (G_max − G_min) / (G_max + G_min) × 100%

对于N片晶圆的批次，进气端第1片的沉积速率最高，排气端第N片最低。将指数衰减模型代入： G_max/G_min = exp（k_s·S·N·d/Q）。要控制均匀性，可以采用三种策略。

第一种是增大气体流量Q。增大流量可以减小指数衰减率 k_s·S/Q，但会增加前驱体消耗和废气处理负担。第二种是温度梯度补偿。沿炉管设置逐渐升高的温度分布，使排气端的 k_s 较高，部分补偿浓度下降的效果。第三种是使用注入管。通过沿炉管分布的注入管在多个位置补充前驱体，维持轴向浓度的近似均匀。

四、台阶覆盖率的几何数学

在先进半导体器件中，薄膜需要沉积在包含沟槽、通孔和台阶等复杂三维结构的表面上。薄膜在这些结构上的覆盖均匀程度用台阶覆盖率（step coverage, SC）来衡量。台阶覆盖率直接影响器件的电气连接可靠性和绝缘完整性，是CVD工艺最重要的质量指标之一。

台阶覆盖率的定义

台阶覆盖率定义为沟槽内部薄膜厚度与顶部平面薄膜厚度的比值：

SC_side = t_side / t_top

SC_bottom = t_bottom / t_top

SC = 1（即100%）表示完全保角沉积。对于典型的CVD工艺，SC值介于0.1至1之间，具体取决于沟槽的深宽比（aspect ratio, AR = depth/width）以及前驱体的反应活性。深宽比定义为 AR = H/W，其中 H 为沟槽深度，W 为沟槽宽度。随着半导体技术节点的推进，深宽比不断增大，先进DRAM存储器的电容沟槽AR可以超过50:1，对台阶覆盖率的要求也越来越苛刻。

反应概率与粘附系数

台阶覆盖率的物理本质在于： 前驱体分子进入沟槽后，能否在到达底部或侧壁之前存活（不被消耗）。决定这一行为的关键参数是表面反应概率（surface reaction probability），也称为粘附系数（sticking coefficient），记为 β。β = 1 表示分子每次碰撞表面都会发生反应（一碰即粘），这会导致沟槽入口处薄膜迅速增厚，而底部几乎没有沉积，形成所谓的视线沉积（line-of-sight deposition）。β ≪ 1 表示分子在表面多次反弹后才有可能反应，因此有更大概率穿透到沟槽深处，提高底部和侧壁的覆盖率。

不同CVD工艺的反应概率与台阶覆盖率

Knudsen扩散-反应模型

在窄沟槽中（Knudsen数 Kn = λ/W ≫ 1，其中 λ 为分子平均自由程），气体分子在壁面间的传输服从Knudsen扩散。设沟槽宽度为 W、深度为 H，沿深度方向坐标为 y（y = 0 为开口，y = H 为底部）。前驱体在沟槽内的浓度 C(y) 满足扩散-反应方程：

D_K · d²C/dy² = (2β/W) · v̄ · C(y) / 4

其中 D_K 为Knudsen扩散系数（D_K = W·v̄/3,v̄ 为分子平均热速度）。将 D_K 代入并整理，可以定义一个无量纲Thiele模数：

Φ = H · √(3β / (2W²))

对于窄长沟槽（AR ≫ 1），底部的台阶覆盖率可以近似表示为：

SC_bottom ≈ 1 / cosh(Φ)

当 Φ ≪ 1（反应概率低或沟槽浅）时，SC ≈ 1，接近完全保角。当 Φ ≫ 1 时，SC ≈ 2·exp（−Φ），台阶覆盖率迅速下降。一个广泛使用的工程简化近似是：

SC ≈ 1 / (1 + β·AR / 2)

这个公式虽然是近似的，但很好地抓住了台阶覆盖率随深宽比和反应概率变化的基本趋势： AR 增大或 β 增大都会使覆盖率下降。

设计指导  

要实现 SC > 90% 的保角沉积，需要满足 β·AR < 0.22。对于 AR = 10 的沟槽，要求 β < 0.022。这就解释了为何LPCVD TEOS（β ≈ 0.04）在高深宽比结构中逐渐被ALD替代。 

图3 台阶覆盖率与深宽比/反应概率的关系

台阶覆盖率计算结果

计算结果清楚地表明，低反应概率的LPCVD Poly-Si在极高深宽比下仍能保持近乎完美的保角性。LPCVD TEOS在中等深宽比下表现良好，但在AR = 20时覆盖率降至约71%。PECVD SiO₂由于反应概率较高，在AR = 20时覆盖率仅33%，基本不适用于高深宽比结构。

总结

本文系统建立了CVD工艺中的核心数学模型，涵盖反应动力学、传质模型、反应器设计和台阶覆盖率四个层面。

CVD沉积速率遵循Arrhenius关系，表观活化能是判断工艺机制的关键参数。活化能高于1 eV对应表面反应控制区域，低于0.5 eV通常表示传质控制或等离子体增强效应。传质系数 h_g = D/δ 与扩散系数和边界层厚度相关，Damköhler数 Da = k_s/h_g 是判断速率控制区域的无量纲判据。

LPCVD炉管中的前驱体浓度沿轴向呈指数衰减： C(z) = C₀·exp（−k_s·S·z/Q）。片间均匀性可以通过增大流量、温度梯度补偿或分布式注入来改善。台阶覆盖率由反应概率 β 和深宽比 AR 共同决定，简化公式 SC ≈ 1/(1 + β·AR/2) 提供了快速工程估算。实现高深宽比的保角沉积需要低反应概率的前驱体或采用ALD技术。

CVD工艺的数学模型看上去分散在不同的物理层面，但背后的逻辑是统一的： 所有模型都围绕反应与传输的竞争关系展开。Damköhler数把这种竞争压缩成一个无量纲判据，Arrhenius方程和Fick定律分别量化了竞争的两端，耗尽模型和台阶覆盖率模型则把这种竞争延伸到空间维度上，分别处理轴向的浓度梯度和径向的几何约束。

掌握这些模型的工程师，可以在调整配方参数之前就预判沉积速率、均匀性和覆盖率的变化方向和幅度。这种预判能力，是工艺开发从试错走向设计的分水岭。