3. 平板与条形波导的解析求解

在硅光波导设计中，三层平板波导是最基础、也最值得先吃透的解析模型。它由芯层（折射率 n₁，厚度 d）、上包层（n₂）和下包层（n₃）构成，光沿 z 方向传播，横向约束发生在 x 方向。虽然实际硅光波导多为条形结构，但平板波导的色散方程、截止条件和有效折射率求解方法，是理解条形波导模式行为的起点。本节围绕这三件事展开：写出色散方程、判断模式是否截止、用数值方法求出 n_eff 并画出模式图。

一、三层平板波导的色散方程

设光在芯层内为振荡场，在上下包层内为指数衰减场。定义横向波数 κ、上下包层衰减常数 γ₂、γ₃，以及传播常数 β = k₀n_eff，其中 k₀ = 2π/λ。芯层内的场可写成 cos(κx) 或 sin(κx) 形式，对应偶模与奇模。由芯层与包层界面处电场（TE）或磁场（TM）切向分量连续，可导出 TE 模色散方程：

TE 偶模：κ·tan(κd/2) = γ
TE 奇模：-κ·cot(κd/2) = γ
其中 γ 由上下包层共同决定：
γ = (γ₂ + γ₃)/2  （对称波导时 γ₂ = γ₃ = γ）

TM 模的色散方程形式相同，但每个界面处多一个折射率平方比因子：

TM 偶模：κ·tan(κd/2) = (n₁²/n₂²)·γ₂ 与 (n₁²/n₃²)·γ₃ 的等效组合
对称波导下：κ·tan(κd/2) = (n₁²/n₂²)·γ

这里需要强调一个关键约定：κ 与 γ 并非独立变量，它们通过 k₀ 和折射率耦合：

κ² = k₀²n₁² - β²
γ₂² = β² - k₀²n₂²
γ₃² = β² - k₀²n₃²

因此色散方程本质上是关于 β（或 n_eff）的超越方程，无法直接解析求解，必须借助数值方法。这也是为什么“有效折射率数值求解”是本节的核心技能。

二、TE/TM 模式截止条件

模式截止意味着该模式不再被芯层约束，能量向包层泄漏。对三层平板波导，截止发生在 β = k₀n₂（假设 n₂ ≥ n₃，即上包层折射率较高的一侧先截止）。此时 γ₂ = 0，色散方程退化为：

TE 偶模截止：κd/2 = 0, π, 2π, ...
TE 奇模截止：κd/2 = π/2, 3π/2, ...
即：κd = mπ（m = 0, 1, 2, ...）

由此可得截止条件：

V = k₀d√(n₁² - n₂²) = mπ
其中 V 为归一化频率，m 为模式阶数

对 TE₀ 模，m = 0，V = 0 即截止，说明基模 TE₀ 没有截止频率，只要 d > 0 就能存在。对 TE₁ 模，需要 V ≥ π 才能传播。TM 模的截止条件类似，但截止 V 值略大，因为 TM 模在界面处的场分布更“敏感”，需要更强的约束。实际设计中，若希望单模工作，应使 V < π，即：

d < λ / (2√(n₁² - n₂²))

例如 λ = 1550 nm，n₁ = 3.48（Si），n₂ = 1.44（SiO₂），则 √(n₁² - n₂²) ≈ 3.17，单模条件 d < 1550/(2×3.17) ≈ 244 nm。这解释了为什么硅光波导的典型截面尺寸在 220 nm × 500 nm 量级——它同时满足单模和强约束。

三、有效折射率数值求解与模式图绘制

数值求解 n_eff 的常用方法是二分法或牛顿法。以对称平板波导 TE 偶模为例，构造关于 β 的误差函数：

f(β) = κ·tan(κd/2) - γ
κ = √(k₀²n₁² - β²)
γ = √(β² - k₀²n₂²)

β 的搜索区间为 [k₀n₂, k₀n₁]。在该区间内 f(β) 单调变化，二分法可稳定收敛。具体步骤：取中点 β_mid，计算 f(β_mid)，若 f > 0 则 β 偏小，调整下界；若 f < 0 则 β 偏大，调整上界。迭代至 |f| < 1e-6 或区间宽度小于 1e-4 k₀。得到 β 后，n_eff = β/k₀。

模式图绘制则是在求得 n_eff 后，回代计算 κ 和 γ，写出芯层与包层内的场分布：

芯层：E(x) ∝ cos(κx)  （偶模）
上包层：E(x) ∝ cos(κd/2)·exp[-γ₂(x - d/2)]
下包层：E(x) ∝ cos(κd/2)·exp[γ₃(x + d/2)]

将 x 从 -3d 到 +3d 采样，即可画出归一化电场分布。观察模式图时，重点看三个特征：芯层内振荡周期数（对应模式阶数）、包层内衰减快慢（对应约束强弱）、以及场是否在包层中拖尾过长（对应接近截止）。

上图左侧给出三层平板波导的折射率剖面，右侧上方是色散方程与二分法求解流程，右侧下方是 TE₀ 与 TE₁ 的模式场分布对比。可以看到 TE₀ 在芯层内只有一个波峰，包层内衰减迅速；TE₁ 在芯层内出现两个波峰，且包层内拖尾更长，说明其约束更弱、更接近截止。实际编写求解程序时，建议先扫描 V 值从 0 到 10，记录每个 V 下存在的模式数量，再针对目标 V 精确求解 n_eff，这样既能画出完整的模式图，也能直观看到模式截止的临界点。

从平板到条形波导，解析求解的思路可以迁移：条形波导可视为两个方向上的平板约束叠加，有效折射率法（EIM）先在一个方向求 n_eff，再将其作为另一个方向平板波导的芯层折射率，从而近似得到二维约束下的模式。掌握平板波导的色散方程、截止条件和数值求解，是后续理解条形波导模式工程的基础。