3. 光场仿真方法与工具链

硅光器件的横向尺寸通常在数百纳米量级，而纵向传播距离可以从几微米（如定向耦合器）到数毫米（如长螺旋延迟线）。这种“横小纵大”的几何特征，决定了没有任何一种光场仿真方法能同时兼顾所有场景。工程上真正需要掌握的，是四种主流方法各自的数学内核、计算代价来源，以及它们分别“吃”哪一类问题。

四种方法的数学内核与代价来源

FDTD（时域有限差分）直接离散麦克斯韦旋度方程，在 Yee 网格上交替更新 E 和 H 分量，时间步长受 CFL 条件约束：Δt ≤ 1/(c·√(1/Δx²+1/Δy²+1/Δz²))。它的优势是宽带一次求解、天然处理非线性与色散、能直接给出时域脉冲响应；代价是三维网格量随尺寸立方增长，且仿真时间与最大传播距离成正比。硅光中典型设置：波导截面网格 Δx=Δy=20 nm，传播方向 Δz=20~40 nm，PML 层 8~12 层，仿真窗口在传播方向留出 2~3 个波长余量。

FDFD（频域有限差分）把场写成 e^{-jβz} 形式，对每个频率求解大型稀疏线性方程组。它天然适合单频或窄带问题，能直接得到模式场与复传播常数，但矩阵规模随三维网格量增长，通常需要迭代求解器配合预条件。硅光中常用于波导模式、微环本征频率、光子晶体带隙的单频精确求解。

EME（本征模展开）把结构沿传播方向切成若干段，每段内用本地本征模展开，段间通过模式重叠积分（S 矩阵）连接。它的计算量主要取决于段数与每段保留的模式数，而不是总长度，因此对长距离、缓变结构极其高效。硅光中典型设置：每段保留 20~40 个模式，段长 0.5~2 μm，适用于锥形波导、MMI、定向耦合器、绝热耦合器。

BPM（光束传播法）基于傍轴近似，假设场主要沿 +z 方向传播且包络缓变，用 Crank-Nicolson 格式逐步推进。它计算极快，但无法处理大角度散射、背向反射和强折射率突变。硅光中主要用于弱耦合、长距离传播的快速估算，如弯曲波导损耗趋势、渐变结构初扫。

硅光器件的典型方法匹配

把方法映射到具体器件，选择逻辑会清晰很多。以 220 nm SOI 平台、1550 nm 波长为背景：

实操：一次 EME 仿真锥形波导的完整设置

以宽度从 0.5 μm 线性过渡到 3 μm、长度 20 μm 的锥形波导为例，说明 EME 的关键输入与结果判读。

# EME 锥形波导仿真参数（示例，非实测）
wavelength     = 1.55 um
taper_length   = 20 um
w_in           = 0.5 um
w_out          = 3.0 um
n_core         = 3.48   # Si @1550nm
n_clad         = 1.44   # SiO2
cell_x         = 4.0 um
cell_y         = 3.0 um
mesh_dx        = 20 nm
mesh_dy        = 20 nm
num_modes      = 30     # 每段保留模式数
segment_length = 0.5 um # 共 40 段
solver         = "EME"
polarization   = "TE"   # Ex 主导

设置完成后，先检查每段的模式有效折射率是否随宽度单调变化。若某段出现模式交叉或简并，说明该段宽度接近模式截止，需要增加模式数或细化段长。随后查看 S 矩阵的透射系数：TE 基模从输入段到输出段的 |S21|² 应接近 1，若低于 0.95，通常是段长过大导致模式失配，把 segment_length 降到 0.25 μm 再试。

方法之间的衔接与交叉验证

实际工程中很少只用一种方法。常见组合是：用 FDFD 提取波导模式作为 EME 或 FDTD 的输入源；用 EME 快速扫描锥形长度和 MMI 宽度，锁定候选尺寸；再用 FDTD 对候选结构做全波验证，确认没有背向反射或高阶模激发。BPM 则适合在参数扫描初期快速排除明显不合理的几何。

交叉验证时，关注三个量是否一致：有效折射率、透射率、模式重叠积分。若 FDFD 与 EME 给出的有效折射率相差超过 1e-3，优先检查网格尺寸和边界条件；若 FDTD 与 EME 的透射率相差超过 0.5 dB，优先检查 EME 的模式数是否足够、FDTD 的 PML 是否反射。

工具链的输入输出约定

四种方法在工具链中并非孤立。一个可复用的约定是：所有方法共享同一套材料折射率数据（n_Si=3.48、n_SiO2=1.44 @1550 nm）和同一坐标系（+z 为传播方向，x 为横向，y 为纵向）。FDFD 输出的模式场保存为 .npz 或 .mat，作为 EME 和 FDTD 的源；EME 输出的 S 矩阵保存为 .dat，供后续链路预算使用；FDTD 输出的时域监视器数据用于提取频谱响应。统一约定后，方法切换时只需替换求解器，几何和材料定义无需重写。

掌握这四种方法的关键，不是记住公式，而是建立“结构特征 → 方法选择 → 参数设置 → 收敛判据”的条件反射。看到长锥形，先想 EME；看到光栅，先想 FDTD；看到模式求解，先想 FDFD；看到快速趋势，先想 BPM。这个判断链条建立起来后，工具链的搭建和切换就会变得自然。