2. 电磁仿真核心方法

在集成量子光子器件的设计中，电磁仿真的任务是求解麦克斯韦方程组在具体几何结构与材料分布下的解，从而得到器件的模式分布、有效折射率、传输损耗、耦合效率、色散与场局域等关键参数。不同仿真方法本质上是对麦克斯韦方程组做了不同的离散化与近似处理，因此它们在适用结构、计算资源、精度与可扩展性上各有取舍。理解这三种方法的数学出发点，是正确选择工具并判断结果可信度的前提。

三种方法的物理与数学出发点

时域有限差分法（FDTD）直接在时域离散麦克斯韦旋度方程。它将空间划分为 Yee 元胞，电场分量位于元胞棱边中点，磁场分量位于面心，通过交替更新 E 与 H 实现时间推进。其核心优势是：一次时域仿真配合傅里叶变换即可获得宽频带响应，天然适合处理色散材料、非线性效应与超短脉冲传播。稳定性受 Courant 条件约束，即时间步长需满足 Δt ≤ 1/(c√(1/Δx²+1/Δy²+1/Δz²))，网格越细时间步越小，三维仿真代价迅速上升。FDTD 适合光子晶体微腔、波导色散、超表面、瞬态响应等宽带问题。

有限元法（FEM）从频域波动方程出发，将计算域剖分为四面体或三角形单元，在每个单元内用基函数展开场量，通过伽辽金加权残量法形成稀疏矩阵方程求解。FEM 的最大优势是能贴合任意复杂曲面与不规则几何，并可通过自适应网格在场变化剧烈处加密，因此特别适合高折射率对比、弯曲波导、微环谐振腔、铌酸锂调制器等需要精确模式与频域响应的结构。其代价是求解大型稀疏矩阵需要较多内存，三维频域扫描多个波长时计算量较大。

光束传播法（BPM）基于傍轴近似，假设光场主要沿传播方向 z 缓慢变化，将波动方程化为沿 z 的一阶抛物型方程，用步进方式逐段推进。BPM 计算极快、内存占用低，非常适合长距离、弱折射率对比、缓慢变化的波导与定向耦合器。但它忽略背向反射，对高折射率对比、强弯曲、谐振结构或需要精确反射谱的场景误差明显，通常用于初步设计与模式演化快速评估。

如何根据器件类型选择方法

选择方法时，先问三个问题：结构是否具有高折射率对比与强谐振？是否需要宽频带或瞬态信息？计算域在传播方向是否远大于横向尺寸？若结构是微环、光子晶体微腔、超表面单元，且需要精确谐振波长与 Q 值，优先 FEM 或 FDTD；若需要一次获得宽谱透射/反射，FDTD 更直接；若结构是数毫米长的弱限制波导或耦合器，且只需快速评估模式演化与耦合长度，BPM 是效率最高的选择。实际设计中常采用“BPM 初筛 + FEM/FDTD 精算”的组合策略。

Lumerical 与 COMSOL 的初步操作

Lumerical 以 FDTD Solutions 和 MODE Solutions 为核心。以硅波导定向耦合器为例，典型流程为：设定材料折射率（Si 3.48、SiO₂ 1.44 @1550 nm），绘制两条宽 500 nm、间距 200 nm 的波导，设置 FDTD 区域与网格（dx=dy=20 nm，dz=20 nm），在输入端放置模式光源，在输出端放置模式扩展监视器，运行后通过 S 参数提取耦合长度与串扰。MODE Solutions 则用 EME 或 FDE 求解器快速获得波导模式与有效折射率。

COMSOL 以 Wave Optics 模块的频域有限元为主。典型流程为：选择电磁波频域接口，绘制波导截面或三维结构，设置材料折射率，添加端口边界（Port）与散射边界（Scattering Boundary Condition），用自由四面体网格剖分并在波导芯区加密，设置参数化扫描波长（如 1500–1600 nm，步长 1 nm），求解后通过 S 参数与场分布图判断传输与反射。COMSOL 的优势在于多物理场耦合，例如电光调制器可同时求解电场与光场。

结果判读与常见误区

仿真结果的可信度取决于边界条件、网格收敛与材料模型。FDTD 中若监视器未完全覆盖模式，S 参数会偏低；FEM 中若端口模式数不足，高阶级数会被截断；BPM 中若忽略反射，谐振结构结果不可用。判断收敛的实用做法是：固定其他条件，仅加密网格或缩小时间步，观察有效折射率、Q 值或耦合效率的变化是否趋于稳定。对于量子光子器件，还需关注自发辐射因子、模式体积与色散，这些参数对网格与材料色散模型尤为敏感。

掌握这三种方法的边界，比熟练点击软件按钮更重要。实际项目中，先用解析或 BPM 估算耦合长度与模式尺寸，再用 FEM 精算谐振与色散，最后用 FDTD 验证宽带响应，是兼顾效率与精度的常见路径。