在结构光场这个圈子里Ince-Gaussian 光束这几年热度一直不低。它不像 Hermite-Gaussian 那样方方正正也不像 Laguerre-Gaussian 那样自带圆周运动但它偏偏夹在两者之间成了研究涡旋阵列的一个特殊入口。配合 VirtualLab Fusion 这类物理光学仿真平台我们不用搭实验台就能把 Ince-Gaussian 光束产生的涡旋阵列预处理、传播过程和后焦面观测一次性跑通。这篇文章就基于我用 VirtualLab Fusion 做 Ince-Gaussian 光束涡旋阵列观测的完整经历来写适合正在做结构光场仿真、涡旋光通信或者微纳光学设计的同行参考。搞物理光学的人都知道所谓“涡旋阵列”简单说就是一束光里同时带着多个相位奇点每个奇点对应一个轨道角动量通道。这种光束在光镊、光通信、并行激光加工里都有用。而 Ince-Gaussian 光束的特殊之处在于它在一个椭圆坐标系里长出来高阶模式天然能往外延伸出多个相位分支叠加之后能非常自然地形成排列规整的涡旋点阵。VirtualLab Fusion 的核心价值是“场追踪”它不把光当成理想光线而是把整个复振幅场从源面一路推到观测面相位、偏振、角谱全保留所以非常适合用来研究这种拓扑结构复杂的场。1. 为什么 Ince-Gaussian 光束适合用来搭建涡旋阵列1.1 Ince-Gaussian 的出身决定了它与涡旋的关系Ince-Gaussian 光束是在椭圆坐标系下求解傍轴波动方程得到的正交解。椭圆坐标里有径向坐标 ξ 和角向坐标 η随之而来的解分奇模和偶模两类模式阶数用 p 和 m 标记。这里还有一个关键参数椭圆度 ε。正是这个 ε 决定了 Ince-Gaussian 到底更接近直角坐标系下的 Hermite-Gaussian还是更接近极坐标系下的 Laguerre-Gaussian。当 ε 趋近于 0 的时候Ince-Gaussian 会退化成 Laguerre-Gaussian所以很多携带涡旋的 LG 模式其实可以看作 Ince-Gaussian 的一种特例。当 ε 变得很大的时候它又会慢慢朝 Hermite-Gaussian 靠拢。这就给了我们一个很有意思的调节旋钮在保持模式阶数不变的情况下通过改变 ε可以让相位分布从“规则螺旋”连续过渡到“直角节点”。这种过渡不是简单的缝隙而是中间态本身就带有多个局域相位奇点。多个 Ince 模式叠加时这些奇点会重新排列组合在横截面上形成整洁的阵列结构。举个例子当 p 和 m 都比较高的时候Ince-Gaussian 的相位分布会在椭圆轨迹上分裂出多个暗核每个暗核附近都有 2π 相位旋绕。这些暗核的位置跟随椭圆坐标的取向和 ε 变化这就给了我们控制阵列旋转角度和疏密的可能。相比直接用多个 LG 光束拼一个晶格Ince-Gaussian 方案在同一个本征模式下就能呈现出多个奇点不用额外考虑光束之间相干叠加的相位匹配问题物理上更干净。1.2 从“单涡旋”到“涡旋阵列”的位置可控性以前大家做涡旋阵列常见做法是取多个 LG 模式把它们按一定几何位置排开再做相干叠加。这样做的问题在于每个 LG 光束需要单独对准阵列的整体对称性高度依赖初始位置参数。Ince-Gaussian 的模式结构自带椭圆坐标族高阶模式的奇点分布自然会沿着某个椭圆轨道展开。我们人为调节 ε 和旋转角就能让奇点分布在圆环、椭圆环甚至近似直线的轨道上。在 VirtualLab Fusion 里这种位置可控性可以直接通过修改初始场数据实现。比如我在脚本里生成一个 p2、m2 的 Ince-Gaussian 偶模再把对应的奇模叠加进去就会形成成对出现的涡旋每对涡旋的拓扑荷符号相反整体阵列呈现中心对称分布。这个效果如果用空间光调制器实验做你得先经过一轮复杂的全息图计算再担心调制器的像素分辨率够不够。而在 VirtualLab Fusion 里我们只需要把复振幅场导入光源然后让物理光学传播算子把场推到目标面观测强度分布和相位分布就能直接确认阵列形态。所以“观测”这个词在仿真里的意义是可以直接打开相位图看暗核周围的等相位线旋转方向。这比实验干涉法测量要直观太多。2. VirtualLab Fusion 里的两种建模路线我在实际做这个课题时发现进入 VirtualLab Fusion 有两条路线一条是把 Ince-Gaussian 的完整复振幅场直接灌进光源另一条是用相位掩模近似。这两条路线各有适用场景下面分开说。2.1 精确路线把 IG 场分布直接灌进光源这条路线最贴近“真正产生 Ince-Gaussian 光束”的物理定义。我们先在外部脚本里把 Ince-Gaussian 电场分布算出来保存为一个二维复数数组再用 VirtualLab Fusion 的场光源导入功能加载它。在 Python 或 Mathematica 里算 Ince 多项式并不复杂。Ince 方程的形式可以写成椭圆坐标下的二阶微分方程解 InceC 和 InceS 可以通过傅里叶级数和特征值方法求。很多开源的 Mathematica 程序包里都有现成的 Ince 函数实现我之前也照着 Bandres 和 Gutiérrez-Vega 在 2004 年那篇经典论文里的递推公式写过一段脚本几百毫秒就能算出 512×512 的场。计算时要特别注意坐标系映射因为场是定义在笛卡尔网格上的而 Ince 多项式是在椭圆坐标里定义的。代码大致是这样的逻辑import numpy as np N 512 x np.linspace(-5e-3, 5e-3, N) y np.linspace(-5e-3, 5e-3, N) X, Y np.meshgrid(x, y) # 椭圆坐标转换 f0 1e-3 # 半焦距由 ε 和束腰决定 u np.sqrt((X**2 Y**2 f0**2 np.sqrt((X**2 Y**2 f0**2)**2 - 4*f0**2*X**2)) / (2*f0**2)) v np.arctan2(Y, X) # 近似角度坐标实际需要完整公式 # 计算 Ince 多项式 C_p^m 和 S_p^m C compute_ince_C(p, m, epsilon, np.cos(v)) S compute_ince_S(p, m, epsilon, np.cosh(u)) field C * S * np.exp(-(X**2 Y**2) / w0**2) # 保存振幅和相位VirtualLab 导入用 np.savetxt(ig_field_amp.dat, np.abs(field)) np.savetxt(ig_field_phase.dat, np.angle(field))这里我特意强调ξ 方向和 η 方向都要使用真正的椭圆坐标转换公式如果直接用极坐标近似高阶模式会有明显变形。保存数据时建议同时输出振幅和相位两个文件或者直接输出复数实部和虚部两个通道这样导入后 VirtualLab Fusion 能重建完整的复振幅场。然后在 VirtualLab Fusion 里添加“自定义光源”把振幅文件和相位文件分别加载进去再设置波长、束腰和能量。这里有个容易踩的坑VirtualLab 内部坐标单位可能默认是毫米脚本计算时如果用了米导入后尺度全部错掉传播几毫米就散得不成样子。我建议在脚本阶段就把坐标单位统一成毫米波长单位用微米这样和 VirtualLab Fusion 的输入窗口保持一致。2.2 快速路线用相位掩模近似这个方法适合只想快速看涡旋阵列外观不想纠缠 Ince 多项式计算的场景。本质上就是做一个相位型衍射元件把目标涡旋阵列的相位分布画成灰度图再导入 VirtualLab Fusion 作为相位掩模。对于多个涡旋组成的阵列相位分布可以直接写成Φ(x, y) arg(Σ_j exp(i l_j θ_j) · exp(-((x-x_j)² (y-y_j)²) / w_j²))其中每个涡旋有自己的位置、拓扑荷 l_j 和束腰 w_j。这种相位掩模加载到 VirtualLab Fusion 里之后再用一个透镜做傅里叶变换焦平面上就会出现对应的涡旋阵列。这种方法非常快但它绕过了 Ince-Gaussian 本身的模式结构更接近“人工构造涡旋阵列”的路子。2.3 方案对比与适用场景对比项精确导入 IG 场相位掩模近似计算量中等需要解 Ince 多项式很小叠加螺旋相位即可物理一致性高保留 Ince 模式的本征结构低近似成立但非 Ince 模式阵列形态控制通过 ε、p、m 控制通过位置、拓扑荷控制适合场景研究 IG 本征模式特性快速原型验证、实验对照我个人的建议是如果课题核心是“Ince-Gaussian 光束的涡旋阵列观测”不要偷懒走精确路线。因为相位掩模近似出来的阵列虽然样子像但涡旋与涡旋之间缺少 IG 模式特有的椭圆相关性之后做模态分解或者纯度分析时会发现和理论预期差很多。如果是工程应用比如想快速估计涡旋阵列在透镜聚焦后的强度分布那相位掩模路线足够用。3. 实操复现从 IG 光源到涡旋阵列观测的完整过程下面这部分我按自己完整跑通一遍的顺序写包含了参数选择、界面操作流程和一些细节处理。版本上我用的是 VirtualLab Fusion 近年比较主流的版本界面大同小异。3.1 准备场数据第一步确定物理参数。我的仿真用的是 1550 nm 波长束腰 w0 取 2 mm计算窗口 8 mm × 8 mm采样点 512×512。选择 1550 nm 是因为涡旋光通信里这个波段很常用而且波长相对较长对采样要求友好一些。第二步确定 Ince-Gaussian 模式阶数。我测试过的组合里p2、m2 的偶模和奇模叠加产生的涡旋阵列最清晰。ε 取 2.0 时阵列呈现出椭圆环状排列中心是一个暗核周围均匀分布四个子涡旋。如果你希望阵列更接近圆环把 ε 调小到 0.5 左右如果希望阵列拉得更开就加大 ε。这个调节非常直观也是 Ince-Gaussian 相对 LG 模式的一大优势。第三步在 Python 里完成计算并保存文件。这一步有个细节就是在初始场里加入适当的相位曲率项。很多第一次做的人直接在束腰位置生成平面波前的 IG 场然后导入 VirtualLab 传播。这在数学上没错但实际物理上从激光器出来的模式通常带有球面波前。如果目标面是透镜前表面你还得在 IG 场里额外乘一个 exp(-i k r²/(2R))。不然后续传播后的阵列大小会和理论焦斑尺寸对不上。我一般会在脚本里预留 R 参数默认 R 很大代表束腰处的平面波前。3.2 VirtualLab Fusion 里搭建系统并配置传播打开 VirtualLab Fusion 后我习惯先新建一个空白系统然后在左侧的 System Explorer 里依次添加下面几类组件光源选择“自定义光源”在属性窗口里导入上一步保存的振幅文件和相位文件。如果数据是复数格式选择相应的实部虚部导入选项。透镜先用理想透镜做初步验证。焦距 f 取 200 mm口径设为 6 mm避免边缘截断影响。传播路径光源到透镜之间设为自由空间传播距离取 50 mm透镜到探测器之间的传播距离设为 200 mm也就是正好一个焦距这样探测器在焦平面上。探测器使用相机探测器窗口 8 mm × 8 mm像素数 512。配置传播算子时我通常选择“菲涅尔传播”或“瑞利-索末菲传播”。对于 200 mm 这种尺度菲涅尔传播已经足够精确计算速度快一半以上。但如果你研究的系统里有强聚焦透镜或者数值孔径超过 0.5那就必须换瑞利-索末菲算子否则高频分量会被低估阵列周围会出现不真实的光强振荡。VirtualLab Fusion 里还有一个“采样密度”设置项这个非常重要。我强烈建议在传播算子里开启自动采样让软件在每一步传播时检查局部相位梯度密度。如果软件版本不支持自动采样就手动把采样点数设成源场网格数的两倍以上。别舍不得计算资源涡旋阵列的相位结构非常细采样不够直接导致暗核消失。3.3 聚焦面上的探测与相位奇点分析仿真跑完后双击相机探测器查看结果。默认显示强度分布这时候应该能看到一个椭圆排布的暗点阵列。每个暗点就是涡旋核的位置中心光强接近零但周围强度高低不同。只看强度还不能确认拓扑荷必须切到相位视图。在 VirtualLab Fusion 的显示选项里把色彩图从强度切到相位或者直接输出相位矩阵到外部软件处理。展开相位后你会看到暗点周围一圈 0 到 2π 的连续变化这个就是涡旋存在的直接证据。拓扑荷的正负看相位旋转方向顺时针旋转是正拓扑荷逆时针是负拓扑荷。多个涡旋同时存在的时候可以通过数相位奇点的个数来确定阵列里涡旋的数量。这一步我自己的心得是把焦点位置稍微偏移几毫米再看一次相位分布非常有助于判断阵列的三维演化。涡旋在传播过程中会发生旋转特别是在 Ince 模式里阵列的整体取向会随着传播距离变化。VirtualLab 允许我们直接扫一系列传播距离这个功能比实验上逐面扫描要方便太多。我建议至少看三个位置透镜前 50 mm、焦点、焦点后 50 mm。这样能得到一个完整的“阵列演化”印象而不是只看到焦平面上的静态图案。4. 我踩过的坑采样、窗口和模式纯度做这类仿真的过程中我遇到了几个非常典型的问题写下来给后来者当参考。尤其是那些只看强度图就觉得“仿真成功了”的同学我的建议是永远要同时看相位图很多问题都藏在相位里。4.1 采样密度不够相位奇点会“凭空消失”这是我碰到最多的问题。当探测器的像素分辨率较低或者传播路径中的网格采样数不足时单个涡旋核周围可能只有一个像素或两个像素。这时候相位展开算法无法正确判断 2π 跳变视觉效果就是暗点还在但相位图变成一个平滑斜坡拓扑荷信息完全丢失。解决方法有两个层面。第一在光源导入时把初始场空间分辨率提高。512×512 对于单个涡旋足够但如果你的阵列有 9 个涡旋甚至更多一个涡旋核只占 5 个像素那就得改用 1024×1024。第二在传播算子里提高默认网格数。VirtualLab Fusion 的自动采样能缓解问题但自动采样有时会漏掉相位的快速变化区域所以最好手动检查一下传播链路里每个组件输出的网格尺寸。我自己的经验准则是每个涡旋核至少占据 5×5 个像素低于这个阈值宁可把计算窗口缩小或者像素数翻倍也不要贪图计算速度。毕竟仿真相对于实验来说最大的优势就是可以随便跑多等几分钟不是问题。4.2 孔径截断与窗口效应另一个常见问题是硬边孔径导致的衍射振荡。当 Ince-Gaussian 场的尺寸超过了透镜口径透镜边缘会对光束进行硬截断。这种截断会在焦平面产生一圈明显的同心环有时候环的强度比涡旋阵列本身还高差点让观测结果直接报废。解决思路有两个方向。一个是在透镜前插入一个软边光阑比如超高斯窗。另一个更干净的方法是反算束腰和透镜口径的比例确保 99% 以上的光功率都能通过透镜。以 w02 mm 为例透镜口径最好不小于 10 mm也就是束腰的五倍以上。如果因为系统限制只能使用小口径透镜那就在 VirtualLab 里把透镜换成“带孔径的理想透镜”并在孔径边缘设置一定的吸收系数。此外在初始场导入时也要注意计算窗口的边界截断。Ince-Gaussian 的尾部在椭圆坐标方向延伸较远如果计算窗口刚好切在强度较高区域会在边缘产生非物理的反射分量。我一般会在脚本里把窗口扩大到束腰的 4 到 5 倍并在边缘加一个余弦衰减窗。4.3 传播距离与 IG 模式的退化Ince-Gaussian 模式不是所有传播距离都能保持稳定形状的。当传播距离超过瑞利长度之后模式之间的相对相位差会导致场重新分布阵列里的涡旋位置会发生明显漂移。这在物理上是合理的但如果观测面设定不合理会把这种演化误认为仿真结果出错。解决办法是提前计算瑞利长度 z_R π w0²/λ。以 w02 mm、λ1550 nm 为例z_R 大约是 8.1 米远远大于实验室里 200 mm 的传播距离所以在短距离内模式保持得很好。如果你用 w00.5 mm 的细光束瑞利长度只有 500 mm 左右这时候选 200 mm 的传播距离基本还在安全范围内但再远就不好说了。我建议在做参数扫描之前先把瑞利长度算好并把它作为传播距离的上限参考。如果研究目的就是观察模式退化那就另说比如想看看 IG 阵列如何过渡成 LG 或 HG 模式这时候可以刻意使用长传播距离和较大 ε观测相位阵列如何逐渐失去椭圆特征。5. 参数、指标与实际体会5.1 常用参数速查表下面这个表是我在多次仿真中沉淀下来比较稳定的起始参数新手可以直接抄作业。参数推荐值备注波长1550 nm / 632.8 nm根据实验条件选束腰 w01~2 mm决定涡旋核尺寸计算窗口6~10 mm至少 3 倍束腰采样数512×512 ~ 1024×1024涡旋数量多时选大值椭圆度 ε0.5~2.0ε 越小越接近圆环阵列模式阶数p2, m2 适合入门高阶模式奇点更多透镜焦距200 mm 典型根据实验布局调整传播距离1~2 倍焦距观测焦点附近阵列探测器像素512×512跟采样数一致5.2 几个评估指标判断涡旋阵列生成质量光靠肉眼不够。我在做数据分析时通常会算四个指标。第一个是相位奇点计数。展开相位后统计相位从 -π 跳到 π 的孤立点数这个数应该和理论涡旋数量一致。如果数量偏少大概率是采样问题如果数量偏多可能是边界伪影。第二个是强度对比度。定义为中心暗核的最小强度与周围亮环平均强度的比值比值越低说明涡旋越“干净”。我通常要求这个比值低于 10%。第三个是阵列几何规则度。将所有暗核坐标提取出来拟合成椭圆看拟合残差。这个指标用于判断 Ince-Gaussian 是否真的保持了椭圆对称性。拟合残差如果超过一个像素就说明传播过程中可能有像差干扰。第四个是模式纯度。VirtualLab Fusion 可以把输出场分解到一组正交基上例如 LG 基或者 HG 基。如果输入的 Ince-Gaussian 纯度很高那么输出场分解到高阶 LG 模式的占比应该非常小。这种分解在实验上很难直接做但在 VirtualLab 里是内置功能配合“模式分析”工具很方便。如果输出场中出现了预期之外的 LG 分量多半是传播算子精度不够或者透镜产生了额外相位畸变。5.3 我自己的几个心得最后聊几句个人感受。VirtualLab Fusion 这个工具界面不算友好菜单层级也比较深但它的“场追踪”理念在结构光场研究里确实能省下大量验证时间。用实验手段测涡旋阵列相位你得搭马赫-曾德尔干涉仪还要特定空间光调制器校准一次不顺利可能就要折腾一周。在 VirtualLab 里只要把复振幅场推进去后焦面的强度和相位随手可得。不过仿真终归是仿真它不会替你解决“怎么生成真正的 Ince-Gaussian 光源”这个工程问题。我在做完这个课题后最大的体会是VirtualLab 最好的用法是作为实验前的快速验证工具用来确定参数窗口和诊断潜在问题。真正要产生稳定的 Ince-Gaussian 光束还是得回到空间光调制器或几何相位光学元件上但这时候你已经知道该用什么参数、该避开什么坑了。如果你也准备做类似方向我建议先从简单的两步开始先在 Python 脚本里算一个 p1、m1 的低阶 Ince-Gaussian导入 VirtualLab 看单涡旋形态确认单涡旋没问题后再升到 p2、m2 做阵列。不要一上来就往 16 个涡旋的复杂阵列使劲那个参数空间里的坑比你想的要多。