3. 光波导模式理论与仿真
光波导是集成量子光子器件的“骨架”，它把光场约束在微纳尺度内传输，并决定模式数量、模式形状、有效折射率与色散。理解波导模式，核心是回答三个问题：允许哪些模式存在（离散本征值）、每个模式的场分布是什么（本征函数）、不同模式之间耦合有多强（重叠积分）。本节从平板波导的解析解出发，过渡到矩形波导的有效折射率法，再给出用 FDTD/FEM 计算单模条件与模式场分布的实际做法。
一、平板波导：三层结构的解析模式求解
平板波导由芯层（折射率 n₁、厚度 d）与上下包层（n₂、n₃）构成，光在 y 方向无限延伸，只在 x 方向受限。设传播方向为 z，场按 exp[i(βz−ωt)] 变化，β 为传播常数，定义有效折射率 n_eff = β/k₀，k₀ = 2π/λ。在弱导近似下，TE 模的横向电场 E_y 满足一维亥姆霍兹方程，芯层内为振荡解，包层内为衰减解：
芯层: d²E_y/dx² + κ²E_y = 0, κ = k₀√(n₁² − n_eff²)
包层: d²E_y/dx² − γ²E_y = 0, γ = k₀√(n_eff² − n₂²) (上包层同理)由芯层—包层界面处 E_y 与 dE_y/dx 连续，得到 TE 模的特征方程：
tan(κd) = κ(γ₂ + γ₃) / (κ² − γ₂γ₃)该方程是超越方程，需数值求解。工程上更常用归一化变量：V = k₀d√(n₁²−n₂²)，b = (n_eff²−n₂²)/(n₁²−n₂²)。V 是归一化频率，b 是归一化传播常数。对对称平板波导（n₂=n₃），TE_m 模的截止条件为 V = mπ/2（m=0,1,2…），因此 V < π/2 时仅存在 TE₀ 基模，这就是平板波导的单模条件。TM 模特征方程形式类似，但界面条件涉及 n₁²、n₂² 的比值，其截止 V 值略大。
二、矩形波导与有效折射率法
矩形波导在 x、y 两个方向都受限，没有严格解析解。有效折射率法（EIM）把二维问题拆成两个一维问题：先沿 x 方向求解一个等效平板波导，得到各阶模式的等效折射率 N_xm；再把 N_xm 当作“芯层折射率”，沿 y 方向求解第二个平板波导，得到最终 n_eff。以宽度 w、高度 h 的条波导为例：
- 沿高度方向（x）取宽度无限大的三层平板，芯层 n₁、厚度 h，包层 n₂、n₃，解出各阶 TE 模的 N_xm。
- 沿宽度方向（y）构造平板：芯层折射率取 N_xm，厚度 w，两侧包层取 n₂（或空气），再解一次特征方程，得到 n_eff,mn。
EIM 的精度在弱导、远离截止时较好，误差通常在 1% 量级；对高折射率差、深亚波长波导误差会增大，此时应以 FEM/FDTD 结果为准。EIM 的价值在于快速给出模式数量与 n_eff 初值，为后续全波仿真提供扫描起点和模式识别依据。
三、模式重叠积分与耦合判据
两个模式场之间的重叠积分定义为：
Γ = |∫∫ E₁*(x,y)·E₂(x,y) dxdy|² / (∫∫|E₁|²dxdy · ∫∫|E₂|²dxdy)Γ 取值 0～1，越接近 1 表示模式越匹配。它直接决定端面耦合效率、定向耦合器耦合系数、光栅耦合器效率以及量子器件中自发辐射因子 β 的估算。例如，光纤基模（模场直径约 10.4 μm）与 220 nm×500 nm SOI 条波导基模的重叠积分通常只有 0.02～0.05，这就是为什么需要倒锥形模斑转换器——通过把波导宽度从 500 nm 渐变到约 100 nm，使模场扩展，Γ 可提升到 0.7 以上。
四、用 FDTD/FEM 计算单模条件与模式场分布
FEM 适合求本征模式：在波导横截面建立二维网格，设置完美匹配层（PML）吸收边界，求解 Maxwell 本征值问题，直接得到 n_eff 与 (E_x,E_y,E_z,H_x,H_y,H_z)。FDTD 则适合宽带激励与传播仿真，通过模式光源注入、监视器记录场分布，再用模式展开得到各模式功率占比。
实际操作流程（以 FEM 为例）：
- 建立横截面几何：Si 芯 500 nm×220 nm，埋氧 2 μm，上包层空气或 SiO₂。
- 材料折射率：Si 用 Sellmeier 或固定 3.48@1550 nm，SiO₂ 取 1.44。
- 网格：芯层内最大单元 ≤ 20 nm，包层可放宽到 100 nm，PML 厚度约 1 μm。
- 求解前若干阶模式，记录 n_eff 与场图。若第二阶模式 n_eff 与基模差小于 1e-3，或第二阶模式被截止（n_eff < 包层折射率），则判定单模。
- 扫描宽度 w 从 300 nm 到 700 nm，观察 TE₀ 与 TE₁ 的 n_eff 曲线交点，交点即单模—多模边界。
判断依据：对 SOI 条波导，TE 单模条件通常为 w < 约 450 nm（h=220 nm）；TM 单模条件更宽。若 n_eff(TE₁) 低于 SiO₂ 折射率 1.44，说明该模式已泄漏，波导实际只导 TE₀。FDTD 中可用模式重叠积分验证：把仿真场与解析基模做 Γ 计算，Γ>0.95 说明激励以基模为主。
五、关键参数与常见问题
| 参数/现象 | 典型值或表现 | 处理建议 |
|---|---|---|
| V 数 | SOI 220 nm 高，约 0.6～1.5 | V<π/2 平板单模；矩形需二维确认 |
| n_eff 基模 | TE₀ 约 2.4～2.8@1550 nm | 与解析/EIM 对比，偏差大查网格与材料 |
| 模式泄漏 | n_eff 接近包层折射率 | 加宽波导或增大限制，检查 PML 距离 |
| 重叠积分低 | Γ<0.1 | 引入倒锥、光栅或透镜模斑转换 |
| FDTD 模式不纯 | 监视器场含高阶分量 | 延长传播距离，用模式展开分离 |
小结：平板波导给出解析本征方程与单模判据 V<π/2；矩形波导用 EIM 把二维拆成两个一维问题快速估算 n_eff；FEM/FDTD 提供精确模式场与单模边界；模式重叠积分 Γ 则把“模式形状”转化为可量化的耦合效率指标。四者配合，才能从理论走到可复现的仿真与优化。