什么是泽尼克多项式

泽尼克多项式的历史

泽尼克多项式以荷兰物理学家弗里茨·泽尼克（Frits Zernike，1888-1966）命名。泽尼克因发明相位对比显微镜技术而于 1953 年荣获诺贝尔物理学奖。他在研究光学系统中的相位问题时，提出了一组定义在单位圆上的正交多项式，用于展开和分析光瞳面上的波前分布。

这组多项式因其优良的数学性质和与经典光学像差的天然对应关系，迅速成为光学像差理论和检测领域的标准工具。如今，泽尼克多项式的应用已远远超出传统光学设计的范畴，扩展到天文学、眼科学、半导体光刻、大气物理等诸多学科。

为什么需要泽尼克多项式

光学系统的性能最终由波前的形状决定。理想光学系统产生完美的球面波前（会聚到一个几何点），而实际系统的波前偏离理想球面——这种偏离就是"波前像差"。要定量描述波前像差，就需要一种数学表示方法。

原则上可以用任何正交函数系来展开波前，但泽尼克多项式因以下几个特殊优势脱颖而出：

• 定义域匹配：大多数光学系统的光瞳（Pupil）为圆形，泽尼克多项式恰好定义在单位圆上，与光瞳的几何形状完美匹配。
• 与经典像差一一对应：低阶泽尼克项直接对应人们熟悉的像差类型（离焦、像散、彗差、球差等），物理含义直观。
• 正交完备：各项之间正交，展开系数相互独立，便于分离和量化每种像差的贡献。
• 逼近效率高：通常只需前 15~36 项就能精确表示绝大多数实际光学系统的波前像差。

数学定义

泽尼克多项式定义在单位圆（半径 ρ ≤ 1）上，采用极坐标 (ρ, θ) 表示。每个泽尼克项由两个指标确定：径向阶数 n（非负整数）和角频率 m（满足 |m| ≤ n，且 n - |m| 为偶数）。

泽尼克多项式的一般形式为：

其中径向多项式 R_n^|m|(ρ) 定义为：

角频率 m 决定了方位方向的变化周期：m = 0 为旋转对称项（无角度依赖），|m| = 1 为单周期变化（如倾斜），|m| = 2 为两周期变化（如像散），以此类推。径向阶数 n 决定了径向方向的复杂程度：n 越大，径向振荡越多。

归一化

为使各项在单位圆上的均方根值归一化为 1，通常采用归一化系数 N_n^m。归一化后的泽尼克多项式使得每一项的系数直接代表该像差的 RMS 值（单位与波前误差相同，通常为波长 λ 的分数或纳米），这在实际应用中非常方便。

低阶项与经典像差对应

泽尼克多项式的最大魅力之一是其低阶项与经典赛德尔像差之间存在清晰的对应关系。下表列出了前 15 项的对应关系：

在实际使用中，Piston（Z1）和 Tilt（Z2、Z3）通常被视为"非像差项"——Piston 是整体相位偏移，对成像无影响；Tilt 只引起像的平移，不影响像质。因此，评价光学系统波前质量时通常从 Z4（离焦）开始计算。

正交性与完备性

泽尼克多项式在单位圆上满足正交性条件：

这一正交性质带来了几个重要的实际优势：

• 系数独立性：增加或去除某一项不会影响其他项的系数值。例如，在波前拟合中增加高阶项不会改变已经确定的低阶系数。
• 最小二乘拟合的唯一性：对于给定的波前数据，最小二乘拟合得到的泽尼克系数是唯一的。
• RMS 分解：波前的总 RMS 值等于各项 RMS 贡献的均方根和，即 σ_total² = Σ a_i²，便于评估各像差的相对重要性。

完备性保证了泽尼克展开可以在单位圆上逼近任意连续波前函数，只要展开到足够高的阶数。

Noll / OSA / Fringe 排序方案

由于泽尼克多项式的二重指标 (n, m) 需要映射为单一序号以便于使用，历史上出现了多种排序方案：

Noll 排序

由 Robert Noll 在 1976 年提出，是天文和大气光学领域最广泛使用的排序方式。按径向阶数 n 递增排列，同一 n 值内按 |m| 递增排列，正余弦项交替。

OSA/ANSI 排序

由美国国家标准协会 (ANSI) 和美国光学学会 (OSA) 联合规范（ANSI Z80.28 标准），主要用于眼科光学。按 n 递增排列，同一 n 值内按 m 从负到正递增。

Fringe（University of Arizona）排序

由亚利桑那大学光学科学中心提出，被 Zygo 等干涉仪厂商广泛采用。其特点是只包含 37 项（最多到 n = 12 的部分项），按"对角线"方式排列，在工业光学检测中最为常用。

使用泽尼克多项式时，最重要的是明确所用的排序方案和归一化约定，否则不同来源的系数无法直接比较。

应用：干涉检测

泽尼克多项式在光学干涉检测中的应用是其最经典的用途之一。

干涉仪波前分析

当使用 Fizeau 或 Twyman-Green 干涉仪测量光学元件的面形误差时，干涉条纹图包含了被测面与参考面之间的波前差信息。通过相移干涉技术提取波前相位后，对波前数据进行泽尼克拟合，即可将复杂的二维波前分解为各阶像差分量。

例如，Zygo 干涉仪的标准软件输出中通常包含 Fringe 排序的前 36 项泽尼克系数。工程师可以直接从系数中识别出：

• Z4（Defocus）：透镜的焦距偏差或检测时的焦点未对准
• Z5/Z6（Astigmatism）：柱面误差，可能源于装夹应力或加工不均匀
• Z7/Z8（Coma）：偏心误差或倾斜
• Z11（Spherical）：球差，反映曲率半径偏差或高阶面形误差

去除低阶项（Piston、Tilt、Defocus 等装调相关项）后的残余 RMS 即为元件的实际面形精度，通常以 λ/N（波长的几分之一）表示。高精度光学元件要求面形精度达到 λ/20 甚至 λ/100。

应用：自适应光学

在天文自适应光学（AO）系统中，泽尼克多项式提供了描述大气湍流引起的波前畸变的天然框架。

大气湍流的泽尼克描述

Kolmogorov 大气湍流模型下，波前各泽尼克项的统计特性（方差和相关性）有解析表达式。Robert Noll 在 1976 年的经典论文中计算了各阶泽尼克模式的方差——低阶模式（Tilt、Defocus）的方差最大，贡献了大部分波前误差；高阶模式方差按幂律递减。这意味着，自适应光学系统只需校正有限数目的低阶模式即可显著改善成像质量。

可变形镜控制

自适应光学系统中的可变形镜（DM）通常用泽尼克模式来定义其校正能力。波前传感器（如 Shack-Hartmann 传感器）测量入射波前后，控制系统将波前分解为泽尼克模式，然后驱动 DM 上的促动器产生相应的补偿面形。低阶 AO（如 Tip-Tilt 系统）只校正 Z2 和 Z3，高阶 AO 可校正数百甚至数千个模式。

应用：眼科光学

泽尼克多项式在眼科学中有着至关重要的应用，用于描述和分析人眼的光学像差。

Hartmann-Shack 波前传感器

现代眼科波前像差仪（如 Zywave、iDesign 等）使用 Hartmann-Shack 波前传感器测量瞳孔平面上的波前。传感器由微透镜阵列和 CCD 组成，每个微透镜将局部波前聚焦为一个光点，通过分析各光点相对于参考位置的偏移，可以重建整个瞳孔上的波前，再用泽尼克多项式展开。

屈光矫正

传统的眼镜和隐形眼镜只能矫正低阶像差：Z4（离焦，即近视/远视）和 Z5/Z6（像散，即散光）。而人眼还存在高阶像差——彗差、球差、三叶草等。泽尼克分析可以精确量化每种像差的大小，为个性化屈光手术（如波前引导的 LASIK）提供数据基础。

波前引导 LASIK

在波前引导的 LASIK 手术中，准分子激光根据泽尼克像差数据进行个性化角膜切削。激光切削的面形被设计为刚好补偿患者的高阶像差，使术后的残余波前误差降至最低。这种技术相比传统 LASIK 可以获得更好的视觉质量，尤其在暗光环境（瞳孔放大、高阶像差影响增大）下的视力表现更为明显。

应用：光学设计与光刻

光学设计软件

在 Zemax OpticStudio、Code V 等光学设计软件中，波前分析功能直接以泽尼克系数形式输出结果。设计师可以通过查看各阶泽尼克系数来快速识别主导像差类型，并针对性地调整光学参数。例如，如果 Z11（球差）系数过大，可以通过调整非球面参数或增加透镜来校正；如果 Z7/Z8（彗差）过大，可能需要检查系统的对称性或光阑位置。

半导体光刻

在半导体光刻领域，投影物镜的波前质量直接决定了可以达到的最小特征尺寸。先进的 DUV/EUV 光刻机投影物镜要求波前精度达到 mλ（毫波长）级别。泽尼克分析在光刻物镜的设计、制造和装调全过程中被广泛使用：

• 设计阶段：以泽尼克系数作为优化目标，确保全视场的波前质量均匀。
• 制造检测：对每个光学元件进行干涉检测，以泽尼克系数衡量面形精度。
• 装调补偿：通过调整元件间距和偏心来最小化系统的总泽尼克像差。
• 在线监控：在光刻机运行过程中，实时测量和补偿因温度变化等引起的波前漂移。

ASML 等光刻机厂商使用泽尼克多项式（通常采用 Fringe 排序的前 64 或更多项）作为描述投影物镜光学性能的标准语言，各芯片制造商根据泽尼克系数来验收光刻机性能。