Open Access
24 September 2012 Excitation of multiple surface-plasmon-polariton waves using a compound surface-relief grating
Author Affiliations +
Abstract
The excitation of multiple surface-plasmon-polariton waves, all of the same frequency but different polarization states, phase speeds, spatial profiles and degrees of localization, by a compound surface-relief grating formed by a metal and a rugate filter, both of finite thickness, was studied using the rigorous coupled-wave approach. Each period of the compound surface-relief grating was chosen to have an integral number of periods of two different simple surface-relief gratings. The excitation of different SPP waves was inferred from the absorptance peaks that were independent of the thickness of the rugate filter. The excitation of each SPP wave could be attributed to either a simple surface-relief grating present in the compound surface-relief grating or to the compound surface-relief grating itself. However, the excitation of SPP waves was found to be less efficient with the compound surface-relief grating than with a simple surface-relief grating.

1.

Introduction

The solution of a canonical boundary-value problem demonstrates that the planar interface of (1) a semi-infinite metal, and (2) a semi-infinite homogeneous, isotropic, dielectric material can guide one surface plasmon-polariton (SPP) wave, at a specific frequency and in a specific direction lying in the interface plane.1,2 This SPP wave is of the p-polarization state. In practice, SPP waves are excited in the prism-coupled and the grating-coupled configurations. The prism-coupled configuration is widely used for chemical sensing.2 Whereas the grating-coupled configuration has been known for at least three decades to be useful in enhancing the absorptance of solar cells,3,4 this enhancement has only recently been interpreted in terms of SPP waves.5,6

In the grating-coupled configuration, the planar interface is replaced by one with periodic corrugations. This interface is a surface-relief grating. Light is then incident on the other face of the partnering dielectric material. The electromagnetic field phasors are expressed in terms of Floquet harmonics. If the period of the surface-relief grating is appropriate, a part of the incident light may couple to the SPP wave at a certain value of the incidence angle. This happens when one of the Floquet harmonics has a wavevector with a component parallel to the mean plane of the grating the same as the real part of the possible SPP wave.2,7

Recently, to enhance the possibility of this match occurring, Dolev et al.8 made the surface-relief grating a compound grating, called a quasiperiodic grating by them. Each period of the compound grating comprises several periods each of more than one simple surface-relief gratings. Dolev et al.8 showed that such a compound grating can be advantageously used to excite the same p-polarized SPP wave at different angles of incidence.

Multiple SPP waves of the same frequency and direction of propagation, but with different phase speeds, spatial profiles, and degrees of localizations, can be guided by the planar interface of a metal and a dielectric material that is periodically nonhomogeneous normal to the interface, whether that interface is either planar9 or a simple surface-relief grating.10 The multiplicity offers exciting prospects for enhancing the scope of applications of SPP waves.11,12 Chief among these are multianalyte chemical sensors13 and thin-film solar cells.14 The higher the number of SPP waves that can be excited, the more attractive are the prospects.

Although multiple SPP waves can be excited in the grating-coupled configuration,10 in line with predictions from the solution of the canonical boundary-value problem,9 a simple surface-relief grating may not suffice for the excitation of all possible SPP waves. With simple surface-relief gratings of different periods, different subsets of the possible SPP waves may be excited,10 which limits the benefits resulting from multiple SPP waves. Therefore, we decided to investigate whether a compound surface-relief grating would be better than a simple surface-relief grating.

Let us note the important difference between the work presented here and that of Dolev et al.:8 the partnering dielectric material in our work is periodically nonhomogeneous normal to the mean metal/dielectric interface, allowing the excitation of multiple SPP waves of the same frequency but different linear polarization states, phase speeds, and spatial profiles, thereby in contrast to Dolev et al.8 whose partnering dielectric material is homogeneous.

Because the theoretical formulation of the boundary-value problem discussed in this paper is the same as in Ref. 10, only a brief overview of the formulation is provided in Sec. 2. Illustrative numerical results are presented and discussed in Sec. 3, and the concluding remarks are presented in Sec. 4. An exp(iωt) time-dependence is implicit, with ω denoting the angular frequency. The free-space wavenumber, the free-space wavelength, and the intrinsic impedance of free space are denoted by k0=ωϵ0μ0, λ0=2π/k0, and η0=μ0/ϵ0, respectively, with μ0 and ϵ0 being the permeability and permittivity of free space. Vectors are in boldface, the asterisk denotes the complex conjugate, and the Cartesian unit vectors are identified as u^x, u^y, and u^z. The real part of a complex number ζ is denoted by Re(ζ).

2.

Boundary-Value Problem

Let us consider the boundary-value problem shown schematically in Fig. 1. The regions z<0 and z>d3 are vacuous. The region 0zd1 is occupied by a rugate filter15,16 with relative permittivity

Eq. (1)

ϵd(z)=[(nb+na2)+(nbna2)sin(πd2zΩ)]2,z>0,
where 2Ω is the period, and na and nb are the lowest and the highest indexes of refraction, respectively. Such filters can be deposited on top of a metallic film with a corrugated surface with a variety of well-established techniques.1719 The region d2zd3 is occupied by the metallic partnering material with spatially uniform relative permittivity ϵm.

Fig. 1

Schematic of the grating-coupled configuration used for the excitation of multiple SPP waves by an incident plane wave. Specular reflection and transmission are denoted by order n=0, whereas nonspecular modes are denoted by orders n0. The SPP waves are guided by the periodically corrugated metal/rugate-filter interface.

JNP_6_1_061701_f001.png

The region d1<z<d2 contains a compound surface-relief grating of period L along the x axis. The relative permittivity ϵg(x,z)=ϵg(x±L,z) in this region is taken to be as

Eq. (2)

ϵg(x,z)=ϵm[ϵmϵd(z)]U[d2zg(x)],x(0,L),
for z(d1,d2), with

Eq. (3)

g(x)={12(d2d1)[1+sin(2πxΛ1)],x(0,L1),12(d2d1)[1+sin(2πx-L1Λ2)],x(L1,L),
and

Eq. (4)

U(ζ)={1,ζ0,0,ζ<0.
The chosen compound surface-relief grating is made of two simple surface-relief gratings, one of period Λ1 and the other of period Λ2. Both of the simple gratings have a sinusoidal profile. For simplicity, both simple gratings have the same crest-to-trough distance d2d1, and the depth Lg of the compound surface-relief grating therefore is also the same: Lg=d2d1. Each period of the compound grating comprises N1 consecutive periods of the simple grating with period Λ1 and N2 consecutive periods of the simple grating with period Λ2; hence, L=L1+L2, L1=N1Λ1 (N1{1,2,3,}), and L2=N2Λ2 (N2{1,2,3,}).

We chose this particular grating-shape function g(x) because sinusoidal surface-relief gratings are the easiest type of grating for theoretical studies of the excitation of SPP waves.7 Simple gratings with nonsinusoidal profiles can be investigated with the theoretical formulation adopted here.10 Furthermore, more than two simple gratings can be incorporated in a compound grating.

In the vacuous half-space z0, let a plane wave propagating in the xz plane at an angle θ to the z axis, be incident on the compound surface-relief grating coated with the rugate filter. Hence, the incident, reflected, and transmitted electric field phasors may be written in terms of Floquet harmonics as follows:

Eq. (5)

Einc(r)=n=Ntn=Nt(snas(n)+pn+ap(n))exp[i(kx(n)x+kz(n)z)],z0,

Eq. (6)

Eref(r)=n=Ntn=Nt(snrs(n)+pnrp(n))exp[i(kx(n)xkz(n)z)],z0,

Eq. (7)

Etr(r)=n=Ntn=Nt(snts(n)+pn+tp(n))exp{i[kx(n)x+kz(n)(zd3)]},zd3.
Here,

Eq. (8)

kx(n)=k0sinθ+nκx
is the wavenumber of the Floquet harmonic of order n, with

Eq. (9)

κx=2πL=2πN1L1+N2L2,
and

Eq. (10)

kz(n)={+k02(kx(n))2,k02>(kx(n))2+i(kx(n))2k02,k02<(kx(n))2.
The s and p-polarization states are represented, respectively, by

Eq. (11)

sn=u^y
and

Eq. (12)

pn±=kz(n)k0u^x+kx(n)k0u^z.
The coefficients as(n)=ap(n)=0n[Nt,Nt], except that as(0)0 when the incident plane wave is s-polarized and ap(0)0 when the incident plane wave is p-polarized. The coefficients rs,p(n) and ts,p(n) have to be determined by solving a boundary-value problem.

That boundary-value problem was formulated using the rigorous coupled-wave approach (RCWA)20,21 and was implemented using a stable algorithm.2224 The accuracy of the solution using RCWA depends on the value of Nt, the number of terms of Floquet harmonics in the representation of field phasors being 2Nt+1. The formulation of the boundary-value problem and the numerical algorithm have been explained elsewhere in detail10 and we skip both to go directly to present representative numerical results in Sec. 3.

3.

Numerical Results and Discussion

For illustrative numerical results, we chose a rugate filter with Ω=λ0, na=1.45 and nb=2.32 (Ref. 10). We set the free-space wavelength λ0=633nm and chose the metal to be bulk aluminum (ϵmet=56+21i).

The plane wave absorptance of the entire structure depicted in Fig. 1 is given by

Eq. (13)

A=1n=NtNt|rs(n)|2+|rp(n)|2+|ts(n)|2+|tp(n)|2|as(0)|2+|ap(0)|2Re(kz(n)kz(0)).
The absorptance for p-polarized incidence is denoted by Ap=A|as(0)=0,ap(0)0, and that for s-polarized incidence by As=A|as(0)0,ap(0)=0. Let us note that rs(n)=0 and ts(n)=0n[Nt,Nt] for p-polarized incident plane waves, and vice versa for the s-polarized case, because the partnering dielectric material is isotropic.

We implemented the RCWA in Mathematica and calculated the absorptances Ap and As as functions of θ. We fixed Nt=40 after ascertaining that the absorptances for Nt=40 converged to within ±1% of the absorptances calculated with Nt=41. Also, the electric field phasor E(x,z), the magnetic field phasor H(x,z), and the time-averaged Poynting vector

Eq. (14)

P(x,z)=12Re[E(x,z)×H*(x,z)],
were calculated as appropriate. For this purpose, the region 0<z<d1 was divided into 3-nm-thick slices parallel to the plane z=0, the grating region (d1<z<d2) was divided into 1-nm-thick slices, but the metallic layer (d2<z<d3) was kept as one slice.

3.1.

Canonical Boundary-Value Problem

The corresponding canonical boundary-value problem, when both the rugate filter and the metal are semi-infinite in thickness and their interface is planar, has been solved elsewhere.9 With the assumption that an SPP wave propagates along the x axis with an exp(iκx) variation and is independent of y, the solution of the canonical problem yields values of the wavenumber κ. We found that for the chosen constitutive parameters of both partnering materials, five p-polarized and four s-polarized SPP waves are allowed. Their relative wavenumbers κ/k0 are provided in Table 1.

Table 1

Relative wavenumbers κ/k0 at λ0=633  nm of possible SPP waves guided by the planar interface of aluminum (ϵmet=−56+21i) and a rugate filter (Ω=λ0, na=1.45 and nb=2.32) from the solution of the canonical boundary-value problem.9

s-pol1.48639+0.00132i1.7324+0.00 4i*1.9836+0.0006i*
2.2128+9.6×105i*
p-pol1.36479+0.00169i1.61782+0.00548i1.87437+0.00998i
2.06995+0.01526i2.21456+0.00246i

*These solutions had been missed when solutions for Fig. 1 of Ref. 9 were numerically searched.

3.2.

Grating-Coupled Configuration

The compound surface-relief grating was chosen with N1=N2=3, Λ1=0.75λ0, Λ2=λ0, and Lg=50nm; furthermore, d3d2=30nm. These specific parameters are consistent with experience acquired from working with simple surface-relief gratings.10

Plots of Ap and As as functions of θ are presented in Figs. 2 and 3, respectively, for d1{4Ω,6Ω}. These plots contain several peaks. The peaks that are independent of the thickness of the partnering dielectric material beyond a certain threshold thickness indicate the excitation of SPP waves,10 whereas the remaining peaks could represent the excitation of waveguide modes25,26 that must depend on the thickness of the partnering dielectric material.

Fig. 2

As versus θ, when N1=N2=3, Λ1=0.75λ0, Λ2=λ0, Lg=50nm, d3d2=30nm, and λ0=633nm. The vertical arrows and the stars identify the peaks that represent the excitation of s-polarized SPP waves: the arrows with single line (double line) represent the peaks that are also present in the absorptance curves when the compound grating is replaced with a simple grating with N2=0 (N1=0), but the stars represent the peaks that are not present when either N1=0 or N2=0.

JNP_6_1_061701_f002.png

Fig. 3

Same as Fig. 2, except that Ap is plotted instead of As.

JNP_6_1_061701_f003.png

3.2.1.

Excitation of s-polarized SPP waves

In the plots of As versus θ in Fig. 2, nine peaks, identified by vertical arrows and stars, are present independent of the thickness of the partnering dielectric material (i.e., the rugate filter). The angular positions of these As-peaks are θ=19.0deg, 22.2 deg, 26.7 deg, 31.0 deg, 34.6 deg, 39.8 deg, 49.3 deg, 56.1 deg, and 61.4 deg; all angular positions in this paper were read from absorptance plots with the angle of incidence varying in steps of 0.1 deg. The conclusion that the identified peaks represent the excitation of s-polarized SPP waves is upheld by the spatial profiles of the x-component of the time-averaged Poynting vector in one unit cell (0<x<L,0<z<d3).

In Fig. 4, the variation of Px(x,z)=u^xP˙(x,z) is provided for two values of θ:19.0deg and 22.2 deg. The figure confirms that the power density of the s-polarized SPP waves is maximal in the region near the plane z=d1 and decays away from that region. Moreover, most of the energy of each SPP wave resides in the rugate filter rather than in the metal. As the spatial profiles in the top and bottom panels of Fig. 4 are different for the different values of θ, each of the two values of θ is associated with the excitation of different SPP waves.

Fig. 4

Variation of the x-component Px (in Wm2 ) of the time-averaged Poynting vector P(x,z) when (top) θ=19.0deg and (bottom) θ=22.2deg, the incident plane wave is s polarized (as(0)=1Vm1, ap(0)=0), and d1=4Ω. Other parameters are the same as for Fig. 2.

JNP_6_1_061701_f004.png

In Fig. 5, the spatial variation of Px is provided for two other values of θ:39.8deg and 56.1 deg. This figure again shows that an s-polarized SPP wave is excited at each of these two values of θ; however, the same s-polarized SPP wave is excited at both of these angles of incidence. This observation is not surprising as it has already been shown elsewhere10 that the same SPP wave can be excited in the grating-coupled configuration at different angles of incidence.

Fig. 5

Same as Fig. 4, except that (top) θ=39.8deg and (bottom) θ=56.1deg.

JNP_6_1_061701_f005.png

Let us also note that Px>0 in all the spatial profiles presented in Figs. 4 and 5 implying that these SPP waves propagate parallel to u^x. This was found true for all s-polarized SPP waves excited at the peaks identified in Fig. 2 by vertical arrows and stars.

In the grating-coupled configuration, an SPP wave is excited whenever the real part of the wavenumber κ from the canonical boundary-value problem matches the wavenumber kx(n) of the Floquet harmonic of order n for some n[Nt,Nt]. These wavenumbers of the Floquet harmonics depend on the period L of the surface-relief grating. Also, if the SPP wave is propagating parallel to u^x (u^x) in our formulation of the problem, it is excited as a Floquet harmonic of positive (negative) order.

In the light of the foregoing discussion, let us now explain the roles of the two simple surface-relief gratings constituting the compound surface-relief grating in the excitation of multiple SPP waves. The relative wavenumbers kx(n)/k0 of several Floquet harmonics are provided in Table 2 for:

  • Case 1: N1=1 and N2=0 (i.e., L=Λ1) and

  • Case 2: N1=0 and N2=1 (i.e., L=Λ2),

at the values of θ identified by vertical arrows and stars in Fig. 2. Plots of As versus θ for these two cases are provided in Fig. 6 to compare the As-peaks for the compound surface-relief grating with those for the individual simple surface-relief gratings. For the computation of absorptances when either N1=0 or N2=0, Nt=8 was sufficient because the period L is much smaller than for the compound grating and a converged solution can be obtained with a smaller value of Nt. We also identify

  • Case 3: N1=N2=3,

which is the case of the chosen compound surface-relief grating, for ease of discussion. The relative wavenumbers of several Floquet harmonics are at the θ-positions of the As-peaks identified by stars are provided in Table 3 for case 3 only.

Fig. 6

As versus θ when Lg=50nm, d3d2=30nm, λ0=633nm, for case 1 and case 2. The vertical arrows identify the peaks that represent the excitation of SPP waves.

JNP_6_1_061701_f006.png

Table 2

Relative wavenumbers kx(n)/k0 of Floquet harmonics, at the θ-values of the peaks identified in Fig. 2 by vertical arrows and stars, when the compound surface-relief grating is replaced by just one of its two constituent simple surface-relief gratings. Case 1: N1=1 and N2=1. Case 2: N1=1 and N2=1. Boldface entries signify SPP waves.

Casen=2n=1n=0n=1n=2
θ=19.0deg12.34111.00780.325571.65892.9922
21.67440.674430.325571.32562.3256
θ=22.2deg12.28880.955490.377841.71123.0445
21.62220.622160.377841.37782.3778
θ=26.7deg12.21730.884010.449321.78273.116
21.55070.550680.449321.44932.4493
θ=31.0deg12.15160.81830.515041.84843.1817
21.4850.484960.51501.51502.5150
θ=34.6deg12.09880.765490.567841.90123.2345
21.43220.432160.567841.56782.5678
θ=39.8deg12.02660.693220.640111.97343.3068
21.35990.359890.640111.64012.6401
θ=49.3deg11.90850.57520.758132.09153.4248
21.24190.241870.758131.75812.7581
θ=56.1deg11.83670.503320.830012.16333.4967
21.170.169990.830011.83002.8300
θ=61.4deg11.78870.455350.877982.21133.5446
21.1220.122020.877981.8782.878

Table 3

Relative wavenumbers kx(n)/k0 of Floquet harmonics, at the θ-values of the peaks identified in Fig. 2 by stars when the s-polarized SPP waves are excited by the compound grating. Case 3: N1=N2=3. Boldface entries signify SPP waves.

Casen=3n=4n=5n=6n=7n=8
θ=19.0deg30.8971.08751.27791.46841.65891.8494
θ=26.7deg31.02071.21121.40171.59221.78271.9731
θ=34.6deg31.13931.32971.52021.71071.90122.0917
θ=56.1deg31.40141.59191.78241.97292.16332.3538

The solution of the canonical boundary-value problem predicts the existence of an s-polarized SPP wave with κ/k0=1.7324+0.0014i. According to Table 2, the simple surface-relief grating of period Λ1 alone is responsible for its excitation in the grating-coupled configuration at θ=22.2deg in Fig. 2, because kx(1)/k0|Case1=1.7112 is quite close to 1.7324. Furthermore, the leftmost As-peak identified in Fig. 6(a) by an arrow is at the same angular position as the As-peak of interest in Fig. 2.

Similarly, the solution of the canonical boundary-value problem predicts the existences of s-polarized SPP waves with κ/k0=1.9836+0.0006i and 2.2128+10-5i (Table 1). As-peaks are present at θ=39.8deg and 61.4 deg in Fig. 2. Both of these SPP waves are excited in the grating-coupled configuration due to the presence of the simple grating with period Λ1 because kx(1)/k0|Case1=1.9734 for θ=39.8deg, and kx(1)/k0|Case1=2.2113 for θ=61.4deg in Table 2. Moreover, these As-peaks are also present at θ=39.9deg and 61.4 deg in Fig. 6(a).

The two As-peaks identified by double arrows, at θ{31.0deg,49.3deg}, in Fig. 2 represent the excitation of s-polarized SPP waves due to the presence of simple grating of period Λ2. A comparison of Tables 1 and 2 shows that kx(1)/k0|Case2=1.5150 for θ=31.0deg, and kx(1)/k0|Case2=1.7581 for θ=49.3deg are close to Re(1.4864+0.0013i) and Re(1.7324+0.0014i), respectively. These As-peaks identified by double arrows in Fig. 2 (at θ=31.0deg and 49.3 deg) are also present in Fig. 6(b) at θ=27.9deg and 45.5 deg, respectively. The difference in the angular positions of each peak in the two figures is entirely due to the presence of the simple surface-relief grating of period Λ1 in the compound surface-relief grating.

The As-peaks identified by stars in Fig. 2 at θ{19.0deg,26.7deg,34.6deg,56.1deg} indicate the SPP waves that cannot be explained by either case 1 or case 2, because corresponding peaks are not present in Fig. 6(a) and 6(b). Instead, the entire compound surface-relief grating is involved as follows: At θ=19.0deg, the SPP wave is excited as a Floquet harmonic of order n=6 because kx(6)/k0|Case3=1.4684 is close to Re(1.4864+0.0013i) in Table 1; at θ=26.7deg, the SPP wave is excited as a Floquet harmonic of order n=8 because kx(8)/k0|Case3=1.9731 is close to Re(1.9836+0.0006i) in Table 1; at θ=34.6deg, the SPP wave is excited as a Floquet harmonic of order n=6 because kx(6)/k0|Case3=1.7107 is close to Re(1.7324+0.0014i) in Table 1; and at θ=56.1deg, the SPP wave is excited as a Floquet harmonic of order n=6 because kx(6)/k0|Case3=1.9729 is close to Re(1.9836+0.0006i) likewise. Thus, these SPP waves arise from a cooperative phenomenon which merges the distinct identities of the two constituent simple surface-relief gratings into a single identity. Let us also note that the excitation of the SPP wave at θ=26.7deg is not as efficient as at the other peaks identified by stars in Fig. 2 because the SPP wave is being excited as a Floquet harmonic of higher order when θ=26.7deg than for the other three values of θ (Refs. 7, 10).

A comparison of the magnitudes of As at the absorptance peaks in Figs. 2 and 6 shows that an As-peak for an individual simple surface-relief grating (cases 1 and 2) is generally higher than the corresponding As-peak for the compound surface-relief grating (case 3).

In summary, nine As-peaks in Fig. 2 that represent the excitation of s-polarized SPP waves can be divided into three groups: (1) three peaks due to the presence of the simple surface-relief grating with period Λ1, (2) two peaks due to the presence of the simple surface-relief grating with period Λ2, and (3) four peaks due to the compound surface-relief grating itself with period L=3Λ1+3Λ2. The relative wavenumbers κ/k0 of s-polarized SPP waves in the canonical problem (Table 1) are correlated in Table 4 with the θ-values of the As-peaks in Fig. 2 that represent the excitation of these SPP waves in case 3. Correlation of some of those peaks with the As-peaks in Fig. 6(a) and 6(b), for cases 1 and 2, respectively, is also presented in the same table. All s-polarized SPP waves predicted by the solution of the canonical boundary-value problem can be excited by the compound surface-relief grating, but neither of the two constituent simple surface-relief gratings is, by itself, able to do the same. The As-peaks that are attributable to the simple surface-relief grating with period Λ2 are shifted to somewhat higher angular positions, and the shift becomes more pronounced as the incidence becomes more oblique, whereas the peaks attributable to the other simple grating with period Λ1 remain almost at the same angular position except for one peak that shifts to lower angular position.

Table 4

The relative wavenumbers κ/k0 of s-polarized SPP waves in the canonical problem (Table 1) and the θ-values of the As-peaks in Figs. 6(a), 6(b), and 2 that represent the excitation of these SPP waves in cases 1, 2, and 3, respectively.

κ/k0θ
Fig. 6(a)Fig. 6(b)Fig. 2
1.48639+0.00132i19.0 deg
27.9 deg31.0 deg
1.7324+0.0014i22.2 deg22.2 deg
34.6 deg
45.5 deg49.3 deg
1.9836+0.0006i26.7 deg
39.9 deg39.8 deg
56.1 deg
2.2128+9.6×105i61.4 deg61.4 deg

3.2.2.

Excitation of p-polarized SPP waves

In the plots of Ap versus θ provided in Fig. 3 for d1{4Ω,6Ω}, the Ap-peaks that are independent of the thickness of the rugate filter are identified by vertical arrows and stars. Such Ap-peaks are present for θ{12.0deg,19.5deg,23.5deg,27.6deg,41.0deg,47,5deg,56.3deg} and represent the excitation of p-polarized SPP waves. Examination of the spatial profiles of Px(x,z) at these θ-values confirmed the localization of the power density to the plane z=d1.

Representative profiles of Px in a unit cell of the compound surface-relief grating are provided for two identified Ap-peaks in Fig. 7. The spatial profiles for θ=23.5deg and θ=27.5deg show that the power density of the p-polarized SPP wave is maximal in the region near the plane z=d1 and decays away from that region. Moreover, Px is predominantly positive implying that the p-polarized SPP waves are propagating along u^x and are excited as Floquet harmonics of positive order. This was found to be true for all other p-polarized SPP waves as well.

Fig. 7

Variation of the x-component Px (in Wm2) of the time-averaged Poynting vector P(x,z) when (a) θ=23.5deg and (b) θ=27.6deg, the incident plane wave is p polarized (as(0)=0, ap(0)=1Vm1), and d1=4Ω. Other parameters are the same as for Fig. 3.

JNP_6_1_061701_f007.png

The relative wavenumbers kx(n)/k0 of several Floquet harmonics are provided in Table 5 for case 1 and case 2 at the θ-values of the peaks identified in Fig. 3. Plots of Ap versus θ for these two cases are provided in Fig. 8 to compare the Ap-peaks for the compound surface-relief grating with those for the individual simple surface-relief gratings.

Fig. 8

Same as Fig. 6, except that Ap is plotted against θ.

JNP_6_1_061701_f008.png

Table 5

Same as Table 2, except that the θ-values are of the Ap-peaks identified by vertical arrows and stars in Fig. 3.

Casen=2n=1n=0n=1n=2
θ=12deg12.45881.12540.207911.54122.8746
21.79210.792090.207911.20792.2079
θ=19.5deg12.33290.999530.333811.66713.0005
21.66620.666190.333811.33382.3338
θ=23.5deg12.26790.934580.398751.73213.0654
21.60130.601250.398751.39872.3987
θ=27.6deg12.20340.870040.46331.79663.13
21.53670.53670.46331.46332.4633
θ=41.0deg12.01060.677270.656061.98943.3227
21.34390.343940.656061.65612.6561
θ=47.5deg11.92940.596060.737282.07063.4039
21.26270.262720.737281.73732.7373
θ=56.3deg11.83470.501380.831952.16533.4986
21.1680.168050.831951.8322.832

The Ap-peaks identified by single arrows in Fig. 3 are also present at approximately the same angular positions in Fig. 8(a). Therefore, the p-polarized SPP waves excited by the compound surface-relief grating for θ{red12.0deg,27.6deg,47.5deg,56.3deg} can be attributed solely to the simple surface-relief grating with period Λ1. This is also evident from a comparison of Tables 1 and 5 because:

  • kx(1)/k0|Case1=1.5412 is close to Re(1.61782+0.00548i) for θ=12.0deg,

  • kx(1)/k0|Case1=1.7966 is close to Re(1.87437+0.00998i) for θ=27.6deg,

  • kx(1)/k0|Case1=2.0706 is close to Re(2.06995+0.01526i) for θ=47.5deg, and

  • kx(1)/k0|Case1=2.1653 is close to Re(2.21456+0.00246i) for θ=56.3deg.

A comparison of Figs. 3 and 8(b) shows that the Ap-peaks due to the independent simple surface-relief grating with period Λ1 (case 1) are slightly shifted to lower angular positions when that simple grating becomes part of the compound surface-relief grating (case 3), the shift being larger at the higher angular locations.

Similarly, the Ap-peaks at θ=23.5deg and 41.0 deg, identified by double arrows in Fig. 3, are also present at approximately the same angular positions in Fig. 8(b) signifying that the p-polarized SPP waves at these peaks are excited in effect by the presence of the simple surface-relief grating of period Λ2. Moreover,

  • kx(1)/k0|Case2=1.3987 is close to Re(1.36479+0.00169i) for θ=23.5deg, and

  • kx(1)/k0|Case2=1.6561 is close to Re(1.61782+0.00548i) for θ=41.0deg.

However, the Ap-peaks at θ=20.5deg and θ=37.2deg in Fig. 8(b) for the simple surface-relief grating with period Λ2 shift to θ=23.5deg and θ=41.0deg, respectively, in Fig. 3 when this simple surface-relief grating becomes part of the compound surface-relief grating.

The Ap-peak identified by a star in Fig. 3 also represents the excitation of a p-polarized SPP waves; however, this peak is not present in the plots of Ap versus θ in Fig. 8 for either case 1 or case 2. Also, the Ap-peak at this value of θ is very low, suggesting that the SPP wave is being excited as a Floquet harmonic of a high order at this peak by the compound surface-relief grating with period L=3Λ1+3Λ2. Indeed, the Ap-peak at θ=19.5deg represents the excitation of a p-polarized SPP wave as a Floquet harmonic of order n=7 by the compound grating because kx(7)/k0|Case3=1.6671 is close to Re(1.61782+0.00548i).

Let us note that an Ap-peak present in Fig. 8(b) at θ=60.6deg for the simple surface-relief grating with period Λ2 is not present in Fig. 3 for the compound surface-relief grating. Furthermore, as was shown for As-peaks in Sec. 3.2.1, an Ap-peak for an individual simple surface-relief grating (cases 1 and 2) is generally higher than the corresponding Ap-peak for the compound surface-relief grating (case 3). Thus, although the use of compound surface-relief gratings appears to be efficacious to excite SPP waves, careful selection of the constituent surface-relief gratings is needed for maximal effect.

To recapitulate, out of seven Ap-peaks in Fig. 3 that represent the excitation of p-polarized SPP waves, four peaks are due to the simple surface-relief grating with period Λ1, two to the simple surface-relief grating period Λ2 alone, and one to the compound surface-relief grating itself. The angular positions of the Ap-peaks representing the excitation of p-polarized SPP waves in Fig. 3 are presented in Table 6 are cross-referenced to values of κ/k0 obtained from the solution of the canonical boundary-value problem. Correlation of some of those peaks with the Ap-peaks in Fig. 8(a) and 8(b), for cases 1 and 2, respectively, can also be deduced from the same table. A comparison of Tables 1 and 6 upholds the conclusions drawn in Sec. 3.2.1 for s-polarized SPP waves: All p-polarized SPP waves predicted by the solution of the canonical boundary-value problem can be excited by the compound surface-relief grating even though either of the individual surface-relief grating alone cannot excite all possible SPP waves. Furthermore, the Ap-peaks that are attributable to the simple surface-relief grating with period Λ1 (Λ2) shift to somewhat lower (higher) angular positions in Fig. 3 in relation to Fig. 8.

Table 6

The relative wavenumerbs κ/k0 of p-polarized SPP waves in the canonical problem (Table 1) and the θ-values of the Ap-peaks in Figs. 8(a), 8(b), and 3 that represent the excitation of these SPP waves in cases 1, 2, and 3, respectively.

κ/k0θ
Fig. 8(a)Fig. 8(b)Fig. 3
1.36479+0.00169i20.5 deg23.5 deg
1.61782+0.00548i12.8 deg12.0 
19.5 deg
37.2 deg41.0 deg
1.87437+0.00998i32.3 deg27.6 deg
60.6 deg*
2.06995+0.01526i51.6 deg47.5 deg
2.21456+0.00246i64.8 deg56.3 deg

*This Ap-peak in case 2 does not have a counterpart peak in case 3.

4.

Concluding Remarks

The excitation of multiple SPP waves of the same frequency and the direction of propagation but of different linear polarization states, phase speeds, degrees of localization, and spatial profiles was studied theoretically in the grating-coupled configuration, when the surface-relief grating is compounded from two simpler surface-relief gratings and the partnering dielectric material is a rugate filter. The absorptances for s and p-polarized incident plane waves were calculated as functions of the angle of incidence for two sufficiently high values of the thickness of the partnering dielectric material to rule out waveguide modes and identify SPP waves.

Our numerical results allowed us to draw the following conclusions:

  • (1) With one fixed compound surface-relief grating, all the SPP waves predicted by the canonical boundary-value problem were excited. The extension of this conclusion to the compounding of two or more simple surface-relief gratings in different configurations requires further investigation.

  • (2) The excitation of each SPP wave could be attributed to either a simple surface-relief grating present in the compound surface-relief grating or to the compound surface-relief grating itself or to both.

  • (3) The same SPP wave could be excited by either one or both of the constituting simple gratings of the compound surface-relief grating and/or the compound grating itself.

  • (4) The angular positions of the absorptance peaks (signifying the excitation of SPP waves) for a simple surface-relief grating shifted somewhat when that grating was made a constituent of a compound surface-relief grating.

  • (5) Not every absorptance peak (indicative of the excitation of an SPP wave) for a simple surface-relief grating by itself had a counterpart peak for the compound surface-relief grating, but this conclusion does not nullify conclusion (1).

However, the excitation of any specific SPP wave turned out to be less efficient with the compound surface-relief grating than with a simple surface-relief grating. This observation is in agreement with the work of Dolev et al.8 for whom the partnering dielectric material was homogeneous.

The foregoing conclusions are expected to be useful in enhancing the quantum efficiency of thin-film solar cells because of the greater possibility of exciting SPP waves of both linear polarization states. For this purpose, the rugate filter will have to be replaced by a multilayered semiconductor material and an easily manufacturable compound surface-relief grating will have to be found to excite all possible SPP waves efficiently over a broad range of visible and infrared frequencies. The introduction of periodic nonhomogeneity in the partnering semiconductor would result in the possibility of guiding multiple SPP waves of both linear polarization states.14 The possibility of exciting a given SPP wave by either a constituent simple surface-relief grating or the compound surface-relief grating at more than one angles of incidence (combined with the possibility of multiple SPP waves of the same frequency) will increase the energy coupled to SPP waves by the incident light not only in diffuse insolation but also in direct insolation when the solar cell is either unable or not designed to track the sun. The theoretical formulation of RCWA used in this paper for continuously nonhomogeneous partnering dielectric material has been recently extended for periodically multilayered dielectric materials.27

Acknowledgments

This work was partly supported by Grant No. DMR-1125591 from US National Science Foundation. A.L. is also grateful to Charles Godfrey Binder Endowment at the Pennsylvania State University for partial support of this work.

References

1. 

S. A. Maier, Plasmonics: Fundamentals and Applications, Springer, New York (2007). Google Scholar

2. 

Surface Plasmon Resonance Based Sensors, Springer, Heidelberg, Germany (2006). Google Scholar

3. 

P. ShengA. N. BlochR. S. Stepleman, “Wavelength-selective absorption enhancement in thin-film solar cells,” Appl. Phys. Lett., 43 (6), 579 –581 (1983). http://dx.doi.org/10.1063/1.94432 APPLAB 0003-6951 Google Scholar

4. 

C. HeineR. H. Morf, “Submicrometer gratings for solar energy applications,” Appl. Opt., 34 (14), 2476 –2482 (1995). http://dx.doi.org/10.1364/AO.34.002476 APOPAI 0003-6935 Google Scholar

5. 

S. XiaoE. StassenN. A. Mortensen, “Ultrathin silicon solar cells with enhanced photocurrents assisted by plasmonic nanostructures,” J. Nanophoton., 6 (1), 061503 (2012). http://dx.doi.org/10.1117/1.JNP.6.061503 JNOACQ 1934-2608 Google Scholar

6. 

M. A. GreenS. Pillai, “Harnessing plasmonics for solar cells,” Nature Photon., 6 (3), 130 –132 (2012). http://dx.doi.org/10.1038/nphoton.2012.30 1749-4885 Google Scholar

7. 

H. Raether, Surface Plasmons on Smooth and Rough Surfaces and on Gratings, Springer, Heidelberg, Germany (1988). Google Scholar

8. 

I. Dolevet al., “Multiple coupling of surface plasmons in quasiperiodic gratings,” Opt. Lett., 36 (9), 1584 –1586 (2011). http://dx.doi.org/10.1364/OL.36.001584 OPLEDP 0146-9592 Google Scholar

9. 

M. FaryadA. Lakhtakia, “On surface plasmon-polariton waves guided by the interface of a metal and a rugate filter with sinusoidal refractive-index profile,” J. Opt. Soc. Am. B, 27 (11), 2218 –2223 (2010). http://dx.doi.org/10.1364/JOSAB.27.002218 JOBPDE 0740-3224 Google Scholar

10. 

M. FaryadA. Lakhtakia, “Grating-coupled excitation of multiple surface-plasmon-polariton waves,” Phys. Rev. A, 84 033852 (2011). http://dx.doi.org/10.1103/PhysRevA.84.033852 PLRAAN 1050-2947 Google Scholar

11. 

J. A. Polo Jr.A. Lakhtakia, “Surface electromagnetic waves: a review,” Laser Photon. Rev., 5 (2), 234 –246 (2011). http://dx.doi.org/10.1002/lpor.v5.2 1863-8880 Google Scholar

12. 

A. Lakhtakia, “Surface multiplasmonics,” Proc. SPIE, 8104 810403 (2011). http://dx.doi.org/10.1117/12.893129 PSISDG 0277-786X Google Scholar

13. 

T. G. MackayA. Lakhtakia, “Modeling chiral sculptured thin films as platforms for surface-plasmonic-polaritonic optical sensing,” IEEE Sensors J., 12 (2), 273 –280 (2012). http://dx.doi.org/10.1109/JSEN.2010.2067448 ISJEAZ 1530-437X Google Scholar

14. 

M. FaryadA. Lakhtakia, “Enhanced absorption of light due to multiple surface-plasmon-polariton waves,” Proc. SPIE, 8110 81100F (2011). http://dx.doi.org/10.1117/12.893492 PSISDG 0277-786X Google Scholar

15. 

B. G. Bovard, “Rugate filter theory: an overview,” Appl. Opt., 32 (28), 5427 –5442 (1993). http://dx.doi.org/10.1364/AO.32.005427 APOPAI 0003-6935 Google Scholar

16. 

P. W. Baumeister, Optical Coating Technology, SPIE Press, Bellingham, WA (2004). Google Scholar

17. 

R. Overendet al., “Rugate filter fabrication using neutral cluster beam deposition,” Vacuum, 43 (1–2), 51 –54 (1992). http://dx.doi.org/10.1016/0042-207X(92)90184-X VACUAV 0042-207X Google Scholar

18. 

E. Lorenzoet al., “Fabrication and optimization of rugate filters based on porous silicon,” Phys. Status Solidi C, 2 (9), 3227 –3231 (2005). http://dx.doi.org/10.1002/(ISSN)1610-1642 PSSCGL 1610-1634 Google Scholar

19. 

C.-C. LeeC.-J. TangJ.-Y. Wu, “Rugate filter made with composite thin films by ion-beam sputtering,” Appl. Opt., 45 (7), 1333 –1337 (2006). http://dx.doi.org/10.1364/AO.45.001333 APOPAI 0003-6935 Google Scholar

20. 

M. G. MoharamT. K. Gaylord, “Diffraction analysis of dielectric surface-relief gratings,” J. Opt. Soc. Am., 72 (10), 1385 –1392 (1982). http://dx.doi.org/10.1364/JOSA.72.001385 JOSAAH 0030-3941 Google Scholar

21. 

L. Li, “Multilayer modal method for diffraction gratings of arbitrary profile, depth, and permittivity,” J. Opt. Soc. Am. A, 10 (12), 2581 –2591 (1993). http://dx.doi.org/10.1364/JOSAA.10.002581 JOAOD6 0740-3232 Google Scholar

22. 

M. G. MoharamE. B. GrannD. A. Pommet, “Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings,” J. Opt. Soc. Am. A, 12 (5), 1068 –1076 (1995). http://dx.doi.org/10.1364/JOSAA.12.001068 JOAOD6 0740-3232 Google Scholar

23. 

F. WangM. W. HornA. Lakhtakia, “Rigorous electromagnetic modeling of near-field phase-shifting contact lithography,” Microelectron. Eng., 71 (1), 34 –53 (2004). http://dx.doi.org/10.1016/j.mee.2003.09.003 MIENEF 0167-9317 Google Scholar

24. 

N. ChateauJ.-P. Hugonin, “Algorithm for the rigorous coupled-wave analysis of grating diffraction,” J. Opt. Soc. Am. A, 11 (4), 1321 –1331 (1994). http://dx.doi.org/10.1364/JOSAA.11.001321 JOAOD6 0740-3232 Google Scholar

25. 

N. S. KapanyJ. J. Burke, Optical Waveguides, Academic Press, New York, NY (1972). Google Scholar

26. 

D. Marcuse, Theory of Dielectric Optical Waveguides, Academic Press, San Diego, CA (1991). Google Scholar

27. 

M. Faryadet al., “Excitation of multiple surface-plasmon-polariton waves guided by the periodically corrugated interface of a metal and a periodic multilayered isotropic dielectric material,” J. Opt. Soc. Am. B, 29 (4), 704 –713 (2012). http://dx.doi.org/10.1364/JOSAB.29.000704 JOBPDE 0740-3224 Google Scholar

Biography

JNP_6_1_061701_d001.png

Muhammad Faryad received his BSc degree in mathematics and physics from University of Punjab, Pakistan, in 2002; MSc and MPhil degrees in electronics from Quaid-i-Azam University, Pakistan in 2006 and 2008, respectively; and PhD degree in engineering science and mechanics from the Pennsylvania State University, in 2012. He is a postdoctoral scholar at the Pennsylvania State University. His research experience includes analysis of high-frequency fields reflected from cylindrical reflectors in an isotropic chiral medium and the fractional curl operator in electromagnetics. Currently, his research interests are surface multiplasmonics, sculptured thin films, and thin-film solar cells.

JNP_6_1_061701_d002.png

Akhlesh Lakhtakia received degrees from the Banaras Hindu University (BTech and DSc) and the University of Utah (MS and PhD), in electronics engineering and electrical engineering, respectively. He is the Charles Godfrey Binder (Endowed) professor of Engineering Science and Mechanics at the Pennsylvania State University, and presently serves as the Editor-in-Chief of the Journal of Nanophotonics. His current research interests include nanotechnology, bioreplication, surface multiplasmonics, complex materials, metamaterials, and sculptured thin films. He is a Fellow of SPIE, Optical Society of America, Institute of Physics (UK), and American Association for the Advancement of Science.

© 2012 Society of Photo-Optical Instrumentation Engineers (SPIE) 0091-3286/2012/$25.00 © 2012 SPIE
Muhammad Faryad and Akhlesh Lakhtakia "Excitation of multiple surface-plasmon-polariton waves using a compound surface-relief grating," Journal of Nanophotonics 6(1), 061701 (24 September 2012). https://doi.org/10.1117/1.JNP.6.061701
Published: 24 September 2012
Lens.org Logo
CITATIONS
Cited by 19 scholarly publications.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Optical coatings

Dielectrics

Stars

Interfaces

Metals

Optical filters

Wave propagation


CHORUS Article. This article was made freely available starting 24 September 2013

RELATED CONTENT

Grating resonant filters: an overview
Proceedings of SPIE (September 22 1995)
Surface multiplasmonics
Proceedings of SPIE (September 26 2011)
On The Numerical Study Of Deep Lamellar Gratings In The...
Proceedings of SPIE (December 12 1984)
Airy beam manipulation based on the metallic slit array
Proceedings of SPIE (September 07 2011)

Back to Top