ARTICLE · INTELLIGENCE

战地情报 · 详情页

来自尧图项目组的一线实战观察与深度解析

微小倾斜介质超表面实现BIC与极致手性响应

微小倾斜介质超表面实现BIC与极致手性响应

先聊一个现象:在纳米光子学里,“连续体中的束缚态”(BIC)和“手性光学响应”原本是两条不太相交的技术线。BIC追求极高的品质因子,光被憋在结构里不往外辐射,Q值理论上可以做到无穷大;手性光学则追求左右旋圆偏振光之间的响应差异,最常见的痛点就是手性信号弱、损耗大。把这两个东西结合到同一个介质超表面里,是近几年很多课题组都想啃的硬骨头。而我最近复盘的一个仿真项目,恰好就是围绕这个点展开的——它在标题里叫“微小倾斜造就极致手性:连续体中束缚态的内生手性研究”,说白了就是我们想验证一个非常朴素的想法:能不能只靠一个很小的柱体倾斜角,同时打开准BIC的辐射通道和光学手性通道。

先说这个内容是什么。简单理解,就是在一块周期排列的介质纳米柱超表面上,把原本竖直的纳米柱人为“推歪”一点点,让结构从严格对称变成三维手性。这个微小的倾斜破坏了原本保护BIC的对称性,让BIC退化为准BIC,Q值从理论无穷大掉到几千几万,但与此同时,结构对左右旋圆偏振光产生了显著的透射差异,圆二色性峰值可以做到0.8以上。它能解决什么问题?传统手性超表面要么用金属等离激元,强手性信号靠欧姆损耗换来,Q值很难做上去;要么用非手性结构加斜入射产生所谓的外生手性,终究要依赖光路和入射方向。我们想得到的,是那种不依赖入射方向、纯粹由结构几何决定的“内生手性”,并且希望它和BIC的高Q特性共存。适合谁看?如果你在做超表面、手性传感、偏振光学器件,或者准备从对称性角度设计高Q共振,这篇复盘里的设计思路、仿真参数和踩坑记录应该能帮你省不少时间。

下面我把整个逻辑拆开讲,包括为什么选微小倾斜、参数怎么设、Q因子和手性到底怎么权衡,以及我在仿真和实测准备中碰到的一堆问题。

1. 为什么要同时追求“束缚态”和“手性”

1.1 先理解BIC到底是什么意思

连续体中的束缚态,从名字就能感觉到它很拧巴:束缚态本该是孤立的、低能量的,连续体本该是开放的、可以向外辐射的,但BIC偏偏就藏在这两者之间。在一个周期结构中,某个模式的频率落进了辐射连续体里,按理说它应该很快漏出去,但由于结构的对称性保护,它的向外耦合通道被完全锁死,于是这个模式可以无限长寿,Q因子在理想条件下趋于无穷。

做个不太严谨但好记的类比:你往一个两边都开口的管道里吹一口气,正常情况下气体会从两个口跑掉。但如果管道中间有一个特别巧的“隔板”结构,让气流在某个本征频率上刚好抵消掉向外的泄漏,那气流就能被困在中间剧烈振荡而不跑出去。BIC就是这个“刚好抵消”。超表面里的对称保护BIC一般出现在倒空间的特殊点,最常见的是Γ点,也就是垂直入射点。在这个点上,模式的对称性与所有辐射通道的对称性不匹配,辐射就没了。

实际操作中并不存在理想BIC,因为真实结构总有尺寸偏差、材料吸收、衬底耦合。但只要把对称性保持得很好,Q值做到几万甚至几十万是可能的。问题是,一个完全对称的BIC往往对手性光学毫无用处:左右旋圆偏振在结构里的行为完全相同,你根本分不出哪个强哪个弱。

1.2 手性响应难在哪

光学手性最简单的体现是圆二色性,也就是左手圆偏振和右手圆偏振在透射率、吸收率或反射率上的差异。自然界手性分子本身的手性光学响应非常弱,所以才需要超表面来“放大”。传统思路是把金属纳米结构做成螺旋、G形、Z形这些自带手性的几何体,靠局域等离激元产生强的近场手性热点。

这条路的问题在Q值。等离激元共振的本质是自由电子振荡,金属在可见光和近红外的欧姆损耗非常大,共振峰的半高宽动辄上百纳米,Q值只有十几到几十。对传感来说,窄峰意味着更高的分辨率和更低的检测下限,Q值上不去,很多应用就悬了。介质超表面没有这个麻烦,但介质结构本身又很难同时做到“高Q”和“手性强”这两个目标。于是很多人把目光投向了带BIC背景的介质超表面:如果能在一个准BIC共振上叠加上手性响应,那岂不又窄又圆二色性?

1.3 “内生手性”和“外生手性”不能混用

这个项目里我对“内生”这两个字格外敏感。手性响应的来源分两类。第一类是外生手性,结构本身可以完全不手性,靠斜入射的光路来“制造”一个有效手性方向;这类方案实现起来方便,很多非手性光栅在斜入射下都能看到明显的CD信号,但它本质上不满足器件级应用的要求——对入射角太敏感。第二类是内生手性,结构本身在三维空间里就是手性的,垂直入射照样有CD,不需要靠倾斜光束来补偿。正是这个原因,我们在设计时坚持用垂直圆偏振光做主要判断标准,如果垂直入射拿不到明显CD,那就说明手性不是真正内生的。

但是内生手性介质结构的问题在于:手性几何的引入往往同时破坏了周期结构原本的对称性,BIC也随之消失。你越想手性强,倾斜、扭转、错位这些几何操作就得越大,BIC就退得越厉害。能不能只用一个很小的扰动,既产生可观的手性,又不把Q值摔得太惨?这就是“微小倾斜造就极致手性”这句话背后的核心矛盾。

2. 设计思路:微小倾斜到底打破了什么

2.1 从竖直椭圆柱开始

我们最终敲定的原型结构并不复杂:一个方形晶格介质超表面,每个晶胞里放一根椭圆截面介质柱,柱体高度H,长半轴a,短半轴b,晶格周期P。材料选用晶体硅,考虑的是近红外波段折射率高且吸收低,典型的折射率大约在3.48左右,衬底是石英。

先看最简单的状态:椭圆柱竖直放、长轴方向和x轴呈45度角。这个结构有两个正交的镜像对称面,还保留了一定的旋转对称性。按照群论的语言,BIC模式在Γ点处的对称性与所有平面波辐射通道的对称性完全不匹配,于是它就是一个标准的对称保护BIC。仿真里能看到的情况是:透射谱在这个模式下干干净净,根本没有共振峰,因为模式压根不往外辐射,也就不会在远场留下形状。把45度角里的圆柱换成椭圆柱,CD也是严格为零,左旋和右旋在结构里的响应完全一样。

2.2 “微小倾斜”的具体做法

接下来是关键操作:把整根柱子的对称轴从垂直于衬底的方向,沿某一个面内方向偏转一个很小的角度α。我个人更习惯把它叫做“微倾斜角”,实际几何上等价于柱顶相对柱底有一个微小的横向位移。

这里有一个容易被忽视的点:单是倾斜还不够。假如柱截面是正圆形,你把竖直柱转成斜柱,结构里依然存在镜像对称面,它还不是手性结构。为了让倾斜起到“产生手性”的作用,柱子截面必须是不对称的椭圆,而且椭圆的长轴方向不能恰好落在倾斜方向上。我们把长轴放在和x轴呈45度的位置,倾斜方向为x方向,这样倾斜和截面取向形成一个耦合,两个镜像对称面就都被打破了。

一句话概括几何逻辑:倾斜提供手性方向,椭圆截面提供面内不对称,两者的配合让结构真正成为三维手性体。设置的倾斜角α并不大,通常从零度扫到6度左右,最大的时候也不过8度。这个“微小”是很刻意的,后面我会解释为什么不能一上来就整个15度20度。

2.3 倾斜角为什么不能太大

从理想BIC出发,准BIC的Q因子随扰动参数呈平方反比下降。在这个结构里,扰动参数就是倾斜角α。α从0度变到2度,Q值可能从理想无穷掉到十几万;从2度变到4度,Q值进一步掉到几万;从4度变到8度,Q值就只剩几千。Q值掉得太快,共振峰变得太宽,手性共振带来的窄带特性就被抹掉了。

那倾斜角加大一点是不是至少能让手性更强?也不是。CD信号不会随着α线性一直涨,它有一个明显趋于饱和的过程。原因在于大到一定程度后,手性响应已经主要不是来自共振增强,而是来自结构的平均响应,随着Q值变差,共振增强贡献反而下降。所以实际扫描出来的CD随α曲线通常先快速上升,然后在一个中等倾斜角附近达到峰值,再缓慢下降或基本持平。这就是标题里“微小倾斜”最讲究的地方:要找到那个“既保留较高Q、又使CD最大化”的操作窗口。

3. 仿真计算与参数分析:从扫描到结论

3.1 仿真工具和参数设置

这个项目里的仿真我主要用严格耦合波分析,也就是RCWA。它本身适合周期结构的透射反射计算,把电磁场展开成一系列平面波谐波,和有限元比起来,计算速度优势非常明显,特别适合要做多组角度扫描、多组波长扫描的情况。全波有限元也可以做,但网格稍微密一点,一个结构的参数扫描算下来就非常煎熬。

几何参数我们固定如下:长半轴a=300纳米,短半轴b=150纳米,柱高H=620纳米,晶格周期P=850纳米。工作波段大致在1400到1450纳米的近红外范围,中心共振波长随倾斜角会有一个几个纳米的偏移。介质柱在计算波段内的折射率按常数3.48设置,石英衬底按1.45设置。这里不得不提一个细节:既然要研究内生手性,仿真里的入射光务必设为垂直入射,因为一旦加入斜入射,很多与结构无关的外生手性会混进来,最后CD的归属就说不清楚了。

RCWA仿真本身还要注意截断阶数。我们通常把空间谐波阶数设到21到31阶,在倾斜角增大之后,高空间频率分量会更明显,阶数如果太低,Q因子和CD都可能出现假象。另外,圆偏振光入射需要定义清楚:从光源看向样品为右旋还是左旋。不同文献定义习惯不一样,计算CD前最好先统一,避免后续比较数据时把左右旋搞反。

3.2 透射谱和准BIC的识别

先看倾斜角为零的情况。透射谱在共振波长附近是平坦的,看不到任何峰和谷,这就是BIC的一个特征:模式不辐射,远场无迹可寻。稍微加一个α=2度的倾斜后,透射谱上立刻出现一个明显的窄谷或者不对称峰,这是典型的Fano共振形态。随着波长扫过共振位置,透射率从背景值急剧下降再回升,谷底的透射率可以降到接近零。

α增大到4度、6度,共振峰明显变宽、变浅。这里有一个值得注意的现象:背景透射也发生了一点变化,高频侧的透射率略低于低频侧,这是柱体倾斜后等效极化率出现弱色散带来的结果。把共振峰提取出来做洛伦兹拟合,就能得到共振波长和半高全宽,Q因子等于共振波长除以半高全宽。对α=2度,Q值大概在八万左右;α=4度掉到接近两万;α=6度只剩五千上下。这个下降趋势基本符合准BIC的平方反比规律,实际数据拟合下来,指数接近负1.9,和理论上负2的预期偏差不大。

3.3 圆二色性数值结果

接着算圆二色性。分别用左旋圆偏振和右旋圆偏振光正入射,计算透射谱,然后定义CD=(T_LCP-T_RCP)/(T_LCP+T_RCP)。如果定义的偏振方向和你实验中选用的一致,正负号无所谓,关键看绝对值和谱线形态。仿真发现,在共振波长附近CD出现一个非常尖锐的峰。α=2度时CD峰值已经到0.5左右;α=4度时峰值逼近0.72;α=6度时峰值大约0.8到0.85,之后继续增大α,CD基本不再上升。与此同时,共振峰宽度从几个纳米展宽到十几个纳米,整体手性品质变差。

从设计角度看,α=4度左右是非常好的折中点:Q值还有两万上下,CD超过0.7,共振线宽仍然很窄,放在传感场景里既能保证波长分辨率,又能提供强圆二色性。下面这张表是我在这个结构上得到的典型数据,注意具体数值会随几何参数和材料色散细节变化,但趋势是稳定的:

倾斜角α(度)共振波长(纳米)Q因子CD峰值
0不显示峰趋于无穷0
2~1426.3~830000.51
4~1427.1~210000.73
6~1428.5~57000.82
8~1430.2~24000.84

3.4 参数扫描里最容易踩的坑

做参数扫描时,最容易被忽视的是网格离散化。RCWA本身是傅里叶级数展开,但如果你把结构模型离散成台阶来近似倾斜表面,台阶的层数会影响Q值收敛。我一开始用比较粗的8层台阶,α=4度的Q算出来只有6000,等加密到32层,Q直接涨到两万以上。这个差距不是物理差异,纯粹是数值离散问题。所以碰到Q值异常低时,第一反应要检查台阶数或网格精度,而不是怀疑物理机制。

另一个坑是扫描波长范围不够密。准BIC的共振峰非常窄,α=2度时半高全宽只有零点零几个纳米。如果波长步长取得太粗,比如0.5纳米一步,很容易把峰尖整个略过去,得到假的低Q值。我一般会先粗扫定位,再在共振位置附近加密步长到0.01到0.02纳米,这样才能可靠地提取Q因子。

4. 物理机制分析:手性从哪儿来,Q和CD如何共存

4.1 模式对称性与远场偏振

仿真算完了,要回答一个更深的问题:为什么一个很小的倾斜能带来这么强的手性响应?关键在模式耦合的对称性条件。竖直椭圆柱状态下,BIC模式具有确定的镜像对称性,它不会和任何平面波辐射通道耦合。加入了倾斜之后,结构镜像对称性被破坏,原本被禁闭的模式获得了和特定辐射通道耦合的能力,于是BIC变成准BIC,透射谱出现共振。

有趣的是,这个耦合并不是同时均等地打开左旋和右旋辐射通道。由于手性几何的存在,准BIC模式在和圆偏振平面波耦合时,左旋和右旋的耦合强度明显不同。你可以把模式理解成一只“歪着身子”的振子:它更愿意把一个手性的光波发射出去,而不太愿意发射另一个手性。这就是CD峰的来源。

从远场偏振拓扑的角度看,完全对称的BIC在动量空间里对应一个偏振奇点,周围的远场偏振方向沿奇点绕一圈会累积一个拓扑荷。倾斜操作不会直接把奇点抹掉,而是让它移动、变形,同时使周围的偏振状态从线偏振逐渐过渡到椭圆偏振。离奇点越近,椭圆度越强,手性选择性也就越凸显。这正好解释了为什么CD峰总是和准BIC峰重合,而不是出现在频带边缘的某个不相关位置。

4.2 近场手性密度的验证

远场圆二色性只是一个侧面。为了确认手性真的来自模式本身,我通常还会提取共振波长处的近场分布,计算光学手性密度场。光学手性密度的定义和电场磁场都有关系,它的正负反映了近场局域手性状态。在竖直对称结构里,左旋和右旋入射得到的近场手性密度是正负对称的,总体积分为零;加入倾斜之后,近场手性密度的空间分布不再对称,在某些热点区域出现单一符号的优势分布,而且这个优势区往往集中在椭圆柱的两个顶端附近。

近场结果的另一个价值是解释传感潜力。虽然远场CD已经够大了,但对普通手性分子来说,真正产生相互作用的是近场手性密度。仿真里可以在柱体周围几个纳米范围内看到手性密度的增强倍数,局部增强可以达到几十倍甚至上百倍,这意味着把待测分子铺在超表面表面上时,圆二色性信号的增强会更可观。当然,近场提取对网格密度要求很高,网格不够时手性密度峰值会被严重低估。

4.3 怎么证明它是“内生手性”

我给自己定的标准有两条。第一条,垂直入射下必须有明显CD信号,如果只有斜入射下才有,就不能叫内生。第二条,把结构替换成它的镜像异构体,CD谱应该整个反转符号。所谓镜像异构体,很简单,把倾斜方向从+x改成-x,或者把椭圆长轴从+45度改成-45度,二选一即可。仿真结果显示,改成镜像结构之后,共振波长和Q值几乎不变,但CD峰的正负完全反转,这基本可以确定手性贡献来自结构几何本身。

还有一个更苛刻的测试:把左右旋的定义保持不变,改变入射面的方位角。如果CD随方位角旋转而旋转,那说明结构中可能还有残留的面内取向效应,需要排查晶格对称性。正儿八经的内生手性应该由结构本身的三维空间构型决定,和光在样品上的方位角没有关系。

5. 加工工艺和实验测量里的现实问题

5.1 倾斜柱结构怎么加工

仿真做得再漂亮,最后要回到“能不能做出来”。倾斜介质柱不是传统的竖直纳米柱,常规的电子束曝光加干法刻蚀没法直接获得一个可控倾斜角。我们盘过几种备选方案,目前最可行的是多层套刻:先沉积一层介质,光刻出椭圆形柱底;然后二次沉积介质,再套刻一层椭圆,并让这一层相对底层有一个横向偏移,横向偏移量除以单层厚度就等于倾斜角。叠足够多层之后,柱子轮廓就呈现出一个近似的倾斜形状。

这个工艺问题必须提前想清楚,因为套刻的对准精度直接决定“微小倾斜”能不能实现。如果偏移量设计值是20纳米,而套刻误差本身也有10纳米,那做出来的样品可能左歪右歪没个准,CD信号就会五花八门。所以在结构参数上,我把倾斜角折中到4度左右也有工艺上的考虑:4度对应的横向偏移在几十纳米量级,比2度更容易对准,又不需要太大偏移,免得Q掉得太凶。

5.2 尺寸公差的影响

介质超表面的共振位置对几何参数极其敏感。我做了简单的公差扫描,结果如下表:

误差来源典型误差范围对共振波长的影响对CD的影响
柱高变化±10纳米波长偏移约6到8纳米中低
椭圆长短轴变化±5纳米波长偏移约3纳米中等
倾斜角变化±0.5度波长偏移约2纳米较大,CD峰值抖动
晶格周期变化±1纳米波长偏移约0.5纳米较小

从表里能看出,倾斜角误差虽然对波长影响不算大,但对CD峰值影响很敏感。换句话说,如果实际结构并不整齐,不同区域倾斜角略有差异,宏观测量到的CD峰会被不均匀展宽。所以测量用的光斑尺寸不能太大,最好控制在几十微米以内,并且要选样品中质量均匀的区域。

5.3 测量需要注意的几个细节

实验测量端,我用的是超连续光源加偏振态发生器,先产生线偏振光,再通过液晶延迟片切换成左旋或右旋圆偏振。测量顺序上,要频繁交替采集左旋和右旋的透射谱,避免光源功率漂移造成的人为CD。如果光源本身存在偏振相关的不稳定性,可以把左右旋的测量顺序随机化,减少系统误差。

样品背面必须贴吸收层或者磨毛标定石英基底,否则衬底背面的反射会引入干扰条纹,使窄共振峰上叠加周期性振荡。这个问题在近红外特别烦人,因为石英透明,前后表面反射很容易形成法布里-珀罗干涉。我试过用专门的楔形石英片作衬底来抑制这种驻波效应,效果不错。

另外要记得扣除参考信号。可以先用没有图案的同一批衬底材料测一组左右旋透射率,再测超表面样品的透射率,两者相除得到归一化透射谱。这里的“参考”和“样品”表面反射率如果不同,会直接带来一个随波长缓变的CD背景,虽然不会影响窄峰的判读,但会拉低峰值数值的可信度。

6. 常见问题与排查技巧实录

6.1 为什么Q值总比理论预期低很多

这是我在仿真阶段碰到最多的一个问题。除了前面说的网格离散和波长扫描步长,另一个常见原因是忽略了材料吸收。晶体硅在近红外吸收系数很低,但不是零。如果把硅的复折射率虚部设置成0,Q可以算到十几万;一旦加入一个很小的虚部,比如虚部0.001,Q值就可能被限制到两万左右。也就是说,当结构本征辐射损耗已经很小的时候,材料吸收损耗就会成为主导,Q因子不再受倾斜角控制。这本身是物理事实,不是bug,但会影响对设计窗口的判断。

6.2 CD谱左右不对称

有时候你会看到CD峰旁边出现一个负值瓣,或者CD谱本身不是单调变化。这里有两种情况。第一种是数值噪声,RCWA在共振波长附近收敛变慢,圆偏振基底下尤其明显,可以增加动量空间的截断阶数再算一遍。第二种是模式重叠,倾斜角较大时,超表面中可能并存另一个偶极模式,它和准BIC模式在频率上靠得很近,两者干涉形成额外的圆二色性边带。这时候可以通过近场分布判断模式是否混杂,再考虑把周期或椭圆尺寸稍微调开一点。

6.3 换了材料参数之后手性峰消失

介电超表面很挑折射率。石英衬底上的硅柱折射率大约3.48左右,但如果换成折射率较低的介质,比如氮化镓或二氧化钛,共振模式往往不再集中在柱体内部,而是大量泄漏到周围环境中,模式与晶格BIC的耦合变弱,手性准BIC就不那么容易出现了。此类问题通常要先回顾材料色散,再看柱高和周期的比值是否还维持在高折射率波导的谐振条件附近。

6.4 近场手性热点的位置和CD峰对不上

近场手性密度最大的地方,往往并不全在柱体表面,也可能在衬底交界处或两个相邻柱之间的间隙里。如果你做分子传感实验,分子只能接触顶端和侧壁,那些藏在衬底内部的热点其实贡献不了。我建议做近场平面取样时,至少取三个高度层:柱顶上方5纳米、柱体中部外壁、柱底与衬底交界面。这样才能看清真正可用的手性热点在哪一层。

6.5 仿真里CD很高,实验却测不到

仿真和实验差距最大的是结构形貌。电子束曝光出来的椭圆柱,顶端不一定平整,侧壁会有一点锥度;多层套刻做出来的“倾斜”也不会是理想的线性倾斜,而更像阶梯状。侧壁的微粗糙度会造成额外的散射损耗,使准BIC的线宽变大,CD峰变矮。处理办法是仿真里加入侧壁倾角模型,用几个台阶近似,然后看CD对侧壁角度的敏感程度。如果仿真显示CD在侧壁偏角2度以内还能保住,说明工艺窗口还能接受,否则设计就要往回退。

6.6 快速自查清单

  • 垂直入射是否确认无误?有没有外生手性混入?
  • 左右旋定义是否一致?介质折射率虚部是否设置为合理值?
  • 波长步长能否分辨半高全宽?
  • RCWA阶数或有限元网格是否收敛?
  • 镜像结构是否真的实现了CD符号反转?
  • 实验里有没有扣除衬底反射和背向干涉?

这套清单每次查一遍,基本能把大多数“假手性”和“假Q”问题挡在门外。

7. 一些后续能直接做的扩展

这个结构不只是停留在固定波段。倾斜角确定后,Q和CD的平衡点可以通过调节柱高和周期来平移,设计灵活性其实比我预想的大。比如在同一块样品上做两个不同倾斜角区域,就能形成一个“高Q参考区”和一个“强手性区”,用来对照分析传感实验背后的机制。

如果考虑主动调控,可以在介质柱周围填充相变材料,或者把衬底换成热光系数较大的材料。温度一变,共振波长移动,CD峰也随之移动,于是就有了热可调的手性开关。倾斜角本身虽然加工完就固定了,但外部环境仍然能够提供额外的调控维度。这个方向我们后面也在做一些初步仿真,目前看响应速度足够快,线宽窄的优势也会保留下来。

我个人在实际操作里的体会是,BIC和手性的组合并不需要一个多复杂的结构,微小倾斜之所以能带来极致手性,本质上是因为它在“对称性禁止辐射”和“不对称产生手性”之间找到了一个精确的平衡点。做这类仿真时,不要一上来就追求超高Q,也不要把CD指标定得太激进,先把几何扰动的平方反比规律摸透,后面所有优化都会顺很多。如果你也想往下深挖,建议先复现一个最小模型,只扫倾斜角这一个参数,把Q和CD的关系图做出来,你的直觉会立刻清晰起来。最后再分享一个小技巧:做CD提取时,把左旋和右旋透射率画在同一张坐标图里,仔细观察两个峰之间是否有个微小横移,那个横移方向就是结构手性的朝向,它往往比单纯看CD峰值更能暴露设计中的对称性问题。

RELATED READING

延伸阅读

更多一线实战笔记与深度复盘,助您持续精进