Metal-insulator-metal (MIM) surface plasmon polaritons (SPPs) waveguides with side-coupled resonators have been widely studied through various approaches. However, few methods are both physically transparent and complete. Here, an analytical approach, which is based on the Green’s function method, is developed in order to investigate electromagnetic wave transmission across SPPs MIM waveguide networks. The proposed method is applied in order to model different MIM-waveguide geometries with weakly-coupled side stubs, comparing to the geometries with strongly-coupled stubs. The weak coupling between the backbone and stubs is taken into account by the electromagnetic field leakage at metal-insulator interface. Analytical expressions for transmittance in cases of single stub and cavity are obtained straightforwardly. Our method shows excellent computational efficiency in contrast with solving Maxwell equations numerically.
© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
Surface plasmon polaritons (SPPs) waveguides with coupled resonators have been widely studied in recent years because of its potential application in highly integrated optical circuits. For example, a designed plasmonic waveguide system with double resonant nano-disks shows a powerful electromagnetically induced transparency (EIT) like effect . Other optical effects such as wavelength splitter , total absorption , slow light effect , and plasmon-induced transparency [5,6] have been realized in plasmonic waveguide with nanoscale side-coupled resonators. A diverse geometry of resonators side-coupled to metal-insulator-metal (MIM) SPPs waveguide have been reported, such as ring resonator [6,7], serial stubs [8,9], rectangular cavity [10,11], as well as irregular cavities . Among them, rectangular side-coupled resonators have the advantage of taking up less physical space. Simplicity in geometry makes rectangular resonators suffer less challenge of fabrications. Recently, various potential applications have been exhibited by MIM SPPs networks with rectangular resonators, including optical switching [13,14], wavelength filter [15–17], optical pressure sensor , microfluidic sensor , and all-optical logic gates .
Methods of treating plasmonic metallic MIM waveguides are either numerical simulation or analytical algorithm. The numerical tools, like Finite Element Method (FEM) and Finite Difference Time Domain (FDTD), are widely used to deal with various SPPs networks. However, generally, they not only lack physical transparency but also are consuming in time and memory. Besides, the use of Cartesian-grid in FDTD may produce artificial surface plasmon around the metallic structures for subwavelength cells , resulting in obscures of investigation. In contrast to simulation method, analytical approaches are efficient and rapid in predicting the behavior of light modes in the waveguides . Conventional analysis methods, mainly containing Coupled-Mode Theory (CMT), and Transfer Matrix Method (TMM), have been successfully applied to model various guided-wave devices, such as optical directional couplers , nanodisks coupled waveguides [1,24], Kretschmann prism , channel fiber , etc. But the usage of TMM is limited to 1D model up to now. In prior studies [1,27], however, perfect conductor condition is adopted at the terminal of stubs, which failed to reflect the real leakage of the electromagnetic field at metal-insulator interface. Further, it is attributed to the zero value of the penetration depth of the perfect conductor that the unconnected cavity structure is always troublesome. As an engineering approach, Transmission Line (TL) theory has enjoyed large popularity in modeling varieties of plasmonic devices, including wavelength demultiplexers [28–30], optical filters [31–33], all-optical switching [34,35], sensor , resonator [13,36], etc. Importantly, the analogy of waveguide impedance is the ratio of voltage and travelling current,, based on the circuit approach . However, the TL theory fails to demonstrate field distribution of arbitrary configuration. Therefore, we find that the analysis approach of electromagnetic wave propagation remains to be improved.
In this paper, we propose a consistent analytical approach based on Green’s function method to model the MIM SPPs waveguides networks containing rectangle side-coupled antennas. In our method, propagating behaviors of waves within the considered system are determined by two key physical quantities, wave vector and effective impedance, which are depended on material parameters and magnitude of structure. We map out the complete theoretical framework of our method in details. Algorithms for constructing Green’s function matrix in the cases of arbitrary stubs, cavities and junction networks respectively are also developed. Thereby, we derive analytical expressions for the transmission spectra of the waveguide coupled to single stub and single cavity respectively. Magnetic-field distribution for the two specifications at given frequency are demonstrated both analytically and numerically. It shows that the analytical model agree with the FDTD method very well, which validates the feasibility of the theoretical analysis. Moreover, a multi-stub antenna structure is studied to show the high efficiency of the analytical model on the premise of high accuracy. Our work is expected to give deep insight into the optical properties of the MIM waveguide networks and can provide useful guidance for the design of integrated optical circuits.
Consider an infinite two-dimensional MIM waveguide, as shown in Fig. 1(a), which is composed of two semi-infinite metallic media and a thin gap insulator layer with thickness. We know that only the TM modes can be excited in the MIM waveguide under the boundary condition for the dielectric-metal interface, and their propagation equation can be expressed as :37], the impedance has been derived and given approximately by the first-order Taylor expansion as : One can easily confirm that is given by:
Before we employ the Green’s function of the magnetic field, it is necessary to seek the propagation wave vector of the MIM waveguide. For the only existing symmetric mode TM0 in the waveguide we focus in this paper, when dielectric thickness is small enough (), the SPP propagation wave vector is determined by the dispersion relation :
In order to characterize the leakage of the field which is unavoidable in SPP MIM networks, here we introduce the conception of metallic waveguides. This is particularly important when solving the problem containing side-coupled antennas. Since two key physical quantities, wave vector and effective impedance, have been obtained in MIM waveguides, we simply consider an infinite metal media, as shown in Fig. 1(b). The thickness of the metallic layer is assumed to be, which is arbitrary but has a strong relation to the structure when particular geometry is concerned. Thus we get:Fig. 1(c), which is composed of an end-terminated MIM waveguide and an infinite metallic block media. As we see from the schematic, metallic media acts as obstacle occurring at the end of the semi-infinite MIM waveguide. The leakage caused by the block metallic material is represented by a semi-infinite metallic waveguide we assumed. In this way, such stub-like structure can be decomposed into a truncated MIM waveguide and a semi-infinite metallic waveguide. By letting, the inverse of Green’s function is:
In what follows, we present the basic idea of the Green’s function method in the MIM network waveguide system. Any composite structure can be decomposed into the following categories, as shown in Fig. 2. They are (a) input and output waveguides, (b) short-circuit waveguide with finite length, (c) open-circuit stub like structure waveguide and (d) unconnected cavity side-coupled structure. Input and output waveguides act as exterior ports, which can be taking as semi-infinite waveguide. Let indicate the domain and denote all the existing interfaces for each constituent. Based on the interface response theory , one can calculate the Green’s function of each composited system which is related to the Green’s function by:Fig. 2(b)], bounded by two free surfaces at and such that electric field at these boundaries, in this case, the cleavage operator is:Fig. 2(a)], electric field E vanishing on its extremity, the cleavage operator is:Fig. 2(c), we take penetration depth of metal into account to represent the insulator-metal interface across the propagation direction. Thus we get the inverse of the surface Green’s function:
Suppose an incident wave is launched in one homogeneous piece of composite material, the magnetic-field distribution along the wave propagation direction in the waveguide can be calculated by :Eq. (18) and obtain the distribution of an arbitrarily transmit domain with finite length as where is the impedance of the input waveguide. Finally, transversal magnetic-field distribution can be expressed in a generalized form as:Eq. (20), which is the superposition of optical field magnitude of four diverse propagation processes during the 1D material (0→1). Among them, and denote the time harmonic wave field transmitting through the first interface (0), which is indicated by . While and describe the propagation process occurring reflection at the second interface (1), which is indicated by. By a further observation, oscillation occurs just after wave propagating through interface (0) in the process of , without reaching at the other interface yet. While, represents the wave oscillation with a reflection at the second interface (1). In addition, the propagation process described by is a reflecting oscillation at the second interface (1). Finally shows the propagation process that optical field is reflected back and forth twice, the first occurs at interface (1) while the second at the interface (0).
As the transversal field distribution can be calculated by Eq. (19), the longitudinal distribution of the magnetic field should be also considered thus magnetic-field spatial distribution in the waveguide can be obtained. The infinite MIM waveguide structure is divided into three parts denoted by I, II and III, as shown in Fig. 1(a). For the TM0 plasmonic mode, by using the boundary connection condition, the general formula of magnetic-field spatial distributions in the MIM waveguide is given byFig. 1(b, c), the magnetic-field spatial distributions in the propagating range is:
Suppose the wave has been launched in the port 1 of the SPP network, the transmittance of the port and reflectance can be obtained directly in terms of the green function:
Thereby we get the two deterministic physical quantities in wave propagation of both MIM and metallic waveguides. In implementing of side-coupled antennas structure, we assume that stubs are terminated by semi-infinite metallic output waveguides. Similarly for cavities, metallic waveguides are supposed within the interconnecting direction with a finite length between the cavity and the backbone waveguide. Armed with these results, we can now set up network analogs for any system including interested-guides, stubs, cavities, etc.-interconnected in arbitrary ways and excited either externally or internally. Since precise numerical values are assigned to all the quantities involved, we can then calculate frequency response, radiated power, and other quantities of concern.
3. Main constructions in SPPs networks
A network structure consisted of four X-junctions whose corresponding can be written as a (12 × 12) matrix, has already been discussed in details . However, the analysis on such problem is limited to short-circuit waveguides which is directly connected to each other composing the network structure with several junctions. The previous framework lacked of skills to model the side-coupled antennas structures, such as T-shaped waveguides [9,35] and unconnected cavities [10,36]. In order to mature the Green’s function method into a complete way for various SPPs network compositions, we introduce an explicit implemental algorithm of the Green function to deal with three typical constructions, namely (i) open circuit stubs, (ii) side-coupled cavities, and (iii) short circuit junctions.
3.1 Open-circuit stub model
By assuming metallic waveguides connected to the ends of the stub structure, arbitrarily side-coupled stubs system can be analyzed by the Green’s function method. In this manner, the side stub-like open circuit structure is considered as two short-circuit junction with distinct impedance for MIM waveguides and for metallic waveguides. A general case is that one backbone MIM waveguide is coupled with open circuit stubs. As the schematic illustrated in Fig. 3(a), first, we label the open circuit stubs by an array of bivariate functions, where refers to theT-junction on the backbone SPP waveguide and the binary variable with denotes the width and length of the stubs respectively. is the distance between and label of the bus waveguide, where .
Then we focus on the Green’s function, which is considered to be a matrix. Two kinds of interfaces occur in such configuration. They are representing the T-junctions on the backbone waveguide, and representing the semi-infinite metallic waveguides. The diagonal matrix elements representing self-energy of each labeled stub open circuits can be expressed as:
3.2 Arbitrary cavity model
In this part, we emphasize on the side-coupled cavities. Owning the view of metallic waveguide, we consider the gap between backbone waveguide and cavities to be metallic waveguides with loss. In detail, one backbone waveguide is coupled with cavities, which is shown in Fig. 3(b). Similar labeling method, i.e., is utilized to characterize the geometry of the proposed structure. To deal with such structure, the Green’s function should be a matrix. Additional geometry parameter is the gap distance between the cavity and the bus waveguide. Because of the appearance of the gap which constitutes a lossy metallic junction, the interface occurs in such structure. represents the position of each cavity on the backbone waveguide, while is the semi-infinite metallic waveguide we supposed, which is equivalent to in the case of general side stubs configuration. Following the matrix elements expression above, we get the Green’s function straightforward.
with The lossy metallic waveguide gap is denoted as where And where And where
- (2) Non-diagonal matrix elements:
For two neighbor junctions on the bus waveguide, we obtain:
For each gap, we present the correlation betweenand as:
For each cavity, we have the correlation between and as:
3.3 Short-circuit junction networks
Such networks contain orthogonally connected MIM-waveguide forming junctions separately at their terminals as shown in Fig. 3(c). Clearly, four short-circuit MIM waveguides with length and width, constitute a junction structure. Our networks constructed by short circuits junctions are compatible to any number of arms, junctions and ports. The inverse of Green’s function concerning short-circuit junction is a matrix. Junction dots are labeled by a pair of numbers, and can be divided into three types. The first is the four junction dots at the corner composed with two semi-infinite waveguides and two finite waveguides (red colored), second type is thejunction dots on the side way with one semi-infinite waveguide and three finite waveguides (blue colored), and the last one is the rest inside junctions (black colored) with four finite waveguides. For the junctions with two semi-infinite waveguide terms, i.e. the junctions on the corner, their inverse of the surface Green’s function are as follows:
Though, of course, it is necessary to state the limitation of the analytical method here. The distance between side antennas or two paralleled waveguides cannot be smaller than the width of the waveguide. Further, the gap distance between the side cavity antennas cannot be too large. Thus far, we have derived the general expression for each element of the green’s function. In this way, transmission properties of plasmonic waveguide system contained by any combination of junctions, stubs or cavities can be characterized analytically. The theoretical analysis mentioned above may provide a guideline for the accurate adjustment of spectral responses.
4. Validity of method
4.1 Single-side antenna
In order to validate our analytical method and highlight the improvement of the coupling efficiency by the metallic waveguides, we calculate transmittance spectra and the field distribution for (i) single stub [Fig. 4] and (ii) single cavity [Fig. 5], both analytically and numerically. For the sake of simplicity, vacuum is used as the insulating media, thus. The metal is chosen by silver. Its permittivity can be expressed by the well-known Drude model :41]. As the simplest example schematically shown in Fig. 4(a), a single side-coupled stub structure is analyzed by the proposed model. The widths of the backbone waveguide and the side-coupled stub are denoted by and, respectively. The length of stub and the backbone waveguide are denoted by and, respectively. By applying Eq. (29), the transmission at the waveguide out port is obtained by:Fig. 4(b), where the transmission coefficient is plotted versus the free space wavelength (solid) and is compared with the FDTD result (dotted). For reference, we also show the transmittance spectrum reported in  with the same geometric and material parameters. Obviously, our transmission spectrum does a better match with the numerical result, especially at low frequency range. Besides, the magnetic field in the structure at the incident wavelength is shown in Fig. 4(c). Comparing with the insert figure calculated by numerical method, we see close agreement between them.
Transmission properties of the MIM waveguide side-coupled to single unconnected cavity structure are also studied, shown as Fig. 5(a). The additional geometry parameter is the gap distance. With Eq. (29) we get the transmission as:Fig. 5(b), the magnetic-field distributions are shown, and the simulation (insert figure) is compared with the result from our theory. Both of them are of perfect agreement with each other. To measure the effect of the length modification, we study the transmission as a function of the gap distancein Fig. 5(c). The red dashed curves show the simulation data of the transmittance, while the blue block lines result from our improved analytical model. For comparison, the black dash lines show the transmission properties without length modification. It is obvious that the analytical results using the proposed method agree well with the FDTD results, particularly at the transmission-dip wavelength where the proposed length modification takes great effect.
4.2 Multi-stub antenna
To show the efficiency in computation, we further consider a MIM waveguide with four periodical stubs perpendicular to it as shown in Fig. 6(a). The insulator medium of the waveguide is assumed to be air and the metal medium is chosen by silver. The length of the stubs and the interval between two stubs are denoted by and, respectively. Figure 6(b) gives the transmission spectra for such serial stubs structure withand by FDTD (dotted, red) and Green function method (solid, blue). It shows that the analytical method provides more detailed features in transmission on the premise of high accuracy. Compared to the case of single stub, taking 1.544s, the computing time of the transmission of the four periodical stubs structure is 2.260s. Swept range is from 500nm to 2000nm with step length 1nm. Resulting from the structural complexity, the increase in computing time is acceptable. The only difference is the increase in matrix dimensions, from 4 × 4 (single stub) to 10 × 10 (four stubs). In comparison, it takes 138s (2 minutes, 18 seconds) and 303s (5 minutes, 3 seconds) for Comsol to calculate the transmission of single stub structure and four periodical stubs structure respectively, with the same swept range and step length. Thus it can be seen that simulation tool is far more time-consulting than the method of Green’s function, which is very efficient and rapid to predict the behavior of light modes in waveguides.
We have developed a complete analytical method based on Green function to calculate the propagation problem in the MIM SPPs waveguide networks containing rectangle side-coupled antennas. Our method provides analytical solutions for arbitrary composition of side-coupled stubs or cavities. In proposed method, phase shift effect caused by metals-insulator interface has been taken into account by utilizing length modification of the stubs or cavities in propagating direction. The leakage of the electromagnetic field at metal-insulator interface described by metallic waveguides reveals the weak coupling between the backbone and antennas structure. Using the proposed theory, we demonstrated the algorithms for constructing Green’s function matrix in the cases of arbitrary stubs, cavities as well as junction networks, which are important in practical. Thus the effectiveness of our method stemming from superposition arithmetic in Green function formalism is estimated. Through the calculation of the transmission spectra and magnetic field distribution in cases of single stub and cavity, it is shown that our method is capable of providing results in high accuracy. The proposed analytical method can significantly reduce the computational time and allow us to design geometrically flexible couplers for plasmonic waveguides. Our work may help to understand the fundamental physics of the MIM waveguide networks and to guide the corresponding applications.
National Key R&D Program of China (Grant No. 2017YFA0303400) and National Natural Science Foundation of China (NSFC) (Grant No. 91630313).
1. G. Lai, R. Liang, Y. Zhang, Z. Bian, L. Yi, G. Zhan, and R. Zhao, “Double plasmonic nanodisks design for electromagnetically induced transparency and slow light,” Opt. Express 23(5), 6554–6561 (2015). [CrossRef] [PubMed]
2. K. Wen, Y. Hu, L. Chen, J. Zhou, L. Lei, and Z. Guo, “Design of an Optical Power and Wavelength Splitter Based on Subwavelength Waveguides,” J. Lightwave Technol. 32(17), 3020–3026 (2014). [CrossRef]
4. C. Li, R. Su, Y. Wang, and X. Zhang, “Theoretical study of ultra-wideband slow light in dual-stub-coupled plasmonic waveguide,” Opt. Commun. 377, 10–13 (2016). [CrossRef]
6. M. R. Rakhshani and M. A. Mansouri-Birjandi, “Dual wavelength demultiplexer based on metal–insulator–metal plasmonic circular ring resonators,” J. Mod. Opt. 63(11), 1078 (2016). [CrossRef]
7. S. Zhan, H. Li, G. Cao, Z. He, B. Li, and H. Yang, “Slow light based on plasmon-induced transparency in dual-ring resonator-coupled MDM waveguide system,” J. Phys. D Appl. Phys. 47(20), 205101 (2014). [CrossRef]
9. R. Zafar and M. Salim, “Wideband Slow Surface Plasmons in Double Resonator Plasmonic Grating Waveguide,” IEEE Photonic. Tech. L 26(22), 2221–2224 (2014). [CrossRef]
10. Z. He, H. Li, S. Zhan, B. Li, Z. Chen, and H. Xu, “π-Network transmission line model for plasmonic waveguides with cavity structures,” Plasmonics 10(6), 1581–1585 (2015). [CrossRef]
11. M. A. Swillam and A. S. Helmy, “Feedback effects in plasmonic slot waveguides examined using a closed form model,” IEEE Photonic. Tech. L 24(6), 497–499 (2012). [CrossRef]
12. X. Zhou, M. Ouyang, B. Tang, Z. Wang, and J. He, “Transparency windows of the plasmonic nanostructure composed of C-shaped and U-shaped resonators,” Opt. Commun. 384, 65–70 (2017). [CrossRef]
13. A. N. Taheri and H. Kaatuzian, “Numerical investigation of a nano-scale electro-plasmonic switch based on metal-insulator-metal stub filter,” Opt. Quantum Electron. 47(2), 159–168 (2015). [CrossRef]
14. A. Noual, O. E. Abouti, E. H. El Boudouti, A. Akjouj, Y. Pennec, and B. Djafari-Rouhani, “Plasmonic-induced transparency in a MIM waveguide with two side-coupled cavities,” Appl. Phys., A Mater. Sci. Process. 123(1), 49 (2017). [CrossRef]
15. S. M. Ebadi, M. S. Bayati, S. Bonyadi Ram, S. M. Poursajadi, and M. Jamili, “High-efficiency nanoplasmonic wavelength filters based on MIM waveguides,” IEEE Photonic. Tech. L 28(22), 2605–2608 (2016). [CrossRef]
16. S. Elbialy, B. Yousif, and A. Samra, “Modeling of active plasmonic coupler and filter based on metal-dielectric-metal waveguide,” Opt. Quantum Electron. 49(4), 145 (2017). [CrossRef]
17. X. Lin and X. Huang, “Numerical modeling of a teeth-shaped nanoplasmonic waveguide filter,” J. Opt. Soc. Am. B 26(7), 1263–1268 (2009). [CrossRef]
18. J. Wu, P. Lang, X. Chen, and R. Zhang, “A novel optical pressure sensor based on surface plasmon polariton resonator,” J. Mod. Opt. 63(3), 219–223 (2016). [CrossRef]
19. X. Li, J. Song, and J. X. J. Zhang, “Design of terahertz metal-dielectric-metal waveguide with microfluidic sensing stub,” Opt. Commun. 361, 130–137 (2016). [CrossRef]
20. H. H. Wu, B. H. Cheng, and Y. C. Lan, “Coherent-controlled all-optical devices based on plasmonic resonant tunneling waveguides,” Plasmonics 12(6), 1–7 (2016).
21. R. Fikri and J. P. Vigneron, “Discrete plasmonic nano-waveguide: numerical and theoretical studies,” Proc. SPIE 6343, Photonics North 2006, 63432Q (2006).
22. I. Haddouche and L. Cherbi, “Comparison of finite element and transfer matrix methods for numerical investigation of surface plasmon waveguides,” Opt. Commun. 382, 132–137 (2017). [CrossRef]
23. L. Y. He, T. J. Wang, Y. P. Gao, C. Cao, and C. Wang, “Discerning electromagnetically induced transparency from Autler-Townes splitting in plasmonic waveguide and coupled resonators system,” Opt. Express 23(18), 23817–23826 (2015). [CrossRef] [PubMed]
24. H. Lu, X. Liu, and D. Mao, “Plasmonic analog of electromagnetically induced transparency in multi-nanoresonator-coupled waveguide systems,” Phys. Rev. A 85(5), 053803 (2012). [CrossRef]
25. H. R. Gwon and S. H. Lee, “Spectral and angular responses of surface plasmon resonance based on the Kretschmann prism configuration,” Mater. Trans. 51(6), 1150–1155 (2010). [CrossRef]
26. S. Zhan, H. Li, G. Cao, Z. He, B. Li, and H. Xu, “Theoretical analysis of plasmon-induced transparency in ring-resonators coupled channel drop filter systems,” Plasmonics 9(6), 1431–1437 (2014). [CrossRef]
29. H. Lu, X. Liu, L. Wang, D. Mao, and Y. Gong, “Nanoplasmonic triple-wavelength demultiplexers in two-dimensional metallic waveguides,” Appl. Phys. B 103(4), 877–881 (2011). [CrossRef]
30. G. Wang, H. Lu, X. Liu, D. Mao, and L. Duan, “Tunable multi-channel wavelength demultiplexer based on MIM plasmonic nanodisk resonators at telecommunication regime,” Opt. Express 19(4), 3513–3518 (2011). [CrossRef] [PubMed]
31. A. N. Taheri and H. Kaatuzian, “Numerical investigation of a nano-scale electro-plasmonic switch based on metal-insulator-metal stub filter,” Opt. Quantum Electron. 47(2), 159–168 (2015). [CrossRef]
32. S. Wu and D. Wu, “Adjusting Spectrum of MIM Optical Filters by Stub Inclination,” in Frontiers in Optics 2016, OSA Technical Digest (online) (Optical Society of America, 2016), paper JTh2A.157.
33. S. M. Ebadi, M. S. Bayati, S. Bonyadi Ram, S. M. Poursajadi, and M. Jamili, “High-Efficiency Nanoplasmonic Wavelength Filters Based on MIM Waveguides,” IEEE Photonic. Tech. L 28(22), 2605–2608 (2016). [CrossRef]
34. S. Paul and M. Ray, “Analysis of plasmonic subwavelength waveguide-coupled nanostub and its application in optical switching,” Appl. Phys., A Mater. Sci. Process. 122(1), 21 (2016). [CrossRef]
35. S. Paul and M. Ray, “Plasmonic switching and bistability at telecom wavelength using the subwavelength nonlinear cavity coupled to a dielectric waveguide: A theoretical approach,” J. Appl. Phys. 120(20), 203102 (2016). [CrossRef]
36. H. Yang, G. Li, X. Su, W. Zhao, Z. Zhao, X. Chen, and W. Lu, “A novel transmission model for plasmon-induced transparency in plasmonic waveguide system with a single resonator,” RSC Advances 6(56), 51480–51484 (2016). [CrossRef]
37. K. Zhang and D. Li, Electromagnetic Theory for Microwaves and Optoelectronics, 2nd ed. (Springer, Berlin, 2008).
38. J. Vasseur, A. Akjouj, L. Dobrzynski, B. Djafarirouhani, and E. Elboudouti, “Photon, electron, magnon, phonon and plasmon mono-mode circuits,” Surf. Sci. Rep. 54(1), 1–156 (2004). [CrossRef]
39. Q. Zhu and Z. Wang, “The Green’s function method for metal-dielectric-metal SPP waveguide network,” EPL 103(1), 17004 (2013). [CrossRef]
40. S. A. Maier, Plasmonics: fundamentals and applications (Springer Science & Business Media 2007).
41. J. Tao, X. G. Huang, X. Lin, Q. Zhang, and X. Jin, “A narrow-band subwavelength plasmonic waveguide filter with asymmetrical multiple-teeth-shaped structure,” Opt. Express 17(16), 13989–13994 (2009). [CrossRef] [PubMed]
42. A. Pannipitiya, I. D. Rukhlenko, M. Premaratne, H. T. Hattori, and G. P. Agrawal, “Improved transmission model for metal-dielectric-metal plasmonic waveguides with stub structure,” Opt. Express 18(6), 6191–6204 (2010). [CrossRef] [PubMed]