做结构光场的人，近几年大概都会听到一个绕不开的名字——Ince-Gaussian光束。如果你正好拿到这个课题，要看明白涡旋阵列是怎么从这类光束里“长”出来的，我的建议是别急着上手调光路，先花两天把Ince-Gaussian的理论框架和仿真吃透，再做实验就是按图索骥的事。这篇文章我就从理论建模、光路搭建到涡旋阵列的观测判读，完整走一遍流程，最后把我踩过的坑和排查思路一并抖出来。
1. 为什么偏偏是Ince-Gaussian：涡旋阵列的思路来源
1.1 Ince-Gaussian光束到底是什么
Ince-Gaussian光束是傍轴波动方程在椭圆坐标系下的一组正交完备解。你要理解它，可以把它想象成高斯光束家族里最“通吃”的一类：当椭圆坐标退化到极限情况，它可以连续过渡到厄米-高斯光束和拉盖尔-高斯光束。换句话说，Ince-Gaussian不是和它们并列的一种新光束，而是能包住它们的一种更一般的解族。这就像四边形里的矩形，矩形里的正方形——Ince-Gaussian是那个大集合，厄米-高斯和拉盖尔-高斯是它两端的特例。
它的场分布一般写成椭圆坐标下的模式叠加形式，核心是两个因子：一个是长椭球方向的Ince多项式，另一个是短轴方向的Ince多项式，两者由一组模式阶数(p, m)和椭圆参数ε共同决定。这里有个地方特别值得注意：Ince多项式分偶模和奇模，也就是IG^e_{p,m}和IG^o_{p,m}，它们在椭圆坐标下本身就有复杂的相位结构。
1.2 为什么它能自然形成涡旋阵列
涡旋光束的核心是相位奇点，也就是某个位置上光场的相位无法定义，强度必须为零。单个涡旋好办，螺旋相位板就能解决。但要做涡旋阵列，也就是在一束光里同时产生十几个甚至几十个涡旋，常用的方法要么是多个光束干涉，要么用达曼光栅做分束，要么叠加多个不同拓扑荷的拉盖尔-高斯模式。这些方法都有个共同缺点：系统复杂，对对准精度极其敏感。
Ince-Gaussian的好处在于，它的涡旋阵列是“内生”的。研究已经证明，相同阶数、相同椭圆参数的一对偶模和奇模Ince-Gaussian模式，以合适的相位差做相干叠加时，会在横截面上产生一串离散的涡旋点，这些涡旋点规则地分布在一条椭圆弧线或闭环上。涡旋的数目和位置由模式阶数m直接决定，不需要额外引入任何干涉器件。这是非常漂亮的生成方式——一个空间光调制器，两页全息图，就能看到干净的涡旋环。
1.3 这个方案对比传统涡旋阵列生成方法的优势
我做这个课题之前，首先对比了一轮已知的涡旋阵列生成方案。达曼光栅方案的好处是均匀性好，但设计门槛偏高，相位台阶的加工精度直接决定阵列均匀度，实验室条件很难做到理想；多光束干涉方案概念直观，但需要极其稳定的干涉光路，对环境振动极其敏感，稍微有点气流条纹就会模糊；拉盖尔-高斯叠加方案理论简单，但实验上要精确控制多个模式的相对权重和相位，调起来很费时间。
Ince-Gaussian方案在这种对比里显得特别实用：你只需要用SLM加载偶模和奇模对应的全息图，调整它们之间的相移量，涡旋阵列就出现了。而且因为它是单光束方案，不涉及多臂干涉，抗振能力天然比干涉方案强很多。我实测下来，用最普通的台式光学平台，曝光时间压到10ms以下，条纹照样清晰锐利。
2. 理论建模与仿真：先让光束在电脑里“跑”一遍
2.1 数值模型怎么写
我一直强调，做光学实验前必须先有可靠的仿真，否则你连“对的实验结果长什么样”都不知道。Ince-Gaussian在近轴条件下其实没有闭式的简单解析表达式，标准做法是用Ince多项式的递推关系来数值构造场分布。具体实现上，可以用Python配合NumPy写一个专门的计算模块，或者直接用MATLAB的符号工具包。
Ince多项式本身的递推关系并不复杂。椭圆坐标下的微分方程经过分离变量后，可以得到一组带椭圆参数ε的常微分方程，它们的解就是Ince多项式。计算时先解一个三对角矩阵的本征值问题，得到展开系数，然后代入Ince多项式的级数表达式，就能算出横截面上的复振幅分布。这个过程我建议直接参考文献里的递推公式，自己推导一遍，基本不会花太多时间，但能帮你彻底理解每个参数在控制什么。
下面是我常用的仿真核心代码骨架，这段代码能生成任意阶数的偶模Ince-Gaussian场分布：
import numpy as np from scipy.special import eval_hermite import scipy.linalg as la def ince_polynomial_coeffs(p, m, eps, parity='e'): # 构造Ince多项式递推所需的矩阵 # 返回展开系数c[0..p] n = (p - m) // 2 N = n + 1 M = np.zeros((N, N)) if parity == 'e': # 偶Ince多项式，对应cos项 for k in range(N): M[k, k] = (2*k + m)**2 if k > 0: M[k, k-1] = eps * (2*k + m) * (2*k + m - 1) / 2 if k < N - 1: M[k, k+1] = eps * (2*k + m + 1) * (2*k + m + 2) / 2 else: # 奇Ince多项式，对应sin项 for k in range(N): M[k, k] = (2*k + m)**2 if k > 0: M[k, k-1] = eps * (2*k + m) * (2*k + m - 1) / 2 if k < N - 1: M[k, k+1] = eps * (2*k + m + 1) * (2*k + m + 2) / 2 # 求解本征值问题，取最小本征值对应的向量 eigvals, eigvecs = la.eigh(M) idx = np.argmin(np.abs(eigvals - 0)) coeffs = eigvecs[:, idx] return coeffs / np.max(np.abs(coeffs))这段代码的关键点在于，Ince多项式的系数本质上是一个三对角矩阵本征问题的解，系数向量取对应最小本征值的特征向量。当然这里为了演示做了简化，实际做研究时建议交叉验证一下，和已发表的模式图对比。
2.2 三个核心参数怎么选
Ince-Gaussian的模式形态由三个参数共同决定：模式阶数(p, m)和椭圆参数ε。我建议一次只扫一个参数，先固定p=7, m=3，然后观察ε从0.5变化到8时的模式演化。
这里有一个规律要记住：ε越大，模式越“像”厄米-高斯；ε趋近于0时，模式退化到拉盖尔-高斯。对涡旋阵列来说，ε的中等取值（1~4）是最理想的，因为这个区间内涡旋点的间距比较均匀，阵列的椭圆环形态最突出，也最容易在实验上被SLM忠实再现。ε取值太大，涡旋点会挤到椭圆短轴两端去；ε太小，阵列又接近圆环，做不出“椭圆弧”的效果。
束腰的选择同样重要。我的经验是束腰w0取1mm到1.5mm之间，既保证了在SLM有效区域内的强度分布足够完整，又不会因为边缘强度太低导致信号噪声变大。束腰太小，高阶模的旁瓣会被SLM孔径切掉，模式纯度直接崩；束腰太大，SLM的像素分辨率又不够，相位量化的误差会被放大。
2.3 偶模叠加奇模：涡旋阵列出现的判据
涡旋阵列不是随便叠加就有的。理论上的做法是：取相同(p, m)的一对偶模IG^e_{p,m}和奇模IG^o_{p,m}，以相对相位差φ叠加，复振幅写为：
E = IG^e_{p,m} + i·IG^o_{p,m}
这个叠加模式下，横截面里会出现2m个拓扑荷为+1的涡旋（或者与相位差相关的拓扑荷）。p的取值主要决定径向节线数，m决定涡旋的个数和排列。我用仿真验证过很多种组合，p=7, m=3时，叠加场里能看到清晰的6个涡旋点，沿一个椭圆环均匀分布；p=5, m=2时是4个涡旋，分布在椭圆内弧上，间距更大一些。
这里有个实操技巧：仿真时不要只看强度图，一定要把相位图也画出来。强度图上的暗点看起来可能只是强度极小值，只有相位图才能确认那里确实存在拓扑荷为1的相位奇点。我通常的做法是计算环绕每个候选涡旋点的相位卷绕变化量，如果一圈下来总变化量为±2π，才认定它是真正的涡旋。
2.4 仿真结果怎么解读
当你拿到第一张像样的仿真强度图时，你会发现一个有趣的现象：涡旋点并不是像拉盖尔-高斯涡旋那样中心一个完美圆环，而是椭圆形的暗核，周围有对称分布的旁瓣。这是Ince-Gaussian坐标系本身带来的“遗传特征”，不要慌张，这正是Ince-Gaussian涡旋阵列区别于其他方案的指纹。
我在仿真阶段还会顺便预测一个关键判据：阵列中相邻两个涡旋之间的相位差应该是一致的。这样实验上如果看到干涉条纹出现等间距的叉形结构，就说明涡旋阵列的相位拓扑是稳定有序的。
3. 实验光路搭建与关键器件标定
3.1 整体光路配置
实验光路我采用的是一套标准的单光束SLM方案，结构简单，复现性好。整体链路如下：
激光器（氦氖激光器，波长632.8nm）→ 半波片+偏振片 → 扩束镜组（20倍）→ 反射式空间光调制器SLM → 透镜L（f = 200mm）→ CMOS相机
SLM对入射偏振极其敏感，反射式LCOS型SLM通常只有s偏振分量被高效调制。所以半波片和偏振片不是可有可无的装饰，而是决定调制效率的关键。我在第一次实验中偷懒跳过偏振控制，结果衍射效率掉了接近40%，模式混叠很严重，后来加上偏振片才算正常。
扩束的目的是让激光束在SLM面板上的光斑直径充分覆盖全息图的主要区域。我的面板是1920×1080像素，像素间距8μm，有效区域约15.36mm×8.64mm，束腰扩到1.2mm后，在面板上能看到完整的高阶模节线结构。如果扩束倍数不够，相当于用一个小窗口去看大图案，截断效应会让模式变形。
3.2 空间光调制器：相位标定必须自己做
市面上常见的SLM出厂时会附带相位标定文件，但我强烈建议不要完全相信它——温度、波长、入射角都会影响灰度-相位映射曲线。我拿到SLM后的第一件事就是做一次Gamma校正。
做法不复杂：在SLM上加载均匀灰度图，用迈克尔逊干涉光路观察干涉条纹位移，逐步改变灰度值从0到255，记录每个灰度下的相位变化。把数据拟合成一条单调曲线，然后取反函数作为查找表。这个步骤做完，SLM的相位调制精度基本能控制在±0.1π以内，对涡旋实验来说完全够用。
如果想偷懒，也可以直接用SLM厂商自带的软件做校正，但一定要确认校正波长为632.8nm，不能拿着别的波长文件直接套用。
3.3 全息图设计：把Ince-Gaussian写到SLM上
这是整个实验流程里技术含量最高的一步。SLM只能直接调制相位，不能直接调制复振幅，所以需要把复振幅编码成纯相位全息图。最常用的方法是闪耀光栅相位叠加法：
把目标复振幅E(x,y)写成振幅A(x,y)和相位φ(x,y)的乘积，然后构造：
H(x,y) = mod(φ(x,y) + 2π·(u0·x + v0·y), 2π)
其中u0、v0是光栅的空间频率，目的是把+1级衍射光和零级光分离开，后续用光阑只让+1级通过。灰度值由相位值换算：
G(x,y) = H(x,y) / (2π) × 255
这里有一个重要细节：纯相位编码会让振幅信息部分丢失。如果直接按原始IG模式去编码，阵列的均匀性会受影响。我采用的方法是先对振幅做归一化处理，再对振幅取一个阈值下限，低于阈值的相位直接随机化，这样能把模式纯度控制在可接受范围。更严格的做法是用复振幅编码算法，比如GS算法或者双相位编码，但对涡旋阵列来说，闪耀光栅叠加法已经足够好。
一个额外提醒：全息图的光栅周期不能太密，否则+1级和零级的分离角度太大，透镜后焦面上可能装不下衍射级次；太疏则分离不开。我常用的光栅周期是8个像素到12个像素，效果都比较稳定。
3.4 透镜选型与傅里叶变换位置
透镜L的作用是把SLM面上的相位分布做傅里叶变换，在后焦面上形成模式的远场强度分布。焦距的选择直接影响观测到的光斑大小和分辨率。
焦距200mm是个稳妥的起步值。焦距太长，光斑大但亮度低，相机动态范围容易吃紧；焦距太短，涡旋点之间的间距变小，容易挨在一起分辨不清。如果阵列的涡旋点间距普遍偏小，可以换更长焦距的透镜，把空间尺度拉开。
放置SLM时还需要考虑一个几何关系：SLM面板本身有一定倾斜角度，反射光和入射光之间存在夹角，这个夹角带来的像散会让椭圆参数ε偏离设定值。我在实验中发现，当入射角大于10度时，模式的椭圆取向会发生明显旋转，涡旋环从正椭圆变成了斜椭圆。解决办法是尽量减小入射角，或者在数据处理时把实测的椭圆取向和仿真一致化。
4. 涡旋阵列的观测与拓扑荷识别
4.1 强度分布观测：别满足于“看到暗点”
把加载了叠加模式全息图的SLM放进光路，调节相机曝光到不饱和的合适范围，你会在屏幕上看到一组由暗点组成的椭圆环。这是最直观的涡旋阵列证据。但必须提醒：强度暗点本身并不严格等于涡旋。有些旁瓣之间的极小值也可能表现为暗点，但它们不是相位奇点。
我在处理这一步时，会结合仿真图做对比。把实验强度图和仿真密度图叠加，看暗点位置是否一一对应，间距是否匹配。如果位置偏差超过一个光斑直径，基本可以判断光路中有像散或者SLM相位校正不准确。
4.2 干涉法观测相位结构：叉形条纹的解读
验证涡旋最可靠的方法是干涉。标准做法是用马赫-曾德尔干涉仪，让信号光束和参考平面波干涉。涡旋点位置会出现典型的叉形条纹，叉的朝向和数量直接对应拓扑荷的正负和大小。
如果你觉得搭马赫-曾德尔干涉仪太复杂，还有一个非常实用的替代方案：把全息图的零级光和+1级光干涉。因为零级光可以近似看作平面波，+1级光就是我们想要的涡旋阵列，两者在像面上自然干涉，同样能看到叉形条纹。这个方案的优点是无需额外增加光路，光路稳定性好，对环境振动不敏感。
在干涉图上数叉形条纹时，要沿着涡旋点所在位置横截过去数分叉数。单个拓扑荷为+1的涡旋对应一个叉形分裂；如果看到两个叉叠加，那可能是有两个涡旋靠得非常近，或者拓扑荷为2的模式混入。判断拓扑荷的更严谨方法是在相位图上沿闭合路径积分相位梯度，总变化量为2πl，l就是拓扑荷。这个积分可以用MATLAB或Python直接对重构的相位数据做。
4.3 用自参考干涉提升鲁棒性
马赫-曾德尔干涉仪对精密调节的要求比较高，尤其在实验室没有气浮台的情况下，空气扰动会让条纹不断漂移。如果你和我一样不愿意在干涉仪上耗时间，强烈推荐自参考干涉。
自参考干涉的原理是将信号光束分出一小部分作为参考，其余部分通过一个涡旋相位板（或再经过一次SLM）后重新合束干涉。由于两束光走的是几乎相同的光路，共模噪声被极大抑制，即使桌面有轻微振动，条纹依然能保持稳定。我实测下来，自参考干涉的条纹对比度比马赫-曾德尔好很多，而且调起来只需要一到两个反射镜，十分钟基本能对准。
4.4 从叉形条纹推算涡旋阵列的全局拓扑结构
涡旋阵列的观测最后要落到一个“全局拓扑”的判断上。我可以用前面仿真时预测的一个判据来做验证：阵列中相邻涡旋之间的相位差应该是均匀的。
实际操作中我会这样做：将干涉图导入图像分析软件，提取叉形条纹的坐标和方向，然后把所有涡旋点坐标拟合到一个椭圆上，检查相邻涡旋方向是否连续变化。如果方向沿椭圆弧逐渐旋转，说明阵列整体的相位梯度是光滑的，这正是Ince-Gaussian叠加模式的典型表现。如果方向出现跳变，可能是叠加相位差设置不对，或者有杂散模式混入。
5. 常见问题与排查技巧实录
5.1 SLM衍射效率低、杂散光明显
这个问题几乎每个人都逃不掉。排查顺序先把偏振片转到最佳位置，观察衍射亮度变化；再检查全息图的灰度范围是否覆盖了整个0-255区间，如果相位编码被压缩在很小的灰度范围内，衍射效率会大幅下降。我用的是8μm像素的SLM，最佳工作灰度范围大约在60~200，超过这个范围相位响应容易饱和。
如果零级光始终去不掉，大概率是光栅周期太密，导致+1级和零级分离不足。把光栅周期从6像素改为10像素，分离效果会显著改善。
5.2 阵列数目对不上、涡旋位置偏移
这是最让人头疼的问题。第一次做的时候，我预期6个涡旋，实际只看到4个，仔细看干涉图，发现有两个涡旋点被杂散光“淹没”了。处理方法是降低激光功率或者缩减曝光时间，让暗点更明显；同时也检查全息图的振幅阈值是否设置得太高，把弱涡旋的相位信息抹掉了。
位置偏移则往往来自入射角和SLM面板的倾斜像散。如果实验台上SLM的入射角固定无法调整，可以尝试在仿真阶段把椭圆参数ε适当调大或调小，补偿实验系统的像散偏差。
5.3 干涉条纹抖动、对比度差
抖动主要是气流和振动。实验室里不要开空调直吹光路，可以在光路附近加一个挡风罩；曝光时间压到10ms以内，可以“冻结”大部分高频抖动。对比度差则优先检查参考光和信号光的强度比，理想情况是参考光强度略低于信号光，这样条纹的调制度最高。
5.4 相机拍摄的几个细节
相机设置上，我习惯关闭自动增益，手动设曝光时间。CCD的暗噪声在长曝光下会明显增大，影响叉形条纹的清晰度，所以能短曝光就短曝光；动态范围有限时，优先保证涡旋暗点不饱和，旁瓣的强弱可以后续通过图像处理增强。拍摄时尽量使用RAW格式，避免有损压缩在涡旋位置附近产生伪像。
最后说一个我自己的习惯：实验记录里一定要把SLM的灰度查找表版本、入射角和透镜焦距都写清楚，不然隔两周回来想复现，会发现参数对不上了。做这类光学观测课题，细节的重复性往往比某一次的“漂亮结果”重要得多。
