微环谐振腔仿真选型：波动光学与波束包络模块的适用边界

做微环谐振腔与环形波导耦合的仿真，有一个绕不开的选择题：到底用波动光学模块老老实实全波求解，还是上波束包络模块省计算资源？我在光子集成相关的项目里两条路都跑过，也看到不少同事在这两个模块之间反复横跳。一根直波导搭一个环形波导，耦合间隙只有一两百纳米，谐振波长对结构尺寸的敏感度可以到亚纳米级别——这种结构画起来极简，真正把传输谱算准，问题全藏在模块的默认假设和网格策略里。

先说结论：如果你的目标是微环谐振腔这种强反馈、准回音壁模式的器件，波动光学（全波）是稳定靠谱的默认选择；波束包络主要用于场包络沿某一大方向平滑演进的结构，强行套到环形波导上往往要在相位参考上做很多手脚，算出来还容易偏离谐振频率。下面把我的完整对比过程、设置细节和踩坑记录写出来，给正在两条路上纠结的朋友一个参考。

1. 微环耦合仿真真正难在什么地方：相位本征条件与耦合间隙

1.1 谐振条件与耦合区域决定仿真范围

微环谐振腔的基本结构不用多说：一个总线波导（bus waveguide）贴近一个环形波导，两者之间隔着几十到几百纳米的间隙。光从总线端口注入，一部分能量通过倏逝场耦合进环，环里转一圈回到耦合区，再和总线里的直通光干涉。系统满足谐振条件时，环内功率不断累加，总线输出谱就出现一个向下凹陷的谷，对应从“通端（THRU）”切换到“下载端（DROP）”的过程。

谐振条件本身很简单：2πR·neff = m·λ0，R是环半径，neff是环内模式有效折射率，m是整数。但这个条件对仿真有两个很麻烦的连带要求。第一，neff必须足够精确，因为谐振波长对有效折射率的变化几乎是等比放大的——微环的Q值动辄几千几万，neff差个10⁻⁴，谐振峰就可能漂出去几个自由谱范围（FSR）。第二，环是非闭合反馈结构，光场要满足周期性边界条件，这意味着任何全局相位近似都必须和环的周长严格匹配。

耦合间隙则是第二个精度敏感点。耦合系数κ由倏逝场在间隙中的指数衰减决定，衰减长度通常只有几十到一两百纳米。如果网格在这个区域不够密，或者边界条件引入了不自然的反射，算出来的耦合强度会明显偏大或偏小，直接反映在谐振峰的深度、半宽和插入损耗上。

1.2 从模块角度看两种方法对“解的形式”的不同假设

COMSOL里的波动光学模块和波束包络模块，底层都是从麦克斯韦方程组出发，差别在于怎么处理“快速振荡”的分量。

波动光学模块直接以电场E为未知量，在频域求解矢量亥姆霍兹方程。网格必须解析真实光场在一个波长尺度内的变化，通常每个有效波长至少要放6～10个网格单元，才能把相位累积和倏逝场形状抓准。这个要求听起来苛刻，但换来的好处是几乎不依赖任何先验假设——直波导、弯波导、耦合器、谐振腔，统统可以算。

波束包络模块做的是变量代换：把电场写成慢变包络E₀(r)乘以一个已知的快变相位项，比如exp(-i k·r)。理论上网格只需要分辨包络的变化尺度，可以比波动光学粗得多，计算量能下降一个数量级。但这个代换成立的前提是：你选的相位函数必须足够接近真实场的快速相位。换句话说，用户必须预判波在哪里、往哪个方向传、传播常数大概是多少。这个预判对直线波导很容易，对环形波导就非常尴尬——光在环里是绕着圆心转的，传播方向一直在变，不能用单一的全局波矢去近似。

这就是微环仿真真正的难点：它不是网格一颗就完事的问题，而是“解的形式假设”是否成立的问题。后面两个模块的全部区别，几乎都从这里派生出来。

2. 波动光学模块的全波路线：设置、网格与边界条件实测

2.1 界面选择与二维模型参数

我做微环仿真时用的还是最常见的一套：二维模型，射频模块或波动光学模块下的“电磁波，频域”物理场接口。不要被模块名字误导，射频模块的电磁波频域接口做光学波长完全没问题，COMSOL只是按常规应用场景给出了不同材料库和默认表达式。波动光学模块在光学材料色散、各向异性晶体、非线性极化等场景更方便，但对无源微环来说，两者等效。

几何尺寸按硅光常见参数来：芯层宽度w=500 nm，环半径R=5 μm，耦合间隙g=200 nm，芯层折射率3.45，包层1.45，工作波长1.55 μm附近。二维模型里把面外方向设为TE偏振（电场垂直纸面）或TM偏振，取决于你要仿真哪种模式。很多初学者会在这里栽跟头：同一个几何，TE和TM的有效折射率差别不小，耦合系数可以差出一倍，谐振峰位置也跟着变，必须先想清楚自己要哪一支模式。

研究类型选“频域”做参数化扫描，频率范围覆盖谐振峰附近，例如150 THz到162 THz。也可以先做“特征频率”研究，快速看环内存在哪些本征模式，然后再转到频域做精确的传输谱。

2.2 端口边界与PML的细节偏好

端口设置是耦合仿真里最容易出错的地方。总线波导两端各放一个端口。输入端选“端口”类型，勾选“激活端口”，模式类型用“数值”——让COMSOL根据总线波导的截面自动计算本征模式。通端和下载端同样设置为“端口”，但不激活激发，只用来读透射功率。这样能在后处理里直接得到S参数，比自己在边界上积分电场省事得多。

端口方向要特别注意：数值端口会假定波导截面是一段规则区域，如果端口边界距离耦合区太近，会截断高阶倏逝场，造成耦合系数的数值误差。我通常让端口离耦合区至少一个波长以上，给场的建立和衰减留出盈余。

外层边界用完美匹配层（PML）。有人认为小模型可以用散射边界条件顶替，经验是：高Q微环的辐射和泄漏场非常弱但又真实存在，散射边界在斜入射下的吸收效果不如PML平缓，容易带来虚假反射，在谐振谷附近制造出肉眼看不出来的小波纹。PML厚度取一个有效波长左右，外层再用“映射”网格划分，保证PML内部网格一层层推进而不是无规则三角形。

2.3 特征频率与频域扫描的取舍

先做特征频率研究往往能少走很多弯路。特征频率求解给出环内的本征模场分布和复数频率，实部对应谐振频率，虚部对应损耗或辐射导致的线宽。它不需要激励端口，直接把环当作封闭结构看，非常适合验证几何、折射率设置是否合理。

频域扫描则给出连接外部总线后的完整响应，能看到谐振谷的深度和下载端的上升，这是特征频率看不到的信息。具体做法是在“研究”设置里用参数化扫描，频率步长按共振线宽决定。Q值预期5000以上时，1.55 μm处线宽约0.3 nm，换算成频率约25 GHz，频率扫描步长必须小于这个量级，我通常取2～5 GHz，把100 GHz范围扫开。

求解器方面，二维模型优先选MUMPS直接求解器。迭代求解器（如GMRES）在大规模三维模型里有用，但在二维微环问题里反而可能出现收敛慢或残差残留的问题，直接求解器的代价完全可以接受。我实测过：连PML在内约6万个三角形网格元素，自由度在30万量级，单个频点用MUMPS求解只需要2～4秒，全程扫100个频点也就是几分钟。这个成本根本不需要省。

2.4 网格分配的实操经验

网格分配的原则是“局部加密，整体宽松”。总线波导和环波导的芯层内部，最大单元尺寸设为中心波长有效值的1/8以内，即约60～80 nm。耦合间隙区域必须加密到10～20 nm，才能把倏逝场的指数衰减尾部抓准。弯曲波导内侧和外侧光场分布不对称，曲率变化越大的地方单元要越小。

我这里有个笨办法：先把全局网格放宽跑一版，再把耦合区加密一倍跑一版，对比谐振波长的变化。如果移动超过0.1 nm，就继续加密；如果几乎不动，就可以认为网格收敛了。微环仿真的误差通常不是网格数量不够，而是网格分布不合理——大量三角形白白放在了远离波导的大片包层区域，真正需要分辨的耦合区反而很稀疏。

3. 波束包络模块解析：它真正省力在“哪一步”

3.1 慢变包络假设和网格收益的来源

波束包络界面的核心思路，是把电场改写为：E(r) = E₀(r)·exp(-iφ(r))。其中φ是用户指定的快速相位函数，E₀是慢变包络。代入亥姆霍兹方程后，方程主体变成对E₀的约化形式，网格只需要分辨E₀的空间变化，不再需要分辨一个波长内的快速振荡。如果φ选得准确，包络可能只是一个缓慢扩散的场型，网格尺寸可以放大到波长量级甚至更大，自由度随之暴降。

COMSOL官方资料里，波束包络最经典的案例是长距离直波导、倾斜入射光束、规则周期结构这类场景。在这些场景中，光的传播方向几乎不动，相位函数是全局线性项，比如φ=kx·x+kz·z，用户只需要告诉软件一个合理的波矢分量，剩下的事情全部交给包络求解。

3.2 为什么微环的弯曲波导会戳破包络假设

把同样的思路移到微环上，麻烦立刻浮现。环内光场沿角向传播，相位实际是角向坐标的函数：φ≈βRθ，其中β是环内模式的传播常数。这个相位不是空间中的单一直线，而是绕一圈后需要满足周期性条件。当你尝试用全局平面波相位去近似它，包络项会自动吸纳剩余的快速振荡，结果就是包络不再“慢”，网格收益被吃掉大半。

有人会说，那我可以直接用角向相位mθ来定义啊，其中m是整数环模式阶数。听起来可行，细想还是有坑：mθ这个相位只对“纯环回音壁模式”成立，而微环仿真通常还要包含总线波导的直通场。总线里的光用平面波相位，环里的光用角向相位，两种相位参考在同一个物理接口里叠加，界面上必然出现不可忽略的相位间断。再加上耦合区存在双向传输和强近场相互作用，慢变包络近似赖以成立的“单方向、缓变包络”前提已经崩坏。

算出来的结果也印证了这一点。我在类似结构的2D模型里试过用波束包络接口跑频域扫描，谐振峰的位置和全波结果相差往往超过1 nm，峰深也明显偏浅。这种误差很难通过网格加密修复——根子不在网格，在于相位函数的选取和真实场不匹配，属于模型层面的误差。

3.3 波束包络仍然擅长耦合类问题的几个场景

公平地讲，波束包络并非对耦合类问题一律无能为力，只是适用边界很清楚。

一是长直定向耦合器。两根直波导平行放置，光沿同一方向传播，耦合长度可能达到几百微米甚至毫米。全波扫描在这种结构上计算量非常大，波束包络却能用一个粗糙得多的网格得到收敛结果，因为相位函数可以很好地近似为βx。

二是周期性波导或光栅结构。如果能提前定义光栅矢量和占空比，慢变包络可以清晰表现出耦合模在空间上的缓慢交换，适合研究布拉格反射和超结构波导。

三是大曲率半径的弯曲波导。当环半径远大于弯曲波导模式的有效横向尺寸时，局部场分布和直波导差别很小，波束包络可以当作“直波导的准全局近似”来用。但此时你要清楚地知道自己在牺牲对谐振精确位置的预测能力。

我的建议很直白：如果模型的归一化曲率（有效半径/波长）不够大，或者目标指标涉及谐振波长、Q值、消光比的定量误差，那就别在波束包络上耗时间。

4. 同一微环模型两边同时跑的实测数据对比

4.1 对比条件与参考结果

为了把两个模块的差别量化，我搭了一个统一基准：硅芯微环，R=5 μm，波导宽度500 nm，耦合间隙200 nm，包层折射率1.45，中心波长1.55 μm，二维TE模式。频域扫描范围覆盖两个FSR，步长5 GHz。

全波模型在波动光学接口下跑通，网格约6万三角形单元，自由度约35万，总计算时间在普通八核工作站上约5分钟完成全部扫描。得到的透射谱在1548.2 nm附近出现明显谐振谷，谷深约22 dB，Q值约5600。

波束包络模型把相位参考设为沿x方向的全局波矢，网格放松到全波模型的约1/8，自由度降到约5万。结果却不乐观：谐振峰移到1549.6 nm附近，整体红移超过1 nm，谷深只有约12 dB。误差来源非常清晰：全局波矢对环内角向模式完全失配，总线和环的相位接地不一致，导致耦合干涉强度被系统性削弱。

4.2 横向对比表

这张表不是为了说明波束包络一无是处，而是提醒大家：省下来的计算时间，本质上是把“对相位先验的建模”偷偷推给了用户。当你对结构理解不深时，这个先验反而才是最难的部分。

4.3 计算资源的取舍心态

有人说微环无源仿真用3D也不怕，其实大可不必。一个二维微环模型在很多情况下已经能给出足够好的初步设计参考，尤其是TE/TM明确的片上集成方案。3D全波虽然更能反映侧壁粗糙度、衬底泄漏等真实效应，但网格量瞬间翻几十倍，参数化扫描变得极慢。我的做法是：先用2D全波把尺寸、间隙、折射率窗口扫出来，确定工作点后，再针对选定的一组几何做3D验证。

这种“2D全波为主、3D验证为辅”的策略，比在单一模块里死磕更高效。波束包络则留给真正需要模拟超大结构（毫米级）且场方向可预判的场景，不要在微环谐振上强行平替。

5. 传输谱与Q值提取：不管用哪个模块都要做对的三件事

5.1 频域扫描步长与共振峰拟合

从COMSOL里导出S参数后，谐振峰是一条密集的离散点序列。直接取最小值附近的点当“峰点”会带来误差，正确做法是先用洛伦兹函数拟合谱线，再从中提取中心频率f₀和半高宽Δf，Q=f₀/Δf。

扫描步长必须明显小于预期线宽。Q=5000的环在1550 nm附近的线宽约0.31 nm，如果参数化扫描步长是1 nm，谐振峰可能只被两三个点捕获，拟合出来Q值误差能到百分之几十。我一般先粗扫一个FSR，找到谐振区后再缩小范围做精细步长扫描。

5.2 网格收敛检验与求解器选择

全波仿真收敛性检验是必须写的作业。我第一次做微环时偷懒，用了默认网格，算出来的谐振谷奇奇怪怪，后来发现是耦合间隙只有两个网格单元，倏逝场尾部被强行做平了。正确的检验方法：固定所有物理设置，只把耦合区最大单元尺寸从20 nm变成10 nm再变成5 nm，对比谐振波长和Q值变化。变化小于2%才算过了网格收敛这一关。

求解器选择也值得多提一句。二维模型用MUMPS直接求解器最省心，数值稳定性好；三维模型才需要换迭代求解器。千万不要在二维小模型里为了省内存强行用迭代求解器，很容易出现残差项收敛不足导致的“隐形”误差，谱线看起来光滑，其实每个点都漂了一点点，拟合出的Q值假性偏高。

5.3 从仿真数据反推耦合系数与后续外推

拿到透射谱之后，可以根据谐振谷深度和背景损耗反推耦合系数κ和环内损耗α。这一点对工程很有用：COMSOL全波仿真本质上帮你准确计算了“电磁结构”，但系统级设计（比如多个微环级联、调制器带宽估计）不适合每次都跑全波。

我习惯的做法：利用2D全波模型在几个不同耦合间隙下得到S参数，提取κ随间隙的变化曲线，再把这个κ放进解析的微环传递函数模型里。后面的系统级扫描全部用解析模型，一分钟能把整个参数空间跑完。这也算是“仿真资源要给在刀刃上”的一个好例子。

5.4 偏振与维度检查

二维模型的“面内/面外”设置直接决定结果是TE还是TM。硅波导中两者有效折射率差0.3～1%很正常，耦合系数的差别可能更大。如果仿真结果和实验对不上，第一步就检查偏振设置，不要急着怀疑模块选错。

另外，二维模型天然忽略侧壁粗糙度和垂直方向的模场挤压效应。对于紧束缚的硅线波导，二维近似相当可靠；对于弱束缚的氮化硅或聚合物波导，模场垂直分量大，二维结果的绝对值要保守看待，最好做一次三维全波锚定。

6. 最终实操建议：什么时候坚定选波动光学，什么时候可以考虑波束包络

把两种模块跑过一遍后，我对这两个接口的定位有了比较清晰的判断。

微环谐振腔与环形波导耦合这类问题，只要定量要求到达谐振波长的亚纳米级别，直接选波动光学模块。它的成本固然高一些，但二维模型完全可承受，换来的却是可靠性和直观性：不需要刻意定义相位参考，不需要担心全局波矢和弯曲波导失配，所有物理细节都在电场解里自然呈现。尤其是设计关键指标围绕Q值、消光比、耦合间隙公差时，全波结果能给你一个清晰的benchmark，后续所有解析近似都可以以它为准。

波束包络模块则建议用在“大而不弯”的结构里。长距离直波导耦合、周期光栅辅助耦合、大面积超表面单元阵列，这些场景中相位参考容易定义、慢变包络假设基本成立，能省下一个数量级的自由度。如果你要做的是毫米级长度的模斑转换器或绝热耦合器，波束包络几乎是最优选择。

如果你实在想在微环里尝试波束包络，请先回答自己三个问题：环内模式阶数m是多少？总线波导的传播常数是多少？两个相位参考在耦合边界上如何平滑过渡？这三个问题都能答清楚，说明你对结构理解已足够深，此时用波束包络可能也能得到可用的结果；答不清楚，就老老实实回全波。

最后分享一个我每次都会检查的小细节：无论用哪个模块，记得在结果里看场分布图而不是只看谱线。谐振波长对得上但场型不对，多半是相位假设或边界条件出了隐患；谱线漂亮但场图混乱，说明后处理里端口归一化有问题。先看场、再看谱，这两个都通过，仿真结果才算真正可交付。