做结构光场的人近几年大概都会听到一个绕不开的名字——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, paritye): # 构造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, k1] 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, k1] 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)和椭圆参数ε。我建议一次只扫一个参数先固定p7, m3然后观察ε从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决定涡旋的个数和排列。我用仿真验证过很多种组合p7, m3时叠加场里能看到清晰的6个涡旋点沿一个椭圆环均匀分布p5, m2时是4个涡旋分布在椭圆内弧上间距更大一些。这里有个实操技巧仿真时不要只看强度图一定要把相位图也画出来。强度图上的暗点看起来可能只是强度极小值只有相位图才能确认那里确实存在拓扑荷为1的相位奇点。我通常的做法是计算环绕每个候选涡旋点的相位卷绕变化量如果一圈下来总变化量为±2π才认定它是真正的涡旋。2.4 仿真结果怎么解读当你拿到第一张像样的仿真强度图时你会发现一个有趣的现象涡旋点并不是像拉盖尔-高斯涡旋那样中心一个完美圆环而是椭圆形的暗核周围有对称分布的旁瓣。这是Ince-Gaussian坐标系本身带来的“遗传特征”不要慌张这正是Ince-Gaussian涡旋阵列区别于其他方案的指纹。我在仿真阶段还会顺便预测一个关键判据阵列中相邻两个涡旋之间的相位差应该是一致的。这样实验上如果看到干涉条纹出现等间距的叉形结构就说明涡旋阵列的相位拓扑是稳定有序的。3. 实验光路搭建与关键器件标定3.1 整体光路配置实验光路我采用的是一套标准的单光束SLM方案结构简单复现性好。整体链路如下激光器氦氖激光器波长632.8nm→ 半波片偏振片 → 扩束镜组20倍→ 反射式空间光调制器SLM → 透镜Lf 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πll就是拓扑荷。这个积分可以用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的灰度查找表版本、入射角和透镜焦距都写清楚不然隔两周回来想复现会发现参数对不上了。做这类光学观测课题细节的重复性往往比某一次的“漂亮结果”重要得多。
