原子层沉积（ALD）气体输送技术因其优异的保形性和亚纳米级厚度控制，已成为先进半导体制造、光学镀膜及新能源器件的核心技术。然而，随着器件特征尺寸微缩至3nm节点以下，高深宽比（HAR，High Aspect Ratio）结构（如3D NAND沟道孔、DRAM电容孔）对前驱体输运提出了严峻挑战。本文从流体力学的基本流态出发，深入分析了粘性流、分子流及过渡流下前驱体分子的输运机制，探讨了输运过程对表面饱和吸附动力学的影响。在此基础上，本文提出了一种基于分区耦合的多尺度输运模型，将宏观气路系统的连续流与微观结构内的稀薄气体动力学相结合，为ALD工艺的均匀性优化提供了理论基础。

1. 引言

ALD的核心在于将气相前驱体交替脉冲送入反应腔室，并利用前驱体分子与基底表面活性位点的自限制化学反应（Self-Limiting Reaction）实现单层生长。理想情况下，只要脉冲时间足够长，整个基底表面（包括深孔内部）均能达到饱和吸附。然而，在实际工艺中，随着深宽比（AR，通常>50:1）的增加，前驱体分子从腔室主体传输到孔底部的过程不再是无阻力的。

这种输运限制导致了一个关键问题：输运受限的非均匀沉积。前驱体分子在进入狭窄结构时，会因与壁面的碰撞或被表面吸附而消耗，导致孔底部的有效剂量远低于孔口，从而引发“夹断效应”（Pinch-off）。要解决这一问题，必须深刻理解不同流态下的输运物理机制，并建立能够跨越从米级（气路管道）到纳米级（反应表面）的多尺度模型。

2. 流态划分及其物理特征

在ALD气体输运中，流态的判别依据是克努森数（Knudsen number, $Kn$），定义为气体分子平均自由程（$\lambda$）与特征尺寸（$d$，如管道直径或孔隙宽度）之比：$Kn=\lambda /d$。

2.1 粘性流（$Kn<0.01$）

在气路总管、质量流量控制器（MFC）以及腔室大空间区域，由于压力较高（通常>100 Pa）且特征尺寸较大，气体处于粘性流（连续流）状态。

• 输运机制：分子间碰撞占主导地位，气体被视为连续介质。流动遵循纳维-斯托克斯方程（Navier-Stokes equations），主要由压力梯度驱动。

• 特点：在粘性流下，前驱体输运速率快，气体分子表现出集体行为，流动阻力相对较小。

2.2 分子流（$Kn>1$）

当气体进入高深宽比结构的内部（如孔径<50 nm）或腔室处于低压（<0.1 Pa）状态时，分子平均自由程大于结构尺寸，进入分子流状态。

• 输运机制：分子-壁面碰撞远多于分子-分子碰撞。此时，输运不再由压力梯度驱动，而是由浓度梯度或分子热运动决定。每个分子独立运动，遵循余弦定律（Knudsen 扩散）。

• 特点：气相传质阻力极大增加，前驱体分子需要经过多次反射才能进入深孔底部。

2.3 过渡流（$0.01\le Kn\le 1$）

在ALD脉冲的中间阶段（如从腔室主体进入通孔入口处），往往处于过渡流状态。该区域分子碰撞与壁面碰撞同等重要，是粘性流与分子流之间的桥梁，也是最难精确建模的区域。

在高深宽比ALD工艺中，分子流和过渡流是限制沉积速率和均匀性的主要物理瓶颈。

 

3. 输运机制对表面饱和吸附动力学的影响

表面饱和吸附动力学通常用Langmuir吸附模型描述，即吸附速率取决于前驱体分压和表面空位浓度。然而，在HAR结构中，宏观的Langmuir模型需要与微观的输运方程耦合。

3.1 反应物消耗引起的梯度效应

在粘性流主导的腔室主体中，前驱体浓度是均匀的。但当气体进入HAR结构内部（分子流或过渡流）时，输运过程与表面反应过程形成竞争关系。

设孔深度为$L$，半径为$R$。在孔内微元 $dz$ 处，前驱体通量 $J(z)$ 的减少量等于侧壁的吸附消耗量：

dJ(z)dz=−2R⋅S⋅F(z)dzdJ(z)​=−R2​⋅S⋅F(z)

其中 $S$ 为粘附系数（Sticking Coefficient），$F(z)$ 为表面通量。

在分子流条件下，通量 $J$ 与浓度梯度遵循Fick第一定律，但扩散系数为Knudsen扩散系数 DkDk​：

Dk=23R8RTπMDk​=32​RπM8RT​​

其中 $R$ 为孔半径，$M$ 为摩尔质量。对于高深宽比结构，有效扩散时间 τ∝L2/Dkτ∝L2/Dk​。这意味着，随着孔深增加，前驱体到达底部的时间呈二次方增长。

3.2 粘附系数的调制作用

粘附系数 $S$ 是连接输运与吸附的关键参数。如果前驱体具有高粘附系数（如 $S>0.1$），分子在第一次碰撞时即被吸附。在分子流状态下，这意味着孔口附近的前驱体分子被迅速消耗，形成一个“遮蔽效应”，导致孔内前驱体通量呈指数衰减。

相反，如果前驱体具有低粘附系数（如 $S<0.01$），分子需要在壁面进行多次反射才能被吸附。这种多次反射增加了分子进入深孔的概率，提高了深宽比容忍度。因此，前驱体分子设计（如使用较大分子量或低粘附系数的前驱体）是改善高深宽比ALD保形性的重要手段。

3.3 瞬态饱和行为

在脉冲时间较短的工业应用中，输运机制决定了饱和特征时间。在粘性流区，饱和时间受限于质量流量控制器的响应速度和腔体体积；在分子流区，饱和时间受限于分子的平均停留时间（Mean Residence Time）。

对于HAR结构，表面达到饱和并非同时发生。通常孔口先饱和，随后饱和前沿向孔底移动。若脉冲时间小于孔底的输运特征时间，则会导致孔底欠饱和沉积。这一动力学过程可以通过扩散-反应方程（Diffusion-Reaction Equation）进行数值求解，揭示出“剂量-深宽比-饱和时间”的三维工艺窗口。

 

4. 多尺度输运模型的建立

由于ALD气体输送设备涉及从宏观（米级）到纳米（埃级）的跨尺度问题，单一尺度的模型无法准确预测工艺结果。我们需要建立一种多尺度耦合模型，将宏观气路、腔室流体力学与微观孔内输运连接起来。

4.1 尺度划分与接口定义

通常将模型划分为三个尺度：

1. 宏观尺度（气路及总管）：尺度为0.1 m ~ 10 m。包括气体钢瓶、阀门、质量流量控制器及传输管道。此区域压力较高（通常1~10 Torr），流态为粘性流（$Kn<0.01$）。建模采用集总参数模型或一维可压缩流，主要关注压力波动的传播速度和脉冲前沿的锐利度。

2. 介观尺度（反应腔室及晶圆表面）：尺度为1 mm ~ 0.5 m。腔室内气体流动可能涉及粘性流或过渡流。此处采用计算流体动力学，但需结合滑移边界条件或直接模拟蒙特卡洛方法。重点在于获取晶圆表面上方（或HAR孔口）的入射通量边界条件。

3. 微观尺度（高深宽比结构内部）：尺度为1 nm ~ 100 μm。流态为分子流或过渡流。此尺度下，连续介质假设失效，需采用Knudsen扩散模型或蒙特卡洛（MC）模拟。

4.2 耦合策略：分区顺序耦合与迭代耦合

建立多尺度模型的核心在于边界条件的传递。

步骤一：宏观输运模拟
首先，建立ALD设备的气路网络模型。输入前驱体蒸汽压、载气流量及阀门时序。通过求解质量守恒和动量守恒方程，计算进入反应腔室的瞬时气体通量 Qin(t)Qin​(t)。输出结果为腔室入口处的压力 Pchamber(t)Pchamber​(t) 和气体组成。

步骤二：腔室流体模拟（介观）
利用CFD软件（如ANSYS Fluent或COMSOL）模拟腔室内的流场。对于低压ALD（如低于1 Torr），需开启稀薄气体模型（如滑移流或DSMC）。该步骤提取晶圆表面（特别是芯片区域）的前驱体通量分布 Jwafer(x,y,t)Jwafer​(x,y,t) 以及孔口处的压力 Pmouth(t)Pmouth​(t)。

步骤三：微观结构内输运模拟（关键）
将孔口压力 PmouthPmouth​ 作为微观模型的入口边界条件。在微观尺度，对于深宽比小于100:1的结构，通常采用一维扩散-反应方程（耦合Knudsen扩散）进行快速工程计算。

然而，对于超高深宽比（>100:1）或复杂几何形状（如弯曲孔道），必须采用颗粒法，如：

• 直接模拟蒙特卡洛（DSMC）：直接求解玻尔兹曼方程，精确模拟分子在结构内部的碰撞与运动轨迹，但计算量巨大。

• 测试粒子蒙特卡洛（TPMC）：忽略分子间碰撞，仅考虑分子-壁面碰撞，适用于纯分子流区域。通过追踪大量粒子的轨迹，统计孔底部的入射通量 JbottomJbottom​ 和累积吸附量。

步骤四：耦合闭环
微观模型计算出的吸附速率或穿透率（Transmission Probability）会反作用于介观模型。例如，如果孔口处消耗了大量前驱体，介观模型中晶圆表面的边界层需要重新校准。因此，通常采用迭代耦合：先假设均匀吸附，计算孔口通量；再根据微观模拟结果修正局部消耗速率，重新运行腔室CFD，直至收敛。

4.3 模型验证与参数校准

多尺度模型必须通过实验校准。通常采用硅通孔或高深宽比测试结构，在ALD沉积后通过扫描电镜（SEM）或透射电镜（TEM）测量膜厚分布。将实测的台阶覆盖率（Step Coverage，定义为孔底厚度与孔口厚度之比）与模型预测对比，反推实际工艺条件下的有效粘附系数。

研究表明，由于表面温度不均匀或表面配体位阻效应，前驱体的有效粘附系数往往低于实验室单晶表面的测量值。因此，模型必须引入一个经验修正因子，以匹配工业实际。

 

5. 结论与展望

本文系统阐述了ALD气体输运过程中粘性流、分子流及过渡流的物理本质，揭示了在高深宽比结构内，输运机制通过控制前驱体分子到达速率与表面吸附反应的竞争关系，直接决定了ALD工艺的饱和动力学特征。

为了精确描述从气路总管到纳米孔底的复杂过程，建立多尺度输运模型是必由之路。该模型通过耦合宏观气路网络、介观腔室流体力学及微观稀薄气体动力学，能够有效预测工艺窗口，并为前驱体化学设计、腔室结构优化以及脉冲时序设定提供理论依据。

随着器件架构向三维堆叠和原子尺度精准制造演进，未来的ALD输运模型将面临更高的精度要求。挑战主要集中在以下两个方面：

1. 瞬态效应：目前多数模型假设稳态或准稳态，而实际ALD脉冲持续时间极短（毫秒级），压力波动的瞬态传播效应不可忽略。

2. 化学反应与输运的强耦合：当前模型往往将化学反应简化为线性粘附系数，而真实情况下，表面配体的覆盖度变化会动态改变粘附系数，形成“输运-反应”的实时反馈回路。

未来的研究需结合原位表征技术（如石英晶体微天平QCM与质谱RGA联用）与机器学习加速的分子动力学，构建数字孪生模型，以实现ALD工艺从“试错法”向“预测性制造”的跨越。

想了解更多原子层沉积 (ALD) 气体输送的基础原理的信息，请登录深圳捷微洛克官网 https://www.jewellok.cn/product-category/ultra-high-purity-regulators了解更多资讯。