在电磁波仿真中，我们经常遇到的一种情况是：计算均匀介电基板表面刻蚀结构的光散射。本文将介绍一种在 COMSOL Multiphysics® 软件中模拟这种情况的方法。
考虑一个小型光散射体：介电基板表面的金制半圆柱体。我们研究仅存在单个散射体的情况。假设一束准直光线，其直径远大于散射体，所以几乎可以近似为一个平面波。平面波以任意的角度和偏振态入射到半圆柱体结构。假设介电基板在负 Z 轴方向无限延伸。
由于模拟的是单个散射体，可以将其表征为对背景场的一个微扰，因此使用散射场公式是合理的。这个公式需要输入一个背景电场，代表没有散射体时麦克斯韦方程的解。也就是说，以频域的形式求解麦克斯韦方程组：
但是，将总电场写成以下形式：
我们求解的不是总场，而是相对场（常被称为散射场）。背景电场是在没有散射体的情况下，有关域上麦克斯韦方程组的解。对于自由空间中的物体，背景场只是一个平面波，比如 RF 模块的基准示例：计算理想导体球的雷达散射截面对比了使用这种方法计算的数值解与解析解。波动光学模块中也有一个类似的示例，它计算了金纳米球的光散射。可能有人会尝试用同样的方法求解半无限空间上的介电质，但这会导致背景场出错。
对于半无限无损介电基板上的散射体，背景场方程必须包括界面的反射和折射。一种方法是基于单独的分析计算场，并将该场用作背景场。COMSOL 案例库中的基板上的颗粒散射教程模型对这种方法进行了演示。但在本文的示例中，我们直接将解析解作为入射到半无限介电基板的背景场输入。
菲涅耳方程与斯涅耳定律共同描述了平面光波入射到具有不同折射率的两种材料界面时的反射和传输。这些方程首先定义了入射平面，即由垂直于两种材料交界面的法向量和入射平面波的波向量描述的平面。入射电场被分解为一个完全垂直于该平面的分量，称为 s 极化或 TE 波，以及一个完全平行于该平面的分量，称为 p 极化或 TM 波。例如，圆偏振光是相同量级的 TE 波和 TM 波的总和，但彼此相位相差 90°。
在入射平面内，可以使用斯涅耳定律，将入射光()和透射光(
)相对于表面法线的角度与两种材料的折射率
和
相关联:
菲涅耳方程用于计算 TE 和 TM 极化的反射和透射系数：
接下来，必须将这些极化分量从入射平面转换回全局直角坐标，以得到入射和透射光束的 k 向量。按照惯例，交界平面是 xy 平面，波长为 的入射光束在负 z 方向传播，将入射角
定义为与 z 轴的夹角，角度
定义为从负 x 轴开始绕 z 轴旋转的角度，如上图所示。因此，入射光束的 k 向量可定义为：
透射光束的 k 向量定义为：
反射光束的 k 向量 与入射光束相似，但 z 分量的符号相反。
已知由 TE 和 TM 极化分量组成的入射光束 和
，入射光束的电场分量为：
对于反射分量 c：
因此，入射侧的总背景场为 ，另一侧的背景场为
，其分量为：
在 COMSOL Multiphysics® 软件中，这些表达式可以在不同域作为一组变量输入，并作为背景场定义。
点击底部“阅读原文”，进入 COMSOL 案例库下载应用上述方法建立的模型。该模型计算并绘制了一个金制半圆柱散射体中的电场大小和损耗。
虽然这种输入分析背景场的方法需要在模型设置方面多做一点工作，但它比基板上的颗粒散射案例模型运行得更快，后者需要先计算背景场。后一个模型方法适用于解析解求解比较困难，甚至不可能求解的情况。因此，对于均匀介电基板，前一种方法可能较好。
