# Analytic Modeling of Trench Waveguides for Channel Plasmon-Polariton Using Full-Hybrid Trial-Field Functions

Authors: Sujit Chattopadhyay; Pradip Kumar Saha
Journal: ISRN Optics
 (2012)

---

## Abstract

A full-hybrid trial-field modeling, which was proposed earlier by the authors for plasmonic structures of rectangular cross-section, has been applied to trench waveguides for channel plasmon polariton. Most of the computed results of complex propagation constant agree well with those obtained by method of lines and effective index method, though the latter two results themselves show some discrepancy. The particular aspect in which there is varying degree of disagreement among the results obtained in the above three   methods is the imaginary part of the effective index seen by the plasmon-polariton.


---

## Body

## 1. Introduction

In the field of nanophotonics, the ability of surface plasmon polariton (SPP) to beat the diffraction limit is a tremendous boon to technology and its usefulness towards achieving higher integration density in photonic integrated circuits is now well established. As such, well-confined wave guidance and improved propagation length are the key aspects being addressed by the researchers. These goals in turn have led to a number of innovations in the SPP guiding structures. Reviews of SPP guides can be found, for example, in [1–3]. Apart from the basic SPP guiding structure of metal strip embedded in homogeneous or inhomogeneous dielectrics and its variants, gap-SPP and channel-SPP have been shown to have several desirable features [4–7]. Channel SPPs, commonly referred to as channel plasmon polariton or CPP, are guided by two basic structures. (a) The channel is a V-shaped groove in a metal. The supported CPP is generally referred to as V-groove channel plasmon polariton or V-groove CPP. (b) The channel is a narrow rectangular groove in a metal. The supported CPP is referred to as trench channel plasmon polariton or trench-CPP. Simultaneous availability of subwavelength lateral confinement and relatively long propagation length is the most desirable feature of CPP. Photonic components based on both types of CPPs have been demonstrated in practice [8–10]. Together with different simulation techniques, the effective index method (EIM), traditionally employed in dielectric waveguide problems, has been used for determining the propagation characteristics of channel-SPP guides and in modeling SP waveguide components [1, 11, 12]. The main attractive feature of EIM is that the characteristics of 2D guides can be predicted by combining results obtained from analysis of 1D guides. Another method, the semianalytic method of lines, has proved to be a very useful tool for numerical modeling of plasmonic waveguides [13–16]. Minh and others [17] recently have applied MoL to calculate the propagation characteristics of trench-CPP waveguides and have compared their results with those obtained by EIM [12] and FDTD [18]. The present authors have carried out a simple analytic modeling of trench-CPP guides using full-hybrid trial-field functions. This approximate yet quite accurate technique was earlier successfully applied to metallic SPP guides of rectangular cross-section and to rectangular nanoholes in metals [19]. The computed values of complex effective index and propagation length for various trench dimensions and dielectric filling are compared with those in [12, 17].


## 2. Formulation of the Problem

The geometry of the trench SPP guide is shown in Figure1(a). It is a rectangular groove of depth d and width 2w in a semi-infinite metal block. It can be bisected longitudinally by a symmetry plane at y=0. The groove is usually air-filled but can be filled with other dielectrics. All modal fields are assumed to have exp(-jβz) dependence in the direction of propagation.(a) Geometry of trench-CPP waveguide. (b) One half of trench-CPP guide segmented into six regions.
(a)(b)The parameters are chosen as in [12, 17]: the wavelengths of operation λ=1550 nm, 1033 nm, 775 nm; the metal is gold with dielectric constant of −131.95–j 12.65 at 1550 nm; trench width 2w=300 nm, 500 nm; depth d is variable.


### 2.1. Trial Field Modeling

In this approach, we would construct the trial transverse spatial functions forEyTE and HyTM fields corresponding to TEx (Ex=0) and TMx  (Hx=0) modes, respectively, and derive all six field components of the hybrid plasmon modes of the structure from superposition of the TEx and TMx mode families [13]. While writing the trial field functions, first we consider the following points: (a) only one-half of the structure, shown in Figure 1(b), that takes into account the symmetry plane at y=0 is used for formulation, (b) the reduced structure is segmented into six regions, R1 to R6, as shown in Figure 1(b), (c) trial fields for each region i (i=1 to 6) are constructed separately such that each trial field Ψi satisfies the Helmholtz equation in respective region i, and (d) continuity of the transverse wave-vector at any interface between two media is also assumed.  Secondly, since the success of trial-field modeling depends on how well the trial functions approximate the true fields, in writing the trial field functions inside the trench (region R1), the nature of field variation, particularly with depth, as reported in [12, 17] was taken into consideration. Clearly, the symmetry plane at y=0 is electric for the dominant field configuration. Provision for magnetic symmetry, however, is retained in the formulation. For x-dependence of the field inside the trench, with respect to the chosen coordinates, one has the option of writing the functional form as sinhkxx or coshkxx, both growing toward the top of the trench and being almost identical beyond the argument value of about two. Choice of sinh function forces the main electric field Ey to be zero at the trench bottom irrespective of the trench depth. The reported variation [12, 17], however, shows that the field is always very small (but not zero) at the bottom, and this value depends on the depth d. For short depths, the field at the bottom is not negligibly small and penetration into the metal is also significant. For deep trenches, both are very small. In fact, our computation with sinh function in the trial field for different trench widths always yielded the propagation constant of the gap-SPP mode [12]. Thus cosh dependence on x inside the trench appears to be the logical choice for the trial field function. On inspection and with the above considerations in mind, the trial-field functions for different regions are written as follows. Although regions R5 and R6 are shown in the segmented structure, the trial fields in these two regions are not written. Obviously, field matching at the R2–R6 and R5-R6 interfaces are not meaningful. Similar observations can also be made about field matching at R1–R3 and R3-R4 interfaces. Furthermore, field matching at the trench bottom (R1–R5 interface) is not necessary because it does not yield any useful information. This is to be expected from the reported variation of field with depth [12], which shows that this interface is not plasmonic; that is, it does not support a surface electromagnetic wave. So finally we write the trial field functions only for three regions R1, R2, and R4:(1)E1yTE=∑n=1∞A1ncosh⁡⁡(kn1xx)(kn1yy)cosh⁡⁡sinh⁡⁡,H1yTM=∑n=1∞B1ncosh⁡⁡(kn1xx)(kn1yy)sinh⁡⁡cosh⁡⁡,E2yTE=∑n=1∞A2ncosh⁡⁡(kn1xx)e-kn2y(y-w),H2yTM=∑n=1∞B2ncosh⁡⁡(kn1xx)e-kn2y(y-w),E4yTE=∑n=1∞A4ne-kn4x(x-d)e-kn2y(y-w),H4yTM=∑n=1∞B4ne-kn4x(x-d)e-kn2y(y-w).The transverse wave numbers are related toβ by the following equation:(2)knix2+kniy2=β2-k02εi.The upper and lower functions correspond, respectively, to magnetic and electric walls at the symmetry planey=0. The expressions (1) for the field components need some clarifications. The summation on n in expressions (1) is for the purpose of writing the trial functions in generalized form and do not represent expansion in terms of orthogonal basis functions. Only the nature of a physical solution for a prospective plasmon mode is envisaged while constructing the trial functions. However, the eigenvalue equations are transcendental in nature and can be solved for a series of complex roots (knix, kniy) which would give the same value of β for a particular mode and different β for different modes. Therefore, it is sufficient to compute, say, k11x and k11y, for computation of β of the lowest mode. However, the root set (knix, kniy) can be computed to obtain the field profiles in different regions.


### 2.2. Eigenvalue Equations

Once all six field components in each region (R1, R2, R4) are derived from the trial fields, the next task would be to apply the continuity conditions. It is not necessary to impose all the boundary conditions at all the interfaces which only overdefines the problem. We found that application of continuity of tangential field components at only the R1-R2 and R2–R4 interfaces (between two different media) yields the desired eigenvalue equations, given by (3), for all values of n:(3)Det[cosh⁡⁡(kn1xd)-kn1xkn2yωβε0ε2sinh⁡⁡(kn1xd)-1-kn4xkn2yωβε0ε4kn1xkn2yωβμ0sinh⁡⁡(kn1xd)cosh⁡⁡(kn1xd)kn4xkn2yωβμ0-1kn2yβcosh⁡⁡(kn1xd)-kn1xωε0ε2sinh⁡⁡(kn1xd)-kn2yβ-kn4xωε0ε4kn1xωμ0sinh⁡⁡(kn1xd)kn2yβcosh⁡⁡(kn1xd)kn4xωμ0-kn2yβ]=0,X12coth⁡⁡tanh⁡(kn1yw)+kn1ykn2y=0,where X12=(kn1x2+k02ε1)/(kn1x2+k02ε2).


## 2.1. Trial Field Modeling

In this approach, we would construct the trial transverse spatial functions forEyTE and HyTM fields corresponding to TEx (Ex=0) and TMx  (Hx=0) modes, respectively, and derive all six field components of the hybrid plasmon modes of the structure from superposition of the TEx and TMx mode families [13]. While writing the trial field functions, first we consider the following points: (a) only one-half of the structure, shown in Figure 1(b), that takes into account the symmetry plane at y=0 is used for formulation, (b) the reduced structure is segmented into six regions, R1 to R6, as shown in Figure 1(b), (c) trial fields for each region i (i=1 to 6) are constructed separately such that each trial field Ψi satisfies the Helmholtz equation in respective region i, and (d) continuity of the transverse wave-vector at any interface between two media is also assumed.  Secondly, since the success of trial-field modeling depends on how well the trial functions approximate the true fields, in writing the trial field functions inside the trench (region R1), the nature of field variation, particularly with depth, as reported in [12, 17] was taken into consideration. Clearly, the symmetry plane at y=0 is electric for the dominant field configuration. Provision for magnetic symmetry, however, is retained in the formulation. For x-dependence of the field inside the trench, with respect to the chosen coordinates, one has the option of writing the functional form as sinhkxx or coshkxx, both growing toward the top of the trench and being almost identical beyond the argument value of about two. Choice of sinh function forces the main electric field Ey to be zero at the trench bottom irrespective of the trench depth. The reported variation [12, 17], however, shows that the field is always very small (but not zero) at the bottom, and this value depends on the depth d. For short depths, the field at the bottom is not negligibly small and penetration into the metal is also significant. For deep trenches, both are very small. In fact, our computation with sinh function in the trial field for different trench widths always yielded the propagation constant of the gap-SPP mode [12]. Thus cosh dependence on x inside the trench appears to be the logical choice for the trial field function. On inspection and with the above considerations in mind, the trial-field functions for different regions are written as follows. Although regions R5 and R6 are shown in the segmented structure, the trial fields in these two regions are not written. Obviously, field matching at the R2–R6 and R5-R6 interfaces are not meaningful. Similar observations can also be made about field matching at R1–R3 and R3-R4 interfaces. Furthermore, field matching at the trench bottom (R1–R5 interface) is not necessary because it does not yield any useful information. This is to be expected from the reported variation of field with depth [12], which shows that this interface is not plasmonic; that is, it does not support a surface electromagnetic wave. So finally we write the trial field functions only for three regions R1, R2, and R4:(1)E1yTE=∑n=1∞A1ncosh⁡⁡(kn1xx)(kn1yy)cosh⁡⁡sinh⁡⁡,H1yTM=∑n=1∞B1ncosh⁡⁡(kn1xx)(kn1yy)sinh⁡⁡cosh⁡⁡,E2yTE=∑n=1∞A2ncosh⁡⁡(kn1xx)e-kn2y(y-w),H2yTM=∑n=1∞B2ncosh⁡⁡(kn1xx)e-kn2y(y-w),E4yTE=∑n=1∞A4ne-kn4x(x-d)e-kn2y(y-w),H4yTM=∑n=1∞B4ne-kn4x(x-d)e-kn2y(y-w).The transverse wave numbers are related toβ by the following equation:(2)knix2+kniy2=β2-k02εi.The upper and lower functions correspond, respectively, to magnetic and electric walls at the symmetry planey=0. The expressions (1) for the field components need some clarifications. The summation on n in expressions (1) is for the purpose of writing the trial functions in generalized form and do not represent expansion in terms of orthogonal basis functions. Only the nature of a physical solution for a prospective plasmon mode is envisaged while constructing the trial functions. However, the eigenvalue equations are transcendental in nature and can be solved for a series of complex roots (knix, kniy) which would give the same value of β for a particular mode and different β for different modes. Therefore, it is sufficient to compute, say, k11x and k11y, for computation of β of the lowest mode. However, the root set (knix, kniy) can be computed to obtain the field profiles in different regions.


## 2.2. Eigenvalue Equations

Once all six field components in each region (R1, R2, R4) are derived from the trial fields, the next task would be to apply the continuity conditions. It is not necessary to impose all the boundary conditions at all the interfaces which only overdefines the problem. We found that application of continuity of tangential field components at only the R1-R2 and R2–R4 interfaces (between two different media) yields the desired eigenvalue equations, given by (3), for all values of n:(3)Det[cosh⁡⁡(kn1xd)-kn1xkn2yωβε0ε2sinh⁡⁡(kn1xd)-1-kn4xkn2yωβε0ε4kn1xkn2yωβμ0sinh⁡⁡(kn1xd)cosh⁡⁡(kn1xd)kn4xkn2yωβμ0-1kn2yβcosh⁡⁡(kn1xd)-kn1xωε0ε2sinh⁡⁡(kn1xd)-kn2yβ-kn4xωε0ε4kn1xωμ0sinh⁡⁡(kn1xd)kn2yβcosh⁡⁡(kn1xd)kn4xωμ0-kn2yβ]=0,X12coth⁡⁡tanh⁡(kn1yw)+kn1ykn2y=0,where X12=(kn1x2+k02ε1)/(kn1x2+k02ε2).


## 3. Results of Computation

With the chosen parameters as mentioned in Section2, complex effective index and propagation length from its imaginary part were computed and are presented in Figures 2(a) to 4(b).(a) Dispersion of the real part of the fundamental trench-CPP mode effective index with depth of the air-filled trench. (b) Dispersion of the propagation length of the fundamental trench-CPP mode with depth of the air-filled trench.
(a)(b)


### 3.1. Effect of Trench Depth

The variation ofRe(neff) with trench depth d, as shown in Figure 2(a), is very much similar to the EIM results [12], though our computed values are always slightly higher. The computed propagation length L (Figure 2(b)) for 2w=300 nm is almost same as that in [12]. For 2w=500 nm our values, however, are slightly lower, particularly at smaller depths. Though we have not shown comparison with MoL data in Figures 2(a) and 2(b); an examination of Figures  3(a) and 3(b) of [17] indicate that our data of Re(neff) agree better with the MoL data than with the EIM data. However, our computed values of propagation length agree much better with EIM data. The MoL value of propagation length hardly changes with trench depth and shows a small increase at small trench depth of about 1 μm. The discrepancy in the nature of variation of the MoL values of L with respect to EIM data for small trench depths has been noted in [17]. All the previously reported data and our results confirm that subwavelength lateral confinement as well as relatively long propagation length can be achieved with trench-CPP waveguides with deep grooves. Specifically, the fundamental trench-CPP mode which shows these features require large depth-to-width ratio.(a) Dispersion of the real part of the fundamental trench-CPP mode effective index with trench depth for different dielectric inserts inside the trench (n=3, 2, 1.44). (b) Dispersion of the imaginary part of the fundamental trench-CPP mode effective index with trench depth for different dielectric inserts inside the trench (n=3, 2, 1.44).
(a)(b)(a) Dispersion of the real part of the fundamental trench-CPP mode effective index for different wavelengths of operation with depth of the air-filled trench. (b) Dispersion of the imaginary part of the fundamental trench-CPP mode effective index for different wavelengths of operation with depth of the air-filled trench.
(a)(b)


### 3.2. Effect of Dielectric Filling of Trench

It has been shown [12, 17] that, in a trench waveguide with air-filled groove (which can be looked upon as two right-angled wedges hybridized by a finite-depth rectangular groove), the field is mainly confined at the two corners at the top of the trench. In [17] it was proposed that dielectric filling (ε1>1) of the groove can be utilized for stronger confinement of field inside the trench. Results of MoL computation have established the projected improvement in normalized phase constant. Unfortunately, this is accompanied by increase in normalized attenuation constant. However, there is still provision for tradeoff between the two characteristics obtained with dielectric filled trench.The trial-field modeling with segmented problem area is well suited for computing the complex effective index of dielectric-filled trench configuration. The results of our computation are presented in Figures3(a) and 3(b), where comparison with MoL data also is shown. The agreement of Re(neff) is excellent. However, there is discrepancy in the imaginary part particularly for the lower trench depths (<1.5 μm), which is just the opposite of the variation obtained from MoL for air-filled trench. In [17], Im(neff) increases with decreasing trench depth which is just the opposite of the variation obtained by MoL for air-filled trench. It is difficult to understand why dielectric filling should cause this reversal in the nature of variation. The authors in [17] also have not offered any plausible explanation. The results obtained from trial-field modeling indicate the Im(neff) decreases with decreasing trench depth for both air- and dielectric-filled grooves. The decrease, however, is not exactly monotonic when ε1>1. For small depths of both air- and dielectric-filled grooves, the confinement is poor, the loss is lower, and hence the propagation length would be longer.


### 3.3. Effect of Wavelength

Figures4(a) and 4(b) show the computed effective index for air-filled trench as function of depth at two other (smaller) wavelengths 775 nm and 1033 nm. These are compared with the MoL data [17]. The agreement of Re(neff) data is reasonably good except at small trench depths where our results decrease faster. Clearly, the wave confinement improves at lower wavelengths. Our computed values of Im(neff) at these wavelengths predict longer propagation lengths than the MoL results suggest. In fact our computed values of L at 1550 nm, which we have already mentioned to be in better agreement with the EIM data, also indicated higher propagation length than the MoL prediction.


## 3.1. Effect of Trench Depth

The variation ofRe(neff) with trench depth d, as shown in Figure 2(a), is very much similar to the EIM results [12], though our computed values are always slightly higher. The computed propagation length L (Figure 2(b)) for 2w=300 nm is almost same as that in [12]. For 2w=500 nm our values, however, are slightly lower, particularly at smaller depths. Though we have not shown comparison with MoL data in Figures 2(a) and 2(b); an examination of Figures  3(a) and 3(b) of [17] indicate that our data of Re(neff) agree better with the MoL data than with the EIM data. However, our computed values of propagation length agree much better with EIM data. The MoL value of propagation length hardly changes with trench depth and shows a small increase at small trench depth of about 1 μm. The discrepancy in the nature of variation of the MoL values of L with respect to EIM data for small trench depths has been noted in [17]. All the previously reported data and our results confirm that subwavelength lateral confinement as well as relatively long propagation length can be achieved with trench-CPP waveguides with deep grooves. Specifically, the fundamental trench-CPP mode which shows these features require large depth-to-width ratio.(a) Dispersion of the real part of the fundamental trench-CPP mode effective index with trench depth for different dielectric inserts inside the trench (n=3, 2, 1.44). (b) Dispersion of the imaginary part of the fundamental trench-CPP mode effective index with trench depth for different dielectric inserts inside the trench (n=3, 2, 1.44).
(a)(b)(a) Dispersion of the real part of the fundamental trench-CPP mode effective index for different wavelengths of operation with depth of the air-filled trench. (b) Dispersion of the imaginary part of the fundamental trench-CPP mode effective index for different wavelengths of operation with depth of the air-filled trench.
(a)(b)


## 3.2. Effect of Dielectric Filling of Trench

It has been shown [12, 17] that, in a trench waveguide with air-filled groove (which can be looked upon as two right-angled wedges hybridized by a finite-depth rectangular groove), the field is mainly confined at the two corners at the top of the trench. In [17] it was proposed that dielectric filling (ε1>1) of the groove can be utilized for stronger confinement of field inside the trench. Results of MoL computation have established the projected improvement in normalized phase constant. Unfortunately, this is accompanied by increase in normalized attenuation constant. However, there is still provision for tradeoff between the two characteristics obtained with dielectric filled trench.The trial-field modeling with segmented problem area is well suited for computing the complex effective index of dielectric-filled trench configuration. The results of our computation are presented in Figures3(a) and 3(b), where comparison with MoL data also is shown. The agreement of Re(neff) is excellent. However, there is discrepancy in the imaginary part particularly for the lower trench depths (<1.5 μm), which is just the opposite of the variation obtained from MoL for air-filled trench. In [17], Im(neff) increases with decreasing trench depth which is just the opposite of the variation obtained by MoL for air-filled trench. It is difficult to understand why dielectric filling should cause this reversal in the nature of variation. The authors in [17] also have not offered any plausible explanation. The results obtained from trial-field modeling indicate the Im(neff) decreases with decreasing trench depth for both air- and dielectric-filled grooves. The decrease, however, is not exactly monotonic when ε1>1. For small depths of both air- and dielectric-filled grooves, the confinement is poor, the loss is lower, and hence the propagation length would be longer.


## 3.3. Effect of Wavelength

Figures4(a) and 4(b) show the computed effective index for air-filled trench as function of depth at two other (smaller) wavelengths 775 nm and 1033 nm. These are compared with the MoL data [17]. The agreement of Re(neff) data is reasonably good except at small trench depths where our results decrease faster. Clearly, the wave confinement improves at lower wavelengths. Our computed values of Im(neff) at these wavelengths predict longer propagation lengths than the MoL results suggest. In fact our computed values of L at 1550 nm, which we have already mentioned to be in better agreement with the EIM data, also indicated higher propagation length than the MoL prediction.


## 4. Conclusions

A simple and approximate yet reasonably accurate method of modeling based on full-hybrid trial field functions has been applied to trench-CPP waveguide. The technique was earlier validated by successfully applying it to metal-strip SPP guide and nanoholes of rectangular cross-section in metals [19]. The literature on trench-CPP is very much limited at this time. The reported comparison of EIM and MoL results showed only partial agreement [12, 17]. And now, our results, obtained with trial-field modeling, agree well in parts with EIM and partly with MoL data. The problem of discrepancy appears to be mainly in the imaginary part of the effective index or the propagation length, particularly for small values of trench depth. As such any claim about one particular method being more accurate than the others may not be justified at this point of time. Perhaps more detailed study with more rigorous formulation will resolve this issue. Meanwhile we put forward our approach as a very much useful method.



---

*Source: 219281-2012-01-11.xml*