4、矩形波导模式分析：矩形波导的 TE/TM 模式、模式的正交性与归一化、模式求解方法（FEM、FDTD 简介）

矩形波导，说白了就是硅光芯片上最基础的传输线。你想想看，光在这么一个小矩形里跑，它可不是随便跑的，得遵循特定的“姿势”——这就是模式。今天咱们就聊聊矩形波导里的TE/TM模式，以及怎么把它们算清楚。

4.1 矩形波导中的TE/TM模式

先说说基本概念。矩形波导的横截面是矩形，光沿着长度方向（z轴）传播。根据电场和磁场在传播方向上的分量，模式分为两类：

• TE模式（横电模）：电场在传播方向（z方向）分量为0，磁场有z分量。说白了，电场是横着的，磁场是竖着带点“拧巴”的。
• TM模式（横磁模）：磁场在传播方向分量为0，电场有z分量。正好反过来。

我个人习惯用下标m和n来标记模式阶数，比如TE_mn或TM_mn。m和n分别代表x和y方向上半波长的个数。举个例子，TE₁₀模式就是x方向有1个半波，y方向是均匀的。我在项目中遇到过，很多新手搞不清这个下标，其实你想想看，它就像驻波在矩形腔里形成的“斑马条纹”，条纹数就是m和n。

这里有个关键点：矩形波导不能传输TEM模（横电磁模），因为TEM模需要双导体结构，而矩形波导是单导体。嗯，这个坑我踩过，刚开始做设计时总想着用同轴线的思路去套波导，结果发现根本对不上。

4.2 模式的正交性与归一化

模式的正交性，说白了就是不同模式之间“互不干扰”。数学上可以写成：

∫∫ (E_m × H_n*) · dS = 0  (当 m ≠ n 时)

这个性质非常重要。为什么？因为光在波导里传播时，不同模式之间不会交换能量。你想想看，这就像不同频率的电台信号，各自走各自的，不会串台。

归一化呢，就是把模式的功率归一化到1瓦。我建议在仿真时一定要做归一化，否则不同仿真工具算出来的场分布没法直接比较。归一化的公式是：

∫∫ (E × H*) · dS = 1

我曾经在项目里吃过亏：两个团队分别用不同工具仿真同一个波导，结果耦合效率对不上。后来发现是归一化方式不同——一个归一化到1W，另一个归一化到峰值场强为1。从那以后，我要求所有仿真必须统一归一化标准。

4.3 模式求解方法：FEM与FDTD简介

模式求解，说白了就是解麦克斯韦方程组。但解析解只适用于简单结构（比如矩形波导），复杂结构就得靠数值方法。这里介绍两种最常用的：FEM和FDTD。

4.3.1 有限元法（FEM）

FEM的核心思想是把连续区域离散化成小单元（网格），在每个单元上近似求解。我个人习惯用FEM做模式分析，因为它对复杂几何形状适应性强，而且精度高。

FEM的步骤大致是：

1. 把波导横截面划分成三角形或四边形网格
2. 在每个网格上假设场分布的多项式近似
3. 组装全局矩阵，求解特征值问题
4. 得到模式的有效折射率和场分布

举个例子，用COMSOL求解矩形波导TE₁₀模式的设置：

// COMSOL中设置矩形波导模式分析
1. 选择“电磁波，频域”物理场
2. 设置波导横截面尺寸：宽a=500nm，高b=220nm
3. 边界条件：完美电导体（PEC）
4. 网格：最大单元尺寸 λ/10
5. 求解：特征值求解，搜索有效折射率范围 1.5~3.0
6. 后处理：查看TE10模式的场分布

嗯，这里要注意：网格密度直接影响精度。我建议至少保证每个波长有10个网格，否则高阶模式的计算误差会很大。

4.3.2 时域有限差分法（FDTD）

FDTD是另一种思路——直接在时域上离散化麦克斯韦方程组。它不求解特征值，而是通过时域脉冲激励，然后做傅里叶变换得到频域响应。

FDTD的优势在于：

• 一次仿真可以得到宽频带结果
• 适合分析瞬态效应和非线性
• 实现简单，容易并行化

但FDTD也有缺点：

• 需要满足CFL稳定性条件（时间步长不能太大）
• 色散误差会累积
• 对精细结构需要大量网格，计算量大

我曾经用FDTD分析一个带光栅的波导，网格设得太粗，结果反射谱完全不对。后来把网格加密到λ/20，结果才收敛。所以我的经验是：FDTD的网格密度至少λ/15起步，关键区域要λ/20以上。

4.4 知识体系总览

为了让你更直观地理解本章的知识结构，我画了一张图：

这张图把本章的核心逻辑串起来了：从矩形波导出发，先分清TE/TM模式，再理解正交性和归一化这两个关键性质，最后掌握FEM和FDTD两种求解方法。你想想看，有了这个框架，后续学习硅光器件设计就会清晰很多。

---

公众号:蓝海资料掘金营，微信deep3321