粗糙化Kapton H表面上氩和氧原子的非弹性散射动力学：实验与建模

《The Journal of Physical Chemistry C》：Inelastic Scattering Dynamics of Argon and Oxygen Atoms on Roughened Kapton H Surfaces: Experiments and Modeling

　　

表面粗糙度与气体-表面散射动力学之间的关系对多种应用至关重要。然而，关于在良好表征的粗糙表面上高能（hyperthermal）气体-表面散射动力学的详细实验数据仍然缺乏。在本工作中，研究人员利用实验和计算工具研究了表面粗糙度与气体-表面散射动力学之间的联系。研究人员进行了分子束-表面散射实验，以评估高能氩（7.1和5.3 eV）和氧（5.4 eV）原子束在六种不同聚酰亚胺聚合物（Kapton H）表面上的散射动力学，这些表面的均方根粗糙度（root-mean-squared roughness, S_q）跨越两个数量级。低粗糙度表面产生准镜面角通量分布，主要由脉冲散射（impulsive scattering, IS）主导，而粗糙表面则表现出宽角分布、增强的热解吸（thermal desorption, TD）以及更大的面外散射。能量转移取决于入射原子及其能量，但在固定偏转角下基本与粗糙度无关。能量适应系数（energy accommodation coefficient, α_E）和热解吸概率均随入射角减小和表面粗糙度增加而增加，在最粗糙表面接近1。为了模拟这些结果，研究人员使用γ分布表面扩展了高斯（Gaussian）洗板粗糙度模型（washboard roughness model）。γ模型相比高斯描述改进了预测，准确再现了角通量和能量转移的实验趋势，以及从实验数据估计的适应系数。拟合的粗糙度参数与原子力显微镜（atomic force microscopy, AFM）测量的形貌参数呈正相关，尽管并未直接匹配这些参数。当结合气体特定的局部气体-表面散射模型时，一组单一的粗糙度参数能够捕捉不同能量和入射原子下的散射动力学。这些结果确立了该模型在稀薄气体条件下模拟气体-表面散射动力学的实际价值，例如在极低地球轨道（very low Earth orbit, VLEO）中卫星与残余大气的相互作用。

**论文解读：粗糙化Kapton H表面上氩和氧原子的非弹性散射动力学：实验与建模**

**研究背景与问题**

表面粗糙度与气体-表面散射动力学之间的关系在多个领域具有重要意义，包括卫星轨道维持、大气再入、材料加工以及微机电系统等。在稀薄气体流动中，气体粒子与表面之间的质量、动量和能量传递由散射动力学决定，直接影响卫星阻力等关键参数。然而，现有研究缺乏针对良好表征粗糙表面的高能（hyperthermal）气体-表面散射动力学的详细实验数据。尽管已有多种现象学模型（如Maxwell模型、Cercignani-Lampis (CL) 模型）和基于统计描述的洗板模型（washboard model）被提出，但先前的研究表明，这些模型在描述高能气体在粗糙表面上的散射动力学时存在明显不足。例如，Liang等人提出的三维互易性满足洗板模型（3D-RW）在模拟各向同性高斯（Gaussian）随机表面时，未能准确再现实验观测到的角通量分布，尤其是在小最终极角（θ_f）和向后半球（backward hemisphere）区域的散射通量。此外，几何粗糙度与散射动力学之间的直接联系仍缺乏深入理解，尤其是当基底材料与表面粗糙度同时变化时，难以分离各自的影响。因此，本研究旨在通过实验和计算手段，系统评估表面粗糙度如何影响气体-表面散射动力学，并发展一种改进的模型来准确描述这一过程。

**研究内容与结论**

研究人员对六种具有不同粗糙度的Kapton H聚酰亚胺聚合物表面进行了分子束-表面散射实验。使用高能氩（Ar）原子束（7.1和5.3 eV）和氧（O）原子束（5.4 eV）作为入射粒子，测量了不同入射角下的角通量分布、能量转移以及热解吸概率。同时，利用原子力显微镜（AFM）表征了表面形貌，获取了均方根粗糙度（S_q）和极角分布。基于实验数据，研究人员扩展了传统高斯洗板模型，引入了γ分布表面模型（3D-RW-Γ），以描述高度非高斯分布的表面形貌。通过对比实验与模型结果，发现γ模型在预测角通量分布、能量转移以及能量适应系数（α_E）和热解吸概率（f_TD）方面优于高斯模型。模型拟合的粗糙度参数与AFM测量的形貌参数呈正相关，但并非直接对应，表明模型中的粗糙度不仅反映了几何形貌，还包含了气体-表面相互作用势的波纹度。此外，同一组模型参数能够结合气体特定的局部散射模型，再现不同能量和入射原子（Ar和O）的散射动力学。该研究发表在《The Journal of Physical Chemistry C》。

**技术方法**

本研究主要采用以下关键技术方法：（1）分子束-表面散射实验：使用激光驱动击穿源产生高能脉冲原子束，通过旋转样品架和可旋转质谱仪检测散射产物的角通量和飞行时间（TOF）分布，获取非弹性散射（IS）和热解吸（TD）组分。（2）原子力显微镜（AFM）表征：在轻敲模式下获取5×5 μm²的形貌图，计算均方根粗糙度（S_q）和极角分布。（3）混合模型（hybrid model）：结合CL+扩散通量模型和软球碰撞模型，从实验数据中估计半球上的能量适应系数和热解吸概率。（4）三维互易性满足洗板模型（3D-RW-Γ）：采用γ分布随机表面描述各向同性粗糙度，耦合CL局部散射模型，通过马尔可夫链蒙特卡洛（MCMC）采样或手动优化拟合实验数据。样本来源为杜邦公司125 μm厚的Kapton H薄膜，通过台面激波管（TTST）中不同距离的原子氧暴露获得不同粗糙度表面。

**研究结果**

**4.1 表面形貌（Surface Topography）**
通过AFM形貌图，六种表面可分为两组：低粗糙度组（S_q = 1.3–5.6 nm，极角<30°）和高粗糙度组（S_q = 16–162.6 nm，极角>45°）。γ分布模型能很好地拟合AFM极角分布，多数表面接近高斯统计（N>16），仅样品V需要较低N值。

**4.2 散射实验（Scattering Experiments）**
**4.2.1 7.1 eV Ar原子的非弹性散射动力学**：低粗糙度表面（如原始表面、样品II、III）的脉冲散射（IS）通量呈准镜面分布，峰值在θ_f=65–70°；高粗糙度表面（样品IV–VI）的IS通量显著降低，热解吸（TD）通量增强，且出现面外散射增多。能量转移随偏转角（χ）变化，与表面粗糙度无关，但受入射原子和能量影响。采用能量分辨方法（ERM）处理的数据显示，相同偏转角下不同表面的平均能量转移曲线重合，可用软球模型拟合（有效表面质量m_s = 92.3 amu，内部能量激发比?E_int?/?E_i? = 0.29）。

**4.2.2 入射能量和气体种类对散射动力学的影响**：在样品II上，5.3 eV Ar和7.1 eV Ar的IS角通量分布相似，但5.4 eV O的分布更宽，TD通量更低。能量转移曲线显示，5.4 eV O的能量转移最小（m_s = 48.1 amu，?E_int?/?E_i? = 0.20），而7.1 eV Ar最大（m_s = 89.2 amu，?E_int?/?E_i? = 0.28）。较低的动量与较短的相互作用时间解释了O原子更低的能量转移。

**4.3 积分量的估计（Estimation of Integral Quantities）**
利用混合模型，从实验数据中估计了能量适应系数（α_E）和热解吸概率（f_TD）。对于7.1 eV Ar，α_E和f_TD随粗糙度增加而显著升高：原始表面在θ_i=60°时α_E=0.81，f_TD=0.20；样品VI时α_E≈0.99，f_TD≈0.91。不同入射角和气体种类下，α_E和f_TD随入射角减小而线性增加，O的f_TD显著低于Ar。

**4.4 3D-RW-Γ模型的适用性（Suitability of the 3D-RW-Γ Model）**
直接使用AFM形貌参数输入3D-RW-Γ模型时，预测的角通量分布远窄于实验观测；允许粗糙度参数自由优化后，模型能较好地再现角通量和能量转移趋势，且γ模型优于高斯模型。模型粗糙度参数A与AFM形貌参数呈正相关，但数值更大（例如原始表面A=0.7 vs AFM的A=0.07），表明模型包含了势能波纹度。对于不同入射能量和气体种类，同一表面模型可通过调整局部散射模型参数（α_n, α_t）拟合数据，但α_n, α_t需随能量变化，表明CL模型在宽能量范围存在局限性。模型预测的平动能分布比实验更窄，但平均值一致。

**总结讨论与结论翻译**

**总结讨论**：尽管3D-RW-Γ模型存在一些局限性，例如当前仅能描述各向同性表面、需要能量依赖的局部适应系数来覆盖宽能量范围、以及预测的平动能分布较窄，但该模型成功地再现了角通量和能量转移分布的主要趋势，以及宏观能量适应系数和热解吸概率。模型能够使用同一组表面参数结合气体特定的局部散射模型，描述不同气体种类和入射能量下的散射动力学，这对于稀薄高能流动模拟（如卫星在极低地球轨道（VLEO）中的气动阻力）具有重要意义。模型参数与AFM形貌参数呈正相关，但并非直接对应，表明模型中的粗糙度是有效粗糙度，综合了几何形貌和势能波纹度。未来工作将扩展γ模型以描述各向异性表面，并改进局部碰撞模型（如采用速度相关波纹度模型）以解决上述局限性。

**研究结论**（翻译自论文结论部分）：
研究人员在六种具有不同粗糙度的Kapton H样品上进行了分子束-表面散射实验，以考察气体-表面散射如何依赖于表面形貌。原子力显微镜用于表征样品形貌，揭示了两个不同的表面组：低粗糙度样品的S_q值在1.3至5.6 nm之间，局部表面法线极角低于30°；粗糙样品的S_q值在16至160 nm之间，局部表面法线极角超过45°。使用不同高能Ar和O束测量了散射动力学。低粗糙度表面产生以脉冲散射为主的叶状角通量分布，而粗糙表面产生宽分布，脉冲散射极小，热解吸增强。对于最粗糙表面，仅观察到热解吸。面外测量表明，更粗糙的表面和更高的入射能量增加了面外散射通量的比例。氧原子产生的角通量分布比氩原子更宽，但热解吸贡献更低。脉冲散射粒子的平动能分布覆盖从热能量到入射束能量的整个范围；这些分布随气体种类和能量变化，但在固定偏转角下基本不受表面粗糙度影响。分数能量转移曲线作为偏转角的函数塌缩为单一曲线，并由软球模型很好地描述，符合结构区散射。能量转移取决于入射能量和气体种类，5.3 eV Ar的能量转移低于7.1 eV Ar，5.4 eV O原子的能量转移最低。研究人员使用混合模型估计宏观气体-表面参数。能量适应系数和热解吸概率随入射角减小而近似线性增加。此外，随着表面粗糙度增加，观察到更高的能量适应系数和更高的热解吸概率。对于光滑表面和60°入射角，最低测量能量适应系数为5.4 eV O的0.68和7.1 eV Ar的0.79，而粗糙表面促进了接近1的能量适应系数和热解吸概率。研究人员开发了一种在洗板框架下的表面粗糙度模型（3D-RW-Γ），将表面表示为各向同性γ随机表面。从AFM数据导出的表面极角分布由γ粗糙度模型很好地描述，大多数表面呈现高斯统计，一个表面呈现γ统计。使用与AFM测量统计相似的模型表面对束-表面散射数据进行3D-RW-Γ模型拟合，产生的角通量分布比实验观测狭窄得多。允许粗糙度参数变化后，对角通量和能量转移趋势达到了相当好的一致性，表明表面模型不仅描述了几何粗糙度，还描述了气体-表面势的有效波纹度。γ表面模型相比高斯模型改进了预测，提供了相当准确的能量适应和热解吸概率估计。尽管模型的粗糙度参数未直接匹配AFM测量，但测量与模型表面粗糙度之间存在清晰的相关性，更粗糙表面的散射由具有更高波纹度参数的模型表面拟合。此外，研究人员的结果表明，同一表面模型能够捕捉不同入射能量和不同气体种类下的实验散射动力学，入射气体性质的调整由局部碰撞模型处理。3D-RW-Γ模型仍存在一些局限性：首先，γ表面模型目前仅能描述各向同性表面，而工程兴趣的表面可能具有各向异性形貌，当前模型无法捕捉；其次，准确再现实验数据需要为5.3 eV Ar使用比7.1 eV Ar更低的局部碰撞模型能量适应系数，表明单一系数集无法覆盖宽能量范围；最后，预测的平动能分布比测量值窄，尽管其平均值相似。这些局限性均指向局部碰撞模型是3D-RW框架中的主要限制因素，改进该组件将是未来工作的重点。尽管存在这些局限性，采用γ表面的3D-RW模型成功再现了角通量和能量转移分布的主要趋势，以及宏观能量适应和热解吸概率。跨不同气体种类使用单一表面模型的能力对于稀薄高能流动模拟是有利的。在这些流动中，每种气体的入射能量范围有限，因此使用恒定适应系数可能是一种准确描述散射动力学的方法。该模型易于在粒子模拟中实现，用于在真实轨道条件下卫星周围的稀薄流动。鉴于上述特征，当使用相关表面和流动条件下的分子束-表面散射实验或高保真分子动力学模拟进行参数化后，该模型预计能比传统模型（如CL或Maxwell模型）更准确地描述气体-表面动量交换（气动阻力的基本驱动因素），从而改进VLEO中卫星的阻力预测。