arxiv:physics/ v2 [physics.optics] 30 Nov 2005

Transcription

1 Orthogonal Phase Velocity Propagation of Electromagnetic Plane Waves arxiv:physics/5v [physics.optics 3 Nov 5 Tom G. Mackay School of Mathematics, University of Edinburgh, Edinburgh EH9 3JZ, UK Akhlesh Lakhtakia CATMAS Computational & Theoretical Materials Sciences Group Department of Engineering Science and Mechanics Pennsylvania State University, University Park, PA 68 68, USA Abstract In an isotropic, homogeneous, nondissipative, dielectric magnetic medium that is simply moving with respect to an inertial reference frame, planewave solutions of the Maxwell curl postulates can be such that the phase velocity and the time averaged Poynting vector are mutually orthogonal. Orthogonal phase velocity propagation thus adds to the conventional positive phase velocity propagation and the recently discovered negative phase velocity propagation that is associated with the phenomenon of negative refraction. Keywords: Minkowski constitutive relations, Poynting vector, Special theory of relativity Introduction Until quite recently, it was commonly held that the time averaged Poynting vector and the phase velocity of a plane wave in any homogeneous and isotropic medium are coparallel [, ; furthermore, the angle between the two vectors was commonly thought to be acute in homogeneous anisotropic mediums [3. The successful fabrication of what are called negative phase velocity ( materials in contrast to their conventional counterparts, the positive phase velocity ( materials led to a radical re think: the angle between the two vectors can be obtuse and even 8. Much has been written on materials during the past five years, and we refer the interested reader to three pedagogical reviews [4, 5, 6 and two Focus Issues of electronic journals [7, 8. Could there be materials that are neither of the nor of the types? We have already provided an affirmative answer to this question: there can be isotropic dielectric magnetic materials wherein the phase velocity is infinite at some frequency [9. In addition Corresponding Author. Fax: ; e mail: T.Mackay@ed.ac.uk. Some researchers prefer the infelicitous term left handed materials.

2 to the infinite phase velocity materials, we provide here a second affirmative answer: materials in which the phase velocity is orthogonal to the time averaged Poynting vector. Such materials may be called orthogonal phase velocity (OPV materials. The OPV materials emerge here as a gift of the special theory of relativity [,,. Relative to a co moving observer, these materials are supposed to be homogeneous, isotropic, and with dielectric and/or magnetic properties. Relative to an inertial observer that is not co moving, planewave propagation may occur such that the OPV condition is satisfied. The OPV condition, however, will never be satisfied for the co moving observer. The plan of this paper is as follows: Section is devoted to the Minkowski constitutive relations, planewave propagation, and the OPV condition. Section 3 contains numerical results for an OPV material that is non dissipative. Conclusions are presented in Section 4. In the notation adopted, Re{Q} and Im{Q} represent the real and imaginary parts, respectively, of a complex valued Q. The complex conjugate is written as Q. Vectors are identified by bold typeface and 3 3 dyadics are double underlined; ˆv is a unit vector co directional with v; the unit dyadic is I; and r denotes the spatial coordinate vector. The permittivity and permeability of free space (i.e., vacuum are ǫ and µ, respectively; c = (ǫ µ / is the speed of light in free space; ω is the angular frequency; and k = ω/c. Analysis. Minkowski constitutive relations As in a predecessor paper [3, two inertial reference frames, Σ and Σ, provide the setting for our analysis. The inertial reference frame Σ moves at constant velocity v = vˆv relative to the inertial reference frame Σ. The electromagnetic field phasors in frames Σ and Σ are represented as {D,E,H,B} and {D,E,H,B }, respectively. The two sets of phasors are related via the Lorentz transformation as [ E = (E ˆvˆv+ B = (B ˆvˆv+ H = (H ˆvˆv+ D = (D ˆvˆv+ β β β β [( I ˆvˆv [ (I ˆvˆv [( I ˆvˆv [ (I ˆvˆv E+v B, ( v E B, ( c H v D, (3 v H D+, (4 c where β = v/c. We consider a homogeneous medium which is characterized by the frequency domain dielectric magnetic constitutive relations } D = ǫ ǫ r E, (5 B = µ µ r B

3 intheco movingreference frameσ, withǫ r astherelative permittivityandµ r astherelative permeability. By virtue of ( (4, the constitutive relations (5 transform to the Minkowski constitutive relations mˆv H D = ǫ ǫ r α E+ c B = mˆv E +µ µ r α c H in the non co moving reference frame Σ; here,, (6 α = αi +( αˆvˆv, (7 α = β ǫ r µ r β, (8 m = β ǫ rµ r ǫ r µ r β. (9 The constitutive relations (6 degenerate to (5 as v. In the remainder of this paper, the propagation of plane waves in a medium described by the Minkowski constitutive relations (6 is investigated, with particular emphasis on the unusual behaviour which occurs when ǫ r µ r <. For further general background details on planewave propagation within the context of special relativity, the reader is referred to standard works [ [.. Planewave propagation Our focus lies on plane waves with field phasors E = E exp(ik r H = H exp(ik r }. ( The plane waves are taken to be uniform; i.e., the wavevector k = k ˆk where the unit vector ˆk R 3 and k = k R +ik I, ( where k R = Re{k} and k I = Im{k}. As a result of combining the frequency domain Maxwell curl postulates with the Minkowski constitutive relations and substituting ( therein, we find that [3 βξ k = k ˆk ˆv± ( ξ β ˆk, ( ˆv 3

4 wherein ξ = ǫ rµ r β, (3 = +(ǫ r µ r δ, (4 ( β ˆk ˆv δ =. (5 β Since the two wavenumbers represented by ( are not independent, the medium with constitutive relations (6 is unirefringent. The choice of sign for the square root term in ( is dictated by the direction of planewave propagation. The amplitudes of the electromagnetic field phasors turn out to be [3 where the orthogonal eigenvectors e and e are given by E = C e +C e, (6 C H = e ωǫ ǫ r C e, ωµ µ r (7 e = k ˆv, (8 e = α [(k+m ωc ˆv e, (9 whereas C and C are arbitrary constants. Exploiting the relation e = a e, ( where we find the time averaged Poynting vector e exp ( k Iˆk r P =.3 Nondissipative mediums a = k+ ξ(ω k v v, ( c { } C Re a + C ωµ µ ωǫ ǫ r a. ( r For nondissipative mediums (i.e., ǫ r R and µ r R, the following three scenarios may be identified: (i Ifǫ r µ r then > andtheexpression(representstworeal-valuedwavenumbers. (ii If < ǫ r µ r < then can be either positive or negative valued depending upon the magnitude of β and the relative orientations of ˆk and ˆv. Hence, either two real valued wavenumbers or a conjugate pair of complex valued wavenumbers are represented by (. In the instance where v =, we have = ǫ r µ r > and the wavenumbers are accordingly real valued. 4

5 (iii If ǫ r µ r < then < and ( provides two complex valued wavenumbers. In the case where v =, we have k R = which does not correspond to a propagating wave. We exclude this scenario from further consideration..4 Positive, negative and orthogonal phase velocity Equation ( for the time averaged Poynting vector simplifies to e exp ( k Iˆk r ( C ξ (ω k Rˆk v P = + C ωǫ ǫ r k Rˆk+ v, (3 ωµ µ r when ǫ r R, µ r R, and k C. The phase velocity v p = c k k Rˆk, (4 is described as being positive, negative or orthogonal depending upon its orientation relative to the time averaged Poynting vector. That is, positive phase velocity( is characterized by v p P >, negative phase velocity ( is characterized by vp P <, and orthogonal phase velocity (OPV is characterized by v p P =. In light of (3, the phase velocity for nondissipative mediums is positive, negative or orthogonal depending upon the sign of the scalar ( W = k R [k R +ξβ k k R β ˆk ˆv ˆk ˆv. (5 We distinguish between the scenarios and < for nondissipative mediums, as follows. (a. Here the wavenumber is real valued and with θ being the angle between ˆk and ˆv. W = k[k +ξβ (k kβcosθcosθ, (6 (b <. Here the wavenumber is complex valued, with real and imaginary parts given by βξ cosθ k R = k ξ(β cosθ ±, (7 k I = k ξ(β cosθ and W is null valued. c 5

6 3 Numerical illustrations Let us now illustrate the theory presented in the preceding section by means of some numerical results. In Figure, the real and imaginary wavenumbers are plotted as functions of the real valued product ǫ r µ r for β {.3,.6,.9} with θ {π/6,π/3,3π/4}, wherein θ denotes the angle between ˆk and ˆv. The wavenumber generally has a nonzero imaginary part, despite ǫ and µ being real valued. Furthermore, the positive valued upper bound on the range of ǫ r µ r values for which k I increases as the relative speed β increases. The magnitudes of the phase velocities, corresponding to the wavenumbers presented in Figure, are plotted in Figure. By definition, the phase velocity is aligned with the direction of ˆk. The positive, negative and orthogonal nature of the phase velocity, as deduced from the sign of the scalar W in (5, are mapped in Figure as functions of β and θ for ǫ r µ r {.,.4,.6,.8,.,.,.5,,3,4,5,}. The OPV region occupies a substantial portion of the βθ plane for small values of ǫ r µ r. As ǫ r µ r increases towards unity, the OPV region decreases and vanishes in the limit ǫ r µ r. The OPV region is symmetric about the line θ = π/, whereas the region is restricted to those regions of the βθ plane for which π/ < θ π [3 which shows that arises only when the velocity vector v casts a negative projection onto the wavevector k. When ǫ r µ r =, the phase velocity is positive for all values of β and θ. For ǫ r µ r >, OPV does not arise: the phase velocity is either positive or negative, with the region expanding as the product ǫ r µ r increases. 4 Concluding remarks Consider a homogeneous medium which is an isotropic dielectric magnetic medium in a co moving reference frame Σ, characterized by the relative permittivity scalar ǫ r > and the relative permeability scalar µ r >. From this perspective, the medium supports planewave propagation with real valued wavenumbers, and the phase velocity is directed parallel to the time averaged Poynting vector. The nature of planewave propagation is different when viewed in an inertial reference frame Σ, which moves uniformly at velocity v relative to Σ. The following two regimes are then distinguished: (a If ǫ r µ r then all wavenumbers are real valued, regardless of v. The projection of the phase velocity onto the time averaged Poynting vector can be either positive or negative, depending upon the magnitude and direction of v. (b If < ǫ r µ r < then planewave propagation characterized by both k R and k C (with k I is supported. Whether k I or not depends upon the magnitude of v and its orientation relative to k. For plane waves with k I, the time averaged Poynting vector is directed orthogonal to the phase velocity. 6

7 The regime < ǫ r µ r < lies between either vacuum (ǫ r = µ r = or anti vacuum (ǫ r = µ r = and nihility (ǫ r = µ r = [4, 5, 6. Anti vacuum is the vaunted material for the so called perfect lenses [7. Many schemes for materials satisfying the condition < ǫ r µ r < have been formulated [8, 9,, although we must note that dissipation then is hard to avoid. However, our calculations (not reported here show that the OPV condition can be very nearly satisfied by weakly dissipative materials. References [ Durney CH and Johnson CC 969 Introduction to Modern Electromagnetics (New York: McGraw Hill Chap 6 [ Lakhtakia A 994 Beltrami Fields in Chiral Media (Singapore: World Scientific Chap 7 [3 Born M and Wolf E 98 Principle of Optics, 6th ed (Oxford: Pergamon Chap 4 [4 Lakhtakia A, McCall M W and Weiglhofer W S 3 Negative phase velocity mediums Introduction to Complex Mediums for Optics and Electromagnetics ed WS Weiglhofer and A Lakhtakia (Bellingham, WA, USA: SPIE Press pp [5 Boardman AD, King N and Velasco L 5 Negative refraction in perspective Electromagnetics [6 Ramakrishna SA 5 Physics of negative refractive index materials Rep. Prog. Phys [7 Pendry JB (ed 3 Focus Issue: Negative Refraction and Metamaterials Opt. Exp. id=86 [8 Lakhtakia A and McCall M (eds 5 Focus on Negative Refraction New J. Phys. 7 [9 Lakhtakia A and Mackay TG 4 Infinite phase velocity as the boundary between positive and negative phase velocities Microw. Opt. Technol. Lett [ Chen HC 983 Theory of Electromagnetic Waves (New York: McGraw Hill Chap 8 [ Pappas CH 988 Theory of Electromagnetic Wave Propagation (New York: Dover Press Chap 7 [ Chawla BR and Unz H 97 Electromagnetic Waves in Moving Magneto plasmas (Lawrence, KS, USA: University Press of Kansas [3 Mackay TG and Lakhtakia A 4 Negative phase velocity in a uniformly moving, homogeneous, isotropic, dielectric magnetic medium J. Phys. A: Math. Gen

8 [4 Lakhtakia A An electromagnetic trinity from negative permittivity and negative permeability Int. J. Infrared Millim. Waves ; ibid [5 Ziolkowski RW 4 Propagation in and scattering from a matched metamaterial having a zero index of refraction Phys. Rev. E [6 Lakhtakia A and Mackay TG 6 Fresnel coefficients for a permittivity permeability phase space encompassing vacuum, anti vacuum, and nihility Microw. Opt. Technol. Lett. 48 at press [7 Pendry JB 4 Negative refraction Contemp. Phys [8 Rachford FJ, Smith DL, Loschialpo PF and Forester DW Calculations and measurements of wire and/or split ring negative index media Phys. Rev. E [9 Simovski CR and He S 3 Frequency range and explicit expressions for negative permittivity and permeability for an isotropic medium formed by a lattice of perfectly conducting Ω particles Phys. Lett. A [ Paniou NC and Osgood RM 3 Numerical investigation of negative refractive index metamaterials at infrared and optical frequencies Opt. Commun

9 5 4 3 = =.3 k R k I k R 5-5 =.6 k I = =.9.8 = k R k I Figure : The real (left and imaginary (right parts of the wavenumbers k, normalized with respect to k, plotted against ǫ r µ r for β {.3,.6,.9}. The values of θ are: 3 (solid curves, 6 (dashed curves and 35 (broken dashed curves. 9

10 8 6 =.3 4 v p =.6 v p =.9 v p Figure : The phase velocity, normalized with respect to c, plotted against ǫ r µ r for β {.3,.6,.9}. The values of θ are: 3 (solid curves, 6 (dashed curves and 35 (broken dashed curves. The phase velocity is directed along ˆk.

11 ΕrΜr=. OPV.6.4 ΕrΜr=.4 OPV.6.4 ΕrΜr=.6 OPV ΕrΜr=.8 OPV ΕrΜr= ΕrΜr= Figure 3: The distribution of positive phase velocity (, negative phase velocity ( and orthogonal phase velocity (OPV, as functions of β (, and θ (,8 for ǫ r µ r =.,.4,.6,.8,,.,.5,, 3, 4, 5,. The labels and shown hold for ǫ r > and µ r >, but must be reversed for ǫ r < and µ r < ΕrΜr= ΕrΜr=..6.4 ΕrΜr= ΕrΜr= ΕrΜr= ΕrΜr= Figure 3: continued