搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

含双挡板金属-电介质-金属波导耦合方形腔的独立调谐双重Fano共振特性

陈颖 曹景刚 谢进朝 高新贝 许扬眉 李少华

Wang Xiao-Xiao, Zhang Mei-Ling, Han Yue-Lin, Jia Li-Qun. Mei symmetry and Mei conserved quantity of Nielsen equation in a dynamical system of the relative motion with nonholonomic constraint of Chetaev's type. Acta Phys. Sin., 2012, 61(20): 200203. doi: 10.7498/aps.61.200203
Citation: Wang Xiao-Xiao, Zhang Mei-Ling, Han Yue-Lin, Jia Li-Qun. Mei symmetry and Mei conserved quantity of Nielsen equation in a dynamical system of the relative motion with nonholonomic constraint of Chetaev's type. Acta Phys. Sin., 2012, 61(20): 200203. doi: 10.7498/aps.61.200203

含双挡板金属-电介质-金属波导耦合方形腔的独立调谐双重Fano共振特性

陈颖, 曹景刚, 谢进朝, 高新贝, 许扬眉, 李少华

Resonance characteristics of independently tuned dual Fano of metal-dielectric-metal waveguide coupling square cavity with double baffles

Chen Ying, Cao Jing-Gang, Xie Jin-Chao, Gao Xin-Bei, Xu Yang-Mei, Li Shao-Hua
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 基于表面等离子激元在亚波长结构的传输特性, 设计了一种含双挡板金属-电介质-金属波导耦合两个方形腔的结构. 由F-P谐振腔产生的宽谱模式与两个方形谐振腔产生的两个窄谱模式发生干涉作用, 形成了独立调谐的双重Fano共振, 而且可以通过改变两个方形腔的大小及填充介质实现双重Fano共振的独立调谐. 基于耦合模理论, 定性分析了该结构产生双重Fano共振的机理. 利用有限元仿真的方法, 定量分析了结构参数对可独立调谐双重Fano共振和折射率传感特性的影响. 结果表明, 优化参数后该结构的灵敏度分别高达1020和1120 nm/RIU, FOM值分别高达3.59 × 105和1.17 × 106. 该结构可为超快光开关、多功能高灵敏度传感器和慢光器件的光学集成提供有效的理论参考.
    A metal-dielectric-metal (MDM) waveguide coupling two square cavities with double baffles is designed in this paper based on the transmission characteristics of surface plasmon polaritons in subwavelength structure. The independent tuning of the dual Fano resonance is implemented by the interference between the wide-spectrum mode generated by the F-P (Fabry Perot) cavity and the two narrow-spectrum modes generated by the two square cavities. Moreover, the independent tuning of the dual Fano resonance can be achieved by changing the sizes of the two square cavities and filling medium. The coupled-mode theory (CMT) is adopted to analyze the transmission characteristics of the dual Fano resonance. The structure is simulated by the finite element method to quantitatively analyze the influence of structural parameters on the independent tuning of the dual Fano resonance and the refractive index sensing characteristics. The proposed sensor yields respectively sensitivity higher than 1020 nm/RIU and 1120 nm/RIU and a figure of merit of 3.29 × 105 and 1.17 × 106 by optimizing the geometry parameters. This structure provides an effective theoretical reference in the optical integration of ultra-fast optical switches, multi-function high-sensitivity sensors and slow-light devices.
      PACS:
      81.07.-b(Nanoscale materials and structures: fabrication and characterization)
      81.15.Cd(Deposition by sputtering)
      68.55.-a(Thin film structure and morphology)
      74.25.nd(Raman and optical spectroscopy)
      通信作者: 陈颖, chenying@ysu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 61201112, 61475133)、河北省自然科学基金(批准号: F2016203188, F2016203245)、中国博士后基金(批准号: 2018M630279)、河北省高等学校科学技术研究项目(批准号: ZD2018243)资助的课题.
      Corresponding author: Chen Ying, chenying@ysu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 61201112, 61475133), the Natural Science Foundation of Hebei Province, China (Grant Nos. F2016203188, F2016203245), the China Postdoctoral Science Foundation (Grant No. 2018M630279), and the Scientific Research Foundation of the Higher Education Institutions of Hebei Province, China (Grant No. ZD2018243).

    由于在空间光通信、量子纠缠、粒子操控等领域的应用, 光学涡旋已成为近年来的研究热点. 与普通光束相比, 涡旋光束可以携带轨道角动量. 拉盖尔-高斯(L-G)光束就是一种典型的涡旋光束, 其具有螺旋相位结构和环形场强分布, 光束中心的强度为零[1]. L-G光束沿其传播方向携带每光子的轨道角动量, 表示L-G光场的拓扑荷数[2].

    实验表明, 通过交换轨道角动量, L-G光束可以对物体施加扭矩, 包括微观吸附颗粒[3]、介观玻色-爱因斯坦凝聚[4]及亚微米布朗粒子[5]. 2007年, 基于此性质, 拉盖尔-高斯旋转腔系统被提出[6], 该系统由两个用作腔镜的螺旋相位元件构成, 分别为一个固定腔镜和一个围绕腔轴旋转的旋转腔镜(可视为扭摆). 在这种光旋转系统中, 腔模是L-G光, 它携带轨道角动量, 可以与螺旋相位元件交换轨道角动量. 近期人们对此类系统的兴趣逐渐增加, 并做了一些相关研究, 包括旋转腔镜的冷却[6,7]、光场轨道角动量的检测[8,9]、二阶边带效应[10]和纠缠[11,12]. 然而, 据我们所知, 关于L-G光束在拉盖尔-高斯旋转腔中的非互易传输现象至今还未被报道过. 由于光旋转相互作用[13]和光机械相互作用[14,15]之间存在类似的哈密顿量, 光机械相互作用产生的许多效应[16-29], 在拉盖尔-高斯旋转腔中通过与L-G光场轨道角动量的交换也会发生[6-12]. 这为利用光的轨道角动量而不是动量提供了可能性.

    光的传播通常是互易的, 不过实现光学非互易对于构建非对称量子网络极为重要. 光学非互易性描述的是光场从一个方向通过光学系统但是不能沿原路返回的特性. 光学非互易是实现定向放大器、隔离器、循环器的基础[30], 并且需要打破时间反演对称性. 传统上, 非互易传输依赖于施加的偏置磁场 [31], 但这些传统器件通常体积庞大, 与超低损耗超导电路不兼容, 并且需要相当大的磁场. 近年来, 有人利用辐射-压力诱导光力耦合来打破时间反转对称性, 实现光的非互易效应[32-46]. 这一领域已经取得了重要成果, 包括非互易传输和放大[32-34]、非互易单光子效应[35]、非互易慢光[37]等. 文献[38]提出了一种基于具有光增益的光机械系统来实现光场和微波场之间光定向放大的方案, 发现放大方向可以通过调节有效光机耦合之间的相对相位(非互易相位差[46])来控制. 文献[43]给出了由两个光学膜和一个机械模组成的三模光机械系统中的光学非互易响应, 并证明了光学非互易响应是通过调节光机械耦合速率之间的相位差来打破系统的时间反转对称性来实现的. 然而, 通过光旋转耦合来实现携带轨道角动量的涡旋光束的非互易传输还有待进一步研究, 且之前工作在利用光力耦合来实现普通光束的非互易传输时, 都是利用光学模式与机械模式交换动量来实现, 而非交换角动量.

    本文通过在单个拉盖尔-高斯旋转腔中额外加一个固定腔镜来构建双拉盖尔-高斯旋转腔, 并通过光纤将两个腔模线性耦合, 并且用两个强驱动场驱动不同的腔模, 进而研究涡旋光束的非互易传输现象. 将以往通过光与机械振子交换动量来实现普通光束的非互易传输推广到光与旋转腔镜交换角动量来实现涡旋光束的非互易传输. 研究发现, 非互易相位差可以决定涡旋光束非互易传输的产生及方向; 在一定的拓扑荷比值下, 该系统可以实现完美的非互易性; 系统耗散对传输振幅的变化有比较大的影响; 此外, 通过调节非互易相位差, 在该系统中可以实现非互易的慢光效应. 在此系统中的涡旋光束传输的非互易性, 来源于光旋转耦合相互作用和腔模耦合相互作用形成的两条路径之间发生的量子干涉效应.

    图1所示, 本方案考虑了双拉盖尔-高斯旋转腔光旋转系统, 其中3个腔镜均为螺旋相位元件, 其中两个输入腔镜部分透明且刚性固定, 不改变透过它们光束的轨道角动量而会给反射光束增加2122的拓扑荷; 处于中间位置的旋转腔镜是完全反射的, 可以给反射光增加+21+22的拓扑荷, 并可以在支架上绕腔轴z旋转(角平衡位置ϕ0=0); 同时两个腔场通过光纤线性耦合在一起. 两个频率均为ωc, 振幅分别为εcεd的强驱动场分别从两侧注入系统驱动频率为ω0的L-G腔模. 利用一束频率为ωp、振幅为εL(εR)的弱探测场从系统左侧(右侧)入射, 从而检查该光旋转系统相对于左侧(右侧)探测光的响应特性. 在驱动场频率ωc的旋转框架下, 系统的哈密顿量( = 1)可写为

    图 1 双拉盖尔-高斯旋转腔示意图, 两L-G腔场通过光旋转相互作用与中间的旋转腔镜耦合, 振幅为$ {\varepsilon _{\text{c}}} $和$ {\varepsilon _{\text{d}}} $($ {\varepsilon _{\text{L}}} $和$ {\varepsilon _{\text{R}}} $)的强驱动场(弱探测场)分别从两侧入射系统来驱动L-G腔模$ {c_1} $和$ {c_{\text{2}}} $, 同时用光纤将两L-G腔模线性耦合. 旋转腔镜的平衡位置为$ {\phi _{\text{0}}} $, 在扭力作用下的角位移用$ \phi $角表示\r\nFig. 1. Schematic diagram of double Laguerre-Gaussian (L-G) rotational-cavity. The two L-G cavity modes are coupled with a rotating cavity mirror in the middle via the optical rotation interaction. Two strong pump fields (weak probe fields) with amplitudes $ {\varepsilon _{\text{c}}} $ and $ {\varepsilon _{\text{d}}} $ ($ {\varepsilon _{\text{L}}} $ and $ {\varepsilon _{\text{R}}} $) are incident on the system from both sides to drive the L-G cavity modes $ {c_1} $ and $ {c_{\text{2}}} $, and the two L-G cavity modes are linearly coupled with an optical fiber. The equilibrium position of the rotational mirror is $ {\phi _{\text{0}}} $, and the angular displacement is indicated by angle $ \phi $ under the action of the torsion.
    图 1  双拉盖尔-高斯旋转腔示意图, 两L-G腔场通过光旋转相互作用与中间的旋转腔镜耦合, 振幅为εcεd(εLεR)的强驱动场(弱探测场)分别从两侧入射系统来驱动L-G腔模c1c2, 同时用光纤将两L-G腔模线性耦合. 旋转腔镜的平衡位置为ϕ0, 在扭力作用下的角位移用ϕ角表示
    Fig. 1.  Schematic diagram of double Laguerre-Gaussian (L-G) rotational-cavity. The two L-G cavity modes are coupled with a rotating cavity mirror in the middle via the optical rotation interaction. Two strong pump fields (weak probe fields) with amplitudes εc and εd (εL and εR) are incident on the system from both sides to drive the L-G cavity modes c1 and c2, and the two L-G cavity modes are linearly coupled with an optical fiber. The equilibrium position of the rotational mirror is ϕ0, and the angular displacement is indicated by angle ϕ under the action of the torsion.
    H=Δc(c1c1+c2c2)+L2z2I+12Iω2ϕϕ2  +gϕ1c1c1ϕ+gϕ2c2c2ϕ+J(c1c2+c2c1)  +i(εcc1εcc1)+i(εdc2εdc2)  +iεL(c1eiδtc1eiδt)+iεR(c2eiδtc2eiδt)
    (1)

    式中, 第1项表示L-G腔模的自由哈密顿量, 其中c1(c1)c2(c2)为两个腔的湮灭(产生)算符, Δc=ω0ωc为驱动场与腔模的失谐; 第2项和第3项给出了旋转镜的自由哈密顿量, ωϕ, Lzϕ分别为旋转镜的角频率、角动量和角位移, I=MR2/2为旋转镜的转动惯量, 其中MR分别为旋转镜的质量和半径; 第4项和第5项表示L-G腔模与旋转镜的光旋转耦合, 单光子耦合强度分别为gϕ1=c1/Lgϕ2=c2/L, 其中L为两个腔的长度; 第6项表示两个腔模的线性耦合, 其中J为线性耦合强度; 最后4项分别表示强驱动场和弱探测场对腔模的驱动, δ=ωpωc为探测场与驱动场之间的失谐.

    为了方便计算, 设定b=Iωϕ2(ϕ+iIωϕLz)b=Iωϕ2(ϕiIωϕLz)分别表示转动腔镜的湮灭(产生)算符, 系统的哈密顿量可以写为

    H=Δc(c1c1+c2c2)+ωϕbb+(gϕ1c1c1+gϕ2c2c2)(b+b)+J(c1c2+c2c1)+i(εcc1εcc1)+i(εdc2εdc2)+iεL(c1eiδtc1eiδt)+iεR(c2eiδtc2eiδt).
    (2)

    由系统的哈密顿量, 可以写出描述系统的动力学的相关算符的海森伯-郎之万方程:

    ˙c1=[iΔc+κ12+igϕ1(b+b)]c1+εc+εLeiδtiJc2˙c2=[iΔc+κ22+igϕ2(b+b)]c2+εd+εReiδtiJc1˙b=iωmbγϕ2bi(gϕ1c1c1+gϕ2c2c2)
    (3)

    其中κjγϕ分别为两腔和转动腔镜的弛豫速率. 在乘积算符的平均值满足分解假设bcj=bcj且忽略探测场的情况下, 可得到算符的稳态解:

    bs=i(gϕ1|c1s|2+gϕ2|c2s|2)γϕ2+iωϕc1s=(κ22+iΔ2)εciJεdJ2+(κ12+iΔ1)(κ22+iΔ2)c2s=(κ12+iΔ1)εdiJεcJ2+(κ12+iΔ1)(κ22+iΔ2)
    (4)

    其中Δ1,2=Δc+gϕ1,ϕ2(bs+b*s)为L-G腔模和驱动场之间的有效失谐.

    在强驱动场的驱动下, 每个算符可以由其平均值和量子涨落的和表示, 即cj=cjs+δcj, b=bs+δb. 将cjb的表达式代入方程(3), 可以得到线性化的海森伯-郎之万方程:

    δ˙c1=(κ12+iΔ1)δc1igϕ1c1s(δb+δb)+εLeiδtiJδc2δ˙c2=(κ22+iΔ2)δc2igϕ2c2s(δb+δb)+εReiδtiJδc1δ˙b=(γϕ2+iωϕ)δbigϕ1(c1sδc1+c1sδc1)igϕ2(c2sδc2+c2sδc2)
    (5)

    采用变换δbδbeiωϕt, δc1δc1eiΔ1t, δc2δc2eiΔ2t, 则方程(5)为

    δ˙c1=κ12δc1iG1[δbei(ωϕ+Δ1)t+δbei(ωϕΔ1)t]+εLei(δΔ1)tiJδc2ei(Δ1Δ2)tδ˙c2=κ22δc2iG2eiθ[δbei(ωϕ+Δ2)t+δbei(ωϕΔ2)t]+εRei(δΔ2)tiJδc1ei(Δ2Δ1)tδ˙b=γϕ2δbiG1[δc1ei(ωϕΔ1)t+δc1ei(ωϕ+Δ1)t]iG2[eiθδc2ei(ωϕΔ2)t+eiθδc2ei(ωϕ+Δ2)t]
    (6)

    其中G1, G2为两个涡旋腔场与旋转腔镜之间的有效光旋转耦合强度, G1=gϕ1c1s, G2=gϕ2c2seiθ, gϕ1c1sgϕ2c2s之间的非互易相位差θ可以通过驱动场εcεd的相对相位来调节.

    如果每个强驱动场都处于力学红边带(Δ1Δ2ωϕ), 同时旋转镜频率ωϕ远大于有效耦合强度Gj, 方程(6)可以化简为

    δ˙c1=κ12δc1iG1δbiJδc2+εLeixtδ˙c2=κ22δc2iG2eiθδbiJδc1+εReixtδ˙b=γϕ2δbiG1δc1iG2eiθδc2
    (7)

    其中x=δωϕ. 由方程(7)的具体形式, 可以假设方程(7)的解具有下面的形式:

    δc1=δc1+eixt+δc1eixtδc2=δc2+eixt+δc2eixtδb=δb+eixt+δbeixt.
    (8)

    将方程(8)代入方程(7)得

    δc1+=(8G22+2γxκ2x)εL(8G1G2eiθ+4iJγx)εRΩ16iG1G2Jcosθδc2+=(8G21+2γxκ1x)εR(8G1G2eiθ+4iJγx)εLΩ16iG1G2Jcosθb+=(4iG1κ2x+8G2Jeiθ)εL+(8G1J+4iG2κ1xeiθ)εR16iG1G2JcosθΩ
    (9)

    其中γx=γ2ix, κjx=κj2ix, Ω=4G22κ1x+4G21κ2x+4J2γx+κ1xκ2xγx.

    为了研究涡旋光束传输的非互易性, 需要得到εoutLεoutR. 根据输入输出关系[47,48]

    εoutL+εinLeixt=κ1δc1εoutR+εinReixt=κ2δc2
    (10)

    其中εinL, R=εL, R/κ1κ2, 输出场具有δc1=δc1+eixt+δc1eixtδc2=δc2+eixt+δc2eixt的形式, 由方程(9)和方程(10)得

    εoutL + =κ1δc1+εL/κ1εoutR + =κ2δc2+εR/κ2εoutL=0  εout=0.
    (11)

    本文定义T12(T21)为入射腔c1(c2)的探测场通过系统后输出时的传输振幅, 根据(9)式和(11)式计算得

    T12=TLR=|εoutR+εinL|εinR=0=κ1κ2(8G1G2eiθ+4iJγx)Ω16iG1G2JcosθT21=TRL=|εoutL+εinR|εinL=0=κ1κ2(8G1G2eiθ+4iJγx)Ω16iG1G2Jcosθ.
    (12)

    本节将详细说明在双拉盖尔-高斯旋转腔系统中非互易相位差θ、涡旋光束所携带拓扑荷和系统耗散对传输振幅的影响. 首先讨论非互易相位差θ对传输振幅的影响.

    图2给出传输振幅T12(红色实线)和T21(蓝色虚线)在不同非互易相位差θ下随标准化失谐x/κ1的演化. 可以看出, 相位差θ对涡旋光束传输有比较大的影响, 在θ = 0时, 探测场在两个方向上的传输是互易的, 但当θ0探测场的传输会呈现明显的非互易性. 观察图2(b)图2(e)可以发现, 非互易现象会出现在共振点及附近, 还可以发现相位差的正负可决定非互易的方向(θ<0时, T21>T12; 而θ>0时, T21<T12). 当θ=±π /4时, 非互易性在x=0.4κ处呈现的最明显. 当θ=±π /2时, 系统可以在共振处(x=0)呈现完美的涡旋光束非互易传输(T12=0T21=1, 或T12=1T21=0). 由(12)式也可以看出, 当θ=nπ 时, T12=T21, 涡旋光束传输呈现互易性, 而当θnπ 时, T12T21, 即涡旋光束的传输是非互易的. 实际上, 这里的涡旋光束的非互易传输主要源于光子不同传输路径之间的量子干涉作用. 为了更清楚地展示更多非互易相位差对传输振幅的影响, 图2(f)给出了传输振幅T12(红色实线)和T21(蓝色虚线)随非互易相位差θ的演化.

    图 2 传输振幅$ {T_{12}} $(红色实线)和$ {T_{21}} $(蓝色虚线)在不同非互易相位差$ \theta $下随标准化失谐$ x/{\kappa _1} $的演化 (a) $ \theta  = 0 $; (b) $ \theta  =  - {\text{π }}/4 $; (c) $ \theta  = {\text{π }}/4 $; (d) $ \theta  =  - {\text{π }}/2 $; (e) $ \theta  = {\text{π }}/2 $. (f)在标准化失谐为零(x = 0)时, 传输振幅$ {T_{12}} $(红色实线)和$ {T_{21}} $(蓝色虚线)随非互易相位差$ \theta $的演化. 其他参数为$ {\kappa _1} = 1 $, $ {\kappa _2} = 1 $, $ \gamma  = 1 $,  $ {G_i} = \sqrt {{\kappa _i}\gamma } /2 $, $ J = \sqrt {{\kappa _1}{\kappa _2}} /2 $.\r\nFig. 2. Transmission amplitudes $ {T_{12}} $ (red solid line) and $ {T_{21}} $ (blue dotted line) versus normalized detuning $ x/{\kappa _1} $ under different nonreciprocal phase difference: (a) $ \theta  = 0 $; (b) $ \theta  =  - {\text{π }}/4 $; (c) $ \theta  = {\text{π }}/4 $; (d) $ \theta  =  - {\text{π }}/2 $; (e) $ \theta  = {\text{π }}/2 $. (f) Transmission amplitudes $ {T_{12}} $ (red solid line) and $ {T_{21}} $ (blue dotted line) versus nonreciprocal phase difference with $ x = 0 $. Other parameters are $ {\kappa _1} = 1 $, $ {\kappa _2} = 1 $, $ \gamma  = 1 $, $ {G_i} = \sqrt {{\kappa _i}\gamma } /2 $, $ J = \sqrt {{\kappa _1}{\kappa _2}} /2 $.
    图 2  传输振幅T12(红色实线)和T21(蓝色虚线)在不同非互易相位差θ下随标准化失谐x/κ1的演化 (a) θ=0; (b) θ=π /4; (c) θ=π /4; (d) θ=π /2; (e) θ=π /2. (f)在标准化失谐为零(x = 0)时, 传输振幅T12(红色实线)和T21(蓝色虚线)随非互易相位差θ的演化. 其他参数为κ1=1, κ2=1, γ=1, Gi=κiγ/2, J=κ1κ2/2.
    Fig. 2.  Transmission amplitudes T12 (red solid line) and T21 (blue dotted line) versus normalized detuning x/κ1 under different nonreciprocal phase difference: (a) θ=0; (b) θ=π /4; (c) θ=π /4; (d) θ=π /2; (e) θ=π /2. (f) Transmission amplitudes T12 (red solid line) and T21 (blue dotted line) versus nonreciprocal phase difference with x=0. Other parameters are κ1=1, κ2=1, γ=1, Gi=κiγ/2, J=κ1κ2/2.

    接下来讨论系统耗散对涡旋光束传输非互易性的影响. 为不失一般性, 取非互易相位差θ=π /2. 图3(a)图3(b)给出了传输振幅T12T21在旋转腔镜的不同耗散率γ下随标准化失谐x/κ1的演化. 可以看出在共振点(x = 0)附近, 随着旋转镜耗散率γ增加, 从腔c1到腔c2的传输振幅T12 = 1保持不变, 而从腔c2到腔c1的传输振幅T21呈现出明显的变化, 显示出先降低后增加的趋势. 当旋转镜耗散率γ = κ1 = 1时, 涡旋光束的传输非互易性达到最佳(T12=1, T21=0, 在共振点处).

    图 3 (a)传输振幅$ {T_{12}} $和(b)传输振幅$ {T_{{\text{21}}}} $作为标准化失谐$ x/{\kappa _1} $和旋转镜与腔$ {c_1} $耗散比$ \gamma /{\kappa _1} $的函数; (c)传输振幅$ {T_{12}} $和(d)传输振幅$ {T_{{\text{21}}}} $作为标准化失谐$ x/{\kappa _1} $和腔$ {c_2} $与腔$ {c_1} $耗散比$ {\kappa _2}/{\kappa _1} $的函数. 其他参数为$ \theta  = {\text{π }}/2 $, $ {\kappa _1} = 1 $, $ {G_1} = {G_2} = J = 1/2 $\r\nFig. 3. (a) Transmission amplitude $ {T_{12}} $ and (b) transmission amplitude $ {T_{21}} $ as a function of the normalized detuning $ x/{\kappa _1} $ and the dissipation ratio $ \gamma /{\kappa _1} $ of rotating mirror to cavity $ {c_1} $; (c) transmission amplitudes $ {T_{12}} $ and (d) transmission amplitudes $ {T_{21}} $ as a function of the normalized detuning $ x/{\kappa _1} $ and the dissipation ratio $ {\kappa _2}/{\kappa _1} $ of cavity $ {c_2} $ to cavity $ {c_1} $. Other parameters are $ \theta  = {\text{π }}/2 $, $ {\kappa _1} = 1 $, $ {G_1} = {G_2} = J = 1/2 $.
    图 3  (a)传输振幅T12和(b)传输振幅T21作为标准化失谐x/κ1和旋转镜与腔c1耗散比γ/κ1的函数; (c)传输振幅T12和(d)传输振幅T21作为标准化失谐x/κ1和腔c2与腔c1耗散比κ2/κ1的函数. 其他参数为θ=π /2, κ1=1, G1=G2=J=1/2
    Fig. 3.  (a) Transmission amplitude T12 and (b) transmission amplitude T21 as a function of the normalized detuning x/κ1 and the dissipation ratio γ/κ1 of rotating mirror to cavity c1; (c) transmission amplitudes T12 and (d) transmission amplitudes T21 as a function of the normalized detuning x/κ1 and the dissipation ratio κ2/κ1 of cavity c2 to cavity c1. Other parameters are θ=π /2, κ1=1, G1=G2=J=1/2.

    图3(c)图3(d)给出了传输振幅T12T21在腔c2的不同耗散率下随标准化失谐x/κ1的演化. 可以看出, 在共振点(x = 0)附近, 随着腔c2耗散率κ2增加, 从腔c2到腔c1的传输振幅T21 = 0保持不变, 而从腔c1到腔c2的传输振幅T12呈现逐渐减小的趋势. 当腔c2耗散率κ2 = κ1 = 1时, 涡旋光束传输的非互易性达到最佳(T12=1, T21=0, 在共振点处). 所以对于涡旋光束的非互易传输来说, 腔场c2的耗散不是越低越好.

    接下来讨论两个涡旋光场所携带拓扑荷比值2/1对传输振幅的影响. 通常, 涡旋光场所携带拓扑荷实际上会影响旋转腔镜与L-G腔模的耦合强度. 为了方便, 这里通过调节驱动场振幅εcεd的大小从而使得拓扑荷比2/1与有效耦合强度比G2/G1相等. 图4(a)图4(f)给出了在标准化失谐x = 0时传输振幅T12T21在不同的非互易相位差下随两个涡旋光场所携带拓扑荷比2/1的演化. 可以发现, 当θ=0θ=π 时, 拓扑荷的比值对传输振幅T12T21的影响是相同的, 即系统对涡旋光束的传输呈现互易性. 当θ=±π /4θ=±π /2时, 系统对涡旋光束传输呈现明显的非互易性. 当θ=±π /4时, 随着拓扑荷比值的增加, 传输振幅T12T21的变化趋势是相同的, 且是偶对称, 但相同拓扑荷比值所对应的传输振幅T12T21是不同的, 在拓扑荷比值为±1处的非互易性是最好的. 当θ=±π /2时, 随着拓扑荷比值的增加, 传输振幅T12仍呈现关于拓扑荷比值等于0处的对称性, 但此时传输振幅T21保持为0不变. 在拓扑荷比值为±1处, 系统对涡旋光束传输可以呈现完美的非互易性(T12=0T21=1, 或T12=1T21=0). 此外, 拓扑荷的正负代表涡旋光矢量的旋转方向, 根据涡旋光矢量的旋转方向可以将涡旋光束分为左旋涡旋光束和右旋涡旋光束. 从图4(d)图4(e)可以了解到涡旋光束想要呈现完美的非互易性与涡旋光束的左、右旋无关, 只要保证|2/1| = 1.

    图 4 传输振幅$ {T_{12}} $(红色圆圈)和$ {T_{21}} $(蓝色圆圈)在不同非互易相位差$ \theta $下随两个涡旋光场所携带拓扑荷之比$ {\ell _2}/{\ell _1} $的演化 (a) $\theta  =   0$; (b) $ \theta  =  - {\text{π }}/4 $; (c) $ \theta  = {\text{π }}/4 $; (d) $ \theta  =  - {\text{π }}/2 $; (e) $ \theta  = {\text{π }}/2 $; (f) $ \theta  = {\text{π }} $. 其他参数为$ {\kappa _1} = 1 $, $ {\kappa _2} = 1 $, $ \gamma  = 1 $, $ J = 2{G_1}{G_2}/\gamma $\r\nFig. 4. Transmission amplitudes $ {T_{12}} $ (red circle) and $ {T_{21}} $ (blue circle) with ratio of topological charges carried by two vortex optical fields $ {\ell _2}/{\ell _1} $ under different nonreciprocal phase difference: (a) $ \theta  = 0 $; (b)$ \theta  =  - {\text{π }}/4 $; (c)$ \theta  = {\text{π }}/4 $; (d) $ \theta  =  - {\text{π }}/2 $; (e) $ \theta  = {\text{π }}/2 $; (f) $ \theta  = {\text{π }} $. Other parameters are $ {\kappa _1} = 1 $, $ {\kappa _2} = 1 $, $ \gamma  = 1 $, $ J = 2{G_1}{G_2}/\gamma $.
    图 4  传输振幅T12(红色圆圈)和T21(蓝色圆圈)在不同非互易相位差θ下随两个涡旋光场所携带拓扑荷之比2/1的演化 (a) θ=0; (b) θ=π /4; (c) θ=π /4; (d) θ=π /2; (e) θ=π /2; (f) θ=π . 其他参数为κ1=1, κ2=1, γ=1, J=2G1G2/γ
    Fig. 4.  Transmission amplitudes T12 (red circle) and T21 (blue circle) with ratio of topological charges carried by two vortex optical fields 2/1 under different nonreciprocal phase difference: (a) θ=0; (b)θ=π /4; (c)θ=π /4; (d) θ=π /2; (e) θ=π /2; (f) θ=π . Other parameters are κ1=1, κ2=1, γ=1, J=2G1G2/γ.

    在腔系统中, 光的群延迟是描述光群速度的一个重要指标, 首先介绍群延迟的定义, 即[49,50]

    τij=dΘ21(12)dωp,
    (13)

    其中Θ21(12)是在频率为ωp下输出场c1c2(c2c1)的相位, 群延迟τij<0对应快光, τij>0对应慢光.

    图5(a)图5(c)给出了针对涡旋光束所对应的群延迟τ12τ21在不同的非互易相位差下随两个涡旋光场所携带拓扑荷比2/1的演化. 在θ=0时群速度是互易的且关于拓扑荷比值偶对称. 在θ=±π /2时, 系统可以呈现明显的非互易慢光效应, 且相位的变化可以影响非互易慢光的群延迟的变换. 随着拓扑荷比值的增加, τ12τ21均呈现先增大后减小的趋势且关于拓扑荷比值偶对称, 但群延迟的差值(τ12τ21τ21τ12)一直保持不变.

    图 5 群延迟$ {\tau _{12}} $(红色圆圈)和$ {\tau _{21}} $(蓝色圆圈)在不同非互易相位差$ \theta $下随两个涡旋光场所携带拓扑荷比$ {\ell _2}/{\ell _1} $的演化 (a) $\theta  =  0$; (b) $ \theta  =  - {\text{π }}/2 $; (c) $ \theta  = {\text{π }}/2 $. 其他参数为$ {\kappa _1} = 1 $, $ {\kappa _2} = 1 $, $ \gamma  = 1 $, $ J = 2{G_1}{G_2}/\gamma $\r\nFig. 5. Group delay $ {\tau _{12}} $ (red circle) and $ {\tau _{21}} $ (blue circle) with ratio of topological charges carried by two vortex optical fields $ {\ell _2}/{\ell _1} $ under different nonreciprocal phase difference: (a) $ \theta  = 0 $; (b) $ \theta  =  - {\text{π }}/2 $; (c) $ \theta  = {\text{π }}/2 $. Other parameters are $ {\kappa _1} = 1 $, $ {\kappa _2} = 1 $, $ \gamma  = 1 $, $ J = 2{G_1}{G_2}/\gamma $.
    图 5  群延迟τ12(红色圆圈)和τ21(蓝色圆圈)在不同非互易相位差θ下随两个涡旋光场所携带拓扑荷比2/1的演化 (a) θ=0; (b) θ=π /2; (c) θ=π /2. 其他参数为κ1=1, κ2=1, γ=1, J=2G1G2/γ
    Fig. 5.  Group delay τ12 (red circle) and τ21 (blue circle) with ratio of topological charges carried by two vortex optical fields 2/1 under different nonreciprocal phase difference: (a) θ=0; (b) θ=π /2; (c) θ=π /2. Other parameters are κ1=1, κ2=1, γ=1, J=2G1G2/γ.

    最后, 根据目前的实验进展, 对本文使用的系统方案的可行性进行了讨论. 本文讨论的双拉盖尔-高斯旋转腔系统是在文献[6]提出的由两个螺旋相位元件所构成的拉盖尔-高斯旋转腔的基础上, 添加一个固定的螺旋相位元件所实现的. 通过三个螺旋相位元件作为腔镜来组成两个光腔, 并且用光纤将两个光腔线性耦合. 其中, 两个输入腔镜FM1和 FM2为部分透明且被刚性固定, 不会改变透射光束的轨道角动量但会给反射光束增加2122的拓扑荷, 这一效应已经被实验观测到[51]. 另外, 处于中间位置的旋转腔镜RM是完全反射的, 可以给反射光增加+21+22的拓扑荷[6]. 随着纳米技术的发展, 螺旋相位元件可以利用超精密点金刚石车削机床直接加工铝盘表面得到[6,52]. 利用螺旋相位元件的反射和透射可以改变激光光束的方位结构, 进而改变其携带的拓扑荷的值. 拓扑荷依赖于螺旋相位元件的结构, 通过将螺旋相位元件的方位角坡道划分为离散的阶跃并控制阶跃高度和旋向, 可以设计出具有特定拓扑荷的涡旋光场[6,53]. 最近实验已经可以通过螺旋相位元件产生携带拓扑荷高达1000的L-G光束[54]. 对于本工作, 若想实现理想的涡旋光隔离器, 就需要在本系统中实现涡旋光束的完美非互易性, 即涡旋光束的单向传输. 在理论分析中得到, 若要呈现完美的非互易涡旋光束传输, 需要保证|2/1| = 1, 因此这里的涡旋光束所携带的拓扑荷需是可以调节的. 在实验上, 可以通过调整旋转腔镜两侧的阶跃高度和旋向来实现. 另外, 理论分析中得到非互易相位差对涡旋光束的非互易传输有着重要影响, 而非互易相位差的改变在实验上可以通过调节两个强驱动激光场εcεd的相对相位来实现[45]. 值得强调的是, 本文所构建的系统在实验实现时需要较高精细度的光腔, 在光腔的精细度比较低时, 在该系统所要探究相关现象会不明显[53].

    本文研究了涡旋光束在双拉盖尔-高斯旋转腔系统中的非互易传输. 将以往大多在光机械系统中研究的普通光束的非互易性推广到在光旋转系统中研究携带轨道角动量的涡旋光束传输的非互易性. 系统中, 光旋转耦合作用以及光纤的线性耦合作用提供的两条路径, 从而可以产生量子干涉. 首先, 研究了非互易相位差对传输振幅的影响, 发现相位差可以决定涡旋光束非互易性的发生及方向. 然后, 在相位差θ=π /2时, 研究了系统耗散对传输振幅的影响, 分析了涡旋光束实现完美的非互易传输时系统的耗散应该满足的条件. 接下来, 分析了涡旋光束携带的拓扑荷比值对系统非互易性的影响, 探究了实现涡旋光束完美的非互易传输时拓扑荷比值应该满足的条件, 发现涡旋光的左、右旋不会对传输产生影响. 拓扑荷对光传输产生影响, 主要是因为涡旋光束携带的拓扑荷数会影响旋转腔镜与腔模的耦合强度. 最后, 分析了该系统产生的非互易的慢光效应. 这些研究成果可用于设计针对携带轨道角动量的涡旋光束的理想光隔离器, 有望应用于光通信等领域.

    [1]

    Yankovich A B, Verre R, Olsén E, Persson A E O, Trinh V, Dovner G, Käll M, Olsson E 2017 ACS Nano 11 4265Google Scholar

    [2]

    Zeng C, Cui Y D 2013 Opt. Commun. 290 188Google Scholar

    [3]

    Huang L L, Chen X Z, Bai B F, Tan Q F, Jin G F, Zentgraf Z, Zhang S 2013 Light-Sci. Appl. 2 e70Google Scholar

    [4]

    Jankovic N, Cselyuszka N 2018 Sensors 18 1Google Scholar

    [5]

    Yan Z D, Wen X M, Gu P, Zhong H, Zhan P, Chen Z, Wang Z L 2017 Nanotechnol. 28 475203Google Scholar

    [6]

    Khatir M, Granpayeh N 2013 J. Lightwave Technol. 31 1045Google Scholar

    [7]

    陈颖, 罗佩, 田亚宁, 刘晓飞, 赵志勇, 朱奇光 2017 光学学报 37 0924002

    Chen Y, Luo P, Tian Y N, Liu X F, Zhao Z Y, Zhu Q G 2017 Acta Opt. Sin. 37 0924002

    [8]

    Fu H X, Li S L, Wang Y L, Song G, Zhang P F, Wang L L, Yu L 2018 IEEE Photonics J. 10 1

    [9]

    Chen Z, Yu L 2014 IEEE Photonics J. 6 1

    [10]

    Li C, Li S L, Wang Y L, Jiao R Z, Wang L L, Yu L 2017 IEEE Photonics J. 99 1

    [11]

    Wang D Q, Yu X L, Yu Q M 2013 Appl. Phys. Lett. 103 824

    [12]

    Artar A, Yanik A A 2011 Nano Lett. 11 3694Google Scholar

    [13]

    Wu C H, Khanikaev A, Shvets G 2011 Phys. Rev. Lett. 106 107403Google Scholar

    [14]

    Zhang Z D, Wang H Y, Zhang Z Y 2013 Plasmonics 8 797Google Scholar

    [15]

    Rakhshani M R, Mansouri-Birjandi M A 2016 IEEE Sens. J. 16 3041Google Scholar

    [16]

    Kim J, Soref R, Buchwald W R 2010 Opt. Express 18 17997Google Scholar

    [17]

    Zheng S, Ruan Z S, Gao S Q, Long Y, Li S M, He M G, Zhou N, Du J, Shen L, Cai X L, Wang J 2017 Opt. Express 25 25655Google Scholar

    [18]

    Guo Z C, Wen K H, Hu Q Y, Lai W H, Lin J Y, Fang Y H 2018 Sensors 18 1348Google Scholar

    [19]

    Lu H, Liu X, Mao D 2012 Phys. Rev. A 85 53803Google Scholar

    [20]

    Piao X J, Yu S, Koo S, Lee K 2011 Opt. Express 19 10907Google Scholar

    [21]

    Wu C, Ding H F, Huang T Y, Wu X, Chen B W, Ren K X, Fu S N 2017 Plasmonics 13 251

    [22]

    Wen K H, Hu Y H, Chen L, Zhou J Y, Liang L, Meng Z M 2016 Plasmonics 11 315Google Scholar

  • 图 1  结构示意图和透射光谱图 (a)含F-P腔的MDM波导耦合方形腔结构;(b)透射光谱;(c)相位图

    Fig. 1.  Schematic diagram and transmission spectrum: (a) MDM waveguide coupled square cavity structure with F-P cavity;(b) transmission spectrum; (c) phase diagram.

    图 2  Hz场分布 (a) FR1波峰处的Hz场分布;(b) FR1波谷处的Hz场分布;(c) FR2波峰处的Hz场分布;(d) FR2波谷处的Hz场分布

    Fig. 2.  The Hz field distribution: (a) The Hz field distribution at the peak of FR1; (b) the Hz field distribution at the dip of FR1; (c) the Hz field distribution at the peak of FR2; (d) the Hz field distribution at the dip of FR2.

    图 3  参数l1l2对传感特性的影响 (a)参数l1对FR1的影响;(b)参数l2对FR2的影响;(c)参数l1 = l2l1/l2对Fano共振线型的影响

    Fig. 3.  Influence of parameters l1 and l2 on sensing characteristics: (a) Influence of parameters l1 on the FR1; (b) influence of parameters l2 on the FR2; (c) influence of parameters l1 = l2 or l1/l2 on the Fano resonance.

    图 4  参数L1对传感特性的影响 (a)参数L1对FR1的影响;(b)参数L1对FR1的FOM值的影响;(c)参数L1对FR2的影响;(b)参数L1对FR2的FOM值的影响

    Fig. 4.  Influence of parameters L1 on sensing characteristics: (a) Influence of parameters L1 on the FR1; (b) influence of parameters L1 on the FOM value of FR1; (c) influence of parameters L1 on the FR2; (d) influence of parameters L1 on the FOM value of FR2.

    图 5  参数g1, g2和折射率n1, n2对传感特性的影响 (a)参数g1对FR1的影响;(b)参数g2对FR2的影响;(c)折射率n1对FR1的影响;(d)折射率n2对FR2的影响

    Fig. 5.  Influence of parameters g1, g2 and refractive index n1,n2 on sensing characteristics: (a) Influence of parameters g1 on the FR1; (b) influence of parameters g2 on the FR2; (c) influence of refractive index n1 on the FR1; (d) influence of refractive index n2 on the FR2.

  • [1]

    Yankovich A B, Verre R, Olsén E, Persson A E O, Trinh V, Dovner G, Käll M, Olsson E 2017 ACS Nano 11 4265Google Scholar

    [2]

    Zeng C, Cui Y D 2013 Opt. Commun. 290 188Google Scholar

    [3]

    Huang L L, Chen X Z, Bai B F, Tan Q F, Jin G F, Zentgraf Z, Zhang S 2013 Light-Sci. Appl. 2 e70Google Scholar

    [4]

    Jankovic N, Cselyuszka N 2018 Sensors 18 1Google Scholar

    [5]

    Yan Z D, Wen X M, Gu P, Zhong H, Zhan P, Chen Z, Wang Z L 2017 Nanotechnol. 28 475203Google Scholar

    [6]

    Khatir M, Granpayeh N 2013 J. Lightwave Technol. 31 1045Google Scholar

    [7]

    陈颖, 罗佩, 田亚宁, 刘晓飞, 赵志勇, 朱奇光 2017 光学学报 37 0924002

    Chen Y, Luo P, Tian Y N, Liu X F, Zhao Z Y, Zhu Q G 2017 Acta Opt. Sin. 37 0924002

    [8]

    Fu H X, Li S L, Wang Y L, Song G, Zhang P F, Wang L L, Yu L 2018 IEEE Photonics J. 10 1

    [9]

    Chen Z, Yu L 2014 IEEE Photonics J. 6 1

    [10]

    Li C, Li S L, Wang Y L, Jiao R Z, Wang L L, Yu L 2017 IEEE Photonics J. 99 1

    [11]

    Wang D Q, Yu X L, Yu Q M 2013 Appl. Phys. Lett. 103 824

    [12]

    Artar A, Yanik A A 2011 Nano Lett. 11 3694Google Scholar

    [13]

    Wu C H, Khanikaev A, Shvets G 2011 Phys. Rev. Lett. 106 107403Google Scholar

    [14]

    Zhang Z D, Wang H Y, Zhang Z Y 2013 Plasmonics 8 797Google Scholar

    [15]

    Rakhshani M R, Mansouri-Birjandi M A 2016 IEEE Sens. J. 16 3041Google Scholar

    [16]

    Kim J, Soref R, Buchwald W R 2010 Opt. Express 18 17997Google Scholar

    [17]

    Zheng S, Ruan Z S, Gao S Q, Long Y, Li S M, He M G, Zhou N, Du J, Shen L, Cai X L, Wang J 2017 Opt. Express 25 25655Google Scholar

    [18]

    Guo Z C, Wen K H, Hu Q Y, Lai W H, Lin J Y, Fang Y H 2018 Sensors 18 1348Google Scholar

    [19]

    Lu H, Liu X, Mao D 2012 Phys. Rev. A 85 53803Google Scholar

    [20]

    Piao X J, Yu S, Koo S, Lee K 2011 Opt. Express 19 10907Google Scholar

    [21]

    Wu C, Ding H F, Huang T Y, Wu X, Chen B W, Ren K X, Fu S N 2017 Plasmonics 13 251

    [22]

    Wen K H, Hu Y H, Chen L, Zhou J Y, Liang L, Meng Z M 2016 Plasmonics 11 315Google Scholar

  • [1] 陈召, 马昕新, 李童, 王艺霖. 耦合谐振系统中基于Fano共振的光学压力传感器. 物理学报, 2024, 73(8): 084205. doi: 10.7498/aps.73.20232025
    [2] 杨其利, 张兴坊, 刘凤收, 闫昕, 梁兰菊. 劈裂环-盘二聚体结构的多重Fano共振. 物理学报, 2022, 71(2): 027802. doi: 10.7498/aps.71.20210855
    [3] 陈颖, 周健, 丁志欣, 张敏, 朱奇光. 亚波长介质光栅/MDM波导/周期性光子晶体中双重Fano共振的形成及演变规律分析. 物理学报, 2022, 71(3): 034202. doi: 10.7498/aps.71.20211491
    [4] 杨其利, 张兴坊. 劈裂环-盘二聚体结构的多重Fano共振研究. 物理学报, 2021, (): . doi: 10.7498/aps.70.20210855
    [5] 鹿利单, 祝连庆, 曾周末, 崔一平, 张东亮, 袁配. 基于硅基光子器件的Fano共振研究进展. 物理学报, 2021, 70(3): 034204. doi: 10.7498/aps.70.20200550
    [6] 耿逸飞, 王铸宁, 马耀光, 高飞. 拓扑表面等离激元. 物理学报, 2019, 68(22): 224101. doi: 10.7498/aps.68.20191085
    [7] 张兴坊, 刘凤收, 闫昕, 梁兰菊, 韦德全. 同心椭圆柱-纳米管结构的双重Fano共振研究. 物理学报, 2019, 68(6): 067301. doi: 10.7498/aps.68.20182249
    [8] 李爱云, 张兴坊, 刘凤收, 闫昕, 梁兰菊. 对称纳米棒三聚体结构的Fano共振特性研究. 物理学报, 2019, 68(19): 197801. doi: 10.7498/aps.68.20190978
    [9] 刘仿, 李云翔, 黄翊东. 基于双表面等离子激元吸收的纳米光刻. 物理学报, 2017, 66(14): 148101. doi: 10.7498/aps.66.148101
    [10] 杨傅子. 从plasmon到nanoplasmonics——近代光子学前沿及液晶在其动态调制中的应用. 物理学报, 2015, 64(12): 124214. doi: 10.7498/aps.64.124214
    [11] 张志东, 高思敏, 王辉, 王红艳. 三角缺口正三角形纳米结构的共振模式. 物理学报, 2014, 63(12): 127301. doi: 10.7498/aps.63.127301
    [12] 吴青峻, 吴凡, 孙理斌, 胡晓琳, 叶鸣, 徐越, 史斌, 谢昊, 夏娟, 蒋建中, 张冬仙. 基于表面等离子激元的超薄金属减色滤波器的研究. 物理学报, 2014, 63(20): 207801. doi: 10.7498/aps.63.207801
    [13] 王五松, 张利伟, 冉佳, 张冶文. 微波频段表面等离子激元波导滤波器的实验研究. 物理学报, 2013, 62(18): 184203. doi: 10.7498/aps.62.184203
    [14] 黄茜, 张德坤, 熊绍珍, 赵颖, 张晓丹. 降低表面等离子激元寄生吸收损失的途径研究. 物理学报, 2012, 61(21): 217301. doi: 10.7498/aps.61.217301
    [15] 佟建波, 黄茜, 张晓丹, 张存善, 赵颖. 纳米Ag颗粒表面等离子激元对上转换材料光致发光性能影响的研究. 物理学报, 2012, 61(4): 047801. doi: 10.7498/aps.61.047801
    [16] 黄茜, 熊绍珍, 赵颖, 张晓丹. 表面等离子激元非线性表面增强拉曼散射效应. 物理学报, 2012, 61(15): 157801. doi: 10.7498/aps.61.157801
    [17] 郑俊娟, 孙刚. 嵌入电介质小球的金属薄片的电磁波透射特性. 物理学报, 2010, 59(6): 4008-4013. doi: 10.7498/aps.59.4008
    [18] 张存喜, 王瑞, 孔令民. 太赫兹场辅助的单量子阱自旋共振输运. 物理学报, 2010, 59(7): 4980-4984. doi: 10.7498/aps.59.4980
    [19] 宋国峰, 汪卫敏, 蔡利康, 郭宝山, 王青, 徐云, 韦欣, 刘运涛. 表面等离子激元调制的亚波长束斑半导体激光器. 物理学报, 2010, 59(7): 5105-5109. doi: 10.7498/aps.59.5105
    [20] 郑俊娟, 孙 刚. 周期排列的电介质小球所诱发的金属-电介质表面上的表面等离子激元的光学性质. 物理学报, 2005, 54(11): 5210-5217. doi: 10.7498/aps.54.5210
计量
  • 文章访问数:  8454
  • PDF下载量:  86
出版历程
  • 收稿日期:  2018-11-07
  • 修回日期:  2019-02-27
  • 上网日期:  2019-05-01
  • 刊出日期:  2019-05-20

/

返回文章
返回