企业官网建设流程全解析

1. 先把这个方向看明白：PT对称不是玄学，是“损耗换来的对称”

PT对称系统这几个字一听像是量子力学理论里的概念，但实际上这几年在光子学领域已经被玩出花来了。不管是波导、微盘还是光纤体系，PT对称的核心思路都一样：让系统在一个特殊的“增益-损耗配平”条件下表现出反直觉的光学行为，其中最著名的就是奇异点（Exceptional Point，EP）附近的本征态简并和自发对称性破缺。用一句稍微工程化的话说：如果一个光学体系的复折射率满足实部空间偶对称、虚部空间奇对称，那它就是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.45（SiO2），波导芯层折射率实部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(ω)趋近于0，Q值会发散，这也是理论文章里常说的“奇异点处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对称扫描曲线跑出来。希望这篇实操整理能让你少走我当年走过的弯路。