Simulations and realizations of active right-handed metamaterials with negative refractive index

The theory of determining the sign of the refractive index in active materials is discussed. Animations of numerical simulations are presented, supporting the claim that negative refractive index may occur in right-handed media. An example of such a medium, in the form of a lumped circuit model with active and passive resonances, is presented. © 2007 Optical Society of America OCIS codes: (160.4670) Optical materials; (160.4760) Optical properties; (260.2065) Effective medium theory; (260.2110) Electromagnetic optics; (350.4010) Microwaves; (260.2030) Dispersion. References and links 1. V. G. Veselago, “The electrodynamics of substances with simultaneously negative ε and μ ,” Sov. Phys. Usp. 10(4), 509–514 (1968). 2. J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, “Low frequency plasmons in thin-wire structures,” J. Phys.: Condens. Matter 10(22), 4785–4809 (1998). 3. J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, “Magnetism from conductors and enhanced nonlinear phenomena,” IEEE Trans. Microwave Theory Tech. 47(11), 2075–2084 (1999). 4. D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, “Composite medium with simultaneously negative permeability and permittivity,” Phys. Rev. Lett. 84(18), 4184–4187 (2000). 5. G. V. Eleftheriades, A. K. Iyer, and P. C. Kremer, “Planar negative refractive index media using periodically L-C loaded transmission lines,” IEEE Trans. Microwave Theory Tech. 50(12), 2702–2712 (2002). 6. 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(20), 201101 (2003). 7. 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, 3022 (2006). 8. A. K. Popov and V. M. Shalaev, “Compensating losses in negative-index metamaterials by optical parametric amplification,” Opt. Lett. 31, 2169 (2006). 9. V. M. Shalaev, “Optical negative-index metamaterials,” Nat. Photonics 1, 41–48 (2007). 10. Y.-F. Chen, P. Fischer, and F. W. Wise, “Negative refraction at optical frequencies in nonmagnetic two-component molecular media,” Phys. Rev. Lett. 95(6), 067402 (2005). 11. J. Skaar, “Fresnel equations and the refractive index of active media,” Phys. Rev. E 73, 026605 (2006). 12. Y.-F. Chen, P. Fischer, and F. W. Wise, “Sign of the refractive index in a gain medium with negative permittivity and permeability,” J. Opt. Soc. Am. B 23, 45–50 (2006). 13. T. G. Mackay and A. Lakhtakia, “Comment on “Negative refraction at optical frequencies in nonmagnetic twocomponent molecular media”,” Phys. Rev. Lett. 96(15), 159701 (2006). 14. Y.-F. Chen, P. Fischer, and F. W. Wise, “Chen, Fischer, and Wise reply,” Phys. Rev. Lett. 96, 159702 (2006). 15. S. A. Ramakrishna, “Comment on “Negative refraction at optical frequencies in nonmagnetic two-component molecular media”,” Phys. Rev. Lett. 98(5), 059701 (2007). 16. Y.-F. Chen, P. Fischer, and F. W. Wise, “Chen, Fischer, and Wise reply,” Phys. Rev. Lett. 98(5), 059702 (2007). 17. S. A. Ramakrishna and O. J. F. Martin, “Resolving the wave vector in negative refractive index media,” Opt. Lett. 30, 2626 (2005). #85000 $15.00 USD Received 9 Jul 2007; revised 14 Aug 2007; accepted 14 Aug 2007; published 15 Aug 2007 (C) 2007 OSA 20 August 2007 / Vol. 15, No. 17 / OPTICS EXPRESS 10935 18. A. N. Grigorenko, “Negative refractive index in artificial metamaterials,” Opt. Lett. 31, 2483 (2006). 19. J. Skaar, “On resolving the refractive index and the wave vector,” Opt. Lett. 31, 3372 (2006). 20. L. Brillouin, Wave propagation and group velocity (Academic Press, New York and London, 1960). 21. L. D. Landau and E. M. Lifshits, Electrodynamics of continuous media, chap. IX.62 (Pergamon Press, 1960). 22. In general, we can set γ = 0+ provided the refractive index is analytic in the upper half-plane and the denominator in the Fresnel equations is nonzero in the upper half-plane. Excluding media with absolute instabilities, the refractive index can always be identified as an analytic function in the upper half-plane. 23. Poles in the upper half-plane mean that the slab will start lasing. 24. T. Koschny, P. Markos, 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(24), 245105 (2005). 25. B. Nistad and J. Skaar, in Photonic metamaterials: From random to periodic, V. M. Shalaev and A. Genack, eds. (OSA, 2006). ISBN:1-55752-808-X. 26. J.-S. Lee and Y.-S. Kwon, “Negative resistance circuit for monolithic resonators using gate-to-source resistive feedback,” Electron. Lett. 34(18), 1758–1760 (1998). 27. L. O. Chua, J. Yu, and Y. Yu, “Bipolar-JFET-MOSFET negative resistance devices,” IEEE Trans. Circuits Syst. 31(1), 46–61 (1985). 28. D. M. Pozar, Microwave engineering, chap. 4, 2nd ed. (Wiley, 1998). 29. M. Born and E. Wolf, Principles of Optics, chap. 1.6.5, pp. 66–67, 3rd ed. (Pergamon Press, 1965).


Introduction
In 1968 Veselago predicted that isotropic materials with simultaneously negative permittivity ε and permeability μ may yield negative refraction [1].Veselago referred to such materials as left-handed due to the fact that the electric field, the magnetic field, and the wave vector constitute a left-handed set of vectors.Thus, the phase velocity and Poynting vector point in opposite directions, in contrast to the situation in right-handed materials.Left-handed materials do not exist naturally; however, it is possible to make artificial metamaterials with negative ε and μ in the same frequency range [2,3,4,5].
Passive metamaterials often suffer from limited performance due to large losses.It has therefore been suggested to combine them with active gain media [6], or even create active metamaterials [7].Nonlinear processes such as parametric amplification have also been suggested as a method for overcoming inherent losses in metamaterials [8].As the interest in active metamaterials seem to be increasing [9], it is crucial to understand fundamentally how they behave electromagnetically.
Of particular interest are the sign of the refractive index (direction of phase velocity), and the direction of the energy flow.The sign of the refractive index depends on the location of zeros and poles of εμ in the complex ω-plane.This indicates that it is possible to obtain negative refractive index without any magnetic resonances, μ = 1, but instead two electric resonances.Indeed, it has been suggested that certain nonmagnetic media, with active and passive dielectric resonances, can exhibit negative refraction [10].In such right-handed materials the wave vector and Poynting vector both point towards the source [11].
Recently, these results have been discussed in the literature [11,12,13,14,15,16].While the conclusions differ, all authors seem to argue that causality must determine the sign of the refractive index.However, the interpretations of causality differ and in fact, give different results for certain active media [10,11,12,13,14,15,16,17,18,19].To identify the correct interpretation and resolve this controversy, it is necessary to go back to first principles: Causality means that if an excitation starts at time t = 0, no effect can be detected a distance z away before t = z/c.In other words, the front of an electromagnetic wave cannot travel faster than c, the vacuum velocity of light [20].Using this principle, it has been shown analytically how to determine the sign of the refractive index in active media [11].
In this work, we will first present animations of numerical simulations, to demonstrate negative refractive index in right-handed, active media [10,11].The electromagnetic solution in a slab of thickness d can be found unambiguously as it only dependents on ε, μ, and d, and is independent of the refractive index.Moreover, by causality, the time-domain solution for t < d/c must be identical to that of a semi-infinite medium when the slab is excited at one side, starting at t = 0. Letting d be sufficiently large, the solutions can be compared in an arbitrarily large time window.The movies clearly support the claim that negative refractive index may be obtained in active, right-handed media.
To concretize this class of media, we propose a 1-D transmission line model with lumped circuit elements, implementing a right-handed medium with negative refractive index.This model can be used as a starting point for analysis and possible implementation.Although we consider a one-dimensional structure, generalization to two dimensions is straightforward.

Determining the sign of the refractive index
When a plane wave is normally incident from vacuum to a semi-infinite medium, the reflection R and transmission S are given by the Fresnel equations: Here k = ωn/c, η = μ/n, and n 2 = εμ; ε and μ are the permittivity and permeability of the semi-infinite medium.The interface between vacuum and the semi-infinite medium is the plane z = 0; the region z > 0 is the location of the semi-infinite medium.For convenience, the propagation factor exp(iωnz/c) is included in S. The sign of the refractive index n must be identified to ensure causality; for passive materials this can be achieved simply by requiring that the Poynting vector point in +z-direction, or that the wave decays in the +z-direction.As an example of an active medium, we consider a material with ε(ω) = (1 + f (ω)) 2 and μ(ω) = 1, where We take F = 2 and Γ = 0.005ω 0 .This material is causal and realizable as ε is analytic in the upper half-plane of complex frequency and ε(ω Using the method of Refs.[20,10,11,12,19], we find the refractive index n 1 (ω) = 1 + f (ω).On the other hand, using the approaches in Refs.[13,15,17,18], we find the refractive index The real and imaginary parts of the refractive index for the two different solutions are plotted in Fig. 1.Note that there is a frequency band where Re n 1 (ω) < 0; thus the first solution predicts that both phase velocity and steadystate energy flow may point towards the source.For the other solution, Re n 2 (ω) > 0 for all frequencies; thus this solution means that the phase velocity and energy flow are directed away from the source.Note the two points, ω = ω 0 and ω = 1.73ω 0 , where n 2 (ω) is discontinuous.
To resolve the above controversy, we first compute the time-domain solutions associated with the choices n 1 (ω) and n 2 (ω).This is achieved with the inverse Laplace transform.We excite the medium by a plane wave, normally incident at z = 0 − , starting at t = 0.For example, if the excitation at z = 0 − is u(t) cos(ω 1 t), where u(t) is the unit step function and ω 1 is the excitation frequency, the time-domain field in the material can be expressed as Here γ is a sufficiently large, positive number, such that Im ω < γ for all poles or nonanalytic points.For S given by Eq. (1b), assuming that the refractive index has meaning for real frequencies, the integral can be evaluated just above the real frequency axis (γ = 0 + ) [11,22].The excitation frequency is taken to be ω 1 = 1.4ω 0 ; at this frequency ε(ω 1 ) ≈ +1.13 − 0.32i.Now we can examine the solutions associated with n(ω) = n 1 (ω) (Fig. 2) and n(ω) = n 2 (ω) (Fig. 3), and judge whether they satisfy causality.Note that causality should be interpreted in its most fundamental form; no field can arise at z before the time t = z/c.We clearly see that the wave in Fig. 2 satisfies causality, while the wave in Fig. 3 does not.In fact, in Fig. 3 a nonzero field exists everywhere, even at t = 0.In Fig. 2, the wave front propagates at exactly the speed c [20].Similarly, when the excitation is turned off, the high-frequency components induced by the abrupt envelope propagate at c.Note that the field eventually dies out after the excitation has been turned off.Some ripple can be seen in front of the red line in Fig. 2.This is a numerical error, which relate to Gibbs' effect, and can be reduced if a broader frequency representation is used.
To see that n 1 (ω) not only gives a causal solution, but also the correct one, we compare the solution in Fig. 2 to the electromagnetic field in a slab of finite thickness d.Such a comparison makes sense since, by causality, the time-domain fields in a slab and in a semi-infinite medium must coincide for t < d/c.In an active slab, the fields may blow up with time and thus Fourier transformed fields may not exist.A natural remedy is to use Laplace transformed fields, i.e., to introduce complex frequencies with Im ω > 0. For a normally incident electromagnetic field, the reflected field R, the field in the slab S = S + exp(ikz) + S − exp(−ikz), and the transmitted field T can be found as [11 Note that R, S, and T are unchanged if n → −n; thus the sign of n is irrelevant in Eq. ( 4).The time-domain field in the slab is computed by an inverse Laplace transform.The parameter γ must be chosen such that the integration path is located above all poles of S in the upper halfplane [23].For the present choice of medium [ε(ω) = (1 + f (ω)) 2 and μ(ω) = 1], by choosing a sufficiently large d < ∞, it can be shown that there are no such poles.Then we can set γ = 0 + .We take the thickness to be ω 0 d/c = 35, which in this case means that S does not have poles in the upper half-plane.The resulting time-domain field is given in Fig. 4. We clearly see that the solution in Fig. 4 equals the solution in Fig. 2 when t < d/c.Actually, the solutions are rather similar also for later times, which can be attributed to the fact that the field is small at z = d.Again we note that the field dies out after the excitation has been turned off; thus the excitation of the amplified "backward" wave in Fig. 4 is fundamentally different to conventional lasing.
We conclude that n 1 (ω) is the correct refractive index function.The steady-state (monochromatic) solution, approached using the excitation u(t) cos(ω 1 t) in the limit t → ∞, demonstrates that the phase velocity points towards the source.Since μ = 1, also the energy flow points towards the source.This backward wave draws energy from the active medium.
By analytic arguments, one can arrive at the same conclusion [11,19].Provided εμ is continuous for real frequencies, one finds that the refractive index of (possibly) active media can be related to the permittivity and permeability using the formula where ϕ ε + ϕ μ is the complex argument of εμ, unwrapped so that it is continuous and tends to zero as ω → ∞.More generally, ε and μ may not be continuous for real frequencies (the Kramers-Kronig relations only imply analyticity in the upper half-plane Im ω > 0, not for Im ω = 0).Moreover, εμ may contain odd-order zeros in the upper half-plane, which means that n cannot be identified as an analytic function there.Excluding the latter possibility (which corresponds to materials with so-called absolute instabilities [11]), we find the refractive index from the analytic branch of √ εμ that tends to +1 as Im ω → ∞.The refractive index for real frequencies is the limit as Im ω → 0 + of this analytic branch.However, if εμ contains only a single discontinuity or zero for real frequencies (at ω = 0), we can still use Eq. ( 5), unwrapping in the frequency interval (0, ∞), and ensuring the limit ϕ ε + ϕ μ → 0 as ω → +∞.
The necessary phase unwrapping procedure means that the sign of the refractive index cannot be determined at a single frequency without knowing the global behavior of ε(ω)μ(ω) (at all frequencies).In other words, two media with ε = 1 − iα, α > 0 at a single frequency may have different refractive indices there.Thus the treatment of the example medium above does not contradict the fact that conventional gain media, with e.g.inverted Lorentzian responses, have positive real part of the refractive index.We note that any direct identification of n from ε and μ at a single frequency, such as in Refs.[13,15,17,18], are incorrect in general.
It may seem that we have used causality as an extra principle in addition to Maxwell's equations, to calculate the field.However, causality is in fact built into the time-domain versions of the Maxwell equations: Provided ε and μ describe a microscopically causal medium (i.e., the polarization and magnetization do not precede the electric and magnetic fields), an electromagnetic excitation cannot give any response before the time z/c, where z is the distance between the excitation point and the observation point.Thus, one can actually verify Fig. 2 by solving Maxwell's equations in the time-domain, with no extra assumptions.

Effective medium based on a transmission line structure
The example in the previous section represents a novel class of active metamaterials [10,11].To realize this type of material both passive and active resonances are needed.An example of a typical passive resonance that fits this model is the serial RLC-resonance circuit.The inverted Lorentz function has similar behavior as a parallel RLC-resonance including an active element, which points us towards a transmission line with lumped element inclusions.
In 2002 Eleftheriades et.al. [5] proposed a negative index medium based on periodically L-C loaded transmission lines.As long as the periodicity of the structure remains less than roughly λ 0 /30, where λ 0 is the vacuum wavelength [24], the structure can be described electromagnetically by an effective dielectric permittivity εε 0 [ F/ m] and an effective magnetic permeability μμ 0 [ H/ m], where ε 0 and μ 0 are the vacuum permittivity and permeability.Consider a transmission line as shown in Fig. 5, with series impedance per length unit Z, and shunt admittance per length unit Y .The telegrapher's equation for a general 1-D transmission line can be expressed as: Combining Eqs.(6a) and (6b) yields where β is the propagation constant.Mapping the voltage V to E y and the current I to −H x [5], Eqs. ( 6a) and (6b) can be rewritten as field equations: which yield the effective material parameters as The propagation constant satisfies It can easily be seen that interchanging the inductance and the capacitance of a conventional transmission line leads to negative ε and μ; thus the transmission line becomes left-handed.Since our goal is a right-handed transmission line medium with μ = 1, we will leave the impedance inductive, Z = −iωμ 0 , and rather change the admittance Y so that ε becomes a sum of one passive and one active resonance [10,25].An example circuit is given in Fig. 6, where the total admittance is Thus the effective relative permittivity and permeability become Here Λ, C l , and L l are the length, capacitance, and inductance for a section of conventional transmission line; L p , C p , and R p are the inductance, capacitance and resistance of the passive resonator; L a , C a , and R a are the inductance, capacitance and resistance of the parallel resonator part of the circuit; R n is the negative resistor element which provides gain.Negative resistance can be realized in several ways [26,27]; transistors and Gunn-diodes are feasible alternatives.The parallel resonance has zero admittance at the resonance frequency.This acts together with the gain element to give negative admittance and negative imaginary part of the permittivity.The inductance L n is included to remove a zero of ε from the positive imaginary frequency axis.
One must be careful when choosing the circuit parameters; it is not sufficient to have negative real part of the refractive index to obtain a medium similar to that in Section 2. It is important that ε does not have any odd-order zeros in the upper half plane, as this would create branch cuts in n.When there are such branch cuts, the refractive index loses its usual interpretation for real frequencies; moreover, the media are electromagnetically unstable (that is, with absolute instabilities, see [11]).
We have found that the following parameters will give negative refractive index in a small frequency window: , and C a = 10 pF. Figure 7 shows the refractive index versus frequency calculated by Eq. ( 5). Figure 8 shows the reflection coefficient calculated by Eq. (4a) using d = 20Λ.Negative real part of the refractive index exists in a small bandwidth from 1.11 GHz to 1.30 GHz, with the minimum value Re n = −0.48 at 1.20 GHz.A maximum 2.1 of the reflection coefficient is attained at a slightly larger frequency.At 1.22 GHz the permittivity has positive real part, ε = 0.0063 − 0.23i, corresponding to n = −0.34+ 0.33i and |R| = 1.9.The bandwidth where Re n < 0 coincides with the bandwidth where |R| > 1.

Periodic structure
The effective medium theory yields a useful design guide for realization of the effective parameters as a transmission line with discrete elements.To investigate such periodic structures exactly, we use the transmission matrix method.The transmission matrix of a transmission line where Λ, β 0 = ω √ L l C l , and Z 0 = L l /C l are the length, propagation constant, and characteristic impedance, respectively.For a general shunt admittance Y s , the transmission matrix is The transmission matrix for a unit cell is T uc = T tl T s , and the transmission matrix for N unit cells becomes Since det T uc = 1, the Nth power of T uc can be calculated analytically with Chebyshev's identity [29].T can now be converted into a scattering matrix [28] The scattering matrix element S 11 is the reflection coefficient and S 21 is the transmission coefficient.Note that S is symmetric due to reciprocity.
In the limit N → ∞ while d ≡ NΛ and Y s /Λ are fixed, the transmission matrix model should correspond exactly to the effective medium theory.Indeed, we may expand T uc to first order in Λ to obtain T uc = I + ϒΛ + O(Λ 2 ) = exp(ϒΛ) + O(Λ 2 ), where I is the identity matrix and Thus, in the above limit, T = exp(ϒd), and analytical expressions for the scattering parameters S 11 and S 21 can be found and verified to be equal to Eqs. (4a) and (4d).
To analyze the behavior of the periodic structure for finite Λ, the fields are calculated for each unit cells of the periodic structure.This is achieved by first calculating the scattering matrix for the total structure and then utilizing the transmission matrix method to calculate the fields within the slab.Let the field (or voltage) incident to the structure be V + (ω).The reflected voltage is V − (ω) = V + (ω)S 11 (ω), and the voltage and current at the material interface are calculated as V 1 (ω) = V + (ω) +V − (ω) and I 1 (ω) = (V + (ω) −V − (ω))/Z 0 , respectively.Then the voltage and current in the next unit cell are given by and so forth.(Since det T uc = 1, T uc is invertible.)By standard inverse Laplace transforms, we finally find the time-domain fields at each unit cell of the structure.Figure 9 shows the difference between the scattering parameter S 11 for the periodic structure and the reflection coefficient R for the effective medium slab.The length scale Λ is varied while Y /Λ and d = ΛN are fixed.To enable fair comparison, the reference plane of the periodic structure is shifted by an amount Λ/2 to the right, i.e., the periodic structure is symmetrized with respect to the point z = d/2.We see that there is a good agreement between the two calculations for Λ = λ 0 /30 and even better when Λ = λ 0 /50, as predicted by [24].Figure 10 shows the timedomain field when exciting the periodic structure with the frequency f 1 = 1.22 GHz.At this frequency the effective medium has ε = 0.063 − 0.23i.The excitation, cos(2π f 1 t), is initiated at t = 0 and shut off at t = 20 ns.We observe that the forerunner propagates straight through the material at the vacuum light velocity.The backward wave builds up before the reflection from the far end has returned to the region near z = 0 + .After 20 ns the excitation is turned off.Then the associated transient propagates at c, and the backward wave dies out, demonstrating stability.This behavior is consistent with the example given in Sec. 2. We stress that we have analyzed a discrete lumped model with finite d; the scattering parameters are calculated by the transmission matrix method so that the refractive index is not considered.
The periodic circuit model is based on ideal passive components along with an ideal, constant negative resistor.For a more realistic model of the components it is possible to include resistive, capacitive and inductive parasitics into the components in use.Nevertheless, preliminary sensitivity studies show that the circuit is not particularly sensitive to perturbations.The implementation of the negative resistor element should yield a broad-banded response compared to the two resonances.

Conclusion
We have presented theory and simulations which demonstrate how to determine the refractive index in active metamaterials.We have also proposed and analyzed a lumped circuit transmission line model for an active right-handed metamaterial with negative refractive index.This model demonstrates a backward propagating wave with phase velocity and Poynting's vector towards the source, and can be used as a starting point for implementation.