# Design of Matched Absorbing Layers for Surface Plasmon-Polaritons

Authors: Sergio de la Cruz; Eugenio R. Méndez; Alexei A. Maradudin
Journal: Advances in OptoElectronics
 (2012)

---

## Abstract

We describe a procedure for designing metal-metal boundaries for the strong attenuation of surface plasmon-polaritons without the introduction of reflections or scattering effects. Solutions associated with different sets of matching materials are found. To illustrate the results and the consequences of adopting different solutions, we present calculations based on an integral equation formulation for the scattering problem and the use of a nonlocal impedance boundary condition.


---

## Body

## 1. Introduction

In the numerical solution of spatially unbounded electromagnetic problems it is often necessary to truncate the computational domain. The use of a region of finite spatial size can introduce reflections and other spurious effects in the calculations. Several techniques have been proposed to overcome this problem. An established one, used commonly in finite-difference time-domain (FDTD) calculations [1], is the perfectly matched layer (PML) technique proposed by Berenger [2]. Such a layer can absorb electromagnetic waves without reflections at the vacuum-PML interfaces. Methods of calculation based on Green’s theorem (see, e.g., [3]) do not present such complications with volume waves, but the truncation of the interfaces can produce spurious effects in the presence of surface waves, like surface plasmon-polaritons (SPPs).The basic properties of SPPs have been known for some time, but their importance for nanophotonic applications has produced a renewed interest on the subject [4, 5]. In studies of the interactions of SPPs with objects or surface structures, the computational problem grows as a function of the physical size of the sample and, thus, it is desirable to reduce the computational domain as much as possible. The PML techniques known to us were not designed to handle truncation effects involving SPPs, and are not well-adapted for situations involving metallic structures and evanescent waves [6].In this paper, we present a procedure for determining the optical constants of absorbing materials for the attenuation of SPPs without introducing, or minimizing at least, spurious reflections and/or radiative scattering effects. Although SPPs are already lossy traveling waves, a reduction of their propagation length in the matched medium permits an important reduction in the dimensions of the region over which the computational domain extends. Although our approach is related in spirit to the usual PML techniques, we point out that it addresses a different problem, namely, the termination of surfaces over which surface waves propagate without the introduction of artificial scattering or reflection effects.The absorbing materials that we propose are homogeneous and isotropic. The flexibility for choosing their optical properties permits the selection of materials that are appropriate for the use of impedance boundary conditions. All these facts simplify considerably the electromagnetic problem. The results presented here can be applied not only for electromagnetic calculations based on integral equations, like the one used for our examples, but also for those based on differential equations, like the FDTD.


## 2. Theory

To describe the technique, let us consider the generic geometry illustrated in Figure1. The figure shows a flat interface involving three media: a semi-infinite dielectric and two metals, defined through their permittivities and permeabilities, with a vertical boundary. We assume that SPPs traveling from left to right, impinge on the vertical boundary between media 1 and 2. The optical properties of medium 2 are to be chosen in such a way that the SPPs are strongly attenuated without the introduction of significant reflections or coupling into radiative modes.Figure 1
Schematic diagram of the three media considered. The upper medium is a dielectric, characterized by its electric permittivityϵd and magnetic permeability μd in the frequency region of interest. Similarly, the metals 1 and 2 are characterized by ϵ1, μ1 and ϵ2, μ2, respectively.Let us consider the magnetic field associated with an SPP traveling on the flat surface of medium 2, in the direction+x1. This field may be conveniently written in terms of the magnetic field along the x2-direction. We write(1a)H2>(x1,x3)=H0eiksp(2)x1+iαd(ksp(2))x3x3≥0,(1b)H2<(x1,x3)=H0eiksp(2)x1-iα2(ksp(2))x3x3≤0,where H0 is a constant,
(2)ksp(2)=ωcϵdϵ2(ϵ2μd-ϵdμ2)ϵ22-ϵd2,αd(ksp(2))=ϵdμd(ωc)2-(ksp(2))2,α2(ksp(2))=ϵ2μ2(ωc)2-(ksp(2))2.
The signs of the square roots are taken in such a way that the imaginary parts of ksp(2), αd(ksp(2)), and α2(ksp(2)), are positive. In the dielectric, the attenuation constant along x3 is given by
(3)κ2>=ℑm{αd(ksp(2))}=ωcℑm{ϵd2(ϵ2μ2-ϵdμd)ϵ22-ϵd2}.
Similar expressions can be written for the interface between the dielectric and medium 1.To minimize the scattering losses due to the vertical metallic interface of Figure1, the transverse profile of the SPP wave must be the same on the two sides of the interface [7–9]. That is, one must have that κ1>=κ2>. On the other hand, the reflection of the SPP at the boundary should be minimized when ℜe{ksp(2)}=ℜe{ksp(1)}. We will refer to these conditions as the no scattering and no reflection conditions.At this stage, we focus our attention on the case in whichϵd=1, μd=1, and assume that the metal on the left is nonmagnetic at optical frequencies (i.e., μ1=1), while that on the right has arbitrary properties. As a particular example, we consider the design of absorbing layers that match the properties of gold at the vacuum wavelength λ=980nm. Then, ϵ1=-40.44+i2.97 [10] and μ1=1.0. Since ϵ2 and μ2 are complex quantities, one must seek solutions in a four-dimensional space. To simplify matters, we have decided to fix the imaginary parts of ϵ2 and μ2, and explore the space (ℜe{ϵ2}-ℜe{μ2}).For physical reasons, we chooseℑm{ϵ2}>0. Although there seems to be no problem with the choice of a negative ℑm{μ2} [11], we also took it as positive. With these choices, the search space is restricted to the quadrant defined by the conditions ℜe{ϵ2}<0 and ℜe{μ2}>0.For our first example we have setℑm{ϵ2}=3.0 and ℑm{μ2}=0.0. The solutions to the equations for no scattering and no reflection are shown in Figure 2(a). We observe that the two curves practically coincide when ℜe{ϵ2}<-10 and that they gradually separate as ℜe{ϵ2} approaches zero. It is thus advisable to seek solutions in the region ℜe{ϵ2}≤-10.(a) Solution curves corresponding to the no reflection and no scattering conditions assumingℑm{ϵ2}=3.0 and ℑm{μ2}=0.0. (b) Behavior of ℑm{ksp(2)} along the solution curves. The value of ℜe{ϵ2} chosen for the example is denoted by the vertical dashed line.
(a)
(b)One can see in Figure2(a) that the pair of x,y values (-40.44,1), which corresponds to the properties of medium 1, lie on the solution curve. This is due to the closeness between the chosen value of ℑm{ϵ2} and that of ℑm{ϵ1}.As one moves along these solution curves from left to right, the absorption coefficient of the SPPs,ℑm{ksp(2)}, increases. The behavior is illustrated in Figure 2(b). On the basis of these results, and since we are looking for solutions with a strong attenuation coefficient, we have chosen a value of ℜe{ϵ2}=-10.0, which leads to ℜe{μ2}=0.19 and an attenuation constant ℑm{ksp(2)}=0.033μm-1. We will refer to this example as case I.Values of the real parts ofϵ2 and μ2 that result in higher attenuation coefficients are not appropriate for our purposes, as they produce a mismatch between the no scattering and no reflection conditions. Moreover, as can be inferred from the discussion in the following section, they lead to solutions that are more difficult to handle computationally.For our second example, we choseℑm{ϵ2}=3.0 and ℑm{μ2}=10.0. The solution map is shown in Figure 3(a), and the attenuation coefficient for SPPs is shown in Figure 3(b). We see that, with these parameters, the curves corresponding to the no scattering and no reflection conditions are fairly close, but differ by an increasing amount as the magnitude of ℜe{ϵ2} decreases. The attenuation coefficient curves, on the other hand, are practically indistinguishable in this case. For this second example, we choose ℜe{ϵ2}=-531.62, and based on the no-scattering-condition curve, one has that ℜe{μ2}=13.44 and ℑm{ksp(2)}=0.06μm-1. We will refer to this example as case II. It is worth mentioning that the results obtained with these parameters are practically the same as those obtained with the parameters corresponding to the no reflection curve.(a) Solution curves corresponding to the no reflection and no scattering conditions assumingℑm{ϵ2}=3.0 and ℑm{μ2}=10.0. (b) Behavior of ℑm{ksp(2)} along the solution curves. The value of ℜe{ϵ2} chosen for the example is denoted by the vertical dashed line.
(a)
(b)Scattering calculations involving gold and the materials proposed in this section will be presented and discussed in Section4.


## 3. A Nonlocal Impedance Boundary Condition

The attenuation of SPPs in the proposed matching media can be visualized through rigorous electromagnetic simulations of the problem. We base our calculations on an integral equation formulation that has been used for rough-surface-scattering studies [3]. Geometries like the one depicted in Figure 1 pose some difficulties for this kind of approach, but the problem can be simplified with the use of an impedance boundary condition [12–14]. The use of an impedance boundary condition in the present problem leads to a formulation of the scattering problem that does not require knowledge of the field below the interface.Impedance boundary conditions have been used in the past for scattering calculations and are known to provide accurate results for good conductors, like gold and silver, in the near infrared [15]. We shall see however that, due to the rapid decay of the SPPs in the absorbing region, a local impedance boundary condition is insufficient to deal with some of the media considered here. So, before discussing the scattering calculations, we present a brief derivation of the impedance boundary condition for flat surfaces, which includes the first nonlocal term.With reference to Figure1, we consider the field below a flat interface, inside a medium with the properties of medium 2. It can be written in the form
(4)H2<(x1,x3)=∫-∞∞‍dq2πA(q)eiqx1-iα2(q)x3,
where A(q) is the scattering amplitude, or angular spectrum of the field below the interface.Here,α2(q)=ϵ2μ2(w/c)2-q2, with ℑm{α2(q)}>0.Since(5)A(q)=∫-∞∞‍dx1′H2<(x1′,0)e-iqx1′,
one can write
(6)H2<(x1,x3)=∫-∞∞‍dq2π∫-∞∞‍dx1′×H2<(x1′,0)eiq(x1-x1′)-iα2(q)x3.
From this expression, one can establish the following integral equation involving the field and its normal derivative on the interface:
(7)∂H2<(x1,x3)∂x3|x3=0=∫-∞∞‍dx1′K(x1,x1′)H2<(x1′,0),
where
(8)K(x1,x1′)=1i∫-∞∞‍dq2πα2(q)eiq(x1-x1′).We now seek the solution of (7) by expanding α2(q) as follows:
(9)α2(q)≈ϵ2μ2(ωc)[1-q22ϵ2μ2(ω/c)2+⋯],
which leads to
(10)K(x1,x1′)=1iϵ2μ2(ωc)δ(x1-x1′)+12i(ω/c)ϵ2μ2δ′′(x1-x1′)+⋯.
In this expression, δ(x1) is a delta function and δ′′(x1) its second derivative.Keeping only the first two terms, we find the relation(11)∂H2<(x1,x3)∂x3|x3=0=1iϵ2μ2(ωc)H2<(x1,x3)|x3=0+12iϵ2μ2(ω/c)∂2H2<(x1,x3)∂x12|x3=0,
and using the continuity of the tangential components of the fields across the interface we find that
(12)1ϵdL>(x1)=1ϵ2[12iϵ2μ2(ω/c)d2H>(x1)dx121iϵ2μ2(ωc)H>(x1)1ε2+12iϵ2μ2(ω/c)d2H>(x1)dx12],
where
(13)H>(x1)=H2>(x1,x3)|x3=0,L>(x1)=∂H2>(x1,x3)∂x3|x3=0.
Finally, we write (12) in the form
(14)L>(x1)=Kp(0)H>(x1)+Kp(2)d2H>(x1)dx12,
where
(15)Kp(0)=ωicϵdμ2ϵ2,Kp(2)=ϵd/ϵ22iϵ2μ2(ω/c).
The first term in the impedance expansion (14) represents a local relation between L>(x1) and H>(x1) (or between the E and H fields). The second term, represents the first nonlocal correction.We see that as the attenuation of the SPPs becomes stronger, the second derivative of the field appearing in (14) becomes larger. Looking at the quantities entering the expression for Kp(2), we conclude that it is desirable to have media with a high refractive index n2=ϵ2μ2. For a more quantitative evaluation of this issue, we consider the ratio between the first two terms of the expansion. This ratio ℛ constitutes an incomplete but, nevertheless, essential knowledge for establishing the validity of the local relation and the convergence of the expansion. Using the general form of the field associated with an SPP on a flat surface given by (1a), we find that
(16)ℛ=|Kp(2)(d2H>(x1)/dx12)Kp(0)H>(x1)|=12|ksp(2)n2(ω/c)|2,
where n2 is the complex refractive index of medium 2.In normal circumstances (i.e., when μ2=1), for a good conductor one has that n2 is large and that |ksp(2)|≈(ω/c). Thus, ℛ is small, and the nonlocal term can be neglected. For instance, for a gold-vacuum interface and the wavelength considered, this ratio is ℛ=0.0125.Turning our attention to the absorbing materials proposed in Section2, for our first example (case I) this ratio turns out to be ℛ=0.258. It is clear that a local impedance boundary condition would not be accurate on the absorbing side of the boundary and that nonlocal correction terms are needed. This conclusion is supported by the numerical calculations presented in Section 4.For case II, the ratio between the first two terms of expansion (14) is ℛ=1.67×10-5. Not surprisingly, the local impedance boundary condition is very accurate in this case. This case leads to parameters that are more convenient for our type of calculation.


## 4. Results and Discussion

In this section, we present scattering calculations corresponding to the two examples of matching media proposed in Section2, starting with a summary of their properties: For case I we have ϵ2=-10.0+3.0i and μ2=0.19, while for case II, ϵ2=-531.62+3.0i and μ2=13.44+10.0i. In the first case, the attenuation constant for SPPs is ℑm{ksp(2)}=0.033μm-1, while for the second ℑm{ksp(2)}=0.06μm-1.The geometry employed for the scattering calculations is shown in Figure4. Based on previous results [16], to excite SPPs traveling to the right, we use an array of 5 rectangular grooves of period T=0.863λ, illuminated by a Gaussian beam of width g=2λ with an angle of incidence θ0≈-6.5∘. The coupling efficiency for the SPPs travelling to the right is about 45%. The distance between the start of the buffer material and the end of the surface is 45μm.Figure 4
Schematic diagram of the scattering geometry considered for the calculations. A Gaussian beam illuminates a short gold grating that couples a fraction of the incident light to SPPs traveling to the right. The absorbing material is represented by the darker region.The magnitude of the surface field associated with the SPP propagating along the surface for case I is shown in Figure5. To simplify the visualization of the region of interest, the section with the grating was omitted from the figure. The curve shown with a continuous line was calculated with a local impedance boundary condition, while for the dashed-line curve we included also the first nonlocal term. One observes that the attenuation constant of SPP changes as the wave enters the absorbing material. As expected, the nonlocal correction does make a substantial difference in the calculations. The small ripples observed in the upper curve are due to standing waves caused by the reflection of the SPP at the end of the surface. With this material it is important to use at least the first nonlocal correction in expansion (14).Figure 5
Magnitude of the surface field calculated using the local impedance boundary condition and the first nonlocal correction for case I. The boundary between the two metals is atx1=0 and, to facilitate the visualization of the interaction of the SPP with the metal-metal boundary, the region of the grating is not shown.In Figure6, we show the calculated magnetic near-field intensity map (i.e., |H|2) for case I. As in the previous figure, only the rightmost section of the sample is shown. It can be observed that the SPP arrives from the left and decays quite steeply as it enters the absorbing medium. The results are encouraging. The incident SPP is not perturbed much by the interface, and it is clear that the reflection and scattering effects are low.Figure 6
Squared modulus of the magnetic field in the vacuum region, above the flat metallic surfaces for case I. The position of the vertical border between the two metals atx1=0 is denoted by the vertical dashed white line.The near-field intensity map corresponding to case II is shown in Figure7. The attenuation of the SPP is stronger in this case, and the local impedance boundary condition gives accurate results. As desired, the SPP wave practically vanishes after only a few microns. As in case I, the incident SPP does not appear to be perturbed much by the interface.Figure 7
Squared modulus of the magnetic field in the vacuum region, above the flat metallic surfaces for case II. The position of the vertical border between the two metals atx1=0 is denoted by the vertical dashed white line.To investigate the possible leakage of the SPPs guided waves due to the absorbing boundary, we present, in Figure8, far-field calculations corresponding to a complete gold sample and to samples terminated with the absorbing materials that we have called cases I and II (see Figure 4).The differential reflection coefficient for a gold surface (a), a gold surface terminated with the material of case I (b), and a gold surface terminated with the material of case II (c).
(a)
(b)
(c)The scattering curves show that there is a relatively strong reflected beam in the specular direction (θs≈-6.5∘) as well as a broad, grazing diffraction order. The rapid oscillations observed in the diffraction order for the case of the gold surface are due to the leakage of the SPP at the end of the surface, which interfere with the grazing order. Moreover, the SPP reflected at this end propagate back and leak into the specular direction upon interaction with the grating. This gives rise to the small-scale oscillations over the whole curve.The differential reflection coefficient curves corresponding to cases I and II are much cleaner. In particular, the fine-scale oscillations that can be observed over the whole curve associated with the gold surface have disappeared. Also, the strong and rapid oscillations in the region of the diffraction order have been replaced by smoother and less important oscillations that seem to be due to residual scattering at the metal-metal interface.


## 5. Summary and Conclusions

In this paper, we have presented a procedure for designing metal-metal boundaries for the strong attenuation of surface plasmon-polaritons, minimizing the introduction of reflections or scattering effects.We have illustrated the results by means of computer simulations based on an integral equation formulation of the scattering problem and an impedance boundary condition. Depending on the parameters chosen for the matching medium, nonlocal corrections might be needed for the calculations.The results show that with the method described here, SPPs can be attenuated strongly without introducing significant effects in near field calculations. The propagation of the SPP is practically unperturbed along the original interface and spurious reflection and scattering effects are minimal. The implementation of the technique is fairly simple and can be used effectively in plasmonic and nanophotonic calculations by different methods.The far-field results, however, seem to indicate the presence of some residual scattering originating from the metal-metal boundary. It is worth pointing out that this leakage is not due to the use of an impedance boundary condition or to the small difference between the values of the dielectric constants corresponding to the no-scattering and no-reflection conditions.



---

*Source: 598213-2012-08-29.xml*