Electromagnetic time-harmonic and static field polygonal rotator with homogeneous materials

We propose a scheme of designing polygonal rotator with homogenous materials by using linear coordinate transformation. Our strategy is available for both time-harmonic electromagnetic field case and static field case. In particular, we found that only one anisotropic material is needed in static field case, and the density of field in the central region can be altered to be denser or sparser, or stay the same. The magnetostatic field rotator can be realized by multilayered structure composed of ferromagnetic materials and superconductor, and the direct current rotator can be realized by metals with different conductivity. Numerical results verify the effectiveness of our strategy in both time-harmonic field case and static field case.

We propose a scheme of designing polygonal rotator with homogenous materials by using linear coordinate transformation. Our strategy is available for both time-harmonic electromagnetic field case and static field case. In particular, we found that only one anisotropic material is needed in static field case, and the density of field in the central region can be altered to be denser or sparser, or stay the same. The magnetostatic field rotator can be realized by multilayered structure composed of ferromagnetic materials and superconductor, and the direct current rotator can be realized by metals with different conductivity. Numerical results verify the effectiveness of our strategy in both timeharmonic field case and static field case. Pendry et al. proposed the transformation optics method to control electromagnetic field 1 . This new method pave a new way to design electromagnetic devices such as invisibility cloaks [1][2][3][4][5][6][7][8][9][10][11] , field concentrators 12 , electromagnetic wormholes 13 , optic black holes 14,15 , bending wave guides [16][17][18] , metalens [19][20][21] , field rotators 4,[22][23][24][25] , and so on. The field rotator was first proposed by Chen et al. 4 in 2007. This rotator can be divided into two parts, rotator shell and inner cylinder. The inner cylinder rotates a certain angle around the symmetry axis. When the radius approaches external radius of the rotator shell, the rotational angle is reduced to zero. However, the electromagnetic parameters of this kind of field rotator are inhomogeneous and anisotropic, which are difficult to realize. A rotator with reduced parameters were experimentally proposed to simplify the design and fabrication 22 , but imperfect rotation field performance was introduced due to impedance mismatch. Inspired by Chen's work, versions of heat fluxes were proposed 23,25 , which are composed of inhomogeneous and anisotropic conductivity materials. As the parameters in the rotators involve inhomogeneous parameters, they are not easy to implement.
One approach to simplify the constructive parameters of rotator is the linear homogeneous coordinate transformation, which can eliminate the inhomogeneity of parameters. The linear homogeneous coordinate transformation was proposed by Xi et al. 26 . By transforming a large triangle to a smaller one, a homogeneous one-dimensional cloak was designed. This approach was later extended to omnidirectional polygonal cloaks 27 . The whole polygonal cloak is divided into several triangles, and in each segment a linear homogenous coordinate transformation is applied, and the cloaked region is transformed from a large area to a smaller one. A twofold spatial compression method was also proposed to design a polygonal cloak by expanding several lines or a line into a concealed region 28,29 . The linear homogenous coordinate transformation is further applied to concentrator 30 , waveguides 31 and plasmonics 32,33 to overcome the spatial variation.
In this paper, we use the linear coordinate transformation to design an electromagnetic field rotator which is composed of homogenous materials. We extend this strategy to static field case which can be implemented by realizable materials. We show that two kinds of anisotropic materials are needed to design a time-harmonic electromagnetic field rotator. While in static field case, we show that only one anisotropic material is needed. This anisotropic material can be realized with multilayered structure composed of two isotropic materials based on the effective medium theory. Moreover, for static field cases, the field in the central region can be rotated and at the same time be denser, sparser, or unchanged.
Electromagnetic field rotator. We firstly consider the linear transformations in two-dimensional space from a triangle in virtual space with local coordinate axes (x, y) to another triangle in physical space with local coordinate axes (x′, y′) ( Fig. 1 Substituting coordinate of the three points in (1), we can get where J is the Jacobian matrix. By using transformation optics method, we can obtain the electromagnetic parameters of the triangle in physical space, which are Our strategy is dividing a polygon into several triangle segments, and then applying a linear coordinate transformation in each triangle segment. We take the quadrate rotator shown in Fig. 2 as an example to demonstrate our scheme. The coordinate system is built with the origin at the center of rotator. The rotator shell can be divided into 8 segments, which can be grouped into two types: segment I and segment II marked in blue and red, respectively. Each triangle segment in virtual space is linearly transformed to its corresponding segment in physical space. The central region is only rotated by certain degree (θ) without changing its sizes, which means that the constructive parameters of this area is unchanged after coordinate transformation. Compared with ref. 24 , the transformation function in our paper is different, as a result, the rotation angle is very flexible. From Eq. (2), with a linear transformation, the Jacob matrix will be independent of the position (x, y), i.e., it will be a constant  www.nature.com/scientificreports www.nature.com/scientificreports/ matrix. Therefore, from Eqs (3) and (4), one can see that the constitute parameters of the rotators in all regions will be independent of the position and become homogeneous. This facilitate the practical realizations.
Let's consider the TM wave case and that all segments in virtual space are vacuum. According to the procedure mentioned above, we can get the parameters of each segment, which are , for segment I, , for segment II.
In the following, we make full wave simulations by employing the commercial finite element method software, COMSOL Multiphysics, to demonstrate the performance of the rotator. We schematically illustrate the simulation in Fig. 3(a) The electromagnetic field in central region is rotated anticlockwise by 45°. Figure 3(a,b) show the magnetic field distribution in and near rotator, when an Hz polarized plane wave is incident along x direction and +x y 2 2 2 2 direction, respectively. It is obvious that rotator shell rotates the field in central region 45° anticlockwise, and the field out of the rotator is undistorted. Note that when the rotational angle is larger than the critical angle, for example, −69.3° in our case, both permeability and permittivity become negative, which increase the difficulty of practical implementation. However, this limitation can be easily circumvented by adding another rotator shell out or in the rotator, as shown in Fig. 3(c,d), in which there is a second smaller rotator in the central region and field in the inner core is rotated anticlockwise by 90°.
Magnetostatic field rotator. We go a step further and extend the scheme of rotator to magnetostatic field case. The difference between static field and time-harmonic electromagnetic field is that the magnetic field and electric field are decoupled in static field. Thus we only consider the permeability of materials in magnetostatic field, which provides more degree of freedom in designing different kinds of devices and will simplify the parameters of materials. The scheme of magnetostatic field rotator is same as that of electromagnetic rotator illustrated in Fig. 2. By applying the equations in (5) and (6), we can get the parameters of the rotator, , for segment I, 1 2 and according to Eq.(4), the eigenvalue of µ ′ 1 and µ ′ 2 will be the same, which means that the segment I and II can be realized by the same material. As an illusion, the geometries of magnetostatic field rotator are similar to that of electromagnetic field rotator in Fig. 3(a) . This material can be realized by a multilayered structure composed of two isotropic materials. µ = .
0 2903 a can be achieved by using ferromagnetic materials and µ = .
3 4445 b can be achieved by using dc metamaterials 34,35 . Thus our strategy manifests that only one anisotropic material is used to design a perfect magnetostatic field rotator. Besides, if the size of central square in physical space is smaller or larger than that in virtual space, the magnetic field in this area will be denser or sparser. We can still use the same anisotropic material mentioned above to design a concentrator or a sparse device.
We firstly set = . ′ = ′ = + a x y t x y 0 5, 0, 2 (if let t > 1, then we obtain a concentrator; or if let < < t 0 1 , then we obtain a sparse device. To get an exact concentrator, here we set t = 2), and set x 1 and y 1 as unknowns. If let , we get the following parameters of the concentrator: It is interesting to see that we only use one anisotropic material to design a rotator, a concentrator and a sparse device.
We use COMSOL Multiphysics to demonstrate the effectiveness of our strategy. In the simulation, we set the left boundary as magnetic potential with = A 1 Wb/m z , right boundary as magnetic potential with = A 0 Wb/m z , top and bottom boundaries as periodic conditions. To verify the omnidirectional performance of the devices, we also show the cases when the magnetic field source is rotated anticlockwise by 45°. Figure 4(a,d) show the contour of magnetic potential in and near rotator. One can see that the contour of the magnetic potential in the central region is rotated anticlockwise 39.85°, and that out of the rotator is undistorted. Figure 4(b,e) show the contour of magnetic potential in and near concentrator. It is obvious that the contour of magnetic potential in the central region is rotated anticlockwise 27.20°, and the density of contour of magnetic potential is double that out of the concentrator. Moreover, the contour out of the concentrator is undistorted. Figure 4(c,f) show the contour of magnetic potential in and near sparse device. One can see that the contour of magnetic potential in the central region is rotated anticlockwise 24.22°, and the density of the contour of magnetic potential is half that out of the sparse device. Moreover, the contour out of the sparse device is undistorted. Note that all these devices are composed of only one anisotropic material, which can be realized by the multilayered structure consisting of ferromagnetic materials and superconductor.
This strategy can also be applied to electrostatic field case and dc case. In dc case, we can use metals with different conductivity to design these devices mentioned above, for example, stainless iron, copper, and iron, which are used by Han et al. to design a cloak and a concentrator 36 .

conclusion
In this paper, a rotator consisting of homogeneous anisotropic materials is proposed, by applying the linear coordinate transformation. Our strategy is available in both time-harmonic electromagnetic field and static field. We further demonstrate that, for static field, only one anisotropic material is needed to design a rotator. In its central region, the density of the field can be altered to be denser or sparser, or stay the same. We show that the rotator for magnetostaic field can be realized by ferromagnetic materials and superconductor. For direct current, the rotator can be easily realized by metals with different conductivity, such as stainless iron, direction, respectively. All these devices above are composed of only one anisotropic material.