Resonances in complementary metamaterials and ... - OSA Publishing

Resonances in complementary metamaterials and nanoapertures Carsten Rockstuhl,1,∗ Thomas Zentgraf,2,3 Todd P. Meyrath,3 Harald Giessen,3 and Falk Lederer1 1

Institute of Condensed Matter Theory and Solid State Optics, Friedrich Schiller University Jena, Max-Wien-Platz 1, 07743 Jena, Germany

2

Max-Planck-Institut für Festkörperforschung, Heisenbergstr. 1, 70569 Stuttgart, Germany 3

4. Physikalisches Institut, Universität Stuttgart, Pfaffenwaldring 57, 70550 Stuttgart, Germany ∗ Corresponding

author: [email protected]

Abstract: We theoretically analyze the properties of metamaterials which have been designed by taking advantage of Babinet’s principle. It is shown that the complementary structure exhibits both a complementary spectral response and field distribution of the respective eigenmodes. For complementary split-ring resonators, we show that the spectral resonance features can be explained from two different perspectives. On one hand they can be explained as plasmon polariton resonances in dielectric nanostructures surrounded by metal, on the other hand they can be understood as guided mode resonances with vanishing propagation constant. The physical origin of these modes and differences to the conventional split-ring geometry are discussed. © 2008 Optical Society of America OCIS codes: (160.4670) Optical materials; (240.6680) Surface plasmons; (260.3910) Metal optics; (260.5740) Resonance

References and links 1. V. M. Shalaev, “Optical negative-index metamaterials,” Nature Photonics 1, 41-48 (2007). 2. C. M. Soukoulis, S. Linden, and M. Wegener, “Negative refractive index at optical wavelengths,” Science 315, 47-49 (2007). 3. J. D. Joannopoulos Photonic crystals: Molding the flow of light (Princeton University Press, Princeton, New Jersey, 1995). 4. D. R. Smith, S. Schultz, P. Markoˇs, and C. M. Soukoulis, “Determination of effective permittivity and permeability of metamaterials from reflection and transmission coefficients,” Phys. Rev. B 65, 195104 (2002). 5. D. Seetharamdoo, R. Sauleau, K. Mahdjoubi, and Anne-Claude Tarot, “Effective parameters of resonant negative refractive index metamaterials: Interpretation and validity,” J. Appl. Phys. 98, 063505 (2005). 6. C. Rockstuhl, T. Paul, F. Lederer, T. Pertsch, T. Zentgraf, T. P. Meyrath, and H. Giessen, ”Investigating the transition from thin film to bulk properties of metamaterials” Phys. Rev. B 77, 035126 (2008). 7. T. C. Choy, Effective Medium Theory, Principles and Applications, (Oxford University Press, 1999). 8. R. D. Grober, R. J. Schoelkopf, and D. E. Prober, “Optical antenna: Towards a unity efficiency near-field optical probe,” Appl. Phys. Lett. 70, 1354-1356 (1997). 9. L. Lewin, “The electrical constants of a material loaded with spherical particles,” Proc. Inst. Elec. Eng., Part 3, 94, 65-68 (1947). 10. C. Rockstuhl, F. Lederer, C. Etrich, T. Pertsch, and T. Scharf, “ Design of an Artificial Three-Dimensional Composite Metamaterial with Magnetic Resonances in the Visible Range of the Electromagnetic Spectrum,” Phys. Rev. Lett. 99, 017401 (2007).

#90887 - $15.00 USD

(C) 2008 OSA

Received 14 Dec 2007; revised 18 Jan 2008; accepted 22 Jan 2008; published 30 Jan 2008

4 February 2008 / Vol. 16, No. 3 / OPTICS EXPRESS 2080

11. J. P. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, “Magnetism from conductors, and enhanced nonlinear phenomena,” IEEE Trans. Microwave Theory Tech. 47, 2075-2084 (1999). 12. F. Falcone, T. Lopetegi, M. A. G. Laso, J. D. Baena, J. Bonache, M. Beruete, R. Marqués, F. Mart´ın, and M. Sorolla, “Babinet principle applied to the design of metasurfaces and metamaterials,” Phys. Rev. Lett. 93, 197401 (2004). 13. H.-T. Chen, J. F. O’Hara, A. J. Taylor, R. D. Averitt, C. Highstrete, M. Lee, and W. J. Padilla, “Complementary planar terahertz metamaterials,” Opt. Express 15, 1084-1095 (2007). 14. M. Born and E. Wolf, Principles of optics (Cambridge University Press, Cambridge, 1999). 15. R. Ullrich, “Far-infrared properties of metallic mesh and its complementary structure,” Infrared Phys. 7, 37-55 (1967). 16. L. C. Botton, R. C. McPhedran, and G. W. Milton, “Perfectly conducting lamellar gratings: Babinet’s principle and circuit models,” J. Mod. Opt. 42, 2453-2473 (1995). 17. T. Zentgraf, C. Rockstuhl, T. P. Meyrath, A. Seidel, S. Kaiser, F. Lederer, and H. Giessen, “Babinet’s principle for optical frequency metamaterials and nanoantennas,” Phys. Rev. B 76, 033407 (2007). 18. C. Rockstuhl and F. Lederer, “Negative-index metamaterials from nanoapertures,” Phys. Rev. B 76, 125426 (2007). 19. J. A. Porto, F. J. Garca-Vidal, and J. B. Pendry, “Transmission Resonances on Metallic Gratings with Very Narrow Slits,” Phys. Rev. Lett. 83, 2845-2848 (1999). 20. X. Shi, L. Hesselink, and R. L. Thornton, “Ultrahigh light transmission through a C-shaped nanoaperture,” Opt. Lett. 28, 1320-1322 (2003). 21. C. Rockstuhl, F. Lederer, T. Zentgraf, and H. Giessen, “Enhanced transmission of periodic, quasi-periodic, and random nanoaperture arrays,” Appl. Phys. Lett. 91, 151109 (2007). 22. F. J. Garc´ıa de Abajo, R. Gómez-Medina, and J. J. Sáenz, “Full transmission through perfect-conductor subwavelength hole arrays,” Phys. Rev. E 72, 016608 (2005). 23. F. J. Garc´ıa de Abajo, “Light scattering by particle and hole arrays,” Rev. Mod. Phys. 79, 1267-1290 (2007). 24. P. B. Johnson and R. W. Christy, “Optical constants of the noble metals,” Phys. Rev. B 6, 4370-4379 (1972). 25. D. H. Dawes, R. C. McPhedran, and L. B. Whitbourn, “Thin capacitive meshes on a dielectric boundary: theory and experiment,” Appl. Opt. 26, 3498-3510 (1989). 26. L. Li, “New formulation of the Fourier modal method for crossed surface-relief gratings,” J. Opt. Soc. Am. A 14, 2758-2767 (1997). 27. C. Rockstuhl, T. Zentgraf, C. Etrich, J. Kuhl, F. Lederer, and H. Giessen, “On the reinterpretation of resonances in split-ring-resonators at normal incidence,” Opt. Express 14, 8827-8836 (2006). 28. C. Rockstuhl, T. Zentgraf, H. Guo, N. Liu, C. Etrich, I. Loa, K. Syassen, J. Kuhl, F. Lederer, and H. Giessen, “Resonances of split-ring resonator metamaterials in the near infrared,” Appl. Phys. B. 84, 219-227 (2006). 29. A. Farjadpour, D. Roundy, A. Rodriguez, M. Ibanescu, P. Bermel, J. D. Joannopoulos, S. G. Johnson, and G. Burr, “Improving accuracy by subpixel smoothing in the finite-difference time domain,” Opt. Lett. 31, 29722974 (2006). 30. C. Rockstuhl, T. Zentgraf, E. Pshenay-Severin, J. Petschulat, A. Chipouline, J. Kuhl, H. Giessen, T. Pertsch, and F. Lederer, “The origin of magnetic polarizability in metamaterials at optical frequencies - an electrodynamic approach,” Opt. Express 15, 8871-8883 (2007). 31. Z. Zhu and T. Brown, “Full-vectorial finite-difference analysis of microstructured optical fibers,” Opt. Express 10, 853-864 (2002). 32. V. V. Varadan and A. R. Tellakula, “Effective properties of split-ring resonator metamaterials using measured scattering parameters: Effect of gap orientation,” J. Appl. Phys. 100, 034910 (2006). 33. M. Kafesaki, T. Koschny, R. S. Penciu, T. F. Gundogdu, E. N. Economou, and C. M. Soukoulis, “Left-handed metamaterials: detailed numerical studies of the transmission properties,” J. Opt. A: Pure Appl. Opt. 7, S12-S22 (2005).

1. Introduction Metamaterials (MMs) are media that permit control over the propagation characteristics of light [1, 2]. They derive their properties primarily from their geometry rather than from intrinsic material properties of the constituents. MMs usually consist of periodically arranged unit cells so that the light propagation is entirely governed by the dispersion relation of the respective Bloch modes [3]. Frequently, effective material parameters for MMs can be derived if a single Bloch mode dominates the propagation characteristics of the wave field [4]. If metals are incorporated into the unit cell, this dominating Bloch mode has a significantly smaller imaginary part of its complex wave vector when compared to other Bloch modes [5, 6]. Additionally, strong dispersion of the effective material parameters can be encountered if resonances are employed. #90887 - $15.00 USD

(C) 2008 OSA



These resonances must be independent from the period of the structure and are determined by the response of the unit cell. Such an approach resembles the concept of a polarizable atom where an effective polarization can be determined upon spatial averaging [7]. In designing a MM, one should therefore aim for a nanostructure that shows a resonantly enhanced scattering response similar to an optical nanoantenna [8]. If the scattering response is dominated by an electric/magnetic dipole, the effective permittivity/permeability is altered. Presently, two resonance mechanisms are employed to control the MMs response. The first is based on Mie resonances in dielectric spheres [9, 10] and the second on plasmon polariton resonances in appropriately shaped metallic nanoparticles [11]. By applying Babinet’s principle to the latter class of MMs, a resonant scattering response has been encountered in the microwave domain [12, 13]. Basically, Babinet’s principle states that a complementary structure illuminated by a complementary wave field (for a plane wave this implies that the polarization is rotated by 90 o in the normal direction) causes a complementary scattering response with reflectance and transmittance being interchanged [14]. Babinet’s principle in its original formulation requires the scattering structures to be made of an infinitely thin perfect conductor. Babinet’s principle is a fundamental principle in optics with a variety of applications. In the past, it was exploited to understand light propagation in conductive and inductive wire grids [15, 16]. Recently, experiments with the complementary split-ring resonator (c-SRR) have probed the limits of the applicability of Babinet’s principle to design MMs with resonance frequencies in the optical domain [17]. It was shown that the principle remains valid, although the strict requirements, as imposed in its original formulation, are violated to a certain extent. For c-SRR structures at normal incidence the complementary behavior comprises the scattering response, the classification of observable resonances, and the effective material parameters. As an open question, it remains to clarify to which extent the physical origin of the resonances in complementary MMs is likewise complementary to the excitation of plasmon polaritons in MMs. The understanding of the resonance mechanism will certainly lead to more efficient design approaches for this MM type, which is mandatory for the transfer of the concept into practical devices. With the present work, we aim at bridging this gap with a detailed analysis of the origin of resonances in c-SRRs. However, this physical analysis applies to nanoaperture MMs in general, in particular to apertures having a rectangular or circular shape [18]. Our observations are connected to the work related to enhanced light transmission through sub-wavelength apertures evoked by vertical resonances [19, 20, 21] which are essentially Fabry-Pérot resonances in the nanoapertures. Additionally, there are also horizontal resonances causing an enhanced transmission by the excitation of surface-plasmon-polaritons (SPPs) on the periodically patterned metal film surface[22]. A comprehensive review about the scattering of light at particles and hole arrays can be found in Ref. 23. Enhanced transmission evoked by the latter effect depends sensitively on the nanoaperture period and can be disregarded in the present work where the resonance spectrum is independent of the period. We show that Babinet’s principle is an excellent means of applying a unified point of view on various fields of modern optics, e.g., light propagation in MMs and through sub-wavelength apertures as well as localized plasmons in metallic nanoparticles and optical nanoantennas. Although not all implications can be revealed in their entirety in the present work, we shall put forward some essential ideas. 2. Properties of the resonant modes First, we compare the scattering response of a split-ring resonator (SRR) with that of its complementary structure (c-SRRs) and summarize the essential features. Most of these properties #90887 - $15.00 USD

(C) 2008 OSA



SRR

Normalized intensity Normalied intensity

0.6

T R

0.4 0.2

(a)

0

2000

4000

6000

8000

10000

Wavenumber in cm-1 1

Normalizedintensity intensity Normalied

1

(3)

(1)

0.8

(2)

0.6 0.4 0.2 0

2000

4000

6000

(2c)

0.6 0.4

EInc

0.2 0

2000

8000 -1

Wavenumber in cm

10000

4000

6000

8000

10000

BInc

Wavenumber in cm-1 1

0.8

(b)

0.8

(c) Normalizedintensity intensity Normalied

Normalied intensity Normalized intensity

1

c-SRR

(3c)

(1c)

0.8 0.6 0.4

BInc

0.2

(d)

0

2000

4000

6000

8000

10000

EInc

-1

Wavenumber in cm

Fig. 1. Normalized reflected (dashed red lines) and transmitted intensity (solid blue lines) for SRRs (a, b) and c-SRRs (c, d), as shown in the geometry on top. The structures are illuminated by a plane wave at normal incidence with its electric field component polarized parallel (a,c) or perpendicularly (b,d) to the gap. The numbers indicate the order of the resonance which can be deduced from the amplitude nodes of the electric/magnetic field component normal to the SRR/c-SRR plane. The SRR/c-SRR is characterized by an arm length of l = 300 nm, a width of w = 40 nm, and a height for the SRR of h = 15 nm which is the thickness of the metal film for the c-SRR. They are assumed to be arranged in a periodic lattice with the period being Λ = 400 nm.

have already been reported in the literature [12, 13, 17]. For our investigations, the geometry of both structures (SRR and c-SRR) is characterized by an arm length of l = 300 nm, a width of w = 40 nm, and a height of h = 15 nm. The gap is completely open on one side of the ring. They are assumed to be arranged in a periodic lattice with the period being Λ = 400 nm. In the simulations, the SRR is assumed to be made of gold and the c-SRR is assumed to be the respective air nanoaperture in a thin gold film. The dielectric function of gold has been taken from Ref. 24. The effect of an experimentally necessary substrate is omitted in our work as it has solely a weak quantitative impact on the optical response and, strictly speaking, Babinet’s principle to be applied requires a symmetric cladding [25]. Figure 1 shows the normalized reflected and transmitted intensity for normal incidence and the two polarization directions. All spectral calculations have been performed by using the Fourier modal method [26]. In the simulation 31 × 31 plane waves were retained, sufficiently enough to ensure convergence of all quantities in the order of a percent. We observe resonances appearing as reflection/transmission peaks for the SRR/c-SRR at approximately the same wavenumbers, if the polarization of the illuminating wave field is interchanged. The resonances are labeled in ascending order of the wavenumber. For SRR resonances, that number equals the number of amplitude nodes of the electric field component normal to the SRR surface [27]. The origin of these resonances is the excitation of plasmonic eigenmodes in the gold

#90887 - $15.00 USD

(C) 2008 OSA



nanostructures [28]. The larger line width of the c-SRR resonances in Fig. 1 can be potentially attributed to the larger part of the electric mode volume in gold, which increases the absorption losses. Furthermore, the spectral shift of the higher-order resonances is attributed to a partial violation of the validity conditions of Babinet’s principle. Particularly, the conductivity of gold degrades at optical frequencies and the ratio of structure height to wavelength becomes larger - contradicting the assumption of nanoapertures in an infinitely thin and perfectly conducting film. It has to be mentioned that the skin depth of gold in the spectral region of interest takes approximately 24 nm, which is already comparable to the thickness of the metal film. -0.2

-0.2

EInc

1

0.1

-0.2

0.8

BInc

0

0.6 0.4

-0.1

0

x in !m

0.1

0.2

0.1

0.2

(2c)

(b)

0.2 -0.2

-0.1

0

x in !m

0.1

0.2

0

(c)

(3c) 0.2 -0.2

2

BInc

-0.1

0.1

0.5

(1c) 0.2

(a)

yy in !m !m

y in !m

1.5 0

EInc

-0.1

1

in !m !m yy in

BInc

-0.1

-0.2

2

1.5

EInc

1 0.5

-0.1

0

x in !m

0.1

0.2

Fig. 2. Amplitude of the magnetic field component normal to the c-SRR surface at the spectral positions of the first three lowest-order eigenmodes. The electric polarization of the incident plane wave is normal to the c-SRR gap for odd (a, c) and parallel for even eigenmodes (b).

The close relation between resonances of the SRR and the c-SRR can not only be deduced from the spectral response, but also from the field distribution of the eigenmodes. Babinet’s principle states that E − cBc = EInc , therefore at resonance, where the strength of the incident field is negligible (EInc E, Bc ), the total electric field of the SRR eigenmode must resemble the total magnetic field of the c-SRR eigenmode [17]. Here, E is the total electric field behind the structure, Bc is the total magnetic field behind the complementary structure, and E Inc is the illuminating electric field component. This analogy can be easily perceived from Fig. 2. There, the amplitude of the magnetic field component, normal to and shortly (20 nm) behind the surface of the c-SRR, is shown for the three lowest-order resonance frequencies of Fig. 1(c,d). The fields were calculated by using the finite-difference time-domain (FDTD) method [29]. As already outlined above, a striking similarity with the eigenmodes of a SRR is observed (see Ref. [27]). The c-SRR higher order modes have an increasing number of nodes in the magnetic field component normal to the c-SRR plane. This similarity of field shapes is not restricted to the normal component, but holds for the other components too. 3. Origin of resonances in complementary metamaterials After having outlined the basic properties of c-SRR resonances, we will proceed in clarifying their physical origin. In the optical domain SRR resonances are associated with plasmonic eigenmodes [27, 30]. For c-SRRs the mechanism is very similar but even richer in some details. The origin of these resonances can be explained by using two different pictures, namely on the one hand plasmon polariton resonances in a dielectric (air) nanostructure embedded into metal, and on the other hand Fabry-Pérot resonances of guided waves propagating perpendicular to the structure. Figure 3 shows the transmittance of the c-SRR as a function of the metal film thickness for the two possible polarizations. The spectral positions of the resonances at ν¯ 1c = 4000 cm−1 , ν¯ 2c = 8500 cm−1 , and ν¯ 3c = 10000 cm −1 are almost independent of the film thickness, although a minor shift to lower frequencies is observed for decreasing thickness, particularly for #90887 - $15.00 USD

(C) 2008 OSA



0.6

0.5

0.5

3c

0.4 0.3 0.2

1.5

0.1 2

(a)

0.2

0.4

Height in !m

0.6

0.8

-1

Wavenumber in cm

-1

1c

BInc

4

x 10

Wavenumber in cm

4

x 10

1

EInc

BInc EInc 0.5

2c

0.6

1

0.4

4c 1.5

(b)

2

0.2 0.2

0.4

Height in !m

0.6

Fig. 3. Normalized transmitted intensity through the c-SRR as a function of metal thickness for different polarizations of the incident electric field. The lowest order Fabry-Pérot resonance (kz d ≈ 0) of the four lowest order guided modes are indicated by the dashed lines. For an increasing thickness higher order Fabry-Pérot resonances of these guided modes move asymptotically towards the spectral position of the lowest order resonance .

the higher-order resonances. An additional weak resonance (4 th order eigenmode) is observed at ν¯ 4c = 15000 cm −1 in Fig. 3(b) with a resonance frequency already in the spectral domain where the conductivity of gold is abated. However, this eigenmode behaves similarly as the 2 nd order eigenmode which was discussed in detail. Most notably, we further observe that for all resonances the transmittance degrades only weakly with increasing thickness. In addition to these thickness independent transmission resonances, a further set of (weaker) resonances appears in the spectrum. With increasing thickness their spectral positions tend asymptotically towards the corresponding fundamental transmission resonance. This response resembles that observed for enhanced transmission through nanoapertures in the presence of vertically propagating guided modes [19]. The transmission spectrum can be qualitatively understood in terms of the plane wave response of a symmetric Fabry-Pérot resonator, given by T =

t 2 eikz h , 1 + r2ei2kz h

(1)

which also holds in first approximation for guided modes. Here, T is the amplitude transmission coefficient, t and r are the complex amplitude transmission and reflection coefficients of the guided mode, k z is their propagation constant, and h is the height of the structure. If a transmission resonance is independent of the height of the structure (see Fig. 3), the guided mode has to be close to its cut-off frequency (all modes in a dielectric waveguide embedded into metal have a zero cut-off propagation constant). Thus the√real part of the wave vector has to be close ω 2 −ω 2

c , where ωc is the cut-off frequency. to zero. Hence, k z can be approximated by k z ≈ c Higher-order transmission peaks at larger height occur because the guided mode is strongly dispersive (see Fig. 5). For a given height of the structure, there exist wavenumbers which minimize the denominator in Eqn. 1, thereby enhancing transmission. Therefore, these resonances may be understood as higher order Fabry-Pérot resonances of the guided eigenmodes.

#90887 - $15.00 USD

(C) 2008 OSA



-0.2

0

10

EInc

5

0.1 x

-0.1

-0.2

y in !m

0 0.1 x in !m

4

0.1

2 x

-0.1

0

x in !m

0.1

0.2

-0.2 15

8

EInc

-0.1 10

EInc

0

5

0.1

Ey

1c

0.2 -0.2

BInc

0

(c) 0.2-0.2

0.2

BInc

-0.1

(b)

6

E2c

E1c

0.2 -0.2

y in !m

BInc

8

EInc

-0.1

y in !m

y in !m

-0.1

(a)

-0.2

15

6

BInc

0

4

0.1

2

E2y c

-0.1

0

x in !m

0.1

0.2

(d) 0.2-0.2

-0.1

0

x in !m

0.1

0.2

Fig. 4. Ex and Ey component of the first-(a,b) and second-order (c,d) guided mode. The mode profiles were calculated by solving the eigenvalue equation in an infinite structure normal to the c-SRR surface for the eigenvalues (kz ≈ 0). These mode fields are qualitatively identical to those of plasmon resonances obtained by calculating the transmission of a h = 15 nm film by using the Fourier Modal Method.

For the very small thickness (h = 15 nm) considered here, these higher-order transmission resonances are therefore not excited. For larger thicknesses, they move towards the mode cutoff frequency, as observed in Fig. 3. It is now clear that for very small thicknesses, the system adopts to the very first resonance at (ℜkz h = 0) which entails (ℜkz = 0), hence c-SRR resonances will appear at the cut-off frequencies of the modes. The relationship between c-SRR and nanoaperture cavity resonances can be understood by looking at the field distribution and the dispersion relation of the lowest-order guided modes for a nanoaperture of infinite thickness. For this purpose, the wave equation has been solved with the ansatz E(x, y, z) = E0 (x, y)eikz z using a finite-difference scheme [31]. The amplitudes of the x- and y-component of the electric field for the two lowest-order guided modes (1 c and 2c ) are shown in Fig. 4. They are almost identical (by visual inspection) to the corresponding amplitudes of the electric field of a c-SRR with a finite height at normal light incidence, when the frequency of the illumination corresponds to the first- or second-order plasmon resonance, respectively (see Fig. 1). The x-components for both guided modes are nearly identical [Fig. 4(a) and (c)], whereas the y-components differ in the number of nodes. We conclude that the observed spectral resonance originates in a coupling of the incident field to a guided eigenmode inside the nanoaperture. As the guided modes have a continuous spectrum, it is not a priori clear at which frequency the transmission resonance will appear. However, this question can be addressed in analyzing the entire dispersion relation of the nanoaperture modes. Figure 5 shows the effective refractive index n eff of the first three eigenmodes of the infinitely thick c-SRR as a function of the frequency. The refractive index was derived from the calculated propagation constant k z by neff = kz /k0 with k0 = 2π /λ being the wavenumber at a #90887 - $15.00 USD

(C) 2008 OSA



1 0.5 0

(a)

1st mode 2nd mode 3rd mode

8

3rd mode

eff

Real(n ) Re

1.5

10

1st mode 2nd mode

Imag(n Im eff)

2

6 4 2

0.5

1

1.5 -1

Wavenumber in cm

2 4

x 10

0

(b)

0.5

1

1.5

Wavenumber in cm-1

2 4

x 10

Fig. 5. Effective refractive index of the first three eigenmodes of the c-SRR nanoparture separated in the real (a) and the imaginary (b) part.

particular wavelength. A large imaginary part of n eff implies that the guided mode is strongly attenuated upon propagation, therefore being nearly evanescent. The mode changes its character at a certain frequency (cut-off frequency). There, the imaginary part of the mode tends towards zero. As it remains close to zero beyond this frequency, the guided mode propagates in the nanoaperture with low losses. In contrast to the imaginary part of n eff , the real part remains close to zero below the cut-off frequency. At frequencies above the cut-off frequency, the modes have a propagating character with a strong dispersion in the real part of n eff . The cut-off frequencies for the three guided modes displayed in Fig. 5 coincide precisely with the height-independent transmission resonances of the c-SRR (shown in Fig. 3). At cutoff frequency, the calculated propagation constant shows exactly the properties required for a transmission resonance which is independent of the thickness of the structure. Both, the real and the imaginary parts of the propagation constant are close to zero. The slight increase in the imaginary part of n eff for all modes at frequencies beyond 15000 cm −1 is associated with the interband absorption of gold in this spectral domain. 4. Eigenmodes at parallel light incidence An essential issue for SRR MMs is the comparison of resonant modes excited for normal and in-plane incidence [33]. In this section, we will address this issue for c-SRRs. In the following, we restrict ourselves to the case where the electric field is also in-plane. Then two c-SRR orientations have to be discriminated. The c-SRR gap can be situated either normal or parallel with respect to the illuminating wave. Figure 6(a) shows the reflected and transmitted intensity for in-plane illumination of normal oriented c-SRRs where it does not matter whether the gap is at the front or the back side except a phase difference for asymmetric structures [32]. The single resonance observed at ν¯ 2c = 7000 cm−1 is close to the resonance of the second-order guided mode. To verify that this guided mode is indeed excited, the z-component of the magnetic field 20 nm above the c-SRR plane is shown in Fig. 6(b). Furthermore, the electric field in the central plane of the c-SRR, split into its x- and y-component, is shown in (c) and (d), respectively. Note that the height of the c-SRR is again only 20 nm. Indeed, the field distribution corresponds to that of the secondorder mode. The magnetic field component perpendicular to the c-SRR exhibits two amplitude minima which are not exactly zero because this field is a superposition of the incident and the mode field. Nevertheless, the expected amplitude modulation of the field along the c-SRR is clearly seen. Furthermore, the electric field equals the second-order guided mode as calculated #90887 - $15.00 USD

(C) 2008 OSA



xxin in !!m m

-0.2

T R

0.8

1.2 1

kInc

0 0.1

(b)

0.2

0.8 0.6

EInc -0.2

Bz2c

0.4 0.2 -0.1

0

y in !m

0.1

y in !m

0.2

0.3

-0.2

0.6

xx in in !!m m

Normalied intensity Normalized intensity

1

1.4

-0.1

0.4

0 0.1

(c)

0.2

6

-0.1

0.2

EInc -0.2

4

kInc

2

-0.1

0

0.1

in !!m m yyin

0.2

0.3

-0.2

(a)

2000

4000

6000

8000 -1

Wavenumber in cm

10000

x xinin !m !m

0

(d)

-0.1

8

0.2

Air

Air

kInc

0 0.1

E2yc

Au

EInc -0.2

Air

4

0

0.1

in !!m m yyin

E2xc

2

Air -0.1

6

0.2

0.3

Fig. 6. (a) Normalized reflected (red dotted line) and transmitted intensity (blue bold line) for an in-plane illuminated c-SRR with its gap at the front side. The amplitude of the magnetic field component perpendicular to the c-SRR and 20 nm above its surface at the resonance frequency of ν¯ 2c = 7000 cm−1 is shown in (b). The amplitude of the x- and ycomponent of the electric field amplitude at the same frequency are displayed in (c) and (d), respectively.

for the infinite extended system and for normal incidence (see Fig. 4). In contrast to a SRR, which supports in this configuration odd modes, only even modes are excited in the c-SRR. This response can be explained by the symmetry constraints. The illuminating magnetic field has an even symmetry, hence it can only couple to eigenmodes with the same symmetry. Therefore, the first- and the third-order modes are not excited. This adds another complementary aspect, as for SRRs, only odd modes are excitable in this configuration [30]. The spectral response of a c-SRR for parallel gap orientation is shown in Fig. 7(a). Two pronounced resonances are observed at frequencies corresponding to even modes for excitation normal to the c-SRR plane (see Fig. 1). The calculation of the magnetic field component B z perpendicular to the c-SRR plane for the resonance frequencies of ν¯ 1c = 3500 cm−1 and ν¯ 3c = 8000 cm −1 reveals that the resonances correspond indeed to the first- and third-order guided modes [Fig. 7(b) and (c)]. The second-order guided mode may also be excited because no restrictions on the symmetry are imposed by the structure. It appears as an extremely weak dip in reflectance at ν¯ 2c = 6000 cm −1 . The weakness of this excited resonance can be explained by analyzing the y-component of the electric mode field. This component is provided by the incident field. The electric field is displayed for the first two lowest modes in Fig. 4. The ycomponent of the first-order (and likewise for the third-order) modes shows no node along the bottom of the c-SRR (opposite to the gap). Therefore, an illuminating plane wave - showing a phase variation much larger than the c-SRR size - can excite these modes easily. In contrast, the amplitude has a single node for the second-order resonance, hence reducing significantly a potential coupling.

#90887 - $15.00 USD

(C) 2008 OSA



-0.2

yyin in !m !m

T R

0.8

EInc

1.5 1

0.2

0

x in m

0.1

x in !!m

0.2

2

-0.1

2000

4000

6000

8000

10000

Wavenumber in cm-1

(c)

EInc

0 0.1 0.2

Bz1c

0.5 -0.2 -0.1

-0.2

0.2

(a)

kInc

0.1

(b)

0.4

2

0

0.6

0

2.5

-0.1

in ! m yy in !m

Normalied Normalizedintensity intensity

1

1

Au Air

Air -0.2 -0.1

1.5

kInc

0

Air 0.1

m xx in !!m

Bz3c

0.5

0.2

Fig. 7. Reflected and transmitted normalized intensity for an in-plane illuminated c-SRR with parallel orientation (a). The amplitude of the magnetic field component perpendicular to the c-SRR and 20 nm above its surface at the resonance frequency of ν¯ 1c = 3500 cm−1 are shown in (b) and at the resonance frequency of ν¯ 3c = 8000 cm−1 are shown in (c).

5. Conclusion In conclusion, we have analyzed the resonance mechanisms which govern the spectral response of complementary MMs. The geometry of such MMs is derived by applying Babinet’s principle to geometries of MMs that employ localized plasmon polaritons as resonance mechanisms. For complementary MMs the resonances have been shown to be Fabry-Pérot resonances guided modes in nanocavities. These guided modes are resonances of the infinitely extended nanostructure. It has been shown that for complementary MMs with a near zero thickness, resonances appear only at spectral positions where the eigenvalue disappears (cut-off frequencies). This makes the spectral positions of these resonances independent of the height of the MM. The spatial distribution of the field of the resonances in c-SRRs as compared to the field of resonances in SRRs is shown in the conceptual sketch in Fig. 8. The field distributions are extracted from rigorous simulations using FDTD. It can be seen that the resonances in the SRR are associated with plasmons, with an electric current inside the metallic nanoparticle that oscillates parallel to the structure. The resonances in the c-SRR are associated with the field that oscillates perpendicular to the structure. With increasing resonance number, an additional node appears in both amplitudes and a phase jump of π occurs. At increased height, further resonances appear in the spectra for the c-SRRs, corresponding to higher order Fabry-Pérot transmission resonances of the guided modes in the nanocavity. These resonances are similarly excitable for an in-plane illumination relative to the complementary structure. Although we have analyzed in the present work the example of the complementary split-ring resonator, the origin of the resonance mechanism applies equally to all the other possible shapes of complementary MMs, such as rectangular or annular apertures. By employing Babinet’s principle in the design of MMs, we expect to see in the near future a constant flow of ideas among the different fields of nanophotonics. Acknowledgments Parts of the computations were done on the IBM p690 cluster (JUMP) of the ForschungsZentrum in Jülich, Germany. T.M. and T.Z. thank the Alexander von Humboldt Foundation and #90887 - $15.00 USD

(C) 2008 OSA



(1)

(2)

(3)

Fig. 8. Conceptual sketch of the electric field in the nanoresonators that form the MM for the first three resonances. At a fixed time, the arrows denote the electric field, where their directions resolve the phase relation, and their thickness resolve the excitation strength. Resonances in the SRR are shown in the top row, resonances in the c-SRR are shown in the bottom row.

the Landesstiftung Baden-Württemberg, respectively. Our work was supported by the German Federal Ministry of Education and Research (Grant No. 13N9155).

#90887 - $15.00 USD

(C) 2008 OSA

