In transformation optics, the space transformation is viewed as the deformation of a material. The permittivity and permeability tensors in the transformed space are found to correlate with the deformation field of the material. By solving the Laplace’s equation, which describes how the material will deform during a transformation, we can design electromagnetic cloaks with arbitrary shapes if the boundary conditions of the cloak are considered. As examples, the material parameters of the spherical and elliptical cylindrical cloaks are derived based on the analytical solutions of the Laplace’s equation. For cloaks with irregular shapes, the material parameters of the transformation medium are determined numerically by solving the Laplace’s equation. Full-wave simulations based on the Maxwell’s equations validate the designed cloaks. The proposed method can be easily extended to design other transformation materials for electromagnetic and acoustic wave phenomena.
© 2009 Optical Society of America
CorrectionsJin Hu, Xiaoming Zhou, and Gengkai Hu, "Design method for electromagnetic cloak with arbitrary shapes based on Laplace’s equation: erratum," Opt. Express 17, 13070-13070 (2009)
Coordinate transformation method for designing electromagnetic materials with a prescribed function has received much attention [1-5]. Based on the form invariance of Maxwell’s equations under a coordinate transformation, a topological variation in the transformed space is effectively equivalent to the variation of material parameters in the original space. Through a specific transformation, Pendry et al  firstly design the well-known spherical and cylindrical cloaks capable of guiding electromagnetic waves around embedded obstacles. The theory is immediately verified by a microwave experiment  with help of the metamaterial technology , thus opening a new approach by utilizing materials to control electromagnetic waves at will. The transformation method is a general method whatever the shape and function of the transformed media [2, 7, 8], however the determination of the corresponding material parameters for a given shape and function is a formidable task, which is still an active topics in this field. For cloaks with regular forms, such as square  and elliptical [10-13] shapes, they can be designed without much difficulty by constructing explicitly the transformation tensor. To design cloaks with irregular shapes, there are different techniques recently proposed based on the transformation theory [14-17]. Yan et al.  theoretically show that an arbitrary cloak is possible by analyzing wave property at cloak boundaries. Jiang et al.  use non-uniform B-spline (NURBS) functions to describe the geometrical boundary of a cloak with an irregular shape. Nicolet et al.  employ finite Fourier expansions to characterize a random boundary. After a linear mapping, the material parameters of an arbitrary cloak can be obtained. In a different way, Ma et al.  propose a numerical method by using position surfaces and tracing lines to represent the spatial compression during a linear transformation along the tracing lines. However, the reported works all rely on semi-analytic and semi-numerical approaches to evaluate the transformation matrix A for an arbitrary cloak. So a general and flexible method for designing an arbitrary cloak is still necessary. Especially, the answer to the following question needs to be clarified: is there a governing equation that allows us to evaluate the transformed material parameters for cloaks with arbitrary shapes in a systematic way?
In this paper, we will use the Laplace’s equation to determine the transformation matrix in a unified manner. It is shown that the function of the transformation media corresponds in fact to specific boundary conditions. So the transformed material parameters are determined by the solutions of the Laplace’s equation under proper boundary conditions whatever the shape of a cloak. The paper will be arranged as follows: The design method will be explained in detail in Sec. 2. For regular cloaks, the expressions of the permittivity and permeability tensors can be analytically derived with the proposed method. Two examples are given in Sec. 3. For irregular cloaks, the proposed method is easily used to compute their material parameters. The full-wave simulations based on the Maxwell’s equations validate the effect of the cloak, which will be discussed in Sec. 4. The discussion and conclusion will be given in Sec. 5.
2. Design method
According to the coordinate transformation method [2, 8, 18], under a space transformation from a flat space x to a distorted space x′(x), the permittivity ε′ and permeability μ′ in the transformed space are given by 
where ε and μ are the permittivity and permeability of the original space, respectively. A is the Jacobian transformation tensor with the components Aij = ∂x′i/∂xj, which characterizes the geometrical variations between the original space Ω and the transformed space Ω′. The determination of the matrix A is the crucial point for designing transformation mediums.
In the continuum mechanics , the tensor A is called the deformation gradient tensor for an infinitesimal element d x deformed to d x′ under the space transformation. The transformation can be decomposed into a pure stretch deformation (described by a positive definite symmetric tensor V ) and a rigid body rotation (described by a proper orthogonal tensor R ), so the tensor A can be expressed as A = VR . Suppose that material parameters in the original space are homogeneous and isotropic, and they are expressed by scalar parameters ε 0 and μ 0. Consider the left Cauchy-Green deformation tensor B = V2 = AA T, Eq. (1) can be rewritten as
In the principal system, the tensor B can be expressed in the diagonal form
Equation (4) show that the material parameters of the transformation material can be calculated from the principal stretches of the deformation induced by the space transformation. Note that the transformed material parameter and the deformation tensor have the same principal directions and the rigid body rotation in deformations has no contributions on the material parameters. Thus, the calculation of material parameters is converted to evaluate the deformation field for the flat grid element distorted by the space transformation.
Figure 1 shows the scheme for constructing an arbitrary cloak. Suppose an original region enclosed by the outer boundary represented by b, inside of this region, we define a point denoted by a. An arbitrary cloak can be constructed by enlarging the point a to the inner boundary a′, while keeping the outer boundary of the region fixed (b= b′). This condition can be expressed by U′(a) = a′ and U′(b) = b′, where the operator U′ is the new coordinates for a given point during the transformation. Now the problem is how to determine the deformation field ∂U′ /∂ x within the cloak layer enclosed by the inner and outer boundary a′ and b′ for a specific transformation. The commonly used operator U′ for designing a cloak is a linear transformation [2, 9-13, 15-17]. For example, the radial displacement of a spherical cloak is assumed to be a linear relation r′ = (b′−a′)r/b′ + a′. However, for an arbitrary cloak, it is very difficult to express analytically the boundaries a′ and b′, the calculation of deformation field is usually very complicated.
To ensure a cloak without reflections, the internal deformation field of the cloak layer must be continuous. The deformation tensor is calculated by the partial derivative of displacement with respect to the original coordinate, so the displacement fields must be smooth enough. It is known that the Laplace’s equations with Dirichlet boundary conditions will always give rise to harmonic solutions . This suggests that the displacement field inside of the cloak layer can be calculated by solving Laplace’s equations ∆U′ = 0 with the boundary conditions U′(a) = a′ and U′(b) = b′. In order to keep from the singular solution of the Laplace’s equations, we can use the inverse form of the Laplace’s equations as
where Ui denotes the original coordinates in the original space. The corresponding Dirichlet boundary conditions then become U(a′) = a and U(b′) = b. After solving Eq. (5) with proper boundary conditions, we can get the deformation field inside of the cloak layer that characterizes the distortion of the element. The transformed material parameters can then be calculated from Eq. (4). This method is flexible whatever the shape of cloaks. In the next two sections, we will explain in detail how to obtain analytical expressions of material parameters for regular cloaks and how to design irregular cloaks by solving Eq. (5) numerically.
3. Application to regular cloaks
3.1 Spherical cloak
The Laplace’s equation in the spherical coordinate system is expressed as
For a spherical cloak, we let θ′ = θ and φ′ = φ. Then the radial coordinate satisfies the following equation
With the boundary condition r(r′ = a′) = 0 and r(r′ = b′) = b′ (b′ = b ), the solution of Eq. (7) is given by
where a′ and b′ are respectively the radii of the inner and outer boundaries of the cloak. Equation (9) indicates a nonlinear transformation, which is different from the linear one r′ = (b′−a′)r/b′+a′ commonly used for a spherical cloak . The principal stretches corresponding to the nonlinear transformation (9) are given by
From Eq. (4), the material parameters of the spherical cloak are
Figure 2 shows the electric field distribution of a spherical cloak with the material parameters (11), illuminated by a plane electromagnetic wave. The results validate the invisibility of the designed cloak.
3.2 Elliptical cylindrical cloak
In this section, we will design a 2D elliptical cloak by solving the Laplace’s equation (5). The Cartesian coordinates (x,y,z) can be expressed by the elliptic cylindrical coordinates (ξ, η, z) as
where c is the focal length of the elliptic cylinder. The scale factors of the elliptical cylindrical coordinate system are 
We can construct an elliptic cloak by setting η′= η and z′ = z, then Eq. (14) becomes ∂ 2 ξ/∂ ξ′2 = 0 . With the boundary conditions ξ(ξ′ = a′) = 1 and ξ(ξ′ = b′)= b (b′ = b), we can obtain a linear transformation relation ξ′ = (b′ − a′)(ξ-1)/(b′−1)+ a′, where a′ and b′ are respectively the coordinates of the inner and outer boundaries of the elliptic cloak. The principal stretches for this linear transformation are given by
According to Eq. (4), the material parameters of the elliptic cylindrical cloak are given by
where we have assumed ε 0 = μ 0 = 1. It can be found that Eq. (16) agrees with the results given by Ma et al  from a different method. In fact, the solution of Laplace’s equation gives a certain spatial mapping of transformed coordinates. It is expected that the same result can be found if the obtained mapping is the same as that pre-assumed by the other methods.
4. Application to irregular cloaks
For an arbitrary cloak, its boundary is difficult to be expressed in an analytical form, so the Laplace’s equation must be solved numerically. In the following, we will show how to use the proposed method to design a cloak of irregular shape with help of the commercial software COMSOL Multiphysics. Firstly, we depict two boundaries a′ and b′ for an arbitrary cloak and set the boundary conditions U(a′) = 0 and U(b′)=b (b′=b), then solve the Laplace’s equations to obtain the deformation field ∂U/∂ x′ inside of the cloak layer. The material parameters of the cloak can be calculated numerically by the principle stretches λ 1, λ 2, and λ 3 (evaluated from the deformation field) according to Eq. (4). This can be achieved with the PDE solver provided by the commercial software COMSOL Multiphysics. In order to check the performance of the designed cloak, the cloak device will be illuminated by a plane electromagnetic wave with help of the RF module in the same software. So the design and validation of the cloak can be well integrated in a two-step modeling.
As an example, we define a two-dimensional arbitrary cloak, which boundaries are formed by the circles of radii 1cm and 3cm, as shown in Fig. 3. The cloak is embedded in the air of a square area surrounded by the PML regions. In the PDE modes of Laplace’s equation, we set the boundary conditions U(a′) = 0 and U(b′) = b (b′ = b) for the inner and outer boundaries, then solve the equation ∆′U = 0 to get the coordinates x = U(x′) [x′ ∊ (a′,b′) ], from which we obtain x′ = U′(x) [x∊(0,b)]. Figure 4 gives the coordinate lines of x′ transformed from a flat space x. The line path implies that electromagnetic waves will propagate around the area bounded by a′ and any objects put inside of this area will be invisible. For the verification, we let a plane electromagnetic wave of frequency 20 GHz with the electric field polarized along the z direction illuminate the cloak. The full fields of the system are solved with help of the TE waves mode of COMSOL Multiphysics. Fig. 5 (a) and 5(b) shows the contour plots of the electric field Ez for the waves incident on the cloak horizontally and at an angle of 45° from the x direction, respectively. It can be seen from the figures that the constructed arbitrary cloak doesn’t disturb the incident waves and can shield an irregular obstacle from detection. For the physical realization, one can easily retrieve the material parameters of the arbitrary cloak from the deformation field according to Eq. (4). The values of some components of the permittivity and permeability tensors, εzz, μxx, μxy, and μyy of the cloak are shown in Figs. 6(a), 6(b), 6(c), and 6(d), respectively. It can be found that the cloak are highly anisotropic and must be realized with the metamaterial technology.
The proposed method can also be used to design arbitrary cloaks with many separated embedded obstacles. For example, we use the same system as shown in Fig. 3, but separate the inner circles of the radius 1cm apart, as shown in Fig. 7. The coordinates of central points of the four circles in centimeters are O 1(0,1), O 2(1,0), O 3(0,-1), and O 4(-1,0). The boundary conditions for numerically solving the Laplace’s equation then become U(ai′) = Oi(i =1, 2,3,4) and U(b′) = b. Figures 8 (a) and 8(b) show the contour plots of the electric fields Ez for the waves incident on this cloak horizontally and at an angle of 45° from the x direction, respectively. We can see from the figures that the perfect invisibility is still achieved even if the embedded obstacles are separated apart.
Finally the proposed method is applied to design the cloaks with the boundaries drawn randomly or imported from other tool packages, such as Microsoft Office Visio by DXF files. In the following example, we draw an arbitrary cloak at the user-guide interface in the COMSOL software. After solving the Laplace’s equation with the boundary conditions for the cloak, we calculate the distorted coordinate lines, as shown in Fig. 9. Figures 10(a) and 10(b) show the contour plots of the electric fields Ez for TE electromagnetic waves impinging the cloak horizontally and vertically respectively. It can be seen that the designed cloak does not disturb the outside fields and shields the inner random object from detection.
It must be mentioned that there are small scatterings in the external region of the cloaks shown in Figs. 5, 8, and 10. This imperfection is believed to be induced by the numerical method due to discretization; this problem can be improved by fine mesh with cost of computation time.
5. Discussion and conclusion
The proposed method starts from a new understanding of the transformation tensor A from the mechanics point of view. Through a rigorous derivation, the material parameters of cloaks are related to the deformations attached to the coordinate grids. The design of a cloak is completely converted to the calculation of deformation fields; we then propose to use Laplace’s equation to evaluate this deformation field whatever the shape of the cloak. In this sense, the proposed method is systematic and built completely on the solution of Laplace’s equation with proper boundary condition without any assumption on the transformation, especially for irregular cloaks. Different from other methods where the transformation tensor A is obtained explicitly by a given explicit coordinate transformation, the proposed method derives the tensor A implicitly by solving Laplace’s equation. The material parameter of a cloak with given shapes is not unique and depends on the coordinate transformation we take. The geometric nature of the (linear or nonlinear) mapping implied by the Laplace’s equation still needs further investigation.
In conclusion, based on coordinate transformation method, a systematic method is proposed to evaluate the transformed material parameter for arbitrary cloaks. In this method, the calculation of the material parameters for a cloak becomes the computation of deformation field governed by the Laplace’s equation with proper boundary conditions. The principle stretches induced by the transformation are related to the permittivity and permeability tensors in the principle system. By imposing different boundary conditions, we can design arbitrary cloaks with single or many separated embedded obstacles. For regular cloaks with spherical, cylindrical and elliptical shapes, explicit expressions of the material parameters can be easily obtained by solving the Laplace’s equation analytically. The design and full-wave simulation of a cloak can be well integrated in a two-step modeling. It is worth to note that the proposed method is not just limited to the design of cloaks; it can also be extended to design other transformation materials, such as rotators , concentrators , beam shifter  and beam bender  by the modification of boundary conditions. In addition, we can use this method to design acoustic transformation materials. These issues are currently under investigation.
Dr. Xiaoning Liu is acknowledged for helpful discussions. This work is supported by the National Natural Science Foundation of China (90605001, 10702006, 10832002), and the National Basic Research Program of China (2006CB601204).
References and links
1. A. Greenleaf, M. Lassas, and G. Uhlmann, “On non-uniqueness for Calderón’s inverse problem,” Math. Res.Lett. 10, 685–693 (2003).
4. 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]
5. W. Cai, U. K. Chettiar, A. V. Kildishev, and V. M. Shalaev, “Optical cloaking with metamaterials,” Nat. Photonics 1, 224–227 (2007). [CrossRef]
7. U. Leonhardt, “Notes on conformal invisibility devices,” New J. Phys. 8, 118 (2006). [CrossRef]
8. U. Leonhardt, “General relativity in electrical engineering,” New J. Phys. 8, 247 (2006). [CrossRef]
9. M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, and D. R. Smith, “Design of Electromagnetic Cloaks and Concentrators Using Form-Invariant Coordinate Transformations of Maxwell’s Equations,” Photon. Nanostruct. Fundam. Appl. 6, 87–95 (2008). [CrossRef]
10. H. Ma, S. B. Qu, Z. Xu, J. Q. Zhang, B. W. Chen, and J. F. Wang, “Material parameter equation for elliptical cylindrical cloaks,” Phys. Rev. A 77, 013825 (2008). [CrossRef]
12. D. Kwon and D. H. Werner, “Two-dimensional eccentric elliptic electromagnetic cloaks,” Appl. Phys. Lett. 92, 013505 (2008). [CrossRef]
13. W. X. Jiang, T. J. Cui, G. X. Yu, X. Q. Lin, Q. Cheng, and J. Y. Chin, “Arbitrarily elliptical-cylindrical invisible cloaking,” J. Phys. D: Appl. Phys. 41, 085504 (2008). [CrossRef]
14. W. Yan, M. Yan, Z. Ruan, and M. Qiu, “Coordinate transformations make perfect invisibility cloaks with arbitrary shape,” New J. Phys. 10, 043040 (2006). [CrossRef]
15. W. X. Jiang, J. Y. Chin, Z. Li, Q. Cheng, R. Liu, and T. J. Cui, “Analytical design of conformally invisible cloaks for arbitrarily shaped objects,” Phys. Rev. E 77, 066607 (2008). [CrossRef]
17. H. Ma, S. Qu, Z. Xu, and J. Wang, “Numerical method for designing approximate cloaks with arbitrary shapes,” Phys. Rev. E 78, 036608(2008). [CrossRef]
18. G. W. Milton, M. Briane, and J. R. Willis, “On cloaking for elasticity and physical equations with a transformation invariant form,” New J. Phys. 8, 248(2006). [CrossRef]
19. W. M. Lai, D. Rubin, and E. Krempl, Introduction to Continuum Mechanics, 3 edition. (Butterworth-Heinemann, Burlington, 1995).
20. R. Courant and D. Hilbert, Methods of Mathematical Physics, Vol.2, 1 edition. (Wiley-Interscience ,New York, 1989).
21. G. Arfken, Mathematical Methods for Physicists (Academic Press, Orlando, 1970).
22. H. Chen and C. T. Chan, “Transformation media that rotate electromagnetic fields,” Appl. Phys. Lett. 90, 241105 (2007). [CrossRef]
23. M. Rahm, S. A. Cummer, D. Schurig, J. B. Pendry, and D. R. Smith, “Optical Design of Reflectionless Complex Media by Finite Embedded Coordinate Transformations,” Phys. Rev. Lett. 100, 063903 (2008) [CrossRef] [PubMed]