## Vertical coupling between gap plasmon waveguides

Optics Express, Vol. 16, Issue 17, pp. 12677-12687 (2008)

http://dx.doi.org/10.1364/OE.16.012677

Acrobat PDF (453 KB)

### Abstract

This work examines vertical coupling between gap plasmon waveguides for use in high confinement power transfer and power splitting applications at 1.55 µm free space wavelength. The supermode interference method is used to obtain key coupler performance parameters such as coupling length, extinction ratio, net coupled output power, radiated power, and reflected power as a function of waveguide center-to-center spacing, core refractive index, and gap width. The initial power distribution among the two coupler supermodes is obtained via the mode matching method for a single input waveguide feed. Excellent agreement with three-dimensional finite difference time domain simulations is observed for the case of square 50 nm gaps with core refractive indices of 2.50 and a center-to-center spacing of 112 nm. Local maxima in the net coupled output power are found to coincide with local minima in the coupling length. An increase in the core refractive index from 1.00 to 2.5 increases the local maximum net coupled output power from 6.4% to 49% but decreases the extinction ratio from 12.7 to 6.94. A sweep of the width of the core from 25 to 100 nm increases the net coupled output power from 43.7% to 52.0%, but increases the coupling length from 1.58 to 3.19 µm and decreases the extinction ratio from 7.39 to 6.57.

© 2008 Optical Society of America

## 1. Introduction

1. J. A. Conway, S. Sahni, and T. Szkopek, “Plasmonic interconnects versus conventional interconnects: a comparison of latency, cross-talk and energy costs,” Opt. Express **15**, 4474–4484 (2007). [CrossRef] [PubMed]

2. K. Jung, F. L. Teixeira, and R. M. Reano, “Au/SiO_{2} nanoring plasmon waveguides at optical communication band,” IEEE J. Lightwave Tech. **25**, 2757–2765 (2007). [CrossRef]

3. P. Berini, “Air gaps in metal stripe waveguides supporting long-range surface plasmon polaritons,” J. Appl. Phys. **102**, 33112 (2007). [CrossRef]

4. G. Veronis and S. Fan, “Guided subwavelength plasmonic mode supported by a slot in a thin metal film,” Opt. Lett. **30**, 3359–3361 (2005). [CrossRef]

*L*) is defined as the length over which the component of the time-average power in the direction of propagation of a mode decays to 1/e of its original value. The confinement is less strictly defined and often depends on the waveguide geometry, though usually it is stated to be the width or area of the region in which the power density is greater than some proportion of its maximum value [7

_{p}7. R. Buckley and P. Berini, “Figures of merit for 2D surface plasmon waveguides and application to metal stripes,” Opt. Express **15**, 12174–12182 (2007). [CrossRef] [PubMed]

2. K. Jung, F. L. Teixeira, and R. M. Reano, “Au/SiO_{2} nanoring plasmon waveguides at optical communication band,” IEEE J. Lightwave Tech. **25**, 2757–2765 (2007). [CrossRef]

8. P. Berini, “Plasmon-polariton waves guided by thin lossy metal films of finite width: Bound modes of symmetric structures,” Phys. Rev. B **61**, 10484–10503 (2000). [CrossRef]

9. R. Zia, J. A. Schuller, A. Chandran, and M. L. Brongersma, “Plasmonics: the next chip-scale technology,” Mat. Today **9**, 20–27 (2006). [CrossRef]

10. S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, and T. W. Ebbesen, “Channel plasmon-polariton guiding by subwavelength metal grooves,” Phys. Rev. Lett. **95**, 46802 (2005). [CrossRef]

11. M. Yan and M. Qiu, “Guided plasmon polariton at 2D metal corners,” J. Opt. Soc. Am. B **24**, 2333–2342 (2007). [CrossRef]

4. G. Veronis and S. Fan, “Guided subwavelength plasmonic mode supported by a slot in a thin metal film,” Opt. Lett. **30**, 3359–3361 (2005). [CrossRef]

4. G. Veronis and S. Fan, “Guided subwavelength plasmonic mode supported by a slot in a thin metal film,” Opt. Lett. **30**, 3359–3361 (2005). [CrossRef]

**30**, 3359–3361 (2005). [CrossRef]

**30**, 3359–3361 (2005). [CrossRef]

12. G. Veronis and S. Fan, “Modes of subwavelength plasmonic slot waveguides,” IEEE J. Lightwave Tech. **25**, 2511–2521 (2007). [CrossRef]

13. G. Veronis and S. Fan, “Crosstalk between three-dimensional plasmonic slot waveguides,” Opt. Express **16**, 2129–2140 (2008). [CrossRef] [PubMed]

14. H. S. Won, K. C. Kim, S. H. Song, C. Oh, P. S. Kim, S. Park, and S. Kim, “Vertical coupling of long-range surface plasmon polaritons,” Appl. Phys. Lett. **88**, 11110 (2006). [CrossRef]

## 2. Design approach

### 2.1 Geometry specification

*h*. Propagation of the wave occurs in the +z-direction, in which the waveguide geometry is considered to be invariant for the purpose of mode solving. The refractive indices of the core, metal, and cladding are denoted via

*n*,

_{core}*n*, and

_{metal}*n*, respectively.

_{clad}*s*. The refractive index of the middle spacer layer between the two metal films is denoted

*n*. This structure is referred to as a vertical gap plasmon coupler (VGPC). We take the geometry to be vertically-symmetric. To quantify how the total power is distributed between the two waveguides as a function of propagation distance along the coupler, we define the channel 1 region as being that above the line midway between the two gaps, and the channel 2 region as being below it, and thus this line is the boundary between the two semi-infinite waveguide cross-section domains over which the integration of time average power density is done (dashed line in Fig. 1(b)).

_{mid}### 2.2 Simulation method

*e*

^{-γz}for forward propagating waves, where

*γ*=

*α*+

*iβ*is the complex propagation constant,

*α*is the attenuation constant and

*β*is the phase constant. Thus the propagation length is

*L*=1/(2

_{p}*α*). The real part of the effective propagation index (

*n*) is obtained by dividing

_{eff}*β*by

*k*

_{0}, where

*k*

_{0}is the free space phase constant. The real part of the effective propagation index and the propagation length of the fundamental gap plasmon mode are shown in Fig. 2(a) and Fig. 2(b) calculated via the finite difference frequency domain (FDFD) method [15

15. J. A. Pereda, A. Vegas, and A. Prieto, “An improved compact 2D full-wave FDFD method for general guided wave structures,” Microwave Opt. Technol. Lett. **38**, 331–335 (2003). [CrossRef]

12. G. Veronis and S. Fan, “Modes of subwavelength plasmonic slot waveguides,” IEEE J. Lightwave Tech. **25**, 2511–2521 (2007). [CrossRef]

12. G. Veronis and S. Fan, “Modes of subwavelength plasmonic slot waveguides,” IEEE J. Lightwave Tech. **25**, 2511–2521 (2007). [CrossRef]

*w*=

*h*=50 nm,

*n*=(-143-10i)

_{metal}^{0.5}, refractive index of 1.44 for core, cladding, and substrate) [16]. For both FDFD and FEFD formulations,

*γ*is obtained as the eigenvalue of the respective mode. Our FDFD implementation uses a uniform orthogonal grid whose points sit on the geometry boundaries. We enforce homogeneous Dirichlet boundary conditions at the simulation edges and keep them far enough away from the waveguide so as not to have a numerically significant effect. Both the electric (E-FDFD) and magnetic (H-FDFD) field formulations were implemented as a cross-check. For the COMSOL simulations, we use the scattering boundary condition to truncate the boundary, use cubic Lagrangian elements, and keep the maximum element size along the vertical sidewalls of the GPW equal to 0.1

*h*. There is good agreement among the four approaches. The FEFD method was found to converge on a solution to a given accuracy with fewer unknowns and computational effort compared to our FDFD implementation with a uniform grid. Further numerical computation in this paper is accomplished via FEFD.

### 2.3 Parameter space reductions

17. H. J. Hagemann, W. Gudat, and C. Kunz, “Optical constants from the far infrared to the x-ray region: Mg, Al, Cu, Ag, Au, Bi, C, and Al_{2}O_{3},” J. Opt. Soc. Am. **65**, 742–744 (1975). [CrossRef]

*n*≈(0.434-12.7i). The substrate, cladding and middle are composed of silica with a refractive index of 1.44. The gap cores can have refractive indices of 1.00, 1.44, 2.00, and 2.50. Materials that span this range of refractive indices include air, silicon dioxide, silicon nitride, and Group IV–VI amorphous glasses. The fabrication of air gaps has been feasibly demonstrated using multiple techniques [18

_{metal}18. R. M. Reano and S. W. Pang, “Sealed three-dimensional nanochannels,” J. Vac. Sci. Technol. B **23**, 2995–2999 (2005). [CrossRef]

19. Bruce McConnel, “Using self-assembly to create airgap microprocessors” (IBM, 2007). http://www-03.ibm.com/press/us/en/presskit/21463.wss. Accessed 10-30-07.

**25**, 2511–2521 (2007). [CrossRef]

20. D. F. P. Pile, D. K. Gramotnev, M. Haraguchi, T. Okamoto, and M. Fukui, “Numerical analysis of coupled wedge plasmons in a structure of two metal wedges separated by a gap,” J. Appl. Phys. **100**, 13101 (2006). [CrossRef]

### 2.4 Analysis method

*E*for the coordinate system that we are considering. We denote the symmetric supermode with

_{x}*E*and the asymmetric supermode with

_{s}*E*. The phase constants for the symmetric and asymmetric supermodes are donated

_{a}*β*and

_{s}*β*, respectively.

_{a}*E*. The input plane of the coupler is our phase reference point.

_{i}*E*,

_{i}*E*, and

_{s}*E*for

_{a}*w*=

*h*=50 nm and

*s*=100 nm, with

*n*=1.44. The white lines show the boundaries of the metal and the core, the surface plot shows the time average z-propagating power density (

_{core}*S*) profile of the modes, and the arrow plot shows the transverse electric field vector. All of the quantities are normalized to the peak value. From the power density profile, one can see that the wave’s power density peaks at the four corners of each of the gaps, demonstrating that the dominant guiding mechanism in this particular case is due to the coupled “wedge plasmon” mode [20

_{avgz}20. D. F. P. Pile, D. K. Gramotnev, M. Haraguchi, T. Okamoto, and M. Fukui, “Numerical analysis of coupled wedge plasmons in a structure of two metal wedges separated by a gap,” J. Appl. Phys. **100**, 13101 (2006). [CrossRef]

*E*(Fig. 4(b)) and

_{s}*E*(Fig. 4(c)) are as shown, they approximately superimpose to

_{a}*E*(Fig. 4(a)). When the phase difference between them is 180°, the input mode is approximately replicated at the bottom waveguide position, at y=-50 nm.

_{i}*L*=π/|

_{c}*β*-

_{s}*β*|) and is half the beat length between

_{a}*E*and

_{s}*E*. This length sets the minimum required length of any coupler, and is itself a performance parameter. Given that GPW’s are highly lossy and that their advantage rests with their ability to confine radiation in a subwavelength geometry, it is desirable that this length be as short as possible for compact, efficient devices. We take the effective end of the coupler to be at the plane

_{a}*z*=

*L*and evaluate performance parameters at this location.

_{c}*E*and

_{s}*E*in the coupler, calling them

_{a}*t*and

_{s}*t*, respectively, normalizing each mode such that

_{a}*P*is equal to one.

_{avgz}*E*and

_{s}*E*in the coupler, calling them

_{a}*t*and

_{s}*t*, respectively, normalizing each mode such that

_{a}*P*is equal to one.

_{avgz}*t*superscript denotes the transverse components of the field. The parameter

*r*denotes reflection coefficient. We then apply the unconjugated (applicable for lossy modes) orthonormality condition as the overlap integral using the output modes of the waveguide to reduce (1) and (2) to scalar equations [24]:

*n*and

*m*denote modes of the same waveguide, and

*δ*is the Kronecker delta. Using Eq. (3), Eqs. (1) and (2) are transformed to a set of four equations that can be written in the following matrix form:

_{nm}*i*,

*a*, and

*s*subscripts label the input, asymmetric, and symmetric supermodes, respectively. Since this is an over-determined system, we take the pseudo-inverse of the matrix on the left of (4) to find the scattering coefficients to obtain the best (least squares) approximation [25]. Each

*ξ*is the overlap between the electric field of mode

_{nm}*n*and the magnetic field of mode

*m*, where the modes here are from different waveguides:

### 2.5 Performance parameters

*P*), over from one waveguide to another. One performance parameter of the directional coupler is the extinction ratio (

_{avgz}*ER*), which is the ratio of the output channel power (the time-averaged

*z*-component of the power density,

*S*, integrated over the channel 2 region) to the input channel (channel 1) power and is ideally infinite so that the desired output channel is isolated from the other. Since we are interested in extracting as much power in channel 2 as possible, we define another performance parameter to be the net coupled output power,

_{avgz}*P*, which is the channel 2 power at the output plane, taking into account both the losses at the input of the coupler due to reflected and radiated power, and the absorption losses along the coupler. This intrinsically includes information about

_{net}*ER*as well because a large

*ER*coupler will have most of its power in channel 2. We evaluate both

*P*and

_{net}*ER*at a length along the coupler equal to

*L*. We have always found that the

_{c}*ER*is maximized at this point, though

*P*is actually maximized at a slightly shorter distance. We also examine the

_{net}*guided*reflected power (

*P*) and total estimated radiated power (

_{refl}*P*) in terms of the transmission and reflection coefficients:

_{rad}### 2.6 Design strategy

*s*, requires consideration that the electric field of

*E*between the two metal films (to either side of the gaps) points vertically, hence the electric field resembles that of the fundamental mode (TM

_{a}_{0}) of the two-dimensional MIM slab. Previous work has shown that the TM

_{0}mode has no cutoff and that its phase constant increases without bound as the slab width is reduced [26

26. R. Zia, M. D. Selker, P. B. Catrysse, and M. L. Brongersma, “Geometries and materials for subwavelength surface plasmon modes,” J. Opt. Soc. Am. A. **21**, 2442–2446 (2004). [CrossRef]

27. S. Collin, F. Pardo, and J. Pelouard, “Waveguiding in nanoscale metallic apertures,” Opt. Express **15**, 4310–4320 (2007). [CrossRef] [PubMed]

*E*can therefore couple to the semi-infinite slab mode for small enough

_{a}*s*[13

13. G. Veronis and S. Fan, “Crosstalk between three-dimensional plasmonic slot waveguides,” Opt. Express **16**, 2129–2140 (2008). [CrossRef] [PubMed]

_{0}slab mode and the

*E*,

_{i}*E*, and

_{s}*E*effective indices for a coupler with

_{a}*n*=1.44 and

_{core}*w*=h=50 nm obtained via FEFD simulations. Both

*E*and TM

_{a}_{0}converge on each other for

*s*values below about 70 nm and

*E*is effectively cutoff for

_{a}*s*values below 70 nm. At

*s*=80 nm, with the corners of the gaps rounded with a 10 nm radius (weakening the control over the mode that the corners have), and with a simulation domain that was 32 µm wide, we find that the FWHM of

*S*of

_{avgz}*E*along the

_{a}*x*-direction at

*y*=0 is only about 2 µm (not cutoff). We do not observe any cutoff for

*E*for

_{s}*s*values greater than 55nm (corresponding to a metal-metal spacing of 5 nm), even though the transverse electric field off to the sides of the gaps resembles the TM

_{1}slab. This is consistent with the fact that the TM

_{1}slab mode is cutoff at this slab width. As expected, the phase constants of

*E*and

_{s}*E*approach that of

_{a}*E*(the single gap), for large gap spacings. Furthermore, we find that the real part of the

_{i}*n*’s of

_{eff}*E*and

_{s}*E*cross each other at

_{a}*s*≈110 nm leading to an infinite coupling length, as can be seen in Fig 5(b). Also, putting the gaps either less than 110 nm apart or at about 190 nm leads to large

*n*differences and short

_{eff}*L*’s.

_{c}*S*in Eq. (8), where the “tot” suffix denotes a field quantity that is the weighted sum of the two supermodes:

_{avgz}*α*′ and

*β*≡

_{Δ}*β*-

_{s}*β*. The power densities, S

_{a}_{avgz-nm}, are given by

*S*and

_{avgz-ss}*S*are the mode power densities of each supermode by itself. However,

_{avgz-aa}*S*and

_{avgz-sa}*S*are the cross power densities. The power in each channel is obtained by integrating (8) over the regions specified in Fig. 1(b).

_{avgz-as}## 3. Optimization

*n*, the CTC gap spacing,

_{core}*s*, and the gap width,

*w*. We vary

*s*from 60 to 500 nm, take

*w*=25, 50, or 100 nm. All other parameters are as specified in section 2.3 above.

### 3.1 Square gaps

*w*=h=50 nm. The coupling length in general decreases as the waveguides are brought closer together, but there is a singularity (denoted

*s*) where the supermode phase constants cross each other. There is a local minimum in

_{inf}*L*for

_{c}*s*values larger than

*s*. An infinite coupling length, however, does not necessarily imply complete isolation. Even when the phase constants are the same, power is transferred by the more rapid attenuation of one supermode versus another (

_{inf}*E*has always been observed to have a higher loss constant for these analyses), in which case, at the limit of large coupling length, almost half of the total power would be transferred. The ER is shown in Fig. 6(b). ER also decreases with

_{a}*s*, exhibiting a local minimum near

*s*where a close to unity extinction ratio is obtained. Strictly speaking the

_{inf}*ER*is not defined at

*s*, so only

_{inf}*ER*values away from

*s*are considered. Also shown is

_{inf}*P*in Fig. 6(c), which decreases monotonically as the waveguides are separated further apart, corresponding to a better match of the supermodes to the single input mode. In all cases tested it was found that

_{rad}*P*was less than 0.5%.

_{refl}*s*values greater than

*s*. For lower

_{inf}*s*values, the

*ER*and the modal loss overcome further decreases in the

*L*. It is observed that there is a close correspondence between the center-to-center (CTC) spacing that maximizes the net coupled output power, denoted

_{c}*s*, and the relative minimum in the

_{opt}*L*, for each

_{c}*n*tested. This is due to the high loss of these waveguides dominating the output characteristics for the range of

_{core}*s*values above

*s*, where the

_{inf}*ER*is still appreciable, and causing shorter

*L*’s to be favored over large

_{c}*ER*’s. The

*P*is maximized at

_{net}*s*=112, 136, 188, and 296 nm for

_{opt}*n*=1.00, 1.44, 2.00, and 2.50, respectively. Figure 6(e) is a plot of the average supermode power propagation length,

_{core}*L*.

_{p-avg}*n*produces a smaller minimum

_{core}*L*, larger

_{c}*ER*, smaller

*s*and greater

_{opt}*P*. This can be explained by the fact that the increased

_{net}*ER*and decreased

*L*are able to overcome the decrease in the propagation length. Figure 6(d) shows that the increase in the maximum

_{c}*P*increases at a decreasing rate with

_{net}*n*. Because increasing

_{core}*n*also leads to more highly confined modes,

_{core}*P*and

_{rad}*P*also decrease.

_{refl}*n*=2.5 and

_{core}*s*=

*s*=112 nm. The power in channel 1 and 2 is denoted

_{opt}*P*and

_{ch1}*P*, respectively. The FDTD simulation consists of an input waveguide feeding into the coupled structure, with the input mode field profile generated by COMSOL as the excitation. Our simulation setup was exactly as depicted in Fig. 3, using the input mode as a CW source. We extracted the power values along the coupler length by projecting the total FDTD field onto the weighted sum of the coupler supermodes obtained from COMSOL, using overlap integrals (similar to those used to obtain Eq. (4) from Eqs. (1) and (2)) to obtain the weights and normalizing the weights to the input mode power at coupler input plane. For this case, |

_{ch2}*t*|=0.8024 and |

_{s}*t*|=0.5687 (both with negligible imaginary parts). Our results show that not only does the SI method given here give the correct propagation length, but also that the MMM, using the pseudo-inverse of the over-determined matrix to find the scattering parameters at the discontinuity, correctly calculates initial power distribution of the supermodes at the coupler input. A benefit of the SI-MMM approach is that it transforms a 3D problem into a 2D one, allowing one to analyze the structure in a matter of minutes instead of days (on a desktop computer).

_{a}*P*is actually maximized a little bit before the coupling length. This is due to the high loss of the waveguide causing the total propagating power of the waveguide to decay before the amount transferred to channel 2 has maximized. However, one can see that the difference in powers and coupler lengths is small, so we simply evaluate

_{ch2}*P*at

_{ch2}*L*. In all cases that we studied, we found that the

_{c}*ER*along the length of the coupler was maximized at

*z*=

*L*as well, and this could be advantageous for noise reduction.

_{c}### 3.2 Comparison of square gaps to rectangular gaps

*w*=25 and 100 nm (h=50 nm). The relationships of the performance parameters with respect to

*n*and

_{core}*s*, shown in Fig. 6, mainly carry over to the cases where the gap is no longer square. Because we have found that

*P*is maximized at larger values for larger

_{net}*n*, we summarize the results obtained in Table 1, with

_{core}*n*=2.50.

_{core}*w*causes the modes to be less-highly confined in the gaps, and so

*E*is more easily perturbed by the proximity of another metal layer. Thus, the

_{a}*s*values tend to increase with

_{opt}*w*. The coupling length also increases, doubling from about 1.5 µm to 3 µm with

*w*increasing from 25 to 100 nm, this being due to the magnitude of the propagation constants being smaller. The

*ER*does decrease by a smaller amount, and with the

*L*increasing significantly, this causes

_{p-avg}*P*to increase by 9 %, even with the longer coupling length. This means that a designer may decide to use a wider gap to increase the output power at the small expense of a decreased extinction ratio. The radiated power stays around 3%, and the reflected power is negligible.

_{net}## 4. Conclusions

*ER*and resulting in a maximum in the net coupled output power. This maximum corresponds well to the local minimum seen in the coupling length versus center-to-center spacing curves and also occurs at smaller center-to-center spacing values for higher core refractive indices. Making the width of the core larger from 25 to 100 nm makes the modes less confined, increasing the coupling length and the net coupled output power from 1.58 to 3.19 and from 43.7% to 52.0%, respectively, but decreasing the extinction ratio from 7.39 to 6.57. Supermode propagation length can be increased by making the core refractive index low which can be understood in terms of making the magnitude of the complex propagation constant smaller. This structure could be used as a power splitter at half the coupling length. Even though the losses are high, long propagation lengths are not necessary as the gap plasmon waveguide has a high enough confinement that it could be used to connect micron-sized low-loss dielectric waveguides to deep submicron transistor circuits. The tradeoff for confinement is reduced efficiency.

## References and links

1. | J. A. Conway, S. Sahni, and T. Szkopek, “Plasmonic interconnects versus conventional interconnects: a comparison of latency, cross-talk and energy costs,” Opt. Express |

2. | K. Jung, F. L. Teixeira, and R. M. Reano, “Au/SiO |

3. | P. Berini, “Air gaps in metal stripe waveguides supporting long-range surface plasmon polaritons,” J. Appl. Phys. |

4. | G. Veronis and S. Fan, “Guided subwavelength plasmonic mode supported by a slot in a thin metal film,” Opt. Lett. |

5. | L. Chen, J. Shakya, and M. Lipson, “Subwavelength confinement in an integrated metal slot waveguide on silicon,” Opt. Lett. |

6. | D. Gramotnev, “Adiabatic nanofocusing of plasmons by sharp metallic grooves: Geometrical optics approach,” J. Appl. Phys. |

7. | R. Buckley and P. Berini, “Figures of merit for 2D surface plasmon waveguides and application to metal stripes,” Opt. Express |

8. | P. Berini, “Plasmon-polariton waves guided by thin lossy metal films of finite width: Bound modes of symmetric structures,” Phys. Rev. B |

9. | R. Zia, J. A. Schuller, A. Chandran, and M. L. Brongersma, “Plasmonics: the next chip-scale technology,” Mat. Today |

10. | S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, and T. W. Ebbesen, “Channel plasmon-polariton guiding by subwavelength metal grooves,” Phys. Rev. Lett. |

11. | M. Yan and M. Qiu, “Guided plasmon polariton at 2D metal corners,” J. Opt. Soc. Am. B |

12. | G. Veronis and S. Fan, “Modes of subwavelength plasmonic slot waveguides,” IEEE J. Lightwave Tech. |

13. | G. Veronis and S. Fan, “Crosstalk between three-dimensional plasmonic slot waveguides,” Opt. Express |

14. | H. S. Won, K. C. Kim, S. H. Song, C. Oh, P. S. Kim, S. Park, and S. Kim, “Vertical coupling of long-range surface plasmon polaritons,” Appl. Phys. Lett. |

15. | J. A. Pereda, A. Vegas, and A. Prieto, “An improved compact 2D full-wave FDFD method for general guided wave structures,” Microwave Opt. Technol. Lett. |

16. | |

17. | H. J. Hagemann, W. Gudat, and C. Kunz, “Optical constants from the far infrared to the x-ray region: Mg, Al, Cu, Ag, Au, Bi, C, and Al |

18. | R. M. Reano and S. W. Pang, “Sealed three-dimensional nanochannels,” J. Vac. Sci. Technol. B |

19. | Bruce McConnel, “Using self-assembly to create airgap microprocessors” (IBM, 2007). http://www-03.ibm.com/press/us/en/presskit/21463.wss. Accessed 10-30-07. |

20. | D. F. P. Pile, D. K. Gramotnev, M. Haraguchi, T. Okamoto, and M. Fukui, “Numerical analysis of coupled wedge plasmons in a structure of two metal wedges separated by a gap,” J. Appl. Phys. |

21. | K. Okamoto, |

22. | C. Pollock, |

23. | W. C. Chew, |

24. | D. Marcuse, |

25. | G. Strang, |

26. | R. Zia, M. D. Selker, P. B. Catrysse, and M. L. Brongersma, “Geometries and materials for subwavelength surface plasmon modes,” J. Opt. Soc. Am. A. |

27. | S. Collin, F. Pardo, and J. Pelouard, “Waveguiding in nanoscale metallic apertures,” Opt. Express |

28. |

**OCIS Codes**

(130.1750) Integrated optics : Components

(130.2790) Integrated optics : Guided waves

(130.3120) Integrated optics : Integrated optics devices

(250.5403) Optoelectronics : Plasmonics

**ToC Category:**

Integrated Optics

**History**

Original Manuscript: July 14, 2008

Revised Manuscript: August 4, 2008

Manuscript Accepted: August 5, 2008

Published: August 6, 2008

**Citation**

Galen B. Hoffman and Ronald M. Reano, "Vertical coupling between gap plasmon waveguides," Opt. Express **16**, 12677-12687 (2008)

http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-17-12677

Sort: Year | Journal | Reset

### References

- J. A. Conway, S. Sahni, and T. Szkopek, "Plasmonic interconnects versus conventional interconnects: a comparison of latency, cross-talk and energy costs," Opt. Express 15, 4474 - 4484 (2007). [CrossRef] [PubMed]
- K. Jung, F. L. Teixeira, and R. M. Reano, "Au/SiO2 nanoring plasmon waveguides at optical communication band," IEEE J. Lightwave Tech. 25, 2757 - 2765 (2007). [CrossRef]
- P. Berini, "Air gaps in metal stripe waveguides supporting long-range surface plasmon polaritons," J. Appl. Phys. 102, 33112 (2007). [CrossRef]
- G. Veronis and S. Fan, "Guided subwavelength plasmonic mode supported by a slot in a thin metal film," Opt. Lett. 30, 3359 - 3361 (2005). [CrossRef]
- L. Chen, J. Shakya, and M. Lipson, "Subwavelength confinement in an integrated metal slot waveguide on silicon," Opt. Lett. 31, 2133 - 2135 (2006). [CrossRef] [PubMed]
- D. Gramotnev, "Adiabatic nanofocusing of plasmons by sharp metallic grooves: Geometrical optics approach," J. Appl. Phys. 98, 104302 (2005). [CrossRef]
- R. Buckley and P. Berini, "Figures of merit for 2D surface plasmon waveguides and application to metal stripes," Opt. Express 15, 12174 - 12182 (2007). [CrossRef] [PubMed]
- P. Berini, "Plasmon-polariton waves guided by thin lossy metal films of finite width: Bound modes of symmetric structures," Phys. Rev. B 61, 10484 - 10503 (2000). [CrossRef]
- R. Zia, J. A. Schuller, A. Chandran, and M. L. Brongersma, "Plasmonics: the next chip-scale technology," Mat. Today 9, 20 - 27 (2006). [CrossRef]
- S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, and T. W. Ebbesen, "Channel plasmon-polariton guiding by subwavelength metal grooves," Phys. Rev. Lett. 95, 46802 (2005). [CrossRef]
- M. Yan and M. Qiu, "Guided plasmon polariton at 2D metal corners," J. Opt. Soc. Am. B 24, 2333 - 2342 (2007). [CrossRef]
- G. Veronis and S. Fan, "Modes of subwavelength plasmonic slot waveguides," IEEE J. Lightwave Tech. 25, 2511 - 2521 (2007). [CrossRef]
- G. Veronis and S. Fan, "Crosstalk between three-dimensional plasmonic slot waveguides," Opt. Express 16, 2129 - 2140 (2008). [CrossRef] [PubMed]
- H. S. Won, K. C. Kim, S. H. Song, C. Oh, P. S. Kim, S. Park, and S. Kim, "Vertical coupling of long-range surface plasmon polaritons," Appl. Phys. Lett. 88, 11110 (2006). [CrossRef]
- J. A. Pereda, A. Vegas, and A. Prieto, "An improved compact 2D full-wave FDFD method for general guided wave structures," Microwave Opt. Technol. Lett. 38, 331 - 335 (2003). [CrossRef]
- http://www.comsol.com
- H. J. Hagemann, W. Gudat, and C. Kunz, "Optical constants from the far infrared to the x-ray region: Mg, Al, Cu, Ag, Au, Bi, C, and Al2O3," J. Opt. Soc. Am. 65, 742-744 (1975). [CrossRef]
- R. M. Reano and S. W. Pang, "Sealed three-dimensional nanochannels," J. Vac. Sci. Technol. B 23, 2995 - 2999 (2005). [CrossRef]
- Bruce McConnel, "Using self-assembly to create airgap microprocessors" (IBM, 2007), http://www-03.ibm.com/press/us/en/presskit/21463.wss, Accessed 10-30-07.
- D. F. P. Pile, D. K. Gramotnev, M. Haraguchi, T. Okamoto, and M. Fukui, "Numerical analysis of coupled wedge plasmons in a structure of two metal wedges separated by a gap," J. Appl. Phys. 100, 13101 (2006). [CrossRef]
- K. Okamoto, Fundamentals of Optical Waveguides, (Academic Press, 2006).
- C. Pollock, Fundamentals of Optoelectronics, (Ceramic Book and Literature Service, 2003), Chap. 11.
- W. C. Chew, Waves and Fields in Inhomogeneous Media, (IEEE Press, 1995), Chap. 6.
- D. Marcuse, Light Transmission Optics, (Van Norstrand Reinhold Company, 1972), pp. 322-326.
- G. Strang, Linear Algebra and Its Applications 3rd Ed., (Thomson Learning, 1988), pp. 448 - 449.
- R. Zia, M. D. Selker, P. B. Catrysse, and M. L. Brongersma, "Geometries and materials for subwavelength surface plasmon modes," J. Opt. Soc. Am. A. 21, 2442 - 2446 (2004). [CrossRef]
- S. Collin, F. Pardo, and J. Pelouard, "Waveguiding in nanoscale metallic apertures," Opt. Express 15, 4310 - 4320 (2007). [CrossRef] [PubMed]
- http://www.rsoftdesign.com

## Cited By |
Alert me when this paper is cited |

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

« Previous Article | Next Article »

OSA is a member of CrossRef.