Analysis of Filter and Waveguide Effect Based on the MIM Nanodisk with a Metallic Block

Abstract

A novel and meaningful plasmonic filter based on the metal-insulator-metal (MIM) waveguides directly connected to a nanocavity with a metallic block is proposed and demonstrated numerically. By the effective index method and the resonant theory of disk-shaped nanocavity, we reveal that the resonant wavelengths can be simply tuned by adjusting the height of the block, which is in accordance with the results calculated by finite-difference time-domain (FDTD) simulations. We also rotate the metallic block to achieve two resonant modes. One mode shows a blue shift, and the other mode shows a red shift. It is shown that the proposed structure performs as a bend waveguide not a filter when the width of the block increases to hundreds of nanometers. The proposed structure will have significant potential applications in nanointegrated circuits for optical filtering and processing.

Academic Editor:Christian Brosseau

Key Lab for Micro-Nano Physics and Technology, School of Physics and Electronics, Hunan University, Changsha 410082, China

Received 3 November 2014; Revised 12 January 2015; Accepted 12 January 2015; 3 February 2015

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Surface plasmon polaritons (SPPs) are deriving from coupling of electromagnetic waves with oscillations of free electrons in a metal and supported at the insulator/metal interface [1]. Since SPPs can be confined in a subwavelength-scale area beyond the diffraction limit, they are regarded as information carriers for highly integrated photonic circuits, such as particles [2, 3], grooves [4-6], cavities [7-11], and apertures [12, 13]. In addition to the merit of miniaturization for these devices, their unique characteristics can satisfy various applications. Take the basic waveguide, for instance, directional couplers and splitters, and even modulators utilizing electro/thermoeffect [14, 15]. Typical plasmonic waveguides are the insulator-metal-insulator (IMI) waveguides and the metal-insulator-metal (MIM) waveguides [5, 16-21]. It is well known that the IMI waveguide exhibits less propagation loss, leading to longer propagation length than the MIM waveguide [16]. On the other hand, from a point of confining light, since the skin depth of SPPs fields into the dielectric medium is much longer than that into the metallic medium at the interfaces [19], the MIM waveguide is better than the IMI waveguide. Thus, the MIM waveguide is more suitable for high optical integration.

Specifically, a fine filter-like behavior was achieved in a Bragg reflector alternatively stacking MIM waveguides with two kinds of dielectrics [22]. For the sake of simplifying the fabrication of Bragg reflectors and decreasing the propagation length for SPPs, some simple plasmonic filters have been proposed. For example, the proposed structure based on MIM waveguides with a nanodisk resonator acts as a band-pass filter [23]. It has also been reported that the one-end-sealed MIM waveguide reveals fine filtering effect [24]. Subsequently, unidirectional excitation of SPPs in T-shaped waveguide with nanodisk resonator has been analyzed theoretically and numerically [25]. More recently, comb-like filters have been demonstrated [26]. These previous structures merely perform as filters. Thus, the investigation and design of a structure based on the MIM waveguides functioning not only as a filter but also as a bend waveguide are highly desired for various applications.

In this paper, we propose a novel structure based on the MIM waveguides directly connected to a nanodisk resonator with a metal-block core. In order to make our proposed structure be easily seized, we discuss the structure without the metal-block core at the beginning. The difference between simulation and theoretical results is related to the nanocavity not being sealed. Then we put the metallic block into the nanocavity. With the effective index method (EIM) and the resonant mode theory (RMT), we examine wavelength shifts of the resonant modes, which are also validated by finite-difference time-domain (FDTD) [27] method. Another result obtained by FDTD simulations reveals that regular shifts of resonant wavelengths can be manipulated by the block's rotation degree. Finally, when the width of the block changes from tens of nanometers to hundreds of nanometers, the proposed structure performs as a bend waveguide, not a filter.

2. Model Creation and Theories

Figure 1 illustrates the designed plasmonic filter which consists of two parts: two straight MIM waveguides (I, II) and a nanodisk cavity with a metal-block core. The two straight waveguides are connected to the cavity on the opposite sides, and the metallic block is located at the center of the cavity perpendicular to the straight waveguides. The width of straight waveguides [figure omitted; refer to PDF] is 50 nm, and radius of nanodisk [figure omitted; refer to PDF] is 200 nm. The height and width of the metallic block are labeled as [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , respectively, and the distance between the top of the block and the disk is [figure omitted; refer to PDF] . The medium (the white areas) filled in the disk cavity and two straight waveguides is assumed as air with the refractive index being 1 ( [figure omitted; refer to PDF] ). The metal (the green area) is silver, whose complex relative permittivity modeled by Drude equation can be expressed as follows [4]: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] represents the dielectric constant at infinite angular frequency with a value of 3.7, [figure omitted; refer to PDF] is the angular frequency of incident light, [figure omitted; refer to PDF] ( [figure omitted; refer to PDF] ) stands for the bulk plasma frequency, and [figure omitted; refer to PDF] ( [figure omitted; refer to PDF] ) means the damping frequency of the electron's oscillation [28].

Figure 1: Schematic diagram of the proposed plasmonic filter. [figure omitted; refer to PDF] : radius of nanodisk resonator, [figure omitted; refer to PDF] : height of metallic block, [figure omitted; refer to PDF] : width of metallic block, [figure omitted; refer to PDF] : width of Ag-air-Ag waveguide in the core of nanocavity, and [figure omitted; refer to PDF] : metal slit width.

[figure omitted; refer to PDF]

In this paper, we employ our 2D-FDTD program code to simulate the interaction between metal and incident waves. In our calculation, the current density is induced into Maxwell's equations since the real part of the permittivity of metal is negative. The grids are terminated by ten-cell-thick convolutional perfectly matched layer (i.e., CPML). The spatial mesh cells are set to [figure omitted; refer to PDF] nm, and the time step and the total number of time steps are taken as [figure omitted; refer to PDF] s and 2 × 10⁴ , respectively. The small size of cells and large total number of time steps ensure that the values of electromagnetic field and current density in spatial grids are convergent and stable. The fundamental TM mode of the MIM waveguide is excited by a plane wave source whose wavelength ranges from 600 nm to 1550 nm. One frequency-domain field and power monitor is set to detect the transmission power [figure omitted; refer to PDF] .

According to the Maxwell equations and boundary conditions, the effective refractive index of SPPs ( [figure omitted; refer to PDF] ) in the metal-insulator-metal (MIM) structure can be obtained as follows [17]: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] stands for the dielectric constant of the insulator and [figure omitted; refer to PDF] is the dielectric constant of the metal, which can be obtained from (1). [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are the transverse propagation constants of dielectric and metal, respectively. [figure omitted; refer to PDF] is the free space wave vector of the excitation light. Figure 2 reveals that [figure omitted; refer to PDF] relies on the width of the Ag-air-Ag waveguide in the core of the nanocavity [figure omitted; refer to PDF] and the wavelength [figure omitted; refer to PDF] . [figure omitted; refer to PDF] declines with the increase of [figure omitted; refer to PDF] at the same wavelength, while it will be of no alteration by increasing the incident wavelength with [figure omitted; refer to PDF] fixed. SPPs can propagate through the proposed structure when a stable standing wave is built up constructively in the cavity. In this case, the following resonant condition should be satisfied [29]: [figure omitted; refer to PDF] where positive integer [figure omitted; refer to PDF] is the number of antinodes of the standing SPPs wave. [figure omitted; refer to PDF] is the phase delay per round-trip in the cavity, and we have got [figure omitted; refer to PDF] . So we can get the resonant wavelength: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is incident wavelength and [figure omitted; refer to PDF] is the length per round trip. The resonant wavelength [figure omitted; refer to PDF] increases with the increase of effective index [figure omitted; refer to PDF] . [figure omitted; refer to PDF] increases with the decline of [figure omitted; refer to PDF] (i.e., the increase of [figure omitted; refer to PDF] ) at the same wavelength from Figure 2. Therefore, [figure omitted; refer to PDF] increases with the block's height increasing.

Figure 2: Influence of the width of the waveguide in the core of the nanocavity and wavelength on [figure omitted; refer to PDF] .

[figure omitted; refer to PDF]

3. Simulation Results and Discussion

To start with, the transmission properties of the structure without a metal-block core inside the cavity are studied to be compared with the proposed structure. Such a structure presents a filter-like behavior. The resonant wavelength of the ideal disk-shaped nanocavity [30, 31] can be obtained by the following equation: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the radius of the nanocavity. [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are the first kind Bessel function with the order [figure omitted; refer to PDF] and its derivative. [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are the first kind Hankel function with the order [figure omitted; refer to PDF] and its derivative, respectively. Figure 3 reveals that the wavelength shift of the resonant mode has approximately linear relation with the radius of nanocavity from both (5) and FDTD result. The unideal nanocavity in our proposed structure (two edges in the middle connected to two waveguides and then not being sealed) has its effective radius. The actual radius not effective radius is taken into (5), which leads to the high discrepancy up to 100 nm between simulation and theoretical results.

Figure 3: Resonant wavelength with respect to the radius of the cavity. The dot-dashed line is simulated from FDTD, and the solid line is calculated from (5).

[figure omitted; refer to PDF]

Next, we put a metallic block into the nanodisk. We investigate the influence of the height of the block on the resonant wavelengths firstly by using the FDTD method and then by using the RMT. The height of the block [figure omitted; refer to PDF] is variable and the width [figure omitted; refer to PDF] is set as 50 nm, while other parameters are fixed as above. Figure 4(a) shows the transmission spectra with the different heights of the block. From Figure 4(b) we can find that the resonant wavelengths tend to exhibit a red shift with the increase of the block's height, which is in accordance with (2)-(4). Magnetic-field [figure omitted; refer to PDF] distributions at [figure omitted; refer to PDF] nm with [figure omitted; refer to PDF] , [figure omitted; refer to PDF] nm with [figure omitted; refer to PDF] nm, and [figure omitted; refer to PDF] nm with [figure omitted; refer to PDF] nm are depicted in Figures 4(c)-4(e). We observe that the excitations of resonant modes emerge in the nanocavity under the three cases where the incident light can pass the MIM plasmonic structure. According to above simulations and analysis, it is obvious that we can achieve a fine filter-like behavior by tuning the height of the metallic block.

(a) Transmission response with respect to different heights of the block ( [figure omitted; refer to PDF] nm). (b) Resonant wavelength with respect to the height of the block. Contour profiles of field [figure omitted; refer to PDF] at (c) [figure omitted; refer to PDF] nm with [figure omitted; refer to PDF] nm, (d) [figure omitted; refer to PDF] nm with [figure omitted; refer to PDF] nm, and (e) [figure omitted; refer to PDF] nm with [figure omitted; refer to PDF] nm.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

(e) [figure omitted; refer to PDF]

Another parameter [figure omitted; refer to PDF] , which stands for the block's rotation degree, plays an important role in modulating the resonant wavelength. The rotation of the block leads to a displacement of the dielectric core, which makes the direction of magnetic oscillation variable. In addition, the change of the direction of magnetic oscillation makes [figure omitted; refer to PDF] different, affecting [figure omitted; refer to PDF] directly as shown in Figure 2. Hence, it is clear that the effective index changes when the symmetry of the structure is destroyed by modulating the block's rotation degree [figure omitted; refer to PDF] . The effective index [figure omitted; refer to PDF] is an important factor influencing resonant wavelength [figure omitted; refer to PDF] as shown in (4). Thus, we can manipulate the resonant wavelengths by changing the value of the block's rotation degree [figure omitted; refer to PDF] , which is verified by FDTD simulations. The structural parameters are set as [figure omitted; refer to PDF] nm and [figure omitted; refer to PDF] nm. Then we start to rotate the block counterclockwise. Figure 5(a) shows the transmission response corresponding to [figure omitted; refer to PDF] , 22.5°, 45°, 67.5°, and 90°. There is only one transmission peak at [figure omitted; refer to PDF] nm when [figure omitted; refer to PDF] is zero. When the rotation degree is [figure omitted; refer to PDF] , two transmission peaks appear around [figure omitted; refer to PDF] nm and [figure omitted; refer to PDF] nm, which are called mode 2 and mode 1, respectively. With increase of rotation degree [figure omitted; refer to PDF] , mode 1 exhibits a slight red shift, while mode 2 presents a blue shift (as shown in Figure 5(b)), leading to the wavelength gap between two modes getting larger. When [figure omitted; refer to PDF] is the value of 90°, only one transmission peak is observed at [figure omitted; refer to PDF] nm. Furthermore, Figures 5(c)-5(e) depict the field distributions of [figure omitted; refer to PDF] of transmission peaks at [figure omitted; refer to PDF] nm and [figure omitted; refer to PDF] nm, and transmission valley at [figure omitted; refer to PDF] nm when [figure omitted; refer to PDF] is the value of 45°, which match Figure 5(a) perfectly. Therefore, one can choose special rotation degree to design the filter with demanded channel.

(a) Transmission response with respect to [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] . (b) The wavelength at which the transmission spectrum peaks with respect to the block rotation degree. Contour profiles of field [figure omitted; refer to PDF] at (c) [figure omitted; refer to PDF] nm, (d) [figure omitted; refer to PDF] nm, and (e) [figure omitted; refer to PDF] nm with [figure omitted; refer to PDF] .

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

(e) [figure omitted; refer to PDF]

Finally, we find that the proposed structure changes from a filter to a waveguide under the special parameters. When we set the width of the block to be 250 nm and change the height from 210 nm to 250 nm, the structure performs as a bend waveguide, not as a filter. In Figure 6(a), we show the calculated transmission as a function of wavelength. We can see that the loss is about 10% at longer wavelength. The critical wavelength [figure omitted; refer to PDF] regarded as the original wavelength of the waveguide under which the transmission coefficient is max is 790 nm (823 nm, 845 nm, 903 nm, and 934 nm) for [figure omitted; refer to PDF] nm (220 nm, 230 nm, 240 nm, and 250 nm). Hence, it is clear that the wavelength [figure omitted; refer to PDF] increases monotonously with [figure omitted; refer to PDF] increasing as shown in Figure 6(b). That is, the operating wavelength range widens as [figure omitted; refer to PDF] decreases. The reason why our structure changes from a filter to a waveguide is that the nanocavity is connected to two MIM waveguides directly, leading to lots of pathways along which SPPs travel. Therefore, the resonant wavelengths differ slightly, which contributes to the change from a filter to a waveguide. For the sake of confirming our analysis, field distributions [figure omitted; refer to PDF] under some resonances are shown in Figures 6(c)-6(e) as examples.

(a) Transmission response with respect to different heights of the block ( [figure omitted; refer to PDF] nm). (b) The cut-off wavelengths with respect to the heights of the block. Contour profiles of field [figure omitted; refer to PDF] at (c) [figure omitted; refer to PDF] nm, (d) [figure omitted; refer to PDF] nm, and (e) [figure omitted; refer to PDF] nm (dashed red lines in (a)) with [figure omitted; refer to PDF] nm and [figure omitted; refer to PDF] nm.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

(e) [figure omitted; refer to PDF]

4. Conclusion

In this paper, the novel filter effect can be realized based on the MIM plasmonic nanodisk with a metal-block core. The structure without a block inside the cavity is studied to be compared with the proposed structure. The difference between simulation and theoretical results is related to the cavity not being sealed, while (5) is fit for ideal disk. The proposed structure reveals perfect filter effect by easily modulating the height of the block. The resonant wavelengths tend to exhibit a red shift with the increase of the height of the block. The result achieved by FDTD simulations can be accurately analyzed through the EIM and the RMT. We also find the proposed structure produces two modes (i.e., Mode 1 and Mode 2) through tuning the block's rotation degree by means of rotating the strip counterclockwise. Mode 1 and mode 2 tend to exhibit a red shift and a blue shift, respectively. Furthermore, when the width of the block is assumed as 250 nm and the height of the block changes from 210 nm to 250 nm, the proposed structure performs as a bend waveguide, not a filter. The FDTD results are verified by the theoretical analysis. Undoubtedly, the ways to give rise to filter effect are ingenious and wonderful. We believe that our fundamental study could find promising applications in the highly integrated plasmonic devices for optical processing and filtering.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant nos. 11074069 and 61176116) and the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant no. 20120161130003).

Conflict of Interests

The authors would like to declare that they do not have any commercial or associative interests that represent a conflict of interests in connection with the submitted paper.

References

[1] W. L. Barnes, A. Dereux, T. W. Ebbesen, "Surface plasmon subwavelength optics," Nature , vol. 424, no. 6950, pp. 824-830, 2003.

[2] S. A. Maier, P. G. Kik, H. A. Atwater, S. Meltzer, E. Harel, B. E. Koel, A. A. G. Requicha, "Local detection of electromagnetic energy transport below the diffraction limit in metal nanoparticle plasmon waveguides," Nature Materials , vol. 2, no. 4, pp. 229-232, 2003.

[3] W. H. Weber, G. W. Ford, "Propagation of optical excitations by dipolar interactions in metal nanoparticle chains," Physical Review B-Condensed Matter and Materials Physics , vol. 70, no. 12, 2004.

[4] X.-S. Lin, X.-G. Huang, "Tooth-shaped plasmonic waveguide filters with nanometeric sizes," Optics Letters , vol. 33, no. 23, pp. 2874-2876, 2008.

[5] S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, T. W. Ebbesen, "Channel plasmon-polariton guiding by subwavelength metal grooves," Physical Review Letters , vol. 95, no. 4, 2005.

[6] K. H. Wen, L. S. Yan, W. Pan, B. Luo, Z. Guo, Y. H. Guo, X. G. Luo, "Spectral characteristics of plasmonic metal-insulator-metal waveguideswith a tilted groove," IEEE Photonics Journal , vol. 4, no. 5, pp. 1794-1800, 2012.

[7] P.-H. Lee, Y.-C. Lan, "Plasmonic waveguide filters based on tunneling and cavity effects," Plasmonics , vol. 5, no. 4, pp. 417-422, 2010.

[8] Y. H. Guo, L. S. Yan, W. Pan, B. Luo, K. Wen, Z. Guo, H. Li, X. Luo, "A plasmonic splitter based on slot cavity," Optics Express , vol. 19, no. 15, pp. 13831-13838, 2011.

[9] K. H. Wen, L. S. Yan, W. Pan, B. Luo, Z. Guo, Y. H. Guo, "A four-port plasmonic quasi-circulator based on metal-insulator-metal waveguides," Optics Express , vol. 20, no. 27, pp. 28025-28032, 2012.

[10] F. F. Hu, H. X. Yi, Z. P. Zhou, "Wavelength demultiplexing structure based on arrayed plasmonic slot cavities," Optics Letters , vol. 36, no. 8, pp. 1500-1502, 2011.

[11] M. A. Bavil, L. Gao, X. D. Sun, "A compact nanoplasmonics filter and intersection structure based on utilizing a slot cavity and a Fabry-Perot Resonator," Plasmonics , vol. 8, no. 2, pp. 631-636, 2013.

[12] P. Lalanne, J. P. Hugonin, J. C. Rodier, "Theory of surface plasmon generation at nanoslit apertures," Physical Review Letters , vol. 95, no. 26, 2005.

[13] M. J. Kofke, D. H. Waldeck, Z. Fakhraai, S. Ip, G. C. Walker, "The effect of periodicity on the extraordinary optical transmission of annular aperture arrays," Applied Physics Letters , vol. 94, no. 2, 2009.

[14] T. Nikolajsen, K. Leosson, S. I. Bozhevolnyi, "Surface plasmon polariton based modulators and switches operating at telecom wavelengths," Applied Physics Letters , vol. 85, no. 24, pp. 5833-5835, 2004.

[15] A. Boltasseva, T. Nikolajsen, K. Leosson, K. Kjaer, M. S. Larsen, S. I. Bozhevolnyi, "Integrated optical components utilizing long-range surface plasmon polaritons," Journal of Lightwave Technology , vol. 23, no. 1, pp. 413-422, 2005.

[16] P. Berini, "Long-range surface plasmon polaritons," Advances in Optics and Photonics , vol. 1, no. 3, pp. 484-588, 2009.

[17] J. T. Kim, J. J. Ju, S. Park, M.-S. Kim, S. K. Park, M.-H. Lee, "Chip-to-chip optical interconnect using gold long-range surface plasmon polariton waveguides," Optics Express , vol. 16, no. 17, pp. 13133-13138, 2008.

[18] T. Rosenzveig, P. G. Hermannsson, K. Leosson, "Modelling of polarization-dependent loss in plasmonic nanowire waveguides," Plasmonics , vol. 5, no. 1, pp. 75-77, 2010.

[19] S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, T. W. Ebbesen, "Channel plasmon-polariton guiding by subwavelength metal grooves," Physical Review Letters , vol. 95, no. 4, 2005.

[20] S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, J. Y. Laluet, T. W. Ebbesen, "Channel plasmon subwavelength waveguide components including interferometers and ring resonators," Nature , vol. 440, no. 7083, pp. 508-511, 2006.

[21] Z. Y. Fang, H. Qi, C. Wang, X. Zhu, "Hybrid plasmonic waveguide based on tapered dielectric nanoribbon: excitation and focusing," Plasmonics , vol. 5, no. 2, pp. 207-212, 2010.

[22] J. Park, H. Kim, B. Lee, "High order plasmonic Bragg reflection in the metal-insulator-metal waveguide Bragg grating," Optics Express , vol. 16, no. 1, pp. 413-425, 2008.

[23] H. Lu, X. M. Liu, D. Mao, L. R. Wang, Y. K. Gong, "Tunable band-pass plasmonic waveguide filters with nanodisk resonators," Optics Express , vol. 18, no. 17, pp. 17922-17927, 2010.

[24] Z. Zhang, H. Wang, Y. Zhao, D. Lu, "Transmission properties of the one-end-sealed metal-insulator-metal waveguide," Optik , vol. 124, no. 2, pp. 177-179, 2013.

[25] G. X. Wang, H. Lu, "Unidirectional excitation of surface plasmon polaritons in T-shaped waveguide with nanodisk resonator," Optics Communications , vol. 285, no. 21-22, pp. 4190-4193, 2012.

[26] K. Wen, L. Yan, W. Pan, B. Luo, Z. Guo, Y. Guo, X. Luo, "Design of plasmonic comb-like filters using loop-based resonators," Plasmonics , vol. 8, no. 2, pp. 1017-1022, 2013.

[27] A. Taflove, S. C. Hagness Computational Electrodynamics: The Finite-Difference Time-Domain Method , Artech House, Boston, Mass, USA, 2000., 2nd.

[28] Z. Han, E. Forsberg, "Surface plasmon bragg gratings formed in metal-insulator-metal waveguides," IEEE Photonics Technology Letters , vol. 19, no. 2, pp. 91-93, 2007.

[29] Q. Zhang, X.-G. Huang, X.-S. Lin, J. Tao, X.-P. Jin, "A subwavelength coupler-type MIM optical filter," Optics Express , vol. 17, no. 9, pp. 7549-7554, 2009.

[30] S.-L. Qiu, Y.-P. Li, "Q-factor instability and its explanation in the staircased FDTD simulation of high-Q circular cavity," Journal of the Optical Society of America B: Optical Physics , vol. 26, no. 9, pp. 1664-1674, 2009.

[31] I. Chremmos, "Magnetic field integral equation analysis of interaction between a surface plasmon polariton and a circular dielectric cavity embedded in the metal," Journal of the Optical Society of America A , vol. 26, no. 12, pp. 2623-2633, 2009.