一家精密光学企业正在验收两批镜片。供应商A提供的检测报告Surface Figure\[ \lambda/10 \]RMS18 nm。供应商BSurface Figure同样是\[ \lambda/10 \]RMS17 nm。从采购规格来看两批镜片几乎没有区别。甚至B还稍微好一点。于是工程团队分别装机测试。结果却完全不同。A批镜片装入系统以后PSF正常MTF接近设计值成像对比度符合要求。B批镜片中心Spot并没有明显失焦但PSF外围出现了一圈奇怪的能量裙边MTF在某些空间频率明显下降强光目标附近甚至出现额外背景。工程师重新检查曲率没问题中心厚度没问题偏心没问题装调也没有明显异常。于是一个非常现实的问题出现了两块镜片的RMS几乎一样为什么系统性能会完全不同问题很可能并不在误差有多大。而在误差长什么样。一、RMS只告诉你“有多少误差”没有告诉你“误差在哪里”假设有两块直径100 mm的光学表面。第一块表面的误差是一条非常缓慢的起伏\[ z_1(x)A\sin(2\pi x/100) \]第二块表面则存在周期约2 mm的加工纹\[ z_2(x)A\sin(2\pi x/2) \]如果适当调整振幅两者完全可能拥有相同的\[ RMS(z_1)RMS(z_2) \]但是它们对光学系统的影响并不相同。第一种误差更接近PowerAstigmatismComa低阶Figure Error。它主要改变焦点像差PSF核心。第二种误差属于更高的空间频率。它会把光场能量重新分配到不同角度。于是可能出现PSF haloMTF下降散射增加结构化伪影。这就是为什么现代高性能光学系统越来越不能只使用一个总RMS数字描述表面质量。NASA 的大口径光学规范研究明确把表面误差划分到不同空间频率区域低空间频率误差主要改变PSF核心及其附近结构中频误差会使PSF核心扩散而高频误差和粗糙度则把能量散射到更宽的范围。NASA技术报告服务器因此\[ RMS \]实际上只是一个积分后的结果。它把“误差发生在哪里”这个极其重要的信息丢掉了。二、真正应该问的是误差存在于哪个Spatial Frequency把表面高度写成\[ h(x,y) \]经过二维傅里叶变换以后可以观察它的空间频率组成\[ H(f_x,f_y) \mathcal{F}\{h(x,y)\} \]进一步得到\[ PSD(f_x,f_y) \propto |H(f_x,f_y)|^2 \]PSD也就是Power Spectral Density。它回答的不是这块镜片一共有多少纳米误差而是这些纳米分别分布在哪些空间尺度上例如100 mm周期10 mm周期1 mm周期100 μm周期10 μm周期。这些误差即使拥有相同RMS对系统造成的结果也可能完全不同。这也是为什么 International Organization for Standardization 的 ISO 10110-8 不仅允许使用传统RMS粗糙度描述表面纹理也提供了用于控制中空间频率waviness的PSD、RMS slope等方法。SPIE但这里必须特别强调一点Low / Mid / High Spatial Frequency 并不存在一组适用于所有系统的绝对边界。NASA的大口径空间望远镜工程资料也明确指出不同频段之间并不存在普遍接受的固定分界边界会随科学目标和制造工艺而变化。NASA技术报告服务器所以工程上更正确的思路不是背一个“0.01–1 mm⁻¹就是MSF”这样的数字。而应该问对于我的这个光学系统哪些空间频率最危险三、为什么中频误差尤其容易成为“漏网之鱼”这是现代精密光学制造里一个非常有意思的问题。传统加工时代人们主要关心两个极端。一个是Figure。例如PowerIrregularityPVRMS。另一个是Roughness。例如0.5 nm RMS1 nm RMS。于是规格体系很容易形成“大尺度面形 微观粗糙度”两头管理。但现代制造方式越来越大量采用CNC polishingsmall-tool polishingbonnet polishingMRFrobotic polishingsub-aperture deterministic finishing。这些工艺非常擅长确定性修正低频面形。但“小工具沿着表面移动”本身也可能留下具有特定周期的Tool MarksRaster PatternSpiral PatternRipple。它们恰好落在Figure 与 Roughness之间。也就是Mid-Spatial FrequencyMSF。2017年Applied Optics的实验研究指出仅使用传统微观RMS粗糙度并不足以完整描述先进光学表面毫米量级的空间误差已经可能增加散射背景并降低大型光学系统的成像对比度。研究进一步把测得的PSD与研磨方式、抛光界面和抛光材料等制造参数建立了联系。光学出版集团2024年另一项BK7镜片研究则直接针对bonnet polishing产生的MSF进行PSD分析并通过优化spiral tool path显著降低相关PSD误差。数字对象标识符而这件事并没有因为自由曲面加工成熟而消失。相反自由曲面更加容易暴露这个问题。2023年的自由曲面研究甚至发现MSF并不一定只来自加工和检测阶段高阶XY polynomial表面的设计数据采样和数据处理本身也可能产生可识别的中频结构。施普林格自然链接因此今天一块镜片的误差链可能是\[ Design\ Data \]↓\[ Tool\ Path \]↓\[ Material\ Removal \]↓\[ Surface\ Ripple \]↓\[ Wavefront \]↓\[ PSF \]这已经是一个非常典型的系统级工程问题。四、真正重要的不是PSD好不好看而是它最后把光送到了哪里到这里很容易再次掉进另一个陷阱以前只看RMS。现在改成只看PSD。这仍然不够。因为企业真正关心的不是“PSD曲线是否漂亮”而是这个PSD经过真实光学系统以后会产生什么对于反射面表面高度误差会形成近似\[ \Delta OPD\approx2h \]的光程变化。对于透射面情况则与折射率入射角表面位置光束footprint共同相关。于是需要把真实测得或者统计生成的\[ h(x,y) \]放进系统。得到\[ OPD(x,y) \]进一步形成\[ E_{pupil}(x,y) A(x,y) e^{i\frac{2\pi}{\lambda}OPD(x,y)} \]然后通过物理光学传播得到\[ E_{image}(x,y) \]最终\[ PSF|E_{image}|^2 \]再由PSF得到MTFEncircled EnergyContrastHaloPeak Irradiance。2010年的经典研究已经专门分析结构化MSF如何产生structured image artifacts并强调应在真实光学系统中对MSF进行评价和容差分析而不是脱离系统单独判断表面。光学出版集团更值得注意的是2023年的研究指出传统MSF表面规格经过简化后可能丢失复杂topography的重要信息因此仅靠常规表面规格很难直接建立与系统relative modulation之间的对应关系。光学出版集团这恰好说明Surface Specification和System Performance之间还缺少一座桥。这座桥就是Physical Optical Propagation。五、同样20 nm RMS我们可以做一个Virtual A/B Test现在做一个非常有工程意义的虚拟实验。假设三个表面A20 nm RMS Low-Frequency Error。B20 nm RMS Mid-Frequency Ripple。C20 nm RMS High-Frequency Random Error。注意三者RMS完全一样。然后把三张surface map分别加载到同一套成像系统中。保持光源镜头装调孔径探测器全部不变。只改变Error Spectrum。最终很可能得到完全不同的结果。AStrehl下降PSF core发生变化低阶像差明显。B中心峰值下降PSF周围形成结构化halo中高频MTF受到明显影响。C中心结构变化可能没有B那么明显但宽角背景和散射增加。这时候我们第一次可以真正回答对于这套系统20 nm误差应该出现在哪里甚至进一步做\[ KPIF(Spatial\ Frequency) \]例如扫描一系列不同频率0.1 cycles/mm0.20.512510……每个频率都保持相同surface RMS。然后计算StrehlMTF50Encircled EnergyHalo。最终得到System Sensitivity vs Spatial Frequency。这张曲线可能比传统“Surface RMS 20 nm”规格有价值很多。六、这时候企业才能真正回答供应商应该做到什么程度这其实是这类仿真最有商业价值的一步。假设光学设计部门给供应商的图纸只写Surface Figure\[ \lambda/10 \]Roughness\[ 1\,nm\ RMS \]供应商完全可能同时满足两个指标。但是中间仍然留着一个巨大的空间MSF。于是双方都会说“我的指标是合格的。”而整机却不合格。更合理的方法是从系统KPI反推元件规格。例如产品要求\[ MTF(100lp/mm)0.45 \]或者\[ Strehl0.8 \]或者\[ Halo(5\lambda/D)10^{-4} \]那么系统仿真可以逐步增加某个空间频带的surface error直到系统刚好触碰KPI边界。得到\[ PSD_{limit}(f) \]于是规格不再来自“行业一般做到多少”而来自我的系统到底允许多少NASA/JPL很早就提出使用模拟结果推导大口径光学表面的PSD specification再将实际测得surface profile与预测PSD进行比较这实际上就是“系统需求→PSD规格→制造测量”的闭环思路。NASA技术报告服务器ISO 10110-8提供PSD等表面纹理描述方式则让这种系统需求进一步有可能转化为制造图纸和验收语言。SPIE这时候光学仿真的角色已经改变。它不再只是“检查这个镜头成像怎么样。”而开始参与供应商规格定义。七、而且PSD还能告诉制造部门应该改哪一个工艺参数故事还没有结束。如果检测发现某一频率出现明显PSD peak。例如\[ f0.5\,mm^{-1} \]意味着对应空间周期大约\[ P\frac{1}{f}2\,mm \]那么工程师就可以开始追查什么制造动作具有约2 mm的空间周期Tool sizeRaster spacingFeed ratePolishing pad grooveMachine vibrationPath overlap这时候PSD不再只是Quality Metric。它开始变成Manufacturing Fingerprint。2019年发表于Precision Engineering的研究建立了全口径抛光运动学和pad groove参数与MSF之间的模型并通过实验验证groove spacing会影响MSF周期而调整polishing speed和swing参数可以显著降低PSD-RMS。科学直通车针对MRF的国内研究也观察到类似关系磁流变抛光带的波动会影响高功率激光元件表面的MSF其波动量与特定中频段PSD RMS表现出近似线性关系因此可以反向通过工艺参数控制表面中频误差。PubMed Central (PMC)而NASA 2026年更新的自由曲面制造项目仍然把MSF作为核心制造问题之一并指出sub-aperture polishing本身容易直接形成MSF因此需要额外的full-aperture smoothing过程进行抑制。NASA TechPort于是完整链条变成\[ Machine\ Parameter \]↓\[ Tool\ Path \]↓\[ Surface\ PSD \]↓\[ Wavefront \]↓\[ PSF/MTF \]↓\[ Product\ Performance \]这才是真正有意义的Manufacturing-to-Performance Model。八、对于高功率激光系统问题甚至不只是“成像变差”如果这块光学元件不是用于普通成像而是用于高功率激光EUV高对比度望远镜光刻精密计量MSF的后果会进一步变化。例如高功率激光系统中周期性表面误差会调制传播光场。经过一定距离以后phase modulation可能逐渐转化为intensity modulation。于是某些位置形成Hot Spot。这时系统真正关心的KPI已经不再是MTF。而可能变成\[ I_{max}/I_{avg} \]或者Near-field ModulationFar-field EnergyDamage Risk。公开的高功率激光研究指出元件MSF会增加光束调制度并可能通过非线性自聚焦进一步威胁系统安全。PubMed Central (PMC)在空间望远镜领域同一个问题又会变成另一种形式。NASA的大口径UVO空间望远镜工程研究指出中高空间频率误差会把恒星光从PSF核心重新分配出去这对于需要极高对比度的系外行星成像尤其关键。NASA技术报告服务器而在EUV光刻中MSF相关散射则可能表现为FlareContrast LossProcess Window下降。公开的EUV光学资料将flare描述为PSF核心周围增加的能量“skirt”本质上同样是能量从理想位置被重新分配。半导体实验室所以同一张PSD曲线没有脱离系统的绝对“好”与“坏”。必须问它处在什么系统什么波长什么孔径什么光束footprint最终KPI是什么这再次回到了系统级仿真的核心。九、真正值得建立的是“Surface Error Budget”到了这里我们可以重新定义传统的Tolerance Budget。以前可能写Lens 120 nm RMSLens 215 nm RMSMirror 110 nm RMS。然后做RSS\[ \sigma_{total} \sqrt{ \sigma_1^2 \sigma_2^2 \sigma_3^2 } \]但更先进的系统应该进一步问这些误差分别在哪个空间频带于是Error Budget变成\[ PSD_{system}(f) \]或者更工程化地划分LSF BudgetMSF BudgetHSF Budget。然后再考虑每一个表面的Beam FootprintConjugate PositionMagnificationIncidence Angle。因为同样一个2 mm周期的表面纹理放在系统不同位置对最终像面的映射可能不同。2010年的结构化MSF研究在实际测试镜头中就专门分析了不同表面的beam footprint与可见空间频率范围这正说明MSF必须结合具体系统位置判断而不能只看单独元件。光学出版集团最终得到的不是“每块镜片都做到最好。”而是把制造成本花在系统真正敏感的地方。十、所以“λ/10”并不是故事的终点回到文章开始的两批镜片。A18 nm RMS。B17 nm RMS。为什么17 nm反而更差现在答案已经很清楚因为\[ 17\,nm \]只是一个总量。真正决定系统性能的是\[ Error(f_x,f_y) \]以及它经过整个光学系统之后把能量送到了哪里。因此一个更加完整的工程判断应该是\[ Surface\ Map \]↓\[ PSD \]↓\[ System\ Sensitivity \]↓\[ Wavefront \]↓\[ Physical\ Propagation \]↓\[ PSF/MTF/Scatter \]↓\[ Product\ KPI \]这时候我们才能真正回答这块镜片为什么失败哪个空间频率最危险哪一个表面最敏感供应商规格应该怎么写加工设备应该改哪个参数返修应该优先去掉哪一段频率下一代产品要不要提高这个指标十一、这正是系统级物理光学仿真最容易被低估的价值很多时候我们以为仿真的价值是把一个光学系统设计出来。但真正进入量产以后企业面对的问题会慢慢发生变化。设计部门问我的设计对什么surface frequency最敏感制造部门问这个tool mark到底有没有关系检测部门问为什么RMS合格整机仍然失败采购部门问供应商到底应该做到什么指标质量部门问这块镜片能不能放行这些问题实际上都指向同一个缺口元件检测指标与最终系统性能之间没有建立物理联系。系统级物理光学仿真的价值就是把这条链连接起来\[ Manufacturing \rightarrow Surface \rightarrow Wavefront \rightarrow Field \rightarrow System \rightarrow Product \]于是最终得到的并不是一张更复杂的PSD图。而是一套从产品性能反推光学制造规格再从真实测量数据预测整机表现的工程模型。对于真正进入高性能、高精度和批量制造阶段的光学企业这种模型的价值往往比单纯再把设计MTF提高几个百分点更大。因为它开始回答另一个更现实的问题设计出来的光学系统怎样才能真正被稳定地制造出来。