Spacial Energy Distribution Manipulation with Multi-focus Huygens Metamirror

Huygens metasurface is a planar array of crossed electric and magnetic dipoles, which provide specific surface current to tailor the electromagnetic field distribution. By changing the geometrical parameters of the proposed unit cell, the manipulation range of reflection phase can achieve 2π, while the amplitude of the reflection coefficient can keep above 0.993. Based on the designed Huygens meta-atoms, a novel multi-focus Huygens metamirror is proposed at microwave range in this paper. Utilizing the meta-atoms with the desired reflection phase distribution as calculated, the incident plane wave can be converged to designated points in any desired fashion including focal number, location and intensity distribution, which exhibits outstanding manipulation capability. Our research on Huygens metamirror provides a fascinating design of multi-focus imaging in microwave region, which makes it potential applications in antenna and imaging systems.

indirectly according to the boundary condition. Surface electric and magnetic impedance are defined here to build the relationship between the surface current and the electromagnetic field. It can be modified by changing the structure of meta-atom appropriately. Then, based on the scattered field distribution, the relationship between the reflection and transmission coefficient and the surface impendence can be built. Overall, according to the equivalence principle and boundary conditions, arbitrary reflection and transmission coefficient can be realized with specific surface electromagnetic impedance 25,26 .
m where = φ R re j r and T te j t = φ are the predetermined reflection and transmission coefficients. Z e and Z m denotes surface electric and magnetic impedance respectively. For the design of reflective Huygens metasurface, the reflection coefficient should remain to be 1 while the transmission coefficient equals 0. Therefore, when r = 1 and t = 0, equations (1) and (2) where ϕ r represents the reflection phase, which changes from 0 to 2π. η is set to be 377Ω as the impedance of the free space. From the above two expressions, Z e and Z m are pure imaginary, so the reflection phase can be controlled only by a lossless metasurface theoretically. The desired electric and magnetic impendence can be achieved by the proposed meta-atom, as shown in Fig. 1. Two mirror-symmetry metal split rings are placed on the top layer of the substrate, operated as the magnetic dipoles. By tuning the arm length of the magnetic dipoles, marked as l w and l g , the magnetic response can be manipulated. Meanwhile, the electric dipole is composed of two mirror-symmetry T-shaped metal stripes and the induced electric current can be modified by changing the capacitor width w e . The magnetic and electric current should be mutually perpendicular, so the generated electromagnetic field can be manipulated individually without interference. The proposed unit cell is simulated using the commercial software, CST Microwave Studio, to analysis the feature of the proposed meta-atom. Meanwhile, the incidence is set as normal TM microwave to ensure the electric field parallel to electric dipole. According to simulation results, it works at 10.9 GHz. Figure 1(b) shows the designed meta-atom, which composed of metal stripes integrated on dielectric substrate. Except the variable parameters mentioned above, the other parameters remain to be constant, as d = 2.45 mm, c = 0.1 mm, d e = 0.9 mm, t = 0.3 mm. Compared with previous researches, the proposed Huygens meta-atom composed of double mirror-symmetric magnetic dipoles, instead of the single one. This mirroring effect can support toroidal dipole excitation. As shown in Fig. 2(a), the simulated surface current along two metal split rings circles in opposite directions. Therefore, two induced magnetic dipole (m → ) are head-to-tail in a small circular region, which generates the toroidal dipole 27 . There are two merits of this modification. First, since the quality factor of toroidal resonance is higher than the magnetic resonance, full manipulation of the reflection phase can be achieved by changing the arm length, as shown in Fig. 2(b). Therefore, there is no need to change the shape of metal patch significantly, which simplifies the design of meta-atom. Second, the mirror-symmetric structure guarantees the direction of magnetic dipole opposite with neighboring one. This design can dramatically reduce the radiation loss and improve the efficiency 28 . Figure 1(a) illustrates the simulation results of the proposed Huygens meta-atom. The range of reflection phase can cover the whole range of 360 degrees and the step of phase shift is set to be 30 degrees. Actually, arbitrary phase shift can be realized with accurate geometric design. Meanwhile, the amplitude of the reflection coefficient keeps above 0.993 and the amplitude of the transmission coefficient keeps below 0.154, which means the energy of the incidence can be reflected and controlled totally by our design.

Multi-focus Metamirrors
Multi-focus imaging can be achieved by modulating phase distribution of the metasurface. By introducing desired phase changes to the reflected wave utilizing the proposed meta-atom, the incident wave can be focused to specified positions. One of the fundamental methods of designing hologram is setting virtual point sources at the preset focal locations 29,30 . Then the electric field generated by each virtual source can be superimposed to obtain the virtual electric field distribution. Here the electromagnetic field propagation from each virtual point source is described by Green function. Then the virtual electric field is mapped onto the metasurface and the desired phase delay at the positions of different meta-atoms ϕ(x, y, 0) can be calculated. To perform the desired focusing image, the reflection phase of each meta-atom should be set as −ϕ(x, y, 0) under the incident plane wave 24 , since the calculation of phase conjugation in frequency domain is equivalent to the time reversal for the incident monochromatic wave 31,32 . Accordingly, the reconstructed electric field is converged to the positions where virtual sources locate. However, the electric field intensity decreases linearly with the reciprocal of the distance between meta-atoms and virtual sources. The calculation of phase conjugation only ensures the incident wave focus at preset locations, but can't eliminate the difference of intensity between the focal points. Therefore, the method mentioned above is only applicable to the situation of single focus. For multiple-focus imaging, the intensity of virtual point sources should be modified with different weight based on the specific focal intensity distribution. Hence, the assigned reflection phase of each meta-atom is modified as = − + − + is the distance between the j th meta-atom located at (x j , y j , 0)(j = 1 to M) and the i th focal point at (x i , y i , z i )(i = 1 to N). k is the phase constant. w i is the weight factor of i th virtual source, which presents the intensity ratio of the i th virtual source to the first one. Particularly, w 1 is set to 1. Then, the reflection phase is determined to be −ϕ(x j , y j , 0) for the meta-atom located at (x j , y j , 0). Next, based on the established reflection phase distribution, the reconstructed electric field distribution can be calculated through the superposition of the electric field component generated by each meta-atom. Since the electric field component generated by each meta-atom can be described by the Green function theoretically, the reconstructed electric field distribution under the incident plane wave is described as = − + − + denotes the distance between the meta-atom at (x j , y j , 0) and the observation point at (x, y, z). Here the observation point is set at (x i , y i , z i ) to get the intensity of i th reconstructed focal point, |E(x i , y i , z i )|, denoted by the amplitude of the electric field at the position of the i th preset focal point.
Then, the relationship between w i and the intensity distribution of focal points is analyzed to achieve the desired imaging. To describe the intensity distribution of focal points, the intensity ratio of the i th focal point to the first one is defined as k i . Based on the equation (5), w i is the weight factor of intensity of the i th virtual focal point, which can directly affect the intensity distribution of reconstructed electric field. Particularly, w i doesn't equal to k i here, since the reflection coefficient of each meta-atom on the surface keeps the same, and the only variable parameter is the reflection phase. Therefore, the virtual electric field differs from the reconstructed one slightly, which means it is not a completely reversible process. Hence, w i is modified to guarantee the reconstructed electric field at focal points equal to the desired intensity ratio (k i ). Accordingly, the formulation is described as For application of the proposed algorithm, we design single-focus and double-focus metamirrors with specific energy intensity distribution respectively. For the single-focus metamirror, the first step is to set a virtual source at the position of focal point to get the phase delay (φ 1 , φ 2 …, φ M ) on the metasurface, as shown in Fig. 3(a). Next, as Fig. 3(b) illustrates, each meta-atom is assigned negative reflection phase delay (−φ 1 , −φ 2 …, −φ M ) to ensure the electric field emitted from each meta-atom superimposes in phase at the position of foal point.
When the number of focal points exceeds one, using our proposed algorithm, specific focal intensity ratio of the focal points can be realized. As shown in Fig. 3(c), the virtual electric field generated from double virtual sources with the intensity ratio (w i ) is superposed to calculate the phase delay onto each meta-atom (φ 1 , φ 2 …, φ M ). Similarly, the reflection phase distribution on the metasurface is set to the negative phase delay (−φ 1 , −φ 2 …, −φ M ), as shown in Fig. 3(d). This operation is to divert the energy to the desired focal points. Furthermore, by determining the value of w 2 according to the equation (7), the intensity ratio of the reconstructed electric field at two focal points can equal to desired k 2 . Therefore, the manipulation of the intensity distribution between the focal points can also be achieved utilizing the phase change of each meta-atom.
To demonstrate our theoretical analysis, two types of metamirrors have been simulated with the commercial simulation software FDTD Solutions respectively. Here meta-atoms change along the x-axis to achieve the desired reflection phase distribution, while the meta-atoms keep the same along the y-axis. In our simulations, the boundary condition along y-direction is set to be the perfect electric conductor (PEC) boundary. Hence, the focusing performance remains the same along the y-axis, and xoz plane reflects the focusing features of the designed metamirror, as shown in Fig. 3(e).
Firstly, we design single-focus metamirrors to converge the incident plane wave to desired focal points in xoz plane. Figure 4 shows the electric field distribution in xoz plane and the required spatial phase of each meta-atom. The simulation results show the single focal point can be moved along x and z direction at will. It demonstrates that the incidence can be converged to any single point in xoz plane with calculated phase distribution described by equation (5).
Different from the design of single-focus metamirrors, the weight factor needs to be taken into consider for double-focus imaging. Figure 5 illustrates the design of two focal points with the same intensity placed along x-axis and z-axis respectively. To achieve the uniform intensity ratio, the weight factor is modified to compensate the decrease of electric field intensity. As shown in Fig. 5, the simulation results of the electric field distribution agree very well with the theoretical ones. Furthermore, the corresponding phase distribution of the metasurface is shown in the right column. Besides the uniform intensity ratio, by changing the weight factor, we can also achieve arbitrary focal intensity distribution. Figure 6 illustrates the design of two focal points ranged along z-axis with different focal intensity ratio. In both situations, the incident plane wave is converged to the desired focal location, and the simulation results of focal intensity ratio are in good agreement with the theoretical ones. Hence, it indicates that the weight factor can manipulate the intensity ratio of each focal point respectively, adding a new degree of freedom to the manipulation of electromagnetic wave.

Discussion
In conclusion, lossless reflective Huygens metasurface composed of a series of crossed electric and magnetic dipoles are designed and simulated in this paper. The reflection phase of the proposed meta-atoms can be manipulated in the whole range of 2π with high efficiency. Based on them, multi-focus Huygens metamirrors have been first proposed. Arbitrary focus imaging can be achieved at microwave range in consideration of the factors of focal number, location and intensity ratio. To demonstrate it, the intensity distribution of the electric field on the focal plane is simulated, which agrees very well with the theoretical analysis. The proposed multi-focus Huygens metamirrors can be applied to many fields, which makes our design a potential application in antenna and imaging systems. Furthermore, the design method of reflective Huygens metasurface with the function of multi-focusing proposed in this paper can be applied to other frequency ranges.