陈颖 周健 丁志欣 张敏 朱奇光

1)(燕山大学电气工程学院,测试计量技术与仪器河北省重点实验室,秦皇岛 066004)

2)(燕山大学信息科学与工程学院,河北省特种光纤与光纤传感器重点实验室,秦皇岛 066004)

基于衍射原理和模耦合理论,提出了一种由亚波长介质光栅/金属-电介质-金属(metal-dielectric-metal,MDM)波导/周期性光子晶体组成的复合微纳结构.结合反射角谱深入分析了表面等离子激元的传输特性以及在固定波长下不同入射角时刻形成的双重Fano 共振的产生机理.研究表明,双重Fano 共振是由在亚波长介质光栅/MDM 波导结合的上层结构中产生的独立可调的双离散态分别与在周期性光子晶体中形成的连续态相互耦合形成的.接着定量讨论了结构参数对双重Fano 特性的影响,探究了双重Fano 共振的演变规律.结果表明,改变结构参数可实现双Fano 共振曲线和谐振角度之间的调谐,且在最优条件下,共振A 区FR a和FR b 的品质因数(figure of merit,FOM)可高达460.0 和 4.00×104,共振B 区FR a 和FR b 的FOM 值可高达269.2 和 2.22×104 .该结构可为基于Fano 共振的折射率传感器设计提供有效的理论参考.

1 引言

随着微纳光学的快速发展,光与微纳结构的相互作用增强在机理分析、实际应用等方面的研究越来越深入.而关于亚波长光栅[1,2]、光子晶体[3,4]以及波导[5,6]等光学微纳结构在材料选取以及结构组成等方面的设计也越来越新颖,利用光与不同的微纳结构作用后会产生表面等离子激元(surface plasmon polaritons,SPPs)[7]、导模共振[8,9]、衍射效应[10,11]以及光子禁带[12]等具有实质性差异的光学效应这一特点,在光谱图中动态观察对称或非对称的反射、透射等谱线变化来改进微纳结构的光学性能,以更好地实现对光学器件的调控越来越受到国内外研究人员的关注.

而近年来,大量的研究结果表明,相较于光谱中对称的洛伦兹线型,Fano 共振线型具有尖锐不对称、高敏感性、高品质因子以及相位突变等优势,在提高微纳器件性能的同时,可实用性更强,故如今更偏重于设计可实现Fano 共振的微纳结构.Klimov 等[13]提出了一种由多孔金薄膜与多层膜光子晶体相结合的微纳传感结构,利用光学Tamm态与波导模式相互作用在光谱中形成的非对称Fano 线型探究了结构的性能.Chen 等[14]提出了一种在亚波长光栅和缺陷态一维光子晶体中间夹有石墨烯单层的结构,通过Fano 共振和法布里-珀罗腔(F-P)共振控制石墨烯嵌入结构中光的吸收和透射.Jiang 等[15]提出了一种基于Kretschmann构型且夹有单层MoS2薄膜的多层膜结构.利用平面波导模式之间的耦合产生Fano 共振以应用于生物传感器中,其品质因子(figure of merit,FOM)可达到9630.Zhao 和Yu[16]提出了一种带有凹槽和短节谐振器的金属-绝缘体-金属波导结构,探究了光在传输过程中产生Fano 共振响应光谱的特性,并且实现了波长调谐,其FOM 值可达到3.04×103.以往基于Fano 共振的微纳结构多注重结构设计的新颖以及光学性能的提高,这使得深入研究Fano 共振形成的基础理论以及Fano 共振的演变规律变得尤为重要.

本文提出了一种由亚波长介质光栅/金属-电介质-金属(metal-dielectric-metal,MDM)波导/周期性光子晶体组成的复合微纳结构.结合反射角谱对光的传输特性、双窄带峰以及宽频带的产生机理进行了分析,解释了双重Fano 共振形成的物理机制.在满足波矢匹配的条件下,入射光通过介质光栅-上层金属Ag 薄膜时会产生SPPs 模式,传播到MDM 波导中会产生波导模式共振,SPPs 模式和波导模式共振所提供的不同入射角时刻下的离散态会分别与光子晶体中产生的光子禁带效应所提供的连续态发生近场耦合,从而产生双重Fano共振.接着定量分析了结构参数对双重Fano 共振的影响,进一步探究了双重Fano 共振的演变规律以及最优值.

2 结构模型建立及理论分析

2.1 结构模型建立

本文提出的由亚波长介质光栅/MDM 波导/周期性光子晶体组成的复合结构模型如图1 所示.运用有限元分析软件COMSOL Multiphysics 建立了几何模型,其中亚波长介质光栅层的材料采用折射率为nA=2.58 的 TiO2,光栅凹槽宽度为W=55 nm,高度为dg=46 nm.MDM 波导由上下厚度均为dAg1=dAg2=2 nm 的金属Ag 层和折射率为nb=3.88,高度dz=262 nm 的介质Si 层组成.周期性光子晶体 (A/B)N则由周期性排列的折射率为nA=2.58 的 TiO2和折射率为nB=1.45 的 SiO2组成,周期层数N=3.外界环境介质均为空气,na=1.00.模拟时以一个光栅周期为一个单元,将左右边界设置为周期性边界条件,结构顶端设置入射端口和反射端口,底端设置透射端口,详细探究其角度调制下的光学传输特性.

图1 复合结构模型图Fig.1.Composite structure model diagram.

Ag 的相对介电常数满足Drude-Lorentz 模型[17]:

式中,ω为光的角频率,ε∞=2.4064 为当ω趋于无穷时的介电常数,ωP=2π×2214.6×1012Hz 为等离子体振荡频率,γ=2π×4.8×1012Hz 为碰撞频率,Δ=1.6604 为洛伦兹项的权重系数,Ω=2π×1330.1×1012Hz 为洛伦兹谐振子的强度,Γ=2π×620.7×1012Hz 为振动谱宽.银的复折射率,其中实部nM描述银对光的折射特性,虚部k称为消光系数,描述银对光的吸收特性.

2.2 理论分析

SPPs 是一种由入射光和金属表面自由电子相互耦合形成的高度局域的倏逝波[18],沿着金属表面水平传播且在竖直方向上呈指数衰减.根据电磁场边界条件以及周期性光子晶体的光学特性,选用图2(a)所示的光子禁带的中心波长λ0=632.8 nm的TM 偏振光入射.首先入射光经过光栅衍射效应产生的衍射波会与上层金属Ag 薄膜表面的自由电子耦合,产生等离子体振荡现象.其衍射效应可表示为[19]

图2 结构曲线图 (a)连续态光谱图;(b)连续态角谱图;(c)双离散态角谱图;(d)相位角谱图;(e)双Fano 共振形成角谱图Fig.2.Structural diagram:(a)Continuous spectrum;(b)continuous state angular spectrum diagram;(c)double discrete states angular spectrum diagram;(d)phase angular spectrum diagram;(e)double Fano resonances form angular spectrum diagram.

式中θa为入射角,φa为衍射角,na为光栅凹槽处介质折射率,m为介质光栅衍射级数,λ0为入射光波长.

经过光栅耦合后,入射光波的动量得以增强,补偿了与SPPs 波矢之间的差值,使其在满足波矢匹配条件下可激发SPPs.条件为[20]

式中,k0=2π/λ0为入射光在真空中的波矢,εf为介质光栅的有效介电常数,εAg为金属Ag 薄膜的介电常数.

此时会有大量的磁场能量被束缚在介质光栅-上层金属Ag 薄膜界面附近,使得金属表面的光吸收率增大,反射光能量急剧减少,从而形成窄带峰.与图2(c)反射角谱中入射角θa=40.1°处的peak a相符.

由于MDM 波导中金属Ag 薄膜厚度均满足SPPs 在金属Ag 层中的穿透深度,故首先会有能量传播到MDM 波导中.当介质层厚度可比拟SPPs 的趋肤深度时,M-D 与D-M 交界面上分别激发的SPPs 会相互作用发生耦合共振现象.根据波导理论,可知TM 模式在MDM 波导中的本征方程为[21]

式中,k0为入射光波矢,β为TM 模式的传播常数,εb为电介质层介电常数,dz为电介质层厚度,m0为TM 模的模阶数,φ为SPPs 在电介质与金属交界面发生反射时的相位差.

当入射角θb=48.9°时,φ将满足(4)式中的相位匹配条件,MDM 波导中将产生波导模式共振,形成图2(c)所示反射角谱中的共振峰peak b,磁场主要局域在电介质和下层金属Ag 薄膜交界面处.此时仍会有部分能量渗透到周期性光子晶体中,进而产生如图2(b)所示的光学效应稳定的较宽频带,为Fano 共振提供一个连续态.由图2(d)所示的相位角谱图可知,入射光在θa=40.1°和θb=48.9°时刻传播的相位差均趋近于π,相位相反,满足相位匹配条件,结构间发生了干涉相消.此时独立可调的双离散态会分别与宽连续态发生近场耦合,实现双重Fano 共振,进而导致反射角谱线型发生强烈变化,形成与入射角度θa和θb相对应的具有明显非对称性的双Fano 共振曲线.具体形成过程如图2(e)所示.

3 演变规律分析

探究双重Fano 共振的演变规律时,选取FOM值作为评价的重要指标,其表达式如下[22]:

式中,αFOM为FOM 值;ω为离散态和连续态耦合时的频率;R(ω,n0)和R(ω,n)分别表示为初始状态和折射率变化后的反射率; Δn=n-n0为改变共振区域A,B 处介质折射率值时所对应的折射率差.

3.1 光子晶体周期层数N 对双Fano 特性的影响

基于光子禁带特性,光子晶体可在反射角谱中产生一个较宽的连续态,而光子晶体的周期层数N会直接影响连续态的选频范围和反射率大小,进而影响连续态与在介质光栅/MDM 波导结构中所形成的双窄带的耦合程度.随着层数N的增加,光子禁带反射率会逐渐变高,带隙变窄.光子晶体与介质光栅/MDM 波导结构间的耦合成分变少,耦合能力变弱,双Fano 曲线旁带反射率会逐渐变高,与图3(a)和图3(b)中的角谱曲线变化情况相符.且在图3(b)中,当N=6 时,出现了几乎对称的洛伦兹线型,旁带反射率极高,已不再具有Fano 特性.根据(5)式计算出在不同周期层数的条件下,共振A 区和共振B 区所对应的介质折射率na,nb改变时,FR a 和FR b 曲线的FOM 值.由图3(c)可知,共振A 区和共振B 区处FR a 和FR b 的FOM 值均随N的增加逐渐减少.故为了获得最佳的双Fano 线型及高FOM 值,光子晶体层数被设定为N=3,此时FOM 值最高可以达到4.00×104.

图3 N 对双Fano 特性的影响 (a)N 对连续态曲线的影响;(b)N 对双Fano 曲线的影响;(c)N 对FOM 值的影响Fig.3.Effect of N on characteristics of double Fano resonances:(a)Effect of N on the continuous state curves;(b)effect of N on the double Fano curves;(c)effect of N on FOM values.

3.2 金属Ag 层高度 dAg 对双Fano 特性的影响

在固定入射波长下,金属层的厚度会直接影响双窄带峰与宽频带之间的耦合程度.此时需要考虑SPPs 在金属Ag 层中的趋肤深度,可表示为[23]

式中,k0为入射光波矢,εm为金属相对介电常数,εb为电介质相对介电常数.

根据SPPs 传播特性可知,当金属层厚度越厚时,共振深度会越小,光波渗透能力越弱,与图4(a)所示的双Fano 曲线变化相符.且图4(a)中的双Fano 曲线均发生蓝移,当dAg1=dAg2=12 nm 时,FR a 与FR b 曲线共振深度减小明显.说明渗透的能量极少,与连续态耦合只产生微小的Fano 线型,与理论分析相一致.由(5)式及图4(b)可知,随着dAg的增加,共振A 区和共振B 区的FR a 曲线的FOM 值均逐渐减小,FR b 曲线的FOM 值均先增大后减小.当dAg1=dAg2=2 nm,FOM 值最高可达 4.00×104.此时金属Ag 层厚度正满足(6)式中的趋肤深度δ1.

图4 dAg 对双Fano 特性的影响 (a)dAg 对双Fano 曲线的影响;(b)dAg 对FOM 值的影响Fig.4.Effect of dAg on characteristics of double Fano resonances:(a)Effect of dAg on the double Fano curves;(b)effect of dAg on FOM values.

3.3 光栅凹槽宽度W 对共振A 区Fano 特性的影响

根据等效介质理论,在TM 偏振光的入射下,介质光栅的等效折射率为[24]

式中,nA和na分别为介质光栅高低折射率,f=Λ-W/Λ为介质光栅的占空比,W为介质光栅的凹槽宽度.

由(2)式、(3)式和(7)式可知,凹槽宽度W主要影响θa=40.1°时共振A 区的Fano 线型,故扫描角度范围设定为30.0°—42.5°.当TM 偏振光入射到W不同的介质光栅层时,会具有不同的衍射波、波矢匹配条件以及等效折射率,进而导致窄带共振深度以及位置发生改变.如图5(a)所示,当W增加时,光栅等效折射率会变大,为满足波矢匹配条件,FR a 线型会发生蓝移,且FR a 共振峰波谷反射率有变大趋势,与理论分析相一致.根据(5)式计算出当共振A 区介质折射率na改变时,FR a 曲线的FOM 值,如图5(b)所示.可以看出,随着W的增加,FOM 值逐渐减小,在W=55 nm处FOM 值最高可达460.0.

图5 W 对共振A 区Fano 特性的影响 (a)W 对FR a 曲线的影响;(b)W 对FOM 值的影响Fig.5.Effect of W on Fano characteristics of resonance region A:(a)Effect of W on FR a curve;(b)effect of W on FOM values.

3.4 MDM 波导中介质层高度 dz 对共振B 区Fano 特性的影响

在MDM 波导中,SPPs 在电介质中传播的趋肤深度决定了介质层的最小厚度.可表示为[23]

式中,k0为入射光波矢,εm为金属相对介电常数,εb为电介质相对介电常数.

由(4)式可知,介质层高度dz主要影响在角度扫描区间42.5°—55.0°中存在的共振B 区的Fano曲线.当dz过小时,波导内只存在TM0 模式,共振强度明显减弱,与图6(a)中dz=252 nm 时的FR b曲线相符.当dz过大时,上下金属Ag 层与电介质界面产生的SPPs 的相关系数为0,共振减弱,与图6(a)中dz=280 nm 时的曲线相符.而当dz可与δ2相比拟时,M-D 与D-M 界面所激发的SPPs 将会在向前传播的过程中相互叠加耦合,共振增强,如图6(a)中dz=262 nm 时的红线所示,此时dz已满足(8)式中介质层的最小厚度.当dz由262 nm增加到271 nm 时,FR b 曲线会发生蓝移,波谷反射率逐渐增大,且旁带反射率明显变小,与理论分析相一致.根据(5)式计算出当共振B 区介质折射率nb改变时,FR b 曲线的FOM 值,如图6(b)所示.可以看出,随着dz的增加,FOM 值会逐渐减小,当dz=262 nm 时FOM 值最高可达 2.22×104.

图6 dz 对共振B 区Fano 特性的影响 (a)dz 对FR b 曲线的影响;(b)dz 对FOM 值的影响Fig.6.Effect of dz on Fano characteristics of resonance region B:(a)Effect of dz on FR b curve;(b)effect of dz on FOM values.

综合以上分析,最终选取光子晶体周期层数N=3,凹槽宽度W=55 nm,MDM 波导中金属Ag 层高度dAg1=dAg2=2 nm,电介质层高度dz=262 nm.此时共振A 区FR a 和FR b 曲线的FOM值可高达460.0 和 4.00×104,共振B 区的FR a 和FR b 曲线的FOM 值可达到269.2 和 2.22×104.基于双重Fano 共振对结构参数及介质折射率的高度敏感性,通过观察共振谱线的偏移以及FOM 值的变化趋势,可实现对共振A 区和共振B 区不同折射率范围内待测介质的动态检测.

4 结论

基于角度调制提出了一种由亚波长介质光栅/MDM 波导/周期性光子晶体组成的可实现双重Fano 共振的复合微纳结构.在满足波矢匹配的条件下,入射光经过亚波长介质光栅-上层Ag 薄膜会产生SPPs,传播到MDM 波导中会产生波导模式共振,SPPs 模式和波导模式共振会分别与在光子晶体中产生的光子禁带效应发生近场耦合,从而形成双重Fano 共振.定量分析与双离散态、连续态有关的结构参数对双Fano 共振谱线的影响,阐述了双重Fano 共振的演变规律.研究表明,当周期性光子晶体层数N=3 ,介质光栅凹槽宽度W=55 nm,MDM 波导中金属Ag 层高度dAg1=dAg2=2 nm,电介质层高度dz=262 nm 时,共振A 区FR a 和FR b 的FOM 值可高达460.0 和 4.00×104,共振B 区FR a 和FR b 的FOM 值可高达269.2和 2.22×104.本文的研究结果对设计基于Fano 共振的微纳传感结构有一定的指导意义.