光学仪器  2018, Vol. 40 Issue (4): 47-55   PDF    
太赫兹金属光子晶体结构研究
陈麟1,2, 杨冰洋2, 雍有2, 杨涵2, 杨洁2, 许恩豪2     
1. 上海理工大学 上海市现代光学系统重点实验室, 上海 200093;
2. 上海理工大学 光电信息与计算机工程学院, 上海 200093
摘要: 鉴于金属光子晶体结构在滤波、极化、传感、成像等领域的重要作用,研究了基于三角形晶格的金属光子晶体的超常透射、光子能带和表面波特性。采用三维时域有限差分法对金属光子晶体结构建模,同时解释超常透射机理;推导和讨论光子能带沿着ΓK和ΓM晶格方向的能带图,分析不同的表面波特征,并解释非对称入射下观察到的三角形的两个晶格方向模式分裂效应;建立了描述这种模式分裂效应的精确的表面模式谐振频率公式,并对金属光子晶体的色散进行了定量分析。研究结果可为太赫兹功能器件的研究提供参考。
关键词: 时域有限差分     金属光子晶体     光子能带     表面波     角度依赖性    
Study of terahertz metallic photonic crystals
CHEN Lin1,2, YANG Bingyang2, YONG You2, YANG Han2, YANG Jie2, XU Enhao2     
1. Shanghai Key Laboratory of Modern Optical System, University of Shanghai for Science and Technology, Shanghai 200093, China;
2. School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China
Abstract: Owing to its wide applications in the fields of filtering, polarization property, sensing, and imaging, the metallic photonic crystal structure has become the most attractive research field in the past decade.Hence the questions on how to design and analyze the photonic band structure and surface wave modes in the metallic photonic crystals have become important in the development of such frequency selective surfaces devices.Here, based on our previous work, we present a short review to describe the mode characteristics of metallic photonic crystal with a triangular lattice, including extraordinary transmission, photonic bands, and the interaction between surface wave and metallic photonic crystals.Firstly, we used 3-dimension finite difference time domain method with perfectly matched layer to provide a generation process of the surface wave in order to explain the transmittance mechanism of metallic photonic crystal.Then, we derived and discussed photonic band diagrams for two lattice directions and analyzed the characteristics of surface wave along the planar surface of metallic photonic crystal.Splitting of resonance modes had also been observed by varying the incident polar angle along ΓK and ΓM lattice directions, which is due to the off-normal incidence.In addition, exact formulas of angular-dependent resonance frequencies were derived to quantitatively analysis the dispersion of metallic photonic crystal.Such angle dependence results drive the development of related applications and can be used to analyze and design devices based on metallic photonic crystal as well as its attractive application prospects.
Key words: finite difference time domain     metallic photonic crystal     photonic band     surface waves     angular dependence    
引言

在20世纪70年代末, Yeh和Yariv首次报道了周期性平面结构的表面波, 随后带动了微纳光子学领域相关技术的发展。特别是在Ebbesen等发现在金属膜中的亚波长孔阵列能够实现特定波长的超常透射特性后[1], 这个现象引起了广泛关注[2-3]。Pendry等理论证明了表面等离子体的存在[4], 2002年, Thio等报道在牛眼结构中也有类似的现象[5], 不久之后, Degiron等在实验中对单侧和双侧牛眼结构的传输信号进行了比较分析[6]。由于金属孔阵列的光学厚度是由介电衬底涂覆的, 表面等离子体可以在可见光和太赫兹范围内激发, 并在透射谱中出现Fano共振[7-25]。当金属孔阵列的周期与电磁波在谐振频率上的波长可比拟时, 表面波的作用起主导作用[26-28]。与金属孔阵列相关的内容可以在中找到。

对比金属膜厚度只有亚微米量级的金属孔阵列, 其金属光子晶体因有着能够控制电磁波传播的潜能而受到广泛关注。此时金属层的厚度范围远超过金属孔阵列, 达到波长量级, 表明金属光子晶体有着独有的光子能带与表面等离子激元性质。由于金属光子晶体在太赫兹频段易于制造, 大量研究关注于太赫兹频段内表面波引起的金属光子晶体的传播特性上, 例如, 金属光子晶体表面电场的强约束使我们有机会研究具有高灵敏度和大面积成像特性的传感技术。根据透射谱的角度变化可建立表面等离子模式引起的光子能带结构。近年来, 在太赫兹频段内, 正方形阵列的金属光子晶体的超透射现象已有报道。对于三角形阵列, 作者所在的课题组对这种表面波的色散特性在太赫兹域通过角分辨光谱进行了详细研究。本文基于之前的成果讨论了金属光子晶体中ΓK和ΓM方向的完整太赫兹光子能带图及其角度依赖特性。

1 利用三维时域有限差分法分析超透射

三维时域有限差分法是通过离散化的方式精确求解麦克斯韦方程。通过使用Yee元胞, 如图 1所示, 每个部分的精确坐标如表 1所示, 并可表示为

(1)
图 1 Yee单元上的电磁场分量 Figure 1 Electromagnetic field components on Yee unit

表 1 电磁分量的差分格式 Table 1 Difference of electromagnetic components

式中:ijk为不同坐标轴上的Yee单元序列; n为计算的时间步骤; xyzt分别为空间坐标和时间坐标。

利用表 2所示的对偶关系可导出离散的麦克斯韦方程满足的一系列的递归方程组, 其物理意义是表示了不同传播位置处的电场分布。

表 2 递归方程的对偶关系 Table 2 The dual relation of the recursion equation

图 2所示的吸收边界中, 选择基于吸收边界条件算法的完全匹配层。

图 2 区域的划分和样品, 输入和输出面的位置 Figure 2 The division of domain and the location of sample, input and output plane

我们利用三角排列的金属光子晶体来验证上述三维时域有限差分算法的准确性。如图 3所示, 所给出的金属光子晶体由在铝(Al)板上的三角形孔阵列组成, 参数为t=0.25 mm, d=0.7 mm, p=1.13 mm, 太赫兹波垂直入射到金属光子晶体平面上。金属光子晶体几何尺寸为50 mm, 用以提供足够的周期性延伸。

图 3 尺寸为t=0.25 mm, d=0.7 mm和p=1.13 mm的金属光子晶体的示意图, 以及沿着x(ΓM)或y(ΓK)方向斜入射的示意图 Figure 3 Schematics of a metallic PC with the regarding dimension, t=0.25 mm, d=0.7 mm, and p=1.13 mm, and θ the polar angle of the incident radiation along x(ΓM) or y(ΓK) directions

在太赫兹范围内, 铝的介电常数可以用Drude模型来描述, 即

(2)

式中:ω为入射波的角频率; ε(ω)=1;ωp=2.276×1016 rad/s; γ=7.141×1016 rad/s。

此时, 透射率T可以表示为

(3)

式中Win为入射波的总能量。图 4显示了ExHz在垂直x轴的结构对称平面上的场分布。图 4(a)(b)的场分布表明, 电磁波的方向被结构束缚, 在金属和空气的界面上确实存在表面等离子波[29]

图 4 ExHz的场分布 Figure 4 Field distributions of Ex and Hz

图 5所示为入射波的透射率与入射波频率的关系。很容易看出, 在共振频率0.26 THz时, 透光率高达76.52%。

图 5 基于三维时域有限差分方法的透射率曲线 Figure 5 Transmittance based on 3D-FDTD method

为了验证上述理论结果, 我们利用太赫兹时域光谱仪测量了上述样品的透射谱[30-32]。参考光谱和输出信号光谱如图 6所示, 输出谱的谐振频率为0.27 THz。计算结果与实验结果吻合较好。

图 6 使用太赫兹时域光谱仪的实验结果 Figure 6 Experimental results by using THz-TDS
2 金属光子晶体上的ΓM和ΓK方向上太赫兹光子能带图

我们重点研究三角阵列的金属光子晶体的完整光子能带图。结构示意图如图 3所示, 其中ΓM和ΓK方向已在图中标出。考虑到金属光子晶体中的波矢和入射波波矢之间相互作用及匹配条件, 如果入射太赫兹波在金属光子晶体上与光栅常数p、入射角θ等参数密切相关, 则金属光子晶体平面上的表面波矢量ksur可表示如下[24]:

(4)

式中:mn为整数; kxx方向传播的的波矢分量; GxGy为六边形网格中沿xy方向的倒格矢。kx又可定义为

(5)

式中k//为平行于金属光子晶体表面的波矢。GxGy可以表示为:

(6)
(7)

则式(4)可以改写为

(8)

当电磁波从外部入射时, 金属光子晶体可以支持表面波, 这种表面波在金属和空气的界面处传播, 并耦合到自由空间中[28]。在孔直径缩小至零时, 即在占空比为零的晶格近似中, 表面波波矢可以表示为

(9)

式中:ksp为表面等离子体共振矢量; fsp为表面波的谐振频率; ε1为覆盖金属光子晶体的材料介电常数, ε2为金属的介电常数。Al的介电常数在1 THz波段为-44 900+i511 000[33-34], 这比可见光范围大得多[35]。所以式(9)右边根号内的项可近似为1。

根据式(8)和式(9), 我们画出了金属光子晶体在ΓK和ΓM方向的完整的光子能带图, 如图 7所示。为了方便起见, 图 7中的色散图折叠到第一布里渊区里[36]。在ΓM方向, 有一个强的谐振峰色散线和三个弱的色散线, 强色散线对应的是主表面波的(-1, 0)模式, 我们发现(-1, 0)模式有陡峭的负斜率, 与真空中的光色散线几乎相同。在ΓK方向, 有一个强的谐振峰线和两个弱谐振峰线, 强谐振峰线对应于表面波(0, -1)模式。

图 7 金属光子晶体在ΓM和ΓK的光子能带图 Figure 7 Photonic band of metallic PC for ΓM and ΓK directions
3 入射角对透射的影响以及模式分裂效应

本节讨论表面等离子体激元的共振峰随入射角的变化情况。太赫兹波在ΓM和ΓK两个方向上以一定的入射角入射, 如图 3所示。从公式(8)和(9)可以知道, 表面波的共振频率fsp不仅受孔阵的几何构造影响, 还受入射角θ的影响。显而易见, k//可以表示为

(10)

如果我们定义f0是在垂直入射时的零介等离子体共振频率, 并定义归一化共振频率为

(11)

利用式(8)、式(10)和式(11), 我们可以得到表面等离子体的共振频率为

(12)

为了验证仿真结果, 我们在实验中将金属光子晶体样品沿x轴方向和y轴方向旋转, 以代替不同的太赫兹波入射角。如图 8所示, 当入射角不是垂直入射时, 透射峰会分裂成两个峰, 而峰值共振频率会发生红移。图 9显示了由式(12)得到的金属光子晶体在表面波模式下的色散特性以及不同入射角产生的模式分裂效应。当入射波垂直照射在表面时, 太赫兹波的能量与金属光子晶体的本征模耦合。随着入射波的改变, 本征模会分为若干个模式, 每个模式都有自身的色散曲线(如图 9)。本征模的分裂与布里渊区有关, 由金属光子晶体的几何结构决定的。从图 8图 9可以看出, 当入射角改变时, 金属光子晶体表面波模式的色散和与共振频率对应的真空中的光色散线之间的相互作用点发生了移动(如图 9(a)中的点A, B, C, D图 9(b)中的点E, F, G, H)和分离(如图 9(a)中的点AI)。当极角为0°~40°(ΓM方向)且模式为(-1, 0)、(0, 1)和(-1, 1)时, 金属光子晶体的共振频率如图 10所示, 从图 10可以看到, 实验结果与理论计算得到的结果相吻合。

图 8 太赫兹波在金属光子晶体上透射的特性 Figure 8 The THz transmission property from metallic PC

图 9 金属表面等离子波模式的色散, 光在真空中不同的角度 Figure 9 Dispersions of plasma wave pattern on metal surface and the light in vacuum at different angle

图 10 理论与实验结果的共振频率 Figure 10 Theoretical and experimental results of resonance frequency
4 结论

简要讨论了金属光子晶体的异常传输, 光子能带图和非对称入射下的模式分裂效应。该研究将对从事相关现象研究的人员有着借鉴作用。因此, 本文给出的结果不仅太赫兹波段, 同样在滤波、成像和其他学科具有潜在的应用。未来, 金属光子晶体结构在太赫兹近场显微术、高分辨率太赫兹成像和可调谐太赫兹滤光片中的应用仍有许多未知的重要现象有待研究和发现。

参考文献
[1] EBBESEN T W, LEZEC H J, GHAEMI H F, et al. Extraordinary optical transmission through sub-wavelength hole arrays[J]. Nature, 1998, 391(6668): 667–669. DOI:10.1038/35570
[2] SIRTORI C, GMACHL C, CAPASSO F, et al. Long-wavelength(lambda approximately 8-11.5 microm) semiconductor lasers with waveguides based on surface plasmons[J]. Optics Letters, 1998, 23(17): 1366–1368. DOI:10.1364/OL.23.001366
[3] DEGIRON A, LEZEC H J, BARNES W L, et al. Effects of hole depth on enhanced light transmission through subwavelength hole arrays[J]. Applied Physics Letters, 2002, 81(23): 4327. DOI:10.1063/1.1526162
[4] PENDRY J B, MARTÍN-MORENO L, GARCÍA-VIDAL F J. Mimicking surface plasmons with structured surfaces[J]. Science, 2004, 305(5685): 847–848. DOI:10.1126/science.1098999
[5] THIO T, LEZEC H J, EBBESEN T W, et al. Giant optical transmission of sub-wavelength apertures:physics and applications[J]. Nanotechnology, 2002, 13(3): 429–432. DOI:10.1088/0957-4484/13/3/337
[6] DEGIRON A, EBBESEN T W. Analysis of the transmission process through single apertures surrounded by periodic corrugations[J]. Optics Express, 2004, 12(16): 3694–3700. DOI:10.1364/OPEX.12.003694
[7] AZAD A K, ZHAO Y, ZHANG W. Transmission properties of terahertz pulses through an ultrathin subwavelength silicon hole array[J]. Applied Physics Letters, 2005, 86(14): 141102. DOI:10.1063/1.1897842
[8] HAN J G, AZAD A K, GONG M F, et al. Coupling between surface plasmons and nonresonant transmission in subwavelength holes at terahertz frequencies[J]. Applied Physics Letters, 2007, 91(7): 071122. DOI:10.1063/1.2771093
[9] ZHANG W, AZAD A K, HAN J, et al. Direct observation of a transition of a surface plasmon resonance from a photonic crystal effect[J]. Physical Review Letters, 2007, 98(18): 183901. DOI:10.1103/PhysRevLett.98.183901
[10] AZAD A K, CHEN H T, KASARLA S R, et al. Ultrafast optical control of terahertz surface plasmons in subwavelength hole arrays at room temperature[J]. Applied Physics Letters, 2009, 95(1): 011105. DOI:10.1063/1.3168510
[11] QU D X, GRISCHKOWSKY D, ZHANG W L. Terahertz transmission properties of thin, subwavelength metallic hole arrays[J]. Optics Express, 2004, 29(8): 896–898.
[12] SUCKLING J R, SAMBLES J R, LAWRENCE C R. Resonant transmission of microwaves through a hexagonal array of holes in a thin metal layer[J]. New Journal of Physics, 2007, 9(4): 101. DOI:10.1088/1367-2630/9/4/101
[13] MASSON J B, PODZOROV A, GALLOT G. Anomalies in the disappearance of the extraordinary electromagnetic transmission in subwavelength hole arrays[J]. Optics Express, 2008, 16(7): 4719–4730. DOI:10.1364/OE.16.004719
[14] JIANG Y W, TZUANG L D, YE Y H, et al. Effect of Wood's anomalies on the profile of extraordinary transmission spectra through metal periodic arrays of rectangular subwavelength holes with different aspect ratio[J]. Optics Express, 2009, 17(4): 2631–2637. DOI:10.1364/OE.17.002631
[15] STONE E K, HENDRY E. Dispersion of spoof surface plasmons in open-ended metallic hole arrays[J]. Physical Review B, 2011, 84(3): 035418. DOI:10.1103/PhysRevB.84.035418
[16] CHEN L, WEI Y M, ZANG X F, et al. Excitation of dark multipolar plasmonic resonances at terahertz frequencies[J]. Scientific Reports, 2016, 6: 22027. DOI:10.1038/srep22027
[17] CHEN L, XU N N, SINGH L, et al. Defect-induced fano resonances in corrugated plasmonic metamaterials[J]. Advanced Optical Materials, 2017, 5(8): 1600960. DOI:10.1002/adom.201600960
[18] TANAKA M, MIYAMARU F, HANGYO M, et al. Effect of a thin dielectric layer on terahertz transmission characteristics for metal hole arrays[J]. Optics Letters, 2005, 30(10): 1210–1212. DOI:10.1364/OL.30.001210
[19] PRZYBILLA F, DEGIRON A, GENET C, et al. Efficiency and finite size effects in enhanced transmission through subwavelength apertures[J]. Optics Express, 2008, 16(13): 9571–9579. DOI:10.1364/OE.16.009571
[20] LEE J W, PARK T H, NORDLANDER P, et al. Optimum areal coverage for perfect transmission in a periodic metal hole array[J]. Applied Physics Letters, 2010, 97(26): 261112. DOI:10.1063/1.3533658
[21] GHAEMI H F, THIO T, GRUPP D E, et al. Surface plasmons enhance optical transmission through subwavelength holes[J]. Physical Review B, 1998, 58(11): 6779–6782. DOI:10.1103/PhysRevB.58.6779
[22] LI C F, JIN Z M, XUE X, et al. Terahertz coherent control of surface Plasmon polariton propagation in subwavelength metallic hole arrays[J]. Applied Physics Letters, 2012, 100: 19115.
[23] MARY A, RODRIGO S G, MARTÍN-MORENO L, et al. Theory of light transmission through an array of rectangular holes[J]. Physical Review B, 2007, 76(19): 195414. DOI:10.1103/PhysRevB.76.195414
[24] GARCIA N, BAI M. Theory of transmission of light by sub-wavelength cylindrical holes in metallic films[J]. Optics Express, 2006, 14(21): 10028–10042. DOI:10.1364/OE.14.010028
[25] HENDRY E, LOCKYEAR M J, GÓMEZ RIVAS J, et al. Ultrafast optical switching of the THz transmission through metallic subwavelength hole arrays[J]. Physical Review B, 2007, 75(23): 235305. DOI:10.1103/PhysRevB.75.235305
[26] LALANNE P, HUGONIN J P. Interaction between optical nano-objects at metallo-dielectric interfaces[J]. Nature Physics, 2006, 2: 551–556. DOI:10.1038/nphys364
[27] AIGOUY L, LALANNE P, HUGONIN J P, et al. Near-field analysis of surface waves launched at nanoslit apertures[J]. Physical Review Letters, 2007, 98(15): 153902. DOI:10.1103/PhysRevLett.98.153902
[28] LIU H, LALANNE P. Microscopic theory of the extraordinary optical transmission[J]. Nature, 2008, 452(7188): 728–731. DOI:10.1038/nature06762
[29] MIYAMARU F, KAMIJYO M, TAKANO K, et al. Characteristics and generation process of surface waves excited on a perfect conductor surface[J]. Optics Express, 2010, 18(16): 17576–17583. DOI:10.1364/OE.18.017576
[30] CHEN L, XU J M, GAO C M, et al. Manipulating terahertz electromagnetic induced transparency through parallel plate waveguide cavities[J]. Applied Physics Letters, 2013, 103(25): 251105. DOI:10.1063/1.4852115
[31] CHEN L, GAO C M, XU J M, et al. Observation of electromagnetically induced transparency-like transmission in terahertz asymmetric waveguide-cavities systems[J]. Optics Letters, 2013, 38(9): 1379–1381. DOI:10.1364/OL.38.001379
[32] CHEN L, CHENG Z X, XU J M, et al. Controllable multiband terahertz notch filter based on a parallel plate waveguide with a single deep groove[J]. Optics Letters, 2014, 39(15): 4541–4544. DOI:10.1364/OL.39.004541
[33] CHENG Z X, CHEN L, ZANG X F, et al. Ultrathin dual-mode filtering characteristics of terahertz metamaterials with electrically unconnected and connected U-shaped resonators array[J]. Optics Communications, 2015, 342: 20–25. DOI:10.1016/j.optcom.2014.12.042
[34] CHENG Z X, CHEN L, ZANG X F, et al. Tunable Plasmon-induced transparency effect based on self-asymmetric H-shaped resonators meta-atoms[J]. Journal of Optics, 2015, 17(3): 035103. DOI:10.1088/2040-8978/17/3/035103
[35] CHEN L, CAO Z Q, OU F, et al. Observation of large positive and negative lateral shifts of a reflected beam from symmetrical metal-cladding waveguides[J]. Optics Letters, 2007, 32(11): 1432–1434. DOI:10.1364/OL.32.001432
[36] STONE E K, HENDRY E. Dispersion of spoof surface plasmons in open-ended metallic hole arrays[J]. Physical Review B, 2011, 84(3): 035418. DOI:10.1103/PhysRevB.84.035418