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.
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
[5

7]
, couplers
[8
,
9]
, multidemultiplexers
[10

13]
, 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
,
16

24]
, narrow bandpass filters (NBPFs)
[10
,
17
,
11
,
20
,
25

30]
, 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. Huang
et 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
.
Parameters of the DrudeLorentz model for silver[25]
Parameters of the DrudeLorentz model for silver [25]
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
π
, Δ
ϕ
=2
mπ
, where
m
is an integer. Then, the phase delay inside the cavity is Δ
ϕ
=
k
(
ω
)2L+2
ϕ_{ref}
where
k
(
ω
) 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:
(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 d_{i}=d_{o}=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, d_{i}=d_{o}=100 nm, and W=180 nm.
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 inTable 2, the numerical and analytical results forλras 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λrfor the same length is the effect of change inϕrefand 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 ofFig. 1(a)comply well with those in[10]and[17].
Simulated specifications of the single cavity NBPF ofFIG. 1
Simulated specifications of the single cavity NBPF of FIG. 1
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.
(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 d_{i} = d_{o} = 100 nm, W = 180 nm, G = 16 nm, and S_{i} = S_{o} = 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 S_{i} = S_{o} = 0.3. (e) Quality factor and spectra peak values vs. metallic gap size, G, assuming L = 524 nm, d_{i} = d_{o} = 100 nm, W = 180 nm and S_{i} = S_{o} = 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.
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.
Simulated specifications of the single cavity NBPF with tapered waveguides ofFIG. 2
Simulated specifications of the single cavity NBPF with tapered waveguides of FIG. 2
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.
(a) Schematic view of a double cavity NBPF, (b) its transmission spectra for three different values of the cavity length L, with d_{i} = d_{o} = 100 nm, W = 180 nm, G = 16 nm, and Gb = 44 nm.
(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 d_{i} = d_{o} = 100 nm, W = 180 nm, G = 16 nm, Gb = 44 nm, and S_{i} = S_{o} = 0.3. (c) Comparison of analytical and numerical resonance wavelengths. (d) Quality factor and spectra peak values vs. cavity length L, assuming d_{i} = d_{o} = 100 nm, W = 180 nm, G = 16 nm, Gb = 44 nm, and S_{i} = S_{o} = 0.3. (e) Quality factor and spectra peak values vs. the slope of the taper S, assuming d_{i} =d_{o} = 100 nm, L = 564 nm, W = 180 nm, G = 16 nm, and Gb = 44 nm. (f) Filter specifications, adapted from Matlab Simulink.
Simulated specifications of double cavity NBPF with tapered waveguides
Simulated specifications of double cavity NBPF with tapered waveguides
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.
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).
 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.
Schematic view of the proposed two port power splitter structure, (b) its transmission spectra with L = 556 nm, W = 356 nm, G = 12 nm, d_{i} = d_{o} = 100 nm, S_{i} = 2, and S_{o} = 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.
(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, d_{i} = 100 nm, d_{o} =56 nm, S_{i} = 1, and S_{o} = 0.5. Magnetic field distributions for wavelengths of (c) 1320 nm, demonstrate the splitting performance of the proposed device.
(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, d_{i} = 100 nm, d_{o} = 48 nm, S_{i} = 1, and S_{o} = 0.5. Magnetic field distributions for wavelengths of (c) 1320 nm, demonstrate the splitting performance of the proposed device.
 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.
(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 S_{i}, for L1 = 628 nm, L2 = 524 nm, W1 = W2 = 164 nm, G = 16 nm, S_{o}=1, CD = 88 nm, and DL = 0.88 um. (e) Quality factor and spectrum peak value (SPV) of both channels, vs. input taper slope, S_{i}, assuming waveguide width, d_{i} = d_{o} = 100 nm, L1 = 628 nm, L2 = 524 nm, cavity widths of W1 = W2 = 164 nm, metallic gap size, G = 16 nm, CD= 88 nm, S_{o} = 0.5. (f) Quality factor and spectrum peak value (SPV) of both channels, versus the distance between two cavity, CD, assuming waveguide width, d_{i} = d_{o} = 100 nm, L1 = 628 nm, L2 = 524 nm, cavity widths of W1 = W2 = 164 nm, metallic gap size, G = 16 nm, S_{i} = 0.5 and S_{o} = 1. (g) Effects of L1 on quality factor and spectrum peak value (SPV) of both channels, assuming L2= 524 nm, waveguide width, d_{i} = d_{o} = 100 nm, cavity widths of W1 = W2 = 164 nm, metallic gap size, G = 16 nm, CD = 88 nm, S_{i} = 0.5, S_{o} = 1. (h) Optimized two wavelength multiplexer for plasmonic dense wavelength division demultiplexing, the structure parameters are: W = 96 nm, d_{i} = d_{o} = 100 nm, L2 = 524 nm, L1 = 532 nm, G=16 nm, S_{i} = 0.5, S_{o} = 1, CD = 120 nm.
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
Effect of changing output ports’ width, do, on filter specifications of twoport demultiplexer. Assuming: d_{i} = 100 nm, L1 = 628 nm, L2 = 524 nm, W1 = W2 = 164 nm, G = 16 nm, S_{i} = 0.5 and S_{o} = 1, CD=88 nm
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
Effect of cavities’ width, W=W1=W2, on filter specifications of twoport demultiplexer. Assuming: d_{i} = d_{o} = 100 nm, L1 = 628 nm, L2 = 524 nm, G = 16 nm, S_{i} = 0.5 and S_{o} = 1, CD = 88 nm
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.
(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, d_{i} = 48 nm, WD = 80 nm, S_{i} = 1, and S_{o} = 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.
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.
Gramotnev D. K.
,
Bozhevolnyi S. I.
2010
“Plasmonics beyond the diffraction limit,”
Nature Photon.
4
83 
91
Zhang J.
,
Zhang L.
2012
“Nanostructures for surface plasmons,”
Advances in Opt. and Photon.
4
157 
321
Zia R.
,
Selker M. D.
,
Catrysse P. B.
,
Brongersma M. L.
2004
“Geometries and materials for subwavelength surface plasmon modes,”
J. Opt. Soc. Am. A
21
2442 
2446
Chu H. S.
,
Ahmed I.
,
Ewe W. B.
,
Li E. P.
2008
“Guiding light in different plasmonic nanoslot waveguides for nanointerconnect application,”
in Proc. 2008 AsiaPacific Symposium on Electromagnetic Compatibility & 19th International Zurich Symposium on Electromagnetic Compatibility
590 
593
Guo Y.
,
Yan L.
,
Pan W.
,
Luo B.
,
Wen K.
,
Guo Z.
,
Li H.
,
Luo X.
2011
“A plasmonic splitter based on slot cavity,”
Opt. Express
19
13831 
13838
Tai C. Y.
,
Chang S. H.
,
Chiu T.
2007
“Design and analysis of an ultracompact and ultrawideband polarization beam splitter based on coupled plasmonic waveguide arrays,”
IEEE Photon. Technol. Lett.
19
1448 
1450
Zhu J. H.
,
Huang X. G.
,
Mei X.
2011
“Improved models for plasmonic waveguide splitters and demultiplexers at the telecommunication wavelengths,”
IEEE Trans. Nanotechnol.
10
1166 
1171
Xiao S.
,
Liu L.
,
Qiu M.
2006
“Resonator narrow bandstop filters in a plasmonpolaritons metal,”
Opt. Express
14
2932 
2937
Hosseini A.
,
Massoud Y.
2007
“Nanoscale surface plasmon based resonator using rectangular geometry,”
Appl. Phys. Lett.
90
181102 
1~1811022
Noual A.
,
Akjouj A.
,
Pennec Y.
,
Gillet J. N.
,
DjafariRouhani B.
2009
“Modeling of twodimensional nanoscale Ybent plasmonic waveguides with cavities for demultiplexing of the telecommunication wavelengths,”
New J. of Phys.
11
103020 
1~10302019
Hu F.
,
Zhou Z.
2011
“Wavelength filtering and demultiplexing structure based on aperturecoupled plasmonic slot cavities,”
J. Opt. Soc. Am. B
28
2518 
2523
Lu H.
,
Liu X. M.
,
Wang L. R.
,
Mao D.
,
Gong Y. K.
2011
“Nanoplasmonic triplewavelength demultiplexers in twodimensional metallic waveguides,”
Appl. Phys. B
103
877 
881
Wen K.
,
Yan L.
,
Pan W.
,
Luo B.
,
Guo Z.
,
Gu Y.
2012
“Wavelength demultiplexing structure based on a plasmonic metalinsulatormetal waveguide,”
J. Opt.
14
075001 
1~0750015
Min C.
,
Veronis G.
2009
“Absorption switches in metaldielectricmetal plasmonic waveguides,”
Opt. Express
17
10757 
10766
Dolatabady A.
,
Granpayeh N.
2012
“All optical logic gates based on two dimensional plasmonic waveguides with nanodisk resonators,”
J. Opt. Soc. Korea
http://dx.doi.org/10.3807/JOSK.2012.16.4.432
16
432 
442
Min B.
,
Ostby E.
,
Sorger V.
,
UlinAvila E.
,
Yang L.
,
Zhang X.
,
Vahala K.
2009
“HighQ surfaceplasmonpolariton whisperinggallery microcavity,”
Nature Lett.
457
455 
459
Noual A.
,
Pennec Y.
,
Akjouj A.
,
DjafariRouhani B.
,
Dobrzynski L.
2009
“Nanoscale plasmon waveguide including cavity resonator,”
J. Phys.: Condens. Matter.
21
375301 
1~3753016
Lee P.
,
Lan Y.
2010
“Plasmonic waveguide filters based on tunneling and cavity effects,”
Springer Science
5
417 
422
Zhang Q.
,
Huang X. G.
,
Lin X. S.
,
Tao J.
,
Jin X. P.
2009
“A subwavelength couplertype MIM optical filter,”
Opt. Express
17
7549 
7555
Huang X.
,
Tao J.
,
Lin X.
2010
Research on NanoPlasmonic Waveguide Filters
IEEE Conference Publishing
Guangzhou, China
Xue W.
,
Guo Y. N.
,
Zhang J.
,
Zhang W.
2009
“Propagation properties of a modified slot surface plasmonic waveguide,”
IEEE J. Lightwave Technol.
7
2634 
2641
Jafarian B.
,
Nozhat N.
,
Granpayeh N.
2011
“Analysis of a triangularshaped plasmonic metalinsulatormetal bragg grating waveguide,”
J. Opt. Soc. Korea
http://dx.doi.org/10.3807/JOSK.2011.15.2.118
15
118 
123
Zhu J. H.
,
Wang Q. J.
,
Shum P.
,
Huang X. G.
2011
“A simple nanometeric plasmonic narrowband filter structure based on metalinsulatormetal waveguide,”
IEEE Trans. Nanotechnol.
10
1371 
1376
Zheng G.
,
Su W.
,
Chen Y.
,
Zhang C.
,
Lai M.
,
Liu Y.
2012
“Band stop filters based on a coupled circular ring metalinsulatormetal resonator containing nonlinear material,”
J. Opt.
14
055001 
1~0550016
Setayesh A.
,
Mirnaziry S. R.
,
Abrishamian M. S.
2011
“Numerical investigation of a tunable bandpass\bandstop plasmonic filter with hollowcore circular ring resonator,”
J. Opt. Soc. Korea
http://dx.doi.org/10.3807/JOSK.2011.15.1.082
15
82 
89
Hu F.
,
Yi H.
,
Zhou Z.
2011
“Bandpass plasmonic slot filter with band selection and spectrally splitting capabilities,”
Opt. Express
19
4848 
4855
Xu Y.
,
Miroshnichenko A. E.
,
Lan S.
,
Guo Q.
,
Wu L. J.
2011
“Impedance matching induce high transmission and flat response bandpass plasmonic waveguides,”
Plasmonics
6
337 
343
Lu H.
,
Liu X.
,
Mao D.
,
Wang L.
,
Gong Y.
2010
“Tunable bandpass plasmonic waveguide filters with nanodisk resonators,”
Opt. Express
18
17922 
17927
Park J.
,
Kim H.
,
Lee I.M.
,
Lee B.
2008
“Plasmonic nano cavity using the cut off property in the metalinsulatormetal waveguide,” paper ThF2
in Proc. The 13th Optoelectronics and Communications Conference (OECC)
Wang T. B.
,
Wen X. W.
,
Yin C. P.
,
Wang H. Z.
2009
“The transmission characteristics of surface plasmon polaritons in ring resonator,”
Opt. Express
17
24096 
24101
Zhu J. H.
,
Shum P.
,
Wang Q. J.
,
Huang X. G.
2011
“A nanoplasmonic highpass wavelength filter based on a metalinsulatormetal circuitous waveguide,”
IEEE Trans. Nanotechnol.
10
1357 
1361
Zhu J. H.
,
Huang X. G.
,
Tao J.
,
Mei X.
2011
“A nanoscale plasmonic longwavelength cutoff filter,”
IEEE Trans. Nanotechnol.
10
817 
821
Lu H.
,
Liu X.
,
Gong Y.
,
Mao D.
,
Wang L.
2011
“Enhancement of transmission efficiency of nanoplasmonic wavelength demultiplexer based on narrow bandstop filters and reflection nanocavities,”
Opt. Express
19
12885 
12890
Tao J.
,
Huang X. G.
,
Zhu J. H.
2010
“A wavelength demultiplexing structure based on metaldielectricmetal plasmonic nanocapillary resonators,”
Opt. Express
18
11111 
11116
Barnes W. L.
,
Dereux A.
,
Ebbesen T. W.
2003
“Surface plasmon subwavelength optics,”
Nature
424
824 
830
Maier S. A.
2007
Plasmonics: Fundamentals and Applications
Springer
Cai W.
,
Shalaev V.
2010
Optical Metamaterials: Fundamentals and Applications
Springer
Tropf W. J.
,
Thomas M. E.
,
Harris T. J.
1995
Handbook of Optics: Devices, Measurements, and Properties
Sponsored by the Optical Society of America
Part 4: Optical and Physical Properties of Materials, Chapter 33: Properties of Crystals and Glasses
Han Z.
,
Van V.
,
Herman W. N.
,
Ho P. T.
2009
“Aperturecoupled MIM plasmonic ring resonators with subdiffraction modal volumes,”
Opt. Express
17
12678 
12684
Johnson P. B.
,
Christy R. W.
1972
“Optical constants of noble metals,”
Phys. Rev. B
6
4370 
4379