删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

金属-介质-金属多层结构可调谐Fabry-Perot共振及高灵敏折射率传感

本站小编 Free考研考试/2021-12-29

摘要:顶层透光、底层不透光的金属-介质-金属多层结构可以产生窄带完美吸收共振, 用于测量介质层待测液体的折射率变化. 本文通过构建Fabry-Perot共振解析模型, 准确复现了该结构的响应光谱, 给出了其共振波长、品质因子、半高全宽和灵敏度的解析表达式, 并分析了介质层厚度对光谱的谐振波长和线宽的调谐机制, 明确了其物理机理. 基于8阶Fabry-Perot共振的金属-介质-金属多层结构, 用于折射率传感时的品质因子及优值分别达到2162.8和1648.1 RIU–1. 针对极小的折射率扰动, 通过在奇异点状态叠加Fabry-Perot共振的调谐机制, 提出了通过测量奇异点波长处反射系数的增加量或散射矩阵本征值的分裂量, 实现对待测液体折射率的可调谐式传感的方案. Fabry-Perot模型解析结果显示, 当待测液体的折射率变化为10–4 RIU时, 基于8阶Fabry-Perot共振的金属-介质-金属多层结构的前向反射系数增加量和本征值分裂量分别达到0.319和1.1279.
关键词: 金属-介质-金属谐振器/
Fabry-Perot共振模型/
可调谐式传感/
奇异点

English Abstract


--> --> -->
金属-介质-金属(metal-dielectric-metal, MDM)谐振器是等离激元共振传感的典型结构, 主要包括平面分层结构、精心设计的超构材料或超构表面[1]. 基于等离激元的超构材料或超构表面谐振器通常需要复杂且昂贵的微纳米加工工艺, 而基于Fabry-Perot共振的平面分层结构仅需要成本相对低廉的薄膜工艺, 是吸收器[1-3]、彩色滤光片[4-6]、传感器[7-10]等光学器件的首选. 由于金属材料具有显著的欧姆损耗, 等离激元共振谐振器的响应线宽通常为几十或几百纳米, 严重阻碍了其在窄线宽高灵敏传感检测方面的应用. 为降低欧姆损耗的影响, 可采用增益补偿或金属-介质混合纳米结构等方法[11-13]. 2019年, Qin等[9]报道了一种分布式布拉格反射器-隔离层-金层集成的平层结构折射率传感器, 其中的Tamm等离激元共振(tamm plasmon resonance, TPR)可使电磁场严格限制在待测液体腔内, 吸收谱线的半高全宽(full width at half maxima, FWHM)为2.2 nm, 传感灵敏度(sensitivity, S)为860 nm/RIU, 优值(figure of merit, FOM)达到391 RIU–1. 2020年, Liao等[13]设计了介质阵列和金属平面相结合的完美吸收结构, 结合了介质材料的低损耗和金属材料的高灵敏特性, 反射谱线的FWHM为0.38 nm, S和FOM分别是1052 nm/RIU和2768 RIU–1. 然而, 上述降低欧姆损耗的方法通常伴随着较为复杂的材料配置和工艺流程. 本文设计了MDM多层窄带完美吸收结构, 在不改变平面分层结构简单工艺的前提下, 提出了通过调节介质层厚度产生窄带高阶Fabry-Perot共振, 进而实现高灵敏折射率传感的方法.
通常折射率传感以扫描波长方式进行实验测量, 通过灵敏度S、品质因子(quality factor, Q)和优值FOM (${\rm{FOM}} = {S / {{\rm{FWHM}}}}$)进行性能评估, S是单位折射率变化的谐振波长偏移量, FWHM是响应光谱的半高全宽, Q是谐振波长和FWHM的比值[12]. 2007年, Chen等[14,15]提出通过增加导波层厚度, 可以在入射光入射角较小时激发MDM结构的高阶导模, 进而实现高Q共振, 促进传感灵敏度的提升. 之后, Liao和Zhao[16]以及Li等[17]提出入射光垂直入射到MDM结构时, 增加介质层厚度激发的高阶Fabry-Perot共振在高灵敏传感设计中具有潜在价值. 高阶Fabry-Perot窄带共振可以提供更大的Q因子和FOM值, 但实验中由于人眼或光谱仪的分辨率有限, 可以检测的最小折射率微扰量通常较大. 针对较小折射率扰动的传感, 可选择更易被测量的固定波长下的反射率或透射率的变化量作为观测量, 本文基于此传感测量方式, 借助散射矩阵奇异点理论对MDM多层结构的折射率传感进行了分析. “奇异点”(exceptional point, EP)是指非厄米系统中至少两个本征值和相应本征向量相等的简并点[18,19]. 与厄米系统中本征值相等、本征向量正交的“狄拉克点”(diabolic point, DP)不同, 奇异点处本征向量形成的希尔伯特空间不再完备, 导致非厄米简并系统对外部微扰异常敏感, 为弱微扰超灵敏传感提供了新的途径[20,21]. 由于存在欧姆损耗, MDM谐振器属于非厄米系统. 根据散射矩阵理论[7,8,22]可知, 当前向反射系数或后向反射系数中的一个为0时, 系统处于奇异点简并态. 本文研究表明, 将波长固定在奇异点波长处, 折射率发生微弱变化时, 前向反射系数变化量和本征值分裂量会随着Fabry-Perot共振阶次的增加而增加, 通过对介质层厚度优化取值, 可以设计符合不同灵敏度要求的传感器.
基于上述原理, 本文设计了一种顶部金属层透光、底部金属层不透光的MDM窄带吸收结构, 通过构建Fabry-Perot共振解析模型, 可以精确复现多层结构的响应光谱, 并在理论上阐明了介质层厚度对谐振波长和反射光谱半高全宽实现调谐的物理机理. 当该结构以谐振波长作为观测量进行折射率传感时, 其灵敏度S为789 nm/RIU, 8阶Fabry-Perot共振的Q因子和FOM值分别是2162.8和1648.1 RIU–1. 当该结构的探测波长固定在奇异点波长处, 以反射率作为观测量进行折射率传感时, 灵敏度会随着Fabry-Perot共振阶次的增加而增加. 此时, 若折射率发生10–4 RIU的微弱变化, 在8阶Fabry-Perot共振奇异点波长处, 前向反射系数增加量和本征值分裂量分别达到0.319和1.1279.
2
2.1.数值仿真
-->MDM多层结构的示意图如1(a)所示, x方向和y方向是无限延伸的金属-介质-金属三层结构, 上、下金属层均沉积在折射率为1.5的玻璃衬底上, 上、下金属层之间是待测液体空腔, 可以通过控制上下金属层的距离改变待测液体空腔的厚度, 入射平面波垂直于各层界面入射(沿z轴负方向传播). 金属层银采用Drude模型, 折射率${n_{{\rm{silver}}}} $$ = \sqrt {{\varepsilon _\infty } - {{\omega _{\rm{p}}^2}}/{{\left( {\omega + {\rm{i}}\gamma \omega } \right)}}}$, 其中${\varepsilon _\infty } = 1$, 等离子体频率${\omega _{\rm{p}}} = 1.366 \times {10^{16}} \;{{{\rm{rad}}} / {\rm{s}}}$, 阻尼常数$\gamma = 3.07\; \times $$ {10^{13}}\; {{\rm{s}}^{ - 1}}$[23], 上下银层厚度分别为d1 = 41.54 nm和d2 = 150 nm. 待测液体的初始折射率n设为水的折射率1.312, 厚度d = 350 nm.
图 1 (a) MDM多层结构示意图, 箭头表示均匀平面波的传播方向, 红色、蓝色箭头对应前向、后向入射均匀平面波; (b) 前向反射率R、透射率T和吸收率A的光谱; (c) 不同玻璃层厚度dGlass对应的反射率谱, 蓝色实线表示玻璃厚度无限大时的反射率谱(即图(b)中的红色实线), 插图为上下玻璃层均为有限厚度时的结构示意图; (d) 双层增透膜结构的反射率谱, 插图表示位于空气、玻璃之间的双层增透膜结构; (e) 不同玻璃层厚度dGlass对应的反射率谱, 蓝色实线表示玻璃厚度无限大时的反射率谱, 插图为上下玻璃层均为有限厚度且玻璃表面镀双层增透膜的结构示意图
Figure1. (a) Schematic of MDM multilayer structure. The solid arrows indicate the propagation directions of plane waves. The red and blue arrows correspond to forward and backward travelling incident plane waves, respectively. (b) Spectra of forward reflectance R, transmittance T and absorptance A of the proposed structure. (c) Reflectance spectra for different thicknesses dGlass of the glass layer, where the blue solid curve corresponds to infinite glass thickness [i.e. the red solid curve in Fig. (b)]. The inset shows the structure diagram with upper and lower glass layers set to have a finite thickness. (d) Reflectance spectra of the double-layer antireflection coating. The inset shows the structure of the double-layer antireflection coating between the air and glass regions. (e) Reflectance spectra for different thicknesses dGlass of the glass layer, where the blue solid curve corresponds to infinite glass thickness. The inset shows the structure diagram where the upper and lower glass layers are both set to have a finite thickness and both glass surfaces are coated with the double-layer antireflection coating.

采用严格耦合波分析(rigorous coupled wave analysis, RCWA)[24]方法进行数值仿真, 由于结构沿x, y方向无变化, RCWA中只需选取0级谐波. 图1(b)是该MDM结构的反射率、透射率及吸收率光谱, 共振波长是1035.967 nm, 此时的反射率$R = 1.32 \times 1{0^{ - 8}}$, 透射率$T = 5.31 \times 1{0^{ - 5}}$, 吸收率$A = 1 - R - T$, 表现出近乎完美的窄带吸收特性. 考虑到实际结构的玻璃层厚度有限, 将上下玻璃层设置成相同的有限厚度dGlass (图1(c)插图), 平面波从顶部的空气区域入射. dGlass取值不同时仿真的反射率谱如图1(c)所示, 可以看到有限玻璃厚度与无限玻璃厚度对应的反射率曲线有很大的差别, 即玻璃厚度对结构特性存在较大影响. 针对该问题, 实际应用中可以通过在玻璃表面镀宽带增透膜[25,26], 使空气中选定波长范围的入射平面波全部透射到玻璃层, 实现与玻璃衬底厚度无限时(此时平面波从玻璃中入射)相同的效果. 为了证明上述方法具有可行性, 图1(d)插图设计了针对空气-玻璃界面的双层增透膜, 紧挨玻璃的介质层折射率为n1 = 1.7, 紧挨空气的介质层折射率为n2 = 1.38, 两层介质的光学厚度(即折射率与几何厚度的乘积)均为谐振波长的1/4(谐振波长取值1035.967 nm). 根据菲涅耳公式可知, 光垂直入射到空气-玻璃界面的反射率为0.04, 加入双层增透膜后得到的反射率谱如图1(d)所示, 在所选波长范围内反射率被极大地降低. 将图1(c)对应结构的上下玻璃表面均镀上双层增透膜, 如图1(e)插图所示. 此时dGlass取值不同时仿真的反射率谱如图1(e)所示, 可见加入增透膜后, 玻璃厚度变化时的反射率曲线与玻璃厚度无限时的反射率曲线符合很好. 因此, 通过在有限厚度的玻璃上镀增透膜, 能够实现与玻璃衬底厚度无限大时相同的效果. 为了便于分析, 下文的数值仿真和模型计算均基于玻璃衬底厚度无限的假设.
为了计算此传感器的灵敏度, 将中间层空腔中待测液体的折射率从1.312逐渐增加到1.352, 递增间隔为0.01. 图2(a)是5种不同折射率液体的反射率谱曲线, 随着被测液体折射率的增大, 共振峰发生红移. 图2(b)是共振波长随折射率n的变化曲线, 通过线性拟合可得灵敏度S = 789 nm/RIU.
图 2 (a) 不同折射率n对应的反射率谱线; (b) 谐振波长随被测液体折射率n变化的曲线
Figure2. (a) Reflectance spectra for different refractive indices n; (b) resonance wavelength plotted as a function of the refractive index n of the measured liquid.

2
2.2.Fabry-Perot模型
-->此结构的窄带吸收特性来自Fabry-Perot共振, 上、下不对称银层是形成Fabry-Perot共振腔的两个反射面, 共振时入射光能量被限制在待测液体空腔内. 下面通过建立Fabry-Perot共振解析模型, 解释谐振的形成机制并分析谐振的可调特性. 图3(a1)(a4)是针对MDM结构构建的Fabry-Perot模型的参数定义, 对于本文的平面多层结构, a, b分别表示介质层中的下行和上行平面波模式系数, ${r_{{\rm{FP}}}}$, ${t_{{\rm{FP}}}}$分别表示平面波的反射和透射系数, 满足下列模式耦合方程:
图 3 Fabry-Perot模型 (a1) MDM结构中待求解的反射系数rFP, 透射系数tFP, 介质层中模式系数a, b的定义; (a2)—(a4) Fabry-Perot模型中双金属界面散射系数r1, t1, r2, t2, r3, t3的定义; (b1)—(b4) 对于r1, t1的Fabry-Perot模型, 中间金属层模式系数a', b'的定义, 以及单界面散射系数定义; (c1)—(c4) 对于r2, t2的Fabry-Perot模型, 相应的金属层模式系数和单界面散射系数定义; (d1)—(d4) 对于r3, t3的Fabry-Perot模型, 相应的金属层模式系数和单界面散射系数定义
Figure3. Fabry-Perot model: (a1) Definitions of reflection coefficient rFP, transmission coefficient tFP to be solved in MDM structure, and mode coefficients a, b in dielectric layer; (a2)–(a4) definitions of bimetal interface scattering coefficients r1, t1, r2, t2, r3, t3 in Fabry-Perot model; (b1)–(b4) Fabry-Perot model of r1 and t1, definitions of mode coefficients a', b' in intermediate metal layer, and single interface scattering coefficients; (c1)–(c4) Fabry-Perot model of r2 and t2, the corresponding metal layer mode coefficients and single interface scattering coefficients are defined; (d1)–(d4) Fabry-Perot model of r3 and t3, the corresponding metal layer mode coefficients and single interface scattering coefficients are defined.

$a = {t_1} + bu{r_3}, \tag{1a}$
$b = au{r_2},\tag{1b}$
${r_{{\rm{FP}}}} = {r_1} + bu{t_3},\tag{1c}$
${t_{{\rm{FP}}}} = au{t_2},\tag{1d}$
其中$u = \exp \left( {{\rm{i}}knd} \right)$是介质薄膜中平面波沿z方向的相移因子, 波数$k = {{2{\text{π}}} / \lambda }$, $\lambda $表示波长, ${t_1}$${r_1}$表示入射平面波激励介质层中平面波的透射和反射系数, ${t_2}$${r_2}$表示介质层中平面波在介质层底端的透射和反射系数, ${t_3}$${r_3}$表示介质层中平面波在介质层顶端的透射和反射系数. 求解方程(1a)和(1b)可确定a, b, 然后代入方程(1c), 可得
${r_{{\rm{FP}}}} = {r_1} + {t_3}\frac{{{r_2}{t_1}{u^2}}}{{1 - {r_2}{r_3}{u^2}}},\tag{2}$
于是得到平面波的反射率${R_{{\rm{FP}}}} = {\left| {{r_{{\rm{FP}}}}} \right|^2}$.
平面波在不同介质分界面的反射和透射系数可以根据菲涅耳反射和透射公式计算. 当入射光从介质1垂直入射到介质2时, 单一界面上平面波的反射和透射系数分别为$r = \dfrac{{{n_1} - {n_2}}}{{{n_1} + {n_2}}}$, $t = \dfrac{{2{n_1}}}{{{n_1} + {n_2}}}$, n1n2分别表示介质1和介质2的折射率. 图3(a2)(a4)中定义的反射和透射系数涉及2个金属界面, 本文针对2个金属界面结构进一步建立Fabry-Perot模型, 具体的参数定义如图3(b1)(b4)图3(c1)(c4)图3(d1)(d4)所示, 可以借助菲涅耳公式给出${r_1}$, ${t_1}$, ${r_2}$, ${t_2}$, ${r_3}$, ${t_3}$的解析表达式. 这里以${r_1}$, ${t_1}$为例, 借助和(1)式类似的模式耦合方程可以得到:
${r_1} = {r_1'} + {t_3'} \frac{{{r_2'} {t_1'} {{u'}^2}}}{{1 - {r_2'} {r_3'} {{u'}^2}}},\tag{3a}$
${t_1} = \frac{{{t_1'} {t_2'} u'}}{{1 - {r_2'} {r_3'} {{u'}^2}}},\tag{3b}$
其中$u' = \exp \left( {{\rm{i}}k{n_{{\rm{silver}}}}{d_1}} \right)$是上层金属中平面波沿z方向的相移因子, ${r_1'}$, ${t_1'}$, ${r_2'}$, ${t_2'}$, ${r_3'}$, ${t_3'}$是单一界面的反射和透射系数, 解析表达式由菲涅耳公式给出(如前文所述). 同理可以得到${r_2}$, ${t_2}$, ${r_3}$, ${t_3}$的解析表达式. 至此, 平面波入射MDM结构的反射和透射系数${r_{{\rm{FP}}}}$, ${t_{{\rm{FP}}}}$被完全解析化, 这里不再给出具体的表达式. 值得注意的是, ${r_1}$, ${t_1}$${r_3}$, ${t_3}$仅与分界面折射率及d1有关, ${r_2}$, ${t_2}$仅与分界面折射率及d2有关.
基于上述Fabry-Perot模型, 改变介质层厚度能够调节谐振波长和光谱线宽, 下面具体分析. 首先分析谐振波长${\lambda _0}$的调节机制. 发生Fabry-Perot共振时, 方程(2)的分母取极小值, 由此得到, $\arg \left( {{r_2}} \right) + $$ \arg \left( {{r_3}} \right) + {{2 nd} / {{\lambda _0}}} \times 2{\text{π}} = 2 m{\text{π}}$$( m = 0, 1, 2, \cdots )$, 其中m为共振阶次. 从而可确定谐振波长:
${\lambda _0} = \frac{{4{\rm{\pi }}nd}}{{2m{\rm{\pi }} - \arg \left( {{r_2}} \right) - \arg \left( {{r_3}} \right)}}.$
考虑到银层折射率与近红外波长的缓变关系, 仅改变d时, 可近似认为${r_2}$, ${r_3}$是常数, 根据(4)式可知${\lambda _0}$d呈线性关系. 图4(a)是不同介质层厚度的反射率谱曲线, d增大时${\lambda _0}$红移. 谐振波长与介质层厚度的关系曲线如图4(b)所示, 可见由反射率谱极小值点得到的谐振波长与方程(4)计算的结果一致.
图 4 (a) 不同介质层厚度d对应的反射率谱曲线, 点表示RCWA数值仿真结果, 实线表示Fabry-Perot模型计算结果; (b) 谐振波长随介质层厚度d的变化, 红色和蓝色实线表示Fabry-Perot模型和方程(4)的计算结果; (c) 不同谐振阶次m对应同一谐振波长1035.967 nm时(m越大则d越大), RCWA计算(点)、Fabry-Perot模型(实线)给出的反射率谱曲线; (d) Q, FOM, Sd的变化, 红色、蓝色和灰色实线代表由Fabry-Perot模型计算的反射率谱得到的Q, S, FOM, 粉色和青色虚线代表方程(5)和(6)计算的Q, S
Figure4. (a) Reflectance as a function of wavelength for variable dielectric layer thicknesses d, dots represent RCWA numerical simulation results, solid lines represent Fabry-Perot analytical model calculation results; (b) resonance wavelength as a function of dielectric layer thickness d, red and blue solid lines represent the calculation results of Fabry-Perot model and Eq. (4); (c) reflectance spectrums given by RCWA calculation (dots) and Fabry-Perot model (solid lines), and different resonance orders m correspond to the same resonance wavelength 1035.967 nm (larger m is, larger d is); (d) Q, FOM, S as a function of d, red, blue and gray solid lines represent Q, S, FOM obtained from reflectance spectrum calculated by Fabry-Perot model, pink and cyan dotted lines represent Q, S calculated by Eqs. (5) and (6).

下面考虑光谱线宽FWHM的调谐机制. 由于MDM多层结构是与外界有能量交换的开放光学谐振腔, 其本征模式(即无源麦克斯韦方程组的解)为准简正模式(quasi-normal mode, QNM)[27], 具有有限光子寿命和复数谐振频率${\omega _m} = {\omega _m'} \, + $$ {\rm{i}}{\omega _m ''} $, 相应的复数谐振波长为${\lambda _m} = 2{\text{π}}c/{\omega _m} = {\lambda _m '} + $$ {\rm{i}}{\lambda _m''}$ (c为真空中光速), ${\lambda _m'} $接近(4)式中的${\lambda _0}$. 根据QNM展开定理, 当平面波激励频率趋于复数谐振频率, 即$\omega = {\omega _m}$时, 平面波激发的电磁场趋于无穷大, 于是, 根据方程(2), 可得${\omega _m}$满足[28]$1 - $$ {r_2}{r_3}{u^2} = 0$, 由此可得, ${{{\omega _m}}}/{c} = [{\rm{i}}\ln |{r_2}{r_3}| - \arg ({r_2}) - $$ \arg ({r_3}) + 2 m{\text{π}}]/{{2 nd}}$. 进而, 得到QNM的品质因子[27] $Q = \dfrac{{{\omega _m'} }}{{ - {\rm{2}}{\omega _m''}}} = \dfrac{{2 m{\text{π}} - \arg ({r_2}) - \arg ({r_3})}}{{ - {\rm{2}}\ln |{r_2}{r_3}|}} \approx \dfrac{{2{\text{π}}nd}}{{ - {\lambda _0}\ln |{r_2}{r_3}|}},$ 其中利用了方程(4). 故而, Fabry-Perot共振谱线的线宽可近似表示为
${\rm{FWHM}} = \frac{{{\lambda _m'}}}{Q} = - \frac{{\lambda _0^2\ln \left| {{r_2}{r_3}} \right|}}{{2{\rm{\pi }}nd}}.$
由(5)式可知, FWHM反比于d, 因此通过改变介质层厚度d可以调节响应光谱的线宽. 为了验证上述线宽调谐机制, 图4(c)给出了m取值不同时, RCWA数值仿真和Fabry-Perot模型计算得到的反射率谱线. 为了便于比较, 各个谱线的谐振波长统一调至1035.967 nm. 此时, 方程(4)中m = 0, 2, 4, 6, 8对应的介质层厚度分别为350, 1139.609, 1929.218, 2718.827, 3508.436 nm. 由图4(c)可知, Fabry-Perot模型能精确复现RCWA数值仿真的结果, 所有谱线在共振波长处的反射率都接近0, 且m越大(则介质层厚度d越大), FWHM越小, 与方程(5)一致. 图4(d)给出了反射率谱线的Qd变化的曲线, RCWA数值仿真中(红色实线), Q值定义为谐振波长和反射率光谱半高全宽的比值, Fabry-Perot解析模型中(粉色虚线), Q值定义为谐振波长和(5)式计算得到的半高全宽的比值. 结果表明, 二者一致, 并且都正比于d. 当d = 3508.436 nm(m = 8)时, 对应的FWHM和Q值分别是0.479 nm和2162.8, 约是d = 350 nm(m = 0)时的0.111倍和9倍.
下面根据上述解析模型, 分析MDM结构传感器的S和FOM随介质层厚度的变化. 考虑到反射系数${r_2}$, ${r_3}$与待测液体折射率n的缓变关系, 可近似认为二者是常数. 将(4)式对折射率求偏导得到灵敏度的表达式为
$S = \frac{{\partial {\lambda _0}}}{{\partial n}} = \frac{{4{\rm{\pi }}d}}{{2m{\rm{\pi }} - \arg \left( {{r_2}} \right) - \arg \left( {{r_3}} \right)}} = \frac{{{\lambda _0}}}{n}, $
其中最后一个等号利用了方程(4). 对于不同阶次Fabry-Perot共振, 其用于折射率传感时S的变化曲线如图4(d)中蓝色实线和青色虚线所示. 可见随着Fabry-Perot共振的阶次m增大(同时介质层厚度d增大), 由模型计算的反射率谱得到的灵敏度与(6)式计算结果一致, S均近似保持不变(取值为${{{\lambda _0}} / n}$). 考虑到共振阶次增加时, FWHM线性下降而S近似保持不变, 所以FOM呈线性增加(如图4(d)中灰色实线所示). 其中, 8阶Fabry-Perot共振的FOM值是1648.1 RIU–1, 是目前折射率传感器中的较大数值[29-34].
对于微小的折射率变化, 谐振波长的极小偏移量会受到光谱仪分辨率的限制, 导致传感检测失效, 此时可选择更易被观测的固定波长处反射率变化量作为观测量, 这里借助散射矩阵的奇异点理论对此类折射率传感进行分析. 可采用散射矩阵构造MDM平面分层结构的非厄米等效哈密顿量, 通过求解散射矩阵${{M}}\! =\! \left( {\begin{array}{*{20}{c}} t&{{r_{\rm{b}}}} \\ {{r_{\rm{f}}}}&t \end{array}} \right)$的本征方程${{Mx}} = $$ v{{x}}$(设散射矩阵的本征值为$v$, 对应于本征值$v$的本征向量为${{x}}$), 可得到本征值为${v_{1, 2}} = t \pm \sqrt {{r_{\rm{b}}}{r_{\rm{f}}}} $, 本征向量为${{{x}}_{1, 2}} = \left( {\sqrt {{r_{\rm{b}}}}, \pm \sqrt {{r_{\rm{f}}}} } \right)$. ${r_{\rm{f}}}$, ${r_{\rm{b}}}$, t分别为前向、后向反射系数、透射系数, 如图1(a)所示. 值得注意的是, 由于互易定理[35], 前向和后向透射系数相等, 因此都用t表示. 当${r_{\rm{b}}}$, ${r_{\rm{f}}}$中的一个为0时(即单向零反射), 散射矩阵成为类Jordon块矩阵形式, 此时发生EP简并[7,20,36,37]. 由于完美吸收为单向零反射的特例(即${r_{\rm{f}}} = t = 0$), 因此图1中MDM多层结构在谐振波长处形成了EP. 为了显示该MDM结构的EP简并特性, 图5给出了该结构的M矩阵的元素及其本征值随波长的变化曲线. 图5(a)表明, 在谐振波长处(即EP处), ${r_{\rm{f}}}$几乎为0, $|{r_{\rm{b}}}|$则接近1. 图5(b)表明, t接近0. 图5(c)图5(d)显示, 两个本征值的实部和虚部在谐振波长处同时相等, 在谐振波长附近又分别表现出交叉和排斥的现象[8], 这是奇异点简并的典型特征.
图 5 奇异点简并 (a) 前向(红色实线)和后向(蓝色虚线)反射系数光谱; (b) 前向(红色实线)和后向(蓝色虚线)透射系数光谱; (c), (d) 散射矩阵本征值ν1 (红色实线)、ν2 (蓝色虚线)的实部和虚部随波长变化的曲线
Figure5. Exceptional point degeneracy: (a) Spectra of forward (red solid line) and backward (blue dashed line) reflection; (b) spectra of forward (red solid line) and backward (blue dashed line) transmission; (c), (d) real and imaginary parts of eigenvalue ν1 (red solid line), ν2 (blue dashed line) of the scattering matrix vary with wavelength.

基于上述分析, 下面将观测波长固定为奇异点波长1035.967 nm, 以反射率作为观测量, 实现EP处折射率传感. 在EP波长处, 前向反射系数为0, 且散射矩阵的本征值发生简并, 加入折射率微扰后, 系统远离EP简并点, 前向反射系数及本征值分裂量发生显著增大, 增大数值远高于DP简并的情形, 从而有利于提高传感灵敏度[7,8]. 此时传感灵敏度定义为微弱折射率变化导致的前向反射系数增加量, 或者本征值分裂量. 为了直观显示上述传感机制, 这里以图6(a)d = 3508.436 nm (m = 8)的奇异点结构为例进行说明, 图6(b)给出了其在不同折射率微扰下的前向反射系数光谱曲线. 可见随着折射率的增大, 共振波长红移, 观测波长处的前向反射系数逐渐增大. 当折射率从1.312增加至1.3121时, 谐振波长偏移量为0.079 nm, 前向反射系数增加量$\Delta \left| {{r_{\rm{f}}}} \right| = 0.319$, 本征值分裂量$\left| {\Delta v} \right| = 1.1279$.
图 6 (a) 基于8阶Fabry-Perot共振的EP简并系统; (b) EP简并系统的前向反射系数光谱随折射率n的变化, 竖直黑色实线表示观测波长; (c) DP简并系统; (d) 对称银层系统; (e) 前向反射系数随折射率微扰的变化曲线, 圆圈实线对应基于Fabry-Perot共振(取不同阶次m)的EP简并系统, 三角和方形实线分别对应DP简并系统和对称银层系统; (f) m = 0和m = 8时本征值分裂量随折射率微扰的变化曲线
Figure6. (a) EP degenerate system based on 8 th order Fabry-Perot resonance; (b) forward reflection coefficient spectra of the EP degenerate system versus refractive index n, where the vertical black solid line represents the observation wavelength; (c) DP degenerate system; (d) symmetrical silver layer system; (e) forward reflection coefficient as a function of refractive index perturbation, where circle solid lines correspond to the EP systems based on Fabry-Perot resonance (taking different orders m), and triangular and square solid lines correspond to the DP degenerate system and the symmetrical silver layer system, respectively; (f) eigenvalue splitting amount as a function of the refractive index perturbation for m = 0 and m = 8.

为了验证EP简并系统具有更高的传感灵敏度, 可以设计DP简并系统和对称银层系统进行对比. DP简并系统(图6(c))由厚度为3553.241 nm的待测液体层组成, 观测波长处其前向反射系数和后向反射系数均为0. 对称银层系统 (图6(d))由 厚度为41.54 nm的对称银层和厚度为3507.505 nm的待测液体层组成, 观测波长处其前向反射系数和后向反射系数相等但不等于0, 既不处于DP简并状态又不处于EP简并状态. 图6(e)中紫色三角、紫色方形以及紫色圆圈实线分别是观测波长处的DP简并系统、对称银层系统以及介质层厚度相差不多的EP简并系统的前向反射系数随折射率微扰的变化曲线. 当折射率从1.312增加至1.3121时, DP简并系统、对称银层系统、以及EP简并系统的前向反射系数增加量分别为2.893 × 10–4, 0.0517, 0.319. 比较可知, 针对弱微扰传感时EP简并系统具有更高的灵敏度.
考虑到Fabry-Perot共振的可调谐性, Fabry-Perot共振的阶次m增加时(同时介质层厚度d增大, 以维持谐振波长λ0不变), 由于灵敏度S近似保持不变(基于方程(6)), 因此相同的折射率微扰导致的谐振波长偏移量相同, 而反射率谱曲线的线宽减小(基于方程(5)), 这导致在固定的观测波长处, 高阶Fabry-Perot共振结构的反射系数增加量较大. 故而, 可以采用上述EP结构, 通过增大共振阶次m以获得更高的灵敏度. 基于不同共振阶次m的EP结构前向反射系数随折射率变化的曲线如图6(e)所示, 随着中间层厚度d的增大, Fabry-Perot共振的阶次m逐渐增大, 相同折射率变化导致的前向反射系数增加量越大. 图6(f)给出了m = 0和m = 8时观测波长处本征值分裂量$\left| {\Delta v} \right| = \left| {{v_1} - {v_2}} \right|$随折射率微扰的变化曲线. 折射率增加10-4 RIU时, m = 8对应的前向反射系数增加量($\Delta \left| {{r_{\rm{f}}}} \right|$ = 0.319)和本征值分裂量($\left| {\Delta v} \right|$ = 1.1279)分别是m = 0对应的前向反射系数增加量($\Delta \left| {{r_{\rm{f}}}} \right|$ = 0.0366)和本征值分裂量($\left| {\Delta v} \right|$ = 0.3826)的8.72倍和2.95倍. 由此证实了通过提高Fabry-Perot共振阶次m, 能够提高EP处折射率传感的灵敏度.
本文设计并数值仿真了MDM多层奇异点结构形成的窄带完美吸收器, 该结构可用于超灵敏传感器的可调谐式传感. 通过构建Fabry-Perot解析模型准确复现了其响应光谱, 给出了共振波长、Q因子、FWHM和S的解析表达式, 揭示了共振波长与介质层厚度成正比, FWHM与介质层厚度成反比的谱线调谐机制. 当该结构的8阶Fabry-Perot共振用作折射率传感时, Q因子和FOM值分别达到2162.8和1648.1 RIU–1.
基于MDM多层结构在谐振波长处于奇异点简并状态, 当折射率变化较小时, 可通过固定波长测量反射率的方式进行传感测量. 通过测量奇异点波长处反射系数的增加量或者散射矩阵本征值的分裂量, 再叠加不同阶次Fabry-Perot共振的调谐机制, 可实现高灵敏度可调谐式折射率传感. 当该结构的8阶Fabry-Perot共振用作折射率传感时, 被测液体的折射率发生10–4 RIU的微弱变化导致的前向反射系数增加量是0.319, 本征值分裂量是1.1279. 若继续增大谐振阶次m, 该MDM结构的前向反射系数增加量和本征值分裂量将取更大的数值.
相关话题/结构 金属 系统 观测 光谱

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于深度学习的相位截断傅里叶变换非对称加密系统攻击方法
    摘要:大多数光学加密系统都是对称加密系统,在光学图像加密中明文和密文之间具有线性关系,其系统的安全性有待加强.而基于相位截断傅里叶变换(phase-truncatedFouriertransform,PTFT)的非对称加密系统,其非线性的相位截断操作使加密系统的安全性得到了极大提升.本文提出使用深度 ...
    本站小编 Free考研考试 2021-12-29
  • 疏水表面振动液滴模态演化与流场结构的数值模拟
    摘要:为了研究疏水表面垂直振动液滴的运动规律,本文建立了振动液滴的三维模型,考虑了振动液滴的动态接触角变化过程,通过流体体积函数和连续表面张力(volumeoffluid-continuesurfaceforce,VOF-CSF)方法实现了液滴受迫振动的数值模拟,得到了液滴的四种模态(2,4,6和8 ...
    本站小编 Free考研考试 2021-12-29
  • 单分子表面增强拉曼散射的光谱特性及分析方法
    摘要:单分子检测代表分子检测灵敏度的极限,能够提供传统检测方式无法提供的物理信息,在化学分析、分子动力学机理、蛋白质解析等领域具有广阔的应用前景,具有重要的科学研究价值.具体而言,单分子检测能对复杂体系中的分子进行识别和计数,给出分子的分布信息;也可以对单个分子在吸附、反应等过程中的实时衍变进行追踪 ...
    本站小编 Free考研考试 2021-12-29
  • 光源光谱特性对空间相机调制传递函数检测的影响
    摘要:调制传递函数(modulationtransferfunction,MTF)检测是评价空间相机像质的重要手段.空间相机光学系统透过率和色差、探测器量子效率均与波长相关,采用不同光谱特征的光源所得到的MTF会出现偏差,光源光谱特性的影响不可忽略.针对这一问题,提出了一种分析光源光谱特性对空间相机 ...
    本站小编 Free考研考试 2021-12-29
  • 基于单元辐射叠加法的结构声源声场重建方法
    摘要:为了提高分布式结构声源的声场重建精度,本文提出了基于单元辐射叠加法的结构声源声场重建方法.该方法首先利用声场叠加原理和结构振声传递特性,建立了结构声源表面振动与辐射声场之间的振声传递关系解析表达式,得到便于快速计算的振声传递矩阵,能够解决连续分布、相干结构噪声源的声传播模型精细化表征问题.然后 ...
    本站小编 Free考研考试 2021-12-29
  • 液态金属软表面池沸腾传热的实验研究
    摘要:池沸腾是一种高效的传热方式,目前主要通过刚性固体表面改性强化沸腾传热.本文以乙醇为工质,实验研究了光滑铜表面和液态金属软表面池沸腾传热.发现液态金属软表面可有效降低沸腾起始点(ONB)壁面过热度,饱和沸腾时,ONB壁面过热度从光滑铜表面的约18℃降低到软表面的约6℃,沸腾传热系数最大提高了14 ...
    本站小编 Free考研考试 2021-12-29
  • In<sub>1+</sub><i><sub>x</sub></i>Te化合物的结构及热电性能研究
    摘要:InTe化合物中In+的孤对电子引发的晶格非谐性振动使得其具有本征极低的热导率,因此被认为是一种具有潜力的中温热电材料.然而,较低的电输运性能使得InTe的热电性能不高.在本工作中,采用熔融、退火结合放电等离子活化烧结工艺制备了一系列In1+xTe(x=0,0.001,0.003,0.005, ...
    本站小编 Free考研考试 2021-12-29
  • 第一性原理研究Mg掺杂对LiCoO<sub>2</sub>正极材料结构稳定性及其电子结构的影响
    摘要:LiCoO2作为商业化最早的锂离子电池正极材料,至今仍受到许多研究人员的广泛关注.高电压下LiCoO2面临着严重的容量衰减和性能下降等问题,实验上通常采用体相元素掺杂以稳定LiCoO2在高电压下的晶体结构,从而提高其电化学性能.Mg元素掺杂被认为是一种能够提高LiCoO2高电压循环稳定性的有效 ...
    本站小编 Free考研考试 2021-12-29
  • 溶剂热制备铬掺杂硫化锌和硫化纳米结构和磁性能
    摘要:以乙醇胺和乙二胺为混合溶剂,通过简单溶剂热法,用S,ZnO和CdO为源,成功制备了不同量Cr掺杂ZnS和CdS半导体纳米结构.X-射线衍射测试表明,纳米结构ZnS和CdS具有纤锌矿结构.扫描电子显微镜给出了不同铬含量的ZnS和CdS的形貌.用电子能量散射谱观察到产物的成分为Cr,Zn,Cd和S ...
    本站小编 Free考研考试 2021-12-29
  • 1907号台风“韦帕”登陆期间眼壁区的垂直电场探空观测
    摘要:2019年夏季在海南文昌获取到1907号台风“韦帕”眼壁内的一次电场探空资料,本文利用卫星、雷达、地面电场及海南省地闪定位资料对该台风特征进行详细分析,在此基础上,对探空路径区域内的电场廓线以及电荷区分布进行了分析,结果显示:在海拔高度5.74—9.10km之间共有4个正电荷区和3个负电荷区, ...
    本站小编 Free考研考试 2021-12-29