搜索

x

留言板

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

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

一种低损耗的对称双楔形太赫兹混合表面等离子体波导

王芳 张龙 马涛 王旭 刘玉芳 马春旺

引用本文:
Citation:

一种低损耗的对称双楔形太赫兹混合表面等离子体波导

王芳, 张龙, 马涛, 王旭, 刘玉芳, 马春旺

A symmetrical wedge-to-wedge THz hybrid SPPs waveguidewith low propagation loss

Wang Fang, Zhang Long, Ma Tao, Wang Xu, Liu Yu-Fang, Ma Chun-wang
PDF
HTML
导出引用
  • 本文研究了一种适用于太赫兹频段的对称双楔形混合等离子体波导. 采用有限元法对其混合模式特征进行数值模拟, 阐述了波导的传播长度、归一化有效模场面积和品质因子等特性随波导几何参数的变化规律. 结果表明, 对称双楔形混合等离子体波导在太赫兹频段可获得良好的模场约束能力和低损耗特性, 在有效模场面积达到λ2/10280时, 传播长度达到51 × 103 μm. 波导参数最优时的平行对称楔形波导在忽略波导间串扰时, 波导间距小于16 μm时的耦合长度约为8958 μm. 通过对比发现, 相比于先前的对称微楔形波导结构和对称蝴蝶结波导结构, 对称双楔形混合等离子体波导在传输特性和耦合特性方面都具有更大的优势, 将在光学力捕获、生物分子传输以及高密度集成电路设计等方面具有潜在的应用前景.
    A symmetrical wedge-to-wedge THz hybrid SPPs waveguide (WWTHSW) with low propagation loss is investigated. The WWTHSW consists of two identical dielectric wedge waveguides symmetrically placed on each side of a micro wedge-patterned thin metal film. The mode characteristics of the WWTHSW, such as the propagation length (Lp), the normalized effective mode area (A) and the figure of merit (FOM) are analyzed by using the finite element method (FEM). Firstly, the influences of the height of Si micro wedge waveguide (H) and the gap between two wedges (g) on Lp and A are studied. For the same g, A first decreases and then increases with the increase of H. A achieves a minimum at an H of ~40 μm. However, Lp monotonically increases as H increases. The change of Lp slows down when H is greater than 40 μm. At a fixed H, Lp slightly increases with the increase of g. But A achieves a minimum when g is ~50 nm. Secondly, the dependencies of the mode characteristics of the WWTHSW on Si wedge tip angle (α) and Ag wedge tip angle (θ) are analyzed. At a fixed α, θ has less effect on Lp and A. As α increases at a fixed θ, Lp increases monotonically but A decreases firstly and then increases. A reaches a minimum when α increases to ~100°. Then, the change of Lp and A with the thicknesses of Ag film (d) and Ag wedge (h) are demonstrated. At a fixed h, both Lp and A slightly decrease as d increases. For the same d, Lp and A decrease with the increase of h. A for h = 0 μm is distinctly larger than those for h = 2 μm and h = 5 μm. According to the above optimizations, the parameters of the WWTHSW are chosen as d = 100 nm, g = 50 nm, h = 2 μm, θ = 80°, α = 100°, H = 40 μm. Under the optimal parameters, Lp of ~51 mm is obtained when Am reaches ~λ2/10280. Compared with the previous hybrid THz plasmonic waveguide, Lp of the WWTHSW increases by 3 times, and A decreases by an order of magnitude. This result reveals that the WWTHSW enables low-loss propagation and ultra-deep-subwavelength mode confinement at THz frequencies. At last, the coupling property of the parallel WWTHSW is investigated. The coupling length of ~8958 μm is achieved without the crosstalk between two parallel waveguides. By comparison, the WWTHSW has more advantages in terms of transmission and coupling characteristics than the previous micro wedge waveguide structure and bow-tie waveguide structure. In summary, due to the excellent transmission and coupling characteristics, the WWTHSW has great potential in the fields of optical force in trapping, biomolecules transporting, and in high-density integrated circuits design.
      通信作者: 马涛, matao_bupt79@163.com
    • 基金项目: 国家级-国家自然科学基金(61627818)
      Corresponding author: Ma Tao, matao_bupt79@163.com
    [1]

    Andersen J, Solodukhov V 1978 IEEE T. Antenn. Propag. 26 598Google Scholar

    [2]

    Oulton R F, Sorger V J, Genov D A, Pile D F P, Zhang X J 2008 Nature Photon. 2 496Google Scholar

    [3]

    Li Qing, Pan D, Wei H, Xu H X 2018 Nano Lett. 18 2009Google Scholar

    [4]

    Bozhevolnyi S I, Volkov V S, Devaux, E, Laluet, J Y, Ebbesen, T W 2006 Nature 440 508Google Scholar

    [5]

    Zhu Z H, Garcia-Ortiz C E, Han Z H, Radko I P, Bozhevolnyi S I 2013 Appl. Phys. Lett. 103 061108

    [6]

    Hassan K, Leroy F, Colas-Des-Francs G, Weeber, J C 2014 Opt. Lett. 39 697Google Scholar

    [7]

    Papaioannou S, Giannoulis G, Vyrsokinos K, Leroy F, Zacharatos F, Markey L, Weeber C J, Dereux A, Bozhevolnyi S, Prinzen A, Apostolopoulos D, Avramopoulos H, Pleros N 2015 IEEE Photon. Technol. Lett. 27 963Google Scholar

    [8]

    Hui F, Berini P 2016 J. Lightwave Technol. 34 2631Google Scholar

    [9]

    涂鑫, 陈震旻, 付红岩 2019 物理学报 68 104210Google Scholar

    Tu X, Chen Z M, Fu H Y 2019 Acta Phys. Sin. 68 104210Google Scholar

    [10]

    Wang Y L, Li S L, Yan J Y, Li C, Jiang P, Wang L L, Yu L 2019 Nanophotonics 8 1271Google Scholar

    [11]

    Wang Y L, Li C, Duan G Y, Wang L L, Yu L 2019 Adv. Opt. Mater. 7 1801362Google Scholar

    [12]

    Mai W, Wang Y, Zhang Y, Cui L, Yu L 2017 Chin. Phys. Lett. 34 024204Google Scholar

    [13]

    Gosciniak J, Volkov V S, Bozhevolnyi S I, Markey L, Massenot S, Dereux A 2010 Opt. Express 18 5314Google Scholar

    [14]

    Berini P 2009 Adv. Opt. Photonic 1 484Google Scholar

    [15]

    Ma Y Q, Gerald F, Yuliya S, Wu Q 2014 Opt. Lett. 39 973Google Scholar

    [16]

    Pan M Y, Lin E H, Wang L, Wei P K 2014 Appl. Phys. A. 115 592

    [17]

    Ma Y Q, Farrell G, Semenova Y, Wu Q 2015 J. Lightwave Technol. 33 3827Google Scholar

    [18]

    贾智鑫, 段欣, 吕婷婷, 郭亚楠, 薛文瑞 2011 物理学报 60 057301Google Scholar

    Jia Z X, Du X, Lü T T, Guo Y N, Xue W R 2011 Acta Phys. Sin. 60 057301Google Scholar

    [19]

    Gong Q, Bian Y 2014 J. Lightwave Technol. 32 4504Google Scholar

    [20]

    Jacek G, Volkov V S, Bozhevolnyi S I, Markey L, Massenot S, Dereux A 2010 Optics Express 18 5314

    [21]

    彭滟, 施辰君, 朱亦鸣, 庄松林 2019 中国激光 46 0614002Google Scholar

    Peng Y, Shi C J, Zhu Y M, Zhuang S L 2019 Chin. J. Las. 46 0614002Google Scholar

    [22]

    陈华, 汪力 2009 物理学报 58 4605Google Scholar

    Chen H, Wang L 2009 Acta Phys. Sin. 58 4605Google Scholar

    [23]

    Eldlio M, Ma Y Q, Maeda H, Cada M 2017 Infrared Phys. Technol. 80 93Google Scholar

    [24]

    Cao W, Song C Y, Lanier T E, Singh R, O’Hara J F, Dennis W M, Zhao Y P, Zhang W L 2013 Sci. Rep. 3 1766Google Scholar

    [25]

    Fitch M J, Osiander R 2004 J. Hopkins Apl. Technol. D. 25 348

    [26]

    Berini P 2006 Opt. Express 14 13030Google Scholar

    [27]

    Moreno E, Rodrigo S G, Bozhevolnyi S I, Martín-Moreno L, García-Vidal F J 2008 Phys. Rev. Lett. 100 023901Google Scholar

    [28]

    Ma Y Q, Farrell G, Semenova Y, Chan H P, Zhang H Z, Wu Q 2013 Plasmonics 8 1259Google Scholar

    [29]

    Georgios V, Shanhui F 2008 Opt. Express 16 2129Google Scholar

  • 图 1  WWTHSW示意图 (a)三维图; (b)截面图

    Fig. 1.  Schematic diagram of the proposed WWTHSW: (a) 3D diagram; (b) cross-section.

    图 2  不同Hg时, WWTHSW的模式分析 (a) MC, (b) Lp, (c) A; 模场分布: (d) [H, g] = [10, 0.05] μm, (e) [H, g] = [40, 0.05] μm, (f) [H, g] = [90, 0.05] μm

    Fig. 2.  Modes analysis of the WWTHSW with different H and g: (a) MC, (b) Lp, and (c) A;and normalized EM energy density distributions: (d) [H, g] = [10, 0.05] μm; (e) [H, g] = [40, 0.05] μm; (f) [H, g] = [90, 0.05] μm.

    图 3  不同αθ时, WWTHSW的模式分析 (a) Lp, (b) A; (c)模场分布随α的变化(θ = 80°); (d)模场分布随θ的变化(α = 100°)

    Fig. 3.  Modes analysis of the WWTHSW with different α and θ, (a) Lp, (b) A; and normalized EM energy density distributions: (c) with different α at a fixed θ of 80°, (d) with different θ at a fixed α of 100°.

    图 4  不同dh时, WWTHSW的模式分析 (a) Lp、(b) A; (c)沿x方向的归一化能量密度

    Fig. 4.  Modes analysis of the WWTHSW with different d and h, (a) Lp, (b) A, and (c) normalized EM energy density.

    图 5  不同波导性能比较 (a) WWTHSW, HTMWSPPs和HTBTSPPs波导的截面图; (b) WWTHSW, HTMWSPPs和HTBTSPPs波导的ALp关系图; (c)品质因数

    Fig. 5.  Performance comparison of the WWTHSW, HTMWSPPs and HTBTSPPs wavguide: (a) cross-section views; (b) the relationship between A and Lp; and (c) FOM with different parameters.

    图 6  波导耦合特性分析 (a)平行波导三维结构示意图; (b)耦合长度随D的变化; (c)最大传输功率随D的变化

    Fig. 6.  Coupling characteristic of waveguides: (a) schematic diagram parallel waveguides; (b) Lc versus the separation between the two waveguides; (c) the maximum transfer power (Pmax) as a function of distance D.

  • [1]

    Andersen J, Solodukhov V 1978 IEEE T. Antenn. Propag. 26 598Google Scholar

    [2]

    Oulton R F, Sorger V J, Genov D A, Pile D F P, Zhang X J 2008 Nature Photon. 2 496Google Scholar

    [3]

    Li Qing, Pan D, Wei H, Xu H X 2018 Nano Lett. 18 2009Google Scholar

    [4]

    Bozhevolnyi S I, Volkov V S, Devaux, E, Laluet, J Y, Ebbesen, T W 2006 Nature 440 508Google Scholar

    [5]

    Zhu Z H, Garcia-Ortiz C E, Han Z H, Radko I P, Bozhevolnyi S I 2013 Appl. Phys. Lett. 103 061108

    [6]

    Hassan K, Leroy F, Colas-Des-Francs G, Weeber, J C 2014 Opt. Lett. 39 697Google Scholar

    [7]

    Papaioannou S, Giannoulis G, Vyrsokinos K, Leroy F, Zacharatos F, Markey L, Weeber C J, Dereux A, Bozhevolnyi S, Prinzen A, Apostolopoulos D, Avramopoulos H, Pleros N 2015 IEEE Photon. Technol. Lett. 27 963Google Scholar

    [8]

    Hui F, Berini P 2016 J. Lightwave Technol. 34 2631Google Scholar

    [9]

    涂鑫, 陈震旻, 付红岩 2019 物理学报 68 104210Google Scholar

    Tu X, Chen Z M, Fu H Y 2019 Acta Phys. Sin. 68 104210Google Scholar

    [10]

    Wang Y L, Li S L, Yan J Y, Li C, Jiang P, Wang L L, Yu L 2019 Nanophotonics 8 1271Google Scholar

    [11]

    Wang Y L, Li C, Duan G Y, Wang L L, Yu L 2019 Adv. Opt. Mater. 7 1801362Google Scholar

    [12]

    Mai W, Wang Y, Zhang Y, Cui L, Yu L 2017 Chin. Phys. Lett. 34 024204Google Scholar

    [13]

    Gosciniak J, Volkov V S, Bozhevolnyi S I, Markey L, Massenot S, Dereux A 2010 Opt. Express 18 5314Google Scholar

    [14]

    Berini P 2009 Adv. Opt. Photonic 1 484Google Scholar

    [15]

    Ma Y Q, Gerald F, Yuliya S, Wu Q 2014 Opt. Lett. 39 973Google Scholar

    [16]

    Pan M Y, Lin E H, Wang L, Wei P K 2014 Appl. Phys. A. 115 592

    [17]

    Ma Y Q, Farrell G, Semenova Y, Wu Q 2015 J. Lightwave Technol. 33 3827Google Scholar

    [18]

    贾智鑫, 段欣, 吕婷婷, 郭亚楠, 薛文瑞 2011 物理学报 60 057301Google Scholar

    Jia Z X, Du X, Lü T T, Guo Y N, Xue W R 2011 Acta Phys. Sin. 60 057301Google Scholar

    [19]

    Gong Q, Bian Y 2014 J. Lightwave Technol. 32 4504Google Scholar

    [20]

    Jacek G, Volkov V S, Bozhevolnyi S I, Markey L, Massenot S, Dereux A 2010 Optics Express 18 5314

    [21]

    彭滟, 施辰君, 朱亦鸣, 庄松林 2019 中国激光 46 0614002Google Scholar

    Peng Y, Shi C J, Zhu Y M, Zhuang S L 2019 Chin. J. Las. 46 0614002Google Scholar

    [22]

    陈华, 汪力 2009 物理学报 58 4605Google Scholar

    Chen H, Wang L 2009 Acta Phys. Sin. 58 4605Google Scholar

    [23]

    Eldlio M, Ma Y Q, Maeda H, Cada M 2017 Infrared Phys. Technol. 80 93Google Scholar

    [24]

    Cao W, Song C Y, Lanier T E, Singh R, O’Hara J F, Dennis W M, Zhao Y P, Zhang W L 2013 Sci. Rep. 3 1766Google Scholar

    [25]

    Fitch M J, Osiander R 2004 J. Hopkins Apl. Technol. D. 25 348

    [26]

    Berini P 2006 Opt. Express 14 13030Google Scholar

    [27]

    Moreno E, Rodrigo S G, Bozhevolnyi S I, Martín-Moreno L, García-Vidal F J 2008 Phys. Rev. Lett. 100 023901Google Scholar

    [28]

    Ma Y Q, Farrell G, Semenova Y, Chan H P, Zhang H Z, Wu Q 2013 Plasmonics 8 1259Google Scholar

    [29]

    Georgios V, Shanhui F 2008 Opt. Express 16 2129Google Scholar

  • [1] 马赛群, 邓奥林, 吕博赛, 胡成, 史志文. 低维材料极化激元及其耦合特性. 物理学报, 2022, 71(12): 127104. doi: 10.7498/aps.71.20220272
    [2] 马涛, 马家赫, 刘恒, 田永生, 刘少晖, 王芳. 一种电光可调的铌酸锂/钠基表面等离子体定向耦合器. 物理学报, 2022, 71(5): 054205. doi: 10.7498/aps.71.20211217
    [3] 张利胜. 基于金纳米阵列表面等离子体驱动的光催化特性. 物理学报, 2021, 70(23): 235202. doi: 10.7498/aps.70.20210424
    [4] 王芳, 陈亚珂, 李传强, 马涛, 卢颖慧, 刘恒, 金婵. 非对称银膜多孔硅-氟化钙等离子体波导及其波导灵敏度特性. 物理学报, 2021, 70(22): 224201. doi: 10.7498/aps.70.20210704
    [5] 马昊军, 王国林, 罗杰, 刘丽萍, 潘德贤, 张军, 邢英丽, 唐飞. S-Ka频段电磁波在等离子体中传输特性的实验研究. 物理学报, 2018, 67(2): 025201. doi: 10.7498/aps.67.20170845
    [6] 李丹, 刘勇, 王怀兴, 肖龙胜, 凌福日, 姚建铨. 太赫兹波段石墨烯等离子体的增益特性. 物理学报, 2016, 65(1): 015201. doi: 10.7498/aps.65.015201
    [7] 李志全, 张明, 彭涛, 岳中, 顾而丹, 李文超. 基于导模共振效应提高石墨烯表面等离子体的局域特性. 物理学报, 2016, 65(10): 105201. doi: 10.7498/aps.65.105201
    [8] 孙杰, 杨剑锋, 闫肃, 杨晶晶, 黄铭. 等离子体辅助平板波导的传输特性及应用研究. 物理学报, 2015, 64(7): 078402. doi: 10.7498/aps.64.078402
    [9] 陈文波, 龚学余, 邓贤君, 冯军, 黄国玉. THz电磁波在时变非磁化等离子体中的传播特性研究. 物理学报, 2014, 63(19): 194101. doi: 10.7498/aps.63.194101
    [10] 刘亚青, 张玉萍, 张会云, 吕欢欢, 李彤彤, 任广军. 光抽运多层石墨烯太赫兹表面等离子体增益特性的研究. 物理学报, 2014, 63(7): 075201. doi: 10.7498/aps.63.075201
    [11] 张利伟, 赵玉环, 王勤, 方恺, 李卫彬, 乔文涛. 各向异性特异材料波导中表面等离子体的共振性质. 物理学报, 2012, 61(6): 068401. doi: 10.7498/aps.61.068401
    [12] 胡海峰, 蔡利康, 白文理, 张晶, 王立娜, 宋国峰. 基于表面等离子体的太赫兹光束方向调控的模拟研究. 物理学报, 2011, 60(1): 014220. doi: 10.7498/aps.60.014220
    [13] 邹文康, 陈林, 周良骥, 王勐, 杨礼兵, 谢卫平, 邓建军. Z箍缩驱动器与丝阵负载耦合特性研究. 物理学报, 2011, 60(11): 115204. doi: 10.7498/aps.60.115204
    [14] 程木田. 经典光场相干控制金属纳米线表面等离子体传输. 物理学报, 2011, 60(11): 117301. doi: 10.7498/aps.60.117301
    [15] 李山, 钟明亮, 张礼杰, 熊祖洪, 张中月. 偏振方向及结构间耦合作用对空心方形银纳米结构表面等离子体共振的影响. 物理学报, 2011, 60(8): 087806. doi: 10.7498/aps.60.087806
    [16] 方春易, 张树仁, 卢俊, 汪剑波, 孙连春. 一种圆孔单元厚屏频率选择表面结构的传输特性研究. 物理学报, 2010, 59(7): 5023-5027. doi: 10.7498/aps.59.5023
    [17] 陈华, 汪力. 金属导线偶合THz表面等离子体波. 物理学报, 2009, 58(7): 4605-4609. doi: 10.7498/aps.58.4605
    [18] 周仁龙, 陈效双, 曾 勇, 张建标, 陈洪波, 王少伟, 陆 卫, 李宏建, 夏 辉, 王玲玲. 金属光子晶体平板的超强透射及其表面等离子体共振. 物理学报, 2008, 57(6): 3506-3513. doi: 10.7498/aps.57.3506
    [19] 花 磊, 宋国峰, 郭宝山, 汪卫敏, 张 宇. 中红外下半导体掺杂调制的表面等离子体透射增强效应. 物理学报, 2008, 57(11): 7210-7215. doi: 10.7498/aps.57.7210
    [20] 高建霞, 宋国峰, 郭宝山, 甘巧强, 陈良惠. 表面等离子体调制的纳米孔径垂直腔面发射激光器. 物理学报, 2007, 56(10): 5827-5830. doi: 10.7498/aps.56.5827
计量
  • 文章访问数:  3782
  • PDF下载量:  65
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-31
  • 修回日期:  2019-12-24
  • 刊出日期:  2020-04-05

一种低损耗的对称双楔形太赫兹混合表面等离子体波导

  • 1. 河南师范大学电子与电气工程学院, 新乡 453007
  • 2. 河南省光电传感集成应用重点实验室, 新乡 453007
  • 3. 河南师范大学物理学院, 新乡 453007
  • 通信作者: 马涛, matao_bupt79@163.com
    基金项目: 国家级-国家自然科学基金(61627818)

摘要: 本文研究了一种适用于太赫兹频段的对称双楔形混合等离子体波导. 采用有限元法对其混合模式特征进行数值模拟, 阐述了波导的传播长度、归一化有效模场面积和品质因子等特性随波导几何参数的变化规律. 结果表明, 对称双楔形混合等离子体波导在太赫兹频段可获得良好的模场约束能力和低损耗特性, 在有效模场面积达到λ2/10280时, 传播长度达到51 × 103 μm. 波导参数最优时的平行对称楔形波导在忽略波导间串扰时, 波导间距小于16 μm时的耦合长度约为8958 μm. 通过对比发现, 相比于先前的对称微楔形波导结构和对称蝴蝶结波导结构, 对称双楔形混合等离子体波导在传输特性和耦合特性方面都具有更大的优势, 将在光学力捕获、生物分子传输以及高密度集成电路设计等方面具有潜在的应用前景.

English Abstract

    • 表面等离子体学是纳米光电子领域非常重要的研究方向, 表面等离子体激元(surface plasmonpolaritons, SPPs)因其亚波长尺寸传输和局域场增强效应的优点已成为研究的热点[1-3]. 表面等离子体波导具有能够突破衍射极限、整体尺寸较小且制作方便等优势, 因此备受研究者的关注[4], 可广泛应用于耦合器[5]、偏振转换器[6]、M-Z干涉结构[7,8]、光开关[9]、纳米天线[10,11]以及折射率传感器[12]等光子和光电器件的设计. 为了提高SPPs波的传输距离, 可以用介质波导替代SPPs波导(即减小SPPs波导的长度)来实现低损耗传输[13]. 在介质芯层中增加金属薄层[14]、利用混合SPPs波导[15]和双层介质加载SPPs波导[16]等方法也可以增加SPPs波的传输距离. 近年来, 混合表面等离子体波导(hybrid plasmonicwaveguides, HPWs)以其独特的突破衍射极限和在亚波长尺度下提供紧密光约束的能力, 被认为是大规模光子集成器件最有前途的候选者之一[4]. 目前, HPWs波导有圆柱形[2]、楔形[17]、领结形[18]及蝴蝶结(bow-tie, BT)[19]等多种形状. HPWs波导的设计不仅降低光在其中的传播损耗, 还进一步提高对光场的限制能力, 在超高密度光子集成器件和高灵敏度传感等领域具有很大的潜在应用.

      太赫兹(terahertz, THz)波具有许多优异特性, 比如: 能量需求低、穿透衰减小、信噪比高、对生物大分子拥有强吸收和谐振特性等[20]. 与其他波段的电磁波相比, THz波具有直接探测生物分子信息的能力, 它在药品生产的监控、生物大分子的光谱识别[21]等方面具有巨大的应用潜力. 目前, 对支持混合等离子体模式的混合等离子体太赫兹波导(hybrid plasmonic terahertz waveguide, HPTWG)的研究具有良好的应用前景. 由于THz波的频率占据微波和红外之间, HPTWG可以实现低损耗微米级传播及较高的模场限制[22]. 虽然长程太赫兹混合等离子体(long-rang hybrid THz SPPs, LR-HTSPPs)波导减少了传播损耗, 实现了厘米级的传播长度, 但其束缚能力较弱[23].

      为了更好地权衡模式束缚和传播长度之间的相互关系, 本文研究了一种对称双楔形太赫兹混合表面等离子体波导(wedge-to-wedge THz hybrid SPPs waveguide, WWTHSW)结构. 利用有限元方法(finite element method, FEM)研究了WWTHSW的模式特性和传输特性及其随波导结构尺寸变化的规律, 并分析了两条平行WWTHSW间的耦合特性. 最后, 我们对比了WWTHSW与对称微楔形太赫兹混合等离子体(hybrid THz microwedgeSPPs, HTMWSPPs)波导和对称型蝴蝶结太赫兹混合等离子体(hybrid THz bow-tie SPPs, HTBTSPPs)波导的性能.

    • WWTHSW的三维结构示意图和截面图如图1(a)图1(b)所示, 由两个对称放置的完全相同的硅楔形波导(silicon wedge waveguides, SWWs)和夹在中间的Ag层构成, 两者间隔为g. 波导的包层材料选用氮化硅(Si3N4). SWWs的高度为H, 楔形尖端采用圆角化处理, 曲率半径为r, 圆角为α. Ag层的厚度为d = 100 nm, 楔形尖端曲率半径也为r, 圆角为$\theta $. 曲率半径设定为r = 2 μm. Si, Si3N4折射率分别为3.455和2. 根据Drude模型[24,25], Ag的介电常数可表示为

      图  1  WWTHSW示意图 (a)三维图; (b)截面图

      Figure 1.  Schematic diagram of the proposed WWTHSW: (a) 3D diagram; (b) cross-section.

      $ {\varepsilon _{\rm{r}}} = {\varepsilon _\infty } - \frac{{\omega _{\rm{p}}^{\rm{2}}}}{{\omega \left( {\omega + {\rm{j}}\gamma } \right)}}, $

      其中, ε为高频介电常数, ω为入射光角频率, γ为阻尼项, ωp为等离子体角频率.

      $ {\omega _{\rm{p}}} = \sqrt {\frac{{n{e^{\rm{2}}}}}{{{\varepsilon _{\rm{0}}}{m^*}}}}, $

      其中, m*, en分别为电子有效质量、电子电荷和载流子密度; ε0为空间的介电常数. 参数ωp = 1.37 × 106 rad/s, γ = 8.20 × 1013 rad/s[24].

      为了定量分析HPTWG的模式和传输特性, 通常用有效折射率(neff)、传播长度(Lp)、归一化模场面积[17](A)、品质因数[26] (figure of merit, FOM)以及模式分布特性[2] (the mode character, MC)等参数衡量其特性.

      传播长度Lp

      ${L_{\rm{p}}} = {\lambda / {(4{\text{π}}{\rm{Im}}({n_{{\rm{eff}}}}))}},$

      其中, λ是真空中的波长, Im(neff)是模式有效折射率的虚部.

      A定义为Am/A0. 其中: A0 = λ2/4表示自由空间中的衍射限制面积; Am为有效模场面积, 反映波导对模式的束缚能力, 定义为

      $\begin{split} {A_{\rm{m}}} \; &= \frac{{{W_{\rm{m}}}}}{{\max \left\{ {W\left( {x,y} \right)} \right\}}} \\ &= \frac{1}{{\max \left\{ {W\left( {x,y} \right)} \right\}}}\iint {W\left( {x,y} \right){\rm{d}}x{\rm{d}}y},\end{split}$

      其中, Wm是模式的总能量, W(x, y)是能量密度.

      为综合衡量波导的传输特性, 定义品质因数为

      ${\rm{FOM}} = {{{\rm{Re}}({n_{{\rm{eff}}}})}/ {{\rm{Im}}({n_{{\rm{eff}}}})}},$

      其中, Re(neff)是模式有效折射率的实部.

      为描述混合等离子体波导模式中SPPs模式和波导模式之间的差异, 引入模式特性(MC):

      ${\rm{MC }}=\frac{{{n_{{\rm{eff:HSPPs}}}} - {n_{{\rm{eff:SPPs}}}}}}{{\left( {{n_{{\rm{eff:HSPPs}}}} - {n_{{\rm{eff:SPPs}}}}} \right) + \left( {{n_{{\rm{eff:HSPPs}}}} - {n_{{\rm{eff:W}}}}} \right)}},$

      其中neff:HSPPs, neff:SPPsneff:W分别表示混合模式、SPPs模式和波导模式的有效折射率. MC可用于评估波导模式特性, 当neff:W < neff:SPPs, 即MC < 0.5, 为类似表面等离子体模式, 当neff:W > neff:SPPs, 即MC > 0.5, 为类似波导模式.

      本文采用FEM对太赫兹频段的模式特性和传播特性进行研究. 使用商用ComsolMultiphysics软件中的波动光学模块, 模型采用散射边界条件, 仿真时的模型外层区域宽度均设置为2λ, 应用非均匀三角网格进行剖分并设置最小单元网格大小为d/2, 由于Ag层内几乎没有场分布, 因此, d/2能够保证模型计算的收敛性的精度.

    • 为了优化WWTHSW的性能, 利用Comsol Multiphysics软件中的模式分析研究不同波导几何参数下波导中存在的HPTWG的模式特性. 以下仿真中, 入射光的频率设置为1 THz. 楔形波导的高度H与间隙g对MC的影响如图2(a)所示, 当H小于40 μm时, MC < 0.5, 波导模式以SPPs模式为主; 当H大于40 μm时, MC > 0.5, 波导模式以SWWs模式为主. 图2(d)(f)所示, H = 10 μm时, 能量主要集中在间隙中Ag与Si3N4的交界面附近; 当H = 40 μm时, 能量依然主要集中在间隙中, 然而有部分能量向Si基楔形波导转移; 当H = 90 μm时, 模场的能量主要分布在Si基楔形波导内部. 为了优化设计参数, 我们研究了波导的Hg对混合对称双楔形波导LpA的影响, 如图2(b)(c)所示. WWTHSW的其他参数设定为: h = 2 μm, θ = 80°, α = 100°和d = 100 nm. 从图2(b)可以看出, 对于同一g, 随着H的增加, Lp呈上升趋势. 当H小于40 μm时, 波导模式以SPPs模式为主, H的增大使得模场能量逐渐远离金属表面, 降低波导内的欧姆损耗, 从而使Lp迅速增加[27]. 当H大于40 μm时, 模式能量开始向Si基楔形波导转移, 模式以SWWs为主, Ag引起的欧姆损耗减小, 使Lp有所增加. 从图2(b)内局部放大图可以看出, 当H小于37 μm时, Lpg的增大而减小, 但变化较小; 当H大于37 μm时, Lpg的增大而增大. 在H等于37 μm时, 不同g下的传播长度均为42.5 × 103 μm. 之后随H的增大Lp逐渐增大; 当H达到90 μm时, Lp从47.7 × 103 μm增大至61.1 × 103 μm.

      图  2  不同Hg时, WWTHSW的模式分析 (a) MC, (b) Lp, (c) A; 模场分布: (d) [H, g] = [10, 0.05] μm, (e) [H, g] = [40, 0.05] μm, (f) [H, g] = [90, 0.05] μm

      Figure 2.  Modes analysis of the WWTHSW with different H and g: (a) MC, (b) Lp, and (c) A;and normalized EM energy density distributions: (d) [H, g] = [10, 0.05] μm; (e) [H, g] = [40, 0.05] μm; (f) [H, g] = [90, 0.05] μm.

      图2(c)可以看出, 对于相同的g, A随着H的增加而先减小再增大, H约40 μm处A达到最小值. 对于较小的SWWs(H小于40 μm), 混合模式主要由SPPs模式决定(MC < 0.5), 观察归一化能量密度分布图2(d)可发现, 该模式约束相对较弱. 当楔形波导的尺寸增大时(H大于40 μm), WWTHSW模式表现为低损耗SWWs模式(MC > 0.5), 大多数模式能量主要集中在高介电常数的SWWs芯部(如图2(f)), 从而导致模式面积增大. 在H = 40 μm时(MC = 0.5), 波导模式介于SWWs和SPPs模式之间, 处于一种临界状态, 同时具有SWWs和SPPs的模式特性, 此时模式面积A有最小值. 当H不变, A随着g的增大而增大, 这是由于g的增大使得模式能量分布逐渐远离Ag表面, 波导模式能量更加分散, A随之增大, 变化范围从3.62 × 10–4到3.69 × 10–3. 因此, 本文以下的研究中, 选择H = 40 μm, g = 50 nm, 此时, Lp为51.5 × 103 μm, A为3.62 × 10–4. 与文献[23]提出的LR-HTSPPs波导相比, Lp增大将近3倍, 而A减小1个数量级.

      SWWs的αθ也是影响WWTHSW模式性能的重要几何参数. 传播长度Lpαθ变化的规律如图3(a)所示, 对于任意θ, Lp随着α的增大而增大. 当α < 70°时, θ的变化对Lp几乎没有影响; 而当α > 70°时, Lp随着θ的增大而减小. Aαθ的变化关系如图3(b)所示, 对于任意θ, A随着α的增大先减小再增加, 在α = 100°附近存在最小值. 如图3(b)的插图所示, 当α = 100°时, A也随着θ的增大先减小再增加. 图3(c)给出在θ = 80°时, 不同α的WWTHSW的归一化能量密度分布图. 当α = 40°和100°时, 波导模式以SPPs模式为主, 电磁场能量主要集中在Si楔形角与Ag楔形角之间, α增大, 模场能量逐渐远离金属表面, 欧姆损耗减小, 从而使Lp增大; 此时, 模式能量更加集中, A反而减小. 当α = 140°时, 模式以SWWs为主, 模式能量主要集中在SWWs内, 导致欧姆损耗减小, Lp增大, 模场能量无法限制在交界面, 从而使A也增大. 图3(d)给出在α = 100°时, 不同θ的WWTHSW的归一化能量密度分布图. 此时, 波导模式以SPPs模式为主, 随着θ从40°增加到80°再到140°, 模场能量与金属楔形的接触面积不断扩大, Ag引入的欧姆损耗增加, Lp减小. 当θ = 80°时, 模场能量更集中于双楔形之间, A达到最小值3.62 × 10–4, 此时传播长度为5.14 × 104 μm. 在本文以下研究中, 选用参数θ = 80°, α = 100°.

      图  3  不同αθ时, WWTHSW的模式分析 (a) Lp, (b) A; (c)模场分布随α的变化(θ = 80°); (d)模场分布随θ的变化(α = 100°)

      Figure 3.  Modes analysis of the WWTHSW with different α and θ, (a) Lp, (b) A; and normalized EM energy density distributions: (c) with different α at a fixed θ of 80°, (d) with different θ at a fixed α of 100°.

      图4(a)(b)给出WWTHSW的LpA随着dh的变化规律. 由图4可知, LpA随着d的增加而略微减小. 当d不变, h = 0时, LpAh = 2 μm和5 μm时大. 结果表明, 增加Ag楔形结构(h = 2 μm和5 μm)后的波导相比于没有增加Ag楔形结构(h = 0)时的Lp略有减小, 但却具有更好的模场限制能力. 从图4(c)中的归一化能量密度分布可以看出, 当h为2 μm和5 μm时, 光场全部被限制在超深亚波长区域内, 从而减小了有效模场面积; 而在h = 2 μm和h = 5 μm时的微米边缘的模式轮廓差几乎可以忽略不计. 对于h = 2 μm, WWTHSW的横向模式宽度W = 0.01 μm(W, 其定义为能量密度衰减到其峰值的1/e的全宽度[17]), 相比于h = 0时横向波导模式的W = 1.2 μm, 减小12倍. 本文以下的研究中, 选用d = 100 nm, h = 2 μm.

      图  4  不同dh时, WWTHSW的模式分析 (a) Lp、(b) A; (c)沿x方向的归一化能量密度

      Figure 4.  Modes analysis of the WWTHSW with different d and h, (a) Lp, (b) A, and (c) normalized EM energy density.

    • 为了比较不同对称性混合表面等离子体波导结构的特性, 利用FEM对WWTHSW、HTMWSPPs波导和HTBTSPPs波导的模式特性和传输特性进行了分析. HTMWSPPs波导、HTBTSPPs波导的截面模型如图5(a)所示.

      图  5  不同波导性能比较 (a) WWTHSW, HTMWSPPs和HTBTSPPs波导的截面图; (b) WWTHSW, HTMWSPPs和HTBTSPPs波导的ALp关系图; (c)品质因数

      Figure 5.  Performance comparison of the WWTHSW, HTMWSPPs and HTBTSPPs wavguide: (a) cross-section views; (b) the relationship between A and Lp; and (c) FOM with different parameters.

      图5(b)LpA的相关性对比图. 在仿真中WWTHSW和HTMWSPPs波导选择最佳参数H = 40 μm, h = 2 μm, θ = 80°, α = 100°, g = 50 nm, d = 100 nm. HTBTSPPs波导选择最佳参数H = 10 μm, W = L = 30 μm. 由图5(b)可知, WWTHSW比HTMWSPPs波导和HTBTSPPs波导的Lp更长. 相比于HTMWSPPs波导, WWTHSW由于尖端场增强效应, 光场主要聚集在楔形波导的顶点附近[1], 具有更强的模场限制能力; 而相比于HTBTSPPs波导, WWTHSW没有矩形Si波导部分, 有利于模场能量的集中, 因此具有更小的模场面积. g在50 nm到2 μm范围内, WWTHSW的有效模场面积(Am = λ2/10280)相比于HTBTSPPs波导(Am = λ2/5405)减小近2倍, 相比于HTMWSPPs波导(Am = λ2/7407)减小近1.5倍. 在相同的有效模场面积Am = λ2/5405的情况下, WWTHSW的Lp为65×103 μm, 是HTMWSPPs波导(Lp = 54 × 103 μm)的1.2倍, 是HTBTSPPs波导(Lp = 33 × 103 μm)的2倍. 而在相同的传播长度Lp = 51 × 103 μm的情况下, WWTHSW的有效模场面积(Am = λ2/10280), 比HTBTSPPs波导(Am = λ2/4422)减小2倍. 由图5(c)可知, 相比于HTBTSPPs波导、HTMWSPPs波导, WWTHSW品质因数更好. 因此, WWTHSW相比HTBTSPPs波导、HTMWSPPs波导在相似的传播长度下, 具有更强的模场限制能力和更好的品质因数.

      为了分析WWTHSW的耦合特性, 构建两根间距为D的平行WWTHSW, 如图6(a)所示. 通常, 波导的耦合特性由波导的耦合长度Lc = π/(βsβa)来衡量, 其中βsβa分别是两个相邻波导的对称和反对称模的传播常数[28].

      图  6  波导耦合特性分析 (a)平行波导三维结构示意图; (b)耦合长度随D的变化; (c)最大传输功率随D的变化

      Figure 6.  Coupling characteristic of waveguides: (a) schematic diagram parallel waveguides; (b) Lc versus the separation between the two waveguides; (c) the maximum transfer power (Pmax) as a function of distance D.

      WWTHSW、HTBTSPPs和HTMWSPPs波导的Lc随两条平行波导中心间距D的对应关系如图6(b)所示. 在仿真中, 采用上述最佳参数, H = 40 μm, h = 2 μm, θ = 80°, α = 100°, g = 50 nm, d = 100 nm. 图6(a)为WWTHSW在D = 12 μm参数下的对称和反对称模式的归一化能量密度分布. 从图(6)中可以看出, WWTHSW的Lc高于HTBTSPPs波导、HTMWSPPs波导的Lc, 这种差异随着间隔D增大变化更加明显. 当D = 12 μm时, HTMWSPPs波导的Lc为2195 μm, 而WWTHSW的Lc为2749 μm, 是HTMWSPPs波导的1.2倍. 当D = 20 μm时, HTMWSPPs波导的Lc为20361 μm, WWTHSW的Lc为28896 μm, 是HTMWSPPs波导的1.5倍. 该结果表明, 在相同的Lc下, 与HTBTSPPs波导和HTMWSPPs波导相比, 本文提出的WWTHSW可显著降低串扰, 这对于高度集成的光子器件来说是一个极大的优势.

      另外, 波导耦合过程中的传播损耗是不可避免的, 因此, 用最大功率传递[29](Pmax)表明波导损耗特性,

      ${P_{\max }} = \frac{{\exp \left( - 2\tau \arctan \dfrac{1}{\tau }\right)}}{{1 + {\tau ^2}}},~~\tau = \dfrac{{2{L_{\rm{c}}}}}{{\text{π}{L_{{\rm{mp}}}}}},$

      其中, Lmp为对称和反对称模式的平均衰减长度,

      ${L_{{\rm{mp}}}} = {2 / {\left[ {{\rm{Re}} \left( {{\beta _{\rm{s}}}} \right) + {\rm{Re}} \left( {{\beta _{\rm{a}}}} \right)} \right]}},$

      其中, Re(βs)和Re(βa)分别是βsβa的实部.

      图6(b)为WWTHSW、HTMWSPPs波导和HTBTSPPs波导的Pmax与间距D的关系图. 从图6(b)可以看出, 对于相同的最大功率传输Pmax, 两个HTWWSPP波导之间的间隔距离D小于两个HTBTSPPs波导和HTMWSPPs波导之间的间隔距离D. 忽略波导间串扰(Pmax < 10–3)[29] 时, WWTHSW最小间隔距离D为15.8 μm, HTBTSPPs波导最小间隔距离D为17.3 μm, HTMWSPPs波导最小间隔距离D为18.6 μm; WWTHSW最小间隔距离相比于HTBTSPPs波导最小间隔距离减小1.5 μm, 相比于HTMWSPPs波导最小间隔距离减小2.8 μm. 结果表明WWTHSW能够实现更高的封装密度.

    • 本文研究了一种可实现低损耗传输的对称双楔形太赫兹混合表面等离子体波导. 利用数值模拟对其支持的表面等离子体模式的传输特性和场分布特性进行研究发现, Ag层的厚度d、双楔形介质波导间隙g及楔形波导几何参数θαH都会对WWTHSW的模式特性产生影响. 在最佳参数H = 40 μm, h = 2 μm, θ = 80°, α = 100°, g = 50 nm, d = 100 nm的条件下, 当Am达到λ2/10280时, Lp接近5.1 × 104 μm. 相比于其他已经提出的混合等离子体波导结构, WWTHSW的归一化模式面积减少一个数量级, Lp增加约3倍. WWTHSW能够在太赫兹频率实现低损耗传播和超深亚波长模式限制. 通过对WWTHSW的耦合特性进行分析可知, 在忽略波导间串扰时, 波导间距小于16 μm时的耦合长度约为8958 μm, 相比于HTMWSPPs波导、HTBTSPPs波导结构, 具有串扰低、封装密度高的优势. 因此, WWTHSW的研究对深度亚波长尺寸下光子和太赫兹器件的高密度集成具有重要意义, 为微米光子器件和太赫兹器件的设计提供一种可行的波导结构.

参考文献 (29)

目录

    /

    返回文章
    返回