模拟电磁波和周期性结构..docxVIP

模拟电磁波和周期性结构我们经常想要模拟入射到周期性结构中的电磁波（光、微波），例如衍射光栅、超材料，或频率选择表面。这可以使用 COMSOL 产品库中的 RF 或波动光学模块来完成。两个模块都提供了 Floquet 周期性边界条件和周期性端口，并将反射和透射衍射级作为入射角和波长的函数进行计算。本博客将介绍这类分析背后的概念，并将介绍这类问题的设定方法。场景首先，让我们来考虑代表周期性重复晶胞的自由空间平行六面体，平面波沿一个角度入射到其上，如下图所示。入射波矢??在全局坐标系中有三个分量，其大小分别是：?,??，和?。该问题可通过在域侧面使用周期性边界条件，并在顶部和底部使用端口边界条件来模拟。该问题设定最复杂的地方是定义入射波和出射波的方向和偏振。定义波方向虽然 COMSOL 软件非常灵活，支持对基矢坐标系的任意定义，但在本博客中，我们将选定一个坐标系并始终使用它。入射光的方向由两个角度定义，??和?；以及两个矢量，??为模拟空间中向外指向的法向，??为入射面的一个矢量。我们这里所选择的约定将??与全局?x?轴对齐，将??与全局?z?轴对齐。因此，入射波的波矢与全局?z?轴之间的夹角为?，即入射仰角，其中??，?，意味着垂直入射。入射波矢和全局?x?轴之间的夹角为入射方位角?，位于范围??之内。 由该定义，??和??的正值表示波沿着 x 轴 和 y 轴正向传播。如要使用入射方向的上述定义，我们需要指定??矢量。这可以通过选定一个周期性端口参考点来完成，它必须是入射端口的角点之一。软件使用从该点出发的面内边来定义两个矢量，和?，且?。在下图中，我们可以看到满足这一条件的四组??和?。因此，俯瞰?z?轴及平面时，入射面端口上的周期性端口参考点应为?x-y?平面左下角的点。通过选择此点，?矢量变得与全局?x?轴对齐。既然选择了周期性端口参考点而在入射面定义了??和??，那么模拟域中出射面的端口也必须定义。法向矢量??指向相反的方向，因此必须调整所选择的周期性端口参考点。四个角点都无法提供与入射面矢量对齐的??和?，因此我们必须选择这四个点之一，并调整??和??的定义。通过在出射面选择与入射面所选定点完全相反的周期性端口参考点，并将??旋转?，?的方向被旋转到?，指向与入射面??相反的方向。由于这种旋转，模拟域中出射面的??和??的正负号发生了转换。下一步，考虑一个代表介电半空间的模拟域，在入射和出射端口面之间存在折射率差异，这会使波方向发生改变，如下图所示。根据斯涅耳定律，我们知道折射角为。这使我们可以计算出射端口处的波矢方向。此外，请注意，即使有额外的介电层夹在两个半空间之间，这种关系也仍然成立。总结一下，要定义通过一个晶胞的平面波方向，我们首先需要选择两个点，即周期性端口参考点，它们在入射面和出射面的位置完全相反。这些点定义了矢量??和?。因此，入射面的??和??可以相对于全局坐标系定义。在出射面上，方向角变为：?和?。定义偏振入射平面波的偏振可以是二者之一，即电场或磁场与?x-y?平面平行。所有其他偏振，例如圆形或椭圆形，都可以由这二者的线性组合建立。下图显示了?，且磁场与?x-y?平面平行的情况。当??时，全局坐标系中的入射和出射端口的磁场大小是?。由于光束旋转使?，磁场大小变为?。对于正交偏振，可以用类似方法定义入射面的电场大小。在出射端口，x-y?平面中的场分量可以用相同的方式定义。到目前为止，我们已经看到了如何定义传播经过介电界面周围晶胞的平面波的方向与偏振。对这类问题，您可以看到案例库中?Fresnel 方程模型的结果与解析解一致。定义衍射级接下来，让我们来检查一下把周期性结构引入模拟域之后会发生什么变化。考虑一个入射到下图周期性结构中、且??的平面波。如果波长相对光栅间距足够短，就可能会存在一个或多个衍射级。要理解这些衍射级，我们必须观察由矢量??和??所定义的平面，以及由矢量??和??所定义的平面。首先，沿着由??和??所定义平面的法向观察，我们能看到零阶透射传输模式的存在，其方向如前所述按照斯涅尔定律定义。零阶反射分量也存在。结构中也可能存在一些光的吸收，不过未在图片中显示。下图仅显示了零阶透射传输模式。间距??是由矢量??和??所定义平面中的周期性。对于足够短的波长，也可以有更高阶的衍射模式。当??时，下图显示了这些情况。这些模式存在的条件是：for:?…其中??时的情况可以简化为斯涅耳定律。当??时，如果光程差等于真空中波长的整数倍，则发生相长干涉，??阶衍射光束的衍射角为?。注意，正负??阶的个数无需相同。下一步，我们沿着由矢量??和??所定义的平面观察。也就是说，我们将视角绕?z?轴旋转，直到入射波矢看上去沿表面的法向。进入此平面的衍射标注为??阶光束。注意该平面上的周期性间距??不同，且正负??阶