Applications of Nano-Scale Plasmonic Structures in Design of Stub Filters — A Step Towards Realization of Plasmonic Switches

As electronic devices and circuits were shrinking to the nano-scale chips, some drawbacks hindered the reaching to the speed higher than tens of Giga Hertz, such as the higher power consumption, delay, and interference of the signals. However, the need for transmitting the huge amount of data over communication networks urged everyone in this area of technology to find an alternative for pure electronic devices. Despite rapidly developing the photonics technology, in practice the implementation of the photonic counterparts of the electronic devices encountered many problems such as inability to become integrated. The size of photonic devices halted on the half the wavelength of the operating signal due to the Diffrac‐ tion limit of light.


Introduction
As electronic devices and circuits were shrinking to the nano-scale chips, some drawbacks hindered the reaching to the speed higher than tens of Giga Hertz, such as the higher power consumption, delay, and interference of the signals. However, the need for transmitting the huge amount of data over communication networks urged everyone in this area of technology to find an alternative for pure electronic devices. Despite rapidly developing the photonics technology, in practice the implementation of the photonic counterparts of the electronic devices encountered many problems such as inability to become integrated. The size of photonic devices halted on the half the wavelength of the operating signal due to the Diffraction limit of light.
In recent years, Surface Plasmon Polaritons (SPPs) have been considered as one of the most promising ways to overcome the diffraction limits of photonic devices. The field of manipulating the SPPs as the carrier of the signal is known as Plasmonics. In plasmonics, the advantages of electronics such as nano-scale component design are joint together with the benefits of photonics. Therefore, plasmonics shapes a key part of the fascinating field of nanophotonics, which discovers how electromagnetic fields can be confined over dimensions on the order of or smaller than the wavelength. It is constructed based on interaction of electromagnetic radiation and conduction electrons at metallic interfaces or in small metallic nanostructures, leading to enhance the confinement of optical field in a sub-wavelength dimension.
The strong interaction between microscopic metal particles and light has been utilized for thousands of years, the Lycurgus Roman Cup dating from 4'th century A.D. [1], brightly colored stained-glass window panels by annealing metallic salts in transparent glass [1]. The practical instructions for this ancient "nanofabrication" technique survive from as early as the 8th century A.D. [1]. From the beginning of the 20 th century many phenomena related to the unusual optical properties of particles and spheres were explained; the blue color of the sky in terms of a simple derivation of the scattering power of spheres by Rayleigh [2], the bright color of metal glass by Garnett [3], presentation of a general formulation for the scattering of light from spherical surface by Gustav Mie [4]. In 1957, Ritchie [5] described surface plasmons in thin films in terms of electron energy loss spectroscopy. Otto [6], Kretschmann and Raether [7] explained the procedures for exciting surface plasmons on thin films optically in 1968. In 1974, the enhancement of Raman scattering from molecules influenced by the enhanced local fields at a rough metal surface was first observed by Fleischmann et al. [8].
In recent years, the applications of the surface plasmons are extended in optical integrated circuits and create a rapid developing field of photonics known as Plasmonics. Many researches are being carried out on developing plasmonic structures. Plasmonic switches and modulators are the most important components used for light routing and switching in the rapidly developing area of broadband optical communications and designing ultra-high speed switching devices, where sub-femtosecond interaction of light with matter is our main design concern [9]. Numerous plasmonic devices and components are introduced and huge of research are published. By exciting of the surface plasmons at the boundary of dielectric and metallic materials, an optical field can be confine into a subwavelength dimension [10]. A plasmonic circuit transforms the light into surface plasmons at the same frequency but reduced dimensions, in addition to guidance that is more flexible and highest speed.
In the next section of this chapter, Metal-Optics Electromagnetism is discussed. In section 3, wave analysis is used for explaining the propagation of surface plasmon polaritons. Section 4 is dedicated to investigate the various schemes of plasmonic guiding and their characteristics. In section 5, the modeling methods of plasmonic waveguides are introduced, and in section 6, a case study of plasmonic application in designing a switch is presented.

Metal-optics electromagnetism
Surface plasmons polaritons (SPPs) are surface electromagnetic modes that propagate at the interface of a dielectric with real electric permittivity ε 1 and a metal with permittivity ε 2 (ω)<0 as shown in figure 1.

Photonic Crystals 94
The complex permittivity or complex dielectric function can be expressed as [11]: ' '' i e e e = + (1) where ε ' is the real part and ε '' is th imaginary part of the permitivitty. In addition, the response of the material to the incoming optical field is expressed as the complex refractive index as [12]: where n and k are the real and imaginary part of the refractive index. The equations (1) and (2) are related to each other by: These parameters are known as optical constants of the material however, in many of the materials they change with the frequency of the incident optical field. Especially in metals, the dielectric parameters are strongly dependent to the optical frequency. In order to use the permittivity of metals in the calculations of the spectrometry some mathematical models are introduced, such as Lorentz model, Lorentz-Drude model andExtended Drude model.

Lorentz model
In plasmonic, the optical properties of material can be expressed as a classical physical model of the microscopic structure [3], in which the charge carriers are assumed as damped harmonic oscillators subjected to incident electromagnetic field as driving forces. In the Lorentz model a charge carrier is supposed to have a mass of m, charge e and displacement from equilibrium x. The force on this particle as a linear spring force is F = K .x, the velocity dependent damping is F = b.ẋ and the driving force by any incident light is F = eE . Therfore, the equation of motion is [12,13]: normalizing the above equation by mass: where ω 0 2 = K / m and γ = b / m. By solving the equation (6) as a time-harmonic and substitution of ẋ ↔ − iωx and ẍ ↔ − ω 2 x : Supposing N particles per volume V and dipole p = ex for each particle, the polarization is Therfore, equation (7) can be expressed: where ω p is the plasma frequency and is defined as: From the Maxwell's equation, we have: Therefore, from the equations (8), (10) and (11) the Lorentz model for the dielectric function is derived:

Multi-oscillator Drude-Lorentz model
Here, the model is the superposition of j+1 individual oscillators [12]: where f j ' s are the oscillator strengths, ω p, j ' s are the plasma frequencies, γ j ' s are damping rates and ω j ' s are the oscillator frequencies. The second term of equation (13) is the Drude term and represents the free electrons, which are not under the restoring force. Using 5 terms of the Lorentz model constructs an excellent approximation of optical properties of metals and many semiconductors in the visible spectrum [12].

Wave approach for the propagation of the surface plasmon polaritons
Surface plasmon polaritons are electromagnetic waves propagating at the interface between a dielectric and a metal evanescently bounded in the perpendicular direction [4]. These electromagnetic surface waves are created from the coupling of the optical fields to oscillations of the metal's electron plasma. Based on the dispersion relation and the spatial field distribution, the surface plasmons are described quantitatively.

Deriving the dispersion relation
By solving the Maxwell's equation at a single metal-insulator interface, the wave equation is produced, in which for the transverse magnetic (TM or p) modes is [14]: and for transverse electric (TE or s) modes is: where k 0 = ω / c 0 is the wave vector of the propagating wave in vacuum. Now, we should consider a simple flat boundary between an insulator (z>0) with a real positive permittivity ε 2 and a metal (z<0) with a complex permittivity ε 1 (ω) (for metals Re{ε 1 (ω) } < 0). Supposing the condition of propagating wave bounded to the interface with evanescent falloff in z-direction, the solutions for TM waves in z > 0 is [14]: and for z < 0 is: ( ) ( ) where k i ≡ k z,i is the component of the wave vector at z-direction. The reciprocal value of the k i , ẑ = 1 / | k z | , is evanescent decay length of the fields in the perpendicular of the propagation [14].
The continuity condition of H y and εE z at the interface implies that A 1 = A 2 and 2 2 Because at the interface of insulator and metal Re{ε 2 } > 0 and Re{ε 1 } < 0 the propagating waves confines at the surface. Based on equation (14) for TM waves at metal and insulator [14], Therefore, by combining the equations (22), (23) and (24), the dispersion relation of the SPPs propagating at the metal-insulator yields

Forbidden modes for surface plasmons
Based on equation (15), the field components for TE modes are In calculations of the equations (26) to (31), the following equations are used: From the continuity of E y and H x at the interface ( ) and because of the confinement condition requires Re{k 1 } > 0 and Re{k 2 } > 0, hence A 1 = 0 and A 2 = 0 ; therefore, there are no TE modes at the surface. That means, the surface plasmon polaritons propagate only in TM polarization [14].

Propagation Length of SPPs
An important parameter in comparing plasmonic waveguide structures is the propagation length of a supported mode. As a surface plasmon propagates along the surface, it quickly loses its energy to the metal due to absorption. The intensity of the surface plasmon decays with the square of the electric field, so at a distance x, the intensity has decreased by a factor of exp(2Im(β z )z) [15]. The propagation length is defined as the distance for the surface plasmon to decay by a factor of 1/e. This condition is satisfied at a length ( ) The typical value for propagation length in plasmonic metal-insulator waveguides is up to 100μm in visible regime [14].

Plasmonic waveguide schemes
The fundamental element in integrated plasmonic circuitry, as in photonics, is the waveguide. Because the necessity of existence of metal layer and various structural possibilities of using metal layer, many of schemes are developed in the recent years [11]. Each of the schemes has some benefits and drawbacks, in which make a trade-off between their characteristics. Two main characteristics of a waveguide are the confinement area of the mode and the propagation length. This features of the plasmonic waveguides are in opposite and there are a trade-off between them. In this section, some of the well-known types of plasmonic waveguides structures will be briefly introduced. In each case, the advantage(s) and weakness(s) of the waveguide are investigated.

Insulator-Metal (IM) waveguide
The simplest structure for guiding the plasmon polaritons is the Insulator-Metal (IM) waveguide ( figure 2). This scheme consists a thin film metal coated on a simple insulator strip waveguide [11]. The silicon core of the waveguide has dimension about 300nm×300nm and the cladding of the waveguide is the air. The metal cap has a 50nm thickness and can be gold or silver. The structure placed on a SiO 2 substrate. The propagation length of this fundamental scheme is 2μm at 1550nm [11]. In this structure, a relative good balance exists between propagation length and confinement. In addition, the material used in fabrication allows feasibility of integrated plasmonic circuitry [11].

Dielectric-loaded SPP waveguide
The dielectric-loaded SPP waveguides (DLSPPWs) are formed by placing a polymethylmethacrylate (PMMA) ridge with a cross section of 500 nm in width by 600nm in height as the dielectric material on top of a 65nm thick and 3μm wide gold strip [16] (figure 3a). This scheme provides a full confinement in the plane perpendicular to the propagation direction. This geometry is one of the most popular for plasmonic waveguides because PMMA can work as both resist and the dielectric core for the DLSPPW [16]. There is a convenient and simple process to fabricate plasmonic devices based on DLSPPW using deep-ultraviolet (UV) lithography for those devices with larger dimensions working at telecom wavelengths or standard E-beam lithography (EBL) for those working in the near-infrared. In addition, due to the physical dimension and therefore the mode size increase, this structure has slightly better propagation length for SPP's going from 5μm to 25μm.
Running Title 1 Fig. 2. Optical field distribution at cross section of insulator-metal waveguide 2 In this section, some of the well-known types of plasmonic waveguides structures will be 3 briefly introduced. In each case, the advantage(s) and weakness(s) of the waveguide are 4 investigated. 5

Insulator-Metal (IM) waveguide 6
The simplest structure for guiding the plasmon polaritons is the Insulator-Metal (IM) 7 waveguide (figure 2). This scheme consists a thin film metal coated on a simple insulator 8 strip waveguide [11]. The silicon core of the waveguide has dimension about 300nm×300nm 9 and the cladding of the waveguide is the air. The metal cap has a 50nm thickness and can be 10 gold or silver. The structure placed on a SiO2 substrate. The propagation length of this 11 fundamental scheme is 2m at 1550nm [11]. In this structure, a relative good balance exists 12 between propagation length and confinement. In addition, the material used in fabrication 13 allows feasibility of integrated plasmonic circuitry [11]. 14

Dielectric-loaded SPP waveguide 15
The dielectric-loaded SPP waveguides (DLSPPWs) are formed by placing a 16 polymethylmethacrylate (PMMA) ridge with a cross section of 500 nm in width by 600nm in 17 height as the dielectric material on top of a 65nm thick and 3μm wide gold strip [16] ( figure  18 3a). This scheme provides a full confinement in the plane perpendicular to the propagation 19 direction. This geometry is one of the most popular for plasmonic waveguides because 20 PMMA can work as both resist and the dielectric core for the DLSPPW [16]. There is a 21 convenient and simple process to fabricate plasmonic devices based on DLSPPW using 22 deep-ultraviolet (UV) lithography for those devices with larger dimensions working at 23 telecom wavelengths or standard E-beam lithography (EBL) for those working in the near-24 infrared. In addition, due to the physical dimension and therefore the mode size increase, 25 Compared with other waveguide schemes, DLSPPWs are well-matched with different dielectrics and have a rather good trade-off between mode confinement and propagation distance [16]. These characteristics make them suitable for realization of dynamic components by utilizing of material (e.g., thermo-and electro-optic) effects, while strong mode confinement and long propagation distances are required for realization of compact and complex plasmonic circuits [16].

Long-Range Dielectric-Loaded Surface Plasmon Polaritons waveguides (LR-DLSPPs)
The long-range plasmonic waveguide is, unlike previous studies of DLSPPWs, based on a metal film with finite width [17], i.e., a metal strip, and can thus be considered a hybrid plasmonic waveguide (discussed later). The LR-DLSPPW configuration consists of a dielectric ridge placed on top of a thin metal strip, which is supported by a dielectric film (figure 4a). The entire structure is supported by a low-index substrate that ensures mode confinement to the dielectric ridge and underlying dielectric film [17].
Photonic Crystals 8 this structure has slightly better propagation length for SPP's going from 5µm to 25µm. 1 Compared with other waveguide schemes, DLSPPWs are well-matched with different 4 dielectrics and have a rather good trade-off between mode confinement and propagation 5 distance [16]. These characteristics make them suitable for realization of dynamic 6 components by utilizing of material (e.g., thermo-and electro-optic) effects, while strong 7 mode confinement and long propagation distances are required for realization of compact 8 and complex plasmonic circuits [16]. The long-range plasmonic waveguide is, unlike previous studies of DLSPPWs, based on a 12 metal film with finite width [17], i.e., a metal strip, and can thus be considered a hybrid 13 plasmonic waveguide (discussed later). The LR-DLSPPW configuration consists of a 14 dielectric ridge placed on top of a thin metal strip, which is supported by a dielectric film 15 (figure 4a). The entire structure is supported by a low-index substrate that ensures mode 16 confinement to the dielectric ridge and underlying dielectric film [17]. 17 propagation length to L = 3100m, while retaining strong mode confinement [17]. 23

Hybrid Plasmonic Waveguide 24
As a solution to the issue of propagation loss of SPPs, a new waveguide known as hybrid 25 plasmonic waveguide (HPWG) is presented [18]. Figure 5a shows the cross section of a two 26 dimensional HPWG. It consists of a high index region (silicon) disjointed from a silver 27 surface by a low index layer (SiO2). The close area of the silver-silica interface and the silicon 28 slab results in coupling of the SPP mode and dielectric waveguide mode supported by these 29 two structures. Figure 5b shows the resulting hybrid mode. The guide also supports a 30 Because of varying the thickness of the PMMA layer and hence, balancing the mode field on either side of the metal strip (figure 4b), this structure minimizes the losses and increase the propagation length to L = 3100μm, while retaining strong mode confinement [17].

Hybrid plasmonic waveguide
As a solution to the issue of propagation loss of SPPs, a new waveguide known as hybrid plasmonic waveguide (HPWG) is presented [18]. Figure 5a shows the cross section of a two dimensional HPWG. It consists of a high index region (silicon) disjointed from a silver surface by a low index layer (SiO 2 ). The close area of the silver-silica interface and the silicon slab results in coupling of the SPP mode and dielectric waveguide mode supported by these two structures. Figure 5b shows the resulting hybrid mode. The guide also supports a conventional TE mode, which is shown in figure 5c. The mode sizes for both TE and TM modes are comparable in this case and are very similar to mode size achievable in case of silicon waveguide. HPWG offers a number of benefits [18]: it has an improved compromise between loss and confinement compared to purely plasmonic waveguides, and is compatible with silicon on insulator technology. The power of the TE and TM modes in the HPWG are concentrated in two different layers, therefore their guiding characteristics can be controlled in different manners by altering the material properties, and waveguide dimensions of the layers [18].

Metal-insulator-metal waveguide
In contrast to the hybrid waveguides, metal-insulator-metal (MIM) structure has a very good modal confinement in cost of a short propagation length. A MIM waveguide is a dielectric slot sandwiched as core between two layers of metal as cladding. The coupling of two SPP's from two metal-dielectric boundaries is the only aspect in common. In the MIM structure, these two SPP's coupled into the central dielectric slot and thus gives rise to a huge field concentration. However, due to the close proximity of the mode with both metal layers, the losses are extremely large. A typical value for propagation length is in the few microns as in the IM case. The most noticeable advantage of the MIM mode is its sub-wavelength size. In addition, in the MIM waveguides one could make any sharp geometries or bending regardless of wavelength of the incident optical field [20]. MIM subwavelength plasmonic waveguide bends and splitters have low loss over a wide frequency range [20]. While metals are naturally lossy, the bounded SP modes of a single insulator-metal interface can propagate over several microns under optical incident fields [21]. In such a geometry, the field skin depth increases exponentially with wavelength in the insulator but is almost constant (~25nm) in the metal for visible and near-infrared excitation frequencies. Not unlike conventional waveguides, these metal-insulator-metal (MIM) structures guide light via the refractive index differential between the core and cladding [21]. However, unlike dielectric slot waveguides, both plasmonic and conventional waveguiding modes can be accessed, depending on transverse core dimensions. MIM waveguides may thus allow optical mode volumes to be reduced to subwavelength scales-with minimal field decay out of the waveguide physical cross section-even for frequencies far from the plasmon resonance.

Dispersion relation for MIM waveguides
The modal analysis of planar multilayer structures can be solved via the vector wave equation under constraint of tangential E and normal D field continuity [21]. For unpolarized waves in a three-layer symmetric structure, the electromagnetic fields take the form [21]: , , with E y , B x and B z identically zero for transverse magnetic (TM) polarization and E x , E z and B y identically zero for transverse electric (TE) polarization. Inside the waveguide, the field components may be written as [21]: for the TM polarization and as: for the TE polarization. Outside the waveguide, the components are given [21]: for TM polarization and as: for TE polarization. The in-plane wave vector k x is defined by the dispersion relations [21]: where k z2 is the momentum conservation: and L + and L − denote the antisymmetric and symmetric tangential electric field with respect to the waveguide medium, respectively. In above analysis the structure is centered at z=0 with core thickness d and wave propagates along the positive x direction (figure 7). The core (cladding) is composed of material with complex dielectric constant ε 1 (ε 2 ). Since continuity of E y forbids charge accumulation at the interface, TE surface plasmon waves do not generally exist in planar insulator-metal structures [21]. If SPPs are excited at an insulator-metal interface, electrons in the metal make a surface polarization that leads to a localized electric field. In insulator-metal-insulator (IMI) structures, electrons of the metallic core screen the charge configuration at each interface and maintain a near-zero (or minimal) field within the waveguide. Therefore, the surface polarizations on both side of the metal film persist in phase and there is no cutoff frequency for any transverse waveguide dimension [21]. In contrast, screening does not occur within the insulator core of MIM waveguides. At each insulator-metal interface, surface polarizations of each side of interface are independent, and therefore SPPs oscillations need not be energy-or wave-vectormatched to each other. As a result, for certain MIM dielectric core thicknesses, interface SPPs may not remain in phase but will exhibit a beating frequency; as transverse core dimensions are increased, "bands" of allowed energies or wave vectors and "gaps" of forbidden energies will be observed [21].
The effective refractive index n eff of an MIM waveguide is complex. Its real part controls the guided wavelength λ MIM and its imaginary part limits the propagation length L SPP of SPPs through the relation [19]: where k = 2π / λ, λ = 2πc 0 / ω and c 0 is the speed of light in vacuum. One can calculate β using the following dispersion relation for the TM-SPP modes [19]: where the signs ± correspond to symmetric and antisymmetric modes. Based on Equation (46), for h ≪ λ, an MIM waveguide supports only a single antisymmetric mode that is similar to the fundamental TEM mode of a parallel-plate waveguide with perfect-electric-conductor (PEC) boundaries: The symmetric mode stops to exist because it experiences a cut-off in the reciprocal space [19]. Thus, in the deep subwavelength regime, an MIM waveguide operates as a single-mode plasmonic waveguide. Since this regime is the most interesting from the standpoint of nanophotonics applications, we assume in this chapter that the condition h ≪ λ is satisfied and focus on a single-mode plasmonic waveguide [19].

Challenges in modeling of plasmonics
Numerical simulation and modeling of plasmonic devices involves several challenges specific to plasmonics [22]. The permittivity of metals at optical wavelengths is complex, i.e. ε(ω) = ε ' (ω) + iε '' (ω), and is a complicated function of frequency [22]. Thus, several simulation techniques which are limited to lossless, non-dispersive materials are not applicable to plasmonic devices [22]. In addition, the dispersion properties of metals have to be approximated by suitable analytical expressions, such as Drude model and Drude-Lorentz model. Furthermore, in surface plasmons propagating along the interface of a metal and an insulator, the field is confined at the interface, and decays exponentially perpendicular to the direction of propagation. Consequently, for numerical methods based on discretization of the fields, a very fine mesh resolution is required at the metal-dielectric interface [22]. Generally, simulation of plasmonic devices necessitates much finer grid resolution than modeling of low-or high-index-contrast dielectric devices, because of the high localization of the field at metaldielectric interfaces of plasmonic components. The required grid size depends on the shape and feature size of the modeled plasmonic device, the metallic material used and the operating frequency.

General simulation methods for plasmonic components
Due to the wave nature of the surface plasmon polaritons, most of the simulation methods of plasmonic devices and circuits are the same as the radio frequency approaches; therefore, the most common methods in simulation and modeling of the plasmonics are Finite-Difference Time Domain (FDTD) [23], Finite Element Method (FEM) [24], Transmission Line Method (TLM) [19] and Coupled Mode Theory (CMT) [25]. Each of the methods has its advantages over a specific plasmonic structure, e.g. in order to simulation of a 3D hybrid structure of a ring resonator FDTD or FEM are more convenient and for modeling of a metal-insulator-metal stub filter the accurate and fast approach is TLM. Transient response of a structure may be modelled by simulating electromagnetic pulses in time domain, by methods such as FDTD.
The discussion for investigation of the best way to simulate plasmonic waveguides is out of the scope of this chapter; however, some of the most widely used methods are briefly explained.

Finite-Difference Time Domain (FDTD)
Finite Difference Time Domain is a common technique in modelling electromagnetics problems [15]. It is considered easy to understand and implement [26]. It is a time-domain method, so depending on the excitation type used, it could cover a wide range of frequencies in a single run, and hence it is usually the method of choice for wideband systems [15]. This method was introduced first by Kane Yee. FDTD does not require the existence of Green's function and directly approximates the differential operators in Maxwell's equations on a grid discretized in time and space. FDTD is an explicit finite difference approach, i.e. no matrix equation is established, stored and solved. The field values at the next time step are given entirely in terms of the field at the current and the previous time steps [15].
FDTD discretizes into unit cells called "Yee cells" [27,28]. In such Yee cells, Electric Fields are represented by the edges of a cube, where the faces of the cubic cell denote the magnetic field. Given the offset (in space) of the magnetic fields from the electric fields, the values of the field with respect to time are also offset. Time is distributed into small steps, which are corresponding to the amount of time needed for the fields to travel from one Yee cell to the next or less. FDTD method solves the Maxwell's equations using the relationship between the partial time and space derivatives [15]. Yee's algorithm solves both E and H in time and space using the Maxwell's curl equations [26]. The biggest advantage of FDTD is its simplicity of use and fast application.

Finite Element Method (FEM)
The finite element method (FEM) is a numerical technique, which approximates the solutions of differential equations [15]. In electromagnetics, FEM generally approaches the solution to the Maxwell's equations in the frequency domain, therefore is usually used with timeharmonic conditions. Furthermore, it is capable of time domain simulations. FEM can be derived using two methods [15]; First, the variational method, which finds a variational functional whose minimum/maximum/stationary point corresponds to the solution of the PDE subject to certain boundary conditions (a brief formulation is given in section 2.7.1). Second method is called the "weak formulation" in the literatures. It works by introducing a weighted residual error to one of the differentials in the PDE form of Maxwell's equations [29] and equating the sum of the error to zero. If the weighting functions are Dirac delta functions, the resulting procedure is similar to finite difference method. If the weighting functions are the basis functions, then the method is called the Galerkin's method [15].
Several commercial software are presented in which simulate the electromagnetic problems based on FEM, such as CST Microwave Studio, COMSOL Multiphysics and ANSYS.

Transmission Line Modelling (TLM) method
Transmission Line modelling is one of the time domain techniques that can solve the electromagnetic problems [15]. The concept of impedance and understanding the effects of waveguide discontinuities in terms of lumped circuit elements were crucial in this respect [22]. Although the properties of metals are quite different at optical wavelengths compared to the microwave, designs that are qualitatively similar to their low frequency counterparts have been demonstrated at optical frequencies.
In this method, a unit cell is formed by conceptually filling space with a network of transmission-lines in such a way that the voltage and current give information on the electric and magnetic fields [15]. A node represents the intersection point of the transmission-lines. At each time step, voltage pulses come to the node from each of the transmission-lines. These pulses are then scattered to produce a new set of pulses that become incident on adjacent nodes at the next time step. The relationship between the incident pulses and the scattered pulses is determined by the scattering matrix, which is set to be consistent with Maxwell's equations. Additional elements, such as transmission-line stubs, can be added to the node so that different material properties can be represented [15].

Stub filters in MIM plasmonic structures
A stub is one of the key elements in microwave engineering and is employed in various microwave devices to reduce their size [30]. Some research groups have proposed a wavelength filter by using a stub structure in a photonic crystal waveguide [30,31]. Such a structure may be employed in a plasmon waveguide to perform as a wavelength selective filter.
(a) 6. A case study: electro-plasmonic switch based on MIM stub filter 13

Stub filters in MIM plasmonic structures 14
A stub is one of the key elements in microwave engineering and is employed in various 15 microwave devices to reduce their size [30]. Some research groups have proposed a 16 wavelength filter by using a stub structure in a photonic crystal waveguide [30,31]. Such a 17 structure may be employed in a plasmon waveguide to perform as a wavelength selective 18 Figure 8. (a) three-dimension and (b) two-dimensional schematic depiction of the proposed Electro-Plasmonic switch [34].
Stub filters, also could be employed in order to make plasmonic switches [32][33][34]. MIM plasmonic waveguide structure helps to shrink the size of a stub filter to nanoscale and the stub structures cause to reduce the threshold voltage of the switch.

MIM Stub filter as a plasmonic switch
Figures 8(a) and 8(b) show the three and two-dimensional schematic of the proposed electroplasmonic switch based on the MIM structure, respectively [34]. This switch is essentially an MIM stub filter, which comprises of a main waveguide and one or more stub(s) vertically connected to it. Several studies investigated the properties of the stub filters and many researches used their filtering characteristics [19,[35][36].
The width of the waveguide and stubs, the length of the stubs and the displacement between them, and the refractive index of the core adjust the properties of the stub filter [19]. Any change in these factors will change the transfer function of the filter; therefore, one can control the transmittance valley or peak of the filter with a change in above parameters. Thus, through altering the refractive index of the core by applying an electric field to an EO material one could have an EO switch. In an MIM structure, the metal cladding boosts the electro optic effect in the core; subsequently, it decreases the threshold voltage of the switch.
Here, it is supposed that the entire of the switch filled with an EO material known as 4dimethyl-amino-Nmethyl-4-stilbazolium tosylate (DAST) as the core, with a linear refractive index of n d = 2.2 and a large EO coefficient (dn/dE = 3.41 nm/V) [37].
Due to the low losses for the surface plasmons propagation, here silver is selected as the metal cladding of the waveguide [38]. In all simulations, the Drude-Lorentz model of the silver is utilized in order to obtain accurate results. A seven-pole Drude-Lorentz model is used in the wavelength range from 0.4μm to 2μm [19]: where, ω p = 2002.6 THz is the bulk plasma frequency of silver and γ = 11.61 THz is a damping constant [19].
where, w is the width of the waveguide and ε d = n d 2 . The real and imaginary parts of the effective index for w =50nm over the wavelength range from 0.4 to 2μm is depicted in figure 9.  [19]. Figure 9. The real and imaginary part of the effective refractive index for MIM waveguide with the linear refractive index n d =2.2 and W wg = 50nm [34].

Modeling of the device using TLM
The Transmission Line Method (TLM) helps to study the operation of the MIM stub filters more accurately and faster [35], in comparison to other methods. In this method, we convert the plasmonic MIM structure in figure 10a to an equivalent circuit diagram (see figure 10b). This equivalent circuit is formed by a parallel connection of the characteristic impedance of an infinite MIM waveguides (Z MIM ) and the characteristic impedance of a finite and L spp is the propagation length of the SPPs (equation (35)). The length of the stubs is the key point of the design of the stub filter and the electro-plasmonic switch [34]. Knowing that small increase in the real part of the refractive index leads to a red-shift in the transmittance spectrum of the filter, the transmission of the filter at wavelength λ=1550nm must be on the edge of a change, a falling edge of the magnitude for normally ON and a rising edge for normally OFF switch. Figure 11 shows the transmission of a typical 4-stub filter for different length of the stubs. The transmittance is defined as T=P out /P in and is calculated using Transmission Line Method [35].
In the computation of the spectrum, the width of the main waveguide and the width of the stubs has the same value W wg = W stub =50nm, the distances between stubs is d=100nm and the thickness of the silver layer (or the depth of the structure) is t Ag =200nm.
MIM waveguide (Z S ) terminated by Z L accounts for reflection of the SPPs from the stub [19]: (50) β(w) is the complex-valued propagation constant and describes the properties of the MIM waveguide. The circuit diagram in figure 10b may be modified to an equivalent form as shown in figure 10c by replacing Z MIM and Z S by an effective impedance (Z stub ). The value of Z stub can be found from transmission-line theory and is given by [19]: Now, we calculate the transmission spectra of an MIM stub filter with N identical stubs set apart by a distance d, the length of the stubs L stub , the width of waveguide and stubs W wg = W stub = w [39]: where   The corresponding transmission-line representation. (c) The simplified circuit model. Z MIM and Z S correspond to the characteristic impedance of the MIM waveguide and the stub, respectively, Z L accounts for reflection of the SPPs from the stub end and Z stub is the effective stub impedance [19]. show an acceptable conformity [34]. 7 and Lspp is the propagation length of the SPPs (equation (35)). The length of the stubs is the 8 key point of the design of the stub filter and the electro-plasmonic switch [34]. Knowing that 9 small increase in the real part of the refractive index leads to a red-shift in the transmittance 10 spectrum of the filter, the transmission of the filter at wavelength λ=1550nm must be on the 11 edge of a change, a falling edge of the magnitude for normally ON and a rising edge for 12 normally OFF switch. Figure 11 shows the transmission of a typical 4-stub filter for different 13 length of the stubs. The transmittance is defined as T=Pout/Pin and is calculated using 14 Transmission Line Method [35]. In the computation of the spectrum, the width of the main 15 waveguide and the width of the stubs has the same value Wwg= Wstub=50nm, the distances 16 between stubs is d=100nm and the thickness of the silver layer (or the depth of the structure) 17 is tAg=200nm. 18 According to figure 11, by selecting the stub Lstub=300nm (point b) or Lstub=410nm (point a) 19 at the communication wavelength λ=1550nm, the normal state of the switch would be ON 20 or OFF, respectively. Slightly increment in the refractive index turns the switch OFF or ON, 21 According to figure 11, by selecting the stub L stub =300nm (point b) or L stub =410nm (point a) at the communication wavelength λ=1550nm, the normal state of the switch would be ON or OFF, respectively. Slightly increment in the refractive index turns the switch OFF or ON, depending on the length of stubs (L stub =300nm or L stub =410nm, respectively). Figures 12a and  12b show the spectrum of the switch versus the refractive index for L stub =300nm and L stub =410nm, respectively. Figures 13a and 13b show the transmission spectrum for L stub = 410nm and L stub = 300nm and shows reasonable accuracy of three-dimensional FEM simulations [34]. These figures admit the red-shift effect on the ON/OFF status of the switch. As shown in figure 8, the silver cladding is connected to the voltage V and therefore the electric field surrounds the main waveguide and all the stubs. The relation between applied voltages and the refractive index of EO material is expressed using equation [37]: 0 ( ) dn V n n dE w = + where n 0 is the linear refractive index of the EO material that equals 2.2, V is the applied voltage, and w is the width of waveguides and stubs [34]. The applied electric field alters the linear refractive index of the switch and as a result causes the red-shift. Figures 14a-d demonstrate the distribution of the optical field in the structure of the electro-plasmonic switch. In figures 14a and 14c the applied voltage is 0V and corresponds to the normal-state of both switches, while figures 14b and 14d show the electro-optic switch when the voltage is applied to the structure. The applied electric field changes the refractive index of the EO material; therefore, the optical field reflected back to the input port. Figure 15 depicts the transmittance of the normally ON and normally OFF switch as a function of the applied voltage at the wavelength 1550nm. As can be seen, the threshold voltage for switching is 10V [34].
Transmission Line Method [35]. In the computation of the spectrum, the width of the main 15 waveguide and the width of the stubs has the same value Wwg= Wstub=50nm, the distances 16 between stubs is d=100nm and the thickness of the silver layer (or the depth of the structure) 17 is tAg=200nm. 18 According to figure 11, by selecting the stub Lstub=300nm (point b) or Lstub=410nm (point a) 19 at the communication wavelength λ=1550nm, the normal state of the switch would be ON 20 or OFF, respectively. Slightly increment in the refractive index turns the switch OFF or ON, 21