Quick viewing(Text Mode)

Differential Form Approach to the Analysis of Electromagnetic Cloaking and Masking F

Differential Form Approach to the Analysis of Electromagnetic Cloaking and Masking F

MOTL DRAFT 1 Approach to the Analysis of Electromagnetic Cloaking and Masking F. L. Teixeira

Abstract We bring attention to the relationship between (1) electromagnetic masking or cloaking of objects produced by some blueprints and (2) some invariances of obviated by means of the differential forms (exterior ) framework.

I. INTRODUCTION It has been recently shown that it is in principle possible to achieve (under certain idealized conditions) electro- magnetic cloaking of objects by using metamaterial coatings with properly designed constitutive parameters [1]. The basic feature underlying this remarkable possibility is the invariance of Maxwell equations under diffeomorphisms of the metric (metric invariance) [1], [2], [3], [4], [5], i.e., the fact that a change on the metric of space can be mimicked by a proper change of the constitutive . In this work, we discuss how this feature of Maxwell equations is obviated using differential forms and the exterior calculus framework [3], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14]. We then illustrate how this feature also allows for generic masking of objects (again, under idealized conditions) via appropriate metamaterial coatings.

II. MAXWELLEQUATIONSANDDIFFERENTIALFORMS In terms of differential forms, Maxwell equations in four-dimensional write in a very concise manner as [6] dF = 0 d F = J (1) ∗ where d is the four-dimensional , F is the electromagnetic strength 2-form, J is the source density 3-form, and is the four-dimensional . By performing a 3 + 1 splitting and treating time as a ∗ parameter, we can write F in terms of the electric field intensity 1-form E and the magnetic flux density 2-form B, and write J in terms of the (electric) 2-form J and the (electric) density 3-form ρ as F = E dt + B, ∧ J = ρ J dt. (2) − ∧ In the above, d is the three-dimensional (spatial) exterior derivative, and is the exterior product. The exterior ∧ derivatives d and d are both metric-free and nilpotent, and related through ∂ d = d + dt . ∧ ∂t It is interesting to note that the nilpotency of d implies dJ = 0, a statement of charge conservation. On each spatial slice, three-dimensional (spatial) Hodge star operators ⋆ǫ and ⋆µ can be defined such that [9], [10]

B = ⋆µH,

D = ⋆ǫE, (3)

Manuscript received MM DD, YYYY; revised MM DD, YYYY. This work was supported in part by NSF under grant ECS-0347502 and by AFOSR under grant FA-9550-04-1-0359. F. L. Teixeira is with ElectroScience Laboratory and Department of Electrical and Computer , The Ohio State University, Columbus, OH 43210, USA. [email protected] MOTL DRAFT 2 where D is the electric flux density (twisted) 2-form, and H is the magnetic field intensity (twisted) 1-form. In terms of ⋆ǫ and ⋆µ, the four-dimensional Hodge can be written as ∗ −1 G = F = ⋆ B dt + ⋆ǫE = H dt + D. ∗ − µ ∧ − ∧ From dF = 0, we obtain ∂ (d + dt )(B + E dt) = 0, ∧ ∂t ∧ and hence ∂B dB + (dE + ) dt = 0. (4) ∂t ∧ Moreover, from dG = d F = J, we have ∗ ∂ (d + dt )(D H dt)= ρ J dt, ∧ ∂t − ∧ − ∧ and hence ∂D (dD ρ) + (J dH + ) dt = 0 (5) − − ∂t ∧ From (4) and (5), Maxwell equations in 3 + 1 can be decomposed as [5] ∂B dE = , − ∂t dB = 0, ∂D dH = + J, ∂t dD = ρ, (6) where the only spatial operator present is d. When applied to 0-, 1-, or 2-forms, respectively, d is equivalent to the more familiar grad, curl, and div operators of vector calculus, distilled from their metric structure. Indeed, since d is metric-free, none of the equations in (6) depend on the metric of space (i.e., they are under diffeomorphisms) [5]. All metric information is incorporated into the spatial-slice constitutive equations (3) that map 1-forms into 2-forms in 3-space [more generally, Hodge operators map p-forms into (n p)-forms in n-space] together with the material properties. −

III. DIFFEOMORPHISMSONTHEMETRICOFSPACE Under a differentiable but otherwise generic transformation (diffeomorphism) on the metric of space, g g˜, 7→ the spatial-slice constitutive equations transform to

B = ⋆˜µH,

D = ⋆˜ǫE. (7) In other words, by using differential forms, a change on the metric is manifestly equivalent to a change on the material properties of the medium, now encoded by ⋆˜ǫ and ⋆˜µ. This is to say that, for a given (doubly anisotropic, for generality) medium with constitutive tensors [ǫ] and [µ], either one of the transformations ([ǫ], [µ], g) ([ǫ], [µ], g˜) 7→ ([ǫ], [µ], g) ([˜ǫ], [˜µ], g) (8) 7→ where [˜ǫ] and [˜µ] are new particular constitutive tensors, produce the same differential form solutions, as illustrated by the diagram in Fig. 1. In terms of vector fields, either of these two transformations produce the same field distributions outside the region affected by the diffeomorphism g g˜ and up to a scale factor inside that region [15]. More specifically, an arbitrary diffeomorphism on the metric g 7→g˜ can be described by a differentiable coordinate 7→ transformation (u1, u2, u3) (˜u1, u˜2, u˜3) with [S] as the Jacobian of partial derivatives, Sij = ∂ui/∂u˜j. 7→ MOTL DRAFT 3

For an original medium with constitutive tensors [ǫ] and [µ], this transformation on the metric produces identical differential form solutions as those resulting from the following transformation on the constitutive tensors [15], [16] −1 [˜ǫ] = (det[S]) [S] [ǫ] [S]T , −1 · · [˜µ] = (det[S]) [S] [µ] [S]T , (9) · · since the same Hodge operators ⋆˜ǫ and ⋆˜µ are recovered. The differentiability of the coordinate transformation implies continuity of [S]. This in turns implies that the intrinsic impedance is not changed, akin to a PML transformation [15].

Fig. 1. Decomposition of 3+1 Maxwell equations into a topological part (top four equations) and a metric part (boxed equations), obviated by using differential forms. Under diffeomorphisms, the metric g gets mapped to g˜. This can be translated to a modification in the Hodge operators ⋆ǫ and ⋆µ, which can be made equivalent to a change on the (spatial-slice) constitutive properties of the underlying medium: [ǫ] [˜ǫ] and [µ] [˜µ]. 7→ 7→

IV. HODGEOPERATORSANDMASKINGBLUEPRINTS Using the above properties, metamaterial blueprints that produce electromagnetic masking of objects can be derived in a systematic fashion. The term masking denotes here the use of a metamaterial blueprint for coating a physical scatterer S0 so that the resulting scattered field (from the coated object) becomes identical to that from a different (non-coated) scatterer S1 (mask). In the case of a zero volume mask, masking recovers cloaking [1]. This is illustrated in Fig. 2. One example of such masking where a closed-form expression for the final metamaterial blueprints has been derived in [17]. In that case, a convex PEC object S0 with surface radii of curvature in the range r [rm, rM ], ∈ with rM rm, is coated by a metamaterial blueprint of uniform thickness δ such that the resulting scattered field ≥ becomes identical to that of a (noncoated) PEC S1 in a isotropic medium with constitutive parameters ǫ and µ having radii of curvature in the range r [0, rM rm]. At any point p at the boundary of S0, one can define ∈ − orthonormal vectors uˆ1, uˆ2, nˆ, where uˆ1 and uˆ2 are tangent along the principal lines of curvature, and nˆ =u ˆ1 uˆ2 × is outward normal. The local tangential coordinates are denoted as ζ1, ζ2 and the normal coordinate is η, with η = 0 parametrizing the surface of S0 itself. If the principal radii of curvature at p are given by rs1 and rs2, then ′ at a point p along η for fixed (ζ1,ζ2), the radii of curvature of parallel surfaces are ri(ζ1,ζ2, η)= rsi(ζ1,ζ2)+ η, ′ i = 1, 2. In particular, the radii of curvature at a point p on the surface of the mask S1 are r0i(ζ1,ζ2), i = 1, 2, with r0i(ζ1,ζ2)= rsi(ζ1,ζ2) rm, and rm = minp[rs1(p), rs2(p)]. The required diffeomorphism on the metric for − this case is represented by the map : η [ rm, δ] η˜ [0, δ], η˜ = (η)= δ(rm + η)/(rm + δ) (for η > δ, M ∈ − 7→ ∈ M MOTL DRAFT 4

Fig. 2. By coating an arbitrary PEC scatterer S0 (at the left) with a properly designed masking metamaterial the scattered field becomes equal to that of a (mask) object S1 (at the right) in the region outside the dashed curve. If the mask is a zero volume object, cloaking is produced. we have η˜ = η) [17]. The metamaterial coating that results is a inhomogeneous and anisotropic material blueprint with constitutive tensors [˜ǫ]= ǫ[Λ], [˜µ]= µ[Λ] where

1 h˜1h2 1 h1h˜2 h1h2 [Λ](ζ1,ζ2, η) =u ˆ1uˆ1 +u ˆ2uˆ2 +n ˆnsˆ η , (10) sη h1h˜2 sη h˜1h2 h˜1h˜2 −1 ˜ −1 where hi = rsi (rsi +η) and hi = rsi (rsi +η ˜), i = 1, 2, and sη = ∂η/∂η˜ = δ/(rm +δ) [17]. The derivation of this metamaterial blueprint has close analogy with the derivation of Maxwellian perfectly matched layers (PMLs) for conformal boundaries [18]. This is because the latter can be derived in the various coordinate systems by applying a complex coordinate stretching (i.e., complexification on the metric of space) and the above metric invariance [10], [18]. As another example, conformal maps [19] are a class of diffeomorphic maps that preserve both orthogonality of coordinates and isotropy of the (transformed) constitutive parameters. In this case, g g˜ = Ω(u)g, such that 7→ Ω(u) is a strictly positive scalar function of position. In an orthogonal system (u1, u2, u3) it then follows that Sij = δij/√Ω (where δij is the ). Using (9), one can then write: [˜ǫ] = √Ω [ǫ], [˜µ] = √Ω [µ]. (11)

V. CONCLUSIONS We have used the differential forms (exterior calculus) framework to provide an unified analysis of the derivation of metamaterial blueprints for cloaking and masking of electromagnetic scatterers. Differential forms are quite suited for this purpose because they factorize Maxwell equations into a metric part and a (purely) tological part, and hence obviate the invariance of topological part under diffeomorphisms on the metric of space. The discussion has been restricted to the material blueprint level only [20]. Issues about the realizability and limitations on manufacturability of such have not been addressed here. This is touched upon elsewhere [1], [20]. The realization and experimental evaluation of approximate cloaks have been considered in [21].

REFERENCES [1] J. B. Pendry, D. Schurig, and D. R. Smith “Controlling electromagnetic fields,” , vol. 312, pp. 1780-1782, Jun. 2006. [2] D. Van Dantzig, “The fundamental equations of , independent of metrical geometry,” Proc. Cambridge Phil. Soc., vol. 30, pp. 421–427, 1934. [3] G. Deschamps, “Electromagnetics and differential forms,” Proc. IEEE, vol. 69, no. 6, pp. 676-696, June 1981. [4] A. J. Ward and J. B. Pendry, “Refraction and geometry in Maxwell’s equations,” J. Mod. Opt, vol. 43, pp. 773-793, 1996. [5] F. L. Teixeira and W. C. Chew, “Lattice electromagnetic theory from a topological viewpoint,” J. Math. Phys., vol. 40, no. 1, pp. 169–187, 1999. [6] K.S.Thorne, C. W. Misner, and J. A. Wheeler, Gravitation, Freeman, San Fransisco, 1973. MOTL DRAFT 5

[7] C. V. Westenholz, Differential Forms in Mathematical , Elsevier Science, North-Holland, 1980. [8] P. R. Kotiuga, “Helicity functionals and metric invariance in three dimensions,” IEEE Trans. Magn., vol. 25, no. 4, pp. 2813-2186, Jul. 1989. [9] K. F. Warnick and D. V. Arnold, “Green forms for anisotropic, inhomogeneous media,” J. Electromagn. Waves Applicat., vol. 11, pp. 1145-1164, 1997. [10] F. L. Teixeira and W. C. Chew, “Differential forms, metrics, and the reflectionless absorption of electromagnetic waves,” J. Electromagn. Waves Appl., vol. 13, no. 5, pp. 665–686, 1999. [11] A. Bossavit, “Generalized finite differences in computational electromagnetics,” in Geometric Methods for Computational Electromag- netics (F. L. Teixeira, ed.), PIER 32, EMW Publishing, Cambridge, Mass., pp. 45-64, 2001. [12] F. W. Hehl and Yu. N. Obukhov, Foundations of Classical Electrodynamics: Charge, Flux, and Metric, Birkhauser, Boston, 2003. [13] P. W. Gross and P. R. Kotiuga, Electromagnetic Theory and Computations: A Topological Apporach, Cambridge University Press, Cambridge, UK, 2004. [14] B. He and F. L. Teixeira, “On the degrees of freedom of lattice electrodynamics,” Phys. Lett. A 336, pp. 1-7, 2005. [15] F. L. Teixeira and W. C. Chew, “General closed-form PML constitutive tensors to match arbitrary bianisotropic and dispersive linear media,” IEEE and Guided Wave Letters, vol. 8, no. 6, pp. 223-225, 1998. [16] G. W. Milton, M. Briane, and J. R. Wills, “On cloaking for elasticity and physical equations with a transformation invariant form,” New J. Phys., vol. 8, 248, 2006. [17] F. L. Teixeira, “Closed-form metamaterial blueprints for electromagnetic masking of arbitrarily shaped convex PEC objects,” IEEE Antennas Propagat. Lett., vol. 6, 2007 (to appear). [18] F. L. Teixeira and W. C. Chew, “Analytical derivation of a conformal perfectly matched absorber for electromagnetic waves,” Microwave Opt. Technol. Lett., vol. 17, no. 4, pp. 231–236, 1998. [19] C. N. Yang , “Conformal mapping of gauge fields,” Phys. Rev. D, vol. 16, no. 2, pp. 330-331, 1977. [20] F. L. Teixeira, “On aspects of the physical realizability of perfectly matched absorbers for electromagnetic waves,” Radio Sci., vol. 38, no. 2, 8014, 2003. [21] 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 ,” Science, vol. 314, no. 5801, pp. 977-980. 2006.