## Abstract

Omni-directional electromagnetic field concentrators have been recently reported by Marco Rahm et al. [Photon. Nanostruct.: Fundam. Applic. 6, 87 (2008)] based on form-invariant coordinate transformations related to its Jacobi transformation matrix. Using transverse-electric wave illumination, we reduced the complex material parameters of the concentrator for future practical implementation. Concentrators with different set of permittivity and permeability tensors are proposed. The electromagnetic concentrating performance and the scattering properties at the inner and outer boundary of these concentrators are theoretically and numerically analyzed. Finally we obtain a set of material tensors for a concentrator that simultaneously has perfect matched interior and exterior interfaces.

©2008 Optical Society of America

## 1. Introduction

Pendry et al. [1] has proposed an interesting idea of using geometrical transformation combined with Maxwell’s equations to yield unique electromagnetic devices, recently termed “invisibility cloaks” with the electrical permittivity and magnetic permeability both anisotropic and spatially varying, in which the electromagnetic wave can be smoothly guided around a “cloaked” region without disturbing the exterior fields. This new method of controlling electromagnetic fields have been theoretically verified soon by D. Schurig [2] using Cartesian tensors directly to calculate material properties of the transformed medium and by U. Leonhardt [3] who transformed the two dimensional Helmholtz equation to produce similar cloaking effects in the geometrical optic limit. Based on coordinate transformation, other electromagnetic structures with special electromagnetic functionality have been reported such as rotation coatings [4], magnifying perfect lens [5, 6], beam refactors [7], square cloaks and concentrators [8, 9].

Although coordinate transformation associated with Jacobian matrix enables us easily to obtain the material parameters and the development of artificially structured metamaterials provides the flexibility of tuning the electrical permittivity and magnetic permeability independently and arbitrarily covering positive and negative values as desired [10–13], obtaining the tensor components that are generally inhomogeneous and spatially dependent is still challenging for practical fabrication. Most papers reported to date focus on theoretical structural design and numerical analysis [14–18]. Fortunately, transformed material parameters can be simplified by applying EM wave with specified polarization direction, which significantly facilitates the implementation of designed structures in reality. A 2D cloak device with reduced permittivity and permeability tensors have been experimentally fabricated by assembling split ring resonators (SRRs) in desired way [19].

In this paper, using the traditional geometrical transformation approach as first proposed by Pendry for its convenience and intuitiveness in handling co-ordinate systems that are orthogonal, such as cylindrical and spherical, we focus on the concentrator first proposed in Reference [8]. It is noted that the concentrator in the ideal case has the capability of perfectly concentrating the incident EM power without any refection and absorption, but the spatially varying material tensor components are difficult for practical realization. So we derive different sets of reduced material specifications for cylindrical concentrators by taking advantage of TE incident wave and systematically analyze their concentrating performance as well as the scattering property. A set of simpler reduced material tensors is determined to make both the inner and the outer interface perfectly matched, 2D numerical simulations are also performed to confirm our theoretical analysis and conclusions.

## 2. The ideal concentrator

In cylindrical coordinate system, the concentrator experiences a two-step space deformation as Marco Rahm suggested in Reference [8]:

Space in the region of 0 ≤ *r* ≤ *R*
_{2} is compressed into the inner circle, and simultaneously space between *R*
_{2} and *R*
_{3} is stretched to fill the *R*
_{1} ≤ *r* ≤ *R*
_{3} area, the transformation can be mathematically denoted in the following piecewise equation:

According to Pendry’s geometrical transformation method, the transformed material tensor components from Cartesian system to the cylindrical have the following form:

Where ${Q}_{i}=\sqrt{{\left(\frac{\partial x}{\partial {q}_{i}}\right)}^{2}+{\left(\frac{\partial y}{\partial {q}_{i}}\right)}^{2}+{\left(\frac{\partial z}{\partial {q}_{i}}\right)}^{2},}\left(i=\mathrm{1,2,3}\right)$, denotes Lamé coefficient with *q _{i}* representing three coordinate curves of cylindrical system. We should calculate the corresponding Lamé coefficients for the two regions in Equation(1) respectively:

for 0 ≤ *r* ≤ *R*
_{1}, we have

and for *R*
_{1} < *r* ≤ *R*
_{3} region,

Combining Eqs. (4) and (5) with Eqs. (2) and (3), the diagonal material tensor of the permittivity and permeability for the two Cartesian regions are:

here, *M* = (*R*
_{2} - *R*
_{1})*R*
_{3}/(*R*
_{3} - *R*
_{2}) and *N* = [(*R*
_{3} - *R*
_{2})/(*R*
_{3} - *R*
_{1})]^{2} .The results are exactly equivalent to those derived by Marco Rahm.

## 3. Concentrators with reduced parameters

It is noted that in 0 ≤ *r* ≤ *R*
_{1} region material tensor components are all constants and only the z-component need to be considered for the field polarization oriented along the z-axis for future experimental implementation, but the tensor for *R*
_{1} ≤ *r* ≤ *R*
_{3} region shows its components that are all dependent on the radial position, indicating tremendous trouble in satisfying them all at once. Here we use TE incidence as first proposed by Pendy et al. for material simplification in practical realization of cylindrical cloaking structure [19] to reduce material specification in the concentrator. Given the incident plane wave polarized along z-direction propagates in positive x axis, electric field has only one component, which can be expressed in terms of (*r*,*ϕ*,*z*) basis as *E*̄ = *z*̂*E _{z}*. Combining with Maxwell’s equations the wave equation determining the propagation property of EM wave inside the concentrating structure can be derived in the form of cylindrical system:

We note from Eq. (7) that only the radial permeability component μ_{ϕ}, azimuthal permeability component μ_{ϕ} and z-component ε_{z} are involved in the wave equation. We choose a set of material parameters for cylindrical concentrators that can maintain the same wave propagation behavior in the region of *R*
_{1} < *r* ≤ *R*
_{3} , e.g., the reduced parameters:

in which we realize that *ε _{z}* and

*μ*are greatly simplified with the radial permeability specification being position-relevant. Parameters presented in Equation (8) minimize the most challenging part of experimental design. But this set of parameter specification leads to the impedance mismatch at the outer boundary of the concentrator. The ideal parameters in the region of

_{ϕ}*R*

_{1}<

*r*≤

*R*

_{3}in Equation (6) have the impedance of ${Z}_{r={R}_{3}}=\sqrt{\frac{{\mu}_{\varphi}}{{\epsilon}_{z}}}=1$ at the outer surface

*r*=

*R*

_{3}surface and

*Z*

_{r=R1}=

*R*

_{1}/

*R*

_{2}= at the inner surface

*r*=

*R*

_{1}that are both matched at these two interfaces, while the simplified specification corresponds to the mismatched impedance of ${Z}_{r={R}_{3}}={Z}_{r={R}_{1}}=\sqrt{\frac{1}{N}}=\frac{\left({R}_{3}-{R}_{1}\right)}{\left({R}_{3}-{R}_{2}\right)}$ at the two interfaces, evidently the energy loss due to these impedance mismatch will degenerate the concentrating performance. To improve the concentrating efficiency, we derive another set of simplified material tensors, which also have simple tensor components and can make the outer boundary perfectly matched to the free space. With TE wave mode, the wave equation still remain unchanged if we choose the following material components for the

*R*

_{1}<

*r*≤

*R*

_{3}region

In this set of parameters, only the radial permeability is spatially varying, and the region of *R*
_{1} < *r* ≤ *R*
_{3} becomes a perfectly matched layer with a constant impedance of ${Z}_{{R}_{1}\sim {R}_{3}}=\sqrt{\frac{{\mu}_{\varphi}}{{\epsilon}_{z}}}=1$. But at the inner boundary the impedance mismatch is not eliminated. So we can further reduce the parameters for the 0 ≤ *r* ≤ *R*
_{1} region to obtain a zero-reflectance inner interface in the TE mode case, the parameters experience the following simplification:

Actually, in the case of z-polarized wave the simplified material parameters indicate an isotropic material within the inner circular cylinder with the permittivity and permeability being both constants, and the inner surface as well as the outer one are both matched interface without reflection. Now we have proposed three kinds of reduced material specification for the concentrator, their simple forms are expected to provide the ease for the future practical realization. The different impendence patterns at the interfaces give rise to the variation of their concentrating capabilities, which will later be demonstrated by 2D simulations.

## 4. Simulation and discussion

In this section, finite element method is used to simulate the designed concentrators. We use circular concentrators that have the same size as used in Reference [8], where *R*
_{1} = 2*cm*, *R*
_{2} = 4*cm* and *R*
_{3} = 6*cm*. TE wave with working frequency of *f* = 15*GHz* is incident along *x* axis throughout the simulation area from left to right. All computational domain boundaries are perfectly matched layers (PML) so that the scattered fields can be absorbed. We first simulate the 2D concentrator with ideal parameters to confirm our previous analyses, then we illustrate the simulated results of those with different set of reduced parameters, their scattering properties and concentrating abilities are discussed.

Although the concentrators with ideal material parameters have excellent concentrating performance, it is actually difficult to realize the complicated tensor components in practice. Here we numerically demonstrate three kinds of concentrators having reduced material specifications. The first kind is the one with material tensors presented in Eq. (8), which has simplest tensor components but non-zero reflection occurs at both the inner and the outer interface, we denote this kind of concentrator as case (1), the second kind, denoted as case (2) is the one possessing the material parameters given by Eq. (9), in this case only the outer boundary is perfectly matched. The third kind of concentrator is case (3), where Equation (9) and (10) is combined to specify the simplified material in region of *R*
_{1} < *r* ≤ *R*
_{3} and 0 ≤ *r* ≤ *R*
_{1} respectively, eliminating the impendence mismatches at the two boundaries. The simulation results are shown in Fig. (2).

The E field distributions in Fig. (2) clearly demonstrate the scattering properties of these reduced concentrators. In case (1), the field outside the concentrator is disturbed by the scattered wave from the outer boundary, and inside the inner circle, the wave front that is flat in the ideal cases deforms in shape due to perturbation of the scattering occurring at the inner interface. The energy coming into the structure from the left side is no longer equals to that flowing out from the right side, which surely decreases the concentrating efficiency of the concentrator. The same situation appears in case (2), but the scattering wave amount decreases as the simplified material parameters have made the outer boundary perfectly matched, we can see from the middle graph in Fig. (2) that the field distributed in the region between *R*
_{1} and *R*
_{3} is only slightly disturbed. This kind of structure is expected to have better performance for its less energy loss. In case (3), where reflections at both interfaces are eliminated, the field almost has the same pattern when compared with that in the ideal concentrators. We can expect great improvement in the concentrating functionality for its least scattering. The following time-averaged energy densities give a more intuitive illustration of the concentrating property for each case.

As shown in Fig. (3), the energy density is inhomogeneously distributed within the interior circular region in case (1) and (2), the combination of the incident wave with the scattering part strengthens the field energy at some places, and weaken the power at other locations, while the energy intensity in the third case is almost invariant with the positions inside the inner circle in spite of slightly energy loss in the area between the two interfaces. This kind of concentrator with simpler material specification and good concentrating ability is believed to be preferable in the future practical implementation.

## 5. Conclusion

In summary, cylindrical concentrators with different sets of simplified material parameters have been proposed using TE polarized wave incidence. Spatial dependence of the material tensor components was eliminated, dramatically reduced the difficulty in future practical implementation. Their abilities to concentrate incident field energy in the inner circle and the scattering properties at the inner and outer boundary have also been theoretically and numerically analyzed. We believe that the reduced concentrator with matched interior and exterior interfaces would provide the convenience for the realization in practice.

## Acknowledgment

This work was supported by 973 Program of China (No.2006CB302900) and National Natural Science Foundation of China (No.60507014).

## References and links

**1. **J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science **312**, 1780-1782 (2006) [CrossRef] [PubMed]

**2. **D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in
transformation media,” Opt. Express **14**, 9794–9804 (2006), http://www.opticsinfobase.org/abstract.cfm?URI=oe-14-21-9794. [CrossRef] [PubMed]

**3. **U. Leonhardt, “Optical conformal mapping,” Science **312**,1777–1780 (2006). [CrossRef] [PubMed]

**4. **H. Y. Chen and C. T. Chan, “Transformation media that rotate electromagnetic fields,” Appl. Phys. Lett. **90**, 241105 (2007). [CrossRef]

**5. **M. Tsang and D. Psaltis, “Magnifying perfect lens and superlens design by coordinate transformation,” Phys. Rev. B **77**, 035112 (2008). [CrossRef]

**6. **W. Wang, L. Lin, X. Yang, J. Cui, C. Du, and X. Luo, “Design of oblate cylindrical perfect lens using coordinate transformation,” Opt. Express **16**, 8094–8105 (2008), http://www.opticsinfobase.org/abstract.cfm?URI=oe-16-11-8094. [CrossRef] [PubMed]

**7. **L. Lin, W Wang, C. Du, and X. Luo, “Design of electromagnetic refractor and phase transformer using coordinate transformation theory,” Opt. Express **16**, 6815–6821 (2006), http://www.opticsinfobase.org/abstract.cfm?URI=oe-16-10-6815. [CrossRef]

**8. **M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, D. R. Smith, and J. B. Pendry, “Design of electromagnetic cloaks and concentrators using form-invariant coordinate transformations of Maxwell’s equations,” Photon. Nanostruct. Fundam. Appl. **6**, 87 (2008). [CrossRef]

**9. **L Lin, W. Wang, C. Wang, C. Du, and X. Luo, “A cone-shaped concentrator with varying performances of concentrating,” Opt. Express **16**, 6809–6814 (2006), http://www.opticsinfobase.org/abstract.cfm?URI=oe-16-10-6809. [CrossRef]

**10. **A. A. Houck, J. B. Brock, and I. L. Chuang, “Experimental observations of a left-handed material that obeys Snell’s Law,” Phys. Rev. Lett. **90**, 137401 (2003). [CrossRef] [PubMed]

**11. **R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental verification of a negative index of refraction,” Science **292**, 77–79 (2001). [CrossRef] [PubMed]

**12. **D. R. Smith and N. Kroll, “Negative refractive index in left-handed materials,” Phys. Rev. Lett. **85**, 2933–2936 (2000). [CrossRef] [PubMed]

**13. **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**, 4184–4187 (2000). [CrossRef] [PubMed]

**14. **Z. Ruan, M. Yan, C. W. Neff, and M. Qiu, “Ideal cylindrical cloak: perfect but sensitive to tiny perturbations,” Phys. Rev. Lett. **99**, 113903 (2007). [CrossRef] [PubMed]

**15. **F. Zolla, S. Guenneau, A. Nicolet, and J. B. Pendry, “Electromagnetic analysis of cylindrical invisibility cloaks and the mirage effect,” Opt. Lett. **32**, 1069–1071 (2007). [CrossRef] [PubMed]

**16. **W. Cai, U. K. Chettiar, A. V. Kildishev, and V. M. Shalaev, “Optical cloaking with metamaterials,” Nature Photonics **1**, 224 (2007). [CrossRef]

**17. **H.Y. Chen, Z. X. Liang, P. J. Yao, X.Y. Jiang, H. G. Ma, and C.T. Chan, “Extending the bandwidth of electromagnetic cloaks,” Phys. Rev. B **76**, 241104(R) (2007). [CrossRef]

**18. **S. A. Cummer, B. I. Popa, D. Schurig, D. R. Smith, and J. B. Pendry, “Full-wave simulations of electromagnetic cloaking structures,” Phys. Rev. E **74**, 036621 (2006). [CrossRef]

**19. **D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science **314**, 977–980 (2006). [CrossRef] [PubMed]