搜索

文章查询

x

留言板

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

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

涂覆石墨烯的椭圆形电介质纳米线光波导的模式特性分析

程鑫 薛文瑞 卫壮志 董慧莹 李昌勇

涂覆石墨烯的椭圆形电介质纳米线光波导的模式特性分析

程鑫, 薛文瑞, 卫壮志, 董慧莹, 李昌勇
PDF
HTML
导出引用
导出核心图
  • 设计了一种涂覆石墨烯的椭圆形电介质纳米线光波导. 采用分离变量法, 在椭圆柱坐标系中, 借助Mathieu函数, 得到了色散方程. 通过数值求解色散方程, 可以得到模式的有效折射率和场分布, 从而得到模式的传播长度. 研究了工作波长、结构参数以及石墨烯的费米能对模式特性的影响, 并给出了前五个模式的品质因数. 计算表明, 当波长从4.3 ${\text{μ}}{\rm{m}}$增加到8.8 ${\text{μ}}{\rm{m}}$, 这5个模式的有效折射率的实部减小, 基模和一阶模的传播长度增大, 二阶模的传播长度先增大后减小. 当改变纳米线结构参数半长轴和半短轴时, 对基模和一阶模的模式特性影响较小, 对二阶模的模式特性影响较大. 当石墨烯的费米能从0.45 eV增加到0.72 eV时, 有效折射率的实部减小, 传播长度可以达到2 ${\text{μ}}{\rm{m}}$左右. 分离变量法得到的结果与有限元方法得到的结果完全一致. 本文工作可以为基于涂覆石墨烯的电介质纳米线的光波导的设计、制作和应用提供理论基础.
      通信作者: 薛文瑞, wrxue@sxu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 61378039, 61575115)和国家自然科学基金国家基础科学人才培养基金(批准号: J1103210)资助的课题.
    [1]

    Gao Y X, Ren G B, Zhu B F, Liu H Q, Lian Y D, Jian S S 2014 Opt. Express 22 24322

    [2]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Zhang Y, Grigorieva I V, Dubonos S V, Firsov A A 2004 Science 306 666

    [3]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Katsnelson M I, Grigorieva I V, Dubonos S V, Firsov A A 2005 Nature 438 197

    [4]

    Ju L, Geng B S, Horng J, Girit C, Martin M, Hao Z, Bechtel H A, Liang X G, Zettl A, Shen Y R, Wang F 2011 Nature Nanotechnol. 6 630

    [5]

    Wang J C, Song C, Hang J, Hu Z D, Zhang F 2017 Opt. Express 25 23880

    [6]

    Gao Y X, Ren G B, Zhu B F, Liu H Q, Wang J, Jian S S 2014 Opt. Lett. 39 5909

    [7]

    Jablan M, Buljan M, Soljačić M 2009 Phys. Rev. B 80 245435

    [8]

    Christensen J, Manjavacas A, Thongrattanasiri S, Koppens F H, de García Abajo F J 2012 ACS Nano 6 431

    [9]

    Liu P H, Zhang X Z, Ma Z H, Cai W, Wang L, Xu J J 2013 Opt. Express 21 32431

    [10]

    Xing R, Jian S S 2016 IEEE Photonics Technol. Lett. 28 2649

    [11]

    Zhou X T, Zhang T, Chen L, Hong W, Li X 2014 J. Lightwave Technol. 32 3597

    [12]

    Liu Y, Zhong R B, Ding H, Liu S G 2017 Eur. Phys. J. D 71 83

    [13]

    He X Q, Ning T G, Lu S H, Zheng J J, Li J, Li R J, Pei L 2018 Opt. Express 26 10109

    [14]

    Xing R, Jian S S 2017 IEEE Photonics Technol. Lett. 29 1643

    [15]

    Liu J P, Zhai X, Xie F, Wang L L, Xia S X, Liu H J, Luo X, Shang X J 2017 J. Lightwave Technol. 35 1971

    [16]

    Xing R, Jian S S 2017 IEEE Photonics Technol. Lett. 29 967

    [17]

    Zhu B F, Ren G B, Yang Y, Gao Y X 2015 Plasmonics 10 839

    [18]

    彭艳玲, 薛文瑞, 卫壮志, 李昌勇 2018 光学学报 38 0223002

    Peng Y L, Xue W R, Wei Z Z, Li C Y 2018 Acta Opt. Sin. 38 0223002

    [19]

    彭艳玲, 薛文瑞, 卫壮志, 李昌勇 2018 物理学报 67 038102

    Peng Y L, Xue W R, Wei Z Z, Li C Y 2018 Acta Phys. Sin. 67 038102

    [20]

    Xing R, Jian S S 2016 IEEE Photonics Techol. Lett. 28 2779

    [21]

    卫壮志, 薛文瑞, 程鑫, 李昌勇 2018 物理学报 67 108101

    Wei Z Z, Xue W R, Cheng X, Li C Y 2018 Acta Phys. Sin. 67 108101

    [22]

    Hossein M B, Alexander B Y 2015 J. Phys. Condens. Matter 27 185304

    [23]

    陈卫东, 刘丰 2009 电子学报 37 1624

    Chen W D, Liu F 2009 Acta Elec. Sin. 37 1624

    [24]

    Zhao J J, Tang M, Oh K H, Feng Z H, Zhao K, Liao R L, Fu S N, Shum P P, Liu D M 2017 Photon. Res. 8 261

    [25]

    Nejad R M, Tavakoli F, Wang L X, Guan X, Larochelle S, Rusch L A 2018 J. Lightwave Technol. 36 3794

    [26]

    Lee M S, Park B G, Cho I H, Lee J H 2012 IEEE Electron Device Lett. 33 1613

    [27]

    de Abajo F J G 2010 Rev. Mod. Phys. 82 209

    [28]

    翟利, 薛文瑞, 杨荣草, 韩丽萍 2015 光学学报 35 1123002

    Zhai L, Xue W R, Yang R C, Han L P 2015 Acta Opt. Sin. 35 1123002

    [29]

    Nikitin A Y, Guinea F, García-Vidal F J, Martín-Moreno L 2011 Phys.Rev. B 84 195446

    [30]

    Yeh C 1962 J. Appl. Phys. 33 3235

    [31]

    D Erricolo, G Carluccio 2013 ACM Trans. Math. Soft. 40 8

    [32]

    He S L, Zhang X Z, He Y R 2013 Opt. Express 21 30664

    [33]

    Ye L F, Sui K H, Liu Y H, Zhang M, Liu Q H 2018 Opt. Express 26 15935

    [34]

    Ye S, Wang Z X, Sun C R, Dong C B, Wei B Z, Wu B L, Jian S S 2018 Opt. Express 26 23854

    [35]

    Chen M, Sheng P C, Sun W, Cai J J 2016 Opt.Commun. 376 41

  • 图 1  涂覆石墨烯的椭圆形电介质纳米线光波导的横截面示意图

    Fig. 1.  Cross section of the optical waveguide based on graphene-coated elliptical dielectric nanowire.

    图 2  a = 110 nm, b = 80 nm, ${E_{\rm{F}}}= 0{\rm{.5\; eV}}$$\lambda = 7\;{\text{μ}}{\rm{m}}$的情况下, 前五个模式对应的z方向电场分量${E_z}$(a)—(e)与电场强度$\left| E \right|$分布图(f)—(j)

    Fig. 2.  The z direction electric field component ${E_z}$ (a)—(e)and electric field intensity $\left| E \right|$(f)—(j) corresponding to the first 5 modes with a = 110 nm, b = 80 nm, ${E_{\rm{F}}} = 0{\rm{.5\; eV}}$ and $\lambda = 7\;{\text{μ}}{\rm{m}}$.

    图 3  a = 110 nm, b = 80 nm和${E_{\rm{F}}} = 0{\rm{.5\; eV}}$的情况下, 有效折射率实部(a), 传播长度(b)和品质因数(c)与波长的关系

    Fig. 3.  Dependence of the real part of the effective refractive index (a), propagation length (b) and FOM (c) on the wavelength with a = 110 nm, b = 80 nm and${E_{\rm{F}}} = 0{\rm{.5\; eV}}$.

    图 4  a = 110 nm, b = 80 nm和${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$的情况下, 不同波长时Mode 5的电场强度$\left| E \right|$分布图 (a)$\lambda = 5{\rm{.0}}\;{\text{μ}}{\rm{m}}$; (b)$\lambda = {\rm{7}}{\rm{.5}}\;{\text{μ}}{\rm{m}}$; (c)$\lambda = {\rm{8}}{\rm{.5}}\;{\text{μ}}{\rm{m}}$

    Fig. 4.  The distribution of the electric field intensity $\left| E \right|$ of the Mode 5 with different wavelength when a = 110 nm, b = 80 nm and ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$: (a)$\lambda = 5{\rm{.0}}\;{\text{μ}}{\rm{m}}$; (b)$\lambda = {\rm{7}}{\rm{.5}}\;{\text{μ}}{\rm{m}}$; (c)$\lambda = {\rm{8}}{\rm{.5}}\;{\text{μ}}{\rm{m}}$.

    图 5  b = 90 nm, ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$$\lambda = 7\;{\text{μ}}{\rm{m}}$的情况下, 有效折射率实部(a), 传播长度(b)和品质因数(c)与半长轴a的关系

    Fig. 5.  The real part of the effective refractive index (a), propagation length (b) and FOM (c) as a function of semi-major axis when b = 90 nm, ${E_{\rm{F}}}= 0{\rm{.5 \;eV}}$and $\lambda = 7\;{\text{μ}}{\rm{m}}$.

    图 6  b = 80 nm, ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$$\lambda = 7\;{\text{μ}}{\rm{m}}$的情况下, 半长轴长度取不同值时Mode 5的电场强度$\left| E \right|$分布图 (a) a = 100 nm; (b) a = 120 nm; (c) a = 140 nm

    Fig. 6.  The distribution of the electric field intensity $\left| E \right|$ of the Mode 5 with different length of semi-major axis when b = 80 nm, ${E_{\rm{F}}}= 0{\rm{.5 \;eV}}$ and $\lambda = 7\;{\text{μ}}{\rm{m}}$: (a) a = 100 nm; (b) a = 120 nm; (c) a = 140 nm.

    图 7  a = 110 nm, ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$$\lambda = 7\;{\text{μ}}{\rm{m}}$时, 有效折射率实部(a), 传播长度(b)和品质因数(c)与半短轴b的关系

    Fig. 7.  The real part of the effective refractive index (a), propagation length (b) and FOM (c) as a function of semi-minor axis when a = 110 nm, ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$and ${\rm{7}}\;{\text{μ}}{\rm{m}}$.

    图 8  a = 110 nm, ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$$\lambda = 7\;{\text{μ}}{\rm{m}}$的情况下, 半短轴长度取不同值时Mode 5的电场强度$\left| E \right|$分布图 (a) b = 65 nm; (b) b = 85 nm; (c) b = 105 nm.

    Fig. 8.  The distribution of the electric field intensity $\left| E \right|$ of the Mode 5 with different length of semi-minor axis when a = 110 nm, ${E_{\rm F}}= 0{\rm{.5 \;eV}}$ and $\lambda = 7\;{\text{μ}}{\rm{m}}$: (a) b = 65 nm; (b) b = 85 nm; (c) b = 105 nm.

    图 9  a = 110 nm, b = 80 nm和$\lambda = 7\;{\text{μ}}{\rm{m}}$时, 有效折射率实部(a), 传播长度(b)和品质因数(c)与费米能的关系

    Fig. 9.  The real part of the effective refractive index (a), propagation length (b) and FOM (c) as a function of Fermi energy when a = 110 nm, b = 80 nm and $\lambda = 7\;{\text{μ}}{\rm{m}}$.

    图 10  a = 110 nm, b = 80 nm和$\lambda = 7\;{\text{μ}}{\rm{m}}$的情况下, 费米能取不同值时Mode 5的电场强度$\left| E \right|$分布图 (a) ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$; (b) ${E_{\rm F}} = 0{\rm{.63 \;eV}}$; (c) ${E_{\rm F}} = 0{\rm{.72 \;eV}}$

    Fig. 10.  The distribution of the electric field intensity $\left| E \right|$ of the Mode 5 with different values of Fermi energy when a = 110 nm, b = 80 nm, $\lambda = 7\;{\text{μ}}{\rm{m}}$: (a) ${E_{\rm F}} = 0{\rm{.5 \;eV}}$; (b) ${E_{\rm F}} = {\rm{0}}{\rm{.63 \;eV}}$; (c) ${E_{\rm F}} = {\rm{0}}{\rm{.72 \;eV}}$.

    图 11  ${E_{\rm{F}}}= 0{\rm{.5 \;eV}}$b = 100 nm的条件下, 当a = 100和140 nm时, Mode 5 的有效折射率实部(a)和传播长度(b)随波长变化的曲线图

    Fig. 11.  When a = 100 and 140 nm, the real part of the effective refractive index (a) , propagation length (b) as a function of wavelength at ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$ and b = 100 nm.

    图 12  b = 100 nm和$\lambda = 7\;{\text{μ}}{\rm{m}}$的条件下, 当a = 100和140 nm时, Mode 5 的有效折射率实部(a)和传播长度(b)随石墨烯费米能变化的曲线图

    Fig. 12.  When a = 100 and 140 nm, the real part of the effective refractive index (a), propagation length (b) as a function of graphene Fermi energy at b = 100 nm and $\lambda = 7\;{\text{μ}}{\rm{m}}$.

  • [1]

    Gao Y X, Ren G B, Zhu B F, Liu H Q, Lian Y D, Jian S S 2014 Opt. Express 22 24322

    [2]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Zhang Y, Grigorieva I V, Dubonos S V, Firsov A A 2004 Science 306 666

    [3]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Katsnelson M I, Grigorieva I V, Dubonos S V, Firsov A A 2005 Nature 438 197

    [4]

    Ju L, Geng B S, Horng J, Girit C, Martin M, Hao Z, Bechtel H A, Liang X G, Zettl A, Shen Y R, Wang F 2011 Nature Nanotechnol. 6 630

    [5]

    Wang J C, Song C, Hang J, Hu Z D, Zhang F 2017 Opt. Express 25 23880

    [6]

    Gao Y X, Ren G B, Zhu B F, Liu H Q, Wang J, Jian S S 2014 Opt. Lett. 39 5909

    [7]

    Jablan M, Buljan M, Soljačić M 2009 Phys. Rev. B 80 245435

    [8]

    Christensen J, Manjavacas A, Thongrattanasiri S, Koppens F H, de García Abajo F J 2012 ACS Nano 6 431

    [9]

    Liu P H, Zhang X Z, Ma Z H, Cai W, Wang L, Xu J J 2013 Opt. Express 21 32431

    [10]

    Xing R, Jian S S 2016 IEEE Photonics Technol. Lett. 28 2649

    [11]

    Zhou X T, Zhang T, Chen L, Hong W, Li X 2014 J. Lightwave Technol. 32 3597

    [12]

    Liu Y, Zhong R B, Ding H, Liu S G 2017 Eur. Phys. J. D 71 83

    [13]

    He X Q, Ning T G, Lu S H, Zheng J J, Li J, Li R J, Pei L 2018 Opt. Express 26 10109

    [14]

    Xing R, Jian S S 2017 IEEE Photonics Technol. Lett. 29 1643

    [15]

    Liu J P, Zhai X, Xie F, Wang L L, Xia S X, Liu H J, Luo X, Shang X J 2017 J. Lightwave Technol. 35 1971

    [16]

    Xing R, Jian S S 2017 IEEE Photonics Technol. Lett. 29 967

    [17]

    Zhu B F, Ren G B, Yang Y, Gao Y X 2015 Plasmonics 10 839

    [18]

    彭艳玲, 薛文瑞, 卫壮志, 李昌勇 2018 光学学报 38 0223002

    Peng Y L, Xue W R, Wei Z Z, Li C Y 2018 Acta Opt. Sin. 38 0223002

    [19]

    彭艳玲, 薛文瑞, 卫壮志, 李昌勇 2018 物理学报 67 038102

    Peng Y L, Xue W R, Wei Z Z, Li C Y 2018 Acta Phys. Sin. 67 038102

    [20]

    Xing R, Jian S S 2016 IEEE Photonics Techol. Lett. 28 2779

    [21]

    卫壮志, 薛文瑞, 程鑫, 李昌勇 2018 物理学报 67 108101

    Wei Z Z, Xue W R, Cheng X, Li C Y 2018 Acta Phys. Sin. 67 108101

    [22]

    Hossein M B, Alexander B Y 2015 J. Phys. Condens. Matter 27 185304

    [23]

    陈卫东, 刘丰 2009 电子学报 37 1624

    Chen W D, Liu F 2009 Acta Elec. Sin. 37 1624

    [24]

    Zhao J J, Tang M, Oh K H, Feng Z H, Zhao K, Liao R L, Fu S N, Shum P P, Liu D M 2017 Photon. Res. 8 261

    [25]

    Nejad R M, Tavakoli F, Wang L X, Guan X, Larochelle S, Rusch L A 2018 J. Lightwave Technol. 36 3794

    [26]

    Lee M S, Park B G, Cho I H, Lee J H 2012 IEEE Electron Device Lett. 33 1613

    [27]

    de Abajo F J G 2010 Rev. Mod. Phys. 82 209

    [28]

    翟利, 薛文瑞, 杨荣草, 韩丽萍 2015 光学学报 35 1123002

    Zhai L, Xue W R, Yang R C, Han L P 2015 Acta Opt. Sin. 35 1123002

    [29]

    Nikitin A Y, Guinea F, García-Vidal F J, Martín-Moreno L 2011 Phys.Rev. B 84 195446

    [30]

    Yeh C 1962 J. Appl. Phys. 33 3235

    [31]

    D Erricolo, G Carluccio 2013 ACM Trans. Math. Soft. 40 8

    [32]

    He S L, Zhang X Z, He Y R 2013 Opt. Express 21 30664

    [33]

    Ye L F, Sui K H, Liu Y H, Zhang M, Liu Q H 2018 Opt. Express 26 15935

    [34]

    Ye S, Wang Z X, Sun C R, Dong C B, Wei B Z, Wu B L, Jian S S 2018 Opt. Express 26 23854

    [35]

    Chen M, Sheng P C, Sun W, Cai J J 2016 Opt.Commun. 376 41

  • [1] 彭艳玲, 薛文瑞, 卫壮志, 李昌勇. 涂覆石墨烯的非对称并行电介质纳米线波导的模式特性分析. 物理学报, 2018, 67(3): 038102. doi: 10.7498/aps.67.20172016
    [2] 卫壮志, 薛文瑞, 彭艳玲, 程鑫, 李昌勇. 基于涂覆石墨烯的三根电介质纳米线的THz波导的模式特性分析. 物理学报, 2018, 67(10): 108101. doi: 10.7498/aps.67.20180036
    [3] 陈园园, 邹仁华, 宋钢, 张恺, 于丽, 赵玉芳, 肖井华. 纳米银线波导中表面等离极化波激发和辐射的偏振特性研究. 物理学报, 2012, 61(24): 247301. doi: 10.7498/aps.61.247301
    [4] 岳嵩, 李智, 陈建军, 龚旗煌. 基于耦合介质纳米线的深亚波长局域波导. 物理学报, 2011, 60(9): 094214. doi: 10.7498/aps.60.094214
    [5] 李丹, 梁君武, 刘华伟, 张学红, 万强, 张清林, 潘安练. CdS/CdS0.48Se0.52轴向异质结纳米线的非对称光波导及双波长激射. 物理学报, 2017, 66(6): 064204. doi: 10.7498/aps.66.064204
    [6] 李德生, 张鸿庆. 关于分离变量法的一个注记. 物理学报, 2004, 53(6): 1639-1642. doi: 10.7498/aps.53.1639
    [7] 赖云锋. 合成参数对气-液-固法低温生长MgO纳米线的影响. 物理学报, 2010, 59(12): 8814-8819. doi: 10.7498/aps.59.8814
    [8] 陈勇, 李瑞. 纳米尺度硼烯与石墨烯的相互作用. 物理学报, 2019, 68(18): 186801. doi: 10.7498/aps.68.20190692
    [9] 王天会, 李昂, 韩柏. 石墨炔/石墨烯异质结纳米共振隧穿晶体管第一原理研究. 物理学报, 2019, 68(18): 187102. doi: 10.7498/aps.68.20190859
    [10] 谷季唯, 王锦程, 王志军, 李俊杰, 郭灿, 唐赛. 不同衬底条件下石墨烯结构形核过程的晶体相场法研究. 物理学报, 2017, 66(21): 216101. doi: 10.7498/aps.66.216101
    [11] 曾春来, 唐东升, 刘星辉, 海 阔, 羊 亿, 袁华军, 解思深. 化学气相沉积法中SnO2一维纳米结构的控制生长. 物理学报, 2007, 56(11): 6531-6536. doi: 10.7498/aps.56.6531
    [12] 张保磊, 王家序, 肖科, 李俊阳. 石墨烯-纳米探针相互作用有限元准静态计算. 物理学报, 2014, 63(15): 154601. doi: 10.7498/aps.63.154601
    [13] 盛世威, 李康, 孔繁敏, 岳庆炀, 庄华伟, 赵佳. 基于石墨烯纳米带的齿形表面等离激元滤波器的研究. 物理学报, 2015, 64(10): 108402. doi: 10.7498/aps.64.108402
    [14] 顾云风, 吴晓莉, 吴宏章. 三终端非对称夹角石墨烯纳米结的弹道热整流. 物理学报, 2016, 65(24): 248104. doi: 10.7498/aps.65.248104
    [15] 白清顺, 沈荣琦, 何欣, 刘顺, 张飞虎, 郭永博. 纳米微结构表面与石墨烯薄膜的界面黏附特性研究. 物理学报, 2018, 67(3): 030201. doi: 10.7498/aps.67.20172153
    [16] 陈浩, 张晓霞, 王鸿, 姬月华. 基于磁激元效应的石墨烯-金属纳米结构近红外吸收研究. 物理学报, 2018, 67(11): 118101. doi: 10.7498/aps.67.20180196
    [17] 陈令修, 王慧山, 姜程鑫, 陈晨, 王浩敏. 六方氮化硼表面石墨烯纳米带生长与物性研究. 物理学报, 2019, 68(16): 168102. doi: 10.7498/aps.68.20191036
    [18] 冯秋燕, 姚佰承, 周金浩, 夏汉定, 范孟秋, 张黎, 吴宇, 饶云江. 基于飞秒激光抽运的石墨烯包裹微光纤波导结构的级联四波混频研究. 物理学报, 2015, 64(18): 184214. doi: 10.7498/aps.64.184214
    [19] 魏晓林, 陈元平, 王如志, 钟建新. 含孔缺陷石墨烯纳米条带的电学特性研究. 物理学报, 2013, 62(5): 057101. doi: 10.7498/aps.62.057101
    [20] 张慧珍, 李金涛, 吕文刚, 杨海方, 唐成春, 顾长志, 李俊杰. 石墨烯纳米结构的制备及带隙调控研究. 物理学报, 2017, 66(21): 217301. doi: 10.7498/aps.66.217301
  • 引用本文:
    Citation:
计量
  • 文章访问数:  479
  • PDF下载量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-11-26
  • 修回日期:  2018-12-27
  • 上网日期:  2019-03-23
  • 刊出日期:  2019-03-01

涂覆石墨烯的椭圆形电介质纳米线光波导的模式特性分析

  • 1. 山西大学物理电子工程学院, 太原 030006
  • 2. 山西大学激光光谱研究所, 量子光学与光量子器件国家重点实验室, 太原 030006
  • 3. 山西大学, 极端光学协同创新中心, 太原 030006
  • 通信作者: 薛文瑞, wrxue@sxu.edu.cn
    基金项目: 国家自然科学基金(批准号: 61378039, 61575115)和国家自然科学基金国家基础科学人才培养基金(批准号: J1103210)资助的课题.

摘要: 设计了一种涂覆石墨烯的椭圆形电介质纳米线光波导. 采用分离变量法, 在椭圆柱坐标系中, 借助Mathieu函数, 得到了色散方程. 通过数值求解色散方程, 可以得到模式的有效折射率和场分布, 从而得到模式的传播长度. 研究了工作波长、结构参数以及石墨烯的费米能对模式特性的影响, 并给出了前五个模式的品质因数. 计算表明, 当波长从4.3 ${\text{μ}}{\rm{m}}$增加到8.8 ${\text{μ}}{\rm{m}}$, 这5个模式的有效折射率的实部减小, 基模和一阶模的传播长度增大, 二阶模的传播长度先增大后减小. 当改变纳米线结构参数半长轴和半短轴时, 对基模和一阶模的模式特性影响较小, 对二阶模的模式特性影响较大. 当石墨烯的费米能从0.45 eV增加到0.72 eV时, 有效折射率的实部减小, 传播长度可以达到2 ${\text{μ}}{\rm{m}}$左右. 分离变量法得到的结果与有限元方法得到的结果完全一致. 本文工作可以为基于涂覆石墨烯的电介质纳米线的光波导的设计、制作和应用提供理论基础.

English Abstract

    • 表面等离子激元(SPPs)是一种电磁表面波. 利用SPPs可以在亚波长范围内约束和控制光波的传输, 因此SPPs在纳米光子学中有着重要的应用前景[1].

      石墨烯是一种由碳原子紧密堆积而成的蜂巢晶格的二维结构[2-4], 在其表面可以激发SPPs[5]. 石墨烯有三大优点: 超强的模式约束、相对较大的传播长度、可调的化学势和载流子迁移率[6]. 因此, 在中红外到太赫兹波段, 石墨烯是优选的材料. 许多学者对基于石墨烯的各种结构的波导进行了研究, 例如基于石墨烯带的平板波导[7,8]、基于石墨烯的楔形/槽形波导[9,10]和基于石墨烯的混合波导[11-13]等. 与传统的基于贵重金属的表面等离子波导相比, 基于石墨烯的表面等离子波导有两个优点: 欧姆损耗低、对电磁场有较强的亚波长约束能力[7].

      近年来, 为了传输中红外到太赫兹波, 涂覆石墨烯的电介质纳米线波导受到了人们的关注. 这类波导是在单根或多根对称、非对称分布的圆柱形电介质纳米线上涂一层或多层石墨烯构成. 与平板石墨烯波导相比, 这类波导不存在边缘损耗. 研究表明, 在涂覆单层石墨烯的单根圆柱形纳米线波导中, 通过调节工作波长, 可以实现单模工作[6]; 在涂覆单层石墨烯的方形纳米线波导中, 可以实现较小的归一化模式面积[14]; 在涂覆双层石墨烯的单根圆柱形纳米线波导中, 可实现深度亚波长约束[15]; 在涂覆多层石墨烯的单根圆柱形纳米线波导中, 可进一步增强模式的束缚性[16]; 在涂覆单层石墨烯的并行纳米线波导的狭缝区域, 有明显的场增强效应及超强的梯度力[17,18]; 在涂覆单层石墨烯的非对称并行纳米线波导中, 可以灵活调节模式的传输性能[19]. 在涂覆双层石墨烯的纳米并行线间狭缝区, 可以得到更高的场增强[20]. 在涂覆单层石墨烯的三根并行纳米线波导中, 可以灵活调节模式的有效折射率及传播长度[21]. 在低太赫兹频率下, 利用石墨烯单层和纳米结构石墨烯超表面可以使电介质和金属椭圆柱隐身[22].

      本文设计了一种涂覆单层石墨烯的单根椭圆形电介质纳米线光波导, 采用分离变量法, 对这种波导所支持的模式进行研究. 就工作波长、结构参数以及石墨烯的费米能对模式特性的影响进行了分析. 为了验证分离变量法结果的正确性, 与有限元法(FEM)的结果进行了对比. 这种波导有望在中红外到太赫兹波的保偏传输方面得到应用[23-25]. 本文设计的涂覆石墨烯的椭圆形电介质纳米线光波导利用了表面等离子共振效应. 由于表面等离子激元(SPPs)的传输距离在波长量级, 适用于短距离传输, 在存储器[26]、光操控[27]和光子器件的高密度集成[28]等领域有着潜在的应用前景.

    • 涂覆石墨烯的椭圆形电介质纳米线光波导的横截面如图1所示. 它是在相对介电常数为${\varepsilon _{\rm{d}}}$的椭圆形电介质纳米线的表面涂覆单层石墨烯材料构成, 背景介质是空气, 介电常数为${\varepsilon _{{\rm{air}}}}$. 假设${\varepsilon _{\rm{d}}} = 2.1025$, ${\varepsilon _{{\rm{air}}}} = 1.0$. 波导结构的半长轴为$a$, 半短轴为$b$, 椭圆的左、右两个焦点分别为${f_1}$${f_2}$, 半焦距为$q$. 石墨烯的表面电导率${\sigma _{\rm{g}}}$与工作角频率、环境温度和石墨烯的费米能有关, 由带内电导率${\sigma _{{\rm{intra}}}}$和带间电导率${\sigma _{{\rm{inter}}}}$两部分构成(${\sigma _{\rm{g}}} = {\sigma _{{\rm{intra}}}} + {\sigma _{{\rm{inter}}}}$), 并可以通过Kubo公式得到[29]:

      ${\sigma _{{\rm{intra}}}} = \frac{{2{\rm{j}}{e^2}{k_{\rm{B}}}T}}{{{\hbar ^2}{\text{π}}\left( {\omega + {\rm{j}}\varGamma } \right)}}{\rm{ln}}\left[ {2{\rm{cosh}}\left( {\frac{{{E_{\rm{F}}}}}{{2{k_{\rm{B}}}T}}} \right)} \right], $

      $\begin{split}{\sigma _{{\rm{inter}}}} = & \frac{{{e^2}}}{{4\hbar }}\bigg[ \frac{1}{2} + \frac{1}{{\text{π}}}{\rm{arctan}}\left( {\frac{{\hbar \omega - 2{E_{\rm{F}}}}}{{2{k_{\rm{B}}}T}}} \right) \\ &- \frac{{\rm{j}}}{{2{\text{π}}}}{\rm{ln}}\frac{{{{\left( {\hbar \omega + 2{E_{\rm{F}}}} \right)}^2}}}{{{{\left( {\hbar \omega - 2{E_{\rm{F}}}} \right)}^2} + {{\left( {2{k_{\rm{B}}}T} \right)}^2}}} \bigg], \end{split}$

      其中$e$为电子的电荷量, ${k_{\rm{B}}}$是玻尔兹曼常量, $T = 300\;{\rm{ K}}$是环境温度, $\hbar $是约化后的普朗克常量, $\omega = 2{\text{π}}f$是工作角频率, $\varGamma = 2 \times {10^{12}}\;{\rm{ rad/s}}$是载流子散射率, ${E_{\rm{F}}}$为石墨烯的费米能.

      以椭圆的中心为原点, 建立如图1所示的直角坐标系$o {\text{-}} xy$和椭圆坐标系$o{\text{-}} \xi \eta $. 两套坐标系之间的关系为$x = q\cosh \xi \cos \eta $, $y = q\sinh \xi \sin \eta $, 其中$\xi $为径向坐标, $\eta $为角向坐标, 取值范围分别为$0 \leqslant \xi \leqslant \infty $, $0 \leqslant \eta \leqslant 2{\text{π}}$. 椭圆形电介质纳米线的表面对应的径向坐标为${\xi _0}$. 在椭圆坐标系中, 借助Mathieu函数, 采用分离变量法求解波动方程. 假设波导支持的模式沿着z方向传播, 则场分量均可表示为$F\left( {\xi ,\eta } \right)\exp [\;{\rm{j}}\left( {\beta z - \omega t} \right)]$, 其中$F\left( {\xi ,\eta } \right)$代表电场或者磁场在横截面上的分布, $\beta $为传播常数. 把模式分为奇模(${H_z}$奇对称, ${E_z}$偶对称)与偶模(${H_z}$偶对称, ${E_z}$奇对称)两类[30].

      图  1  涂覆石墨烯的椭圆形电介质纳米线光波导的横截面示意图

      Figure 1.  Cross section of the optical waveguide based on graphene-coated elliptical dielectric nanowire.

      以奇模为例, 假设在柱内($\xi < {\xi _0}$)的纵向场分量表示为:

      $\left\{ {\begin{split} & {\mathop H\nolimits_z^{\left( {\rm{1}} \right)} = \sum\limits_{m = 1}^\infty {{A_m}{\rm{R}}_{{\rm{o}}m}^{\left( {\rm{1}} \right)}} \left( {\xi , - {\gamma _1}^2} \right){{\rm{S}}_{{\rm{o}}m}}\left( {\eta , - {\gamma _1}^2} \right) = \sum\limits_{m = 1}^\infty {{A_m}{\rm{R}}_{{\rm{o}}m}^{\left( {\rm{1}} \right)}} {{\rm{S}}_{{\rm{o}}m}},} \\ & {\mathop E\nolimits_z^{\left( {\rm{1}} \right)} = \sum\limits_{m = 0}^\infty {{B_m}} {\rm{R}}_{{\rm{e}}m}^{\left( {\rm{1}} \right)}\left( {\xi , - {\gamma _1}^2} \right){{\rm{S}}_{{\rm{e}}m}}\left( {\eta , - {\gamma _1}^2} \right) = \sum\limits_{m = 0}^\infty {{B_m}{\rm{R}}_{{\rm{e}}m}^{\left( {\rm{1}} \right)}} {{\rm{S}}_{{\rm{e}}m}},} \end{split}} \right.$

      在柱外($\xi > {\xi _0}$)的纵向场分量表示为:

      $\left\{ {\begin{split} & {\mathop H\nolimits_z^{\left( 0 \right)} = \sum\limits_{m = 1}^\infty {{L_m}{\rm{R}}_{{\rm{o}}m}^{\left( 3 \right)}} \left( {\xi , - {\gamma _0}^2} \right){{\rm{S}}_{{\rm{o}}m}}\left( {\eta , - {\gamma _0}^2} \right) = \sum\limits_{m = 1}^\infty {{L_m}{\rm{R}}_{{\rm{o}}m}^{\left( 3 \right)}} {{\rm{S}}_{{\rm{o}}m}},} \\ & {\mathop E\nolimits_z^{\left( 0 \right)} = \sum\limits_{m = 0}^\infty {{P_m}{\rm{R}}_{{\rm{e}}m}^{\left( 3 \right)}} \left( {\xi , - {\gamma _0}^2} \right){{\rm{S}}_{{\rm{e}}m}}\left( {\eta , - {\gamma _0}^2} \right) = \sum\limits_{m = 0}^\infty {{P_m}{\rm{R}}_{{\rm{e}}m}^{\left( 3 \right)}} {{\rm{S}}_{{\rm{e}}m}},} \end{split}} \right.$

      其中${\rm{R}}_{{\rm{o}}m}^{\left( 1 \right)}$为修正的m阶第一类径向Mathieu奇函数, ${\rm{R}}_{{\rm{e}}m}^{\left( 1 \right)}$为修正的m阶第一类径向Mathieu偶函数, ${\rm R}_{{\rm{o}}m}^{\left( 3 \right)}$为修正的m阶第三类径向Mathieu奇函数, ${\rm R}_{{\rm{e}}m}^{\left( 3 \right)}$为修正的m阶第三类径向Mathieu偶函数, ${{\rm S}_{{\rm{o}}m}}$m阶角向Mathieu奇函数, ${{\rm S}_{{\rm{e}}m}}$m阶角向Mathieu偶函数[31]. 其中${\gamma _1}^2 = \left( {{\beta ^2} - k_1^2} \right){q^2}/4$, $\gamma _0^2 \!=\! \left( {{\beta ^2} \!-\! k_0^2} \right){q^2}/4$, $k_1^2 \!=\! {\omega ^2}{\varepsilon _{{\rm{r}}1}}{\varepsilon _0}{\mu _0}$, $k_0^2 \!=\! {\omega ^2}{\varepsilon _{{\rm{r}}0}}{\varepsilon _0}{\mu _0}$, ${\varepsilon _0}$为真空介电常数, ${\mu _0}$为真空磁导率, ${A_m}$, ${B_m}$, ${L_m}$${P_m}$为待定系数. 为了书写方便, 文中公式均采用Mathieu函数的简写形式.

      由Maxwell方程组可得, 柱内($\xi < {\xi _0}$)的横向场分量为:

      $ \left\{ {\begin{split} & {\mathop E\nolimits_\xi ^{\left( 1 \right)} = - \frac{1}{{p_1^2d}}\left\{ {\sum\limits_{m = 0}^\infty {\left[ { - {\rm{j}}\beta {B_m}{\rm R}_{{\rm{e}}m}^{'\left( 1 \right)}{{\rm S}_{{\rm{e}}m}}} \right]} + \sum\limits_{m = 1}^\infty {\left[ {{\rm{j}}\omega \mu {A_m}{\rm R}_{{\rm{o}}m}^{\left( 1 \right)}{\rm S}_{{\rm{o}}m}'} \right]} } \right\},} \\ & {\mathop E\nolimits_\eta ^{\left( 1 \right)} = - \frac{1}{{p_1^2d}}\left\{ {\sum\limits_{m = 0}^\infty {\left[ {{\rm{j}}\beta {B_m}{\rm R}_{{\rm{e}}m}^{\left( 1 \right)}{\rm S}_{{\rm{e}}m}'} \right]} + \sum\limits_{m = 1}^\infty {\left[ {{\rm{j}}\omega \mu {A_m}{\rm R}_{{\rm{o}}m}^{'\left( 1 \right)}{{\rm S}_{{\rm{o}}m}}} \right]} } \right\},} \\ & {\mathop H\nolimits_\xi ^{\left( 1 \right)} = - \frac{1}{{p_1^2d}}\left\{ {\sum\limits_{m = 0}^\infty {\left[ { - {\rm{j}}\omega {\varepsilon _1}\beta {B_m}{\rm R}_{{\rm{e}}m}^{\left( 1 \right)}{\rm S}_{{\rm{e}}m}'} \right]} + \sum\limits_{m = 1}^\infty {\left[ {{\rm{-j}}\beta {A_m}{\rm R}_{{\rm{o}}m}^{'\left( 1 \right)}{{\rm S}_{{\rm{o}}m}}} \right]} } \right\},} \\ & {\mathop H\nolimits_\eta ^{\left( 1 \right)} = - \frac{1}{{p_1^2d}}\left\{ {\sum\limits_{m = 0}^\infty {\left[ { - {\rm{j}}\omega {\varepsilon _1}\beta {B_m}{\rm R}_{{\rm{e}}m}^{'\left( 1 \right)}{{\rm S}_{{\rm{e}}m}}} \right]} + \sum\limits_{m = 1}^\infty {\left[ {{\rm{j}}\beta {A_m}{\rm R}_{{\rm{o}}m}^{\left( 1 \right)}{\rm S}_{{\rm{o}}m}'} \right]} } \right\},} \end{split}} \right. $

      柱外($\xi > {\xi _0}$)的横向场分量为:

      $ \left\{ {\begin{split} & {\mathop E\nolimits_\xi ^{\left( 0 \right)} = - \frac{1}{{p_0^2d}}\left\{ {\sum\limits_{m = 0}^\infty {\left[ { - {\rm{j}}\beta {P_m}{\rm R}_{{\rm{e}}m}^{'\left( 3 \right)}{{\rm S}_{{\rm{e}}m}}} \right]} + \sum\limits_{m = 1}^\infty {\left[ {{\rm{j}}\omega \mu {L_m}{\rm R}_{{\rm{o}}m}^{\left( 3 \right)}{\rm S}_{{\rm{o}}m}'} \right]} } \right\}}, \\ & {\mathop E\nolimits_\eta ^{\left( 0 \right)} = - \frac{1}{{p_0^2d}}\left\{ {\sum\limits_{m = 0}^\infty {\left[ {{\rm{j}}\beta {P_m}{\rm R}_{{\rm{e}}m}^{\left( 3 \right)}{\rm S}_{{\rm{e}}m}'} \right]} + \sum\limits_{m = 1}^\infty {\left[ {{\rm{j}}\omega \mu {L_m}{\rm R}_{{\rm{o}}m}^{'\left( 3 \right)}{{\rm S}_{{\rm{o}}m}}} \right]} } \right\}}, \\ & {\mathop H\nolimits_\xi ^{\left( 0 \right)} = - \frac{1}{{p_0^2d}}\left\{ {\sum\limits_{m = 0}^\infty {\left[ { - {\rm{j}}\omega {\varepsilon _0}\beta {P_m}{\rm R}_{{\rm{e}}m}^{\left( 3 \right)}{\rm S}_{{\rm{e}}m}'} \right]} + \sum\limits_{m = 1}^\infty {\left[ {{\rm{-j}}\beta {L_m}{\rm R}_{{\rm{o}}m}^{'\left( 3 \right)}{{\rm S}_{{\rm{o}}m}}} \right]} } \right\}}, \\ & {\mathop H\nolimits_\eta ^{\left( 0 \right)} = - \frac{1}{{p_0^2d}}\left\{ {\sum\limits_{m = 0}^\infty {\left[ { - {\rm{j}}\omega {\varepsilon _0}\beta {P_m}{\rm R}_{{\rm{e}}m}^{'\left( 3 \right)}{{\rm S}_{{\rm{e}}m}}} \right]} + \sum\limits_{m = 1}^\infty {\left[ {{\rm{j}}\beta {P_m}{\rm R}_{{\rm{o}}m}^{\left( 3 \right)}{\rm S}_{{\rm{o}}m}'} \right]} } \right\}}, \end{split}} \right. $

      其中${p_1}^2 = {\beta ^2} - {k_1}^2$, ${p_0}^2 = \beta ^2 - k_0^2$, $d = \big[\cosh ^2 \xi - $ $\cos ^2\eta \big]^{1/2} $.

      把石墨烯看成具有表面电导率${\sigma _g}$的无厚度的边界, 场的切向分量之间的边值关系为:

      $\left\{ {\begin{aligned} & {E_z^{\left( 1 \right)}{|_{\xi = {\xi _0}}} = E_z^{\left( 0 \right)}{|_{\xi = {\xi _0}}},} \\[-1pt] & {E_\eta ^{\left( 1 \right)}{|_{\xi = {\xi _0}}} = E_\eta ^{\left( 0 \right)}{|_{\xi = {\xi _0}}},} \\[-1pt] & {H_z^{\left( 1 \right)}{|_{\xi = {\xi _0}}} - H_z^{\left( 0 \right)}{|_{\xi = {\xi _0}}} = {\sigma _g}E_\eta ^{\left( 1 \right)}{|_{\xi = {\xi _0}}},} \\[-1pt] & {H_\eta ^{\left( 1 \right)}{|_{\xi = {\xi _0}}} - H_\eta ^{\left( 0 \right)}{|_{\xi = {\xi _0}}} = - {\sigma _g}E_z^{\left( 1 \right)}{|_{\xi = {\xi _0}}},} \end{aligned}} \right.$

      由(7)式可以得到一个奇次的线性代数方程组:

      $\left| {\begin{array}{*{20}{c}} {\left({a_m} + \dfrac{{{\rm{j}}\omega \mu {\sigma _{\rm{g}}}}}{{p_1^2d}}a_m'\right){\alpha _{mn}}}&{\dfrac{{{\rm{j}}\beta {\sigma _{\rm{g}}}}}{{p_1^2d}}{b_m}{\gamma _{mn}}}&{ - {l_m}{\alpha _{mn}}}&0 \\ 0&{{b_m}{\beta _{mn}}}&0&{ - {p_m}{\beta _{mn}}} \\ { - \dfrac{{{\rm{j}}\beta }}{{p_1^2d}}{a_m}{\chi _{mn}}}&{\left({\sigma _{\rm{g}}}{b_m} + \dfrac{{{\rm{j}}\omega {\varepsilon _1}}}{{p_1^2d}}b_m'\right){\beta _{mn}}}&{\dfrac{{{\rm{j}}\beta }}{{p_0^2d}}{l_m}{\chi _{mn}}}&{ - \dfrac{{{\rm{j}}\beta {\varepsilon _0}}}{{p_0^2d}}p_m'{\beta _{mn}}} \\ {\dfrac{{\omega \mu }}{\beta }a_m'{\alpha _{mn}}}&{{b_m}{\gamma _{mn}}}&{ - \dfrac{{p_1^2}}{{p_0^2}}\dfrac{{\omega \mu }}{\beta }l_m'{\alpha _{mn}}}&{ - \dfrac{{p_1^2}}{{p_0^2}}{p_m}{\gamma _{mn}}} \end{array}} \right| = 0, $

      其中${a_m} = {\rm R}_{{\rm{o}}m}^{\left( 1 \right)}\left( {{\xi _0}} \right)$, $a_m' = {\rm R}_{{\rm{o}}m}^{'\left( 1 \right)}\left( {{\xi _0}} \right)$, ${b_m} = {\rm R}_{{\rm{e}}m}^{\left( 1 \right)}\left( {{\xi _0}} \right)$, $b_m' = {\rm R}_{{\rm{e}}m}^{'\left( 1 \right)}\left( {{\xi _0}} \right)$, ${l_m} = {\rm R}_{{\rm{o}}m}^{\left( 3 \right)}\left( {{\xi _0}} \right)$, $l_m' = {\rm R}_{{\rm{o}}m}^{'\left( 3 \right)}\left( {{\xi _0}} \right)$, ${p_m} =$$ {\rm R}_{{\rm{e}}m}^{\left( 3 \right)}\left( {{\xi _0}} \right)$, $p_m' = {\rm R}_{{\rm{e}}m}^{'\left( 3 \right)}\left( {{\xi _0}} \right)$, ${\alpha _{mn}} = {{\rm S}_{{\rm{o}}m}}\left( {{\eta _n}} \right)$, ${\beta _{mn}} = $${{\rm S}_{{\rm{e}}m}}\left( {{\eta _n}} \right)$, ${\chi _{mn}} = {\rm S}_{{\rm{o}}m}'\left( {{\eta _n}} \right)$, ${\gamma _{mn}} = {\rm S}_{{\rm{e}}m}'\left( {{\eta _n}} \right)$. 通过求解(8)式, 可以得到模式的场分布和传播常数$\beta $, 从而得到模式的有效折射率${n_{{\rm{eff}}}}$、传播长度${L_{{\rm{prop}}}}$和品质因数FOM (figures of merit). 其中, 传播长度定义为${L_{{\rm{prop}}}} = {\lambda /[{4{\text{π}}{\rm{Im}} \left( {{n_{{\rm{eff}}}}} \right)]}}$[13,16], 品质因数定义为${\rm{FOM}} = {{{\rm{Re}} \left( {{n_{{\rm{eff}}}}} \right)}/{{\rm{Im}} \left( {{n_{{\rm{eff}}}}} \right)}}$[32-34].

    • 图2(a)图2(e)分别给出在a = 110 nm, b = 80 nm, ${E_{\rm{F}}} = 0.{\rm{5\; eV}}$$\lambda\; {\rm{ = 7}}\;{\text{μ}}{\rm{m}}$的情况下, 该波导支持的五个最低阶模式(依次命名为Mode 1, Mode 2, Mode 3, ${\rm{Mode 4}}$和Mode 5)对应的z方向电场分量${E_z}$分布, 图2(f)图2(j)为其对应的电场强度$\left| E \right|$分布. 可以看出, 场集中分布在石墨稀涂层附近. 其中Mode 1为基模, Mode 2和Mode 3为一阶模, Mode 4和Mode 5为二阶模. Mode 1, Mode 3和Mode 5为偶模, Mode 2和Mode 4为奇模.

      图  2  在a = 110 nm, b = 80 nm, ${E_{\rm{F}}}= 0{\rm{.5\; eV}}$$\lambda = 7\;{\text{μ}}{\rm{m}}$的情况下, 前五个模式对应的z方向电场分量${E_z}$(a)—(e)与电场强度$\left| E \right|$分布图(f)—(j)

      Figure 2.  The z direction electric field component ${E_z}$ (a)—(e)and electric field intensity $\left| E \right|$(f)—(j) corresponding to the first 5 modes with a = 110 nm, b = 80 nm, ${E_{\rm{F}}} = 0{\rm{.5\; eV}}$ and $\lambda = 7\;{\text{μ}}{\rm{m}}$.

      图3(a)图3(c)分别给出了在a = 110 nm, b = 80 nm和${E_{\rm{F}}} {\rm{ = 0}}{\rm{.5 eV}}$的情况下, 有效折射率实部${\rm{Re}} \left( {{n_{{\rm{eff}}}}} \right)$、传播长度${L_{{\rm{prop}}}}$和品质因数FOM与波长$\lambda $之间的依赖关系. 图中实线为 FEM 的模拟结果, 点线为分离变量法得到的计算结果, 下文中均用该方法标注. 从图3中可以看出, 随着波长的增大, 前五个模式的有效折射率的实部都单调减小, 其中Mode 1, Mode 2和Mode 3 的有效折射率的实部下降较慢, 而Mode 4和Mode 5的有效折射率的实部下降较快. 随着波长的增大, Mode 1, Mode 2与Mode 3的传播长度单调增加, 而Mode 4和Mode 5的传播长度先增大后减小. 随着波长的增大, 品质因数FOM先增大再逐渐减小, 其中Mode 1, Mode 2和Mode 3的品质因数下降较慢, 而Mode 4和Mode 5的品质因数下降较快. 因为品质因数FOM越大, 波导的传输性能越好[35], 在所选参数范围内, Mode 2的品质因数最大, 所以Mode 2的传输性能最好.

      图  3  在a = 110 nm, b = 80 nm和${E_{\rm{F}}} = 0{\rm{.5\; eV}}$的情况下, 有效折射率实部(a), 传播长度(b)和品质因数(c)与波长的关系

      Figure 3.  Dependence of the real part of the effective refractive index (a), propagation length (b) and FOM (c) on the wavelength with a = 110 nm, b = 80 nm and${E_{\rm{F}}} = 0{\rm{.5\; eV}}$.

      上述现象可以通过场分布来解释. 以Mode 5为例, 图4给出了$\lambda\; {\rm{ = 5}}{\rm{.0 }}$, 7.5和8.5 ${\text{μ}}{\rm{m}}$三种情况下电场强度分布的对比图. 从图中可以看出, 当$\lambda\; {\rm{ = 5}}{\rm{.0}}\;{\text{μ}}{\rm{m}}$时, 场被紧密地约束在石墨烯涂层附近, 场与石墨烯涂层之间的相互作用强, 此时波导对场的约束性非常强, 能量分布集中, 但模式的传输损耗较大, 传播长度较小. 当$\lambda\; {\rm{ = 7}}{\rm{.5}}\;{\text{μ}}{\rm{m}}$时, 场分布扩散至电介质纳米线和空气中, 场与石墨烯涂层之间的相互作用变弱, 波导对场的束缚性变弱, 传输损耗减小, 因而传播长度增大. 当$\lambda\; {\rm{ = 8}}{\rm{.5}}\;{\text{μ}}{\rm{m}}$时, 虽然场分布扩散到电介质纳米线和空气中的越来越多, 波导对场的束缚性变弱, 但是石墨烯涂层附近的场的强度增强, 场与石墨烯之间的相互作用变强, 传输损耗增大, 因而传播长度减小.

      图  4  在a = 110 nm, b = 80 nm和${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$的情况下, 不同波长时Mode 5的电场强度$\left| E \right|$分布图 (a)$\lambda = 5{\rm{.0}}\;{\text{μ}}{\rm{m}}$; (b)$\lambda = {\rm{7}}{\rm{.5}}\;{\text{μ}}{\rm{m}}$; (c)$\lambda = {\rm{8}}{\rm{.5}}\;{\text{μ}}{\rm{m}}$

      Figure 4.  The distribution of the electric field intensity $\left| E \right|$ of the Mode 5 with different wavelength when a = 110 nm, b = 80 nm and ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$: (a)$\lambda = 5{\rm{.0}}\;{\text{μ}}{\rm{m}}$; (b)$\lambda = {\rm{7}}{\rm{.5}}\;{\text{μ}}{\rm{m}}$; (c)$\lambda = {\rm{8}}{\rm{.5}}\;{\text{μ}}{\rm{m}}$.

      图5(a)图5(c)分别给出了在b = 90 nm,${E_{\rm{F}}}\;{\rm{ = 0}}{\rm{.5 \;eV}}$$\lambda\; {\rm{ = 7}}\;{\text{μ}}{\rm{m}}$的情况下, 有效折射率实部${\rm{Re}} \left( {{n_{{\rm{eff}}}}} \right)$、传播长度${L_{{\rm{prop}}}}$和品质因数FOM与半长轴a的关系. 从图中可以看出, 当半短轴b的长度一定时, 随着a的增大, 前5个模式的有效折射率的实部逐渐增大, 其中Mode 3的有效折射率的实部变化极小, Mode 1和 Mode 2的有效折射率的实部增长速度比Mode 4和Mode 5的慢. 随着a的增大, Mode 1的传播长度基本不变, Mode 3的传播长度缓慢减小, Mode 2, Mode 4和Mode 5的传播长度缓慢增大. 随半长轴a的增大, Mode 1和Mode 2的品质因数增加较慢, Mode 3的品质因数缓慢减小, Mode 4和Mode 5的品质因数增加较快. 在所选参数范围内, Mode 2的品质因数最大, 传输性能最好.

      图  5  b = 90 nm, ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$$\lambda = 7\;{\text{μ}}{\rm{m}}$的情况下, 有效折射率实部(a), 传播长度(b)和品质因数(c)与半长轴a的关系

      Figure 5.  The real part of the effective refractive index (a), propagation length (b) and FOM (c) as a function of semi-major axis when b = 90 nm, ${E_{\rm{F}}}= 0{\rm{.5 \;eV}}$and $\lambda = 7\;{\text{μ}}{\rm{m}}$.

      上述现象可以通过场分布来解释. 以Mode 5为例, 图6给出了当b = 80 nm, ${E_{\rm{F}}}\;{\rm{ = 0}}{\rm{.5 \;eV}}$, $\lambda \;{\rm{ = 7}}\;{\text{μ}}{\rm{m}}$时, a = 100, 120和140 nm情况下的电场强度$\left| E \right|$分布对比图. 从图中可以看出, 当a = 100 nm时, 石墨烯椭圆环的周长较小, 场被约束在石墨烯椭圆环附近, 场与石墨烯之间的相互作用强, 模式的传输损耗大, 传播长度小. 当a = 120 nm时, 石墨烯椭圆环的周长变大, 场分布的范围扩大, 在半长轴增大的过程中, 左右花瓣对应的场的强度基本保持不变, 而上下花瓣对应的场的强度减小, 场与石墨烯之间的相互作用变弱, 传输损耗减小, 因而传播长度增大. 当a = 140 nm时, 石墨烯椭圆环的周长进一步变大, 场分布进一步扩散, 在半长轴增大的过程中, 左右花瓣对应的场的强度基本保持不变, 而上下花瓣对应的场的强度进一步减小, 场与石墨烯之间的相互作用进一步变弱, 传输损耗进一步减小, 因而传播长度进一步增大.

      图  6  在b = 80 nm, ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$$\lambda = 7\;{\text{μ}}{\rm{m}}$的情况下, 半长轴长度取不同值时Mode 5的电场强度$\left| E \right|$分布图 (a) a = 100 nm; (b) a = 120 nm; (c) a = 140 nm

      Figure 6.  The distribution of the electric field intensity $\left| E \right|$ of the Mode 5 with different length of semi-major axis when b = 80 nm, ${E_{\rm{F}}}= 0{\rm{.5 \;eV}}$ and $\lambda = 7\;{\text{μ}}{\rm{m}}$: (a) a = 100 nm; (b) a = 120 nm; (c) a = 140 nm.

      图7(a)图7(c)分别给出了在a = 110 nm,${E_{\rm{F}}}\;{\rm{ = 0}}{\rm{.5 \;eV}}$$\lambda \;{\rm{ = 7}}\;{\text{μ}}{\rm{m}}$的情况下, 有效折射率实部$\rm{Re} \left( {{n_{{\rm{eff}}}}} \right)$、传播长度${L_{{\rm{prop}}}}$和品质因数FOM与半短轴$b$之间的关系. 从图7中可以看出, 当半长轴a的长度一定时, 随着半短轴b的增大, Mode 1, Mode 3, Mode 4和Mode 5的有效折射率的实部逐渐增大, 其中Mode 5的增加速度最快, Mode 2的有效折射率的实部缓慢减小. 随着半短轴b的增大, 除Mode 2外其余4个模式的传播长度都在增大, 其中Mode 5的增长速度最快. 随着半短轴b的增大, 除Mode 2外其余4个模式的品质因数FOM都单调增大, 其中Mode 1和Mode 3的品质因数增加较慢, Mode 5的品质因数增加较快. 在所选参数范围内, 当b < 95 nm时, Mode 2的品质因数最大, 传输性能最好, 当$b$> 95 nm时, Mode 1的品质因数最大, 传输性能最好.

      图  7  当a = 110 nm, ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$$\lambda = 7\;{\text{μ}}{\rm{m}}$时, 有效折射率实部(a), 传播长度(b)和品质因数(c)与半短轴b的关系

      Figure 7.  The real part of the effective refractive index (a), propagation length (b) and FOM (c) as a function of semi-minor axis when a = 110 nm, ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$and ${\rm{7}}\;{\text{μ}}{\rm{m}}$.

      上述现象可以通过场分布来解释. 以Mode 5为例, 图8给出了b = 65, 85和105 nm情况下的电场强度$\left| E \right|$分布对比图. 从图中可以看出, 当b = 65 nm时, 石墨烯椭圆环的周长较小, 场被约束在石墨烯椭圆环附近, 场与石墨烯之间的相互作用较强, 模式的传输损耗大, 因此传播长度小. 当b = 85 nm时, 石墨烯椭圆环的周长变大, 场分布逐渐扩散, 左右花瓣对应的场的强度有所减弱, 而上下花瓣对应的场的强度基本保持不变, 场与石墨烯之间的相互作用变弱, 传输损耗减小, 因而传播长度增大. 当半短轴增大至105 nm时, 石墨烯椭圆环的周长进一步增大, 场分布的扩散范围也进一步变大, 左右花瓣对应的场的强度进一步减弱, 而上下花瓣对应的场的强度基本保持不变, 场与石墨烯之间的相互作用进一步减弱, 传输损耗进一步减小, 因而传播长度也进一步增大.

      图  8  在a = 110 nm, ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$$\lambda = 7\;{\text{μ}}{\rm{m}}$的情况下, 半短轴长度取不同值时Mode 5的电场强度$\left| E \right|$分布图 (a) b = 65 nm; (b) b = 85 nm; (c) b = 105 nm.

      Figure 8.  The distribution of the electric field intensity $\left| E \right|$ of the Mode 5 with different length of semi-minor axis when a = 110 nm, ${E_{\rm F}}= 0{\rm{.5 \;eV}}$ and $\lambda = 7\;{\text{μ}}{\rm{m}}$: (a) b = 65 nm; (b) b = 85 nm; (c) b = 105 nm.

      图9(a)图9(c)分别给出了在a = 110 nm, b = 80 nm和$\lambda = 7\;{\text{μ}}{\rm{m}}$情况下, 有效折射率实部${\rm Re} ( n_{\rm eff})$、传播长度${L_{{\rm{prop}}}}$和品质因数FOM与费米能EF之间的依赖关系. 从图9中可以看出, 随着费米能的增大, 前5个模式的有效折射率的实部都单调减小, 其中Mode 1, Mode 2和Mode 3的有效折射率的实部下降较慢, 而Mode 4和Mode 5的有效折射率的实部下降较快. 随着费米能的增大, Mode 1, Mode 2与Mode 3的传播长度单调增加, Mode 4和Mode 5的传播长度先增大随后减小. 随着费米能的增大, Mode 1, Mode 2和Mode 3的品质因数缓慢增大, Mode 4和Mode 5的品质因数逐渐减小.

      图  9  当a = 110 nm, b = 80 nm和$\lambda = 7\;{\text{μ}}{\rm{m}}$时, 有效折射率实部(a), 传播长度(b)和品质因数(c)与费米能的关系

      Figure 9.  The real part of the effective refractive index (a), propagation length (b) and FOM (c) as a function of Fermi energy when a = 110 nm, b = 80 nm and $\lambda = 7\;{\text{μ}}{\rm{m}}$.

      上述现象可以通过场分布来解释. 石墨烯的表面电导率可以通过改变费米能级的大小来改变, 因此波导模式的传输性能还可以通过石墨烯的费米能级来调节. 以Mode 5为例, 图10给出了${E_{\rm{F}}}\;{\rm{ = 0}}{\rm{.5}}$, ${\rm{0}}{\rm{.63}}$${\rm{0}}{\rm{.72 \;eV}}$情况下的电场强度$\left| E \right|$分布对比图. 从图中可以看出, 当${E_{\rm{F}}}\;{\rm{ = 0}}{\rm{.5 \;eV}}$时, 场被紧密地约束在石墨烯涂层附近, 场与石墨烯涂层之间的相互作用强, 此时波导对场的约束性非常强, 能量分布集中, 模式的传输损耗较大, 传播长度较小. 当${E_{\rm{F}}}\;{\rm{ = 0}}{\rm{.63 \;eV}}$时, 场的分布有所扩散, 场与石墨烯涂层之间的相互作用变弱, 波导对场的束缚性变弱, 传输损耗减小, 因而传播长度增大. 当${E_{\rm{F}}}\;{\rm{ = 0}}{\rm{.72 \;eV}}$时, 场的强度增强, 场与石墨烯涂层之间的相互作用变强, 传输损耗增大, 因而传播长度减小.

      图  10  在a = 110 nm, b = 80 nm和$\lambda = 7\;{\text{μ}}{\rm{m}}$的情况下, 费米能取不同值时Mode 5的电场强度$\left| E \right|$分布图 (a) ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$; (b) ${E_{\rm F}} = 0{\rm{.63 \;eV}}$; (c) ${E_{\rm F}} = 0{\rm{.72 \;eV}}$

      Figure 10.  The distribution of the electric field intensity $\left| E \right|$ of the Mode 5 with different values of Fermi energy when a = 110 nm, b = 80 nm, $\lambda = 7\;{\text{μ}}{\rm{m}}$: (a) ${E_{\rm F}} = 0{\rm{.5 \;eV}}$; (b) ${E_{\rm F}} = {\rm{0}}{\rm{.63 \;eV}}$; (c) ${E_{\rm F}} = {\rm{0}}{\rm{.72 \;eV}}$.

      以Mode 5为例, 在${E_{\rm{F}}}= 0.5 {\rm \;eV}$b = 100 nm的条件下, 图11给出了当椭圆形纳米线半长轴a = 100和140 nm时, 有效折射率的实部与传播长度随波长变化的曲线图. 图中黑色线代表半长轴a = 100 nm的情形, 红色线代表半长轴a = 140 nm的情形. 由图11可知, 当椭圆形纳米线的半长轴由a = 100 nm (圆形纳米线)变化为a = 140 nm时, 有效折射率的实部与传播长度都增大.

      图  11  在${E_{\rm{F}}}= 0{\rm{.5 \;eV}}$b = 100 nm的条件下, 当a = 100和140 nm时, Mode 5 的有效折射率实部(a)和传播长度(b)随波长变化的曲线图

      Figure 11.  When a = 100 and 140 nm, the real part of the effective refractive index (a) , propagation length (b) as a function of wavelength at ${E_{\rm{F}}} = 0{\rm{.5 \;eV}}$ and b = 100 nm.

      以Mode 5为例, 在$\lambda \;{\rm{ = 7}}\;{\text{μ}}{\rm{m}}$b = 100 nm的条件下, 图12给出了当椭圆形纳米线半长轴为a = 100 nm (圆形纳米线)和a = 140 nm时有效折射率的实部与传播长度随石墨烯费米能变化的曲线图, 图中黑线代表半长轴为${\rm{100 \;nm}}$的情形, 红线代表半长轴为140 nm的情形. 从图12可知, 当椭圆形纳米线的半长轴由a = 100 nm (圆形纳米线)变化为a = 140 nm时, 有效折射率的实部与传播长度都增大. 以上两个例子说明, 椭圆形纳米线相对于圆形纳米线具有一定的优势.

      图  12  在b = 100 nm和$\lambda = 7\;{\text{μ}}{\rm{m}}$的条件下, 当a = 100和140 nm时, Mode 5 的有效折射率实部(a)和传播长度(b)随石墨烯费米能变化的曲线图

      Figure 12.  When a = 100 and 140 nm, the real part of the effective refractive index (a), propagation length (b) as a function of graphene Fermi energy at b = 100 nm and $\lambda = 7\;{\text{μ}}{\rm{m}}$.

    • 本文设计了一种涂覆石墨烯的椭圆形电介质纳米线波导. 在椭圆柱坐标系中, 借助Mathieu函数, 通过分离变量法, 推导出色散方程, 并通过数值方法求解了这个色散方程. 研究了波导最低阶的5个模式的传输特性, 以及工作波长、半长轴、半短轴、费米能的变化对这5个模式的传输特性的影响, 并对各参数变化时, 各个模式的品质因数进行了对比, 选出了最优的模式. 结果表明: 工作波长从4.3 ${\text{μ}}{\rm{m}}$增大至8.8 ${\text{μ}}{\rm{m}}$时, 5个模式的有效折射率的实部都减小, 基模与一阶模的传播长度增大, 二阶模的传播长度先增大后减小; 当半短轴一定, 半长轴从95 nm增大至140 nm时, 有效折射率的实部逐渐增大, 其中Mode 3变化极小, 除Mode 3的传播长度缓慢减小外其他模式的传播长度都增大; 当半长轴一定, 半短轴从60 nm增大至105 nm时, 除Mode 2的有效折射率的实部与传播长度减小外, 其他模式的有效折射率的实部与传播长度都增大; 费米能从${\rm{0}}{\rm{.45\;eV}}$${\rm{0}}{\rm{.72 \;eV}}$变化时, 前5个模式的有效折射率的实部都减小, 基模与一阶模的传播长度增大, 二阶模的传播长度先增大后减小. 分离变量法的计算结果与有限元法仿真的结果一致. 本文所研究的涂覆石墨烯的椭圆形电介质波导在中红外到太赫兹波的保偏传输方面具有应用前景.

参考文献 (35)

目录

    /

    返回文章
    返回