Ince-Gaussian光束生成涡旋阵列：VirtualLab Fusion仿真全解析

之前一直在VirtualLab Fusion里折腾结构光束仿真总想着用现成的拉盖尔-高斯或厄米-高斯模式拼出涡旋阵列结果不是对称性不理想就是阵列排布太“正”调参调到怀疑人生。后来换到Ince-Gaussian这一类解系才意识到自己一直没找到合适的“基底”这类模式本身就生长在椭圆对称的坐标系里高阶项自带多个相位奇点只要参数选对涡旋阵列几乎是“长”出来的而不是“拼”出来的。这篇就从Ince-Gaussian光束的物理源头开始结合我在VirtualLab Fusion里的完整建模、参数设计、仿真观测和踩坑记录把“如何用Ince-Gaussian产生涡旋阵列激光束并观测”这件事讲透。适合正在做结构光场仿真、光学涡旋、光镊或光通信方向的朋友参考也欢迎刚从激光器方向转过来的同学对照着操作。1. 这次仿真做什么Ince-Gaussian与涡旋阵列的对应关系1.1 Ince-Gaussian光束在光学仿真里的定位光学里最常用的傍轴光束解族有两套直角坐标系下的厄米-高斯Hermite-Gaussian, HG和极坐标系下的拉盖尔-高斯Laguerre-Gaussian, LG。HG光束在横截面上呈现矩形的强度网格LG光束则带着同心圆环和螺旋相位当轨道角动量不为零时光束中心会出现一个相位奇点。问题是这两套解都只适合对称性“规则”的边界条件如果你想要的是一组按椭圆轨迹排列、相对独立又保持整体相干的多涡旋阵列直接拿HG或LG去叠加往往需要在相位分布上做很不自然的拼接。Ince-Gaussian简称IG光束正好补上这块拼图。它是在椭圆坐标系下求解傍轴波动方程得到的一族完整正交解用两个整数指标p和m再加一个连续椭圆参数ε来描述。ε0时IG模式退化为HGε趋于无穷大时趋近于LG所以它实际上是一个连接HG和LG的连续谱系。当你选用某个非零ε时模式的横向强度分布会是椭圆对称的并且会在椭圆环带上产生规律分布的暗核这些暗核就是光学涡旋的所在位置也就是题目里说的“涡旋阵列激光束”的直接来源。在仿真软件里IG模的价值不只是多了一种可用的数学形式。它天然地给了你两个额外的自由度椭圆率和两个模式指标。想调节涡旋个数改m想调节径向分层结构改p想改变涡旋在椭圆环上的疏密和整体偏心程度改ε。这种“参数和物理现象几乎一一对应”的特性在实际仿真中非常省事不需要来回在优化算法里反复试探。1.2 涡旋阵列是怎么从IG模里“长”出来的要理解涡旋阵列的来源先要明确“涡旋”在光场里到底是什么。一个光学涡旋对应的是复振幅相位分布中的一个奇点在奇点处光强为零相位围绕该点绕一圈后累计变化2π的整数倍这个整数就是拓扑荷。强度图上的“暗核”和相位图上的“螺旋结构”是同一个物理对象的两面。IG模在椭圆坐标下的解具有偶模IG^e和奇模IG^o之分。以某个偶模为例当m取较大值时它的横向相位分布在椭圆环上一圈会经历多次符号交替这种交替自然地在椭圆轨迹上形成多个相位奇点。每个奇点附近光场的等相位线是螺旋状结构从强度图上就能看到一个个暗核很有规律地排在椭圆环上看起来就是一个“阵列”。如果把p也调大径向方向上还会出现多个这样的环涡旋阵列就从单环变成多层结构。这个特性对照实验物理来说非常重要。过去想产生涡旋阵列常见的办法是用空间光调制器SLM叠加多个LG模式再靠全息图“拧”出一组涡旋。但这种叠加以强度干涉为前提不同模式之间的相对相位极难控制实际得到的结果往往包含串扰和杂散强度。IG模作为一个单一解不需要额外叠加它自身的场分布就携带了多个涡旋从源头上避开了多模干涉的稳定性问题。这也是我后来在VirtualLab Fusion里坚持用IG模做仿真的根本原因——它不是“逼近”涡旋阵列而是“直接存在”于模式本身。2. VirtualLab Fusion建模思路与核心参数设计2.1 建模路径怎么选自定义光源还是外部场导入在VirtualLab Fusion里要实现一种非内置的激光模式通常有两条路一个是使用软件自带的“可编程光源”Programmable Source通过脚本定义光源平面上的复振幅分布另一个是在外部数学软件里预先算好场分布再把电场数据文件导入到光源元件中。两条路我都走过优缺点很明显。可编程光源胜在“一次建模、参数联动”。在VirtualLab Fusion的光路图里插入一个可编程光源把IG模的复振幅公式写进脚本之后p、m、ε这些参数可以直接挂在光源属性面板上。改一个参数就能马上重新计算并观察到结果变化特别适合做参数扫描。缺点是对编程有一定的要求而且不同版本的软件接口命名有差异网上能直接查到的例程并不算多。外部场导入则比较“傻瓜”。你在MATLAB或Python里把电场算好按软件要求的网格和格式导出成电场文件再以“导入场”的方式加载到光路中。它的好处是计算IG模时可以借用更完善的数学工具保证Ince多项式的数值精度坏处是网格尺寸、采样点数以及相位包裹方式都得和VirtualLab Fusion内部约定一致否则导入后经常出现图像错位或振幅异常。就我个人的经验第一次做IG模仿真时建议先走外部导入路线。原因很简单IG模的解析表达式写起来不轻松尤其是Ince多项式的数值求值与其在软件自带脚本环境里反复调试不如先在自己熟悉的计算环境里把场分布算准确再交给VirtualLab Fusion去处理传播和衍射。把“源生成”和“传播计算”拆开出问题时定位也方便。等你确认了整套参数规律再迁移到可编程光源里做批量扫描。2.2 决定阵列形态的关键参数p、m、ε、束腰与波长很多初学者拿到IG模第一反应是问“参数应该设多大”。这个问题的前提是先搞清楚每个参数在光场里的物理作用。我整理了一套自己常用的参考关系不一定精确但足够让思路清晰参数作用典型取值范围调整趋势p径向结构指数0~5p越大径向层数越多涡旋环从单层变为多层m角向结构指数2~10m越大椭圆环带上的涡旋个数越多ε椭圆离心率/椭圆参数0.5~5ε越大结构越接近LG模阵列越向圆环汇聚ε越小越接近HG模阵列趋近矩形排列λ波长532 nm、632.8 nm、1064 nm影响衍射角度和传播演化速度不改变横向阵列拓扑排列w₀束腰半径0.5 mm~2 mm决定基模高斯包络的尺度影响暗核间距和观察窗口大小举个例子如果我想在束腰平面看到4个涡旋围绕中心等距排列可以先选m4、p0。此时IG模的相位奇点分布在椭圆环上数量和m直接相关。如果想看到两圈涡旋半径不同的“双环阵列”就把p设为1或2径向方向出现第二层结构涡旋总数会明显增加。这里有个容易犯的误区p并不直接等于“环数”但它控制径向函数的节点数径向节点越多光强分布中暗环和暗核的组合就越丰富。ε的选择则需要结合目标光斑形状。ε0时模式退化成HG横向结构接近矩形亮块和暗线ε取中等数值时涡旋排列在椭圆环带上呈现“椭圆项链”状的阵列ε非常大时模式退化成LG类结构涡旋围成近似标准圆环。如果你希望阵列在某些方向上拉开间距就往中等ε方向调如果希望阵列“抱团”成圆就适当加大ε。波长和束腰这两个参数容易被忽略但它们直接影响仿真中观察窗口的设置。在VirtualLab Fusion里探测器的物理尺寸必须覆盖光束的主要能量范围束腰太小而窗口太大会浪费大量采样点束腰太大而窗口太小边界会切掉光场的高斯尾翼造成严重的边缘衍射伪影。我一般先按束腰的3到5倍来设定观察窗口半宽再根据实际强度图微调。关于仿真网格VirtualLab Fusion的采样设置在光源平面和探测器平面都需要独立确认。网格点数太少暗核位置会“糊”成一片网格点数太多计算时间成倍上涨尤其是在做多距离或参数扫描时性价比骤降。我通常以512×512起步先观察阵列形态确认没问题后再把关键参数附近的扫描任务提升到1024×1024这样能兼顾效率和精度。3. 完整实操流程从光源定义到探测器出结果3.1 第一步构造IG模复振幅场不管选哪条建模路径第一个要解决的就是“怎么把IG模的场写出来”。IG模的电场表达式可以写成高斯包络乘以Ince多项式及其连带函数的组合在平面z0上偶模和奇模的形式可以表示为椭圆坐标(ξ,η)下的分离解。椭圆坐标用半焦距f与直角坐标关联其中f正比于ε和束腰w₀的乘积。η方向上的Ince方程解出角向函数ξ方向解出径向函数两项相乘后再乘上exp[-(x²y²)/w₀²]的基模高斯包络就得到了完整的复振幅。实际编码时不需要重新推导整套公式Mathematica和MATLAB的符号/数值工具箱都有办法计算低阶Ince多项式。我自己习惯先在Python里写一个“预检脚本”把场强分布和相位分布画出来确认参数组合对应的阵列形态再决定是否进入VirtualLab Fusion走完整仿真。下面是一个只用于预检的概念框架import numpy as np # 参数设置预检用 wavelength 632.8e-9 # 氦氖红光 w0 1.0e-3 # 束腰半径 p 1 m 4 epsilon 2.0 N 512 window 5e-3 # 观察窗口半宽 # 生成直角坐标网格 x np.linspace(-window, window, N) X, Y np.meshgrid(x, x) R2 X**2 Y**2 # 注意这里只是示意流程真正的IG模需要替换为Ince多项式数值解 # 下面用LG模作为占位用于提前观察涡旋个数和对称性 E ( (X 1j*Y)/w0 )**m * np.exp(-R2/w0**2) # 后续可输出强度图和相位图做预判这段代码不是完整的IG模实现只是给你一个感受在进入VirtualLab Fusion之前先用脚本把光斑形态大致看一遍能避免反复在软件里试错。真正在软件里建模时如果走外部导入路线就把最终计算好的复数电场保存成软件可读的网格文件如果走可编程光源就把解析表达式按软件语法翻译成脚本。一个非常需要注意的细节是相位奇点的数值稳定性。在奇点附近Ince多项式的符号跳变很容易造成数值尖峰。预检脚本里一旦发现某个暗核周围的相位分布出现非螺旋的杂乱条纹就要考虑是网格分辨率不够还是Ince多项式计算精度不够提前修正否则导入VirtualLab Fusion后会更难排查。3.2 第二步搭一条最简单的传输链路在VirtualLab Fusion中新建光路图Light Path Diagram后整条链路可以保持极简自定义IG光源、自由空间传播段或理想透镜、探测器。很多做工程仿真的朋友习惯一上来就堆一堆透镜和光阑其实对于“观测IG模涡旋阵列”这个目标最优先应该做的是在原生状态下确认模式本身的长相。加再多元件之前先看“裸奔”的结果这样一旦后面加了元件后阵列变形你就知道是元件引入的像差而不是模式本身的问题。自由空间传播距离的选择也很有讲究。在束腰平面上IG模的强度分布最“标准”暗核数量、位置都符合理论预期适合用来做模式验证。传播一段距离后由于不同频率分量衍射速度不同阵列会发生整体旋转和径向扩展暗核间距会变化这适合用来观察涡旋的演化行为。我一般会先设z0观测一次再设z1m和z2m各观测一次形成一个纵向演化序列。如果想让涡旋阵列在有限范围内保持清晰而不至于衍射散得太开可以加一个焦距合适的透镜。但这里有一个容易翻车的点透镜的孔径必须足够大至少是光束束腰的3倍以上。否则透镜边框会硬性截断IG模的场分布产生强烈的边缘衍射环这些环会叠加到阵列强度图上造成“暗核之间出现不明亮纹”的假象。3.3 第三步探测器设置与仿真运行在VirtualLab Fusion里添加探测器时我建议同时放置两类观测强度检测器和复场检测器。强度检测器直接输出光强分布用于看暗核阵列排布复场检测器保留振幅和相位数据用于后续的相位奇点确认和拓扑荷计算。只看强度不看相位是初学时最常犯的错误——光强图里的暗斑不一定都是涡旋也可能是振幅为零但不是相位奇点的点必须用相位分布来确认。窗口尺寸方面按照前面说的以束腰的3到5倍为参考。采样点建议至少512×512检测IG模p、m都不大时512已经能看清暗核但如果m超过6涡旋间距变小512网格容易出现相邻暗核连成暗线的情况可以升到1024。仿真运行完成后保存强度图和相位包裹图再通过软件的后处理功能做相位解包裹unwrapping这个过程能在暗核附近更清楚地显示螺旋相位变化。运行过程中如果遇到报错最常见的原因是光源平面和探测器平面的网格不匹配或者传播算法选了不适合大发散角的类型。VirtualLab Fusion的传播算法有好几种常见的有平面波分解PWD和瑞利-索末菲积分。对于准直性好、发散角小的情况PWD又快又稳对于发散很快或包含强相位梯度的场瑞利-索末菲积分更鲁棒。IG模的高阶项带有剧烈相位变化我习惯先用瑞利-索末菲类型做短距离传播确认结果正常后再换算法对比。4. 结果怎么“观测”强度、相位与涡旋阵列的对照4.1 强度图里的暗核对应什么仿真跑完最先看到的必然是强度分布。一组漂亮的IG模涡旋阵列在强度图上表现为若干个暗核沿着一条或几条椭圆环线分布。暗核是光强为零或接近零的点周围亮度呈“凹陷”状看起来像是椭圆项链上串着的一颗颗黑珠。这里要提醒一句不是所有暗斑都是涡旋。高斯光束在透镜焦点处形成的中心暗斑、两个模式干涉形成的暗线、偏振相关的暗带这些都不是涡旋。判断一个暗核是否为光学涡旋的金标准是检查它周围360°的相位是否完成整数倍的2π累积变化。所以强度图只能提供线索结论要交给相位图来下。另外还要注意暗核的“成对性”。在自由空间中光学涡旋的拓扑荷必须整体守恒。一个IG模的涡旋阵列里1拓扑荷和-1拓扑荷的涡旋通常配对出现整体角动量为零或呈现某种对称分布。如果仿真结果里所有暗核都是同一个旋转方向那多半是参数选择错误或者数值计算出了问题。4.2 用相位判据确认涡旋阵列的质量在相位图上涡旋表现为相位包裹条纹的“分叉”结构。沿某个闭合路径绕暗核一圈相位从0到2π跳变的次数就对应拓扑荷的大小和符号。VirtualLab Fusion里显示相位图时默认是包裹模式相位值被折叠在-π到π之间这时候螺旋条纹上会看到清晰的分叉线。每一条分叉线就是一个相位奇点所在的位置。把强度暗核位置和相位分叉位置叠在一张图里是验证IG模阵列质量的常用做法。理想的IG模仿真结果是每个暗核中心都恰好对应一个相位分叉点没有多余的暗核或多余的分叉。如果出现“暗核有相位分叉无”的情况说明那个暗点可能是振幅零点但不是相位奇点如果出现“相位分叉有暗核不明显”的情况多半是网格采样不够需要提高空间分辨率。还有一种常见问题是相位解包裹后整体带有缓慢变化的倾斜项。这通常是因为光源定义时场分布中心没有对准仿真坐标系原点导致传播后整个波前携带了一个线性相位斜坡。解决办法是回到光源定义那里重新检查场数据的中心坐标和网格原点是否对齐。4.3 扫参看演化从单个涡旋到阵列过渡做这种仿真最有趣的部分不是只看一张强度图而是观察“阵列是如何成形的”。我强烈建议在VirtualLab Fusion里做一次ε或m的参数扫描例如固定p1、m4把ε从0.5扫描到5。ε较小时强度图呈现接近HG的矩形网格暗核分布比较“散”还不成环随着ε增大暗核逐渐向椭圆环上聚拢形成清晰的项链状阵列ε再继续增大椭圆趋于圆涡旋阵列变成圆环结构接近LG模的特征。这个过程本质上是在看在同一个连续解族中涡旋如何从“准矩形排布”过渡到“椭圆排列”再到“圆形排布”。它能直观回答一个前期经常困惑的问题为什么IG模能产生“阵列”而HG和LG不能。因为IG模正好处在直角坐标和极坐标两个体系之间它允许涡旋在保持“共轴相干”的前提下沿椭圆轨迹实现空间有序排列而在两个极限端点要么排成规则的网格暗线HG要么集中在个别环上失去明显的多环阵列感LG。扫描结束后可以把不同ε下的强度图和相位图导出做成横向对比。我自己常用做法是把ε0.5、1.5、3、6四组结果并排放在一起能够很直观地向同事展示IG模、HG模、LG模三者之间的连续性也方便确认实际应用时选哪个ε范围能兼顾阵列数量和美观度。5. 仿真过程中的高频坑与排除方法5.1 强度图有暗斑但相位图对不上这个现象我遇到过不止一次。某个暗核在强度图上非常明显但在相位图上看不到对应的螺旋分叉或者分叉位置偏了半个格点。大概率是数值采样的问题。IG模的相位奇点在物理上是一个数学点但在离散网格上它最多被解析到单个网格的尺寸。如果奇点恰好落在网格边界上强度暗斑会被拉成细长暗条相位分叉则被“夹”在两个网格之间导致看起来“对不上”。解决方法是加密网格或者对场数据做一次亚像素插值后再分析。VirtualLab Fusion里可以提升探测器的采样倍数也可以把光源平面的网格先校准到更均匀的密度。还有一种可能是场数据本身存在明显的相位缠绕错误。在导入外部电场文件时如果相位被截断到-π到π之后没有正确保持连续性解包裹算法就会出错产生虚假的分叉线。这种情况下回到源数据检查相位矩阵的相邻差分是否超过π的阈值如果是说明包裹之前的数据就有问题。5.2 阵列在传播后变形严重涡旋阵列在传播一段距离后发生形变是正常的因为IG模不是传播不变的横向结构会随z演化。但如果形变表现为暗核数量变化、阵列对称性被明显破坏、或者出现理论之外的暗核就要排查两个原因。第一个原因是光源平面的观察窗口不够大。初始场的高斯包络虽然衰减但椭圆坐标下Ince多项式在远离原点的区域仍然有非零结构窗口太小等于把一部分真实场截断了截断处会产生额外的衍射边缘场传播后这些边缘场会干扰阵列甚至形成新的暗核。第二个原因是传播算法选择不当。用PWD计算高度发散场时平面波展开的有限窗口会造成周期性混叠让远场分布出现重复像尤其在窗口尺寸设置不合理时阵列看起来就像“多了一套”。我排查这类问题的方法是先加宽初始窗口和探测器窗口对比结果是否稳定然后把传播距离缩短到0.1m量级观察暗核数量和位置是否与束腰处一致最后如果问题依旧切换到瑞利-索末菲积分再跑一次。这三步基本能覆盖80%的形变异常来源。5.3 仿真耗时太长怎么办高阶IG模比如p3、m8加上1024×1024网格再配合多距离扫描仿真时间会非常感人。优化方向有三个。降低非关键平面的采样密度。束腰处做高分辨率精确采样传播中段的探测器可以先用256×256网格看看趋势最后平面再回到1024×1024出图。VirtualLab Fusion允许不同探测器使用不同采样率这比全局统一网格要省时得多。关闭不必要的多线程辅助输出。光路图里如果有多个可视化监视器、实时动画或伯德图之类的辅助显示计算过程中都会持续刷新拖慢主仿真。先关掉只保留最终结果的探测器运算速度往往能提升一倍以上。对参数扫描任务优先用“参数扫描”模式而不是手动一版一版跑。VirtualLab Fusion的参数运行Parameter Run会为不同参数组合复用公共传播段数据比单独跑多个任务效率高很多。我在做ε扫描时实测同样跑10个参数点参数运行模式大约比手动逐个跑快50%到一倍。还有一个“偷懒”技巧先用较低阶模式p0或m较小确认光路、算法和探测器设置都正确再切换到目标高阶模式做正式计算。这样能在调试阶段避开高阶复杂场带来的长耗时把排查范围压缩到光学设置本身。6. 从虚拟到实验后续可扩展的方向6.1 用空间光调制器复现IG涡旋阵列仿真跑通之后一个自然的问题是这套涡旋阵列能不能真在实验台上做出来答案是肯定的主流方案是用空间光调制器SLM加载相位全息图。具体思路是把VirtualLab Fusion里计算得到的IG模复振幅分布编码成纯相位全息图或复杂调制全息图加载到SLM上SLM用激光照射后在傅里叶平面上通过空间滤波取出目标级次的衍射光就可以获得对应的IG模涡旋阵列。这里最关键的一点是仿真时使用的网格尺寸、采样率以及坐标原点决定了全息图编码的分辨率和衍射效率。如果仿真阶段就把网格设置得胡乱到实验阶段会发现SLM上加载的全息图出现明显的像素化噪声。我建议在仿真阶段就按实验SLM的像素数和像素尺寸来设计探测器窗口。比如常见的SLM像素尺寸是8μm左右像素阵列为1920×1080你在VirtualLab Fusion里就尽量按这个比例设置观察窗口这样从仿真场到全息图的映射关系是1:1的避免后续重新采样带来的信息损失。这个细节听起来琐碎但直接影响最终涡旋阵列的对比度和模式纯度。6.2 下一步能做的研究方向IG模涡旋阵列的仿真平台搭好之后扩展方向相当丰富。一个是把IG模和其他功能元件结合例如通过柱透镜或变形镜系统研究涡旋阵列在像散条件下的演化与拓扑荷转移。一个是把阵列作为“多光镊”的候选光场用仿真评估不同ε和m参数下每个涡旋的强度梯度力分布为实验提供参数依据。还有一个方向是利用p和m的自由度做空间模式复用在光通信中把不同IG模当作独立信道用仿真预判模式串扰和大气湍流影响这类选题无论在学术发表还是工程项目里都有实际价值。我个人在实际操作中的一个体会是IG模的仿真难点不在“算出涡旋阵列”而在“稳定复现一套自己可控的阵列”p、m、ε和w₀这组参数的关系值得在最初建模阶段就形成一条清晰的调参逻辑线。不同版本的VirtualLab Fusion在光源定义和探测器设置上略有差异但只要弄清楚了场分布如何定义、传播算法如何选用、相位奇点如何验证换版本也好、换实验平台也好核心迁移成本都很低。希望这次记录能帮你在自己的结构光场仿真里少走一些弯路。