3. 仿真区域与边界条件设置

FDTD把连续麦克斯韦方程离散到有限计算网格上，因此必须回答两个问题：计算窗口取多大，窗口边缘的电磁场如何“离开”而不反射回来。边界条件选错，往往不会立刻报错，却会以谐振峰偏移、透射率虚高、模式泄漏或场图出现驻波纹路等形式污染结果。本节围绕四类边界——PML吸收边界、周期性边界、对称边界、Bloch边界——说明它们的物理含义、适用对象、参数设置与结果判读。

一、PML：把出射波“吸走”而不反射

PML（Perfectly Matched Layer，完美匹配层）是在仿真区域外侧附加的一层各向异性有耗介质。其核心思想是通过坐标拉伸，使该层与内部介质的波阻抗在界面处严格匹配，理论上任意频率、任意入射角的波进入PML后都不产生反射，并在层内被指数衰减吸收。工程上常用UPML或CPML（卷积PML）实现，后者对斜入射和倏逝波吸收更好，是多数商用FDTD求解器的默认选择。

PML的关键参数有三个：层数、反射系数和吸收剖面阶数。层数通常取8～16层，层数不足会导致残余反射；反射系数目标值一般设为1e-5～1e-8，越小吸收越强但数值代价略增；阶数常用3阶多项式剖面。以1550 nm波段的硅波导为例，若网格步长Δx=20 nm，PML取12层即240 nm厚度，目标反射系数1e-6，通常可将边界反射压到-60 dB以下。判断PML是否足够，可做一次“PML收敛测试”：固定结构，仅改变PML层数（如8、12、16），观察透射谱或Q值变化；当结果变化小于1%时，认为PML已收敛。

需要注意，PML不能紧贴强散射体或高折射率对比结构。若PML距离器件过近，倏逝场尚未衰减就被吸收，会引入非物理损耗。一般建议PML内边界距器件至少半个波长，对高Q微腔应留出更多余量。

二、周期性边界：用周期结构代表无限阵列

当器件在某一方向无限重复，例如光子晶体、超表面、光栅或波导阵列，可用周期性边界（Periodic Boundary）只仿真一个周期单元。其数学条件是场在相对边界上满足 E(x+P)=E(x)、H(x+P)=H(x)，P为周期。这样，一个单元的计算结果就等价于无限周期结构。

周期性边界必须成对使用，且方向与周期方向一致。例如一维光栅沿x方向周期排列，则x方向两侧设为Periodic，y方向仍用PML吸收出射波。若周期结构还伴随斜入射，则需改用Bloch边界（见下节）。设置时要注意：周期长度P必须与结构实际周期严格一致，否则会引入相位失配，导致布拉格峰位置偏移。

三、对称边界：利用结构对称性减半计算量

对称边界（Symmetric/Antisymmetric Boundary）利用器件的几何与激励对称性，把计算区域缩减为原来的1/2、1/4甚至1/8。其依据是：在对称面上，场的法向分量与切向分量具有确定的奇偶性。若结构关于x=0对称，且激励为偶模（如TE基模的Ey为偶函数），则对称面上切向磁场为零，可设为对称边界（PMC，理想磁导体）；若激励为奇模，则切向电场为零，可设为反对称边界（PEC，理想电导体）。

使用对称边界的前提是：几何对称、材料对称、激励也对称。若激励偏离对称面或存在不对称扰动，对称边界会强制错误对称性，结果不可信。实际设置时，先判断模式奇偶性，再选择PEC或PMC；对波导模式，可用模式求解器确认主分量的对称性。正确使用后，内存和仿真时间可近似按缩减比例下降，对三维大结构尤其有价值。

四、Bloch边界：周期结构斜入射的通用工具

Bloch边界是周期性边界的推广，适用于斜入射或周期结构中存在相位差的情形。其条件为 E(x+P)=E(x)·exp(-i kx P)，其中kx为布洛赫波矢。当kx=0时退化为普通周期性边界。Bloch边界广泛用于超表面、光子晶体能带计算和光栅衍射效率分析。

设置Bloch边界时，需要给定入射角θ和方位角φ，求解器据此计算kx、ky。例如超表面周期P=600 nm，入射角θ=30°，则kx=(2π/λ)·sinθ。若同时存在两个周期方向，需在两个方向分别设置Bloch相位。判断结果是否合理，可检查：改变周期单元数（如1个与3个单元）时，反射/透射谱是否一致；若不一致，说明Bloch相位或周期设置有问题。

五、四类边界的对比与选用

下表汇总四类边界的物理条件、典型对象与关键参数，便于快速选型。

六、边界设置的实际检查清单

为避免边界误用，建议在提交仿真前逐项核对：

• PML是否远离强散射体？内边界距器件是否≥λ/2？
• 周期性边界是否成对出现，周期P是否与结构严格一致？
• 对称边界是否与激励奇偶性匹配？是否存在破坏对称的扰动？
• Bloch相位是否由入射角和周期正确计算？是否做了单元数收敛？
• 是否完成PML层数收敛测试，结果变化是否小于1%？

下面用一张信号链示意图总结四类边界在FDTD计算窗口中的位置关系与作用方向。图中内部为器件区，四周为边界层；PML负责吸收出射波，Periodic/Bloch负责周期延拓，Symmetric/Antisymmetric负责对称面截断。

最后强调一个常见误区：不要为了省计算量而把PML设得很薄，也不要在非周期方向误用周期性边界。边界条件的正确性应通过收敛测试和物理合理性双重验证——透射率不超过1、能量守恒、Q值随PML层数稳定，才是可信结果。