1. 课程导入与FDTD工程概览

光子器件仿真并不是“把结构画出来点一下运行”这么简单。它更像一条从物理目标出发、经过几何与材料建模、数值离散、边界与激励设置、求解、后处理，再回到设计判断的工程链路。时域有限差分法（Finite-Difference Time-Domain，FDTD）在这条链路中承担的是“直接求解麦克斯韦旋度方程”的角色：它不假设模式、不预设传播方向，而是在时间步进中让电磁场在空间网格上自然演化。因此，它特别适合处理宽带响应、强色散、非线性、复杂几何以及需要观察瞬态过程的光子器件问题。

理解FDTD的工程定位，先要把它和另外两类常见方法区分开。一类是频域方法，例如有限元法（FEM）和本征模展开，它们通常在单一频率或一组频率上求解，适合高Q腔、模式分析和精细频域扫描；另一类是光束传播法（BPM），它沿传播方向做近似，适合长距离、弱折射率对比的波导器件。FDTD的优势在于一次时域仿真可通过傅里叶变换得到宽频带响应，且能直接处理反射、透射、辐射和模式转换同时发生的场景。代价是网格必须足够细以分辨最小波长和最小特征尺寸，时间步长受Courant稳定性条件限制，三维问题对内存和计算量要求较高。

从工程对象看，本课程会反复接触几类典型光子器件：硅基波导与弯曲波导、定向耦合器、多模干涉器（MMI）、微环谐振器、光子晶体微腔、光栅耦合器、超表面单元以及片上模式复用/解复用结构。它们的共同点是尺寸在亚微米到数十微米量级，工作波段多在近红外通信波段（例如1550 nm附近）。工程目标通常不是“看场图好不好看”，而是回答几个可量化的问题：在目标波段内插入损耗是多少？回波损耗是否低于阈值？带宽是否覆盖C波段？模式纯度或串扰是否满足链路预算？谐振波长与Q值是否落在设计窗口？这些指标决定了仿真设置中必须包含哪些监视器、哪些边界条件以及多长的收敛时间。

为了让后续章节有统一的参照，这里先约定本课程采用的坐标系与方向定义：以笛卡尔坐标为准，x为横向，y为纵向（或厚度方向，视器件类型而定），z为传播方向。对于波导类器件，默认光沿+z方向传播，输入端口位于z=0附近，输出端口位于z=L附近。反射率、透射率均以端口平面上的净功率通量计算，方向由坡印廷矢量S = E × H的z分量符号决定。若某节改用其他约定，会在该节明确说明并保持图示一致。

一个典型的FDTD工程流程可以概括为：明确物理目标与指标 → 选择仿真维度（2D/3D）与对称性 → 建立几何与材料模型（含色散模型） → 设置网格与收敛性 → 配置边界条件（PML、周期、对称） → 设置激励源（模式源、偶极子、平面波） → 布置监视器（时域场、频域场、功率、模式展开） → 运行并检查收敛 → 后处理提取S参数、损耗、Q值、模式纯度 → 与设计目标对比并迭代。这个顺序不是死板的，但每一步都直接影响结果可信度。例如，PML层数不足会导致虚假反射，网格过粗会让有效折射率偏移，监视器位置太靠近散射体则无法分离入射与反射。

上图的链路中，最容易被低估的是“收敛性与可信度检查”。FDTD的结果不会自动告诉你它是否算对了。工程上至少要做三件事：第一，把网格步长缩小10%～20%再跑一次，观察关键指标变化是否小于允许误差；第二，检查PML区域内的场是否已衰减到可忽略，若PML内仍有明显场分布，说明反射会污染结果；第三，对无源器件做能量守恒检查，即输入功率应约等于输出功率加损耗，残差过大通常意味着监视器位置、模式源或材料拟合有问题。

以硅基微环谐振器为例，一个合理的教学示例设置是：环半径5 μm，波导截面220 nm × 500 nm，埋氧层2 μm，工作波段1500～1600 nm。网格步长取Δx = Δy = Δz = 20 nm，时间步长由Δt ≤ Δx / (c√3)估算，约3.8×10⁻¹⁷ s。PML层数设为10层，模式源从直波导输入端注入，监视器分别放在输入端口、直通端口和下载端口，另在环内放置一个频域点监视器记录谐振谱。运行约20000步后，通过傅里叶变换得到下载端口频谱，读取谐振峰位置与3 dB带宽，再计算Q值。若谐振峰偏移超过2 nm或Q值明显低于设计预期，优先检查材料色散模型和网格色散误差，而不是直接改几何。