3. 光波导电磁理论基础

设计光波导时，我们真正关心的并不是“光在玻璃里怎么走”这种直觉描述，而是电磁场在介质分界面约束下能形成哪些稳定分布、这些分布以多快速度沿传播方向前进、以及它们对波长和结构尺寸有多敏感。要回答这些问题，必须回到麦克斯韦方程，把它化为波动方程，再引入模式、有效折射率和传播常数这几个核心量。本节的目标是建立一套能直接用于波导设计的物理图像，而不是做纯数学推导。

一、从麦克斯韦方程到波动方程

在无自由电荷、无自由电流的介质中，麦克斯韦方程组的微分形式为：

∇·D = 0
∇·B = 0
∇×E = -∂B/∂t
∇×H = ∂D/∂t

其中 D = εE、B = μH。对第三式取旋度并代入第四式，利用 ∇×(∇×E) = ∇(∇·E) − ∇²E，在折射率缓变或分段均匀的波导中可近似得到矢量亥姆霍兹方程：

∇²E + k₀²n²(r)E = 0
k₀ = 2π/λ₀

这里 k₀ 是真空波数，n(r) 是横截面上的折射率分布。对硅基波导，典型结构是 Si 芯层（n≈3.48）被 SiO₂ 包层（n≈1.44）包围，折射率差很大，因此场被强烈约束在芯层附近。若只关心沿 z 方向传播的导模，可写成 E(x,y,z) = E(x,y)·exp(−jβz)，代入后得到横向本征值问题：

∇ₜ²E(x,y) + (k₀²n²(x,y) − β²)E(x,y) = 0

这个方程说明：给定截面折射率分布和波长，只有某些离散的 β 值能让场在横向“闭合”并满足边界连续条件，这些 β 就对应不同的模式。

二、模式：波导允许的稳定场分布

模式是波导横截面上一种自洽的电磁场分布，它在传播过程中形状不变，只整体乘以相位因子 exp(−jβz)。可以把它理解为“波导能容纳的驻波图案”。对矩形波导，常用两个整数指标标记模式，例如 TE₀₀、TE₁₀、TM₀₀ 等；对光纤则常用 LP₀₁、LP₁₁ 等。模式是否被截止，取决于归一化频率 V：

V = (2πa/λ₀)·√(n_core² − n_clad²)

以 SOI 条形波导为例，若芯层宽 500 nm、高 220 nm，在 1550 nm 波长下，基模 TE₀₀ 被良好约束，而高阶模通常被截止，这正是单模波导设计的基本依据。模式概念的价值在于：它把“光怎么走”转化为“哪些场分布能存在、它们的 β 是多少、约束强不强”。

波导截面折射率分布与模式约束关系 截面结构（x-y 平面） 上包层 SiO₂ n≈1.44 Si 芯层 n≈3.48 下包层 SiO₂ n≈1.44 w = 500 nm h = 220 nm 基模 TE₀₀ 场分布 芯层区域 x E 关键参数关系 β = k₀ · n_eff n_clad < n_eff < n_core 1.44 < n_eff < 3.48 V = (2πa/λ₀)·NA V < 2.405 → 单模 V 越大 → 约束越强 沿 z 方向传播：E(x,y,z) = E(x,y)·exp(−jβz) 传播方向 z 相位因子 exp(−jβz) β = 2π·n_eff / λ₀ λ_guide = λ₀ / n_eff

三、有效折射率与传播常数

传播常数 β 描述相位沿 z 方向的变化率，单位是 rad/m。为了把“波导中的光”与“真空中的光”联系起来，定义有效折射率：

n_eff = β / k₀ = β·λ₀ / (2π)

n_eff 的物理意义是：波导中模式感受到的等效折射率。它必须满足 n_clad < n_eff < n_core，否则场无法被约束。对上述 SOI 波导，TE₀₀ 模在 1550 nm 下 n_eff 通常在 2.2～2.6 之间，具体值取决于宽度和高度。n_eff 越大，模式越靠近芯层中心，约束越强；n_eff 越接近包层折射率，场越向包层扩展，弯曲损耗和串扰都会增大。

另一个常用量是模式传播长度 Lπ = π/β，它表示相位变化 π 所对应的距离。在定向耦合器中，两个波导间的耦合长度就与奇偶超模的 β 差直接相关：

L_c = π / (β_even − β_odd)

这说明 β 不仅是“相位变化率”，还是耦合器、谐振器、延迟线等器件的核心设计参数。

四、把理论用于波导设计

实际设计时，通常按以下顺序建立物理图像并确定参数：

  1. 确定材料与波长：例如 SOI 平台、1550 nm，得到 n_core 和 n_clad。
  2. 选择截面尺寸：以单模条件为目标，估算 V 值，使高阶模截止。对 500 nm × 220 nm 的 SOI 条形波导，TE₀₀ 基模在 1550 nm 附近保持单模传输。
  3. 求解模式：用有限元或有限差分模式求解器计算 E(x,y) 和 β，得到 n_eff。示例中若 n_eff ≈ 2.45，则 β = 2π×2.45/1550 nm ≈ 9.93×10⁶ rad/m。
  4. 判断约束与损耗：观察场分布是否集中在芯层。若包层中场强明显，说明约束弱，弯曲半径需加大。
  5. 用于器件设计：把 β 代入耦合长度、谐振条件、相位匹配条件，完成耦合器、微环等结构初值设计。
设计提示：n_eff 不是材料的固有常数，而是“结构 + 波长 + 模式”共同决定的结果。改变波导宽度 10 nm，n_eff 可能变化 0.01 量级，进而显著改变耦合长度。因此版图设计时必须把工艺线宽偏差纳入考虑。

五、边界条件与模式正交性

在芯层与包层界面，电磁场切向分量连续：E_t 连续、H_t 连续。这一条件决定了模式场的形状和 β 的离散取值。不同模式之间满足正交关系：

∫∫ E_m × H_n* · ẑ dA = 0  (m ≠ n)

正交性意味着模式之间在理想波导中不会交换能量。当波导结构发生扰动（如弯曲、宽度变化、邻近波导耦合）时，正交性被破坏，才会发生模式耦合与能量转移。这也是定向耦合器、布拉格光栅等工作机理的根源。

常见误区：把 n_eff 当成“芯层折射率”或“平均折射率”来粗略估算。对高折射率差波导，场分布强烈依赖截面形状，简单平均会带来 10% 以上的误差，足以使耦合长度偏离设计值一倍以上。必须通过模式求解获得 n_eff。

总结来看，麦克斯韦方程给出波动方程，波动方程在波导截面上的本征解就是模式，每个模式对应一个传播常数 β 和一个有效折射率 n_eff。β 决定相位如何积累，n_eff 决定模式被约束的强弱，二者共同构成波导设计的物理基础。后续的耦合器、弯曲、谐振器等结构设计，本质上都是在操控这些模式的 β 和场分布。

← 上一章 📖 返回目录 下一章 →