搜索

文章查询

x

留言板

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

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

非局域颗粒复合介质的相干完美吸收效应

陈志鹏 於文静 高雷

非局域颗粒复合介质的相干完美吸收效应

陈志鹏, 於文静, 高雷
PDF
HTML
导出引用
导出核心图
  • 研究了两束相干光以相同的入射角从左、右两侧分别入射到Au-SiO2复合介质板时, 在不同的体系参数下该复合材料体系发生相干完美吸收的情形. 运用有效媒质理论推导出了复合介质的有效介电常数以及有效磁导率; 在得到有效电磁参数的基础上进一步推导得到平面波入射复合介质板时的反/透射系数. 通过比较分析非局域和局域情况下颗粒复合介质的相干完美吸收现象, 发现当颗粒尺寸很小时非局域效应的影响会导致复合介质产生相干完美吸收的入射光的频率范围显著变宽. 在进一步的解析计算中, 通过调节复合介质板的厚度、入射光波长、金属颗粒体积分数等参数得到了不同情况下产生的相干完美吸收现象, 并由此分析非局域情形下对于相干完美吸收现象的调控.
      通信作者: 高雷, leigao@suda.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11774252)、江苏省自然科学基金(批准号: BK20161210)、江苏省青蓝工程、“333”工程(批准号: BRA2015353)和江苏省高校优势学科建设工程资助的课题.
    [1]

    Dutta-Gupta S, Martin O J F, Gupta S D, Agarwal G S 2012 Opt. Express 20 001330

    [2]

    Sanjeeb D 2015 Opt. Commun. 356 515

    [3]

    Fu Y Y, Xu Y D, Chen H Y, Cummer S 2017 New J. Phys. 20 013015

    [4]

    Huang S, Xie Z W, Chen W D, Lei J Q, Wang F L, Liu K, Li L 2018 Opt. Express 26 7066

    [5]

    Müllers A, Santra B, Baals C, Jiang J, Benary J, Labouvie R, Zezyulin D A, Konotop V V, Ott H 2018 Sci. Adv. 4 eaat6539

    [6]

    Ruppin R 1973 Phys. Rev. Lett. 31 1434

    [7]

    Fuchs R, Claro F 1987 Phys. Rev. B 35 3722

    [8]

    Rojas R, Claro F, Fuchs R 1988 Phys. Rev. B 37 6799

    [9]

    Chang R, Leung P T 2006 Phys. Rev. B 73 125438

    [10]

    Xie H Y, Chung H Y, Leung P T, Tsai D P 2009 Phys. Rev. B 80 155448

    [11]

    Huang Y, Gao L 2013 Prog. Electromagn. Res. 133 591

    [12]

    Huang Y, Bian X, Ni Y X, Miroshnichenko A E, Gao L 2014 Phys. Rev. A 89 053824

    [13]

    McMahon J M, Gray S K, Schatz G C 2010 Nano Lett. 10 3473

    [14]

    Toscano G, Raza S, Jauho A P, Mortensen N A, Wubs M 2012 Opt. Express 20 4176

    [15]

    Zuloaga J, Prodan E, Nordlander P 2009 Nano Lett. 9 887

    [16]

    Esteban R, Borisov A G, Nordlander P, Aizpurua J 2012 Nat. Commun. 3 825

    [17]

    Dong T Y, Ma X K, Mittra R 2012 Appl. Phys. Lett. 101 233111

    [18]

    Stell L, Zhang P, García-Vidal F J, Rubio A, García-González P 2013 J. Phys. Chem. C 117 8941

    [19]

    Maxwell G J C 1906 Philos. Trans. Roy. Soc. London 205 237

    [20]

    Bruggeman D A G 1935 Ann. Phys. (Leipzig) 24 636

    [21]

    Huang Y, Gao L 2013 J. Phys. Chem. C 117 19203

    [22]

    Dias E J C, Iranzo D A, Gonçalves P A D, Hajati Y, Bludov Y V, Jauho A P, Mortensen N A, Koppens F H L, Peres N M R 2018 Phys. Rev. B 97 245405

    [23]

    Agarwal G S, O'Neil S V 1983 Phys. Rev. B 28 487

    [24]

    Mcmahon J, Gray S, Schatz G 2009 Phys. Rev. Lett. 103 097403

    [25]

    Dasgupta B B, Fuchs R 1981 Phys. Rev. B 24 554

  • 图 1  相干完美吸收的示意图

    Fig. 1.  Schematic diagram of coherent perfect absorption.

    图 2  有效媒质理论模型, 红色为金属颗粒, 蓝色为基底介质, 灰色为有效介质

    Fig. 2.  The model of effective medium. The red part is metal particles, the blue part is base medium, and the grey part is effective medium.

    图 3  (a1) f = 0.1, (b1) f = 0.01, (c1) f = 0.0012时有效介电常数的实部; (a2) f = 0.1, (b2) f = 0.01, (c2) f = 0.0012时有效介电常数的虚部随$\lambda $的变化; 此时d为5 ${\text{μ}}{\rm m}$, a为2 nm

    Fig. 3.  (a1), (b1) and (c1) are the real parts of effective permittivity as function of $\lambda $, for (a1) f = 0.1, (b1) f = 0.01, (c1) f = 0.0012; (a2), (b2), (c2) are the imaginary parts of effective permittivity as function of $\lambda $, for (a2) f = 0.1, (b2) f = 0.01, (c2) f = 0.0012. d = 5 ${\text{μ}}{\rm m}$, a = 2 nm.

    图 4  (a1), (b1), (c1) a = 2, 5, 10 nm时, 局域效应下$\lg|r_1+t_2|^2$$\lambda $f的函数关系; (a2), (b2), (c2)对应情况下考虑非局域效应时的结果; 入射角$\theta$ = 45°

    Fig. 4.  $\lg |r_1 \!+\! t_2|^2$ as functions of $\lambda $ and f with different metallic nanoparticle radius (a) a = 2 nm, (b) a = 5 nm, (c)a = 10 nm: (a1), (b1) and (c1) are within the local description and (a2), (b2) and (c2) are within the nonlocal description. The incident angle is $\theta$=45°.

    图 5  (a) d = 2 ${\text{μ}}{\rm m}$, (b) d = 5 ${\text{μ}}{\rm m}$, (c)、d = 10 ${\text{μ}}{\rm m}$时散射光强对数${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$$\lambda $f的函数关系图, 此时入射角$\theta $为45°

    Fig. 5.  ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$ as functions of $\lambda $ and f with thickness of medium plate (a) d = 2 ${\text{μ}}{\rm m}$, (b) d = 5 ${\text{μ}}{\rm m}$, (c) d = 10 ${\text{μ}}{\rm m}$. The incident angle is $\theta $ = 45°.

    图 6  a = 2 nm, d = 5 ${\text{μm}}$, ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$$\lambda $f的函数关系

    Fig. 6.  Color map of ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$ as functions of $\lambda $ and f for a = 2 nm, d = 5 ${\text{μm}}$.

    图 7  f = 0.0012, $\theta $ = 45°时, (a) $\left| {r_1 } \right|$(蓝色)、$\left| {t_2 } \right|$(红色)与$\lambda $的函数关系, (b) $\left| {\Delta \phi } \right|/{\text{π}}$$\lambda $的函数关系, (c) ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$$\lambda $的函数关系

    Fig. 7.  For f = 0.0012, $\theta $ = 45°, (a) $\left| {r_1 } \right|$ (blue), $\left| {t_2 } \right|$ (red) as function of $\lambda $, (b) $\left| {\Delta \phi } \right|/{\text{π}}$ as function of $\lambda $, (c) ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$as function of $\lambda $.

    图 8  $\lambda $ = 310 nm, $\theta = 45^\circ $时, (a) $\left| {r_1 } \right|$ (蓝色), $\left| {t_2 } \right|$(红色)与f的函数关系; (b) $\left| {\Delta \phi } \right|/{\text{π}}$f的函数关系; (c) ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$f的函数关系

    Fig. 8.  For $\lambda $ = 310 nm, $\theta = 45^\circ $, (a) $\left| {r_1 } \right|$ (blue), $\left| {t_2 } \right|$ (red) as function of f, (b) $\left| {\Delta \phi } \right|/{\text{π}}$ as function of f, (c) ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$ as function of f.

  • [1]

    Dutta-Gupta S, Martin O J F, Gupta S D, Agarwal G S 2012 Opt. Express 20 001330

    [2]

    Sanjeeb D 2015 Opt. Commun. 356 515

    [3]

    Fu Y Y, Xu Y D, Chen H Y, Cummer S 2017 New J. Phys. 20 013015

    [4]

    Huang S, Xie Z W, Chen W D, Lei J Q, Wang F L, Liu K, Li L 2018 Opt. Express 26 7066

    [5]

    Müllers A, Santra B, Baals C, Jiang J, Benary J, Labouvie R, Zezyulin D A, Konotop V V, Ott H 2018 Sci. Adv. 4 eaat6539

    [6]

    Ruppin R 1973 Phys. Rev. Lett. 31 1434

    [7]

    Fuchs R, Claro F 1987 Phys. Rev. B 35 3722

    [8]

    Rojas R, Claro F, Fuchs R 1988 Phys. Rev. B 37 6799

    [9]

    Chang R, Leung P T 2006 Phys. Rev. B 73 125438

    [10]

    Xie H Y, Chung H Y, Leung P T, Tsai D P 2009 Phys. Rev. B 80 155448

    [11]

    Huang Y, Gao L 2013 Prog. Electromagn. Res. 133 591

    [12]

    Huang Y, Bian X, Ni Y X, Miroshnichenko A E, Gao L 2014 Phys. Rev. A 89 053824

    [13]

    McMahon J M, Gray S K, Schatz G C 2010 Nano Lett. 10 3473

    [14]

    Toscano G, Raza S, Jauho A P, Mortensen N A, Wubs M 2012 Opt. Express 20 4176

    [15]

    Zuloaga J, Prodan E, Nordlander P 2009 Nano Lett. 9 887

    [16]

    Esteban R, Borisov A G, Nordlander P, Aizpurua J 2012 Nat. Commun. 3 825

    [17]

    Dong T Y, Ma X K, Mittra R 2012 Appl. Phys. Lett. 101 233111

    [18]

    Stell L, Zhang P, García-Vidal F J, Rubio A, García-González P 2013 J. Phys. Chem. C 117 8941

    [19]

    Maxwell G J C 1906 Philos. Trans. Roy. Soc. London 205 237

    [20]

    Bruggeman D A G 1935 Ann. Phys. (Leipzig) 24 636

    [21]

    Huang Y, Gao L 2013 J. Phys. Chem. C 117 19203

    [22]

    Dias E J C, Iranzo D A, Gonçalves P A D, Hajati Y, Bludov Y V, Jauho A P, Mortensen N A, Koppens F H L, Peres N M R 2018 Phys. Rev. B 97 245405

    [23]

    Agarwal G S, O'Neil S V 1983 Phys. Rev. B 28 487

    [24]

    Mcmahon J, Gray S, Schatz G 2009 Phys. Rev. Lett. 103 097403

    [25]

    Dasgupta B B, Fuchs R 1981 Phys. Rev. B 24 554

  • [1] 张波萍, 焦力实, 张 芸, 李向阳. Au/SiO2纳米复合薄膜的微结构及光吸收特性研究. 物理学报, 2006, 55(4): 2078-2083. doi: 10.7498/aps.55.2078
    [2] 高 雷, 陆兴中. 颗粒复合介质在高温下的光学双稳特性. 物理学报, 2004, 53(12): 4373-4377. doi: 10.7498/aps.53.4373
    [3] 宋彤彤, 罗杰, 赖耘. 赝局域有效介质理论. 物理学报, 2020, (): . doi: 10.7498/aps.69.20200196
    [4] 吴丹丹, 佘卫龙. 线性吸收介质非局域线性电光效应的耦合波理论. 物理学报, 2017, 66(6): 064202. doi: 10.7498/aps.66.064202
    [5] 宋亚舞, 孙 华. 非磁性半导体异常磁电阻效应的有效介质理论. 物理学报, 2008, 57(11): 7178-7184. doi: 10.7498/aps.57.7178
    [6] 包科达. 含椭球包体多相复合介质电导率的有效介质理论. 物理学报, 1992, 41(5): 833-840. doi: 10.7498/aps.41.833
    [7] 黄永平, 赵光普, 肖希, 王藩侯. 部分空间相干光束在非Kolmogorov湍流大气中的有效曲率半径. 物理学报, 2012, 61(14): 144202. doi: 10.7498/aps.61.144202
    [8] 马余强, 李振亚. 二元无规混合系统的有效介质理论. 物理学报, 1990, 39(3): 457-463. doi: 10.7498/aps.39.457
    [9] 孙楠楠, 施展, 丁琪, 许伟伟, 沈洋, 南策文. 基于有效介质理论的物理性能计算模型的软件实现. 物理学报, 2019, 68(15): 157701. doi: 10.7498/aps.68.20182273
    [10] 吴亚敏, 陈国庆. 带壳颗粒复合介质光学双稳的温度效应. 物理学报, 2009, 58(3): 2056-2060. doi: 10.7498/aps.58.2056
    [11] 张永棠. 一种广义三模腔光机械系统的相干完美吸收与透射. 物理学报, 2017, 66(10): 107101. doi: 10.7498/aps.66.107101
    [12] 沈云, 于国萍, 傅继武. 一维反激光器完美相干吸收理论分析. 物理学报, 2012, 61(16): 164204. doi: 10.7498/aps.61.164204
    [13] 李倩倩, 陈小刚, 包曙红, 郭军明, 翟丽丽. 非线性柱形涂层复合介质有效的直流-交流电响应. 物理学报, 2013, 62(5): 057201. doi: 10.7498/aps.62.057201
    [14] 陆兴中, 陈国庆, 吴亚敏. 金属/电介质颗粒复合介质光学双稳的温度效应. 物理学报, 2007, 56(2): 1146-1151. doi: 10.7498/aps.56.1146
    [15] 徐岩, 樊炜, 冀彦君, 宋仁刚, 陈兵, 赵振华, 陈达. 非相对论弱相互作用玻色气体的有效场理论处理. 物理学报, 2014, 63(4): 040501. doi: 10.7498/aps.63.040501
    [16] 余田, 张国华, 孙其诚, 赵雪丹, 马文波. 垂直振动激励下颗粒材料有效质量和耗散功率的研究. 物理学报, 2015, 64(4): 044501. doi: 10.7498/aps.64.044501
    [17] 许聪慧, 张国华, 钱志恒, 赵雪丹. 水平激励下颗粒物质的有效质量及耗散功率的研究. 物理学报, 2016, 65(23): 234501. doi: 10.7498/aps.65.234501
    [18] 赵天恩, 伍瑞新, 杨 帆, 陈 平. 周期性层状铁氧体-电介质复合材料中导模模式的有效负折射率. 物理学报, 2006, 55(1): 179-183. doi: 10.7498/aps.55.179
    [19] 李 卫, 郭永权, 李岫梅, 高汝伟, 冯维存, 王 标, 陈 伟, 韩广兵, 张 鹏, 刘汉强. 纳米复合永磁材料的有效各向异性与矫顽力. 物理学报, 2003, 52(3): 703-707. doi: 10.7498/aps.52.703
    [20] 张珉, 陶瑞宝, 周世勋. 具有自相似结构的非均匀复合媒质质量分布的标度指数. 物理学报, 1988, 37(12): 1987-1992. doi: 10.7498/aps.37.1987
  • 引用本文:
    Citation:
计量
  • 文章访问数:  239
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-11-28
  • 修回日期:  2018-12-30
  • 上网日期:  2019-03-23
  • 刊出日期:  2019-03-01

非局域颗粒复合介质的相干完美吸收效应

  • 1. 苏州大学物理科学与技术学院, 苏州纳米科技协同创新中心, 苏州 215006
  • 2. 苏州大学, 江苏省薄膜重点实验室, 苏州 215006
  • 3. 江苏理工学院数理学院, 常州 213001
  • 通信作者: 高雷, leigao@suda.edu.cn
    基金项目: 国家自然科学基金(批准号: 11774252)、江苏省自然科学基金(批准号: BK20161210)、江苏省青蓝工程、“333”工程(批准号: BRA2015353)和江苏省高校优势学科建设工程资助的课题.

摘要: 研究了两束相干光以相同的入射角从左、右两侧分别入射到Au-SiO2复合介质板时, 在不同的体系参数下该复合材料体系发生相干完美吸收的情形. 运用有效媒质理论推导出了复合介质的有效介电常数以及有效磁导率; 在得到有效电磁参数的基础上进一步推导得到平面波入射复合介质板时的反/透射系数. 通过比较分析非局域和局域情况下颗粒复合介质的相干完美吸收现象, 发现当颗粒尺寸很小时非局域效应的影响会导致复合介质产生相干完美吸收的入射光的频率范围显著变宽. 在进一步的解析计算中, 通过调节复合介质板的厚度、入射光波长、金属颗粒体积分数等参数得到了不同情况下产生的相干完美吸收现象, 并由此分析非局域情形下对于相干完美吸收现象的调控.

English Abstract

    • 近年来, 纳米金属颗粒对电磁辐射的响应一直是许多研究工作的主题. 特别令人感兴趣的是, 含有这些粒子的复合样品在远红外中表现出反常的吸收现象, 并且这种反常的吸收依赖于颗粒的大小、体积分数等因素, 这类反常的吸收现象称之为相干完美吸收(coherent perfect absorption, CPA). 由于在纳米级光学操控和数据处理等方面有很多潜在的应用, 对于相干完美吸收的研究已经越来越热门. 有研究表明, 通过改变这些系统参数, 可以实现可控的相干完美吸收[1]. 现有的CPA的计算是在平面波下近似进行的, 但实验室中的普通激光光源是以高斯光束作为输出光源, 已有研究表明即使是使用高斯光束也可以观察到相干完美吸收现象[2]. 而在共轭超材料构成的二维圆柱结构当中, 同样可以实现相干完美吸收[3]. 最新的研究表明, 在红外波段下可以实现超构表面的多波段相干完美吸收, 而且吸收带宽比单波段吸收提高了3倍[4]. 最近也有科学家研究了非线性波的相干完美吸收现象, 并在玻色-爱因斯坦凝聚体中进行了实验验证[5].

      以往的复合介质相干完美吸收的研究并没有考虑非局域效应的存在[1], 而在小尺度金属颗粒中, 非局域效应的影响是不能忽略的. 在研究非局域方面, Ruppin[6]第一次将Mie全波理论扩充到包含纵波的模式并解决了非局域效应情况下的等离子球体的光学性质问题. 而在准静态近似中, Fuchs等[7,8]提出了一种比较简单的方法并以此来计算考虑介质响应的非局域效应情况下金属纳米颗粒以及核壳颗粒的多极极化度, 并由此来研究这些颗粒的光学响应. Leung等利用准静态下半经典非局域理论, 研究了金属纳米壳与分子相互作用[9]以及非局域等离子体激元增强福斯特能量转移的非局域效应[10]. 在研究单分子附着的双金属核壳纳米颗粒表面增强拉曼散射[11]和纳米颗粒表面等离激元[12]时, 还考虑了金属非局域效应的影响. 时域有限差分法(FDTD)[13]被提出用于计算具有任意形状的纳米结构的非局域效应; 有限元方法(FEM)[14]则同样可以研究任意二维形状的纳米结构的介电非局域效应. 金属传导电子的量子效应引起了金属的非局域效应, 所以研究人员又提出了一种全量子的处理方法并且可以用来很好地研究纳米尺度的结构问题[15-18]. 在大尺寸结构中非局域效应的影响是可忽略的, 但在纳米颗粒中由于颗粒的尺度非常小, 此时非局域效应的影响就变得尤其重要. 当纳米颗粒的大小远小于入射波波长时, 可采用有效媒质理论(EMT)[19,20]来研究其非局域效应, 而且已有研究得到了纳米柱下的等价介电参数[21]. 最近有研究表明, 金属中的非局域效应对石墨烯表面等离激元的光谱有很大影响. 因此, 石墨烯表面等离子体可以成为探测金属纳米结构(包括金属薄膜)中非局域效应的工具[22].

      因此, 在研究纳米颗粒复合介质的相干完美吸收效应过程中, 由于金属量子效应的存在, 金属非局域效应的影响是不能忽略的. 本文运用有效媒质理论研究金属球颗粒复合介质的介电参数, 发展非局域复合体系的相干完美吸收的理论计算, 运用理论解析和数值模拟的方法, 分析非局域效应对相干完美吸收产生的影响并且讨论相干完美吸收产生的条件以及调控.

    • 为了研究非局域颗粒复合介质的相干完美吸收, 我们首先建立了一个模型, 如图1所示. 在x-y平面中厚度为d的复合介质板由基底介质与填充介质构成, 该填充介质是非局域金属颗粒, 基底介质是二氧化硅. 这个复合层由两个相干单色波以相同的入射角$\theta $分别从左右两侧入射. 为了方便起见, 我们分别用下标1(2)标记从左(右)入射的向前(向后)传播的波所产生的反射波和透射波, 分别由图中的蓝色实线(红色虚线)表示.

      图  1  相干完美吸收的示意图

      Figure 1.  Schematic diagram of coherent perfect absorption.

      如果从左侧入射产生的反射波(${r_1}$)与从右侧入射产生的透射波(${t_2}$)的振幅相同, 相位相差为${\text{π}}$, 即$|{r_1}| = |{t_2}|$, $|\Delta \phi | = |{\phi _{r1}} - {\phi _{t2}}| = {\text{π}}$, (${\phi _{r1}}$${\phi _{t2}}$分别是${r_1}$${t_2}$的相位), 此时二者会相消即产生相干完美吸收. 并且由于内在的对称性, 即${r_1} = {r_2}, \; {t_1} = $${t_2}$, 导致在介质两侧的总散射幅度也是相同的. 因此, 在入射介质中相消意味着在出射介质中同样会相消, 这也就导致了相干完美吸收.

      图2的结构当中, 非局域金属纳米颗粒的介电常数由横模介电常数$\varepsilon _{\rm{T}} $和纵模介电常数$\varepsilon _{\rm{L}} $来描述, 分别由下式表示[23,24]:

      图  2  有效媒质理论模型, 红色为金属颗粒, 蓝色为基底介质, 灰色为有效介质

      Figure 2.  The model of effective medium. The red part is metal particles, the blue part is base medium, and the grey part is effective medium.

      $\begin{split} \varepsilon _{\rm{T}}^{}(\omega ) =& {\varepsilon _{\rm{g}}} - \dfrac{{{\omega _{\rm{p}}}^2}}{{{\omega ^2} + {\rm{i}}\omega \gamma }}, \\ \varepsilon _{\rm{L}}^{}(\omega ,k) = &{\varepsilon _{\rm{g}}} - \dfrac{{{\omega _{\rm{p}}}^2}}{{{\omega ^2} + {\rm{i}}\omega \gamma - {\beta ^2}{k^2}}},\end{split}$

      式中${\varepsilon _{\rm{g}}}$是与间带跃迁有关的背景介电常数(如果${\varepsilon _{\rm{g}}} = 1$, 此时该金属将只考虑自由电子气对介电的影响), ${\omega _{\rm{p}}}$是金属的等离子频率, $\gamma $是电子阻尼常数. 非局域参数$\beta $等于$\sqrt {3/5} {v_{\rm{F}}}$, 其中${v_{\rm{F}}}$是金属中电子的费米速度. 纵波的波矢由方程$\varepsilon _{\rm{L}}^{}(\omega ,{k_{\rm{L}}}) = 0$给出, 而横波的波矢则满足常规的色散定律: ${k_{\rm{T}}} =$$ (\omega /c){(\varepsilon _{\rm{T}}^{})^{1/2}}{(\mu _{\rm{T}}^{})^{1/2}}$.

      运用有效媒质理论考虑上述带壳金属纳米球颗粒的核壳模型在有效介质中的总散射为零(即金属球壳对电磁波没有散射), 此时带壳非局域金属纳米球颗粒对电磁波的散射效果与有效介质一致, 两者具有相同的电磁性质. 从光散射的角度来看, 可以用有效介质的介电常数${\varepsilon _{{\rm{eff}}}}$和磁导率${\mu _{{\rm{eff}}}}$来表示核壳纳米球的介电参数, 从而描述整个颗粒复合介质体系的有效介电参数. 接下来研究核壳纳米球在有效介质中的电磁散射问题.

      平面波入射时, 入射电场可用如下公式表示:

      $\begin{split}{E^{\rm{I}}} =\; & {E_0}{{\rm{e}}^{ - {\rm{i}}\omega t}}\sum\limits_{l = 1}^\infty {{{\rm{i}}^l}} \dfrac{{2l + 1}}{{l(l + 1)}}\{ \nabla \times [{{r}} \cdot {{\rm{j}}_l}({k_{{\rm{eff}}}}r) \\ &\times {\rm{P}}_l^{(1)}(\cos \theta )\sin \phi ] - {\rm{i}} \cdot \dfrac{1}{{{k_{{\rm{eff}}}}}}\nabla \times \nabla \\ & \times [{{r}} \cdot {{\rm{j}}_l}({k_{{\rm{eff}}}}r) \cdot {\rm{P}}_l^{(1)}(\cos \theta )\cos \phi ]\} ,\end{split}$

      散射电场表示为

      $\begin{split}{E^{\rm{R}}} =\; & {E_0}{{\rm{e}}^{ - {\rm{i}}\omega t}}\sum\limits_{l = 1}^\infty {{{\rm{i}}^l}} \dfrac{{2l + 1}}{{l(l + 1)}}\{ \nabla \times [{{r}} \cdot {a^{\rm{R}}}{{\rm{h}}_l}({k_{{\rm{eff}}}}r)\\ & \times{\rm{P}}_l^{(1)}(\cos \theta )\sin \phi ] - {\rm{i}} \cdot \dfrac{1}{{{k_{{\rm{eff}}}}}}\nabla \times \nabla \\ & \times [{{r}} \cdot {b^{\rm{R}}}{{\rm{h}}_l}({k_{{\rm{eff}}}}r){\rm{P}}_l^{(1)}(\cos \theta )\cos \phi ]\},\end{split}$

      这里的${{\rm{j}}_l}(x)$${{\rm{h}}_l}(x)$分别是第一类球面贝塞尔函数和球汉克尔函数, $k_{{\rm{eff}}} $是有效介质中横波的波矢, 其中${k_{{\rm{eff}}}} = (\omega /c){({\varepsilon _{{\rm{eff}}}})^{1/2}}{({\mu _{{\rm{eff}}}})^{1/2}}$.

      壳层中的电场表示为

      $\begin{split} {E^{\rm{s}}} =\; & {E_0}{{\rm{e}}^{ - {\rm{i}}\omega t}}\sum\limits_{l = 1}^\infty {{{\rm{i}}^l}} \dfrac{{2l + 1}}{{l(l + 1)}}\{ \nabla \times [{{r}} \cdot [a_{\rm{s}}^{{\rm{TM}}}{{\rm{j}}_l}({k_{\rm{s}}}r) \\ &+ b_{\rm{s}}^{{\rm{TM}}}{{\rm{y}}_l}({k_{\rm{s}}}r)] \cdot {\rm{P}}_l^{(1)}(\cos \theta )\sin \phi ]\\ & - {\rm i} \cdot \dfrac{1}{{{k_{\rm{s}}}}}\nabla \times \nabla \times [{{r}} \cdot [a_{\rm{s}}^{{\rm{TE}}}{{\rm{j}}_l}({k_{\rm{s}}}r) \\ &+ b_{\rm{s}}^{{\rm{TE}}}{{\rm{y}}_l}({k_{\rm{s}}}r)] \cdot {\rm{P}}_l^{(1)}(\cos \theta )\cos \phi ]\} , \end{split}$

      式中${{\rm{y}}_l}(x)$是诺依曼函数; $ a_{\rm s}$, $b_{\rm s} $分别是壳层磁散射系数和电散射系数; ${k_{\rm{s}}} = (\omega /c){({\varepsilon _{\rm{s}}})^{1/2}}{({\mu _{\rm{s}}})^{1/2}}$是壳层中横波的波矢.

      非局域金属纳米球颗粒中的电场既有横向电场, 也有纵向电场, 它们分别表示为:

      $\begin{split}E^{\rm{T}} =\; & {E_0}{{\rm{e}}^{ - {\rm{i}}\omega t}}\sum\limits_{l = 1}^\infty {{{\rm{i}}^l}} \dfrac{{2l + 1}}{{l(l + 1)}}\{ \nabla \times [{{r}} \cdot a_{\rm{c}}^{\rm{T}}{{\rm{j}}_l}({k_{\rm{T}}}r) \\ & \times {\rm{P}}_l^{(1)}(\cos \theta )\sin \phi ] - {\rm{i}} \cdot \dfrac{1}{{{k_{\rm{T}}}}}\nabla \times \nabla \\ & \times [{{r}} \cdot b_{\rm{c}}^{\rm{T}}{{\rm{j}}_l}({k_{\rm{T}}}r) \cdot {\rm{P}}_l^{(1)}(\cos \theta )\cos \phi ]\} ,\end{split}$

      $\begin{split}E^{\rm L} = \; & {E_0}{{\rm{e}}^{ - {\rm{i}}\omega t}}\sum\limits_{l = 1}^\infty {{{\rm{i}}^l}} \dfrac{{2l + 1}}{{l(l + 1)}} \cdot \dfrac{1}{{{k_{\rm{L}}}}}\nabla [a_{\rm{c}}^{\rm{L}}{{\rm{j}}_l}({k_{\rm{L}}}r)\\ & \times {\rm{P}}_l^{(1)}(\cos \theta )\cos \phi ],\end{split}$

      其中$a_{\rm c} $, $b_{\rm {c}} $ 分别为核层磁散射系数和电散射系数.

      结合电场和磁场的边界条件可以计算出各式中的相关系数. 值得注意的是, 由于非局域金属纳米球核内存在附加的纵波, 在金属纳米球核和外壳之间的界面上应该有附加的边界条件. 这里我们选用${{n}} \cdot {{{P}}_{{\rm{ex}}}} = 0$作为附加的边界条件, 其中${{{P}}_{{\rm{ex}}}} = {{D}} - $${\varepsilon _0}{\varepsilon _{\rm{g}}}{{E}}$表示介电材料中极化矢量的极子部分[25].

      经过一系列的推算, 我们得到了散射系数$a_n $$b_n $:

      ${a_n} = \frac{{\left| {\begin{array}{*{20}{c}} { - {{\rm{j}}_l}({k_{{\rm{eff}}}}b)}&{ - {{\rm{j}}_l}({k_{\rm{s}}}b)}&{ - {{\rm{y}}_l}({k_{\rm{s}}}b)}&0 \\ { - \frac{{{{[{k_{{\rm{eff}}}}b \cdot {{\rm{j}}_l}({k_{{\rm{eff}}}}b)]}'}}}{{{\mu _{{\rm{eff}}}}}}}&{ - \frac{{{{[{k_{\rm{s}}}b \cdot {{\rm{j}}_l}({k_{\rm{s}}}b)]}'}}}{{{\mu _{\rm{s}}}}}}&{ - \frac{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_l}({k_{\rm{s}}}b)]}'}}}{{{\mu _{\rm{s}}}}}}&0 \\ 0&{{{\rm{j}}_l}({k_{\rm{s}}}a)}&{{{\rm{y}}_l}({k_{\rm{s}}}a)}&{ - {{\rm{j}}_l}({k_{\rm{T}}}a)} \\ 0&{\frac{{{{[{k_{\rm{s}}}a \cdot {{\rm{j}}_l}({k_{\rm{s}}}a)]}'}}}{{{\mu _s}}}}&{\frac{{{{[{k_{\rm{s}}}a \cdot {{\rm{y}}_l}({k_{\rm{s}}}a)]}'}}}{{{\mu _{\rm{s}}}}}}&{ - \frac{{{{[{k_{\rm{T}}}a \cdot {{\rm{j}}_l}({k_{\rm{T}}}a)]}'}}}{{{\mu _{\rm{T}}}}}} \end{array}} \right|}}{{\left| {\begin{array}{*{20}{c}} {{{\rm{h}}_l}({k_{{\rm{eff}}}}b)}&{ - {{\rm{j}}_l}({k_{\rm{s}}}b)}&{ - {{\rm{y}}_l}({k_{\rm{s}}}b)}&0 \\ {\frac{{{{[{k_{{\rm{eff}}}}b \cdot {{\rm{h}}_l}({k_{{\rm{eff}}}}b)]}'}}}{{{\mu _{{\rm{eff}}}}}}}&{ - \frac{{{{[{k_{\rm{s}}}b \cdot {{\rm{j}}_l}({k_{\rm{s}}}b)]}'}}}{{{\mu _{\rm{s}}}}}}&{ - \frac{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_l}({k_{\rm{s}}}b)]}'}}}{{{\mu _{\rm{s}}}}}}&0 \\ 0&{{{\rm{j}}_l}({k_{\rm{s}}}a)}&{{{\rm{y}}_l}({k_{\rm{s}}}a)}&{ - {{\rm{j}}_l}({k_{\rm{T}}}a)} \\ 0&{\frac{{{{[{k_{\rm{s}}}a \cdot {{\rm{j}}_l}({k_{\rm{s}}}a)]}'}}}{{{\mu _{\rm{s}}}}}}&{\frac{{{{[{k_{\rm{s}}}a \cdot {{\rm{y}}_l}({k_{\rm{s}}}a)]}'}}}{{{\mu _{\rm{s}}}}}}&{ - \frac{{{{[{k_{\rm{T}}}a \cdot {{\rm{j}}_l}({k_{\rm{T}}}a)]}'}}}{{{\mu _{\rm{T}}}}}} \end{array}} \right|}},\tag{7a}$

      $ {b_n} = \frac{{\left| {\begin{array}{*{20}{c}} { - \frac{{{{[{k_{{\rm{eff}}}}b \cdot {{\rm{j}}_l}({k_{{\rm{eff}}}}b)]}'}}}{{{k_{{\rm{eff}}}}b}}}&{ - \frac{{{{[{k_{\rm{s}}}b \cdot {{\rm{j}}_l}({k_{\rm{s}}}b)]}'}}}{{{k_{\rm{s}}}b}}}&{ - \frac{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_l}({k_{\rm{s}}}b)]}'}}}{{{k_{\rm{s}}}b}}}&0&0 \\ { - \frac{{{k_{{\rm{eff}}}}}}{{{\mu _{{\rm{eff}}}}}}{{\rm{j}}_l}({k_{{\rm{eff}}}}b)}&{ - \frac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{j}}_l}({k_{\rm{s}}}b)}&{ - \frac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{y}}_l}({k_{\rm{s}}}b)}&0&0 \\ 0&{\frac{{{{[{k_{\rm{s}}}a \cdot {{\rm{j}}_l}({k_{\rm{s}}}a)]}'}}}{{{k_{\rm{s}}}a}}}&{\frac{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_l}({k_{\rm{s}}}a)]}'}}}{{{k_{\rm{s}}}a}}}&{ - \frac{{{{[{k_{\rm{T}}}a \cdot {{\rm{j}}_l}({k_{\rm{T}}}a)]}'}}}{{{k_{\rm{T}}}a}}}&{\frac{{{{\rm{j}}_l}({k_{\rm{L}}}a)}}{{{k_{\rm{L}}}a}}} \\ 0&{\frac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{j}}_l}({k_{\rm{s}}}a)}&{\frac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{y}}_l}({k_{\rm{s}}}a)}&{ - \frac{{{k_{\rm{T}}}}}{{{\mu _{\rm{T}}}}}{{\rm{j}}_l}({k_{\rm{T}}}a)}&0 \\ 0&0&0&{l(l + 1)({\varepsilon _{\rm{g}}} - {\varepsilon _{\rm{T}}})\frac{{{{\rm{j}}_l}({k_{\rm{T}}}a)}}{{{k_{\rm{T}}}a}}}&{ - {\varepsilon _{\rm{g}}}{{[{{\rm{j}}_l}({k_{\rm{L}}}a)]}'}} \end{array}} \right|}}{{\left| {\begin{array}{*{20}{c}} {\frac{{{{[{k_{{\rm{eff}}}}b \cdot {{\rm{h}}_l}({k_{{\rm{eff}}}}b)]}'}}}{{{k_{{\rm{eff}}}}b}}}&{ - \frac{{{{[{k_{\rm{s}}}b \cdot {{\rm{j}}_l}({k_{\rm{s}}}b)]}'}}}{{{k_{\rm{s}}}b}}}&{ - \frac{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_l}({k_{\rm{s}}}b)]}'}}}{{{k_{\rm{s}}}b}}}&0&0 \\ {\frac{{{k_{{\rm{eff}}}}}}{{{\mu _{{\rm{eff}}}}}}{{\rm{h}}_l}({k_{{\rm{eff}}}}b)}&{ - \frac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{j}}_l}({k_{\rm{s}}}b)}&{ - \frac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{y}}_l}({k_{\rm{s}}}b)}&0&0 \\ 0&{\frac{{{{[{k_{\rm{s}}}a \cdot {{\rm{j}}_l}({k_{\rm{s}}}a)]}'}}}{{{k_{\rm{s}}}a}}}&{\frac{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_l}({k_{\rm{s}}}a)]}'}}}{{{k_{\rm{s}}}a}}}&{ - \frac{{{{[{k_{\rm{T}} }a \cdot {{\rm{j}}_l}({k_{\rm{T}}}a)]}'}}}{{{k_{\rm{T}}}a}}}&{\frac{{{{\rm{j}}_l}({k_{\rm{L}}}a)}}{{{k_{\rm{L}}}a}}} \\ 0&{\frac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{j}}_l}({k_{\rm{s}}}a)}&{\frac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{y}}_l}({k_{\rm{s}}}a)}&{ - \frac{{{k_{\rm{T}}}}}{{{\mu _{\rm{T}}}}}{{\rm{j}}_l}({k_{\rm{T}}}a)}&0 \\ 0&0&0&{l(l + 1)({\varepsilon _{\rm{g}}} - {\varepsilon _{\rm{T}}})\frac{{{{\rm{j}}_l}({k_{\rm{T}}}a)}}{{{k_{\rm{T}}}a}}}&{ - {\varepsilon _{\rm{g}}}{{[{{\rm{j}}_l}({k_{\rm{L}}}a)]}'}} \end{array}} \right|}}.\tag{7b}$

      考虑极限情况${k_{{\rm{eff}}}}b < < 1$时, 带壳非局域金属纳米球颗粒的总散射截面主要由(7)式中级次$n = 1$的项决定. 因此, 有效介电常数${\varepsilon _{{\rm{eff}}}}$和有效磁导率${\mu _{{\rm{eff}}}}$则分别由${a_1} = 0$${b_1} = 0$项决定. 例如, 对于$b_n $有:

      $\left| {\begin{array}{*{20}{c}} { - \dfrac{{{{[{k_{{\rm{eff}}}}b \cdot {{\rm{j}}_l}({k_{{\rm{eff}}}}b)]}'}}}{{{k_{{\rm{eff}}}}b}}}&{ - \dfrac{{{{[{k_{\rm{s}}}b \cdot {{\rm{j}}_l}({k_{\rm{s}}}b)]}'}}}{{{k_{\rm{s}}}b}}}&{ - \dfrac{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_l}({k_{\rm{s}}}b)]}'}}}{{{k_{\rm{s}}}b}}}&0&0 \\ { - \dfrac{{{k_{{\rm{eff}}}}}}{{{\mu _{{\rm{eff}}}}}}{{\rm{j}}_l}({k_{{\rm{eff}}}}b)}&{ - \dfrac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{j}}_l}({k_{\rm{s}}}b)}&{ - \dfrac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{y}}_l}({k_{\rm{s}}}b)}&0&0 \\ 0&{\dfrac{{{{[{k_{\rm{s}}}a \cdot {{\rm{j}}_l}({k_{\rm{s}}}a)]}'}}}{{{k_{\rm{s}}}a}}}&{\dfrac{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_l}({k_{\rm{s}}}a)]}'}}}{{{k_{\rm{s}}}a}}}&{ - \dfrac{{{{[{k_{\rm{T}}}a \cdot {{\rm{j}}_l}({k_{\rm{T}}}a)]}'}}}{{{k_{\rm{T}}}a}}}&{\dfrac{{{{\rm{j}}_l}({k_{\rm{L}}}a)}}{{{k_{\rm{L}}}a}}} \\ 0&{\dfrac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{j}}_l}({k_{\rm{s}}}a)}&{\dfrac{{{k_{\rm{s}}}}}{{{\mu _{\rm{s}}}}}{{\rm{y}}_l}({k_{\rm{s}}}a)}&{ - \dfrac{{{k_{\rm{T}}}}}{{{\mu _{\rm{T}}}}}{{\rm{j}}_l}({k_{\rm{T}}}a)}&0 \\ 0&0&0&{l(l + 1)({\varepsilon _{\rm{g}}} - {\varepsilon _{\rm{T}}})\dfrac{{{{\rm{j}}_l}({k_{\rm{T}}}a)}}{{{k_{\rm{T}}}a}}}&{ - {\varepsilon _{\rm{g}}}{{[{{\rm{j}}_l}({k_{\rm{L}}}a)]}'}} \end{array}} \right| = 0.$

      考虑上述极限的情况下, 令$x = {k_{{\rm{eff}}}}r$(或$x = {k_{\rm{T}}}r$), 函数${{\rm{j}}_l}(x)$, ${{\rm{h}}_l}(x)$, ${{\rm{y}}_l}(x)$可以近似为${{\rm{j}}_1}(x) \cong \dfrac{x}{3}$, ${{\rm{h}}_1}(x) \cong \dfrac{x}{3} - \dfrac{i}{{{x^2}}}$, ${{\rm{y}}_1}(x) \cong - \dfrac{1}{{{x^2}}}$, 则有效介电常数可表示为:

      $\dfrac{{{\varepsilon _{{\rm{eff}}}} - 2{\varepsilon _{\rm{s}}}\dfrac{{{{\rm{j}}_1}({k_{\rm{s}}}b)}}{{{{[{k_{\rm{s}}}b \cdot {{\rm{j}}_1}({k_{\rm{s}}}b)]}'}}}}}{{{\varepsilon _{{\rm{eff}}}} - 2{\varepsilon _{\rm{s}}}\dfrac{{{{\rm{y}}_1}({k_{\rm{s}}}b)}}{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_1}({k_{\rm{s}}}b)]}'}}}}} = \dfrac{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_1}({k_{\rm{s}}}b)]}'}\left\{ {{Q_1} \cdot {{[{k_{\rm{s}}}a \cdot {{\rm{j}}_1}({k_{\rm{s}}}a)]}'}{\varepsilon _{\rm{T}}} - {{\rm{j}}_1}({k_{\rm{s}}}a){\varepsilon _{\rm{s}}}} \right\}}}{{{{[{k_{\rm{s}}}b \cdot {{\rm{j}}_1}({k_{\rm{s}}}b)]}'}\left\{ {{Q_1} \cdot {{[{k_{\rm{s}}}a \cdot {{\rm{y}}_1}({k_{\rm{s}}}a)]}'}{\varepsilon _{\rm{T}}} - {{\rm{y}}_1}({k_{\rm{s}}}a){\varepsilon _{\rm{s}}}} \right\}}},$

      式中

      ${Q_1} = \dfrac{{{\varepsilon _{\rm{g}}}{{\rm{j}}_1}({k_{\rm{T}}}a){\rm{j}}_1'({k_{\rm{L}}}a)}}{{{\varepsilon _{\rm{g}}}{{[{k_{\rm{T}}}a \cdot {{\rm{j}}_1}({k_{\rm{T}}}a)]}'}{\rm{j}}_1'({k_{\rm{L}}}a) - 2({\varepsilon _{\rm{g}}} - {\varepsilon _{\rm{T}}}){{\rm{j}}_1}({k_{\rm{T}}}a){{\rm{j}}_1}({k_{\rm{L}}}a)/({k_{\rm{L}}}a)}}.$

      同样地, 在${a_1} = 0$时可以得到复合介质的有效磁导率:

      $\dfrac{{{\mu _{{\rm{eff}}}} - 2{\mu _{\rm{s}}}\dfrac{{{{\rm{j}}_1}({k_{\rm{s}}}b)}}{{{{[{k_{\rm{s}}}b \cdot {{\rm{j}}_1}({k_{\rm{s}}}b)]}'}}}}}{{{\mu _{{\rm{eff}}}} - 2{\mu _{\rm{s}}}\dfrac{{{{\rm{y}}_1}({k_{\rm{s}}}b)}}{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_1}({k_{\rm{s}}}b)]}'}}}}} = \dfrac{{{{[{k_{\rm{s}}}b \cdot {{\rm{y}}_1}({k_{\rm{s}}}b)]}'}\left\{ {{{[{k_{\rm{s}}}a \cdot {{\rm{j}}_1}({k_s}a)]}'}{{\rm{j}}_1}({k_{\rm{T}}}a){\mu _{\rm{T}}} - {{[{k_{\rm{T}}}a \cdot {{\rm{j}}_1}({k_{\rm{T}}}a)]}'}{{\rm{j}}_1}({k_{\rm{s}}}a){\mu _{\rm{s}}}} \right\}}}{{{{[{k_{\rm{s}}}b \cdot {{\rm{j}}_1}({k_{\rm{s}}}b)]}'}\left\{ {{{[{k_{\rm{s}}}a \cdot {{\rm{y}}_1}({k_{\rm{s}}}a)]}'}{{\rm{j}}_1}({k_{\rm{T}}}a){\mu _{\rm{T}}} - {{[{k_{\rm{T}}}a \cdot {{\rm{j}}_1}({k_{\rm{T}}}a)]}'}{{\rm{y}}_1}({k_{\rm{s}}}a){\mu _{\rm{s}}}} \right\}}}.$

      在有效媒质理论的基础上, 如果不考虑金属纳米颗粒的非局域效应, 此时金属颗粒中不存在纵波, 可以得到局域情况下复合介质的有效电磁参数${\varepsilon _{{\rm{eff}}}}$, ${\mu _{{\rm{eff}}}}$:

      $\dfrac{{{\varepsilon _{{\rm{eff}}}} - {\varepsilon _{\rm{s}}}}}{{{\varepsilon _{{\rm{eff}}}} + 2{\varepsilon _{\rm{s}}}}} = f \cdot \dfrac{{{\varepsilon _{\rm{m}}} - {\varepsilon _{\rm{s}}}}}{{{\varepsilon _{\rm{m}}} + 2{\varepsilon _{\rm{s}}}}}, \;\; \dfrac{{{\mu _{{\rm{eff}}}} - {\mu _{\rm{s}}}}}{{{\mu _{{\rm{eff}}}} + 2{\mu _{\rm{s}}}}} = f \cdot \dfrac{{{\mu _{\rm{m}}} - {\mu _{\rm{s}}}}}{{{\mu _{\rm{m}}} + 2{\mu _{\rm{s}}}}},$

      式中f为非局域金属颗粒的体积分数, $f = {(a/b)^3}$; $\varepsilon _{\rm{m}} $为金属的介电常数, $\varepsilon _{\rm{s}} $为壳层介质的介电常数.

      对于平面波入射, 其波矢分布在x-y平面内, 且在y方向的分量Q是相同的, 在x方向的分量q可以由波动方程${q^2} = {\varepsilon _{{\rm{eff}}}}{\mu _{{\rm{eff}}}}\dfrac{{{\omega ^2}}}{{{c^2}}} - {Q^2}$给出, Q则表示为$Q = \dfrac{\omega }{c}\sqrt {{\varepsilon _1}{\mu _1}} \sin \theta $. 可以推导得出平面波入射复合介质板时的反射系数r和透射系数t:

      $r = \dfrac{{({S_1}^2 - {S_2}^2)({{\rm{e}}^{ - {\rm{i}}qd}} - {{\rm{e}}^{{\rm{i}}qd}})}}{{{{({S_1} + {S_2})}^2}{{\rm{e}}^{ - {\rm{i}}qd}} - {{({S_1} - {S_2})}^2}{{\rm{e}}^{{\rm{i}}qd}}}},\tag{13a}$

      $t = \dfrac{{4{S_1}{S_2}}}{{{{({S_1} + {S_2})}^2}{{\rm{e}}^{ - {\rm{i}}qd}} - {{({S_1} - {S_2})}^2}{{\rm{e}}^{{\rm{i}}qd}}}},\tag{13b}$

      式中$c$是光在真空中的传播速度, $\omega $为入射光的频率, d为复合介质板的厚度; 相关系数${S_1} = \dfrac{{\sqrt {{\varepsilon _1}{\mu _1}} \cos \theta }}{{c{\varepsilon _0}{\varepsilon _1}}}$, ${S_2} = \dfrac{q}{{\omega {\varepsilon _0}{\varepsilon _{{\rm{eff}}}}}}$.

    • 为了研究复合介质的相干完美吸收效应, 我们分别计算了金属颗粒体积分数f为0.1, 0.01, 0.0012时有效介电常数${\varepsilon _{{\rm{eff}}}}$的实部以及虚部随入射光波长$\lambda $的变化情况, 如图3所示.

      图  3  (a1) f = 0.1, (b1) f = 0.01, (c1) f = 0.0012时有效介电常数的实部; (a2) f = 0.1, (b2) f = 0.01, (c2) f = 0.0012时有效介电常数的虚部随$\lambda $的变化; 此时d为5 ${\text{μ}}{\rm m}$, a为2 nm

      Figure 3.  (a1), (b1) and (c1) are the real parts of effective permittivity as function of $\lambda $, for (a1) f = 0.1, (b1) f = 0.01, (c1) f = 0.0012; (a2), (b2), (c2) are the imaginary parts of effective permittivity as function of $\lambda $, for (a2) f = 0.1, (b2) f = 0.01, (c2) f = 0.0012. d = 5 ${\text{μ}}{\rm m}$, a = 2 nm.

      图3中可以发现, 考虑非局域效应对复合介质相干完美吸收的影响, 金属颗粒的体积分数较低时(f = 0.01, 0.0012)${\varepsilon _{{\rm{eff}}}}$的实部会在入射波波长310 nm左右有一个峰值, 当体积分数较高(f为0.1)时这一峰值会出现在320 nm附近, 这表明金属颗粒体积分数会影响非局域效应下复合介质的相干完美吸收, 金属颗粒体积分数越小复合介质的有效介电常数实部的峰值也会减小.

      本文研究了散射光强的对数${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$, 认为当${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2} < - 3$时, 两侧的散射很小, 可以忽略不计, 在同时满足$|{r_1}| = |{t_2}|$, $|\Delta \phi | = |{\phi _{r1}} - {\phi _{t2}}| = $${\text{π}}$时即可产生所谓的相干完美吸收. 图4绘制出了非局域和局域情况下用不同尺寸的金属颗粒填充复合介质时${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$$\lambda $f 的函数关系.

      图  4  (a1), (b1), (c1) a = 2, 5, 10 nm时, 局域效应下$\lg|r_1+t_2|^2$$\lambda $f的函数关系; (a2), (b2), (c2)对应情况下考虑非局域效应时的结果; 入射角$\theta$ = 45°

      Figure 4.  $\lg |r_1 \!+\! t_2|^2$ as functions of $\lambda $ and f with different metallic nanoparticle radius (a) a = 2 nm, (b) a = 5 nm, (c)a = 10 nm: (a1), (b1) and (c1) are within the local description and (a2), (b2) and (c2) are within the nonlocal description. The incident angle is $\theta$=45°.

      图4(a1), (b1), (c1)的对比可以很清楚地看到, 金属颗粒的尺度越小, 各体积分数下产生完美吸收的入射光波长的范围就越宽, 而且波长的极大值也会对应减小, 即频率增加. 从图4(a1)(a2), (b1)(b2)的对比也很容易看出, 金属颗粒的尺寸越小, 在非局域效应的影响下会导致完美吸收产生的入射光的频率显著增加, 当a = 10 nm时, 非局域和局域的结果几乎没有区别, 即金属非局域效应的影响完全可以忽略不计. 这也进一步说明金属非局域效应在小尺度结构中的影响是不能忽略的.

      在研究小尺度颗粒填充情况下金属非局域效应的影响时, 选取a = 2 nm, 讨论非局域效应下不同的系统参数的变化对复合介质的相干完美吸收产生的影响. 首先, 考虑复合介质板的厚度d对相干完美吸收的影响. 图5绘制了不同介质板厚度下, ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$$\lambda $f 的函数关系图像.

      图  5  (a) d = 2 ${\text{μ}}{\rm m}$, (b) d = 5 ${\text{μ}}{\rm m}$, (c)、d = 10 ${\text{μ}}{\rm m}$时散射光强对数${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$$\lambda $f的函数关系图, 此时入射角$\theta $为45°

      Figure 5.  ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$ as functions of $\lambda $ and f with thickness of medium plate (a) d = 2 ${\text{μ}}{\rm m}$, (b) d = 5 ${\text{μ}}{\rm m}$, (c) d = 10 ${\text{μ}}{\rm m}$. The incident angle is $\theta $ = 45°.

      图5的结果表明, 在体积分数相同的情况下, 复合介质板较窄时可以出现完美吸收的入射光的波长的范围也较窄.

      为了得到产生相干完美吸收所需要的f$\lambda $, 我们同样采取了以上的方法, 在相干完美吸收所需条件的基础上, 首先研究散射光强的问题. 图6绘制了${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$$\lambda $f的函数关系图, 此时选取f的范围为0—0.01.

      图  6  a = 2 nm, d = 5 ${\text{μm}}$, ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$$\lambda $f的函数关系

      Figure 6.  Color map of ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$ as functions of $\lambda $ and f for a = 2 nm, d = 5 ${\text{μm}}$.

      考虑到当散射强度对数${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2} < - 3$时可以产生完美吸收, 可以在图中找到一个符合该条件的点, 如图6中圆圈所示. 该点的位置在f为0.0012左右, $\lambda $为310 nm左右. 接下来我们要确定f = 0.0012, $\lambda $ = 310 nm时能否满足产生相干完美吸收所需要的条件.

      首先考虑f = 0.0012时$\left| {{r_1}} \right|$, $\left| {{t_2}} \right|$$\lambda $$\left| {\Delta \phi } \right|/{\text{π}}$$\lambda $以及${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$$\lambda $的函数关系, 如图7所示.

      图  7  f = 0.0012, $\theta $ = 45°时, (a) $\left| {r_1 } \right|$(蓝色)、$\left| {t_2 } \right|$(红色)与$\lambda $的函数关系, (b) $\left| {\Delta \phi } \right|/{\text{π}}$$\lambda $的函数关系, (c) ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$$\lambda $的函数关系

      Figure 7.  For f = 0.0012, $\theta $ = 45°, (a) $\left| {r_1 } \right|$ (blue), $\left| {t_2 } \right|$ (red) as function of $\lambda $, (b) $\left| {\Delta \phi } \right|/{\text{π}}$ as function of $\lambda $, (c) ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$as function of $\lambda $.

      图7(a)中箭头所示为平面波入射复合介质时满足相干完美吸收条件的点, 可以清楚地看到该满足条件的点在$\lambda $为310 nm处, 而此时采用的f值为0.0012.

      接下来考虑$\lambda $为310 nm时$\left| {r_1 } \right|$, $\left| {t_2 } \right|$f, $\left| {\Delta \phi } \right|/{\text{π}}$f以及${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$f的函数关系, 如图8所示.

      图  8  $\lambda $ = 310 nm, $\theta = 45^\circ $时, (a) $\left| {r_1 } \right|$ (蓝色), $\left| {t_2 } \right|$(红色)与f的函数关系; (b) $\left| {\Delta \phi } \right|/{\text{π}}$f的函数关系; (c) ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$f的函数关系

      Figure 8.  For $\lambda $ = 310 nm, $\theta = 45^\circ $, (a) $\left| {r_1 } \right|$ (blue), $\left| {t_2 } \right|$ (red) as function of f, (b) $\left| {\Delta \phi } \right|/{\text{π}}$ as function of f, (c) ${\log _{10}}{\left| {{r_1} + {t_2}} \right|^2}$ as function of f.

      图8(a)箭头所示, 满足相干完美吸收条件时f在0.0012处. 图7图8的结果充分验证了图6中所选的点在f为0.0012, $\lambda $为310 nm, 此时满足相干完美吸收所需要的条件, 而且该波长也正好是图4所示的f为0.0012时复合介质的有效介电常数实部的峰值, 二者结果相吻合.

    • 本文主要研究了非局域金属颗粒填充的复合介质的相干完美吸收效应, 研究了金属非局域效应对复合介质的相干完美产生的影响以及调控. 运用有效媒质理论建立了由非局域金属颗粒组成的复合材料体系的电磁散射理论, 计算得出了非局域金属颗粒复合材料的有效介电常数及有效磁导率.

      通过对局域和非局域下复合介质相干完美吸收的研究和对比, 我们发现由于金属非局域效应的存在, 对于尺寸非常小的金属颗粒, 可以显著提高复合介质相干完美吸收产生的频率, 而小尺寸金属颗粒填充也会导致相干完美吸收在更宽的频率范围内产生. 进一步研究复合介质的厚度、金属颗粒的体积分数等系统参数对复合介质相干完美吸收的影响, 得出了小体积分数下复合介质相干完美吸收时的金属颗粒体积分数以及入射光波长.

      利用非局域金属颗粒的体积分数、非局域金属颗粒的半径以及复合介质板的厚度等系统参数的变化, 结合非局域效应对复合介质相干完美吸收的影响, 实现对非局域颗粒复合介质的相干完美吸收的调控. 本文对今后的纳米复合材料的光学吸收的研究有一定的帮助.

参考文献 (25)

目录

    /

    返回文章
    返回