Overcoming the losses of a split ring resonator array with gain

Overcoming the losses of a split ring resonator array with gain Anan Fang, 1 Zhixiang Huang, 1, Thomas Koschny, 1,3 and Costas M. Soukoulis 1,3, 1 Department of Physics and Astronomy and Ames Laboratory, Iowa State University, Ames, Iowa 511, USA Key Laboratory of Intelligent Computing and Signal Processing, Anhui University, Hefei 339, China 3 Department of Materials Science and Technology and Institute of Electronic Structure and Laser, FORTH, University of Crete, 7111 Heraklion, Crete, Greece *soukoulis@ameslab.gov Abstract: We present a computational approach, allowing for a selfconsistent treatment of a split ring resonator (SRR) array with a gain layer underneath. We apply three different pumping schemes on the gain layer: (1) homogeneously pumped isotropic gain, () homogeneously pumped isotropic gain with a shadow cast by the SRR and (3) anisotropic gain pumped in a selected direction only. We show numerically the magnetic losses of the SRR can be compensated by the gain. The difference on loss compensations among the three pumping schemes is analyzed by the electric field distribution. Studies also show the dielectric background of gain does not affect the loss compensation much for the gain only pumped in the direction parallel to the SRR plane. 11 Optical Society of America OCIS codes: (16.476) Optical properties; (16.3918) Metamaterials; (6.574) Resonance; (35.438) Nanophotonics. References and links 1. J. B. Pendry, Negative refraction, Contemp. Phys. 45, 191 (4).. S. A. Ramakrishna, Physics of negative refractive index materials, Rep. Prog. Phys. 68, 449 51 (5). 3. C. M. Soukoulis, M. Kafesaki and E. N. Economou, Negative-index materials: New frontiers in optics, Adv. Mater. 18, 1941 195 (6). 4. V. M. Shalaev, Optical negative-index metamaterials, Nature Photon. 1, 41 48 (7). 5. C. M. Soukoulis, S. Linden, and M. Wegener, Negative refractive index at optical wavelengths, Science 315, 47 49 (7). 6. C. M. Soukoulis and M. Wegener, Optical metamaterials more bulky and less lossy, Science 33, 1633 1634 (1). 7. F. Capolino,Theory and Phenomena of Metamaterials (CRC Press, Taylor and Francis Group, 9). 8. J. B. Pendry, Negative refraction makes a perfect lens, Phys. Rev. Lett. 85, 3966 3969 (). 9. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, Metamaterial electromagnetic cloak at microwave frequencies, Science 314, 977 98 (6). 1. J. Zhou, Th. Koschny, and C. M. Soukoulis, An efficient way to reduce losses of left-handed metamaterials, Opt. Express 16, 11147 1115 (8). 11. J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, G. Bartal, and X. Zhang, Three-dimensional optical metamaterial with a negative refractive index, Nature 455, 376 379 (8). 1. J. Zhou, T. Koschny, M. Kafesaki, and C. M. Soukoulis, Negative refractive index response of weakly and strongly coupled optical metamaterials, Phys. Rev. B 8, 3519 (9). 13. D. O. Guney, Th. Koschny, and C. M. Soukoulis, Reducing ohmic losses in metamaterials by geometric tailoring, Phys. Rev. B 8, 1519 (9). (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 1688

14. S. A. Ramakrishna and J. B. Pendry, Removal of absorption and increase in resolution in a near-field lens via optical gain, Phys. Rev. B 67, 111 (3). 15. N. M. Lawandy, Localized surface plasmon singularities in amplifying media, Appl. Phys. Lett. 85, 54 54 (4). 16. M. A. Noginov, G. Zhu, M. Bahoura, J. Adegoke, C. E. Small, B. A. Ritzo, V. P. Drachev, and V. M. Shalaev, Enhancement of surface plasmons in an Ag aggregate by optical gain in a dielectric medium, Opt. Lett. 31, 3 34 (6). 17. T. A. Klar, Negative-index metamaterials: Going optical, IEEE J. Sel. Top. Quantum Electron. 1, 116 1115 (6). 18. A. K. Sarychev and G. Tartakovsky, Magnetic plasmonic metamaterials in actively pumped host medium and plasmonic nanolaser, Phys. Rev. B 75, 85436 (7). 19. Y. Sivan, S. Xiao, U. K. Chettiar, A. V. Kildishev, and V. M. Shalaev, Frequency-domain simulations of a negative-index material with embedded gain, Opt. Express 17, 46 474 (9).. A. D. Boardman, Yu. G. Rapoport, N. King, and V. N. Malnev, Creating stable gain in active metamaterials, J. Opt. Soc. Am. B 4, A53 A61 (7). 1. A. N. Lagarkov, V. N. Kisel, and A. K. Sarychev, Loss and gain in metamaterials, J. Opt. Soc. Am. B 7, 648 659 (1).. M. Wegener, J. Luis Garca-Pomar, C. M. Soukoulis, N. Meinzer, M. Ruther, and S. Linden, Toy model for plasmonic metamaterial resonances coupled to two-level system gain, Opt. Express 16, 19785 19798 (8). 3. A. Fang, Th. Koschny, M. Wegener, and C. M. Soukoulis, Self-consistent calculation of metamaterials with gain, Phys. Rev. B 79, 4114 (9). 4. A. Fang, Th. Koschny, and C. M. Soukoulis, Lasing in metamaterial nanostructures, J. Opt. 1, 413 (1). 5. A. Fang, Th. Koschny, and C. M. Soukoulis, Self-consistent calculations of loss-compensated fishnet metamaterials, Phys. Rev. B 8, 111 (1). 6. S. Wuestner, A. Pusch, K. L. Tsakmakidis, J. M. Hamm, and O. Hess, Overcoming Losses with Gain in a Negative Refractive Index Metamaterial, Phys. Rev. Lett. 15, 1741 (1). 7. N. I. Zheludev, S. L. Prosvirnin, N. Papasimakis, and V. A. Fedotov, Lasing spaser, Nat. Photonics, 351 354 (8). 8. D. J. Bergman and M. I. Stockman, Surface plasmon amplification by stimulated emission of radiation: quantum generation of coherent surface plasmons in nanosystems, Phys. Rev. Lett. 9, 74 (3). 9. M. I. Stockman, Spasers explained, Nat. Photonics, 37 39 (8). 3. K. Tanaka, E. Plum, J. Y. Ou, T. Uchino, and N. I. Zheludev, Multifold Enhancement of Quantum Dot Luminescence in Plasmonic Metamaterials, Phys. Rev. Lett. 15, 743 (1). 31. N. Meinzer, M. Ruther, S. Linden, C. M. Soukoulis, G. Khitrova, J. Hendrickson, J. D. Olitsky, H. M. Gibbs, and M. Wegener, Arrays of Ag split-ring resonators coupled to InGaAs single-quantum-well gain, Opt. Express 18, 414 4151 (1). 3. S. Xiao, V. P. Drachev, A. V. Kildishev, X. Ni, U. K. Chettiar, H.-K. Yuan, and V. M. Shalaev, Loss-free and active optical negative-index metamaterials, Nature 466, 735 738 (1). 33. A. E. Siegman, Lasers (Hill Valley, 1986), Chaps., 3, 6, and 13. 34. A. Taflove, Computational Electrodynamics: The Finite Difference Time Domain Method (Artech House, London, 1995). See Chaps. 3, 6, and 7. 35. D. R. Smith, S. Schultz, P. Markoŝ, and C. M. Soukoulis, Determination of effective permittivity and permeability of metamaterials from reflection and transmission coefficients, Phys. Rev. B 65, 19514 (). 36. Th. Koschny, P. Markoŝ, E. N. Economou, D. R. Smith, D. C. Vier, and C. M. Soukoulis, Impact of inherent periodic structure on effective medium description of left-handed and related metamaterials, Phys. Rev. B 71, 4515 (5). 37. N. Katsarakis, T. Koschny, M. Kafesaki, E. N. Economou, and C. M. Soukoulis, Electric coupling to the magnetic resonance of split ring resonators, Appl. Phys. Lett. 84, 943 945 (4). 38. D. Schurig, J. J. Mock, and D. R. Smith, Electric-field-coupled resonators for negative permittivity metamaterials, Appl. Phys. Lett. 88, 4119 (6). 1. Introduction The field of metamaterials has seen spectacular experimental progress in recent years [1 7]. However, huge intrinsic losses in the metal-based structures have become the major obstacle towards real world applications at optical wavelengths. Generally, losses are orders of magnitude too large for the envisioned applications, such as, perfect lenses [8], and invisibility cloaking [9]. Achieving such reduction of losses by geometric tailoring of the metamaterial designs [1 13] appears to be out of reach. So far, the most promising and generic approach is to incorporate gain material into metamaterial designs. One important issue is not to assume (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 1689

the metamaterial structure and the gain medium are independent from one another [14 1]. In fact, increasing the gain in the metamaterial changes the metamaterial properties, which, in turn, changes the coupling to the gain medium until a steady state is reached. So, there is a need for self-consistent calculations [ 6] for incorporating gain materials into realistic metamaterials. Instead of simply forcing negative imaginary parts of the local gain material s response function, which produces strictly linear gain, the self-consistent approach inherently includes the nonlinearity and gain saturation of the gain material. When the signal amplitudes are very weak and in the linear region of gain, the linear model can be obtained in the self-consistent approach by assuming a constant population inversion in Eq. (1), i.e., a population inversion given by an average field. However, for strong signals, it is necessary to have self-consistent calculations. Time-domain self-consistent calculations of gain incorporated into D magnetic metamaterials [3, 4] and 3D realistic fishnet structures [5, 6] have been recently reported. Results have shown the magnetic resonances of the D split-ring resonators (SRRs) and the fishnet structures can be substantially undamped by the gain material. Hence, the losses of the magnetic susceptibility, μ, are compensated. It is demonstrated the gain medium can give an effective gain much larger than its bulk counterpart due to the strong local-field enhancement inside the metamaterial designs [3, 7 9]. Recent experimental works also report loss compensations in metamaterial nanostructures coupled with quantum dots [3], single quantum wells [31] and organic dyes [3]. In this paper, we apply a detailed 3D self-consistent computational scheme to study the optical response of a realistic SRR array with a gain layer underneath. In section, we present the semi-classical theory of lasing and describe in detail the computational approach. In section 3, we present the geometric dimensions of the SRR array with gain. In section 4, we study the loss compensations of the combined SRR-gain system for three different pumping schemes on the gain layer: (1) the gain is isotropic and pumped with a homogeneous pumping rate, () the gain is isotropic but has a shadow cast by the SRR where the gain is away and (3) the gain is anisotropic, i.e., it is only pumped in one selected direction. In addition, we investigate the effect of the gain dielectric background on the loss compensation. In section 5, we present our conclusions.. Theory and model The gain atoms are embedded in the host medium and described by a generic four-level atomic system, which tracks fields and occupation numbers at each point in space, taking into account energy exchange between atoms and fields, electronic pumping, and non-radiative decays [33]. The two-level system is not taken because in reality it can not achieve the population inversion required for gain and lasing due to the de-exciting processes of spontaneous and stimulated emssions. The four-level system is more efficient in achieving the population inversion and most practical gain media can be modeled by the system of this type. An external mechanism pumps electrons from the ground state level, N, to the third level, N 3, at a certain pumping rate, Γ pump, proportional to the optical pumping intensity in an experiment. After a short lifetime, τ 3, electrons transfer non-radiatively into the metastable second level, N. The second level (N ) and the first level (N 1 ) are called the upper and lower lasing levels. Electrons can be transferred from the upper to the lower lasing level by spontaneous and stimulated emissions. At last, electrons transfer quickly and non-radiatively from the first level (N 1 ) to the ground state level (N ). The lifetimes and energies of the upper and lower lasing levels are τ 1, E and τ 1, E 1, respectively. The center frequency of the radiation is ω a =(E E 1 )/ h, chosen to equal π 1 14 rad/s. The parameters, τ 3, τ 1, and τ 1, are chosen 5 1 14,5 1 1, and 5 1 14 s, respectively. The total electron density, N (t = ) =N (t)+n 1 (t)+n (t)+n 3 (t) =5. (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 169

1 3 /m 3, and the pumping rate, Γ pump, is an external parameter. These gain parameters are chosen to overlap with the resonance of the split-ring resonator. The time-dependent Maxwell equations are given by E = B/ t and H = εε o E/ t + P/ t, where B = μμ o H and P is the dispersive electric polarization density from which the amplification and gain can be obtained. Following the single electron case, we can show [33] the polarization density P(r,t) in the presence of an electric field obeys locally the following equation of motion, P(t) t + Γ a P(t) t + ω a P(t) = σ a ΔN(t)E(t), (1) where Γ a is the linewidth of the atomic transition ω a and is equal to π 1 1 rad/s. The factor, ΔN(r,t) =N (r,t) N 1 (r,t), is the population inversion that drives the polarization, and σ a is the coupling strength of P to the external electric field and its value is taken to be 1 4 C /kg. It follows [33] from Eq. 1 that the amplification line shape is Lorentzian and homogeneously broadened. The occupation numbers at each spatial point vary according to N 3 t N t N 1 t N t = Γ pump N N 3 τ 3, = N 3 τ 3 + 1 hω a E P t N τ 1, = N τ 1 1 hω a E P t N 1 τ 1, = N 1 τ 1 Γ pump N, (a) (b) (c) (d) where hω 1 a E P t is the induced radiation rate or excitation rate depending on its sign. To solve the behavior of the active materials in the electromagnetic fields numerically, the finite-difference time-domain (FDTD) technique is utilized [34]. In the FDTD calculations, the discrete time and space steps are chosen to be Δt =. 1 18 s and Δx =.5 1 9 m. The initial condition is that all the electrons are in the ground state, so there is no field, no polarization, and no spontaneous emission. Then, the electrons are pumped from N to N 3 (then relaxing to N ) with a constant pump rate, Γ pump. The system begins to evolve according to the system of equations above. 3. Geometric dimensions of the SRR array As shown in Fig. 1(a), the SRR is fabricated on a GaAs-gain-GaAs sandwich substrate. It is made from silver with its permittivity given by a Drude model, ε(ω) =1 ω p/(ω + iωγ), where ω p = 1.37 1 16 rad/s and γ =.73 1 13 rad/s. The GaAs layer between the SRR and gain is introduced to avoid the quenching effect. The incident wave propagates along the y direction parallel to the SRR plane and has the magnetic field perpendicular to that plane. The unit cell size along the propagation direction is a. Inz direction, the unit cell size is h, which is larger than h 1 + h + h 3 + h s, where h 1, h, h 3 and h s are the thicknesses of the bottom GaAs layer, the gain layer, the GaAs spacing layer, and the SRR, respectively. Along the unit cell boundaries in x and z directions, periodic boundary conditions are enforced to simulate the infinite periodic structure. All the dimensions are chosen to have the magnetic resonance overlap with the emission frequency of 1 THz of the gain material. For comparison, we also introduce another gain configuration (see Fig. 1(b)), where the gain is embedded in the gap of the SRR instead of a layer underneath. The dimensions are kept the same as Fig. 1(a). (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 1691

H h 4 h h 3 h 1 a l h s w d l a k E h z x y (a) (b) Fig. 1. (Color online) (a) One unit cell for the silver-based SRR structure (light blue) with the gain layer underneath. The dielectric constants ε for GaAs (yellow) and gain (red) are 11 and, respectively. The whole structure is in vacuum background (light gray). The dimensions are a = 5nm, l = 16nm, h = 8nm, h 1 = 15nm, h = h 3 = 1nm, h 4 = 45nm, h s = 5nm, w = 4nm and d = nm. (b) same as (a) except the gain is embedded in the SRR gap with ε = 1 and the gain layer in (a) is replaced by a dielectric layer (ε = ) (blue). 4. Numerical simulations and discussions In this section, we apply the three pumping schemes discussed in section 1 on the gain layer. The linewidths of the magnetic resonances for different pumping rates are investigated to see if the gain can effectively reduce the magnetic losses. We also do simulations for different gain background dielectric constants to see how it affects the loss compensation. 4.1. Isotropic gain We first let a wide band Gaussian pulse of a given amplitude go through one layer of the SRR structure shown in Fig. 1(a) and calculate the transmission T, the reflection R, and the absorption A = 1 T R, as a function of frequency in the propagation direction. With the introduction of gain, the absorption near the resonance frequency f = 1 THz decreases and the transmission increases. To investigate the loss reduction of the magnetic resonators, we plot the retrieved effective permeabilities, μ, without and with gain by inverting the scattering amplitudes [35, 36] in Fig. (a). One can see the gain undamps the magnetic resonance of the SRR and the resonant effective permeability μ of the SRR becomes much stronger and narrower compared to the case without gain. In Fig. (b), we plot the effective permeabilities, μ, without (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 169

Permeability μ 8 6 4 Re(μ) w/o gain Im(μ) w/o gain Re(μ) with gain Im(μ) with gain Permeability μ 8 6 4 Re(μ) w/o gain Im(μ) w/o gain Re(μ) with gain Im(μ) with gain 9 95 1 15 11 115 (a) 9 95 1 15 11 115 (b) Fig.. (Color online) The retrieved results for the real and imaginary parts of the effective permeability μ, without and with gain, for two different gain configurations. (a) the gain is underneath the SRR as shown in Fig. 1(a). The pumping rate is Γ pump = 1. 1 9 s 1. (b) the gain is in the SRR gap as shown in Fig. 1(b). The pumping rate is Γ pump = 7. 1 8 s 1. I/(η ω E) (a.u.) 5 x 1 13 No gain Γ 4 pump =1. 1 9 s 1 Γ pump =1.5 1 9 s 1 3 Γ =. 1 9 s 1 pump Γ =.8 1 9 s 1 pump 1 1 6.4 THz.5 THz 3.5 THz 4. THz 4.9 THz I/(η ω E) (a.u.) 5 x 1 13 No gain Γ 4 pump =1. 1 9 s 1 Γ pump =1.5 1 9 s 1 3 1 1 6.4 THz.5 THz 3.7 THz 9 95 1 15 11 (a) 9 95 1 15 11 (b) Fig. 3. (Color online) The real (solid) and imaginary (dashed) parts of I/(ηω E) as a function of frequency for different pumping rates. (a) For the structure with a gain layer below the SRR, shown in Fig. 1(a). (b) For the structure with the gain in the SRR gap as shown in Fig. 1(b). Notice that the resonance is getting stronger and narrower as the pumping rate increases. and with gain for the case the gain is in the SRR gap. Similar to the results for D SRR in Ref. 3, the weak and broad resonant μ becomes strong and narrow with the introduction of gain in the SRR gap. Note that a lower pumping rate (Γ pump = 7. 1 8 s 1 ) leads to a sharper magnetic resonance comparing with the case the gain is underneath the SRR (Γ pump = 1. 1 9 s 1 ) due to the local electric field concentration in the gap. However, the strong magnetic resonances in Figs. (a) and (b) are not symmetric due to the periodicity effect [36]. This asymmetry causes the difficulty in obtaining the linewidth of the magnetic resonance. The periodicity effect itself is inherent in the retrieval procedure. To distinguish the magnetic resonance of the SRR from the periodicity effect of the structure, we directly calculate the resonant current (i.e., the magnetic moment) flowing around the split ring, without going through the retrieval procedure. Consider (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 1693

the SRR as a simple LCR circuit model, we can have the following equation, L di Idt dt + C + IR = ε emf, (3) where L, C and R are the effective inductance, capacitance and resistance of the SRR, respectively, and I is the current flowing in the SRR and ε emf is the induced electromotive force. From Faraday s law, ε emf = dφ/dt = iaμ ωh = ia c ω E.(Φis the magnetic flux through the SRR, A is the area enclosed by SRR, and c is the speed of light in vacuum.) Then we can obtain the expression with Lorentz resonance shape, I ηω E = 1 ω ω (4) + iγω, where η, ω, and γ are A/(cL), 1/ LC, and R/L, respectively. The detailed results are plotted in Fig. 3(a) for the structure with the gain layer underneath. One can see the current resonances have very nice Lorentz line shapes. As the pumping rate increases, the resonance is getting stronger and narrower. The full width at half maximum (FWHM) reaches.5 THz when the pumping rate Γ pump =.8 1 9 s 1, which is a significant loss reduction compared with the FWHM without gain (FWHM = 6.4THz). So the gain compensates the losses. In addition, we also calculate I/(ηω E) vs. frequency for the structure with gain in the SRR gap to compare the efficiency of the loss compensation for these two different gain configurations. The results are shown in Fig. 3(b). One can see the structure with gain in the SRR gap needs less gain (i.e., smaller pumping rate 1.5 1 9 s 1 ) to reach the same FWHM,.5THz, of the resonance than the case with gain underneath the SRR with the pumping rate Γ pump =.8 1 9 s 1. It is easy to understand the difference in pumping rates in the two designs because of the strong local electric field enhancement in the SRR gap. Though the loss compensation for the structure with a gain layer underneath is not so efficient as the case with the gain in the SRR gap, the results in Fig. 3(a) still show that the magnetic losses can be substantially reduced, especially if we push the pumping rate to a high value. Q' Q 3.9 THz No gain Γ pump = 1. 1 9 s 1 E Q' Q' Q' Q k I Permittivity ε 1 1 3. THz Q' Q 94 96 98 1 1 14 16 (a) (b) Fig. 4. (Color online) (a) Schematic of the excitation of the magnetic resonance by the incident electric field E parallel to the gap bearing side of the SRR. A circular current appears due to the different configuration of surface charges on both sides of the SRR. (b) The retrieved results for the real (solid) and imaginary (dashed) parts of the effective permittivity ε, with and without gain, for the normal incidence in Fig. 1(a). For the case with gain, the pumping rate Γ pump = 1. 1 9 s 1. It is experimentally difficult to have the parallel incidence for such a planar structure like the SRR array. In experiments, the SRR plane is oriented perpendicular to the incidenct wave and (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 1694

its gap bearing side is parallel to the incident electric field. Hence, the electric field can couple to the electric dipole in the gap and induce the magnetic resonance [37, 38] (see Fig. 4(a)). Simulations are done for this case to see if the losses can be compensated by the gain layer underneath. With this incidence direction, the unit cell size in the propagation direction is h, which is much smaller than the wavelength λ (λ/h = 37.5), so the resonance is far below the Brillouin zone edge and we can ignore periodicity effect. Figure 4(b) plots the retrieved effective permittivity ε, with and without gain. Both of them have a very nice Lorentz line shape. Without gain, the resonance is broad and weak, and the FWMH is 3THz. With the introduction of gain, the resonance becomes stronger and narrower, and the FWHM reduces to a much smaller value,.9 THz. So the gain compensates the losses of the SRR for perpendicular incidence. 4.. Isotropic gain with a shadow of the SRR y k z x E H Fig. 5. Top view of the gain layer in Fig. 1(a) when a shadow (blue) is cast by the SRR structure. The gain does not exist in the shadow, while in other area (red) it is homogeneously pumped. So far, the gain material in our simulations is isotropic and pumped by a homogeneous pumping rate Γ pump. This is an ideal case. Consider the case in experiments that we incident an external optical pumping wave on the structure (Fig. 1(a)) from the top to optically pump the electrons from level to level 3, there will be a shadow on the gain layer cast by the SRR structure, where the gain is pumped by a much lower rate. As a simplified model, we turn off the gain in the area which lies directly under the SRR to simply emulate the shadow of the SRR structure, while we still keep a homogeneous pumping rate Γ pump in other gain area (see Fig. 5). In Fig. 6(a), we plot I/(ηω E) as a function of frequency in this case. Compared with the case without the shadow on the gain layer (Fig. 3(a)), the resonance gets much weaker and broader (FWHM = 5.7THz and 5.4THz for the pumping rates Γ pump = 1. 1 9 s 1 and 1.5 1 9 s 1, respectively). This shows the gain in the shadow area plays an important part in the loss compensation. 4.3. Anisotropic gain The gain in our simulations discussed above is isotropic, which is equally pumped in all directions. The realistic gain, such as semiconductor quantum dots/wells, can be anisotropic, i.e., it can only couple to the external field in a certain direction. Since the electric fields in the SRR structure are mainly distributed across the gap, we have the active direction of the gain material parallel to y direction in Fig. 1(a), i.e., the gap bearing side of the SRR. So, the gain only couples to the electric field in y direction. The corresponding I/(ηω E) vs. frequency curves for different pumping rates are plotted in Fig. 6(b). One can see the resonances are also much (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 1695

I/(η ω E) (a.u.).5 x 1 13 No gain Γ =1. 1 9 s 1 pump Γ pump =1.5 1 9 s 1 1.5 1.5.5 5.4 THz 6.4THz 5.7 THz I/(η ω E) (a.u.).5 x 1 13 No gain Γ pump =1. 1 9 s 1 Γ =1.5 1 9 s 1 pump 1.5 1.5.5 5.9THz 6.4 THz 6.1 THz 1 1 9 95 1 15 11 (a) 9 95 1 15 11 (b) Fig. 6. (Color online) The real (solid) and imaginary (dashed) parts of I/(ηω E) as a function of frequency for different pumping rates. (a) For the structure with a shadow on the gain layer cast by the SRR. (b) For the structure where the gain is pumped in y direction only. broader than the case with homogeneously pumped isotropic gain. So the loss compensation is less efficient. 4.4. Explanation of the differences among the loss compensations by the three pumping schemes V/m y x (a) (b) (c) Fig. 7. (Color online) The electric field amplitude distribution at the resonance frequency in the cross-section of the gain layer (z = nm from the bottom of the structure) for different components: (a) E x, (b) E y and (c) E z. The area enclosed by the white line is the projection of the SRR on the gain layer. The electric fields are calculated without gain. To see why these three gain pumping schemes are so different on the loss compensation, we have calculated the electric field amplitude distribution in the cross-section of the gain layer (xy plane in Fig. 1(a)). The detailed results are plotted in Figs. 7(a) 7(c). One can see the x component of electric field, E x, is very weak while the other two components, E y and E z, are relatively strong. So we can ignore the gain contribution by E x and focus on the gain from the coupling with E y and E z. Notice that E y is bounded in the area right below the SRR gap (Fig. 7(b)) while E z mainly has a significant value in the projection of the SRR on the gain layer (Fig. 7(c)). This characteristic of the field amplitude distribution leads to almost no contribution (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 1696

by E z when we have a shadow in the gain layer since there is no gain in that area. Similarly, the gain contribution by E z goes away for the anisotropic gain because the gain only couples with the y component of the electric field, E y. This fact explains the big difference between the homogeneously pumped isotropic gain and the other two gain pumping schemes. 4.5. The effect of the dielectric background of gain Im(I/ηω E) (a.u.) 3.5 x 1 13 ε = w/o gain g 3 ε g = with gain ε g =5 w/o gain.5 ε = 5 with gain g ε =11 w/o gain g ε = 11 with gain g 1.5 1.5 75 8 85 9 95 1 15 11 Fig. 8. (Color online) The imaginary parts of I/(ηω E) as a function of frequency for different background dielectric constants of the gain material, which is only pumped in y direction. For the case with gain, the pumping rate is Γ pump = 1.5 1 9 s 1. Note the resonance enhancements by the gain are almost the same. Since there is a high contrast between the dielectric constants of the GaAs (ε = 11) and gain (ε = ) layers, the electromagnetic fields may be bounded in the high dielectric layer. In this section, we will discuss the effect of the dielectric background of gain on the loss compensation. In Fig. 8, we plot the detailed results for the imaginary parts of I/(ηω E) as a function of frequency, with and without gain, for the background dielectric constants of the gain layer ε g =, 5 and 11. The gain is anisotropic and only couples to the electric field in y direction. We can see the resonance frequency shifts down as the dielectric constant increases. This is expected since the effective capacitance increases with the increment of the dielectric constant. To effectively compensate the losses, we scale the emission frequency to overlap with the corresponding resonance frequencies and then pump with the same rate Γ pump = 1.5 1 9 s 1.We can see from Fig. 8 the resonance enhancements are almost the same for different background dielectric constants of the gain. To explain this phenomenon, we plot the electric field amplitude distributions in a plane crossing the middle of the gap bearing side of the SRR (Fig. 9), for ε g =, 5 and 11, respectively. The E x component is ignored since it is very weak as shown in Fig. 7(a). From Fig. 9(a), we can see the field amplitude distribution of E y, the only component which couples to the gain, does not change much in the gain layer as the gain background dielectric constant changes. Although there is a bounding effect on the fields, the y component of the electric field, E y, does not substantially decay in such a very narrow gain layer (1 nm) neighboring to the high dielectric GaAs layer. The main change in the electric field is the z component of the electric field, E z, decreases in the gain layer as the gain background dielectric constant, ε g, increases, as shown in Fig. 9(b). This is due to the continuity of the normal component of the electric displacement across the interface since there is no free charge accumulation. Hence the normal component of the electric field is inversely proportional to the dielectric constant. The change of E z does not affect the loss compensation due to no coupling between the gain and E z. If the gain can couple (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 1697

E y z = E = g z g y V/m g =5 g =5 g =11 g =11 (a) (b) Fig. 9. (Color online) The electric field amplitude distributions for different background dielectric constants of the gain material (ε g =,5 and 11) at their corresponding resonance frequencies, in a yz plane crossing the middle of the gap bearing side of the SRR. (a) E y and (b) E z. The area enclosed by the dashed green line indicates the position of the gain layer. The electric fields are calculated without gain. to E z, such as the isotropic gain, the background dielectric constant of the gain will significantly affect the loss compensation. 5. Conclusions We have numerically studied the loss compensation of the silver-based SRR structure with a gain layer underneath. Numerical results show that the losses of the SRR can be compensated by the gain layer for both the parallel and perpendicular incidences. Three different gain pumping schemes are applied in the simulations and the efficiencies of their corresponding loss compensations are studied by investigating the linewidth of the resonant current. The homogeneously pumped isotropic gain can significantly reduce the magnetic losses, though it is less efficient in the loss compensation compared to the case with the gain in the SRR gap. The other two schemes, (1) a homogeneously pumped isotropic gain with a shadow cast by the SRR and () anisotropic gain only coupled to E y, the electric field component parallel to the gap bearing side of the SRR, are much less efficient in the loss compensation compared to the isotropic gain case, due to no interactions between the electric field perpendicular to the SRR plane, E z, and the gain in these two schemes. We have also studied the effect of the background dielectric of gain. In a very narrow gain layer, the gain dielectric background mainly affects the electric field perpendicular to the GaAs-gain interface due to the continuity of the normal component of the electric displacement across the interface. So, the dielectric background of gain does not make much difference for the gain pumped in the parallel direction only. Acknowledgments Work at Ames Laboratory was supported by the Department of Energy (Basic Energy Sciences) under Contract No. DE-AC-7CH11358. This work was partially supported by the European Community FET project PHOME (Contract No. 1339) and by Laboratory-Directed Research and Development Program at Sandia National Laboratories. The author Z. Huang (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 1698

gratefully ackowledges support of the National Natural Science Foundation of China (Grant No. 6931). (C) 11 OSA June 11 / Vol. 19, No. 13 / OPTICS EXPRESS 1699