Improved Plasmonic Filter, UltraCompact Demultiplexer, and Splitter
 Author: Rahimzadegan Aso, Granpayeh Nosrat, Hosseini Seyyed Poorya
 Organization: Rahimzadegan Aso; Granpayeh Nosrat; Hosseini Seyyed Poorya
 Publish: Journal of the Optical Society of Korea Volume 18, Issue3, p261~273, 25 June 2014

ABSTRACT
In this paper, metal insulator metal (MIM) plasmonic slot cavity narrow bandpass filters (NBPFs) are studied. The metal and dielectric of the structures are silver (Ag) and air, respectively. To improve the quality factor and attenuation range, two novel NBPFs based on tapered structures and double cavity systems are proposed and numerically analyzed by using the twodimensional (2D) finite difference time domain (FDTD) method. The impact of different parameters on the transmission spectrum is scrutinized. We have shown that increasing the cavities’ lengths increases the resonance wavelength in a linear relationship, and also increases the quality factor, and simultaneously the attenuation of the wave transmitted through the cavities. Furthermore, increasing the slope of tapers of the input and output waveguides decreases attenuation of the wave transmitted through the waveguide, but simultaneously decreases the quality factor, hence there should be a tradeoff between loss and quality factor. However, the idea of adding tapers to the waveguides’ discontinuities of the simple structure helps us to improve the device total performance, such as quality factor for the single cavity and attenuation range for the double cavity. According to the proposed NBPFs, two, three, and fourport power splitters functioning at 1320 nm and novel ultracompact twowavelength and triplewavelength demultiplexers in the range of 13001550 nm are proposed and the impacts of different parameters on their performances are numerically investigated. The idea of using tapered waveguides at the structure discontinuities facilitates the design of ultracompact demultiplexers and splitters.

KEYWORD
Metalinsulatormetal (MIM) waveguide , Optical multiplexer , Optical splitter , Optical plasmonic filter , Nanodisk resonator

I. INTRODUCTION
Modern electronic devices are rapidly approaching their fundamental speed and bandwidth limitations, creating a serious impediment for new applications and requirements. Carrying information by light and replacing electronic signals by lightwave is an auspicious solution towards higher speeds. The diffraction limit prevents the localization of light in areas smaller than a wavelength; therefore, miniaturization and integration of conventional photonic circuits are impossible [1]. Excitation of surface plasmons is one of the most feasible ways to guide electromagnetic waves beyond the diffraction limit to control light in the nanometer scales. Estimations show that data rates of 10 Tbit/s are achievable by decreasing the photonic devices to subwavelength scales [2].
A plasmon is a quasiparticle and quantum of plasma oscillation of free electrons in a metal. Plasmons with collective electron oscillations confined to the interface of metal and dielectric are called surface plasmons (SPs). When a photon (a quantum of electromagnetic wave) is coupled to a SP, the quasiparticle known as the surface plasmon polariton (SPP) is produced. SPPs are surface electromagnetic waves propagating along the metaldielectric interface and are evanescently confined in the perpendicular direction of the interface. In a cylindrical metal nanowire, when the diameter is reduced below the light wavelength, a strong monotonic increase in localization of light takes place and shows that the diffraction limit can be overcome [2].
Different metallic plasmonic nanoguiding structures including thin metal films, metal nanorods, metal wedges, nanogrooves, slot waveguides, etc. have been proposed and their properties have been studied [1]. Two of the wellknown multilayer plasmonic waveguides are metalinsulatormetal (MIM) and insulatormetalinsulator (IMI) structures. MIM structure gives more light confinement, but more propagation loss (lower propagation length). Therefore, there should be a tradeoff between propagation length and light confinement. It is shown that MIM structure is a superior choice [3]. It is also experimentally verified that silver, due to its lower guiding loss, is an appropriate metal for plasmonic waveguides [4].
To realize photonic circuitry based on plasmonics, a variety of components such as splitters [57], couplers [8, 9], multidemultiplexers [1013], switches [14], logic gates [15], etc. have recently been studied. Filters as a fundamental block, have been at the center of attention. Filters can be divided into three categories of narrow bandstop filters (NBSFs) including fiber Bragg gratings, FabryPerot and slot filters [2, 8, 10, 1624], narrow bandpass filters (NBPFs) [10, 17, 11, 20, 2530], highpass [31], and lowpass (cutoff) [32] filters, which have been widely studied.
One of the most important characteristics of these filters (NBPF or NBSF), is their quality factor, defined as
Q =λ_{r} /Δλ , whereλ_{r} is the resonance wavelength of the cavity and Δλ is the full width at half maximum (FWHM) of the transmission spectrum. Without taking the loss of SPPs in metal into account, the quality factor is naturally high, the subatomic collisions are the main cause of quality factor reduction. For example, a ring NBSF is presented in which the quality factor without considering loss is as high as 4000, whereas when loss is taken into account, Q is reduced to 30 [8]. Higher values of Q are normally achieved as a result of lower coupling between the waveguide and the cavity resonator, hence for a desirable Q for NBSFs the dip atλ_{r} is not normally damped enough, and for NBPFs the spectrum peak value (atλ_{r} ) is undesirably low. For example, a novel bandstop whisperinggallery resonance cavity is presented with quality factor of up to 1400, but damping for the transmitted wave at the resonance frequencies is 5dB [16], which is not desirable compared to the other NBSFs, with higher damping [17]. The quality factor for a ring NBSF is as high as 166, but the corresponding value of damping level atλ_{r} is not given for comparison in [9]. A slot NBSF is proposed and a quality factor of 50 is achieved [17]. The Q of a coupler type MIM filter is as low as 25 [19]. Most of the proposed structures lack a good quality factor. Huanget al . proposed a toothshaped structure for a NBPF, but the quality factor was as low as 19 [20]. Using a novel NBPF structure [33], slot cavity [30], and plasmonic slot [26] result in the quality factors of 40, 153, and 37, respectively.Using filters as a fundamental block, different splitters such as one in [26], based on slot cavity filters and several plasmonic demultiplexer structures have been proposed. A plasmonic demultiplexer based on a Y bent waveguide is presented, but it needs a minimum space of 2μm between the input and the outputs, which may not be desirable for compact circuitry [10]. In [12], a nanoplasmonic triplewavelength demultiplexer is presented, but due to the weak coupling effect, the resonance frequencies are not transmitted properly and corresponding spectrum peak values are undesirably low. Tao
et al . proposed a novel demultiplexer based on nanocapillary resonators with spectrum peak value of 10 dB, which again is not a desirable transmission spectrum [34].In this paper, we begin with an inline cavity NBPF structure and based on tapered waveguides, we modify it to yield two novel slot cavity NBPF structures with two different perspectives. Both have an acceptable spectrum peak value. We are trying to obtain a relatively high quality factor filter with one of them and a filter with a high attenuation range with the other one. High attenuation range helps to damp the undesired frequencies around the resonance. Then, based on the proposed filters, a novel splitter and compact demultiplexer are proposed. The demultiplexer is around 900nm long and 350nm wide. The structures are numerically studied.
The paper is organized as follows; in Section II some introductory concepts of plasmons and the FDTD method are discussed. In Section III, the structures proposed are simulated. The paper is concluded in Section IV.
II. SURFACE PLASMON POLARITONS AND METHODS FOR THEIR ANALYSIS
The surface plasmon is a classic electromagnetic topic, known from 100 hundred years ago, but, due to the recent advances in nanotechnology, it has become an interesting topic of research for scientists [35]. At optical frequencies, metals can be described by a gas of free electrons moving against a positive background of ion cores, which are known as the plasma model. In response to an applied electromagnetic field, the electrons oscillate and their motion is damped with a characteristic collision frequency of
γ . The permittivity function resulting from this simple explanation, known as Drude model, is [36]:where
ω_{p} is the plasma angular frequency of the free electron gas, andω is the angular frequency of the lightwave. We have assumed the driving field as 𝐄(t )=𝐄_{0} exp(−iωt ) [36]. However, all the electrons in a metal are not free, interband transitions of electrons change the dielectric function [37]. Therefore, Eq. (1) is modified to the DrudeLorentz model [38]:where
ε (∞) is the highfrequency (electronic transition) contribution to the dielectric constant,ω _{i} denotes the oscillation frequency of the ith bound electron under an applied electric potential,S _{i} andγ _{i} are related to the density and damping of those bound electrons, respectively [38].In this paper, we have used the silver parameters of:
ε (∞) = 3.7,ω _{p} = 1.38×10^{16}rad/s andγ = 0.27×10^{14}rad/s [10, 12, 28, 39], which are extracted from experimental results [40]. Other parameters are given in Table 1.Solving the Helmholtz equation and applying boundary conditions, show that SPPs exist only for TM modes. This is due to the negative real part of the dielectric function of metal at optical frequencies [36]. However, lightwaves cannot directly excite the SPPs in normal conditions, because the photons’ momentum is lower than that of the SPPs. Some exciting techniques such as prism coupling, scattering from a topological defect, and using periodic corrugations in the metal’s surface are required to match the momentum. Once light is coupled to SPPs, it will be propagated along the interface and will also be attenuated accordingly [35].
The schematic view of a typical MIM waveguide, with an inline slot cavity is shown in Fig. 1. The incident wave is either reflected from or coupled to the cavity. Constructive reflections inside the slit at a specific wavelength form standing waves and the structure operates as a resonance device. Therefore, the incident wave is passed or blocked and the total structure operates as a NBPF. To find the empirical relation for the resonance wavelength, the condition for constructive interference is required. Based on the principle of a resonant cavity, the necessary condition for the case is that the phase delay per round trip in the cavity to be an even multiple of
π , Δϕ =2mπ , wherem is an integer. Then, the phase delay inside the cavity is Δϕ =k (ω )2L+2ϕ_{ref} wherek (ω ) is the wavenumber in the slit at frequency ofω , L is the cavity length, andϕ_{ref} is the phase shift of the wave reflected from the insulatormetal interface at each side of the cavity, which is low and can usually be neglected [19]. Considering , the resonance wavelength is derived as follows:where
n_{eff} is the real part of the effective refractive index of the resonant cavity [19]. Thus the resonance wavelength for the fundamental mode (m =1) is approximately:The effective refractive index is a function of waveguide width and the operating frequency. Using International Mathematics and Statistics Library (IMSL) Fortran subroutine,
n_{eff} can be calculated for different waveguide width for the wavelength range of 400 to 2000 nm. For d=100 nm,n_{eff} is calculated to be around [1.211.23] in the 5002000nm wavelength range [25].To simulate the plasmonic filters, we have used the FDTD method in which Maxwell’s equations are solved by discretizing time and space. The spatial grid size is Δ
z = Δx = 4nm, which is found to be sufficient for numerical convergence. For the two dimensional simulation box, the convolutional perfectly matched layers (CPMLs) are used in the boundaries to virtually simulate an infinite space for the propagating wave and to prevent reflections of outgoing waves. Due to the Courant condition, the time integration step is Δt = (0.6Δx/c )s .The numerical calculations are convergent for the number of time steps of up to 2^{20}. A very important point to be noticed is the distance of the source from CPML layers, which should be far enough to converge the numerical calculations. We have found that the distance of 1μm is required to avoid instability of the numerical simulations.
III. ANALYSIS AND SIMULATION
3.1. Narrow Bandpass Filter Structure
The first structure to be analytically and numerically investigated is depicted in Fig. 1(a). The inline resonator cavity with length L and width W operates as a narrow bandpass filter, and it only passes the cavity resonance wavelength. We are interested in the impact of the parameters of the cavity on the transmission spectra. We have mostly focused on the wavelengths longer than 1 μm. Though other modes may appear, we have concentrated on the fundamental mode to be appropriate for the second and third telecommunication windows. For more accuracy, we have calculated the input power without metallic gap to eliminate the cavity reflection effects. Therefore, for each graph we have run the codes two times, with and without the metallic gap and resonator cavity. The results are depicted in graphs of Fig. 1 and Table 2. The important data, quality factor, spectrum peak value and resonance wavelength in Table 2 are derived from the simulated bandpass filters for easier quantitative evaluations. The structure of Fig. 1(a) has been studied in some other references [10], [17], but the wavelength spectrum and dielectric medium are not the same, so for an exact analysis and comparison of the presented modified structures with the NBPF structure, we have chosen air as the dielectric and wavelengths above 1 μm as our simulating spectrum. Three important results should be noticed:
a) Increasing the cavity length increases the resonance wavelength in a linear relationship, which complies with the analytical results for resonance wavelength of Eq. (4). As shown in Table 2, the numerical and analytical results for λr as a function of L are approximately the same. It is also observed that increasing L, would increase Q, but decreases the spectrum peak value. b) Increasing the width of the cavity, W, enhances the filter quality factor, but decreases the spectrum peak value. Furthermore, the resonance wavelength is slightly decreased. The change in λr for the same length is the effect of change in ϕref and effective index which slightly changes the resonance conditions. c) Increasing the metallic gap size, G, enhances the filter Q and decreases the spectrum peak value, which is expected, because the thicker the metal is, the more wavelength selective the slit would be and the lightwave passes through it with more loss. If we let G = 48 nm, the quality factor may increase to as high as 245, but the spectrum peak value is 20 dB, which is too low for a filter. It is also observed that the variation of G, slightly changes λr, which is the effect of variation of the phase shift of the wave in passing through or reflecting from the cavity metal ends. Our results for structure of Fig. 1(a) comply well with those in [10] and [17].
It can be seen from the results that as far as the cavity is wider than the waveguide, the wave does not monotonically transmit to the cavity and part of it is reflected and blocked. One way to compensate the mismatch is the novel structure presented in Fig. 2(a), in which the cavity input and output waveguides are tapered. Impacts of each parameter on the transmission spectra are illustrated in graphs of Fig. 3 and some numerical data including
λ_{r} , spectrum peak value and Q are derived and given in Table 3. Again,λ_{r} is linearly related to the cavity length. Figure 2(b) compares the analytical and simulated results, which are in good agreement. Figures 2(c) and (d), show the effect of variation of metallic gap size on the quality factor and spectrum peak value, respectively. Comparison of the structures of Figs. 1 and 2 and Tables 2 and 3, are compared in Fig. 3, which shows that the novel proposed filter gives higher quality factor with lower loss, which represents an improvement compared to the simple filter. As an example, for the simple structure the quality factor of 97 is achievable with spectrum peak value of 5.7 dB (with G=24 nm, W=180 nm, L=564 nm) but in the tapered structure, higher quality factor of 113 and lower spectrum peak value of 5 dB (with G=28 nm, W=180 nm, L=564 nm, S_{i}=S_{o}=0.3) are achievable. Therefore, the single cavity NBPF with tapered waveguides shows better filter specifications. In Fig. 2(e), the impact of the taper slope, S=S_{i}=S_{o}, on the transmission spectra is investigated. Increasing the slope would decrease the loss, but on the other hand decreases the quality factor, and it also adds more difficulty in the fabrication process. Resonance wavelengths are relatively constant, so are not shown here.Another procedure for improving the specifications of the simple filter structure is to add another cavity along the waveguide, so that the filter could increase its selectivity. Figure 4(a) illustrates the proposed double cavity NBPF structure. It is important that the size of the metal between the two cavities, Gb, has to be thicker than a specific value, because if it is thin, the transmission spectrum shows two nearresonance wavelengths, which is the effect of coupling between the two cavities. The two cavities should be exactly the same to guarantee the same resonance wavelength. Figure 4(b), indicates the impact of the variation of the cavities’ length, L, on the transmission spectra. It can be noticed that as the length is increased, the quality factor and the resonance wavelength increase. The quality factor is improved, but the spectrum peak value is too low and the structure does not give desirable frequency response. In other words, when another cavity is added, although the selectivity is increased, the spectrum peak value is seriously decreased. To compensate again, we taper input and output waveguides of the new double cavity NBPF, as shown in Fig. 5(a). The impact of the variations of the cavities’ length and slope of the tapers on the transmission spectra are demonstrated in Fig. 5. Figures 5(bc), show the linear relation of the resonance wavelength to L, for this structure, which resembles the analytical results with a good approximation, and it seems that adding another cavity does not noticeably change the empirical relation on resonance wavelengths of the cavities. Figure 5(d), similarly to the tapered single cavity filter, shows an approximate linear relation of spectrum peak value and quality factor to L. Table 4 gives the quantitative data on the effect of cavity length on the quality factor, resonance wavelength and spectrum peak value. In Fig. 5(e), the effect of variation of taper slopes (S_{i}=S_{o}=S) on the quality factor and spectrum peak value is depicted.
One of the specifications of a bandpass filter is its attenuation range (A_{stop}), which shows the ability of the filter for damping the nonresonance wavelengths. It is shown in Fig. 5(f) as A_{stop}. A good characteristic of the structure is its high attenuation range (approximately 40 dB), which can seriously damp the undesired wavelengths.
Four different structures’ spectra are compared in Fig. 6. There are significant improvements in the tapered NBPF structures. In the tapered single cavity structure, with the same loss in the transmission spectrum, the quality factor is improved and in the tapered double cavity NBPF, with the cost of minor decreases in spectrum peak value, the attenuation range significantly increases.
3.2. Ultracompact Plasmonic Power Splitter
An ultracompact splitter structure is presented based on the novel proposed single cavity NBPF structure, as shown in Fig. 7(a). The numerical result of the designed power splitter for the wavelength of 1310nm is depicted in Fig. 7(b). In Fig. 7(c), the impacts of the input taper slope, S_{i}, on power splitter transmission are demonstrated. Increasing the taper slope could result in transmission with lower loss, but decreasing the quality factor and it also make the fabrication process more difficult, so there should be a tradeoff between these parameters. The structures for plasmonic triplewavelength and quartetwavelength and the results of simulations to show the power splitting performance of the proposed devices are illustrated in Figs. 8 and 9, respectively. The output ports’ widths, d_{o}, can be designed to be narrower; so that more port power splitters can easily be achieved. The effective parameters for design of the structures are given in the figures’ captions. Power splitting with different ratios is also achievable by changing the distance of the output port from the cavity.
3.3. Ultracompact Plasmonic Demultiplexer
Based on the investigated structures, a novel demultiplexer system is presented in Fig. 10(a). The lengths of the cavities are adjusted so that the resonance wavelengths are 1310 nm and 1550 nm. Design parameters are given in the figure caption. Magnetic field distributions for two wavelengths are also demonstrated in Figs. 10(b) and (c) to show the demultiplexing performance of the device. The impacts of the input taper slope, S_{i}, on the transmission spectrum are shown in Figs. 10(de). Steeper slopes cause better transmission with lower loss, but decrease the quality factor. Also, it would make the fabrication more difficult, so there should be a tradeoff. CD is the distance between the cavities, and it should be higher than a certain value to prevent coupling between cavities which would deteriorate the demultiplexing performance. The effect of CD on quality factor and spectrum peak value is investigated in Fig 10(f). We are interested in achieving a dense wavelength division demultiplexing. In Fig. 10(g), this possibility is tested, making one cavity length, L2, constant and slightly decreasing the other cavity length, L1, the resonance wavelength is collected and simultaneously the spectrum peak value is under supervision so that the loss is not a serious practical problem. In Fig. 10(h), we could successfully optimize and simulate a twowavelength demultiplexer with resonance wavelengths of 1396 and 1417 nm, which is very desirable for dense plasmonic wavelength division multiplexing. Without taking footprints of the plasmonic filter, the demultiplexer length, DL is as small as 900 nm, which is a compact structure compared to other conventional demultiplexers, making it very desirable for photonic integrated circuits. In Tables 5 and 6, the effects of cavities’ width, W, and output ports’ width, d_{o}, are investigated. It is seen that changing d_{o}, does not change considerably the parameters of the achieved filter.
The structure can be designed for triplewavelength demultiplexing, as shown in Fig. 11. Again, resonance wavelengths are adjusted by variation of the cavities' lengths. We have designed the structure for the wavelengths of 1310, 1430 and 1550 nm. Magnetic field distributions for the 3 resonance wavelengths are illustrated in Figs. 11(ce) to show the demultiplexing performance of the ultracompact three channel demultiplexer. Owing to the flexibility for variation of the cavities’ lengths, other desired demultiplexing wavelengths, with even more channels which are distant enough, are achievable.
IV. CONCLUSION
Two novel plasmonic slot cavity narrow bandpass filter (NBPF) structures have been proposed. The first one is a single cavity structure with tapered input and output waveguides, which increases the filter quality factor, and the second one consists of two successive cavities with the tapered waveguides, which increases the attenuation range of the filter. The devices are numerically analyzed and simulated by the two dimensional finite difference time domain (FDTD) method. Resonance wavelength and bandwidth can be adjusted by variation of the device parameters. Based on these filters, the novel plasmonic structures for two, three and four port splitter and two and three channel demultiplexers have been proposed, analyzed and simulated to demonstrate their performances. Different designs for power splitting and demultiplexing are achievable by adjusting the cavities' parameters with total length of less than 900 nm. Also, power splitters with different power ratios can be designed by varying the distance of the ports from the resonance cavities. Our results imply potential applications of the proposed structures in the nanoscale photonic integrated circuits on the flat metallic surfaces.

[]

[]

[TABLE 1.] Parameters of the DrudeLorentz model for silver [25]

[FIG. 1.] (a) Schematic view of a single inline slot cavity narrow bandpass filter (NBPF) structure. (b) Comparison of analytical and numerical resonance wavelengths. (c) Quality factor and spectrum peak values vs. the cavity length L, assuming di=do=100 nm, W=180 nm, and G= 16 nm. (d) Variations of the peak of the transmission spectra for three different values of the metallic gap sizes, G, assuming L=564 nm, d=100 nm, and W=180 nm. (e) Quality factor and spectra peak values vs. metallic gap size, G, assuming L=564 nm, di=do=100 nm, and W=180 nm.

[]

[]

[TABLE 2.] Simulated specifications of the single cavity NBPF of FIG. 1

[FIG. 2.] (a) Schematic view of the single cavity NBPF with tapered waveguides. Specifications of the single cavity NBPF with tapered waveguides. (b) Comparison of analytical and numerical resonance wavelengths. (c) Quality factor and spectra peak values vs. cavity length, L, assuming di = do = 100 nm, W = 180 nm, G = 16 nm, and Si = So = 0.3. (d) Variations of the peak of the transmission spectra for four different values of the metallic gap sizes, G, assuming L = 524 nm, d = 100 nm, W = 180 nm and Si = So = 0.3. (e) Quality factor and spectra peak values vs. metallic gap size, G, assuming L = 524 nm, di = do = 100 nm, W = 180 nm and Si = So = 0.3. (f) Quality factor and spectra peak values vs. taper slope, S, assuming L = 524 nm, d = 100 nm, W = 180 nm, G = 16 nm.

[FIG. 3.] Comparison of single cavity NBPF and tapered single cavity NBPF regarding quality factor and spectrum peak value. The straight lines are linear approximations of the data given.

[TABLE 3.] Simulated specifications of the single cavity NBPF with tapered waveguides of FIG. 2

[FIG. 4.] (a) Schematic view of a double cavity NBPF, (b) its transmission spectra for three different values of the cavity length L, with di = do = 100 nm, W = 180 nm, G = 16 nm, and Gb = 44 nm.

[FIG. 5.] (a) Schematic view of the double cavity NBPF with tapered waveguides. Specifications of the double cavity NBPF with tapered waveguides. (b) Variations of the peak of the transmission spectra for three different values of the cavity length L, assuming di = do = 100 nm, W = 180 nm, G = 16 nm, Gb = 44 nm, and Si = So = 0.3. (c) Comparison of analytical and numerical resonance wavelengths. (d) Quality factor and spectra peak values vs. cavity length L, assuming di = do = 100 nm, W = 180 nm, G = 16 nm, Gb = 44 nm, and Si = So = 0.3. (e) Quality factor and spectra peak values vs. the slope of the taper S, assuming di =do = 100 nm, L = 564 nm, W = 180 nm, G = 16 nm, and Gb = 44 nm. (f) Filter specifications, adapted from Matlab Simulink.

[TABLE 4.] Simulated specifications of double cavity NBPF with tapered waveguides

[FIG. 6.] Comparison of the frequency responses of single cavity NBPF (bluedashed line), tapered single cavity NBPF (reddotted line), double cavity NBPF (blackstripped line), and tapered double cavity NBPF (green line).

[FIG. 7.] Schematic view of the proposed two port power splitter structure, (b) its transmission spectra with L = 556 nm, W = 356 nm, G = 12 nm, di = do = 100 nm, Si = 2, and So = 0.5, and (c) transmission spectra for two different slopes of Sm. Magnetic field distributions for wavelengths of (c) 1320nm, demonstrate the splitting performance of the proposed device.

[FIG. 8.] (a) Schematic view of the proposed three port power splitter structure, and (b) its transmission spectra for L = 556 nm, W = 356 nm, G = 16 nm, di = 100 nm, do =56 nm, Si = 1, and So = 0.5. Magnetic field distributions for wavelengths of (c) 1320 nm, demonstrate the splitting performance of the proposed device.

[FIG. 9.] (a) Schematic view of the proposed four port power splitter structure, and (b) its transmission spectra for L = 556 nm, W = 356 nm, G = 16 nm, di = 100 nm, do = 48 nm, Si = 1, and So = 0.5. Magnetic field distributions for wavelengths of (c) 1320 nm, demonstrate the splitting performance of the proposed device.

[FIG. 10.] (a) Schematic view of the proposed demultiplexer structure. Magnetic field distributions for wavelengths of (b) 1550 nm and (c) 1310 nm, demonstrate the demultiplexing performance of the proposed device. (d) its transmission spectra for two values of Si, for L1 = 628 nm, L2 = 524 nm, W1 = W2 = 164 nm, G = 16 nm, So=1, CD = 88 nm, and DL = 0.88 um. (e) Quality factor and spectrum peak value (SPV) of both channels, vs. input taper slope, Si, assuming waveguide width, di = do = 100 nm, L1 = 628 nm, L2 = 524 nm, cavity widths of W1 = W2 = 164 nm, metallic gap size, G = 16 nm, CD= 88 nm, So = 0.5. (f) Quality factor and spectrum peak value (SPV) of both channels, versus the distance between two cavity, CD, assuming waveguide width, di = do = 100 nm, L1 = 628 nm, L2 = 524 nm, cavity widths of W1 = W2 = 164 nm, metallic gap size, G = 16 nm, Si = 0.5 and So = 1. (g) Effects of L1 on quality factor and spectrum peak value (SPV) of both channels, assuming L2= 524 nm, waveguide width, di = do = 100 nm, cavity widths of W1 = W2 = 164 nm, metallic gap size, G = 16 nm, CD = 88 nm, Si = 0.5, So = 1. (h) Optimized two wavelength multiplexer for plasmonic dense wavelength division demultiplexing, the structure parameters are: W = 96 nm, di = do = 100 nm, L2 = 524 nm, L1 = 532 nm, G=16 nm, Si = 0.5, So = 1, CD = 120 nm.

[TABLE 5.] Effect of changing output ports’ width, do, on filter specifications of twoport demultiplexer. Assuming: di = 100 nm, L1 = 628 nm, L2 = 524 nm, W1 = W2 = 164 nm, G = 16 nm, Si = 0.5 and So = 1, CD=88 nm

[TABLE 6.] Effect of cavities’ width, W=W1=W2, on filter specifications of twoport demultiplexer. Assuming: di = do = 100 nm, L1 = 628 nm, L2 = 524 nm, G = 16 nm, Si = 0.5 and So = 1, CD = 88 nm

[FIG. 11.] (a) Schematic view of the proposed triplewavelength demultiplexer structure, (b) its transmission spectra for L1 = 508 nm, L2 = 556 nm, L3 = 460 nm, W1 = W2 = W3 = 80 nm, G = 12 nm, di = 48 nm, WD = 80 nm, Si = 1, and So = 0.3. Magnetic field distribution for wavelengths of (c) 1430 nm (d) 1550 nm, and (e) 1310 nm, demonstrate the demultiplexing performance of the device.