COMSOL-RF模块计算光学器件透射率的方法探讨

COMSOL-RF 模块计算光学器件透射率的方法探讨

透/反射率的计算在电磁波研究中非常常见，计算结果的准确性与材料参数定义，边界条件的选择，网格剖分有十分紧密的关系。以下是关于电磁波透/反射率计算问题的经验整理，如有错漏欢迎指正和补充。需要计算透/反射率的器件通常可分为几种类型：

1. 波导器件

如各类波导分路器，光纤Bragg 光栅，其入射端及出射端都满足波导模式。

当入射及出射端波导满足端口（Port ）内置结构（同轴/矩形），可直接选择内置的波导类型，如RF 案例库中的H 弯波导（h_bend_waveguide）及环形器（lossy_circulator）案例。

当波导结构与内置类型不同时，需要首先通过模场分析计算出波导模式，通过Port 边界的Numeric 类型耦合到频域分析中，作为入射条件。如V3.5a 及V4中的波导适配器

（Wave_adapter）案例，以及在V4.2a 中更新的dielectric_slab_waveguide案例。 波导常常支持多个模式，为了保证作为频域分析边界条件的模场分布是正确的，可以先进行边界模式分析，设定查找多个模式，根据模场分布从结果中找到作为入射条件的模式对应的模指数，然后在进行整个模型分析时，把此模指数作为参考值（Search for modes

around:），查找模式数（Desired number ofmodes:）设定为1。以此保证入射条件正确。

对于以上两种情况，Port 边界内置的S 参数可计算出透/反射率，其中S11对应端口1的振幅反射率，S21对应从端口1至端口2的振幅透射率，以此类推。

2. 周期性散射体

如金属纳米天线阵列，光栅，光子晶体，在一或两个维度上具有周期性。

在RF 模块中，完美电/磁导体（PEC/PMC）是完全反射边界，散射边界（SBC ）、端口（Port ）边界仅对某些角度或分布的光波透明，其他角度的光波均会有一定程度的反射，而PML 如果设置恰当可以保证各角度入射波均被吸收。可以想像，如果散射场在边界上有反射，最终计算出的透射场及反射场会受到影响。边界的选择十分重要。

此类结构，可用周期性边界条件，或是根据电/磁场的对称性用PMC/PEC边界进行简化，仅对重复单元进行模拟。目前的解决方案主要有两种：

a). 入射及出射端采用完美吸收层PML

当入射和出射端均设置为PML 时，怎样定义光源？

在V3.5a 版本中，可以通过Port 边界内部一致对作为入射条件，在入射端和出射端进行能流积分来计算透射率及反射率。典型案例是Grating 。

在V4版本中，内部一致对方法不可行

（http://www.comsol.com/community/forums/general/thread/11030/）

，光源可通过背景场定义。透射功率可通过出射端总场能流积分算出，而反射功率可通过入射端散射场能流积分算出。

如果所研究的结构在入射端和出射端是同一种介质，背景场可直接定义为平面波。但是当入射与出射端处于两种介质中时，比如一个石英板与空气界面上排列着金属颗粒，电磁波从空气入射到界面上，采用一个单独的平面波作为背景场时，会在出射端的PML 边界上出现不合理的反射这时需要根据Fresnel 公式定义出符合界面反射的场分布，或是添加计算背景场的步骤，见

（http://www.comsol.com/community/forums/general/thread/16715/， http://i.cntech.com.cn/forum/post/show.html?tid=8444）。

针对非均一介质的情形，采用首先计算背景场并结合PML 的方法，复现了V4中的

Plasmonic Wire Grating案例，当网格最大尺寸为波长的1/20时，与原模型Port 边界计算结果误差

b). 多重Port 边界方法

V4模型库案例中的Plasmonic Wire Grating，根据一维光栅的衍射光集中于零级及正负一级衍射角度，在同一个边界上设置多个Port 吸收出射光。但是当散射体较为复杂，衍射光可能会在许多角度上存在较强的分布，如晶体在X 射线下的Bragg 衍射，这时设置多重Port 实现吸收并不现实。

3. 单独散射体

如带孔的金属板。与周期性结构不同之处是外围区域都需要设置为PML ，如案例库模型Radar Cross Section中的铝船散射问题。

以上是关于边界条件选择及入射条件的定义，PML 参数通常情况选择默认值，外侧采用SBC 边界，距离散射体应足够远，从计算结果中的场分布可以判断PML 的吸收效果是否充分，如果没有充分吸收，需要修改参数保证反射足够低。

对于包含金属散射体，或是含增益介质的情况，由于COMSOL 软件对于电磁场振幅位相因子的定义方式，当介电常数或折射率虚部为负值时，对应于损耗介质，反之对应于增益介质，如果材料属性设置错误也会对透射率计算产生不必要的误差。

在模拟金属散射体或共振腔结构时的网格剖分十分重要。由于在发生SPP 共振时，金属表面会出现场增强现象，谐振腔处于共振状态同样存在这种情况，那么边界上需要足够的网格以准确的描述场强的指数衰减，例如银材料在可见光波段趋肤深度约为20nm ，这时需要使趋肤深度以内的网格尺寸远小于趋肤深度，约在nm 量级，必要时需采用边界层网格。如同要看清一幅细节丰富的图片，需要足够的分辨率。

本文中提到的计算方法典型模型可通过中仿社区下载，欢迎您的建议与讨论。