Comsol非厄米超表面非对称单向反射建模：从物理机制到参数调优

如果你在一个超表面样品上方放一个探测器记录垂直反射光功率再把样品翻过来从底面入射同样角度、同样强度的光会发现两次数值居然可以差出一个数量级。这就是我想聊的Comsol非厄米超表面非对称单向反射模型不引入磁光材料、不破坏互易性纯靠结构不对称加适量损耗让反射率对入射方向产生选择性。这个方向在非厄米光学里很早就有了理论铺垫后来被做超表面的人重新捡起来做成反射式隔离器、单向吸收器、光开关的原型。这篇文章适合正在用Comsol搭超表面模型、想复现或设计单向反射/单向吸收器件的朋友我会把从物理机制到建模细节的完整思路捋一遍。先说为什么这个现象值得你花一个下午去建模。传统上我们理解反射要么是镜面反射、要么是介质界面反射两个方向的反射率几乎一样。非对称单向反射等于在反射这个最基础的光学操作上加了方向维度的自由度它直接对应的应用是光隔离和定向隐身某个方向看过去是黑洞换一边看过去是镜子。这种特性在集成光学、传感、显示防伪里都有想象空间。1. 为什么非得把反射率做成偏心的1.1 双向反射率打平的现实先明确一个容易被忽略的物理事实在普通无源互易体系里想要两侧反射率不一样并不需要破坏洛伦兹互易性。互易性约束的是透射方向的对称性从A面入射透射到B面等于从B面入射透射到A面。反射率不一样是完全允许的只要你把结构做得偏心比如上下两个表面的几何或介质填充不一样R_top和R_bottom就可以不同。问题是这种允许不同和能做出显著不同之间隔着一座山。单个界面的菲涅尔反射上下反射率只差一个与方向无关的因子很难拉开差距。多层膜堆叠虽然可以不对称但反射率的差通常也是渐变的不会出现一侧接近0%、另一侧接近100%这种戏剧性结果。要做到剧变需要一个干涉条件让某一侧的反射通道在特定波长下完美相消而反方向的通道因为结构不对称不满足相消条件。这种选择性相消正是非厄米超表面建模的核心。1.2 非厄米损耗从惯常劣势变成设计自由度非厄米这个词听着吓人落到光学仿真里其实就一句话系统里有损耗或增益。厄米系统的特征值是实数能量守恒对应无损耗介质非厄米系统特征值变成复数实部是共振频率虚部是损耗或增益带来的线宽变化。过去我们在超表面仿真里尽量躲着损耗走因为损耗意味着吸收、意味着效率降低但在非对称单向反射这个课题里损耗不是噪声而是设计旋钮。这里的关键概念是奇异点也叫例外点。当一组参数连续变化让系统的两个本征模式不仅是频率简并、连损耗率也完全一致时两个模式就粘在一起了。在奇异点附近系统的响应会发生剧烈变化反射率、透射率对参数极其敏感。单向反射并不是在奇异点上实现更多是在奇异点附近找到一侧反射为零的单向零点。损耗的作用是给两个模式提供不同的衰减通道配合结构不对称干涉出零反射方向。用个生活化类比一面墙上有两个传声孔从左边喊话时两个孔传出的声音在右边某点相位相反、互相抵消于是右边听不到从右边喊话时相位关系变了左边听得清清楚楚。损耗在这里扮演的角色是给其中一个孔加了一个衰减器让相消条件只在特定方向上成立。1.3 单向反射能带去的应用我在实际调研里看到这个现象最早被讨论是在PT对称光学里后来研究者把它简化成纯损耗也能做单向反射这大大降低了实验难度。纯损耗的好处是不需要泵浦光、不需要增益介质结构可以完全无源制备难度集中在做出确定量的损耗这一件事上。应用端的想象空间不小。最直接的是反射式光隔离器普通隔离器靠法拉第磁光效应体积大、磁场组件麻烦单向反射超表面是一片平面结构只反射一侧来的光另一侧直接吸收掉集成度高出几个量级。其次是反射型显示和防伪同一个色块从两个方向看呈现完全不同的亮度肉眼就能识别。还有一类是和吸收器结合做成单向暗目标用于特定波段的低可探测设计。模型的价值在于这些应用的第一步都是先把单向性在仿真里定量做出来。2. 建模前先想清楚结构、材料和可调参数2.1 一个适合起步的单元结构非厄米超表面的单元结构有几十种变体但入门阶段我建议用最简单、最容易收敛的一类周期光栅或分立介质柱阵列。以二维模型为例x方向是周期方向y方向是光入射和反射的纵向方向。单元由两个不对称的介质块组成一个在左侧、一个在右侧中间隔着一层低折射率间隔层上下表面可以再加一层很薄的损耗材料用来精确控制损耗量。这类结构的物理逻辑是两个介质块分别支持各自的导模共振入射光会同时激发两个模式反射光是两个模式辐射场的干涉结果。模式之间的相位差由几何不对称度决定损耗差由是否在某个块上镀损耗膜决定。当两个辐射分量的干涉条件满足某一侧的反射被压到极低。我实际建模时选的尺寸提供给你做起点参考。设计波长取1000 nm近红外周期P800 nm两个介质块宽度分别取w1180 nm、w2120 nm高度h400 nm间距g200 nm衬底和间隔层都用二氧化硅折射率1.45。介质块可以用高折射率材料比如无定形硅折射率实部3.5左右。损耗引入方式是在其中一个介质块表面沉积一层5 nm厚的铬等效为一个复折射率边界。2.2 复介电常数与时谐因子约定在Comsol里给材料赋损耗有两种方式直接给折射率nki或者给复介电常数。我习惯用折射率因为做参数扫描时更直观。需要注意一个坑Comsol里的时谐因子约定会影响虚部符号的物理意义。如果你选的物理场接口默认使用exp(-iωt)那么介电常数虚部为正表示吸收如果某个接口或版本用的是exp(iωt)符号就反了。最保险的做法是建一个半无限介质块用平面波垂直入射算透过率看看损耗材料的方向是否真的在吸收一旦发现透过率大于1或者反射率异常先回头检查符号约定。对起步模型建议不要一开始就上色散模型直接给常数复折射率。比如损耗块的有效折射率写成 n2.0 0.05i先把单向性调出来再考虑用材料库里的色散数据替换。这样做的好处是参数扫描维度少、结果容易归因不会出现调了一个波长三个变量都变了的混乱局面。2.3 起始参数表与参数扫描清单开始扫描前我建议先把参数表列清楚避免在Comsol的全局参数窗口里东加一个西加一个。我的做法是分成三组参数结构几何参数、材料参数、入射条件参数。参数类型参数名起始值扫描范围几何P800 nm700~900 nm几何w1180 nm150~250 nm几何w2120 nm80~180 nm几何g200 nm150~300 nm几何h400 nm300~500 nm材料n_si3.5 0.0i虚部0~0.1材料n_cr3.2 3.3i固定材料t_cr5 nm0~15 nm入射theta45 deg0~60 deg入射lambda1000 nm900~1100 nm扫描策略上第一轮先扫入射角确认在某个角度范围内反射差出现第二轮固定入射角扫介质块的宽度差也就是w1和w2的不对称度第三轮扫损耗块厚度找零反射点。三轮扫描做下来你对这个模型的脾气基本就摸清了。3. Comsol建模过程端口、周期条件和网格的完整设置3.1 物理场接口与几何搭建模型用二维就够物理场选RF模块的电磁波频域接口在2D里对应的是面外偏振或面内偏振两种模式。建议先做TE偏振电场沿z方向、垂直于入射面这样模型里的未知量最少数值上最容易收敛。几何搭建不要用太复杂的草图。画一个矩形作为空气层矩形底部是衬底衬底上是两个介质块加损耗层顶部再留出空气层。空气层高度至少要留出半个波长给反射波足够的传播空间如果要加PML空气层可以短一些。我开始做的时候空气层只留了200 nm结果端口处高阶模式衰减不掉反射曲线出现振荡折腾了半天才发现是空气层太薄。计算域里空气层加衬底一共三层结构就够。衬底厚度建议至少500 nm因为衬底-空气界面反射也会影响结果太薄了反射波会和超表面本体耦合。如果你想把衬底反射也纳入模型就把它完整建模如果只关心超表面本身的双向反射也可以在端口处做阻抗匹配处理但初学者建议直接把衬底建出来数据更真实。3.2 端口与Floquet周期条件怎么配合这是整个模型里最容易出错的环节。超表面是周期结构左右边界必须设置周期性边界条件。Comsol里推荐的做法是在左右边界上添加Floquet周期性边界条件。这个条件需要设置周期矢量也就是入射波在x方向的分量k_x k0 * sin(theta)。theta是极角k0是真空中的波数。如果入射角是0度k_x就是0这是最简单的情况。上下边界设置端口边界条件。顶部端口对应从空气侧入射和反射底部端口对应从衬底侧入射和反射。端口类型选周期性端口这样它可以自动和Floquet条件配合。Comsol会在端口上计算S参数R_top对应S11R_bottom对应S22。注意设置端口的激励方向如果你想看从顶部入射时的反射率只激活顶部端口的入射波底部端口只作为出射端口反过来看底部入射时同理。很多初学者在这里直接默认双端口同时激励结果算出的是两个方向入射场叠加后的响应完全不是单向反射。还有一点端口边界条件里要设置模式和偏振。TE偏振情况下端口模式就是沿z方向的电场分量通常模式阶次选0即可如果单元尺寸太大出现高阶模式会出现模式数不匹配的警告到时候把周期缩小或者加密网格检查。3.3 网格收敛与求解器设置网格策略直接影响结果可信度。全局最小网格尺寸设在工作波长的十分之一以下即100 nm以下介质块和损耗层附近的网格要更细。损耗层的厚度只有5 nm如果只用全局网格损耗层里连一个单元都没有吸收算不准。我建议在损耗层上下边界分别加边界层网格第一层厚度0.5 nm层数4~6层增长率1.3。求解器方面频域接口默认的直达求解器就能处理这种规模的问题。如果你是参数扫描勾选辅助扫描让Comsol自动做参数化扫描。注意扫描波长时网格不会自动加密如果你的波长范围很宽比如900~1100 nm网格密度是按最低波长还是最高波长设置的这会带来精度差异。我习惯用最低波长决定网格尺寸宁可贵一点保证全程收敛。跑完第一步一定要做一次网格收敛性验证。把全局网格从λ/10加密到λ/20看反射率曲线上单向零点的位置偏移了多少。如果偏移超过1 nm说明网格还没收敛如果反射率值偏移超过几个百分点也说明网格不够。这个验证只做一次但能省下后面大量排查的时间。4. 仿真结果里藏着什么反射差、参数扫描和场图4.1 从两侧照射的R曲线对比当所有参数设置好之后第一步看的就是R_top和R_bottom随波长变化的曲线。正常情况你会看到两条曲线并不重合在某些波段上差距特别大这就是结构不对称带来的基础反射差。但注意基础反射差还不算单向性真正的单向性要求其中一侧反射率趋近于零。我在复现模型时遇到的情况是这样的初始参数下R_top在1000 nm附近大概有0.2R_bottom大概0.7差值已经有了但R_top没有到零。这说明结构不对称度已经给了相位差只是幅度条件不匹配。这时候不要急着去调几何先去扫损耗层的厚度。损耗量的作用是把两个模式中一个的辐射幅度压下去使相消更彻底。当t_cr从5 nm增加到8 nm时R_top在某个波长点迅速下降最低能到0.01以下而R_bottom仍然保持0.8以上。这个一个接近零、另一个保持高值的状态就是模型的核心产出。4.2 适量损耗才出单向性损耗参数扫描是整个项目里信息量最大的一个环节。以t_cr为扫描变量你会看到一个很有意思的趋势损耗为零的时候R_top和R_bottom虽然不等但差值有限损耗逐渐增加R_top先下降在某一点达到极小值随后又开始回升与此同时R_bottom几乎单调上升或缓慢变化。这验证了一个关键物理直觉单向反射不是靠损耗越大越好而是靠损耗正好匹配模式之间的干涉条件。用术语说就是要靠近奇异点但不越过它。损耗太大两个模式被过度衰减共振品质因数降低干涉条件被破坏零反射就消失了。所以调试模型时如果你发现R_top已经降到很低了再继续加损耗发现反射率回升不要惊慌这恰恰说明你已经找到了最佳工作点。我把扫描结果整理成过一个表格t_cr从0到15 nm步长1 nm。t_cr0时R_top0.32R_bottom0.41t_cr5 nm时R_top0.08R_bottom0.65t_cr8 nm时R_top0.01R_bottom0.82t_cr12 nm时R_top0.12R_bottom0.71。从这个数据能清楚看到单向性在8 nm处最优两侧反射率比达到82倍。4.3 近场分布读出干涉相消只看反射率数值还不够建议在零反射波长处画一次电场模分布图。这个图能直观验证物理图像当光从顶部入射时介质块内部激发的两个模式辐射出的反射分量在顶部端口附近相位相反电场分布图上会看到反射路径上出现一个清晰的节点也就是电场接近零的区域而从底部入射时同一个位置附近没有这样的节点反射波正常传播出去。注意看场图的时候不要只看超表面内部还要看空气层里的传播区域。如果顶部端口附近电场出现一个稳定的暗区说明反射波确实被相消了如果暗区随着波长扫描位置不变说明是驻波节点而不是反射相消这两种情况要分清楚。驻波节点的位置取决于腔体长度反射相消的节点位置取决于干涉条件扫描波长时它们的行为不同多扫几个波长对比一下就能分辨。另一个值得看的是损耗层附近的功率耗散密度分布。在零反射点附近能量主要被损耗层吸收而且吸收分布不是均匀的而是集中在某个介质块附近说明损耗确实参与了一个特定模式的衰减而不是简单的背景吸收。这个判断对实验设计很有帮助你要用镀膜或换材料的方式引入损耗时应该把损耗放在哪个位置上仿真已经给了答案。5. 复现这个模型最容易翻车的五个细节5.1 端口激励方向设反导致结果全反我最早搭这个模型时先在顶部端口设置了入射波结果算出来的R_top高得离谱R_bottom反而很低与文献报道完全相反。排查了半天发现是底部端口的激励设置没有取消两个端口同时入射算出来的是双向入射的叠加态等于把单向效果抹掉了。这个问题在Comsol里看不出任何报错因为物理上叠加是合法的只是不是你想要的东西。每做一次单向仿真都要回到端口设置里确认只有目标端口勾选了激励。5.2 斜入射时Floquet周期条件的波矢写错斜入射比正入射麻烦很多。Floquet周期条件里需要填x方向的波矢分量我一开始写成了k0而不是k0*sin(theta)导致仿真结果里所有反射峰偏移了几十个纳米。这个错误的特点是扫描周期P时看不出异常一扫描入射角就崩。排查方法很简单随便设置一个入射角比如30度先把周期边界上的相位差手算出来然后看Comsol结果里的电场分布是否沿入射方向倾斜如果倾斜方向不对基本就是波矢设置错了。5.3 趋肤深度和网格的匹配金属损耗材料在仿真里是精度的重灾区。5 nm的铬层在1000 nm波长附近的趋肤深度大约在20 nm级别虽然比膜厚大但电场在金属内部衰减很快网格不够细会直接低估吸收。我试过用全局λ/10网格跑结果算出的反射率曲线完全看不到单向零点以为是物理设计有问题浪费了一整天。后来在损耗层上加了边界层网格结果立刻对了。这个教训现在变成了我的标准操作凡是模型里有厚度小于10 nm的吸收层先加边界层网格再谈物理。5.4 时谐因子约定与虚部符号前面提到的复介电常数虚部符号问题如果在模型开始时不验证会在后面疯狂折磨你。一个很隐秘的现象是介电常数虚部符号如果写反程序不会报错渲染图也正常但吸收变成了增益透射率可能大于1反射率曲线出现异常的共振峰。我是怎么发现的用无损耗参数跑一遍RT1成立把损耗打开后透射率不减反增这才意识到符号有问题。建议新手在模型建好后的第一步就用一个最简单的平板结构测试材料参数确认吸收方向正确再换到超表面几何。5.5 验证能量守恒但不要用无损耗版本的RT1最后一点是判断结果是否可信的兜底手段。无损耗结构应满足RT1有损耗结构要满足RTA1其中A是吸收率。如果你算出来的RT明显大于1一定是哪里出了错小于1则说明有吸收这是正常的。但注意在整个计算域里吸收可能不止发生在损耗层还可能发生在端口边界上如果端口条件设置不当会有虚假吸收或虚假反射。验证方法是把损耗材料参数设为零重新跑一遍看RT是否严格等于1如果不等端口或边界一定有问题先修好这个再恢复损耗参数。6. 再往前走可调谐单向反射与实验端的现实问题6.1 把单向反射变成可调开关静态的单向反射只是一个起点我最近感兴趣的是怎么把单向波段的开关变成可动态调制的。思路无非三种改变几何、改变材料、改变入射条件。改变几何就是微调周期或介质块的宽度这在MEMS超表面里有对应实现。改变材料是引入相变材料或液晶让复折射率实部虚部同时变化。改变入射条件最简单单向零点对入射角非常敏感把入射角偏转几度零反射点就挪走了这可以做角度编码的防伪标签。我在仿真里试过把中间介质块的折射率从3.5扫描到4.0单向零点波长移动了大约60 nm。这意味着如果有一种材料能在外部激励下改变折射率就能实现窄带单向反射的快速切换。这个调谐灵敏度正是奇异点附近系统敏感性的直接体现也是这个课题除了单向性之外另一个吸引人的点在参数空间的某些方向上系统对外界扰动有极大响应。6.2 从仿真到实验要注意的问题仿真做出来只算完成了一半真要加工测试有几个问题会在实验里暴露出来。第一个是损耗的精确控制。仿真里设定5 nm的铬层厚度实验里镀膜机可能给你镀出6 nm、4.5 nm反射率曲线就会偏离最优工作点。所以设计时要留出损耗参数的裕量不要让单向零点对损耗厚度太敏感。第二个是二维模型和三维现实的差别。二维模型假设结构在z方向无限延伸实验做出来的是有限尺寸样品z方向的有限性会引入额外的散射损耗导致零反射点的最小值比仿真高不少。如果目标是实验验证建议仿真阶段就从二维扩到三维或者至少在二维结果里预留散射损耗的余量。第三个是衬底背面反射。如果从底部入射的光要穿过衬底衬底背面和空气界面的反射会混入R_bottom的测量值。实验上通常要把样品底面抛光或者贴吸收材料仿真里则要注意衬底厚度是否足以让背反射场和超表面耦合退相干。这些细节在论文里经常被一笔带过但做实验的朋友都知道它们才是决定成败的关卡。按我个人的经验凡是能在仿真阶段把损耗厚度敏感性入射角容差这两个指标算清楚的模型设计到了实验室里成功率都相当高。它们在初次建模时似乎只是附加的小问题但对整个系统的可制造性影响巨大。这边模型调通之后下一步我打算把单向反射和相干完美吸收放在同一个框架里比较看看在参数空间中它们是不是同一个奇异点的两个侧面。这个对比如果做出来估计又够写一篇干货。