写这篇东西的时候我刚把手头一个双层扭转光子晶体的模型从 COMSOL 里导出来顺手把远场偏振数据绘成了庞加莱球上的轨迹。看着那些数据点在球面上走出平滑的弧线觉得这个思路确实值得分享用扭转光子晶体中的连续谱束缚态BIC来调控远场偏振听起来像是个非常高冷的物理概念但落在 COMSOL 的实操层面其实是一条相当清晰的路线——建模、扫参、寻模、提偏振每一步都有章可循。这篇文章的定位是写给两类人一类是刚接触 BIC 和光子晶体仿真、想在自己体系里复现准 BIC 辐射偏振可控的研究生另一类是已经在做微纳光学仿真、但对 COMSOL 里远场偏振提取流程还不熟悉的工程师。我会把我自己踩过的坑、试过的参数、以及最终稳定可复现的操作路径都整理出来尽量不说空话。1. 为什么是扭转光子晶体BIC 从锁死到可调的物理逻辑1.1 BIC 的三个关键特征连续谱束缚态Bound state in the continuumBIC在光子晶体里的直观理解是一个模式的能量位于辐射连续谱之内却因为对称性或相位条件被锁在结构内部不向外辐射。用仿真的语言说它就是一个拥有无穷大品质因子Q 值、没有远场辐射通道的本征模式。做仿真的人最容易混淆的是BIC 高 Q 模式这个说法。严格讲理想 BIC 的 Q 是无穷大但因为实际数值仿真一定是有限维截断、有限网格、有限边界你永远算不出无穷大的 Q只能看到一个足够大、且随着网格加密不断增大的 Q 值。所以我们在 COMSOL 里找 BIC本质上是在找Q 值对求解精度极度敏感、且空间场分布在某个动量点呈现明显对称性锁定的模式。等真正处理准 BICsymmetry-protected quasi-BIC时Q 值会从无穷大掉到有限值这个掉下来的过程就是偏振调控的窗口。1.2 单层光子晶体的偏振困局单层光子晶体在 Gamma 点支持的 BIC通常是由面内镜像对称性保护的。比如一个椭圆孔硅膜结构TE 模式在某个对称面约束下无法向外辐射远场偏振是零——没有辐射就没有偏振可言。就算你通过改变椭圆率把它变成准 BIC远场的偏振态往往也被钉死在某个线偏振方向自由度很低。为什么因为单层结构的对称性约束太强BIC 的拓扑电荷位于动量空间的某个点上周围的远场偏振呈现涡旋状分布你只能在这个涡旋的局部做微调很难做到想要圆偏振就圆偏振、想要任意方位角线偏振就任意方位角。1.3 扭转的双层结构到底改变了什么扭转光子晶体的核心玩法是当两层同周期光子晶体叠放并相对旋转一个角度 θ 后原本单层结构的面内对称性被协同打破。这个对称性破缺是连续可调的——θ 从 0 变到 30 度对称性的退化程度连续变化因此准 BIC 的 Q 值连续变化远场偏振态也随之连续演化。这个思路和直接改椭圆率不一样的地方在于扭转引入的是一种手性几何层面的对称性破缺它不光打破镜像对称还会让上下两层之间的近场耦合产生一个等效的旋光项。在 COMSOL 的数值结果里你会看到同一个准 BIC 模式的输出偏振场对角 θ 的变化非常敏感灵敏度远大于单纯修改某个孔的几何尺寸。所以扭转结构实现任意偏振态的物理基础是对称性破缺的强度θ和形式双层耦合同时成为调控旋钮你可以在庞加莱球上覆盖从线偏振到圆偏振的一大片区域。2. COMSOL 建模从单层到双层扭转的几何与边界完整设置2.1 三明治结构怎么搭我用的是最经典的三明治构型底层为矩形晶格光子晶体平板厚度通常在几百纳米量级中间是一层低折射率间隔层或者直接是空气顶层是同样的光子晶体平板但绕光轴旋转了 θ 角。这个结构在 COMSOL 里建三维模型并不复杂但有几个几何细节直接影响后续求解。首先是晶格矢量的统一。上下两层必须保证晶格常数一致否则扭转之后根本形不成周期性的超胞。我在建模时习惯用参数化几何设一个全局参数 a0 作为晶格常数底层两个基矢分别定义为 (a0, 0, 0) 和 (0, a0, 0)顶层则在输入几何时直接把第二个晶格矢量旋转 θ例如变成 (-a0 sinθ, a0 cosθ, 0)。很多人直接手动旋转顶层实体后面遇到周期边界条件定义不一致时就会非常痛苦。其次是间隔层厚度。这个参数我建议做成扫描变量。间隔太薄比如小于 50nm两层间的倏逝波耦合过强准 BIC 的 Q 值掉得太快间隔太厚比如超过 300nm两层近似退耦扭转效应几乎消失。我通常扫描 100-250nm 的范围能找到偏振调控最灵敏的工作区间。2.2 Floquet 周期条件的正确姿势顶层扭转之后整个双层结构不再是简单的矩形晶格周期结构。这一点是 COMSOL 建模里最容易踩的坑——你仍然可以在某个超胞里建模但需要想清楚这个超胞的周期是什么。我自己的做法是如果扭转角选的是特殊角比如 30 度上下两层的晶格会构成一个更大的复合超胞这时候周期性边界条件必须定义在这个复合超胞上而不是单独某一层的原胞。如果是任意角比如 17.5 度严格意义上结构不是周期性的只能做近似处理取一个足够大的截断超胞并在边界上施加 Floquet 周期条件来近似模拟无穷大结构。这种做法在文献里很常见所谓 moiré 超胞近似但你要心里有数——超胞越大计算量越大而 BIC 这种高 Q 模式恰恰对超胞截断非常敏感。在 COMSOL 的周期性边界条件里我建议把周期性类型选为 Floquet 周期并显式设置 k 向量的两个分量。做 Gamma 点 BIC 时 kxky0但为了看能带色散和 BIC 拓扑结构最好还是把 k 向量参数化后面扫动量和扫扭转角可以放在同一个研究中。2.3 材料与求解域的取舍材料层面我的标准配置是高折射率介质平板用硅折射率约 3.45损耗先设为 0间隔层和背景是空气n1。如果你关心吸收对 Q 值的影响再单独加一个带虚部的材料扫描但找 BIC 阶段不建议加损耗——否则你分不清 Q 值下降是因为对称性破缺还是材料吸收。求解域方面一定要用完美匹配层PML。虽然周期性边界条件已经处理了面内的无穷延展但上下两个方向的辐射泄漏必须靠 PML 吸收。这个 PML 不要贴得太近结构表面我一般留出至少半个波长距离。边界层要是太近高 Q 模式的倏逝尾巴会被 PML 截断算出来的 Q 值会显著偏小。3. BIC 的搜寻与判定为什么不可能直接看到无穷大 Q3.1 有限网格下的 Q 值上限与解读在 COMSOL 的特征频率研究中Q 值通常通过本征频率的实部除以两倍虚部来获得Q Re(f) / (2|Im(f)|)。理想 BIC 的实频率应该对应一个振动不衰减的解但数值上它的虚部受网格精度限制表现为一个很小但非零的伪值。因此你会发现同一几何下越加密网格、Q 虚部越小这是判断此模式是不是 BIC的有效手段。实际操作时我采用的做法是固定扭转角 θ0即两层完全对准用中等网格先做一轮特征频率扫描找到一组频率靠得很近的模式通常是简并模对然后加密网格两倍看其中某个模式的本征频率虚部是否持续下降。如果虚部比其它模式低两到三个数量级这个模式十有八九就是 BIC。3.2 扭转角扫描与准 BIC Q 值曲线确认 θ0 处的 BIC 后开始做参数化扫描。把 θ 作为扫描参数我从 0.5 度开始间隔 1 度一路跑到 30 度在每个角度重新求解特征频率追踪同一个模式的 Q 值。这里出现一个非常常见的现象随着 θ 增大Q 值快速下降。对于对称性保护的准 BICQ 与 θ 的关系一般近似满足 Q ∝ 1/θ²这来自二次型对称性破缺的耦合理论。你可以在后处理里把 Q 对 θ 画成双对数坐标如果数据点落在斜率接近 -2 的直线上说明这个模式的准 BIC 属性非常干净后面做偏振提取会很舒服。如果 Q 的下降远快于 1/θ²那就要小心了——很可能是两层之间形成了额外的辐射通道比如超胞折叠带来的额外衍射级次。这时候要么减小超胞截断误差要么重新审视 θ 的选取。3.3 简并模分裂与模式认亲扭转角不为零时原本 θ0 处的简并模式对会分裂成两个频率不同的模式。频率差正比于层间耦合强度。在 COMSOL 的模式列表里这两个模式通常紧挨着很容易认错。我的经验是不要靠频率认模式要靠电场分布。BIC 相关的模式在两层平板内部有强烈的对称分布特征E 场在上下层的节面位置高度对称且面外分量很小。你在 COMSOL 结果里打开电场模的切面图重点看 z0 中等截面位置的能量分布。如果是干净的准 BIC能量应该集中在平板内部PML 区域几乎看不到电场延伸。更靠谱的方法是利用全局变量探针在每个模式下算总电能与总磁能准 BIC 模式的总电磁能随 θ 的变化应该相对平滑而伪模式会有明显跳变。4. 远场偏振提取从近场复振幅到庞加莱球4.1 远场投影为什么不能直接用COMSOL 的特征频率解给出的是近场复电场分布 E( x, y, z )。但远场偏振是在无穷远处的辐射特性不能直接把近场的某个点的电场分量拿来当偏振方向。正确的路径是先通过近场数据做平面波展开或者傅里叶变换得到各个辐射方向的远场复振幅再对特定衍射级次提取偏振分量。COMSOL 里做这件事我的习惯是在电磁波、频域接口下额外加一个远场域节点指定一个包含结构的虚拟面通常位于结构上方约一个波长处软件会基于 Stratton-Chu 公式自动把该面上的场向外积分得到远场某方向的电场复分量。4.2 Stokes 参数的具体计算流程拿到某个观测方向例如沿 z 方向的远场电场复振幅 Ex 和 Ey 后偏振态的完整描述需要用 Stokes 参数S0 |Ex|² |Ey|²S1 |Ex|² - |Ey|²S2 2 Re(Ex·Ey*)S3 2 Im(Ex·Ey*)在 COMSOL 的派生值计算中我通常先分别定义远场 Ex、Ey 的实部和虚部再按上面公式做全局表达式计算。注意 Ex 和 Ey 的相位参考要一致——都取远场域节点定义的相位中心否则 S2 和 S3 会因为相位基准不一致而完全没有意义。归一化后用 S1/S0、S2/S0、S3/S0 三个坐标就能把偏振态画在庞加莱球上。我在实践中得到的典型结果是θ 从 0 到 15 度的扫描中准 BIC 模式的 Stokes 矢量从接近 S30线偏振连续转向 S3±1圆偏振S1 和 S2 交替变化说明偏振方位角也在旋转。4.3 观测角度与近场相位中心的影响这里有个容易忽略的细节你提取远场偏振的观测角度必须和 BIC 的辐射方向一致。Gamma 点 BIC 经扭转变成准 BIC 后辐射通常仍然沿近法线方向所以取正 z 方向的远场点基本够用。但有些模式会因超胞折叠辐射到非零角度这时候 Polarization 提取需要同时扫描仰角和方位角找到远场强度最大值对应的方向再在那个方向提取 Stokes 参数。我做了一个小验证在 θ5 度时沿天顶角 0-30 度连续采样远场 Stokes 参数发现 S3 随观测角度变化非常快仅仅偏移 2 度就可能从 0.8 掉到 0.3。因此建议在正式分析偏振前先扫描远场辐射方向图锁定辐射主瓣方向。5. 那些真正折磨人的地方网格、模式混叠与假 BIC5.1 网格收敛性会骗人的高 Q 值特征频率求解时COMSOL 默认网格下也可能算出 Q 值几万甚至几十万的模式很多人第一次看到会以为捡到了 BIC。但等你加密网格Q 反而下降——这多半是网格的色散误差在捣乱不是物理上的 BIC。我的判据是在固定 θ 下加密网格到两倍密度如果 Q 值的变化小于 5%并且模式频率的实部位移小于该模式带宽的 1/10才可以认为收敛。对于 θ0 的 BIC 模式Q 值应该随网格加密持续上升并被截断到 PML 反射水平这才是理想 BIC 的特征。如果你没有做网格收敛性测试后面所有的偏振分析都可能是场分布数值噪声的体现。5.2 模式排序跳跃与连续跟踪策略参数化扫描扭转角时COMSOL 的特征频率列表里模式排序经常跳跃在 θ2 度时的第 5 个模式在 θ4 度时可能跑到第 7 个去了。如果你只是导出某个序号的模式数据画出来的 Q-θ 曲线会乱成一团不连续。解决办法有三个我逐个试过都有效一是使用辅助扫描配合特征频率的临近模式追踪。在 COMSOL 6.1 以上的版本里参数化扫描配置中有一个选项可以基于上一参数步的解作为新的初始猜测确保模式连续性。二是收紧特征频率搜索范围比如锁定在某个频率区间上下 10% 内减少无关模式混进来。三是利用模式库/匹配功能先在某一个 θ 值手动识别目标模式记录其频率和场分布然后用它的空间相关度作为后续角度的匹配依据——这个在 LiveLink 配合 MATLAB 做特别顺手。5.3 PML 反射与假 BIC的干扰PML 不是完美吸收体它是有限层数逼近理想吸收的。高 Q 模式的倏逝场如果延伸到 PML 区域会在边界产生非物理反射形成一种数值引入的高 Q 假 BIC。判断方法很简单画出模式的总电场分布看能量密度在 PML 与物理域边界上是否出现一个小的隆起或驻波状尾巴。如果是把 PML 位置的边界条件从默认换成散射边界 低反射多层或者直接增加 PML 层厚。另外一个很隐蔽的问题是 PML 与周期性边界条件的相互作用。当超胞截断方向与 PML 法线方向不垂直时周期边界会迫使场在 PML 侧向配对导致 PML 反射显著增强。我解决这个问题的办法是用完整的矩形超胞四个侧面做周期条件而上下两侧用 PML同时让 PML 区域本身的网格各向异性足够强。5.4 远场偏振提取前的相位一致性检查前面说过了远场变换基于一个虚拟积分面上的场分布。如果这个面选得太靠近结构倏逝分量还没有衰减完远场计算会把近场项当作辐射项引入误差。我建议积分面距离结构表面至少 200nm 或者一个晶格常数以上同时在频域解中保证该面上的场已经呈现典型的行波相位分布。检查方式在积分面上画电场相位切片如果相位沿面内方向呈线性变化说明已经进入远场区如果相位剧烈振荡赶紧把面向外移。6. 给同样在跑这个仿真的你几条实用建议最后按我的习惯说几个不容易在论文里看到、但直接影响仿真效率的点。第一不要一上来就做全参数三维扫描。先用二维有效折射率法或者简化的多层模型估算准 BIC 出现的频率区间再把三个频段放进 COMSOL 做三维特征频率求解。全参数扫描的计算成本很高而近 BIC 模式隐藏在大量低 Q 模式里全靠网格搜索会浪费大量时间。第二尽量用 COMSOL 与外部脚本的接口做批量后处理。我习惯让 COMSOL 把每个扭转角下的近场数据导出成文本文件再用 MATLAB 做 Stokes 参数计算因为批量导出近场时COMSOL 内置的派生值计算在大批量参数下反而没有脚本灵活。第三时刻记住准 BIC 也是 BIC不要把 Q 值降下来就当成共振态损失。我做这个项目初期很担心扭转以后模式不够 BIC后来才意识到准 BIC 的价值恰恰在有限 Q 和高选择性偏振辐射之间取了一个极好的平衡点。这个体系目前我还在继续推进下一步打算把同样的扭转调控思路移到非周期性 moiré 结构中看看能不能在更大自由度上操控偏振轨道角动量。对 COMSOL 仿真的新手我的建议是先把我上面写的第 2 章模型搭出来然后直接跳到第 5 章的收敛性检查——把这两件事做好了后面提取偏振态就是水到渠成的事。