## Abstract

Plasmonic reflectors based on serial stub structure are studied in this paper. A general theory of periodic stub structure using transmission line model is developed. The transmission characteristics, e.g., periodicity and symmetry of the spectra, are closely related to the ratio of structure period to stub length. Investigation reveals that the transmission valleys of the spectra could be divided into two categories, which is quite different from conventional Bragg reflectors. Finite-Difference Time-Domain (FDTD) method is used in numerical analysis in this paper.

©2009 Optical Society of America

## 1. Introduction

In recent years, surface plasmon waveguides have become a subject of intensive research because of their outstanding performance in subwavelength optical waveguiding [1, 2]. A variety of metallic structures have been proposed as waveguides in nano scale [3, 4, 5]. Among them, metal-insulator-metal (MIM) waveguide [5] is considered to possess remarkable advantages, such as strong field localization, simplicity and convenience for fabrication and integration into optical circuits. Many efforts have been made to develop nano scale MIM waveguide based plasmonic devices, such as reflectors [6, 7, 8, 9], filters [10] and ring couplers [11]. From recent literatures, most of the research efforts concerning MIM waveguide reflectors are concentrated on Bragg gratings which are realized by periodic change of the cladding layer material, the insulator width and material, respectively [6, 7, 8] or simultaneously [9]. This idea originates from Bragg reflection of X rays from periodic planes in crystalline solids. Actually, many periodic structures have passband and forbidden band characteristics similar to those of Bragg gratings. In this paper, we investigate a novel reflector constructed with periodic stub structure in MIM waveguide.

Stub structure is one of the common elements in microwave networks. In previous studies, it has been introduced to theMIMwaveguide and functions as a wavelength filter with submicron size [12, 13]. The properties of such filters, which consist of a single- or double-stub structure, can be described using the concept of characteristic impedance and transmission line theory [12]. However, they only have good filtering effect for discrete wavelengths. In this paper, we extend the single-stub structure to a periodic one so that it can function like a Bragg reflector whose forbidden band characteristics can not be realized with a single stub. Lin *et al* [13] provided some numerical results about this multiple-stub structure, but failed to give analytical prediction about the center wavelength and width of the forbidden band. In this paper, we use the transmission line and Bloch wave theory to obtain a general analytical equation for this structure. First an infinite periodic stub structure is studied, in which the heat-loss due to metals is assumed to be negligible. Then, Finite-Difference Time-Domain (FDTD) method [14] is used to validate our theoretic analysis and investigate the transmission characteristics, such as the transmission valleys, the periodicity and symmetry of spectrum.

## 2. Serial stub structure

We consider a MIM waveguide with a series of stubs perpendicular to it as shown in Fig. 1. The width *w* of waveguide and stubs are set to be equal. The length of the stubs and the interval between two stubs are denoted by *L* and *d*, respectively. The insulator medium of the waveguide is assumed to be air whose permittivity is *ε _{i}*=1. The circumambience is silver. The frequency-dependent complex relative permittivity of silver is characterized by the Drude model

*ε*(

_{m}*ω*)=

*ε*∞-

*ω*/(

_{p}*ω*

^{2}+

*iωγ*). Here,

*ω*=1.38×10

_{p}^{16}Hz is the bulk plasma frequency,

*γ*=2.73×10

^{13}Hz is the damping frequency of the oscillations,

*ω*is the angular frequency of the incident electromagnetic radiation, and

*ε*

_{Ȟ}stands for the dielectric constant at infinite angular frequency with a value of 3.7. The Finite-Difference Time-Domain (FDTD) method is used for numerical simulation. The calculated area is divided into uniform Yee cells of Δ

*x*=Δ

*y*=2 nm and surrounded by perfectly matched layer (PML) absorbing boundary.

#### 2.1. Modeling of lossless periodic stub structure

Assuming that silver is lossless, which means that the damping due to the metal is ignored, the properties of systems such as serial stub structure can be described using the concept of a distributed constant circuit including a loss-free transmission line with a characteristic impedance of ${Z}_{0}=\genfrac{}{}{0.1ex}{}{{k}_{\mathrm{sp}}w}{\omega {\epsilon}_{i}}$ [15], where *k*
_{sp} is the propagation constant of the fundamental propagating TM mode in the MIM waveguide. Based on transmission line theory, the stub structure is equivalent to an open-circuited transmission line [12, 16] with characteristic susceptance *Y _{s}*=

*jY*

_{0}tan(

*k*

_{sp}

*L*), where

*j*=√-1 and

*Y*

_{0}=1/

*Z*

_{0}. An equivalent circuit of the serial stub structure is presented in Fig. 2. Each unit cell of the line consists of a length

*d*of transmission line with a shunt susceptance across the midpoint of the line. The voltage

*V*at the edge of the cell can be written as:

_{n}*V*=

_{n}*V*

^{+}

_{n}exp(-

*jk*

_{sp}

*x*)+

*V*

^{-}

_{n}exp(

*jk*

_{sp}

*x*). We can relate the voltages on either side of the

*n*th unit cell using a 2×2 transfer matrix:

where matrix *M*
_{n,n+1} can be written as

$$=\left[\begin{array}{cc}\multicolumn{1}{c}{{e}^{{\mathrm{jk}}_{\mathrm{sp}}d}\left(1+\genfrac{}{}{0.1ex}{}{j}{2}\mathrm{tan}\left({k}_{\mathrm{sp}}L\right)\right)}& \multicolumn{1}{c}{\genfrac{}{}{0.1ex}{}{j}{2}\mathrm{tan}\left({k}_{\mathrm{sp}}L\right)}\\ \multicolumn{1}{c}{-\genfrac{}{}{0.1ex}{}{j}{2}\mathrm{tan}\left({k}_{\mathrm{sp}}L\right)}& \multicolumn{1}{c}{{e}^{{-\mathrm{jk}}_{\mathrm{sp}}d}\left(1-\genfrac{}{}{0.1ex}{}{j}{2}\mathrm{tan}\left({k}_{\mathrm{sp}}L\right)\right)}\end{array}\right]$$

For a serial stub structure with *N* cells, the calculation of the transmission spectrum involves multiplication of *N* transfer matrices: Π^{N}
_{i=1}
*M*
_{i,i+1}.

For an infinite periodic stub structure, the waves at the edges of the unit cells is similar to Bloch waves that propagate through periodic crystal lattices. The Bloch wave dispersion relation for this periodic stub structure can be written as

where *K*=*α*+*jβ* is the Bloch wave number. Since the right-hand side of Eq. (3) is purely real, we must have either *β*=0 or *αd*=*απ* (*m* is a non-negative integer). The case of *β*=0, *α*≠0 corresponds to a nonattenuating propagating wave. The magnitude of the right-hand side of Eq. (3) is less than or equal to unity. For the case of *αd*=*mπ*, *β*≠0, the wave does not propagate, but is reflected completely. Eq. (3) is reduced to

This corresponds to the forbidden band of the structure. Within the forbidden band, *βd* varies from zero at the band edges to its maximum value at the center of the band.

Using the multiplication operation of *N* transfer matrices Eq.(2), we can get the optical transmittance through *N* stubs in the periodic stub structure

where

From Eqs. (5) and (6), it can be seen that there exist some low transmission points when *F* →±∞ which requires

or

The solutions of Eq. (7) are simply

where *m* is a non-negative integer. In the forbidden band, *Kd* is a complex number *Kd*=*mπ* +*jβd*, the left-hand side of Eq. (8) becomes sinh^{2}(*βd*)/sinh^{2}(*Nβd*) which approaches zero exponentially as *e*
^{-2(N-1)βd}.

Strictly speaking, there is no rigorous Bloch wave and forbidden band for a finite serial stub structure which is constructed with lossy metal. To investigate the reflection characteristics of these “real” stub structures, FDTD method is used in the following section.

#### 2.2. Simulation and results

From Eqs. (5) and (6), it is clear that the transmission spectra are closely related to *L* and *d*. For simplicity, we investigate the case of *d/L*=1 first. Eq. (3) is reduced to cos(*Kd*)=[3cos(*k*
_{sp}
*L*)-sec(*k*
_{sp}
*L*)]/2. The period of the transmission spectra is *π* and the symmetric center of the first period is located at *k*
_{sp}
*L*=*π*/2. To demonstrate the validity of our analysis, we construct a 4 stub structure with “real” lossy metal. Other geometric parameters are set as: *w*=50 nm, *d*=L=400 nm. FDTD method is used to calculate the transmission spectrum which covers the first two periods. Transfer matrix method (TMM) deduced from Eq.(1) and Eq. (2) is also used for lossless case as a comparison. The calculation results, which are depicted in Fig.3(a), agree very well with each other. The main deviation only exists at the high transmission regions. The transmission valleys are located around *k*
_{sp}
*L*=0.5*π* and 1.5*π*, and the spectra possess good periodicity and symmetry as expected. For comparison, the dispersion relation Eq. (3) is solved numerically and shown in Fig. 3(c). Although there is no strict forbidden band for finite serial stub structure, the low transmission regions in Fig. 3(a) correspondent to the forbidden bands in Fig. 3(c). The serial stub structure with *L*=*d* can be used as a wavelength selective reflector whose transmission spectrum is analogous to that of the conventional Bragg gratings.

For another case of *d/L*=3/2, Eq. (3) can be rewritten as cos(*Kd*)=cos(3*k*
_{sp}
*L*/2)-(1/2) tan(*k*
_{sp}
*L*) sin(3*k*
_{sp}
*L*/2). The period of *k*
_{sp}
*L* for transmission spectrum is 2*π*, and the symmetric center of the first period is *k*
_{sp}
*L*=*π*. The transmission spectra with *L*=400 nm and *d*=600 nm are shown in Fig. 3(b). The corresponding dispersion relation is depicted in Fig. 3(d). According to the forbidden band condition Eq.(4), centers of the first four forbidden bands are located at *k*
_{sp}
*L*=0.5*π* 0.732*π*, 1.268*π* and 1.5*π*. This agrees well with the simulation results. In Fig.3(b), it is also worth noting that the transmittances at the first and fourth valleys are obviously lower than those of the second and third ones. To have a clearer comparison, more cases with different stub number N are studied. The transmission spectra by FDTD with *N*=1, 2, 4 and 8 are presented together in Fig. 4. As *N* decreases from 8 to 1, the transmission coefficient at the valley points increases gradually. This is analogous to the conventional Bragg reflector. But at *k*
_{sp}
*L*=0.5*π* and 1.5*π*, the transmission coefficients are fixed at zero no matter how much the stub number changes. According to Eq. (5), there are two categories of low transmission points which are described by Eqs. (7) and (8), respectively. The transmission coefficient *T* reaches zero, if Eq. (7) is satisfied, while *T* only approaches zero with an increasing *N* at the situation of Eq. (8). We present the optical field at the first two valley points (λ=1503 nm and λ=2191 nm) in Figs. 4(b) and (c), respectively. At λ=1503 nm (Fig. 4(b)), the optical field vanishes gradually after about 4 stubs, while at λ=2191 nm (Fig. 4(c)), the transmission becomes negligible just after the first stub, which is simply the single-stub case presented in Ref.[12]. The existence of these zero transmission points make the transmission characteristics of the serial stub structure quite different from those of conventional Bragg reflectors.

In more general cases *d/L*=*q/p*, where *q* and *p* are positive integers, it is easy to see that the transmittance *T* in Eq.(5) is a periodic function of *k*
_{sp}
*L* with period of *pπ*. The numerator and denominator of Eq.(6) have different period (*π* and *pπ*, respectively), which means the two categories of transmission valleys have different intervals. If *d/L* can not be expressed as a fraction of two integers, the transmission spectrumlosses its periodicity, or the period is infinity.

It can also be seen from the numerical results that the loss-free approximation is reasonable for micro-dimensional serial stub structure. In Fig. 3, this model can predict accurately the positions of both the forbidden and the pass bands for *L* and *d* with submicron size. In practical photonic circuits, to shift the first forbidden band center to the desired telecommunication wavelengths (1550 nm or 1310 nm), the structure size (involving *d* and *L*) will be reduced further and the heat-loss will decrease consequently. Furthermore, at the incident wavelength 1550 nm, light intensity of the primary propagation TMmode drops at about 0.1dB/micron in straight Silver-Air-Silver waveguide (50 nm wide). Therefore, the model presented here can be used to analyze the transmission characteristics of the serial stub structure in MIM waveguide and facilitate the designs of stub based reflectors with micro-dimensional size in telecommunication circuits.

## 3. Conclusion

In this paper, the plasmonic reflectors based on serial stub structure are studied both theoretically and numerically. A generalmodel of periodic stub structure using transmission line model is developed. The transmission characteristics, e.g., periodicity and symmetry of the spectra, are closely related to the ratio of structure period to stub length. Investigation reveals that the transmission valleys of the spectra could be divided into two categories, which is quite different from conventional Bragg reflectors. It is expected that the model and numerical results may provide useful information in the design of serial stub based reflectors in photonic circuits.

## Acknowledgments

This work was supported by the National Natural Science Foundation of China under grants 10674038 and 10604042 and National Basic Research Program of China under grant 2006CB302901.

## References and links

**1. **W. L. Barnes, A. Dereux, and T. W. Ebbesen, “Surface plasmon subwavelength optics,” Nature **424**824–830 (2003).
[CrossRef] [PubMed]

**2. **E. Ozbay, “Plasmonics: Merging photonics and electronics at nanoscale dimensions,” Science **311**189–193 (2006).
[CrossRef] [PubMed]

**3. **D. F. P. Pile, T. Ogawa, D. K. Gramotnev, T. Okamoto, M. Haraguchi, M. Fukui, and S. Matsuo, “Theoretical and experimental investigation of strongly localized plasmons on triangular metal wedges for subwavelength waveguiding,” Appl. Phys. Lett. **87**, 061106 (2005).
[CrossRef]

**4. **D. F. P. Pile and D. K. Gramotnev, “Channel plasmon-polariton in a triangular groove on a metal surface,” Opt. Lett. **29**, 1069–1071 (2004).
[CrossRef] [PubMed]

**5. **J. A. Dionne, L. A. Sweatlock, and H. A. Atwater, “Plasmon slot waveguides: Towards chip-scale propagation with subwavelength-scale localization,” Phys. Rev. B **73**, 035407 (2006).
[CrossRef]

**6. **B. Wang and G.P. Wang, “Plasmon Bragg reflectors and nanocavities on flat metallic surfaces,” Appl. Phys. Lett. **87**, 013107 (2005).
[CrossRef]

**7. **A. Hossieni and Y. Massoud, “A low-loss metal-insulator-metal plasmonic bragg reflector,” Opt. Express **14**, 11318–11323 (2006).
[CrossRef] [PubMed]

**8. **Z. Han, E. Forsberg, and S. He, “Surface plasmon Bragg gratings formed in metal-insulator-metal waveguides,” IEEE Photon. Technol. Lett. **19**, 91–93 (2007).
[CrossRef]

**9. **J. Liu, L. Wang, M. He, W. Huang, D. Wang, B. Zou, and S. Wen, “A wide bandgap plasmonic Bragg reflector,” Opt. Express **16**, 4888–4894 (2008).
[CrossRef] [PubMed]

**10. **Q. Zhang, X. Huang, X. Lin, J. Tao, and X. Jin, “A subwavelength coupler-type MIM optical filter,” Opt. Express **17**, 7549–7555 (2009).
[CrossRef]

**11. **S. Xiao, L. Liu, and M. Qiu, “Resonator channel drop filters in a plasmon-polaritons metal,” Opt. Express **14**, 2932–2937 (2006).
[CrossRef] [PubMed]

**12. **Y. Matsuzaki, T. Okamoto, M. Haraguchi, M. Fukui, and M. Nakagaki, “Characteristics of gap plasmon waveguide with stub structures,” Opt. Express **16**, 16314–16325 (2008).
[CrossRef] [PubMed]

**13. **X. Lin and X. Huang, “Numerical modeling of a teeth-shaped nanoplasmonic waveguide filter,” J. Opt. Soc. Am. B **26**, 1263–1268 (2009).
[CrossRef]

**14. **A. Taflove and S. C. Hagness, *Computational Electrodynamics: The Finite-Difference Time-Domain Method*, 2nd ed. (MA: Artech House, Boston, 2000).

**15. **G. Veronis and S. Fan, “Bends and splitters in metal-dielectric-metal subwavelength plasmonic waveguides,” Appl. Phys. Lett. **87**, 131102 (2005).
[CrossRef]

**16. **S. E. Kocabas, G. Veronis, D. A. B. Miller, and S. Fan, “Transmission Line and Equivalent Circuit Models for Plasmonic Waveguide Components,” IEEE J. Sel. Top. Quantum Electron. **14**, 1462–1472 (2008).
[CrossRef]