Arxiv:1905.08341V1 [Physics.App-Ph] 20 May 2019 by Discussing the Inﬂuence of the Presented Method in Not Yet Conducted for Cylindrical Metasurfaces

Illusion Mechanisms with Cylindrical Metasurfaces: A General Synthesis Approach

Mahdi Safari1, Hamidreza Kazemi2, Ali Abdolali3, Mohammad Albooyeh2,∗, and Filippo Capolino2 1Department of Electrical and Computer Engineering, University of Toronto, Toronto, Canada 2Department of Electrical Engineering and Computer Science, University of California, Irvine, CA 92617, USA 3Department of Electrical Engineering, Iran University of Science and Technology, Narmak, Tehran, Iran corresponding author: ∗[email protected]

We explore the use of cylindrical metasurfaces in providing several illusion mechanisms including scattering cancellation and creating fictitious line sources. We present the general synthesis approach that leads to such phenomena by modeling the metasurface with effective polarizability tensors and by applying boundary conditions to connect the tangential components of the desired fields to the required surface polarization current densities that generate such fields. We then use these required surface polarizations to obtain the effective polarizabilities for the synthesis of the metasurface. We demonstrate the use of this general method for the synthesis of metasurfaces that lead to scattering cancellation and illusion effects, and discuss practical scenarios by using loaded dipole antennas to realize the discretized polarization current densities. This study is the first fundamental step that may lead to interesting electromagnetic applications, like stealth technology, antenna synthesis, wireless power transfer, sensors, cylindrical absorbers, etc.

I. INTRODUCTION formal metaurfaces with large radial curvatures (at the wavelength scale) [41]. However, that method was based on the analysis of planar structures with open bound- Metasurfaces are surface equivalents of bulk meta- aries while a cylindrical metasurface can be generally materials, usually realized as dense planar arrays of closed in its azimuthal plane and it involves rather differ- subwavelength-sized resonant particles [1–4]. Since the ent interaction mechanisms. Our formulation covers the emergence of this research topic, most studies have most general case of metasurfaces with dipolar (electric been carried out to analyze and synthesize infinitely ex- and magnetic) resonant responses including bianisotropic tended planar metasurfaces [5–31]. The synthesis of cases [42]. By applying our synthesis approach, we fur- some topologies other than planar, namely cylindrical ther conceptually synthesize several examples with inter- and spherical ones, has been studied in the past [32–40] esting practical applications, e.g., scattering cancellation but without a general systematic approach. Cylindrical from conductive and dielectric cylinders and also gen- topologies are among the commonly used structures in erating fictitious line sources, appearing away from the engineering electromagnetics and optics. For instance, original line source. cloaking, radar cross section reduction, obtaining an ar- The paper is structured as follows. We first present the bitrary radiation pattern from cylindrical objects, and theory for the analysis of infinitely extended (along the etc. are only a few problems which involve metasurfaces axial direction) cylindrical metasurfaces in Secs. II by synthesis in cylindrical coordinates. providing the relation between the desired fields at the For the analysis and synthesis of planar metasurfaces, metasurface boundary and the required surface polariza- several techniques have been reported as in Refs [9– tion densities. Next, we exhibit an expansion of the inci- 14, 16–23, 25–28] for example. These works are based on dent (here, plane waves and line sources) and scattered modeling the metasurfaces with impedance tensors [9– fields in cylindrical coordinates in Secs. III and IV. After 14, 16–20], with equivalent conductivities and reac- that, we conceptually synthesize two examples of scatter- tances [21], or with effective surface susceptibility or po- ing cancellation and an illusion device in Secs. V A, V B, larizability tensors [22, 23, 25–28]. Although the analysis and V C. We further present a practical realization ap- and synthesis of planar metasurfaces are well-established proach by discretizing the continuous required polariza- and discussed in the literature, to the best of our knowl- tion currents in Sec. VI. Finally, we conclude the study edge, such a comprehensive analysis and/or synthesis is

arXiv:1905.08341v1 [physics.app-ph] 20 May 2019 by discussing the influence of the presented method in not yet conducted for cylindrical metasurfaces. There- practical applications and new possibilities by cylindri- fore, we believe there is a need to develop an extensive cal metasurfaces in engineering electromagnetic waves. analysis and synthesis method for assigning prescribed physical properties to cylindrical metasurfaces. In particular here we aim at developing some illusion mechanisms using cylindrical metasurfaces. To that II. PROBLEM FORMULATION AND end, we present a comprehensive analysis and synthesis BOUNDARY CONDITIONS method for cylindrical metasurfaces analogous to what has been previously performed for the case of planar An electromagnetic metasurface is commonly under- metasurfaces in Refs. [26, 27]. Note that the analysis of stood as a composite layer composed of infinite number planar metasurfaces has been already extended to con- of subwavelength-sized inclusions which are densely ar- 2 ranged over a surface and is capable of manipulating the sors relate the polarization surface-density vectors P and wave front in a desired fashion (see e.g. Refs. [3]). This M to the incident electric and magnetic fields Ei and Hi definition includes also the case of non-planar conformal through the constitutive relations [42] surfaces. Here, we study a class of metasurfaces when ee em the inclusions are arranged to form a cylindrical surface P = α¯ · Ei + α¯ · Hi, (1) which is closed in the azimuthal plane and is infinitely extended in the axial direction. Note that in the limiting ¯me i ¯mm i case when the cylinder curvature is infinitely large, the M = α¯ · E + α¯ · H . (2) current problem transitions to the case of planar metasur- The superscript “i” denotes “incident” fields, and the ef- faces which is widely studied (see e.g. Refs. [16, 21–28]). fective surface polarizability tensors take into account of However, in sharp contrast with planar metasurfaces, in the effect of multiple reflections inside the surface, i.e., of a cylindrical configuration waves bounce inside the cylin- all the interactions among all the inclusions with them- drical metasurface and multiple reflections form the inner selves. The main goal in our synthesis problem is to find part of the cylindrical boundaries must be considered in the required effective polarizability tensors α¯ for a de- the analysis. sired wave manipulation. We note that in some other Let us consider a metasurface composed of resonant in- studies, the impedance or susceptibility tensors are the clusions that form a cylindrical surface with radius a [see main unknowns of the sysntehis problem [17–20, 25], ver- Fig. 1]. The cylindrical surface is assumed to be infinitely sus the effective polarizability tensors as in this current investigation. (a) z (b) z Similarly to what was done for planar metasurfaces [25– II, a IIem, em 27], we start from the boundary conditions for cylindrical metasurfaces which read (see Eqs. (4.3a)–(4.3d) in Ref. [43] and also Appendix A for a detailed derivation Ei 0 and the constitutive relations) Hi P

M + − 1 ∂Pρ x Ez − Ez = − + jωMφ, (3) x 0 ∂z   , ++, −− + − 1 ∂Pρ Eφ − Eφ = − − jωMz, (4) FIG. 1: (a) A cylindrical metasurface composed of a0 ∂φ nonidentical inclusions which is infinitly extended in the axial direction (here z axis) and is exposed to incidences from inside (i.e., infinite line sources with electric and + − 1 ∂Mρ Hz − Hz = − − jωPφ, (5) magnetic current amplitudes Ie and Im) and outside µ0 ∂z (i.e., a plane wave). (b) The electromagnetic modeling of the problem in (a). The metasurface is replaced by 1 ∂M two surface polarization densities P and M that provide + − ρ Hφ − Hφ = − + jωPz, (6) the same electromagnetic response as the problem in aµ0 ∂φ (a). Here the + and − signs in the superscripts correspond to the total fields outside and inside of the metasurface extended along the axial direction (here z-direction). boundary (i.e. ρ = a), respectively, and the implicit time Moreover, we assume that the inner and outer spaces dependence ejωt is assumed, with ω being the angular of the cylindrical surface are filled with two different frequency. Moreover, Eφ,z, Hφ,z, Pρ,φ,z, and Mρ,φ,z are isotropic materials with the permittivity ± and perme- correspondingly the field and surface polarization density ability µ±, respectively (see Fig. 1, the “+” sign refers to components in cylindrical coordinates. Furthermore, µ0 the material properties of the outside medium whereas and 0 are the free space permeability and permittivity, the “−” sign corresponds to the inside one). We illu- respectively. Indeed, Eqs. (3)–(6) relate the jump of minate the proposed metasurface either from outside or the tangential components of the electric and magnetic from inside. Due to the external/internal illumination, fields across the metasurface sheet to the required sur- the effective surface electric and magnetic polarization face electric and magnetic polarization densities. The densities P (with the unit of Asm−1) and M (with the right hand sides of Eqs. (3) and (4) represent the total unit of Vsm−1) are, respectively, induced on the meta- equivalent tangential magnetic surface polarization den- surface boundary (see e.g. [22–27]). Next, similar to sities contributing to the discontinuity of the tangential Refs. [22, 23, 25–27], we express the metasurface response electric field, whereas the right hand sides of Eqs. (5) in terms of effective electric, magnetic, magnetoelectric, and (6) represent the total equivalent tangential electric and electromagnetic polarizability tensors α¯ee, α¯mm,α¯em, surface polarization densities contributing to the discon- and α¯me, respectively. These surface polarizability ten- tinuity of the tangential magnetic field [44–46]. 3

∂ Next, for simplicity we only consider cases with vanish- assuming no z variation, and considering ∂φ → jn, the ing normal polarization components Pρ, and Mρ in the boundary conditions for each Fourier coefficient due to boundary conditions (3)–(6) to simplify the analysis and the orthogonality of different n-indexed harmonics of the synthesis without loosing any freedom in the manipula- total field at ρ = a read tion of a desired electromagnetic field [47]. Indeed, based on the Huygens principle, or more precisely from the + − equivalence theorem [59], it can be shown that a desired Ez,n − Ez,n = +jωMφ,n, (15) field in a given volume can be fully engineered know- ing the tangential (with respect to the metasurface sur- + − face) surface polarization components on the metasurface Eφ,n − Eφ,n = −jωMz,n, (16) sheet [45, 46]. As a result, the boundary conditions for the total field (3)–(6) simplify to + − Hz,n − Hz,n = −jωPφ,n, (17) + − Ez − Ez = +jωMφ, (7) + − Hφ,n − Hφ,n = +jωPz,n. (18) + − Eφ − Eφ = −jωMz, (8) Equations (15)–(18) are the final forms of the boundary conditions with the assumed simplifications and will be + − Hz − Hz = −jωPφ, (9) exploited hereafter.

+ − H − H = +jωPz. (10) φ φ III. INCIDENT FIELDS In the next step, since cylindrical structures are inher- ently periodic with respect to the azimuthal coordinate φ, We consider two different scenarios: the illumination we express the general form of the electromagnetic fields source either outside or inside the metasurface boundary (inside and outside of the closed metasurface boundary) [see Fig. 1]. Based on our examples in Sec. V, we consider and the surface polarization densities on the metasurface the plane wave excitation for the case of outside illumi- as Fourier series nation whereas we discuss the infinite line source for the n=∞ case of illumination from inside. However, we note that X jnφ E(ρ, φ, z) = En(ρ, z)e , (11) it is possible to consider any other excitation types based n=−∞ on the same analysis and procedure that follows.

n=∞ X jnφ A. Plane wave illumination from outside H(ρ, φ, z) = Hn(ρ, z)e , (12) n=−∞ For generality, we consider the superposition of two

n=∞ plane waves as external illumination: one with a z- X jnφ polarized electric ﬁeld i.e., P(φ, z) = Pn(z)e , (13)

n=−∞ i −jβ+x ETM = zˆE0e , (19)

n=∞ and the other with a z-polarized magnetic field i.e., X jnφ M(φ, z) = Mn(z)e , (14) i −jβ+x n=−∞ HTE = zˆH0e , (20) respectively, where En(ρ, z) and Hn(ρ, z) are the coef- which are traveling along the +x direction (normal to ficients of the electromagnetic fields in the Fourier se- the metasurface axis z) and impinging on the cylin- z ries whereas Pn(z) and Mn(z) are the coefficients corre- drical metasurface. While the former is called TM – sponding to the electric and magnetic surface polariza- transverse magnetic– polarization, the latter is called z tion densities on the metasurface sheet. TE –transverse electric [see Fig. 1]. Here, β+ is the prop- Before the next step and with the goal of practical imple- agation constant in the medium outside of the meta- mentation of the problem, we add one more level of sim- surface boundary, whereas E0 and H0 are the incident plification to our problem. That is, we consider no varia- electric and magnetic field amplitudes, respectively. Any tion of the fields and polarization densities along the axial plane wave which is propagating along the x direction direction (here, z-direction) as we provide two examples is decomposable into these two modes. For example, in this study. It should be noted that the general synthe- let us consider a plane wave hitting the metasurface as sis of three-dimensional problems is a simple generaliza- shown in Fig. 1, with its electric field component mak- tion of the z-invariant case studied here. Therefore, ap- ing an angle θ with the z-axis. In this case, we have 0 p plying (11)–(14) into the boundary conditions (7)–(10), η+H0 = E0 tan θ0, where η+ = µ+/+ is the intrinsic 4

(2) wave impedance of the corresponding medium. where H0 is the 0-th order Hankel function of the second Next, to solve the problem of scattering from cylindri- kind. Therefore the total incident electric field is cal structures, it is convenient to express the fields of i 1 i i Eqs. (19) and (20) in terms of the cylindrical wave func- E = ∇ × HTE + ETM tions. It can be shown that the TMz incident plane wave jω− with the electric field (19) can be expressed as [48] β h 0 i = − − φˆ jI H(2) (β ρ) + zˆ η I H(2)(β ρ) . 4 m 0 − − e 0 − n=∞ (26) i X −n jnφ Note that the sources could be located also elsewhere ETM = zˆ E0 j Jn(β+ρ)e , (21) 0 n=−∞ inside the cylindrical metasurface boundary e.g. at ρ away from the origin, and in this case one would need 0 and that with TMz polarization has magnetic field (20) to apply the substitution ρ = |ρ − ρ |, where ρ is the represented as observation point. In the next steps, we present the general forms of the n=∞ total fields due to the induced sources on the metasurface i X −n jnφ boundary, i.e., due to the secondary sources. HTE = zˆ H0 j Jn(β+ρ)e . (22) n=−∞

IV. FIELDS IN CYLINDRICAL COORDINATES Here Jn is the n-th order Bessel function and ρ is the radial position. As a result, the total incident electric field which is the superposition of both the TE and TM Generally, the fields inside (Ein) and outside (Eout) of incident plane waves reads the cylindrical metasurface boundary generated by the induced currents (scattered fields) can be expressed as

i 1 i i n=∞ E = ∇ × HTE + ETM X Jn(β−ρ) jω+ Ein = ρˆ nbn β−ρ n=∞ n=−∞ (27) nJ (β ρ) X −n n + ˆ 0 i = j η+H0 ρˆ + φ jJn(β+ρ) ˆ 0 jnφ β ρ +φ jbnJn(β−ρ) + zˆ anJn(β−ρ) e , n=−∞ + jnφ +zˆ E0Jn(β+ρ)] e , and (23) n=∞ " (2) ˆ X Hn (β+ρ) where ρˆ, φ, and zˆ are the unit vectors in cylindrical coor- E = ρˆ nd 0 out n β ρ dinates and Jn is the derivative of the n-th order Bessel n=−∞ + (28) function with respect to the argument β+ρ. ˆ (2)0 (2) i jnφ +φ jdnHn (β+ρ) + zˆ cnHn (β+ρ) e ,

respectively. Note that the scattered field expressions in B. Line source illumination from inside Eqs. (27) and (28) do not include the fields produced by either the inner line source or the incident plane wave, and only consider the fields generated by the induced Here we consider the superposition of two infinitely ex- currents on the cylindrical metasurface, i.e., secondary tended electric and magnetic line sources located at the (2) (2)0 center of the cylindrical metasurface as excitation source. sources. Here Hn and Hn are the n-th order Hankel The current amplitude of the electric and magnetic line function of the second kind and its derivative with re- sources are assumed to be Ie and Im, respectively. Simi- spect to the argument, respectively. Moreover, the first larly to the previous case of plane wave illumination from two components i.e., ρ and φ components in Eqs. (28) outside, here the electric line source creates a TMz field and (27) correspond to the TEz polarized fields while while the magnetic one generates a TEz field. The elec- the z component corresponds to the TMz polarized fields. tric and magnetic fields of the two kinds of sources are, Furthermore, due to the nonsingular nature of the scat- respectively, given by [48] tered fields inside and outside we have expanded the fields in terms of the suitable Bessel and Hankel functions, re- 2 spectively. The magnetic fields (Hout and Hin) can be i β− (2) ETM = −zÎe H0 (β−ρ), (24) calculated using Maxwell’s equations. Note that an and 4ω− bn are coefficients of Bessel-Fourier series representing standing waves inside the metasurface boundary for TMz z and and TE waves, respectively, whereas cn and dn are the coefficients which represent the TMz and TEz waves out- 2 i β− (2) side the metasurface, correspondingly. HTE = −zÎm H0 (β−ρ), (25) 4ωµ− At this stage we have all tools to synthesize a specific and 5 desired electromagnetic wave. To summarize, in order to illusory source since an observer outside the metasurface synthesize a desired wave profile: boundary sees only a line source which is displaced with respect to the original source. The realization of 1. We first expand it in terms of the cylindrical har- an illusion device using the available approaches such monics as given by Eqs. (27) and (28) to obtain the as TO requires the source to be located in a bulky unknown coefficients an, bn, cn, and dn. complex medium which makes it mainly impractical. 2. Next, by applying the given incident fields (23) However, in our approach, the real source is simply or (26) and by exploiting the boundary condi- located in free space and we surround it with a thin metasurface. Therefore, the advantages of our approach tions (15)–(18), we find the required coefficients Pn compared to the previous approaches for the design of and Mn that provide the required polarization densities P and M. electromagnetic devices for cylindrical structures is two- fold; the proposed approach is both practical and general. 3. Finally, when the polarization densities are found, the required effective polarizability tensors α¯ee, α¯mm,α¯em, and α¯me are retrieved from the constitutive relations (1) and (2). A. Scattering cancellation from a PEC cylinder As a final remark, note that each polarization density In this first example we consider a PEC cylinder with vector has two tangential components Pϕ,Pz and Mϕ,Mz providing a total of four functional equations. We recall radius a = 5λ (λ is the wavelength of the excitation field) that for simplicity we consider metasurfaces that express located in free space which is extended infinitely along its only tangential polarization densities and exclude those axis (here the z-axis) and is illuminated by a z-polarized plane wave propagating along the x-axis as given by metasurfaces with normal polarization densities Pρ and z −1 Eq. (19) (i.e., TM incidence) with E0 = 1 Vm and Mρ. However, each polarizability tensor in Eqs. (1) and ω/(2π) = 1 GHz. Our goal is to find a proper meta- (2) has four tangential polarizability components αϕϕ, surface which covers the PEC cylinder and suppresses αϕz, αzϕ, αzz. Therefore, the solution to the synthesis of a metasurface is not unique since there are two polariza- the scattered electromagnetic field by the PEC. In order tion vectors and four polarizability tensors, i.e., there are to cancel the scattered field, the total electric and magnetic fields must satisfy E+ = Ei, H+ = Hi outside the four equations with 16 complex unknown polarizability − − components. metasurface, and inside the PEC, E = 0 and H = 0, respectively. Therefore, by applying the boundary conditions (15)–(18) one finds the required surface polariza- V. ILLUSTRATIVE EXAMPLES tion density coefficients Pn and Mn at the boundary. Hence, the desired surface electric and magnetic polarization densities P and M at the boundary are obtained We present three representative examples to demon- from Eqs. (13) and (14) as strate the application of the proposed method and to clarify the synthesis approach. Other manipulations of n=+∞ E0 X electromagnetic waves in cylindrical coordinates shall be P = −zˆ J 0 (β a)j−nejnφ, (29) ωη n 0 possible using a similar approach. 0 n=−∞ In the first example we consider the application of a metasurface as a scattering cancellation device [34, 49, 50] and for a cylindrical perfect electric conductor (PEC) when n=+∞ E0 X the system is excited by a plane wave. In the second M = φˆ J (β a)j−nejnφ. (30) jω n 0 example we investigate an application similar to that of n=−∞ the first example with the only difference that the PEC cylinder is replaced by a dielectric cylinder. The main Here, η = pµ / is the intrinsic impedance of free 0 0 √0 advantage in using a metasurface for the realization of space and β0 = ω µ00 is the propagation constant of a scattering cancellation device compared to previous the incident wave (i.e., β0a = 10π in this example). Note z methods such as transformation optics (TO) [50, 51] that for the proposed TM incidence we have θ0 = 0 and transmission line method [52] is that one does not [see Fig. 1]. Indeed for this special case, the required φ- need to use bulky materials with exotic properties which component of P and the z-component of M vanish, i.e., inevitably take a large space and/or contain high losses. Pz and Mφ are the sufficient surface polarization density Instead, a simple electromagnetically thin composite components at the boundary to achieve the scattering surface which is easier to fabricate and take less space is cancellation. Similarly, for TEz incidence as in Eq. (20), exploited. the surface polarization density components Pφ and Mz In the last example we consider a metasurface that are sufficient to suppress the scattered fields. Note that surrounds a line source that creates a fictitious line based on the reciprocity theorem the scattered fields of source outside the cylindrical metasurface boundary. any impressed surface electric polarization density on the Following Ref. 51, we call this fictitious source the surface of a PEC is zero, therefore, the obtained surface 6 electric polarization density P in Eq. (29) represents the marize, (i) P is the induced surface polarization den- induced current on the PEC surface due to the incident sity on the PEC surface and not a surface polarization plane wave. However, the magnetic polarization density density impressed on the metasurface; and (ii) the sur- M at the boundary is an impressed polarization den- face polarization density Mϕ given by Eq. (30) to pro- sity provided by the metasurface to cancel the scattered vide scattering cancellation is obtained with a realiza- fields by the PEC. The surface polarization densities of tion of the metasurface with effective surface polarizabil- mm Eqs. (29) and (30) that synthesize the scattering cancel- ity component αϕϕ shown in Fig. 2(c). As it is clear, mm lation of a PEC cylinder are depicted in Fig. 2(a) and the effective polarizability component αϕϕ varies with (b). the angle, i.e., we require a spatially varying magnetic current to cancel the scattering from a PEC cylinder by (a) (b) using a metasurface. For a TEz plane wave incidence, the solution of canceling the scattered field is dual to the one just mentioned and the metasurface would have

mm )

) only surface polarizability component α . Figures 3(a) 1

1 zz -

- and (b) demonstrate the electric ﬁeld distributions in the Vsm

Asm absence and presence of the synthesized metasurface, re-

( (

 spectively, by using full-wave simulations based on the

z P

M ﬁnite element method [58]. It is clear from this Figure

(a) (b) metasurface 6 0  (degrees) 6 (c) -0.2 Re

-0.4 i 

Im / 0 PEC 0 PEC y

-0.6 EE

(m) log /

-0.8

mm 

0 -6 -6 -1  -6 0 6 -6 0 6 x / 

FIG. 3: The normalized (to the incident field amplitude  (degrees) Ei) electric field distribution around the proposed PEC FIG. 2: Induced surface (a) electric P and (b) cylinder illuminated by a plane wave in (a) the absence z and (b) the presence of the proposed scattering magnetic Mϕ polarization densities of Eqs. (29) and (30) required for scattering cancellation. (c) The cancellation cylindrical metasurface. The metasurface is corresponding effective magnetic surface polarizability shown with a black dashed line. mm αϕϕ desired for scattering cancellation from a PEC cylinder with radius a = 5λ. The values in plots are that the by placing the synthesized metasurface around similarly repeated for 180◦ ≤ φ ≤ 360◦ the PEC cylinder, the cylinder’s scattered fields have perfectly canceled, in other words the metasurface acts as a The last step is to retrieve the required effective polariz- cloaking device for a PEC cylinder. We next examine the ability tensors from Eqs. (1) and (2) that realize the po- cancellation of the scattering (i.e., cloaking) of a dielec- larization density M obtained in Eq. (30). However, in tric cylinder by covering it by a metasurface to further Eq. (2) the number of equations are less than the number demonstrate the capability of our proposed general ap- of unknown polarizability components. Therefore, the so- proach. lution to this synthesis problem is not unique since there are multiple polarizability tensors’ sets which lead to the same magnetic polarization density M, hence, to the B. Scattering cancellation from a dielectric desired electromagnetic field. Here, as previously men- cylinder tioned we simplify the synthesis procedure by neglect- ing the polarizability components normal to the meta- We consider now the problem of scattering cancella- surface plane (i.e., ρ-components). Moreover, for the tion from an infinitely long dielectric cylinder with rela- sake of simplicity, we avoid considering realizations that tive permittivity r = 10 and radius a = λ located in free possess bianisotropic polarizabilities α¯em and α¯me. Fur- space and illuminated by an electromagnetic plane wave thermore, we consider realizations with vanishing cross- as in Eq. (19) as in the previous case. There is a major component polarizability components, i.e., with zero off- difference between this case and the previous one with a diagonal terms in the polarizability tensor α¯mm. To sum- PEC cylinder. In the previous example, the electromag- 7 netic fields E− and H− inside the cylinder were required as well as the cross-component polarizabilities. Then, to be zero. Accordingly, based on the boundary con- by using the constitutive relation (1) the required effec- ee ditions (15)–(18), providing scattering cancellation en- tive electric polarizability αzz is derived and depicted forces both electric and magnetic polarization densities in Fig. 4(b). The electric field distributions (assuming P and M to be present at the boundary of the PEC [53– plane wave incidence) in the absence and presence of the 57]. In the present case with a dielectric cylinder, there designed metasurface are, respectively, demonstrated in are infinite solutions that solve the scattering cancellation Fig. 5(a) and (b) by using full-wave simulations. As it problem depending on the precise form which the equivalent principle is applied. Among them, we consider the (a) (b) metasurface one which can be realized when we impose that only the 3 3 0 electric polarization density P is non vanishing. There- -0.2 fore, assuming the magnetic polarization density M to be -0.4 i zero in Eqs. (15) and (16), the electromagnetic field in- 

/ 0 0 r EE side the dielectric cylinder can be easily obtained. Based y r -0.6 on the considered incident polarization, the electric field log / inside the metasurface has a z-component only. More- -0.8 over, to have no scattered fields out of the metasurface -3 -3 -1 (i.e., Eout = 0), we apply the boundary condition (15) -3 0 3 -3 0 3 using Eq. (21) to obtain the only surviving coefficient of x /  −n √ Eq. (27) as an = E0j Jn(β0a)/ Jn(β0 ra) . As a result, the required surface polarization density P is ob- FIG. 5: The normalized (to the incident field amplitude tained from the boundary condition (18) and using the Ei) electric field distribusions in (a) the absence and (b) ± ± Maxwell-Faraday equation ∇ × E = −jωµ0H at the the presence of the synthesized scattering calncellation metasurface boundary (i.e., ρ = a), leads to cylindrical metasurface in the case of a dielectric cylinder. The metasurface is shown with a black dashed line. n=+∞ E0 X −n jnφ 0 P = −zˆ j e [Jn(β0a) ωη0 is clear from the figure, when the metasurface is present n=−∞ (31) √ the scattered field is suppressed outside the metasurface rJn(β0a) 0 √ boundary. That is, outside the metasurface we observe − √ Jn(β0 ra) . Jn(β0 ra) only the propagating incident plane wave. In both the last examples, as seen in Figs. 2 and 4, The longitudinal polarization component is plotted in the imaginary part of the required effective polarizabil- Fig. 4(a) as a function of the azimuthal angle ϕ. Next ity tensors is negative at some azimuthal positions and positive for others. Therefore, one may deduce that the Re (a) -12 (b) desired metasurface consist of active elements. However, ×10 ×10-1 Im 3 4 as discussed in [17, 18], the pointing vector of the wave can switch the direction locally at the metasurface plane 2 and when its direction is opposite to the incident wave, -1 the polarizabilities can be negative. This implies the pos- 0 0 0

/ (m) sibility of realization of the synthesized structures using

(Asm )

ee zz

z only passive particles.  P -2

-3 -4 0 90 180 0 90 180 C. An illusory infinite line source  (degrees)  (degrees) In this example the goal is to create an electromag- FIG. 4: (a) The required polarization density for the netic illusion that a source is moved somewhere else by example of scattering cancellation metasurface for the surrounding the actual source with a metasurface. Stud- proposed dielectric cylinder. (b) The corresponding ies on the design of optical devices based on TO to create effective electric polarizability of the designed optical illusions have been already performed in past re- metasurface. The dielectric radius and relative cent years (see e.g. Ref. 51). However, besides proposing permittivity are a = λ and r = 10, respectively. bulky complex media, there was also the restriction that both the radiation source and its illusion were located we derive the polarizabilities that provide the required inside the engineered medium. These requirements make polarization density. For this simple case only one kind the realization of such devices difficult if not impracti- of polarizability is sufficient. Indeed, we assume that cal and also less useful. However, as we propose here, the bianisotropic polarizabilities α¯em and α¯me are absent the introduction of metasurfaces make the realization of 8 such devices more convenient and practical. Here, the such metasurface, and here we consider only realizations goal is to show how a cylindrical metasurface around a with non vanishing electric surface polarizations density line source is able to modify the source fields to make it and M assumed to be zero. As a result of the above con- look like it is generated by a fictitious line source out- siderations, and by using the electric fields (32) and (34) side the metasurface boundary, translated from the origin the boundary condition (15), the only surviving coef- inal source that is not visible anymore. In summary the ficients in Eq. (27) read metasurface cloaks the original line source and generates an illusory line source at a different location.  I β2 H(2)(β a)  e 0 0 0 0 Let us consider an infinitely long electric line source − [J0 (β0ρ ) − 1] n = 0  4ω0 J0(β0a) I in free space, along the z-axis of a cylindrical coordi-  e an = . (35) nate system. The electric field of such current source is  2 (2)  Ieβ0 Hn (β0a) 0 described by Eq. (24). The goal is to synthesize a cylin-  − Jn (β0ρ ) n 6= 0 drical metasurface around this line source that creates a 4ω0 Jn(β0a) total field that looks like the one generated by a fictitious Finally, by applying the boundary condition (18) and line source translated at (ρ0 > a, φ0), for an observer out- by using the Maxwell-Faraday equation ∇ × E± = side the metasurface [see Fig. 7]. Moreover, the distance −jωµ H± at the metasurface boundary (i.e., ρ = a), between the metasurface and the fictitious line source 0 we obtain the required electric polarization density P for must be subwavelength due to energy conservation con- the desired fields given by (32) and (34) as siderations. Indeed, the total flux of the Poynting vector through any closed surface around the illusory source n=+∞ Ieβ0 (2)0 X h 0 must be zero if such surface does not enclose any part P = −zˆ H (β a) − ejnφ J (β ρ0)H(2) (β a) 4ω 0 0 n 0 n 0 of the metasurface. This would be impossible if the fic- n=−∞ titious line source is located at large distance from the metasurface. 4an 0 + Jn(β0a) , We assume the fictitious line source has the same am- ωµ0Ie plitude Ie as the original line source but is located at (36) ρ = ρ0, φ = φ0. By using the addition theorem (see e.g., where an is given by Eq. (35). This surface polarization Refs. [48, 59]) the electric field created by such a source, density is used in the constitutive relation (1) to retrieve which corresponds to the desired total field outside the the required effective surface polarizability tensors. Fig- cylindrical metasurface, reads ures 6(a) and (b), respectively, illustrate the required surface polarization density and the needed effective electric surface polarizability for the proposed problem assuming 2 +∞ 0 + Ieβ0 X 0 (2) jn(φ−φ0) a = λ and ρ = 4λ/3. Similarly to what was done in E = −zˆ Jn(β0ρ )Hn (β0ρ)e . (32) 4ω0 −∞ (a) (b) Re ×10-9 ×102 Im Comparing the above equation with Eq. (28) that repre- 5 6 sents the total scattered electric field outside the meta- Sampling points surface, one realizes that the only surviving coefficient in

-1 3 Eq. (28) is

0 0 / (m)

2 (Asm ) 0 ee

I β zz e 0 0 z

c = − J (β ρ ) . (33) P

n n 0  4ω0 -3 The total electric ﬁeld inside the cylindrical metasur- -5 0 90 180 0 90 180 face reads  (degree)  (degree)

2 +∞ Ieβ (2) X 0 E− = −zˆ 0 H (β ρ) + zˆ a J (β ρ)ejn(φ−φ ), FIG. 6: (a) The required polarization density for the 4ω 0 0 n n 0 metasurface example that creates an infinite illusory 0 −∞ (34) line source. (b) The corresponding effective electric which is the contribution of the field created by the polarizability of the designed metasurface. original line source [see Eq. (24)] plus the field created from the cylindrical metasurface [see Eq. (27)]. Note the previous example, we again have considered a meta- that the metasurface shall not create any cross compo- surface realization that does not involve cross-component nent scattered field, i.e., the ρ- and φ-components of the and bianisotropic polarizabilities. Figures 7(a) and (b), electric field simply vanish in Eq. (27). Now we assume respectively, show the full-wave simulation results for the the fictitious source is located at ρ0 and φ0 = 0. More- electric field distributions in the absence and the pres- over, we assume the observer is located at distances that ence of the engineered metasurface. It is clear from the satisfy ρ > ρ0. There are many possible realizations of field distributions that the original line source is cloaked 9

(a) (b) (a) line source (b) 5 5 0 5 0 line source metasurface -0.5 metasurface

-1 -1



 / / -1.5 PEC 0

0 0 E

y E y z

log / y

-2 -2 log / x original line illusory -2.5 line source -5 -5 source line source -3 −5 illusion -3 -5 0 5 -5 0 5 0 5 x /  ~ x /  2 FIG. 7: (a) The normalized [to EA = Ieβ0 / (4ω0)] FIG. 8: (a) The schematic of an exemplary practical electric field distribution around an infinitely long realization of the discretized metasurface that creates electric line source in vacuum. (b) The normalized an illusory line source next to the original one that not electric field distribution of the original line source 2 visible. (b) The normalized [to EA = Ieβ0 / (4ω0)] when located at the center of a metasurface that is electric field distribution for the structure given in (a). synthesized to translate the perception of the position of the original line source for an observer outside the cylindrical metasurface. The matasurface radius is a = λ, the illusory source is translated to the position ρ0 = 4λ/3 and φ0 = 0, and the field appears arising from such source. the sampling points. In order to find the necessary load impedance values we first create a cylindrical array of 20 identical short dipole antennas located at each sampling and an illusory line source is observed at 4λ/3 away from point. These points are given by the angular positions the original line source for an observer outside the cylin- 0 0 of the extrema values given in Fig. 6(a). Each wire in der with radius ρ , i.e., for ρ > ρ . Such a device may the plates scatters and couples to all the others, modify- find practical applications in stealth technology, wireless ing their currents and voltages. Therefore by looking at power transfer systems [60], and also sensing technolo- the voltage and current at each load we relate all these gies. quantities with a 20 × 20 impedance matrix Z¯ in which the diagonal elements denote the self impedance of each VI. PRACTICAL REALIZATION antenna and the off-diagonal elements represent the mu- tual couplings between the antennas. The goal is to design the loads at each wire to ensure that the current In all the above examples we have performed concep- distribution is similar to the one shown in Fig. 6(a). For tual syntheses by considering ideal continuous current convenience we define the load impedance diagonal ma- distributions, i.e., continuous surface polarization densi- trix Z¯ , which contains the load impedances Z at each ties. However, a metasurface is practically composed of L Li sampling points φ , i = 1, ..., 20 shown in Fig. 9. The discrete elements that mimic such continuous currents. i scattering problem is now posed as an algebraic prob- In this section, we present a practical approach to real- h ¯ ¯ i ize such metasurfaces with discrete elements, and apply lem as V = Z + ZL I, where I is a 20 element vector it to the last example, i.e., the illusory line source. We with its i-th element Ii = jωaP(φ − φi) · zˆ given by the start by taking samples along the azimuthal angular di- effective surface electric polarization density at the sam- rection that corresponds to the extrema values of the de- pling point φi [Fig. 6(a)] and V is a voltage (per meter) sired continuous current distribution given in Fig. 6(a). vector applied to dipoles. Considering identical supply- According to Fig. 6(a) and considering the whole 360◦ ing voltage 1 (V/m) for all the antennas, the required azimuth angle, we obtain 20 sampling points. Note that load impedances are given in Fig. 9. Here we have pro- only half of the whole 360◦ angular span is shown in vided the impedance values for half of the elements since Fig. 6(a) since the other half of the cylinder is symmetric the values for the other half is symmetrically distributed to the one shown and hence it exhibits the same values. [see Fig. 9]. The electric field distribution of such dis- Next, each of these required currents can be realized by cretized problem is depicted in Fig. 8(b) by using finite wires with periodic set of loads located at the azimuthal element method (FEM) implemented in CST. As it is sampling points. For simulation purposes we restrict our clear from this figure, although the field distribution is analysis to a single periodic unit as shown in Fig. 8 and not perfectly cylindrical due to the discretization, the representing all the others by reflections of two parallel actual line source inside the cylindrical metasurface is perfectly conducting plates. It is noteworthy that the cloaked by the metasurface and a fictitious line source coupling between the dipoles in the parallel plate region is inspected by an observer outside the metasurface at are tuned by the loads to obtain the desired currents at electromagnetic large distances. 10

ZL6 ACKNOWLEDGMENT ZL7 ZL5 Z ZL4 L8 illusory The authors would like to thank DS SIMULIA for pro- ZL3 line ZL9 viding CST Studio Suite that was instrumental in this line source Z study, and Prof. C. Simovski and Prof. S. Tretyakov, L10 source ZL2 from Aalto University, Finland, for fruitful discussions. j ZL11 ZL1 Z Z Z L1 L2 L3 Appendix A: Boundary conditions j (deg.) 0 19 41 Z (W) -0.11+j0.53 -0.14+j0.03 0.29+j0.3 L Let us consider the proposed cylindrical metasurface ZL4 ZL5 ZL6 of radius ρ = a whose axis coincides with the z-axis (see j (deg.) 58 74 95 Fig. 1). We start from the Maxwell’s curl equations with electric and magnetic polarization densities Pv = Pδ(ρ− ZL(W) -0.26+j0.22 0.32+j0.5 -0.21+j0.09 a) and Mv = Mδ(ρ − a) located at ρ = a, i.e., ZL7 ZL8 ZL9 j (deg.) 110 127 148 ∇ × E = −jω (µ0H + Mv) , (A1)

ZL(W) -0.25+j0.17 -0.22+j0.06 -0.25+j0.34

ZL10 ZL11 j (deg.) 160 180 ∇ × H = jω (0E + Pv) , (A2)

ZL(W) -0.25+j0.25 0.46+j0.30 In the Eqs. (A1) and (A2), µ0 and 0 are the free FIG. 9: The relative positions and load impedance space magnetic permeability and electric permittivity, values of the dipoles in the proposed cylindrical respectively, and we consider the constitutive relations metasurface that creates an illusory line source. Here, D = 0E + Pv and B = µ0H + Mv between the electric we show 11 elements while we declare that elements displacement D, magnetic induction B, electric ﬁeld E, number 2 to 10 are repeated symmetrically in the lower magnetic ﬁeld H, electric polarization density Pv, and half circle. magnetic polarization density Mv. Next, since the problem has cylindrical symmetry, it is convenient to represent all vectors in cylindrical coordinates as ˆ E = Eρρˆ + Eφφ + Ezzˆ, (A3)

ˆ VII. DISCUSSION AND CONCLUSIONS H = Hρρˆ + Hφφ + Hzzˆ, (A4)

ˆ Pv = Pvρρˆ + Pvφφ + Pvzzˆ, (A5) We have investigated the synthesis problem of cylindrical metasurfaces for general illusion mechanisms, and M = M ρˆ + M φˆ + M zˆ. (A6) have presented a general analysis and synthesis approach v vρ vφ vz by modeling metasurfaces with eﬀective polarizability tensors. We have presented three practical examples in- By using (A3)–(A6) in the Maxwell curl equations (A1) cluding scattering cancellation from PEC and dielectric and (A2) in cylindrical coordinate system we have cylinders as cloaking examples. Furthermore we have also shown how to generate an illusion, a translated ﬁctitious 1 ∂Ez ∂Eφ − = −jω (µ0Hρ + Mvρ) , (A7) line source, from an original line source by covering such ρ ∂φ ∂z a line source with a cylindrical metasurface. Besides the mentioned applications, the proposed ap- ∂Eρ ∂Ez proach is general and can be applied to many other − = −jω (µ0Hφ + Mvφ) , (A8) problems. That is, any manipulation of electromagnetic ∂z ∂ρ waves including but not restricted to focusing lenses, de- signing antennas with arbitrary radiation patterns, beam 1 ∂ (ρEφ) 1 ∂Eρ forming, etc. Our study opens up possibilities for com- − = −jω (µ0Hz + Mvz) , (A9) ρ ∂ρ ρ ∂φ plex electromagnetic wave manipulations using cylindrical topologies including stealth technology, antenna synthesis problems, wireless power transfer systems, ab- 1 ∂H ∂H z − φ = jω ( E + P ) , (A10) sorbers, etc. ρ ∂φ ∂z 0 ρ vρ 11

∂H ∂H ρ − z = jω ( E + P ) , (A11) ∂z ∂ρ 0 φ vφ 1 1 ∂Ez ∂Eφ Mvρ Hρ = − − − . (A14) 1 ∂ (ρH ) 1 ∂H jωµ0 ρ ∂φ ∂z µ0 φ − ρ = jω ( E + P ) . (A12) ρ ∂ρ ρ ∂φ 0 z vz

Then, we express the normal (with respect to the metasurface boundary) components of the electric and One can obtain the boundary conditions for cylindri- magnetic ﬁelds in terms of their tangential counterparts cal metasurfaces given in (3)–(6), by substituting (A13) from (A10) and (A7), respectively, as into Eqs. (A8) and (A9), and (A14) into Eqs. (A11) and (A12), and then integrating the resulting equations over 1 1 ∂Hz ∂Hφ Pvρ Eρ = − − , (A13) the metasurface thickness along the ρ-direction. jω0 ρ ∂φ ∂z 0

[1] E. F. Kuester, M. Mohamed, M. Piket-May, and C. Hol- [18] V.S. Asadchy, M. Albooyeh, S.N. Tcvetkova, A. D´ıaz- loway, IEEE Trans. Antennas Propag. 51, 2641 (2003). Rubio, Y. Ra’di, and S.A. Tretyakov, Phys. Rev. B 94, [2] C. Holloway, A. Dienstfrey, E. F. Kuester, J. F. OHara, 075142 (2016). A. K. Azad, and A. J. Taylor, Metamaterials 3, 100 [19] A. Epstein and G. V. Eleftheriades, J. Opt. Soc. Am. B (2009). 33, A31-A50 (2016). [3] C. Holloway, E. F. Kuester, J. Gordon, J. OHara, J. [20] J. PS Wong, A. Epstein, and G. V. Eleftheriades, IEEE Booth, and D. Smith, IEEE Antennas Propag. Mag. 54, Antennas Wirel. Propag. Lett. 15, 1293 (2016). 10 (2012). [21] M. Danaeifar, N. Granpayeh, N. A. Mortensen, and S. [4] The term “metasurface” is referred to surfaces with Xiao, J. Phys. D 48, 385106 (2015). negligible electrical (optical) [61] thickness [3]. There- [22] T. Niemi, A. Karilainen, and S. Tretyakov, IEEE Trans. fore, cylindrical surfaces as considered here (with neg- Antennas Propag. 61, 3102 (2013). ligible thickness in the radial direction) are considered [23] M. A. Salem and C. Caloz, Opt. Express 22, 14 530 in the category of metasurafces although it is a three- (2014). dimensional object. [24] H. Kazemi, M. Albooyeh, and F. Capolino, in URSI EM [5] H. Kazemi, M. Albooyeh and F. Capolino, Theory Symposium, EMTS 2019, San Diego, CA, 2731 arXiv:1905.04439 (2019). May 2019. [6] W. J. Padilla, M. T. Aronsson, C. Highstrete, M. Lee, A. [25] K. Achouri, M. A. Salem, and C. Caloz, IEEE Trans. J. Taylor, and R. D. Averitt, Phys. Rev. B 75, 041102(R) Antennas Propag. 63, 2977 (2015). (2007). [26] M. Albooyeh, D. Morits, and C. R. Simovski, Metama- [7] M. Caiazzo, S. Maci, and N. Engheta, IEEE Antennas terials 5, 178 (2011). Wirel. Propag. Lett. 3, 261 (2004). [27] M. Albooyeh, S. Tretyakov, and C. Simovski, Ann. der [8] W. L. Chan, H.-T. Chen, A. J. Taylor, I. Brener, M. J. Phys. 528, 721 (2016). Cich, and D. M. Mittleman, Appl. Phys. Lett. 94, 213511 [28] C.M. Roberts, S. Inampudi, and V.A. Podolskiy, Opt. (2009). Express 23, 2764 (2015). [9] C. Pfeiffer and A. Grbic, Phys. Rev. Applied 2, 044011 [29] M. Kiani, A. Abdolali, and M.M. Salary, IEEE Trans. (2014). Antennas Propag. 63, 631118 (2015). [10] M. Selvanayagam and G. Eleftheriades, IEEE Trans. An- [30] H. Oraizi and A. Abdolali, Prog. Electromagn. Res. B 3, tennas Propag. v, 6155 (2014). 227 (2008). [11] F. Capolino, A. Vallecchi, and M. Albani, IEEE Trans. [31] R. Alaee, C. Menzel, C. Rockstuhl, and F. Lederer, Opt. Antennas Propag. 61, 852 (2013). Express 20, 18370 (2012). [12] C. Pfeiffer, C. Zhang, V. Ray, L.J. Guo, and A. Grbic, [32] A. A. Kishk and P-S. Kildal, IEEE Trans. Antennas Phys. Rev. Lett. 113, 023902 (2014). Propag. 45, 51 (1997). [13] Y. Zhao, N. Engheta, and A. Alú,Metamaterials 5, 90 [33] S. Raffaelli, Z. Sipus, and P-S. Kildal, IEEE trans An- (2011). tennas Propag. 53, 1105 (2005). [14] M. Albooyeh, R. Alaee, C. Rockstuhl, and C. Simovski, [34] Y.R. Padooru, A.B. Yakovlev, Pai-Yen Chen, and An- Phys. Rev. B 91, 195304 (2015). drea Alú,J. Appl. Phys. 112, 034907 (2012). [15] M. Albooyeh, V. Asadchy, J. Zeng, H. Kazemi, and F. [35] Y.R. Padooru, A.B. Yakovlev, Pai-Yen Chen, and A. Alú, Capolino, arXiv:1811.04176 (2018). J. Appl. Phys. 112, 104902 (2012). [16] M. Albooyeh, Y. Ra’di, M.Q. Adil, and C.R. Simovski, [36] H.M. Bernety, and A.B. Yakovlev, IEEE Trans. Antennas Phys. Rev. B 88, 085435 (2013). Propag. 63, 1554 (2015). [17] V.S. Asadchy, A. D´ıaz-Rubio, S.N. Tcvetkova, D.H. [37] B. O. Raeker and S. M. Rudolph, IEEE Antennas Wirel. Kwon, A. Elsakka, M. Albooyeh, and S. A. Tretyakov, Propag. Lett. 15, 1101 (2016). Phys. Rev. X 7, 031046 (2017). [38] B. O. Raeker and S. M. Rudolph, IEEE Antennas Wirel. 12

Propag. Lett. 16, 995 (2017). (John Wiley & Sons, Danvers, MA, 2012). [39] J. Xiao, Y. Vahabzadeh, C. Caloz, Y. Fan, IEEE Trans. [49] P. Y. Chen and A. Alú,Phys. Rev. B, 84, 205110 (2011). Antennas Propag. 67, 2542 (2019). [50] D. Schurig, J.J. Mock, B.J. Justice, S.A. Cummer, J.B. [40] M. Dehmollaian, N. Chamanara, C. Caloz, IEEE Trans. Pendry, A.F. Starr, and D.R. Smith, Science 314, 977 Antennas Propag. doi: 10.1109/TAP.2019.2905711 (2006). (2019). [51] J. Yi, P.H. Tichit, S.N. Burokur, and A. de Lustrac, J. [41] S.A. Tretyakov Analytical Modeling in Applied Electro- Appl. Phys. 117, 084903 (2006). magnetics (Artech House, Norwood, MA, 2003). [52] S. Tretyakov, P. Alitalo, O. Luukkonen, and C. Simovski, [42] A. Serdyukov, I. Semchenko, S. Tretyakov, and A. Si- Phys. Rev. Lett. 103, 103905 (2009). hvola, Electromagnetics of Bi-anisotropic Materials: [53] D.H. Kwon and D. M. Pozar, IEEE Trans. Antennas Theory and Applications (Gordon and Breach Science Propag. 57, 3720 (2009). Publishers, Amsterdam, 2001). [54] M. Albooyeh and C. R. Simovski, Opt. Express 20, 21888 [43] M. Idemen, Universal boundary conditions and Cauchy (2012). data for the electromagnetic field in Essays on the formal [55] Y. Ra’di, V.S. Asadchy, and S.A. Tretyakov, IEEE Trans. aspects of electromagnetic theory (World Scientific, New Antennas Propag 61, 4606 (2013). Jersey, 1993), pp. 657–698. [56] Y. Ra’di, C.R. Simovski, and S.A. Tretyakov, Phys. Rev. [44] L. B. Felsen and N. Marcuvitz, Radiation and Scattering Appl. 3, 037001 (2015). of Waves (Wiley, New York, 1994), pp. 185–187. [57] R. Alaee, M. Albooyeh, and C. Rockstuhl, J. Phys. D: [45] M. Albooyeh, D.H. Kwon, F. Capolino, and S.A. Appl. Phys. 50, 503002 (2017). Tretyakov, Phys. Rev. B 95, 115435 (2017). [58] https://www.comsol.com/ [46] M. Albooyeh, H. Kazemi, F. Capolino, D.H. Kwon, and [59] R. F. Harrington, Time-Harmonic Electromagnetic S.A. Tretyakov, in 2017 IEEE International Symposium Fields (McGraw-Hill, New York, 1961). on Antennas and Propagation & USNC/URSI National [60] M. Song, P. Belov, and P. Kapitanova, Appl. Phys. Rev. Radio Science Meeting, pp. 1707-1708 (2017). 4, 021102 (2017). [47] Note that although there is no theoretical limitation to [61] The term ”electrically (optically) thin” refers to the realize any desired electromagnetic field with only tange- thickness as compared to the wavelength. tial (to the surface) polarization densities, there would [62] M. Safari, A. Abdolali, H. Kazemi, M. Albooyeh, M. be extra practical degrees of freedom in realization of a Veysi, and F. Capolino, in 2017 IEEE International Sym- desired field using normal polarization densities as dis- posium on Antennas and Propagation & USNC/URSI cussed in Refs. [45, 63]. National Radio Science Meeting, pp. 14991500 (2017). [48] C.A. Balanis, Advanced Engineering Electromagnetics [63] C. Pfeiffer and A. Grbic, Phys. Rev. Lett. 117, 077401 (2016).