做PT对称光子学研究这几年几乎每周都会有人问我同一个问题这种带增益带损耗的结构到底怎么在Comsol里搭出来尤其最近做波导、微盘、光纤三类体系时我发现很多人的卡点根本不在于物理而在于“建模细节”和“参数设置”。这篇文章就围绕Comsol里PT对称系统的仿真思路展开重点是三类最常用的载体耦合波导、双微盘、以及光纤体系。我会按自己实际跑过的模型顺序来讲把从哈密顿量到仿真参数的映射关系、模式劈裂与奇异点的判断方法、增益损耗怎么用复数折射率加载、边界条件和扫描策略怎么设以及最常翻车的那些坑一次性讲清楚。无论你是刚开始接触PT对称、想复现文献里的结果还是已经建了模型但结果对不上这篇都值得你花十分钟看完。1. 先搞懂PT对称在光学里到底意味着什么1.1 从量子力学到光子学的映射PT对称这个概念最早出自量子力学核心是一个非厄米哈密顿量在空间反演和时间反演的联合操作下保持不变。但光学里我们很少直接解薛定谔方程而是解麦克斯韦方程。好在这两个方程在标量近似下有非常漂亮的类比亥姆霍兹方程里介电常数分布ε(r)如果满足ε(-r) ε*(r)就对应了PT对称条件。用更直白的话说一根波导的折射率实部是偶对称的虚部是奇对称的。实部偶对称意味着几何和材料的基本折射率左右一致虚部奇对称则意味着一边是增益、一边是损耗。比如两条波导并排左边加损耗右边加增益损耗系数和增益系数相等损耗和增益沿中间对称轴呈一个“奇函数”分布就构成了最典型的PT对称双波导单元。这个映射的好处在于量子力学里需要的粒子源和陷阱光学里用现成材料就能实现。增益可以用掺铒波导、半导体光放大介质或者拉曼泵浦来实现损耗就更容易了金属覆盖层、故意掺杂吸收材料、或者加工一点粗糙度都能引入可控损耗。1.2 模式劈裂、EP点与参数空间考虑两条耦合波导各自传播常数设为β0耦合系数为κ。左波导有损耗γ右波导有增益-γ按复数传播常数虚部来记哈密顿量写成2x2矩阵H [[β0 iγ, κ], [κ, β0 - iγ]]本征值解出来是β± β0 ± sqrt(κ² - γ²)当γ κ时平方根是实数两个模式的传播常数实部发生劈裂但损耗/增益的虚部完全相同。当γ κ时平方根变成虚数模式实部重新简并但是一个模式出现净增益、另一个出现净损耗。在γ κ这个临界点两个模式的本征值和本征矢量同时合并这就是奇异点EPExceptional Point。这里我要强调一个新手特别容易犯迷糊的点EP点不是简单的“模式交叉”而是“模式并合”。交叉点是两个模式擦肩而过但它们的本征态还是可以区分的在EP点尔米特体系里那种模式的正交性消失两个模式变成一个带缺陷的模式这正是PT对称系统能实现单模激光、增强传感灵敏度的物理根源。所以做仿真时我们最大的需求就是找到这个EP点对应的参数耦合间距、增益损耗系数、波长、几何尺寸。这正是Comsol能发挥价值的地方。1.3 Comsol在这类仿真里的定位与优势市面上的光学仿真工具有不少比如FDTD类软件和基于有限元的光学求解器。但PT对称问题有一个天然需求需要在复平面上连续跟踪本征值频率或传播常数并且直观地画出模式劈裂和EP轨迹。Comsol的模式分析研究步骤Mode Analysis和特征频率研究步骤Eigenfrequency可以直接求解本征值而复数折射率材料又允许我们给反射介质定义任意的增益或损耗。这一点很关键FDTD的增益模型处理起来比较麻烦容易发散Comsol的频域本征值框架则天然适合这个问题。另外一个优势是参数化扫描。我们可以把“损耗量”或者“波导间距”设成全局参数一次性扫描几十个点把本征值的实部和虚部随扫描参数变化画成曲线EP点就在那条劈裂曲线的汇聚处。这个流程我后面会在每一类体系里分别演示。2. 波导体系最直观的PT对称实验台2.1 双波导模型设计与参数估算我强烈建议第一次接触PT对称的人从双波导模型起步。原因很简单几何简单、物理清晰、计算便宜还能跟解析解相互验证。建模第一步是确定波导尺寸。我在1550nm通信波段做SOISilicon-on-Insulator平台时常用单模硅波导宽度500nm、高度220nm。硅的折射率在1.55μm大约3.45埋氧层和上包层都是二氧化硅折射率约1.44。两条波导之间的间距d是关键参数影响耦合系数κ。间距越小耦合越强PT阈值γ_PT就越大能在更大范围内看到模式劈裂。怎么估算κ最可靠的办法不是查表而是先用无增益无损耗的两个相同波导做一次模式分析计算两条相同波导形成的对称超模和反对称超模的传播常数差κ (β_symmetric - β_antisymmetric) / 2。我习惯在一开始就把这一步做掉用真实几何而不是理论估计来确定模型的“标定”参数。举个例子。两条500nm宽、220nm高的硅波导间距200nm时超模劈裂对应的κ换算成有效折射率差大约在10⁻³量级这意味着EP点对应的γ也在类似量级。当我把损耗设成有效折射率虚部0.001再把增益材料设成虚部-0.001基本就能在EP附近观察模式并合。2.2 Comsol建模步骤与关键设置打开Comsol新建模型时我建议直接用“电磁波频域ewfd”接口研究步骤选“模式分析Mode Analysis”。这里先说明一下模式分析求解的是指定频率下的本征传播常数非常适合周期均匀的波导截面或沿传播方向不变的结构。如果你的波导沿传播方向是缓变的、或者需要看S参数那就改用频域扫描加端口边界条件。两者用途不同别混。几何上画两个矩形作为波导芯层再加上一个大的包层矩形。包层尺寸至少留出2μm以上不然边界条件会影响模式分布。我在实际仿真里会把包层画成直径6μm左右的大方块外围再加一层完美匹配层PML厚度1~2μm。材料设置的关键来了损耗波导折射率写为3.45 i*γγ 0表示损耗。增益波导折射率写为3.45 - i*γγ 0表示增益。注意这里符号约定。按光学惯例折射率虚部为正代表材料吸收损耗虚部为负代表增益介质。很多教程里用 ε n² 直接设置介电常数也可以但要注意 ε 虚部的正负和折射率虚部之间有一个交叉项不如直接设置折射率直观。我建议在“材料”节点里直接把“折射率实部”和“折射率虚部”作为两个参数来设置符号一目了然。模式分析研究中需要设置“搜索本征值附近”的目标。如果你已经知道单根波导的有效折射率大约2.6左右那就把搜索目标设为 复值传播常数 β_target ≈ 2π/λ * 2.6比如在1550nm下大约1.05e7 rad/m。我习惯勾选“所需模式数”为6~10个并锁定在目标附近搜索避免被包层模式或PML模式干扰。最后做参数化扫描。在“全局定义”里设置参数gamma范围从0到0.003步长0.0002然后扫描。求解完成后打开一维绘图组把两个目标模式的有效折射率实部和虚部分别提取出来即可看到劈裂曲线和EP汇聚过程。2.3 结果解读模式演化与EP判定当你跑完参数扫描结果应该长这样gamma较小时两个超模的有效折射率实部明显劈裂虚部基本一致随着gamma增大实部差逐渐缩小到达临界gamma时实部合并虚部分裂一个模式净增益、一个模式净损耗。那个临界值就是EP点。这里有个非常实用的判据我实测下来比眼睛看曲线更可靠在EP点附近两个模式本征值的“间距”对参数变化的敏感度会急剧增大而且特征频率/传播常数的虚部会从一条平滑曲线变成分叉。如果你用严谨一点的数学工具可以计算两个本征模式的“态距离”曲率在EP附近会出现发散特征这在传感灵敏度分析里是核心指标。另外要提醒一个细节模式分析给出的首阶模式不一定是你要找的那两个耦合超模。包层里的高阶模式、PML模式甚至数值模式都可能混进来。我的习惯是先固定一个参数比如gamma 0单独把所有模式画一遍电场模确认模式形状再用模式跟踪工具或手动指定“依据模式编号”的方式锁定目标。千万不要无脑提取第1和第2个模式那样很容易在扫描过程中跟踪到别的模式上导致曲线跳变。波导体系的优势在复现性很强。你甚至可以先用解析的2x2矩阵把EP位置、劈裂量造出来再跟Comsol结果对比。两者偏差一般在百分之几以内这给了我很大的信心去推更复杂结构。3. 微盘体系回音壁模式的增益-损耗博弈3.1 为什么微盘适合研究PT对称波导体系的耦合是沿传播方向连续发生的而微盘体系是另一种典型场景模式被约束在圆形回音壁内两个相邻微盘之间通过倏逝场耦合。这种结构的PT对称性通常这样构造两个尺寸几乎相同的微盘一个覆盖增益材料另一个覆盖损耗材料盘间距决定耦合强度。微盘的好处是模式体积小、Q值高增益阈值低而且可集成度高。PT对称双微盘的一个经典现象是通过调节两个盘的增益差/损耗差可以控制模式劈裂并在EP点附近使其中一个模式的线宽急剧变窄从而实现单模激光或者超高灵敏传感。实验上已经有不少工作验证了这一点我在仿真里跑出来的结果也确实能复现模式合并和线宽窄化。不过要注意微盘系统比波导多一个维度——方位角模数m。同一个波长下存在不同m的回音壁模式两个盘的耦合会发生模式竞争。仿真时需要先确认关心的具体模式否则很容易在特征频率列表里看到一堆无关模式。3.2 双微盘建模网格、边界与参数扫描微盘模型我用的是“电磁波频域”接口加“特征频率”研究步骤。这里跟波导模式分析的区别是直接求解谐振频率本征值而不是传播常数。研究设置里仍然要指定搜索目标比如在1.55μm附近搜200 THz左右的特征频率。几何上画两个半径相等的圆半径我常用5μm厚度220nm二维模型就忽略厚度或者用有效折射率法把厚度等效进去。两个圆盘中心距设为11μm净间隙1μm保证耦合不是太强也不是太弱。外围设置大气包层和PML。PML对微盘仿真特别重要因为微盘辐射损耗很小本征频率虚部线宽主要来自辐射损耗和材料吸收。如果没有PML边界反射会让虚部算得完全不准。材料设置和波导体系类似一个盘的折射率虚部设为γ另一个设为-γ。这里有个操作细节微盘材料通常用SiN或Si基座是SiO2也可以通过“有效折射率”的方式把厚度效应压到一个二维模里。我在实际建模中如果是220nm SOI微盘就直接用硅折射率3.45然后把微盘和包层画在一个平面上不开三维模型。二维近似的计算代价低很多而且回音壁模式的本征频率和劈裂规律都是准的。想看偏振和厚度色散再开三维也不迟。网格划分时回音壁模式沿圆周有大量振荡需要保证圆周切向上至少20个网格点覆盖一个光学周期。我会用“物理场控制网格”先跑一遍如果模式发散或频率偏差大就换成“用户控制网格”对微盘内部加最大单元尺寸限制比如0.05μm。这种网格尺寸在1550nm下已经是相当保守的加密计算起来仍然很快。参数扫描的策略我建议这样先固定增益损耗γ为一个较小值比如γ 50 cm⁻¹这个单位换算后面讲扫描盘间距看模式劈裂随距离的变化再固定间距扫描γ找EP。每次都只扫一个参数结果才容易解释。3.3 从本征模式谱判断奇异点特征频率求解结果是一堆复频率f Re(f) iIm(f)。有损耗时Im(f)为负值场衰减Q值可以用 Q Re(f)/(-2Im(f)) 计算。增益介质存在时Im(f)可以为正值表示模式净增益在激光问题里正Im(f)就是阈值以下/以上的判定。在微盘双盘系统中你会在每个方位角模数m附近找到一对模式。跟波导体系类似当增益损耗较小的时候这对模式频率实部劈裂虚部几乎相等随着γ增大实部差缩小在EP点合并然后虚部裂开。我习惯把每个扫描点的两个模式都画在复频率平面上横坐标Re(f)纵坐标Im(f)看轨迹怎么随参数移动。EP点就是两条轨迹在复平面上回合的地方。这里要特别小心一个求解问题在EP附近矩阵变成亏损矩阵数值上会出现模式合并的奇异性导致求解器要么丢失模式、要么突然跳变。我的经验是扫描步长不要太粗在EP附近步长要加密同时调整“特征频率搜索”的目标值保持两个模式的求解在同一分支上。如果出现模式跳变可以尝试把目标复数频率设得更靠近上一频点或者把“求解器”中的“特征值变换”方式改成“标准”形式。还有一些实用技巧微盘的EP位置对PML厚度和包层尺寸敏感如果发现EP点对应参数特别离谱先检查PML是否足够吸取辐射场。另一个验证方法是在纯损耗无增益下双盘模式劈裂值应该与单盘通过超模计算出的κ一致从而反推EP阈值。4. 光纤体系从分布式耦合到周期结构4.1 双芯光纤与长周期光栅的PT对称实现光纤体系是我最近两年投入比较多的一块因为PT对称在光纤里能实现一些很有意思的器件比如模式选择性分束、损耗补偿、带宽调控。光纤体系主要有两类实现方式。第一类是双芯光纤。两根纤芯沿轴向并行结构天然就是一对耦合波导。只要一根纤芯通过掺杂稀土离子比如铒提供增益另一根纤芯引入一定损耗通过掺钴或者做侧抛处理就能在横向平面上构成PT对称系统。双芯光纤的好处是长度可以很长弱耦合情况下耦合系数κ可以控制得很小因而增益损耗阈值很低很多在集成光波导里难以观察的效应在光纤里会更容易看到。第二类是长周期光纤光栅。沿光纤轴向周期性地写入折射率调制如果这个调制的实部是偶对称、虚部是奇对称就能在正反向包层模式之间形成PT对称耦合。这类光栅常用于特殊滤波器和模式转换器件仿真维度也包括轴向的周期结构。不过它的建模思路跟横向耦合不太一样需要借助“周期边界条件”或者“耦合模方程”这里我主要说说更普适的双芯光纤。4.2 等效折射率法建模与仿真细节双芯光纤在Comsol里最简单高效的处理方式是用二维截面模型加“模式分析”画一个直径125μm的包层圆里面放两个直径8μm的纤芯间距比如15μm。包层折射率设1.444纤芯折射率高出0.005到0.01左右。这个折射率差决定了单根纤芯的模式约束能力Δn越大模场越紧耦合越弱。增益和损耗的设置同样是给两个纤芯加复数折射率虚部。注意光纤问题里虚部的单位通常用dB/m而Comsol需要的是折射率虚部换算关系为折射率虚部 吸收系数[dB/m] * λ / (40π * log10(e))。我一般在全局参数里定义一个单位换算直接输dB/m然后在材料属性里引用换算后的值。跑模式分析后两个基模有效折射率随传播方向的变化会出现在实部和虚部里。还有一个地方容易踩坑光纤包层半径很大自由空间波长只有1.55μm但是包层直径125μm如果不加限制网格量很容易爆炸。我的做法是把包层外部区域挖掉一圈用“散射边界条件”或者“吸收层”模拟无限包层。更省资源的办法是先跑单根纤芯模式把模场大致范围确定后只在这个范围里加密网格包层外圈用普通粗网格不会对模式劈裂产生显著影响。另外双芯光纤沿轴向很长模场在传播中会发生周期性拍频。在Comsol里直接建三维长光纤是极不划算的。正确做法是先在截面上算出两个超模的传播常数β1和β2然后耦合长度 L_c π / |β1 - β2|功率转移周期就出来了。如果要做传输谱或串扰曲线可以直接在脚本中用耦合模方程后处理不必跑三维的完整光纤。4.3 光纤中的实用指标透射率、串扰与带宽对双芯PT对称光纤我最关心的三个指标是透射率、模式串扰和带宽。透射率可以在“频域扫描”下用端口边界条件求S参数获得但更常用的是长光纤的耦合模方程解。你会发现PT对称双芯光纤有一个很有用的特性当增益损耗系数γ接近κ时原本在两个纤芯间周期性转移的光功率会变成“单向模式”——一个超模快速增长另一个被抑制。也就是说输出端的功率不再是两个纤芯轮流占优而是大部分被放大到增益纤芯一侧。这在模式分束或者放大器设计里是很有价值的。串扰指标的定义是一个纤芯模式泄漏到另一个纤芯的功率比。在PT对称配置下EP点附近串扰对增益差特别敏感所以如果把“损耗差”作为传感参量串扰的急变可以做成一个灵敏度很高的传感器。实测下来在EP附近把γ变化0.1dB/m透射串扰能变化几个dB这个斜率比非PT对称光纤要大很多。带宽方面我习惯扫波长从1500nm扫到1600nm观察串扰谱和透射谱的变化。由于材料色散和结构色散共同影响β1 - β2当工作点靠近EP时谱线会出现明显的窄带响应。做器件时得权衡太靠近EP灵敏度高但带宽窄离EP稍远带宽宽但灵敏度下降。这就需要按具体应用来定。我个人的建议是在用Comsol跑这种长距离弱耦合结构前先在脚本里用解析耦合模方程估算所需的κ和γ再回来校准仿真。这样效率最高也不用在三维模型里浪费大量内存。5. 仿真过程中的常见问题与排查技巧5.1 找不到目标模式怎么办这是我被问得最多的问题。尤其是模式分析时本来想看两条波导的基模超模结果求解器返回的全是包层模式或者场分布长得一团乱的模式。第一个方法是缩小搜索范围。把“搜索本征值附近”的目标值准确设为 β ≈ 2π/λ * n_eff其中n_eff是单波导的模式有效折射率可以先跑一根单波导模型来获得。这样一来求解器会优先接近这个目标的模式包层模式因为有效折射率低接近包层折射率会被排除在外。第二个方法是调高“所需模式数”。有时目标模式被排在列表第8个而你以为只会有前3个结果漏掉了。我一般设置所需模式数到10个跑完后再按有效折射率大小把基模挑出来顺便看电场分布验证。第三个方法是调整初始猜测。在特征频率研究中如果求解器从一个特别差的初值出发可能收敛到不想要的模式。把搜索目标设成复数例如 193.5e12 i*1e9可以更精准地锁定目标区域。这在微盘系统里尤其有用。5.2 内存不足与收敛性问题Comsol内存不够几乎是每个把模型从二维升到三维时必然碰到的问题。我的建议是优先用二维模型验证物理和参数规律再决定是否升级三维。波导体系、双芯光纤截面、微盘二维近似都能给出非常可靠的模式劈裂和EP规律二维跑通了三维才有意义。如果必须做三维可以用“物理场控制网格”降低非关键区域的网格密度然后在关心的波导区域添加“尺寸”节点加密。另外参数化扫描时如果每个点都重新剖分网格会很慢。把“扫描时使用当前网格”这个选项勾上只在几何不变的情况下跳过重新剖分速度提升非常明显。我实测一个双波导三维模型扫描30个点用这个设置能从几个小时降到几十分钟。收敛性上最麻烦的还是增益介质。增益太强时电磁场能量会指数增长求解器很容易发散。对策有一是初始参数从纯损耗或γ0开始逐步增大γ把上一步的解作为下一步的初值二是不要一次性把增益设到阈值以上太多先观察线性区的行为三是如果模型里不需要严格饱和模型就始终用等效复数折射率避免引入额外的非线性项导致收敛失败。5.3 从仿真到实验的差别与修正仿真跟实验总是有差距的PT对称系统尤其明显。实验里最难精确控制的就是增益介质的实际增益值。铒掺杂的增益曲线跟泵浦功率有关跟波长有关还跟模式场和掺杂截面重叠有关。仿真里直接设置的复数折射率虚部在实验里其实是一个随泵浦变化的量。我的处理方法是做一个“等效增益标定”先用同批次波导/微盘做一个无增益损耗下的损耗测量再通过带泵浦的透射谱拟合出等效折射率虚部最后把这个实验拟合值带回Comsol模型重新预测EP位置和模式行为。做完这一步之后仿真和实验的吻合度会明显提升。另外实验中的微盘半径、波导宽度都有加工误差。哪怕只有10nm的尺寸偏差也会让模式劈裂量变化百分之几甚至更多。建议在Comsol里直接把尺寸参数化仿真时加一个±5nm、±10nm的偏差扫描看EP点位置的偏移范围。这样做出来的器件鲁棒性分析在写论文或做设计交付时都更有说服力。有一点要特别提醒PT对称实验里EP点对参数极其敏感这是物理本质也是它的应用价值所在。正因为它敏感仿真里稍微偏离EP实验上可能就完全看不到模式简并。所以仿真时必须做容差分析把“加工偏差-EP漂移”曲线给出来别只盯着完美几何跑。最后说一个我自己的习惯每一次仿真结束我会顺手把材料折射率的虚部、参数扫描范围和目标模式编号记在一个表格里方便复现。PT对称系统参数多、模式多随手存一个项目模板能省很多重复劳动。后面如果环境允许我还会专门写一篇怎么在Comsol里用模式跟踪脚本自动截取EP位置的教程配合参数扫描能更高效地挖出整个参数空间的奇异点轨迹。