4、去除函数模型：去除函数定义、线性假设与影响去除率的因素

离子束去除（Ion Beam Figuring, IBF）的核心在于精确控制材料去除的空间分布与时间累积。去除函数（Removal Function）是描述单位时间内离子束在工件表面某一点所去除材料体积或深度的数学表达，它是连接工艺参数（束流能量、束流密度、入射角等）与面形修正结果之间的桥梁。本章将深入剖析去除函数的定义、线性假设的适用边界，以及束流能量、束流密度、入射角和靶材材料对去除率的影响机理，为后续的驻留时间算法与工艺优化奠定理论基础。

4.1 去除函数的定义与物理内涵

去除函数在数学上定义为：在离子束以恒定工艺参数（束流能量、束流密度、入射角等）作用下，单位时间内工件表面各点被去除的材料深度分布，通常记为 R(x, y)，单位为 nm/min 或 μm/min。更严格地，它描述的是离子束在工件表面扫描驻留时，材料去除速率在空间上的脉冲响应。在工程实践中，去除函数通常通过实验标定获得：使用一束固定参数的离子束，在平面样件上轰击一定时间，然后通过干涉仪或轮廓仪测量去除坑的形貌，再除以轰击时间，即得到去除函数。

去除函数的物理本质是离子束与靶材原子碰撞级联过程的宏观统计结果。当高能离子（通常为 Ar⁺、Kr⁺ 或 Xe⁺，能量在 300 eV ~ 1500 eV 范围）入射到靶材表面时，通过弹性碰撞将动量传递给靶材原子，形成碰撞级联。当级联中的反冲原子获得足够能量并运动到表面时，若其动能超过表面结合能，即可逸出表面，形成溅射。去除函数正是这种溅射过程在空间上的积分效应，其形状通常近似为高斯分布，中心去除速率最高，边缘逐渐衰减。

4.2 线性假设：适用条件与边界

离子束去除的线性假设是经典 IBF 工艺的理论基石，它包含三个基本前提：时不变性、空间不变性和线性叠加性。时不变性指去除函数在加工过程中不随时间变化，即束流参数稳定、靶材特性不因加工而改变；空间不变性指去除函数在工件表面任意位置具有相同的形状和去除效率，不随扫描位置变化；线性叠加性指总去除量等于各驻留点去除贡献的线性叠加，即去除深度与驻留时间成正比。

在实际工程中，线性假设的成立需要满足以下条件：束流能量和束流密度在加工过程中保持稳定（波动通常需控制在 ±1% 以内）；靶材温度控制在合理范围内，避免热效应导致去除率漂移；工件表面形貌变化引起的入射角变化足够小（通常要求表面坡度小于 5°~10°）；以及靶材材料均匀、无晶向依赖或缺陷。当这些条件满足时，去除过程可以表示为卷积模型：

E(x,y) = R(x,y) ⊗ T(x,y) = ∫∫ R(x-x', y-y') · T(x',y') dx' dy'

其中 E(x,y) 为面形误差修正量，T(x,y) 为驻留时间分布。该卷积模型是反卷积求解驻留时间的数学基础。

4.3 影响去除率的关键因素

去除率（即去除函数的峰值高度）受多个工艺参数控制，理解其影响规律是优化加工效率与精度的前提。以下从束流能量、束流密度、入射角和靶材材料四个维度展开分析。

4.3.1 束流能量（E）

束流能量决定了入射离子传递给靶材原子的最大动能，直接影响溅射产额（Y，即每个入射离子溅射出的靶材原子数）。在低能区（< 100 eV），溅射产额极低，接近阈值；随着能量升高，溅射产额迅速增加，在 100 eV ~ 1 keV 范围内近似线性增长；当能量超过约 10 keV 时，溅射产额趋于饱和甚至略有下降，这是因为高能离子更多地注入靶材内部，能量沉积在更深处，表面溅射效率反而降低。对于光学加工常用的 Ar⁺ 离子，典型工作能量范围为 300 eV ~ 1000 eV，此时溅射产额约为 0.5 ~ 3 原子/离子。

束流能量还影响去除函数的空间分布。较高能量导致碰撞级联深度增加，去除函数的半高全宽（FWHM）增大，形状更“宽”而“浅”；较低能量则产生更集中、更尖锐的去除峰。因此，在粗加工阶段可选择较高能量以提高效率，在精修阶段降低能量以获得更精细的空间分辨率。

4.3.2 束流密度（J）

束流密度（单位：mA/cm²）直接决定了单位时间内到达工件表面的离子数量，与去除率呈近似线性正比关系。去除率峰值 R_peak 可近似表示为：

R_peak = (Y · J · M_atom) / (ρ · N_A · e)

其中 M_atom 为靶材原子摩尔质量，ρ 为靶材密度，N_A 为阿伏伽德罗常数，e 为元电荷。在实际系统中，束流密度受限于离子源的最大引出电流和束流聚焦能力。增大束流密度可显著提高加工效率，但过高的束流密度可能导致靶材表面温度升高、微裂纹产生或束流不稳定，需权衡选择。

4.3.3 入射角（θ）

入射角定义为离子束方向与靶材表面法线之间的夹角。溅射产额随入射角的变化呈现典型的“角度依赖”特性：在法线入射（θ=0°）时，溅射产额并非最大；随着入射角增大，溅射产额逐渐增加，在约 60°~80° 时达到峰值（通常为法线入射的 2~4 倍）；超过峰值角度后，溅射产额急剧下降，因为离子在表面发生镜面反射的概率增大，能量传递效率降低。

这一角度依赖特性对 IBF 加工有重要影响。对于平面或低陡度曲面，入射角变化较小，可采用固定角度标定的去除函数；对于高陡度曲面，需考虑入射角变化对去除函数形状和峰值的影响，必要时采用角度分段模型。此外，入射角还会影响去除函数的对称性——斜入射时去除函数不再呈轴对称高斯分布，而是向离子束方向偏移，这增加了驻留时间求解的复杂度。

4.3.4 靶材材料特性

靶材材料的物理性质对去除率有决定性影响。主要因素包括：表面结合能（升华热）、原子序数（质量）、晶体结构和密度。表面结合能越低，原子越容易逸出，溅射产额越高；靶材原子质量与入射离子质量比越接近 1，动量传递效率越高，溅射产额越大。例如，对于 Ar⁺ 离子（质量 40 amu），轰击铜（63.5 amu）的溅射产额约为轰击碳（12 amu）的 5~10 倍。

晶体结构的影响体现在单晶靶材的晶向依赖效应——不同晶面的溅射产额差异可达 20% 以上。对于光学玻璃（如熔石英、微晶玻璃）等非晶材料，不存在晶向依赖，但材料内部缺陷、气泡或应力区域可能导致去除率局部波动。此外，靶材的导热性影响加工过程中的温度分布，进而间接影响去除稳定性。

4.4 案例拆解：熔石英样件去除函数标定与参数优化

背景与目标：某光学工厂需对一块直径 100 mm 的熔石英平面镜进行 IBF 精修，目标面形 RMS 从初始 15 nm 收敛至 2 nm 以内。工艺工程师需先标定去除函数，并确定最优工艺参数组合。

初始条件与约束：设备为射频离子源，可调参数范围：束流能量 300~1000 eV，束流密度 0.2~2.0 mA/cm²，入射角固定为 0°（垂直入射）。要求单次加工时间不超过 30 分钟，且去除深度不超过 200 nm 以避免引入亚表面损伤。

标定实验设计：采用正交实验法，选取 3 个能量水平（400、600、800 eV）和 3 个束流密度水平（0.5、1.0、1.5 mA/cm²），共 9 组参数。每组参数在 3 块 25 mm × 25 mm 的熔石英样件上轰击 5 分钟，使用白光干涉仪测量去除坑深度分布。

关键数据与结果：

• 在 600 eV、1.0 mA/cm² 条件下，去除函数峰值速率约为 12.5 nm/min，FWHM 为 8.2 mm，形状接近高斯分布（拟合残差 RMS < 0.3 nm）。
• 束流密度从 0.5 增至 1.5 mA/cm² 时，峰值去除率从 6.3 增至 18.7 nm/min，线性度良好（R² = 0.998）。
• 束流能量从 400 增至 800 eV 时，峰值去除率从 8.1 增至 15.2 nm/min，但 FWHM 从 6.5 mm 增至 9.8 mm，空间分辨率下降。

分析与决策：考虑到目标面形误差以中频误差为主（空间周期 5~15 mm），需兼顾去除效率与空间分辨率。选择 600 eV、1.2 mA/cm² 作为精修参数，此时峰值去除率约 15 nm/min，FWHM 约 8.5 mm，可在 25 分钟内完成预期去除量。同时，为验证线性假设的适用性，在样件上进行了双点叠加实验：两点间距 10 mm，驻留时间比 1:2，实测叠加去除深度与线性预测偏差小于 2%，确认线性假设在 200 nm 去除深度内成立。

结果与复盘：采用标定参数进行全口径加工，经过 3 次迭代，面形 RMS 从 15 nm 收敛至 1.8 nm，达到预期目标。复盘经验：① 去除函数标定应在与工件相同材料、相同表面状态的样件上进行；② 每次加工前需重新标定去除函数，以消除离子源状态漂移的影响；③ 对于高精度要求，建议在正式加工前进行小面积试加工，验证去除函数的稳定性。

4.5 常见误区与工程建议

总结而言，去除函数模型是离子束加工的核心工具，其准确性直接决定了面形修正的收敛速度与极限精度。深入理解线性假设的适用条件，掌握各因素对去除率的影响规律，并通过规范的标定流程获取可靠的去除函数数据，是每一位 IBF 工艺工程师必须掌握的核心技能。