4. 平板与条形波导模式求解

在量子光子芯片中，波导是承载光信号最基本的“导线”。设计的第一步，往往不是直接画版图，而是回答一个更基础的问题：给定材料折射率和目标波长，波导的芯层厚度、宽度取多少，才能只支持一个模式？这一节围绕三层平板波导的解析求解、条形波导的近似方法，以及单模条件与截止条件，给出可操作的尺寸初选思路。

一、三层平板波导的解析模型

三层平板波导由芯层（折射率 n₁，厚度 d）和上下包层（折射率 n₂、n₃）构成。光在芯层中全反射传播，在包层中呈指数衰减。设传播方向为 z，横向为 x，则电场可写成 E(x)e^(−jβz)。在芯层内，横向波数满足 k_x² = k₀²n₁² − β²；在包层中，衰减常数满足 γ² = β² − k₀²n₂²（下包层同理）。这里 k₀ = 2π/λ 为真空波数，β 为传播常数，有效折射率定义为 n_eff = β/k₀。

把边界处电场和磁场切向分量连续的条件代入，可得到 TE 模的特征方程：

tan(k_x·d) = (γ₂ + γ₃)·k_x / (k_x² − γ₂·γ₃)

TM 模形式类似，只是把 γ 换成 γ/n²。这个方程没有显式解，通常用数值方法（如牛顿迭代）或图解法求根。对设计而言，更实用的做法是引入归一化参数：归一化频率 V = k₀·d·√(n₁² − n₂²)，归一化传播常数 b = (n_eff² − n₂²)/(n₁² − n₂²)。对对称平板波导（n₂ = n₃），TE₀ 模的色散关系可近似写成：

b ≈ (1 + V²/ (1 + V²)) 的拟合形式，或直接用
V = (m+1)·π/2 作为第 m 阶模的截止条件

更常用的是截止条件：当 V < π/2 时，对称平板波导只支持 TE₀ 基模；当 V 超过 π/2，TE₁ 模开始出现。因此单模条件可写成 V < π/2，即 d < λ / (2√(n₁² − n₂²))。这就是平板波导单模厚度上限的解析来源。

二、条形波导的近似方法

实际芯片中的波导是条形（矩形截面），横向和纵向都受限，严格求解需要二维模式求解器。工程上常用两种近似：

• 有效折射率法（EIM）：把条形波导拆成两个方向的一维平板问题。先沿宽度方向求一次等效平板，得到等效折射率 n_eff,x；再把这个等效折射率当作新的芯层折射率，沿厚度方向再求一次平板，得到最终 n_eff。EIM 对宽而浅的波导精度较好，计算量小，适合初选。
• Marcatili 近似：假设模式在芯层内近似正弦分布，在四个包层区域指数衰减，把二维场写成两个一维场的乘积。它给出更直接的解析表达式，但要求芯层尺寸远大于波长，对亚微米波导误差偏大。

以硅波导为例：芯层 Si（n₁ ≈ 3.48），下包层 SiO₂（n₂ ≈ 1.44），上包层可为空气或 SiO₂。若目标波长 λ = 1550 nm，用 EIM 初选：先取宽度 w = 500 nm，沿宽度方向求等效平板，得到 n_eff,x 约 2.8～3.0；再沿厚度方向求平板，若厚度 h = 220 nm，可算出 n_eff 约 2.4～2.5，且 V 值落在单模区。这个 500 nm × 220 nm 就是常见的硅单模条波导初选尺寸。

三、单模条件与模式截止

单模条件不是单一公式，而是“V 值 + 纵横比 + 偏振”共同约束。对条形波导，常用经验判据是：

• 宽度 w 和厚度 h 都小于对应平板单模上限的约 1.2 倍；
• TE 基模与 TM 基模的有效折射率差不宜过大，否则偏振相关损耗上升；
• 高阶模截止：当 w 或 h 增大到某一阈值，TE₁ 或 TM₁ 出现，此时 V > π/2。

截止的物理含义是：当芯层尺寸缩小到一定程度，模式无法被束缚，能量泄漏到包层，β 趋近 k₀n₂，n_eff 趋近包层折射率。设计时要留出安全余量，通常让工作点离截止点有 10%～20% 的尺寸裕度，以抵抗工艺偏差。

四、尺寸初选的实际操作

把上述方法落到具体设计，可以按以下顺序推进。第一步，确定材料平台和波长。例如 SOI 平台，Si 芯层 n₁ = 3.48，SiO₂ 包层 n₂ = 1.44，λ = 1550 nm。第二步，用平板公式估算厚度上限：d < λ / (2√(n₁² − n₂²)) = 1550 / (2 × √(3.48² − 1.44²)) ≈ 1550 / (2 × 3.17) ≈ 244 nm。所以厚度先取 220 nm，留出约 10% 裕度。第三步，固定厚度，用 EIM 扫描宽度。宽度从 400 nm 到 600 nm 逐步计算 n_eff，观察 TE₀ 和 TE₁ 的截止点。通常 500 nm 附近 TE₁ 尚未出现，且 TE₀ 的 n_eff 约 2.4～2.5，模式束缚良好。第四步，用二维模式求解器（如 FDTD 或 FEM）复核，确认只有 TE₀ 和 TM₀ 两个基模，且两者 n_eff 差在可接受范围。

需要强调的是，解析法和 EIM 给出的是初选值，不是最终值。实际设计还要考虑弯曲损耗、耦合器接口、金属电极吸收等因素。但有了这个初选，后续优化就有了合理的起点，避免在二维仿真中盲目扫描。单模条件与截止条件，本质上是在“束缚足够强”和“只允许一个模式”之间找平衡，尺寸初选就是把这个平衡量化成具体的 d 和 w。