基于锚形谐振腔的等离子体波导特性研究
作者：文阅期刊网 来源：文阅编辑中心 日期：2022-10-10 09:14人气：
摘 要：基于表面等离子体激元在金属-介质-金属结构中优良的传输特性，设计了一种由直波导和锚形谐振腔组成的波导滤波器。仿真分析了波导滤波器传输特性和电场分布随结构参数和谐振腔内介质折射率变化规律。结果表明，该锚形谐振腔最佳结构时滤波器半峰全宽低至8nm，品质因数高达121.9。利用共振波长与结构参数变化规律，设计了光通信波长窗口的窄带带阻滤波器。根据SPPs对谐振腔介质折射率敏感的特性，发现透射谱线随折射率变大而发生红移。本文提出的基于锚形谐振腔的等离子体MIM波导滤波器为设计特定波长的窄带带阻滤波器提供了新思路，也为基于谱线红移特性设计的介质折射率传感器提供了技术支持。
关键词:锚形谐振腔;表面等离子激元;有限元方法;滤波特性;传感特性;
Characteristics of Plasma Waveguide Based on Anchored Cavity
CHEN Zhiyuan DONG Xingfa SUN Hao
School of Electronic and Information Engineering, Suzhou University of Science and Technology
Abstract：Surface plasmon is a kind of electromagnetic wave propagating along the interface between metal and dielectric. SPPs have the property of confining electromagnetic wave energy to half wavelength region, which can break the diffraction limit in traditional optics, and are widely used in the design of sub wavelength optical devices. Based on the excellent transmission characteristics of surface plasmon in metal dielectric metal structure, a waveguide filter composed of straight waveguide and anchor resonator is designed and the transmission characteristics of light in anchor resonator waveguide are studied by finite element method. Using COMSOL multiphysics software to simulate the waveguide, on the one hand, adjust the structural parameters of the anchor resonator θ, H and R, and the factors causing the red shift of the spectral line are obtained. Using MATLAB to fit the simulation data, the structural parameters affecting the resonance wavelength λr are obtained. To explore the influence of structural parameters θ, H=300nm and R=150nm were fixed. When using mode 2 to design the filter, the included angle θ In the range of 105°-150°, the fitting result is λrmode2=2.8θ+ 220.5; When using mode 3, the included angle θ should be controlled within the range of 30°-75°, and the fitting result is λrmode3=9.8θ+225; When using mode 4, the included angle θ is controlled within the range of 75°-90°, and the fitting result is λrmode4=4θ+1030. To explore the influence of structural parameter H, θ=75°,R=150nm were fixed. When using mode 2 to design the filter, only when H=350nm meets the filtering requirements. In mode 3, the strong electric field is almost all distributed in the arc structure and the change of H does not affect Leff, The shifting of λmode3 is only related to θ. When θ=75°, R=150nm, resonance wavelength λmode3 is always 975nm. In mode 4, the fitting result is: λrmode4=3.42H+298.5. To explore the influence of structural parameter R, θ= 75°, H=150nm were fixed. With the change of R from 120nm to 180nm, various modes have different degrees of red shift of spectral lines, of which the red shift in mode 3 is obvious, which is consistent with the analysis results of electric field distribution. In mode 2, only R=160nm can obtain a band stop filter with a resonance wavelength of 620nm. In mode 3, When R is from 150nm to 160nm, the fitting result is λrmode3=7R-75; In mode 4, When R is from 120nm to 180nm, the fitting result is λrmode4=2.71R+916.4. The function relationship can be used to realize narrow-band filtering at specific wavelengths, which provides a new structure for designing narrow-band bandstop filters with tunable wavelengths. Taking the communication wavelength as an example, the corresponding narrowband bandstop filter is designed and its performance index is analyzed. It is found that the ΔFWHM of the anchor resonator is the lowest at 9nm and the Q value is the highest at 94.5. The performance index of the anchor resonator is better than that of other simple resonator filters. On the other hand, by changing the refractive index of the medium in the resonator, it is found that the transmission spectrum of the waveguide shifts red with the increase of the refractive index. Through simulation and data fitting, the sensor sensitivity of the anchor resonator waveguide is better than that of the L-shaped resonator filter, which is more suitable for manufacturing various micro biological solution concentration sensors. In conclusion, a new type of anchor resonator is proposed in this paper. Based on this resonator, a wavelength tunable narrow band stop filter is designed, which provides a new idea for designing narrow band stop filters with specific wavelengths. Based on the red shift phenomenon of the spectral line caused by the change of refractive index in the resonator, a new structure biological solution concentration sensor is designed, which provides technical support for the design of a small-scale micro structure medium refractive index sensor.
Keyword：
Anchor resonator; Surface plasmon polaritons; Finite element method; Filtering characteristics; Sensing characteristics;
0 引言
表面等离子体激元(Surface Plasmon Polaritons，SPPs)是一种沿着金属和电介质界面传播的电磁波。SPPs具有将电磁波能量约束在半个波长区域的性质[1]，可以突破传统光学中的衍射极限，被广泛应用于亚波长光学器件设计[2,3]。为激发SPPs，波导一般设计为介质-金属-介质(Insulator-Metal-Insulator,IMI)或金属-介质-金属(Metal-Insulator-Metal，MIM)结构。IMI波导损耗较低，但约束光传播的能力较弱。与之对应的，MIM结构具有较宽的频谱范围，在亚波长支持高群速模式，而且可实现长距离传播中精确的光调控[4]，被广泛应用于纳米量级先进光学器件的设计，例如纳米级相干光源中的纳米腔[5]、矢量全息的双原子超表面[6]、T型带支节分束器[7]、定向耦合器[8]与红外窄谱增强传感器[9]等。
滤波器中的谐振腔具有波长选择性，利用SPPs与谐振腔耦合共振可以实现光的波长调控[10]。近年来，很多学者通过设计谐振器结构，制作了各种高性能滤波器，例如L形谐振腔滤波器[11]、方形凹环结构谐振器[12]与拱型谐振腔滤波器[13]等。
本文设计了一种基于锚形谐振腔的MIM波导滤波器，对该结构进行仿真，分析其传输特性。仿真结果表明，随着弧形结构角度θ、矩形结构高度H和弧形结构半径R的变化，共振波长处的滤波效果变化非常明显，且能够通过调整参数拟合数据得到对应关系精确控制滤波效果。若在锚形谐振腔内充入不同介质，由于SPPs对谐振腔内介质折射率敏感，可用于设计折射率传感器。
1 原理与数据仿真
1.1 模型建立与理论分析
本文设计的锚形谐振腔MIM结构滤波器示意图如图1，此滤波器由轴对称的弧形结构和垂直于波导的矩形结构组成，金属层为金属银。为使得谐振腔内只存在TM模式[13]，将波导宽度和谐振腔宽度(包括弧形结构和矩形结构)固定为d=50nm[14]，直波导与锚形谐振腔的耦合距离固定为ω=10nm。矩形结构高度为H，弧形结构半径为R，与波导垂线夹角为θ。波导和谐振腔(假定填充空气)介电常数εwaveguide=εcavity=1，金属银的介电常数采用Drude模型[15]计算
εAg=ε∞−ω2pω2+iωγ(1)
式中，ε∞为无穷频率处的介电常数，ωp为等离子体共振频率，γ为电子碰撞的阻尼频率，ω为真空中入射光频率。计算时ε∞=3.7，ωp=1.38×1016Hz，γ=2.37×1013Hz。
波导左端为光入射端口，光功率为Pin，右端为光出射端口，光功率为Pout。使用透射系数T=Pin/Pout，评价光进入波导后的传输特性。在波导结构中，横向磁场传播模式的电磁场方程为[16]
tanh(κd)=−κ(pαc+qαs)κ2+pαcqαs(2)
式中，κ和d为波导传播的电磁波波矢和波导宽度。本文设计的波导由同一种金属组成，则p=q,αc=αs。p=εair/εAg,αc=[k02(εair-εAg)+κ]1/2。自由空间电磁波波矢k0=2π/λ0。计算得到波矢κ后，计算波导有效折射率neff=[εAg+(κ/k0)2]1/2,波导中传播的SPPs波长为λSPPs=λ0/Re(neff)。
1.2 仿真设置及结果
使用COMSOL Multiphysics仿真时设置上下表面为完美匹配吸收层以减少边界反射的影响。网格剖分时使用物理场控制网格，选择极细化以保证网格最大为波长的1/5。研究时需要在一个较宽的波长域进行遍历，精度和计算量的平衡显得尤为重要。经过反复测试，波长域定为400nm~1800nm，遍历步长为5nm时得到较为清晰的传输特性曲线。
引入品质因数Q的概念来评价滤波器的优劣，品质因数Q=λr/ΔFWHM[12]。其中，ΔFWHM为传输频谱的半峰全宽。当H=300nm，θ=75°时锚形谐振腔滤波器传输特性如图2(a)。此结构参数下，波导存在4种共振模式如表1。
图2(b)是相同尺度下的水平放置矩形谐振腔和L形谐振腔透射谱线。对比图2(a)可以发现锚形滤波器谐振波长处的透射率更低，半峰全宽更窄，品质因数更高，滤波效果更好。为进一步研究锚形谐振腔内的共振模式，计算4种模式下SPPs波在波导中的电场分布如图3。当SPPs耦合到锚形谐振腔，在谐振腔内形成稳定的驻波需要满足共振条件
Δϕ=βm⋅2(H+Rθ)+ϕr=2mπ(3)
式中，φr为锚形谐振腔右侧面和上面的光束反射的反射相位；m为在锚形谐振腔形成驻波的波腹数；βm为腔内形成的m阶共振模式对应的SPPs的传播常数。在模式1中m没有正整数解，可见模式1中并没有在锚形谐振腔内形成稳定的驻波，SPPs与锚形谐振腔耦合后依然保持电场强度周期性变化，出现透射谱低谷是因为出射端口恰好位于电场强度谷值。模式2，3，4均满足共振条件，但电场分布不同。模式2中电场集中在弧形结构和矩形结构中，模式3中电场集中在弧形结构内，模式4中电场集中在矩形结构内。对比拱形谐振腔[13],本文结构在保证滤波效果的同时易于优化滤波器特征参数并且没有增加制造工艺难度。
2 结构参数影响分析
2.1 结构参数θ对传输特性的影响
为探究结构参数对波导传输特性的影响，进一步优化滤波结构，对结构参数θ进行分析，这里引入共振条件方程来解释谱线移动，共振条件为[17]
Leff=Nλspp=N[λrRe(neff)],N=1,2,3...(4)
式中，Leff是锚形谐振腔的有效长度，λr是谐振腔的共振波长，N表示为模式数，neff是锚形谐振腔的有效折射率。当H=300nm，R=150nm时，角度θ对透射谱线的影响如图4。根据电场分布，模式1不属于滤波模式，模式2、3、4随夹角θ的增大，透射光谱发生红移并且透射低谷透射率出现变化。本文将透射低谷透射率大于0.1，滤波器明显失去滤波效果的情况排除，绘制了共振波长随夹角θ变化关系如图5。使用模式2设计滤波器时，夹角θ应该控制在105°~150°之间，拟合结果为λrmode2=2.8θ+220.5；使用模式3时，夹角θ应该控制在30°~75°之间，拟合结果为λrmode3=9.83θ+225；使用模式4时，夹角θ应该控制在75°~90°之间，拟合结果为λrmode4=4θ+1030。拟合结果中函数斜率kmode2<kmode4<kmode3,与电场分布关系一致。
2.2 结构参数H对传输特性的影响
在θ=75°,R=150nm的情况下，研究结构参数H对传输特性的影响，结果如图6所示。由于模式2中强电场分布在弧形结构和矩形结构中，模式4中强电场大部分分布在矩形结构内，这两种模式随着H的增大，Leff随之增大，对应的共振波长λr发生红移。使用模式2设计滤波器时，只有当H=350nm时满足滤波要求。模式3中强电场几乎全部分布在弧形结构内，H的变化不影响Leff，λmode3的偏移仅与θ有关。当θ=75°，R=150nm时，共振波长λmode3始终为975nm。模式4中拟合结果为λrmode4=3.42H+298.5，可以借此设计共振波长为675nm、975nm和1150nm~1665nm的窄带带阻滤波器。
2.3 结构参数R对传输特性的影响
在θ=75°，H=150nm的情况下，探究结构参数R对传输特性的影响，结果如图7(a)所示。发现随着R从120nm变化至180nm，各模式均出现不同程度的谱线红移，其中模式3中红移现象最明显，与电场分布的分析结果一致。剔除明显不满足滤波要求的情况得到图7(b)。可以发现在模式2中，只有R=160nm可以得到共振波长为620nm的带阻滤波器。在模式3中，当R=150nm~160nm时，拟合结果为λrmode3=7R-75；在模式4中，当R=120nm~180nm，拟合结果为λrmode4=2.71R+916.4。根据上述拟合结果，可以设计波长为620nm(模式2)、975~1045nm(模式3)、1235~1400nm(模式4)的窄带带阻滤波器。
本节谈论了特定情况下的结构参数变化导致的传输曲线偏移，阐述了导致谱线红移的因素。综合考虑结构参数θ、H、R，利用拟合结果，可以设计特定波长的高性能窄带带阻滤波器。
3 锚形谐振腔的应用
3.1通信窗口下的窄带带阻设计
将基于锚形谐振腔结构特性设计通信波长下的窄带带阻滤波器，以多模光纤的三个窗口为例进行设计。第一窗口使用模式3进行滤波，由图5得到当θ1=64°，H1=300nm，R1=150nm时满足要求。第二窗口和第三窗口使用模式4进行滤波，由图6(b)得到当θ2=75°，H2=293nm，R2=150nm和θ3=75°，H3=366nm，R3=150nm时达到对应的滤波效果，仿真结果如图8。仿真证明拟合公式精确描述了共振波长λr变化规律如表2。设计的三种通信窗口滤波器性能优越，与方形凹环结构MIM滤波器[9]相比(ΔFWHMmin=60nm，Qmax=17.07)性能有明显提升，并且对一种窗口滤波时不会影响其他两种窗口的正常传输，实现了光通信中的通信窗口选择调控。
3.2 折射率传感特性
研究共振波长λr对锚形谐振腔中介质折射率n的敏感特性，仿真得到不同折射率下透射光谱如图9(a)。图9(a)中发现随着介质折射率n增大，各模式共振波长均发生红移且透射波谷发生改变。模式4下共振波谷透射率T始终小于0.1且共振波长红移现象最为明显。传感器灵敏度是衡量传感效果的依据，灵敏度S=Δλr/Δn，其中Δλr为共振波长变化量，Δn为折射率变化量。9(b)中数据拟合得到`S2=600nm∙RIU-1,λrmode2=600n+7.5,`S3=933.3nm∙RIU-1,λrmode3=935n+41,`S4=1316.7nm∙RIU-1,λrmode4=1315n+14。MIM波导中的SPPs耦合对锚形谐振腔介质折射率敏感，可使用该特性制作传感设备。与L形谐振腔的MIM波导结构(Smax=1166.7nm∙RIU-1)相比，灵敏度有提升，更适合用于设计介质折射率传感器。
4 结论
设计了一种基于MIM波导结构的锚形谐振腔滤波器，仿真分析表示，锚形谐振腔的结构参数θ、H和R会影响其滤波效果，实现了特定波长(通信窗口)下的窄带带阻滤波。在锚形谐振腔内充入不同折射率的介质，折射率增大引起了透射光谱红移，可为利用锚形谐振腔设计简单结构的介质折射率传感器提供理论依据。
参考文献
[1] BARNES W L, DEREUX A, EBBESEN T W. Surface plasmon subwavelength optics[J]. Nature, 2003, 4, 24(6950): 824-830.
[2] LIU Liu, HAN Zhanghua,HE Sailing. Novel surface plasmon waveguide for high integration[J]. Optics Express, 2005, 13(17): 6645-6650.
[3] DENG Zilan, DENG Junhong, ZHUANG Xin, et al. Facile metagrating holograms with broadband and extreme angle tolerance.[J]. Light Science & Applications, 2018, 7:78.
[4] YAN Yunfei, ZHANG Guanmao, QIAO Litao, et al. Design on the convex ring MIM structure filter based on surface plasmon polaritons[J]. Acta Photonica Sinica, 2019, 48(2): 0223002.
[5] DENG Zilan, DONG Jianwen. Lasing in plasmon-induced transparency nanocavity.[J] Optics Express. 2013, 21(17): 20291-20302.
[6] DENG Zilan , DENG Junhong, ZHUANG Xin , et al. Diatomic metasurface for vectorial holography[J]. Nano Letters, 2018,18(5):2885.-2892.
[7] JIANG Yalan, WANG Jicheng, WANG Yueke, et al. A MIM surface plasmon T-splitter based on a stub structure[J]. Acta Photonica Sinica, 2014，43(9):0923002
[8] WANG Jicheng, JIANG Yalan, WANG Yueke, et al. Directional couplers based on MIM plasmonic waveguide stuctures[J]. Chinese Journal of Lasers, 2015, 42(2): 0217001.
[9] WANG Lin, ZHANG Lei. Narrow-spectrum enhanced sensor based on surface plasmon resonator[J]. Acta Optica Sinica, 2021, 41(7): 0724001.
[10] YANG Hongyan, CHEN Yupeng, XIAO Gongli, et al. MIM tunable plasmonic filter embedded with symmertrical sector metal resonator[J]. Acta Optica Sinica, 2020, 40(11): 1124001.
[11] PANG Shaofang, QU Shixian, ZHANG Yongyuan, et al. Filter characteristic research of MIM waveguide based on L shaped resonator[J]. Acta Optica Sinica, 2015, 35(6): 0623001.
[12] LIU Hairui, ZHANG Guanmao, Wang Zhishuang, et al. Design of the square concave ring resonantor MIM filter based on the surface plasmon polaritons[J]. Acta Photonica Sinica, 2018, 47(2); 0223002.
[13] WU Meng, LIANG Xiyin, YAN Changlin, et al. Design of arch-type resonance cavity filter based on surface plasmon polaritons[J]. Laser&Optoelectronics Progress, 2019, 56(20); 199-206.
[14] KEKATPURE R D, HRYCIW A C, BARNARD E S, et al. Solving dielectric and plasmonic waveguide dispersion relations on a pocket calculator[J]. Optics Express, 2009, 17(26): 24112-24129.
[15] YANG Yunru. Numerical study of plasmonic filter based on metal-insulator-metal waveguide[D]. Nanjing:Nanjing University of Posts and Telecommunications, 2016.
[16] KIM K Y, CHO Y K, TAE H S, et al. Light transmission along dispersive plasmonic gap and its subwavelength guidance characteristics[J]. Optics Express, 2006, 14(1): 320-330.
[17] SHI Songsong, WEI Zhongchao, Lu Zhiyang, et al. Enhanced plasmonic band-pass filter with symmetric dual side-coupled nanodisk resonators[J]. Journal of Applied Physics, 2015, 118(14): 143103.
热门排行
