1. 先把这个方向看明白PT对称不是玄学是“损耗换来的对称”PT对称系统这几个字一听像是量子力学理论里的概念但实际上这几年在光子学领域已经被玩出花来了。不管是波导、微盘还是光纤体系PT对称的核心思路都一样让系统在一个特殊的“增益-损耗配平”条件下表现出反直觉的光学行为其中最著名的就是奇异点Exceptional PointEP附近的本征态简并和自发对称性破缺。用一句稍微工程化的话说如果一个光学体系的复折射率满足实部空间偶对称、虚部空间奇对称那它就是PT对称的。我最早接触这个方向是在做集成光子学器件的时候。当时课题需要设计一个单向传输的波导结构翻文献发现PT对称是绕不开的一条路。与其拿耦合模方程从头写脚本求解我选择了COMSOL做全波仿真。原因很直接COMSOL的波动光学模块能真实处理波导、微盘、光纤这些带具体几何边界的结构模式场分布、辐射损耗、边界效应都能显式体现出来同时PT对称研究必然要反复扫描增益/损耗参数COMSOL的参数化扫描和特征值求解器配合得相当省心。这篇文章就把我在波导、微盘、光纤三个体系里的建模思路、实操细节、踩过的坑完整整理一遍适合正在做光子学方向的研究生以及想从纯理论计算转向全波仿真验证的工程师参考。2. 建模前必须想清楚的三件事2.1 物理场接口怎么选COMSOL里做光频段的电磁仿真最常见的两个入口是波动光学模块的“电磁波频域”接口和射频模块的同名接口。对PT对称研究来说绝大多数场景直接用波动光学模块就够了。这个接口求解的是频域矢量亥姆霍兹方程支持特征值研究、模式分析、频域扫描等多种研究类型。射频模块强在S参数、集总端口、电路耦合这些微波工程功能做光学波导和微盘反而用不上界面里还多了很多干扰选项。我自己的做法是结构简单就先用2D模型快速验证物理图像确认EP点位置之后再用3D模型做精细计算。2D模型计算量小一个数量级但有两个前提要注意。第一2D模型默认结构在第三维无限延伸平面波导、双芯光纤这种拉伸型结构完全适用第二如果你想观察的是沿传播方向的功率演化那就还是得老老实实建3D。我见过不少同学一上来直接抄教程里的3D双波导模型算一个大尺寸结构算了几个小时其实用2D截面先扫一遍参数完全能定位EP的大致位置。2.2 复数折射率的正负号约定这一小节我要放到最前面说因为这是我见过翻车率最高的地方。在COMSOL中材料折射率可以写成复数n n_r i * n_i其中虚部为正对应损耗介质虚部为负对应增益介质。如果你的本构关系里直接写介电常数那么先求平方根再定义虚部符号道理相同。很多文献受不同符号约定的影响有的写虚部为负是损耗你照搬过来就会把增益和损耗写反。最稳妥的办法不是去翻手册而是先做一个一维衰减验证建一段长度为L的均匀介质平板n 1.5 i*0.01在边界加端口激励沿传播方向画电场振幅如果振幅衰减说明你的约定没错再把虚部改成负号振幅增长说明增益也没写错。这一步五分钟就能完成但它能避免你在后续整套模型里浪费一个星期排查“为什么EP点永远找不到”。2.3 边界条件和PML的正确打开方式PT对称仿真里大部分体系是开放边界。波导有泄漏模微盘有辐射损耗光纤如果做弯曲也会有泄漏所以边界条件的选择直接影响计算结果。散射边界条件SBC实现简单但对大角度入射光的吸收效果一般完美匹配层PML把出射波“吃掉”得干净代价是增加计算域和网格量。我的建议是模式分析优先用散射边界因为模式场大多数局域在结构内部边界只是兜底选择SBC不会带来大误差而微盘这类高Q值回音壁结构一定要用PML否则算出来的Q值会受边界反射影响出现明显假象。PML的厚度至少要覆盖两个网格层PML内边界与结构之间要留出至少半个波长的空气区不然会激励出边界伪模。设置PML时还要注意把PML区域内网格各向异性拉开让径向尺寸放大吸收效果才够理想。3. 波导体系实操双通道定向耦合器找EP点3.1 几何建模和网格控制双通道波导是最经典的PT对称载体。两条平行矩形波导一条写增益折射率一条写损耗折射率背景是低折射率包层。我在COMSOL里常用的参数是波导宽度500nm高度220nm背景折射率1.45SiO2波导芯层折射率实部2.1两条波导的间距从200nm扫到1μm工作波长1550nm。增益损耗参数用一个全局参数g控制波导1的折射率写成2.1 - ig波导2的折射率写成2.1 ig。几何层面没有太多花活两个长方体加一个包层长方体在物理场中定义两个不同的“折射率”域即可。网格是这里的关键波导芯区和两条波导之间的间隙必须是细网格至少低于介质波长的十分之一。1550nm波长下硅波导内波长约738nm所以芯区网格控制在70nm左右比较稳背景可以放宽到200nm。使用“用户控制网格”而不是物理场控制网格因为物理场控制网格往往整体加密计算量不必要地暴涨。3.2 模式分析先看传播常数的实虚部研究类型选择“模式分析”面外波数设为变量。COMSOL求解得到的是复传播常数β类比有效折射率n_eff β / k0。PT对称的特征就在复β随g的变化曲线上当g小于临界值g_c时两个超模的有效折射率实部发生劈裂虚部保持简并此时系统处于PT对称态当g大于g_c时实部重新简并虚部劈裂成一对相反符号的值系统进入破缺态g g_c处就是奇异点EP。实操上我一般先粗扫g从0到0.02步长0.002把曲线整体轮廓画出来确定EP大致位置再在EP附近加密步长比如0.0002甚至0.00005。扫描结果是典型的“横杠”形结构实部在EP处收拢虚部在EP处张开。用一句话记忆就是对称相里频率劈开、损耗平衡破缺相里频率简并、增益损耗各自归边。3.3 不要只画曲线还要盯住场分布只看本征值曲线容易得出错误结论因为数值上本征值简并并不代表物理上模式被“冻结”。我每次扫描完都会挑三个典型g值远小于g_c、接近g_c、远大于g_c专门做场分布图。小于g_c时模式场同时覆盖两条波导看起来是一个扩展的超模接近g_c时场分布出现明显的空间再分布超过g_c后一个模式集中在增益波导中另一个集中在损耗波导中。这种空间局域化是PT对称破缺的直接证据也是论文里最有说服力的图。画场分布时建议用归一化电场模的平方即能量密度不要只看电场模。能量密度能更清晰显示模式在增益损耗通道中的占比。如果还想看手性或者相位信息可以叠加“表面图箭头”观察坡印廷矢量的方向分布。3.4 参数化扫描的工程习惯COMSOL的参数化扫描是一个成熟功能但有几个习惯值得养成。第一扫描类型不要选“所有组合”单独扫描g这一个参数就行避免无用计算第二每个参数点都保留求解结果用来做后续场分析第三扫描过程中如果出现断点大概率是求解器没收敛先检查EP附近网格再检查初始猜测值。四我一般会先关掉“智能网格细化”这种自动过程手动控制网格密度保证所有扫描点使用同一套网格否则本征值曲线的连续性会被网格变化带来的数值噪声破坏。4. 微盘体系实操回音壁模式与模式劈裂4.1 为什么微盘和波导的建模思路完全不同微盘谐振器依赖的是回音壁模式光在微盘边缘不断全内反射形成闭合环形轨迹。与波导最大的区别在于微盘的模式是离散的共振频率不是连续传播常数。所以仿真研究类型应该选“特征值/本征频率”而不是“模式分析”。微盘体系的PT对称结构常见有两种双微盘耦合结构一个增益微盘、一个损耗微盘以及单微盘内部沿角向做增益/损耗调制。前者建模更简单后者更能展示旋转对称性破缺的物理。我做的比较典型的模型是半径5μm、厚度200nm的微盘材料折射率实部2.0背景为空气。计算域背景给一个半径15μm的柱体外围包PML。微盘的Q值通常在10^4以上边界对计算结果的影响非常大PML几乎不可省略。4.2 增益损耗的角向调制单微盘内实现PT对称的方法是让增益和损耗沿角向交替分布。举个例子微盘上半部分θ在0到π折射率虚部为正下半部分θ在π到2π折射率虚部为负。用COMSOL表达式实现就是在材料折射率节点的虚部里写一个if(θπ, g, -g)θ通过atan2(y,x)得到。更精细的做法是按角向模式阶数m调制n_i g * cos(m * θ)这样既满足PT对称条件又只在特定m阶回音壁模式下起效。这里要特别注意COMSOL默认的坐标角度范围是-π到π用atan2时方向边界可能出现不连续这会导致材料参数在某个方向发生跳变进而激发伪模式。我建议写表达式时加上平滑过渡或者直接用两段条件剔除边界不连续。增益损耗边界是否是突变其实对物理结果也有影响真实器件往往是渐变掺杂所以也可以用cos渐变表达式来实现。4.3 模式手性交换EP点的隐藏证据微盘PT对称研究中EP点附近最漂亮的现象是两个手性模式顺时针和逆时针的频率简并以及手性的自发交换。在COMSOL里提取两个相邻特征频率把它们的复频率随g画到一起你会发现g小于g_c时两个频率实部劈裂、虚部简并g大于g_c时实部简并、虚部劈裂。虚部为正的模式实际是损耗模式虚部为负的模式实际是增益模式——这两个模式在EP处互换物理性质。需要多解释一句COMSOL求解出的复本征频率虚部符号同样遵循“正虚部衰减、负虚部增长”的约定。物理上我们更习惯说“模式Q值”换算关系是Q |Re(ω)| / (2|Im(ω)|)在EP附近Im(ω)趋近于0Q值会发散这也是理论文章里常说的“奇异点处Q值异常增大”的一个侧面。4.4 模式识别的实用技巧微盘的特征值求解会求出一大堆高阶模式、径向模式、甚至PML伪模。我踩过最狠的坑是算出来的最低频率模式表面看是基阶回音壁模式实际上场分布完全落在PML层里纯属伪模。辨认方法很简单把每个特征模式的场分布图打开真实的回音壁模式场是沿微盘边缘一圈一圈分布中心区域基本没有场强伪模则通常在PML区或背景边界上有明显增强。建立一个筛选习惯先求不带增益损耗的微盘本征模式和解析公式对比确认模式的阶数再开启增益损耗参数去追它们的演化轨迹能省下不少时间。5. 光纤体系实操双芯光纤的模间耦合与PT对称5.1 光纤建模的性价比方案光纤结构的特点是纵向尺寸远大于横向尺寸直接在3D里建完整光纤做仿真非常吃力。我对光纤体系做PT对称研究时99%的情况用二维截面模式分析就够了。建一个圆形包层里面放两个圆形纤芯设置不同的折射率虚部研究类型选“模式分析”求解得到两个超模的传播常数就能完整复现PT对称物理。典型的双芯PT对称光纤参数包层直径125μm纤芯直径8μm芯间距20μm工作波长1550nm。两个纤芯的实部折射率相同比如1.47包层折射率1.45一个纤芯虚部设成-ig增益另一个设成ig损耗。这个构型在COMSOL里的网格很好处理纤芯区网格0.5μm包层区网格2μm外围散射边界条件。计算量非常小几秒内完成一次模式扫描。5.2 传播常数曲线的EP特征与前两种体系完全一样扫描g后画出两个模式的复传播常数实虚部曲线依然能看到对称相到破缺相的转换。光纤体系特别适合和解析耦合模理论做对照因为耦合系数κ有解析估计式κ (2π/λ) * (n_core^2 - n_clad^2)/n_eff * overlap integral实际也可以靠COMSOL先算对称模式无增益损耗的频率劈裂反推κ。EP条件在理论上是γ 2κ这里γ对应增益损耗系数与g的关系就是γ 2 * k0 * g按平面波近似具体看定义。我做这块时习惯在同一个图里同时画解析解和COMSOL数值解。两者趋势一致证明建模无误如果形状都不同先回头查符号约定和网格。5.3 三维光纤传播仿真与边界条件如果你想更直观地看到功率沿光纤长度的演化就需要做3D模型。3D双芯光纤仿真的计算量很大我的建议是压缩光纤长度只算几百微米到几个毫米端面加端口激励侧面包PML吸收边界。重点关注两个端面的功率分配以及光纤中沿传播方向的坡印廷矢量分布。这个做法可以验证“单向无反射”这类PT对称特有现象从增益端入射和从损耗端入射的透射率不再相同这是实验上最容易被测量到的特征。5.4 光纤体系里的实际物理限制用COMSOL做光纤仿真有一个容易被忽略的问题真实光纤里增益介质的饱和效应和泵浦分布并非均匀。COMSOL里如果把增益简单写成常数虚部和实验结果对比时会出现偏差。可以后续引入速率方程耦合但这是一个不小的工程。对大多数原理验证型研究常数虚部模型已经够用只要在论文里把“理想化假设”写清楚就行。另外双芯光纤在制造中纤芯间距误差可能高达±1μm而EP位置对间距极其敏感建议仿真时做一下容差扫描看看EP点偏移量有多大——这种参数敏感性分析对实验合作方非常有价值。6. 常见问题与排查技巧实录6.1 特征值曲线不连续或者乱跳ET、伪模、不收敛三大嫌疑。先做网格收敛性验证在两套网格下分别跑同一条扫描曲线如果曲线形状不一致问题在网格。如果曲线是跳跃性地从一条分支跳到另一条分支尝试把特征值求解器的“搜索基准点”锁定在目标模式附近避免求解器“跳到”其它模式上。还有一个常见原因扫描步长太大EP附近变化过于剧烈加密步长即可。6.2 永远找不到EP点先做最简单的解析验证把结构简化为两条平板波导的耦合模型估算理论EP位置然后检查COMSOL模型中两个波导的实部折射率是否完全一致PT对称要求实部对称哪怕差0.001曲线就不会完美简并。再检查增益损耗写入的符号是否相反。如果这些都确认过最后检查模式是否是同一个分支上的两个模式有时候你画的实部是一个模式、虚部是另一个模式两个模式根本没关联曲线自然对不上。6.3 微盘的伪模式怎么快速识别用排除法真实回音壁模式的场强集中于微盘边缘内侧中心区域极小伪模和背景模的场分布散落在PML区或计算区域边界。有一个更准的判断标准如果在不加PML仅仅加大计算域的模型中某个模式的频率变化超过0.1%基本可以判定为边界依赖的伪模。真模式在边界条件变化时频率几乎纹丝不动。6.4 COMSOL结果怎么和实验对比实验上直接测的是透射谱不是本征频率。想和实验对比一个办法是用频域扫描算透射谱另一个办法是通过Q值和劈裂宽度做间接换算。我一般会同时做两件事先用特征值研究找出模式频率和Q值再搭建一个带端口激励的频域模型算透射谱把特征的频率劈裂对应到透射谱的共振谷劈裂上这样文章出来的图可以直接和实验曲线对比审稿人基本挑不出毛病。最后分享一个小经验。PT对称系统研究在COMSOL里并不复杂真正决定成败的往往是一些看似琐碎的基础设置——符号约定、网格收敛性、PML参数、扫描步长。如果你在仿真的第一周就反复卡在“EP点找不到”这个事上别硬刚先从最简单的单通道增益/损耗验证开始一步步搭上来。我后来带的几个学生都是这样先花一天做一维验证再花两天搭双波导模型最后切到微盘和光纤基本一周内就能把完整的PT对称扫描曲线跑出来。希望这篇实操整理能让你少走我当年走过的弯路。