A Low-Cost Multibeam Switching Antenna Using Reconfigurable Hybrid Metasurface for Beamforming Applications

In this paper, we proposed a multibeam switching antenna based on a low-cost reconfigurable hybrid metasurface applied for beamforming systems. The antenna consists of two parts: a microstrip feed antenna and a transmission hybrid metasurface. The latter is composed of three types of units with different amplitude and phase responses to electromagnetic waves so as to control the beams of the feed antenna. Sixteen PIN diodes are arranged in the metasurface with a simple bias network. When two different direct-current voltages are applied to the PIN diodes, the antenna can dynamically switch between two beams and four beams. For demonstration, the proposed antenna is fabricated, and the measured results show that the antenna operates at 9.07–9.42 GHz (−10 dB bandwidth) with a total size of 1.80λ0 × 1.52λ0 × 0.22λ0 (λ0 corresponds to the wavelength of 9.28 GHz in free space). With the merits of a compact structure, low cost and good radiation performance, the proposed design is suitable for beamforming applications.


Introduction
The multibeam antenna can accurately point to multiple users in time and has the advantages of multi-dimensional coverage and a fast communication rate, so it is of great significance in beamforming applications. With the rapid development of the low-orbit satellite communication industry and 5G wireless communication, the Internet-of-Space-Things (IoST) system [1], mainly composed of ground stations, aerial vehicles and small satellites, has attracted more and more interest. The antennas used in the IoST system should provide the best possible performance and maintain the minimum size. The dynamical multibeam switching antenna (MSA) has become a potential candidate for IoST systems due to its ability to cover multiple areas, improve multipath effect and support high data rates with a compact structure [2]. The flat MSA can also be widely applied in and beyond 5G and 6G mobile communication systems and multi-target radar systems [3].
Generally, the categories of MSA mainly include phased array antenna [2,4], multiport antenna [5,6] and coding metasurface antenna [7,8]. As mentioned in [9], the phased array antennas can achieve flexible beam switching, but the configuration is complicated with a complex structure, and they suffer high costs. Compared with digital beamforming methods, MSAs based on metasurfaces have lower energy consumption and complexity. As discussed in [10], four ports in each unit are used to generate various radiation beams. However, when a large-scale antenna array is composed of these units, the feed network of the system is quite complex. In recent years, with the introduction of digital coding metasurfaces [11], there are more choices for the design of an MSA. Active devices, such 2 of 12 as a PIN diode or varactor, are added to the coding metasurface to discretize the phase or amplitude response of the electromagnetic wave so as to obtain flexible beam control capability. For example, a 2-bit programmable metasurface in [12] is adopted to realize flexible multibeam switching. However, these reflective or transmissive metasurfaces require additional horns to irradiate, resulting in a bulky system volume. In addition, as discussed in [13], a reconfigurable metasurface loaded with PIN diodes is designed to achieve beam switching. However, this method uses a large number of PIN diodes, which undoubtedly increases the loss and cost. Of course, the method of mechanically rotating the coding metasurface can also achieve multibeam switching [14]. However, it lacks real-time capability and flexibility with heavy weights.
The hybrid metasurface is composed of units with different amplitudes and phase responses to incident waves, which makes the manipulation of electromagnetic waves more flexible. For compact and low-cost multibeam antennas based on metasurfaces, hybrid metasurfaces are a promising choice. As presented in [15], a hybrid metasurface with a split-ring resonator and U-shaped unit cell is adopted to split a beam. A hybrid digital coding metasurface consisting of a series of I-shaped structures was proposed in [16] to independently manipulate surface waves and spatially propagating waves. However, once this type of metasurface is fabricated, the generated beams are typically static and cannot be adjusted in real time. In the case of beamforming, it is necessary to steer the pattern radiated by the antenna in real time; obviously, the antenna with a fixed beam cannot achieve this.
In this paper, a single-layer reconfigurable hybrid metasurface is innovatively proposed for the realization of a compact, real-time and low-cost MSA. Firstly, in order to realize the compact structure, a planar microstrip feed antenna is designed. Secondly, three kinds of metasurface elements with different phase and amplitude responses to electromagnetic waves are proposed and studied. Then, the three elements are combined to form different hybrid metasurfaces. Under the irradiation of the feed antenna, the manipulation of the far-field radiation beam using the hybrid metasurface is studied. Finally, after optimization, by adding sixteen PIN diodes in the hybrid metasurface, real-time switching between two beams and four beams can be achieved. The proposed multibeam switching antenna features flexible beam switching, a small size and a simple feeding network with low costs. The idea of hybrid units can be applied to the design of a large aperture metasurface.

Design of the Feed Antenna
In order to realize the low-profile characteristic of the antenna, a planar microstrip antenna is adopted instead of the horn as the metasurface feed. Figure 1a shows the front structure of the feed antenna, with the yellow part being copper and the gray part being a 1.5 mm thick F4B dielectric substrate (ε r = 2.2; tanδ = 0.001). A U-shaped slot is etched on the radiation patch to improve the current distribution and achieve a wider bandwidth and higher gain than an ordinary microstrip without a slot. The back of the substrate is covered with copper to ensure that the antenna keeps the broadside radiation mode, as shown in Figure 1b. It can be seen from Figure 2a

Unit Design of the Hybrid Metasurface
The unit cells of the proposed hybrid metasurface are illustrated in Figure 3a-c. F4B (marked gray) with a dielectric constant of 2.2, a loss tangent of 0.001, and a thickness of 1 mm is used as the dielectric substrate. The thickness of the metal copper (marked yellow) is 0.035 mm. Two types of radiation structures are selected: a rectangular patch (unit 1) and two electric-inductive-capacitive (ELC) resonator structures (unit 2 and unit 3). Unit 1 is a common structure that makes up a metasurface, and different phase responses to transmitted electromagnetic waves can be achieved by adjusting the size of the rectangular patch. Unit 2 has resonance characteristics at a certain frequency bandwidth, and if it is considered as a component of the metasurface, it has filtering characteristics for specific polarized electromagnetic waves during resonance [17]. By combining the above two units into a metasurface, the phase and amplitude of the transmitted beam can be adjusted simultaneously. Unit 3 has made the following improvements on the basis of unit 2: two structures of unit 2 are merged to achieve a more compact metasurface, and diodes are loaded on the two wide arms in the middle to achieve real-time manipulation of electromagnetic waves. A commercially available "SMP 1320-079LF" is chosen as the switchable diode as it has low resistance and low capacitance. According to the datasheet [18], the forward voltage of the diode is 0.85 V, and for the convenience of testing, a 1.5 V battery is used for DC bias. We define the state of the diode as "ON" when 1.5 V voltage is loaded; when 0 V voltage is loaded, the state is defined as "OFF".

Unit Design of the Hybrid Metasurface
The unit cells of the proposed hybrid metasurface are illustrated in Figure 3a-c. F4B (marked gray) with a dielectric constant of 2.2, a loss tangent of 0.001, and a thickness of 1 mm is used as the dielectric substrate. The thickness of the metal copper (marked yellow) is 0.035 mm. Two types of radiation structures are selected: a rectangular patch (unit 1) and two electric-inductive-capacitive (ELC) resonator structures (unit 2 and unit 3). Unit 1 is a common structure that makes up a metasurface, and different phase responses to transmitted electromagnetic waves can be achieved by adjusting the size of the rectangular patch. Unit 2 has resonance characteristics at a certain frequency bandwidth, and if it is considered as a component of the metasurface, it has filtering characteristics for specific polarized electromagnetic waves during resonance [17]. By combining the above two units into a metasurface, the phase and amplitude of the transmitted beam can be adjusted simultaneously. Unit 3 has made the following improvements on the basis of unit 2: two structures of unit 2 are merged to achieve a more compact metasurface, and diodes are loaded on the two wide arms in the middle to achieve real-time manipulation of electromagnetic waves. A commercially available "SMP 1320-079LF" is chosen as the switchable diode as it has low resistance and low capacitance. According to the datasheet [18], the forward voltage of the diode is 0.85 V, and for the convenience of testing, a 1.5 V battery is used for DC bias. We define the state of the diode as "ON" when 1.5 V voltage is loaded; when 0 V voltage is loaded, the state is defined as "OFF".

Unit Design of the Hybrid Metasurface
The unit cells of the proposed hybrid metasurface are illustrated in Figure 3a-c. F4B (marked gray) with a dielectric constant of 2.2, a loss tangent of 0.001, and a thickness of 1 mm is used as the dielectric substrate. The thickness of the metal copper (marked yellow) is 0.035 mm. Two types of radiation structures are selected: a rectangular patch (unit 1) and two electric-inductive-capacitive (ELC) resonator structures (unit 2 and unit 3). Unit 1 is a common structure that makes up a metasurface, and different phase responses to transmitted electromagnetic waves can be achieved by adjusting the size of the rectangular patch. Unit 2 has resonance characteristics at a certain frequency bandwidth, and if it is considered as a component of the metasurface, it has filtering characteristics for specific polarized electromagnetic waves during resonance [17]. By combining the above two units into a metasurface, the phase and amplitude of the transmitted beam can be adjusted simultaneously. Unit 3 has made the following improvements on the basis of unit 2: two structures of unit 2 are merged to achieve a more compact metasurface, and diodes are loaded on the two wide arms in the middle to achieve real-time manipulation of electromagnetic waves. A commercially available "SMP 1320-079LF" is chosen as the switchable diode as it has low resistance and low capacitance. According to the datasheet [18], the forward voltage of the diode is 0.85 V, and for the convenience of testing, a 1.5 V battery is used for DC bias. We define the state of the diode as "ON" when 1.5 V voltage is loaded; when 0 V voltage is loaded, the state is defined as "OFF".
The transmission characteristics are simulated using ANSYS Electromagnetics Suite. Each unit is simulated using periodic boundary conditions and Floquet ports, as shown in Figure 3d. The incident uniform plane wave is a TE linearly polarized wave. During the simulation, we used the equivalent circuit shown in Figure 3e to replace the OFF or ON state of the diode. From Figure 4a,b, we know that unit 1 and unit 2 have different amplitude and phase responses to electromagnetic waves, enabling the manipulation of electromagnetic waves from two dimensions. Unit 2 and unit 3 (in OFF state) have almost the same amplitude and phase properties at 11.2 GHz. In terms of amplitude, when all PIN diodes are in the ON state, the transmissivity of the TE wave in unit 3 is almost 0 at 11.2 GHz. When all PIN diodes are in the OFF state, the transmissivity is 0.96. When electromagnetic waves are incident from different angles, the phase and amplitude responses of the units remain relatively stable. It should be noted that due to the similar characteristics of unit 3 and unit 2, Figure 4c,d only show the situation of unit 1 and unit 2. The amplitude difference close to 1 will have a great impact on electromagnetic waves and help the metasurface to manipulate the transmitted waves. The transmission characteristics are simulated using ANSYS Electromagnetics Suite. Each unit is simulated using periodic boundary conditions and Floquet ports, as shown in Figure 3d. The incident uniform plane wave is a TE linearly polarized wave. During the simulation, we used the equivalent circuit shown in Figure 3e to replace the OFF or ON state of the diode. From Figure 4a,b, we know that unit 1 and unit 2 have different amplitude and phase responses to electromagnetic waves, enabling the manipulation of electromagnetic waves from two dimensions. Unit 2 and unit 3 (in OFF state) have almost the same amplitude and phase properties at 11.2 GHz. In terms of amplitude, when all PIN diodes are in the ON state, the transmissivity of the TE wave in unit 3 is almost 0 at 11.2 GHz. When all PIN diodes are in the OFF state, the transmissivity is 0.96. When electromagnetic waves are incident from different angles, the phase and amplitude responses of the units remain relatively stable. It should be noted that due to the similar characteristics of unit 3 and unit 2, Figure 4c,d only show the situation of unit 1 and unit 2. The amplitude difference close to 1 will have a great impact on electromagnetic waves and help the metasurface to manipulate the transmitted waves.

Design and Optimization of Multibeam Metasurface Antenna
It can be seen from Section 2.2 that there are phase differences among the three types of unit cells, and they have different transmittances. Therefore, by arranging these three units, a metasurface with an amplitude and phase gradient can be formed to control the

Design and Optimization of Multibeam Metasurface Antenna
It can be seen from Section 2.2 that there are phase differences among the three types of unit cells, and they have different transmittances. Therefore, by arranging these three units, a metasurface with an amplitude and phase gradient can be formed to control the antenna radiation beam. The effect of the metasurface with different combinations on the far-field radiation pattern is studied by placing the metasurface at the same height above the feed antenna and keeping other conditions unchanged.
As shown in Figure 5a,b, when the metasurface only includes unit 2, it can be seen that the antenna only radiates one beam and has strong scattering. When the metasurface is combined with units 1 and 2, as shown in Figure 5c,d, the antenna exhibits a zero point in the theta = 0 direction, and the directionality is relatively enhanced compared with case 1. As shown in Figure 5e,f, after replacing unit 2 in the middle with unit 3, the beam radiated by the antenna is more symmetrical compared with that in case 2. From the above three cases, it can be seen that when the hybrid metasurface consists of different units and is irradiated with the planar feed, it tends to realize multiple beams due to the difference in phase and amplitude. After heuristic optimization, the metasurface antenna as shown in Figure 6 is finally obtained. The proposed switched-beam antenna consists of two layers, as shown in Figure  6a. The bottom layer is a microstrip antenna (see Figure 1), and the upper layer is the hybrid metasurface. Eight PIN diodes are placed in the cells of the first row of the hybrid metasurface. And the other eight PIN diodes are placed in the last row with simple bias lines on the back (see Figure 6). The bias lines in Figure 6c are connected to the units with After heuristic optimization, the metasurface antenna as shown in Figure 6 is finally obtained. The proposed switched-beam antenna consists of two layers, as shown in Figure 6a. The bottom layer is a microstrip antenna (see Figure 1), and the upper layer is the hybrid metasurface. Eight PIN diodes are placed in the cells of the first row of the hybrid metasurface. And the other eight PIN diodes are placed in the last row with simple bias lines on the back (see Figure 6). The bias lines in Figure 6c are connected to the units with PIN diodes on the front side through a metal via. The shorter wire is connected to the positive pole of the DC circuit. To reduce the impact of DC bias on microwave signals, the width of the bias line is only 0.2 mm for high impedance. At the same time, in order to prevent microwave signals from entering the DC path, the bias line is connected to the corner with the lowest voltage in each unit. In order to be as realistic as possible, the bias lines and the metal vias are taken into the simulation. using electromagnetic waves radiated by a planar microstrip feed. These reasons result in a red shift in the operating frequency of the antenna.
In order to obtain better radiation performance, the parameter h3 is analyzed. The boundary between the induction near-field and Fresnel region is taken as the original value of h3, which is calculated as (1).
where L is the effective length of the source antenna, and λ is the wavelength in free space at 9.28 GHz. According to Equation (1), the initial value of h3 is 7.4 mm. In order to reduce the antenna profile while maintaining good impedance matching and a low side lobe, the radiation patterns with h3 of 3 mm, 5 mm and 7 mm are analyzed. Figure 8 shows the twodimensional pattern of the antenna when all diodes are OFF. When h3 is equal to 3 mm, as shown with the red curve in Figure 8, the gain of the antenna is lower in the YOZ-plane. When h3 is equal to 7 mm, as shown with the green curve in the XOZ-plane, the larger tail lobe of the antenna results in a lower gain. Overall, the antenna radiation performance is the best when h3 is equal to 5 mm. Table 1 shows the final parameter values after optimization. The two layers are separated by an air cavity with the height of h 3 . As shown in Figure 7, when all PIN diodes are turned on, the antenna radiates two beams; when all PIN diodes are in the OFF state, the antenna radiates four beams. It must be noted that the patterns in Figure 7 are observed at 9.28 GHz, which is lower than the results discussed in Figure 4. This is due to the non-periodic arrangement of the units, which is different from Figure 3. Another reason is the difference in incident plane waves. The simulation of the unit uses uniform plane wave incidence, while the proposed metasurface is illuminated using electromagnetic waves radiated by a planar microstrip feed. These reasons result in a red shift in the operating frequency of the antenna.

Mechanism of Multibeam Switching Antenna Based on the Reconfigurable Hybrid Metasurface
As shown in Figure 4, the hybrid metasurface consists of three types of units that have different transmission characteristics. We study it from the aspect of phase and amplitude, respectively. From the perspective of phase, according to the generalized Snell's law [19], i.e., Equation (1) below, when there is a phase gradient dφ/dx, the incident electromagnetic wave will be transferred to anomalous refraction. When the incident angle θi is known, the refraction angle θt can be designed using the phase gradient dφ/dx caused by the interface.

=
(2) where ni and nt are the refractive indices of the two media. For the reconfigurable hybrid metasurface, the units are non-periodically arranged, and the phase gradient introduced is not continuous. From this analysis, we know that the proposed hybrid metasurface can undergo abnormal refraction of incident electromagnetic waves. In terms of amplitude, when all PIN diodes are in the ON state, the transmissivity of the TE wave in Unit 3 is almost 0. This means that when the electromagnetic wave In order to obtain better radiation performance, the parameter h 3 is analyzed. The boundary between the induction near-field and Fresnel region is taken as the original value of h 3 , which is calculated as (1).
where L is the effective length of the source antenna, and λ is the wavelength in free space at 9.28 GHz. According to Equation (1), the initial value of h 3 is 7.4 mm. In order to reduce the antenna profile while maintaining good impedance matching and a low side lobe, the radiation patterns with h 3 of 3 mm, 5 mm and 7 mm are analyzed. Figure 8 shows the two-dimensional pattern of the antenna when all diodes are OFF. When h 3 is equal to 3 mm, as shown with the red curve in Figure 8, the gain of the antenna is lower in the YOZ-plane. When h 3 is equal to 7 mm, as shown with the green curve in the XOZ-plane, the larger tail lobe of the antenna results in a lower gain. Overall, the antenna radiation performance is the best when h 3 is equal to 5 mm. Table 1 shows the final parameter values after optimization.  As shown in Figure 4, the hybrid metasurface consists of three types of units that have different transmission characteristics. We study it from the aspect of phase and amplitude, respectively. From the perspective of phase, according to the generalized Snell's

Mechanism of Multibeam Switching Antenna Based on the Reconfigurable Hybrid Metasurface
As shown in Figure 4, the hybrid metasurface consists of three types of units that have different transmission characteristics. We study it from the aspect of phase and amplitude, respectively. From the perspective of phase, according to the generalized Snell's law [19], i.e., Equation (1) below, when there is a phase gradient dϕ/dx, the incident electromagnetic wave will be transferred to anomalous refraction. When the incident angle θ i is known, the refraction angle θ t can be designed using the phase gradient dϕ/dx caused by the interface.
where n i and n t are the refractive indices of the two media. For the reconfigurable hybrid metasurface, the units are non-periodically arranged, and the phase gradient introduced is not continuous. From this analysis, we know that the proposed hybrid metasurface can undergo abnormal refraction of incident electromagnetic waves. In terms of amplitude, when all PIN diodes are in the ON state, the transmissivity of the TE wave in Unit 3 is almost 0. This means that when the electromagnetic wave irradiates the unit with low transmittance, the electromagnetic wave will be reflected back for secondary radiation. According to Huygens's principle, the radiated electromagnetic wave will be continuously superimposed to form a new radiation pattern. The amplitude difference and phase shift can be used to control the direction and number of the radiation beams. Generally, the transmission wave from the hybrid metasurface under plane-wave illumination can be calculated as the following equation [20]: where A(m, n) and ϕ(m, n) denote the transmission amplitude and phase from the (m, n) unit cell, d is the periodicity of the unit and (φ, θ) represent the azimuth and elevation angles of the radiation pattern, respectively. To better reveal the radiation properties of the multibeam switching antenna, the current distribution on the hybrid metasurface is studied as shown in Figure 9. It can be seen that when all PIN diodes are in the ON state, the current density of unit 1 and unit 2 is low, while the current distribution on unit 3 located on both sides is large. At this time, two beams in the YOZ-plane were obtained. When all PIN diodes are in the OFF state, the current density of unit 2 and unit 1 in the middle row obviously increases, while the current on both sides remains large. In this case, the antenna not only maintains two radiation beams on the YOZ-plane but also radiates two beams on the XOZ-plane. irradiates the unit with low transmittance, the electromagnetic wave will be reflected back for secondary radiation. According to Huygens's principle, the radiated electromagnetic wave will be continuously superimposed to form a new radiation pattern. The amplitude difference and phase shift can be used to control the direction and number of the radiation beams. Generally, the transmission wave from the hybrid metasurface under plane-wave illumination can be calculated as the following equation [20]: where A(m, n) and φ(m, n) denote the transmission amplitude and phase from the (m, n) unit cell, d is the periodicity of the unit and (ϕ, θ) represent the azimuth and elevation angles of the radiation pattern, respectively. To better reveal the radiation properties of the multibeam switching antenna, the current distribution on the hybrid metasurface is studied as shown in Figure 9. It can be seen that when all PIN diodes are in the ON state, the current density of unit 1 and unit 2 is low, while the current distribution on unit 3 located on both sides is large. At this time, two beams in the YOZ-plane were obtained. When all PIN diodes are in the OFF state, the current density of unit 2 and unit 1 in the middle row obviously increases, while the current on both sides remains large. In this case, the antenna not only maintains two radiation beams on the YOZ-plane but also radiates two beams on the XOZ-plane.

Experiments and Results
As illustrated in Figure 10, the hybrid metasurface antenna was fabricated and tested. The total dimension is 58 mm × 49 mm × 7 mm (1.80λ0 × 1.52λ0 × 0.22λ0). The bias lines are placed on the back of the metasurface through the vias. Figure 11a presents the comparison of the reflection coefficient between measurement and simulation, indicating that the measured −10 dB impedance bandwidth is 9.07-9.42 GHz when all PIN diodes are in the ON state. The measured bandwidth is narrower than the simulation because of manufacturing errors.

Experiments and Results
As illustrated in Figure 10, the hybrid metasurface antenna was fabricated and tested. The total dimension is 58 mm × 49 mm × 7 mm (1.80λ 0 × 1.52λ 0 × 0.22λ 0 ). The bias lines are placed on the back of the metasurface through the vias. Figure 11a presents the comparison of the reflection coefficient between measurement and simulation, indicating that the measured −10 dB impedance bandwidth is 9.07-9.42 GHz when all PIN diodes are in the ON state. The measured bandwidth is narrower than the simulation because of manufacturing errors. beams are 5.54 dBi and 5.37 dBi at the H-plane and 5.76 dBi and 5.98 dBi at the E-plane, respectively. The measured gain of cross-polarization is lower than −12 dBi. It can be seen that the measured gain is lower than the simulation, and the side lobe is higher than the simulation results. This is mainly caused by the inserted loss of PIN diodes and the influence of the bias lines. Nevertheless, the measured results are basically consistent with the simulation results, indicating that the proposed method based on a reconfigurable hybrid metasurface is feasible. beams are 5.54 dBi and 5.37 dBi at the H-plane and 5.76 dBi and 5.98 dBi at the E-plane, respectively. The measured gain of cross-polarization is lower than −12 dBi. It can be seen that the measured gain is lower than the simulation, and the side lobe is higher than the simulation results. This is mainly caused by the inserted loss of PIN diodes and the influence of the bias lines. Nevertheless, the measured results are basically consistent with the simulation results, indicating that the proposed method based on a reconfigurable hybrid metasurface is feasible.  Figure 11b gives the simulated radiation efficiency, from which it can be seen that in the working band of the antenna, the radiation efficiency under two states is greater than 80%. Although the amplitude is modulated by the proposed metasurface, the total radiation efficiency of the antenna is not affected since the transmission amplitude of unit 2 and unit 3 (in the OFF state) is close to 1. Figure 12 shows the 2-D realized gain patterns at 9.28 GHz. It is worth noting that only the patterns in upper-half-plane were measured for their directional radiation properties.
When all PIN diodes are in the ON state, the antenna radiates two beams in the E-plane, which is located at theta = 38 • , 314 • in polar coordinate with the measured gain of 6.77 dBi and 7.32 dBi, respectively. When all PIN diodes are in the OFF state, the antenna generates four beams, two of which are located in the H-plane with theta = 48 • , 312 • , and the other two are in the E-plane with theta = 42 • , 314 • . The measured gains of the four beams are 5.54 dBi and 5.37 dBi at the H-plane and 5.76 dBi and 5.98 dBi at the E-plane, respectively. The measured gain of cross-polarization is lower than −12 dBi. It can be seen that the measured gain is lower than the simulation, and the side lobe is higher than the simulation results. This is mainly caused by the inserted loss of PIN diodes and the influence of the bias lines. Nevertheless, the measured results are basically consistent with the simulation results, indicating that the proposed method based on a reconfigurable hybrid metasurface is feasible.

Performance Comparison
To further illustrate the advantages of the proposed multibeam switching antenna, the comparisons with the multibeam antennas based on metasurfaces are listed in Table  2. The antenna presented in [15] with a smaller size achieves two beams of radiation, but it cannot switch beams. Antennas using multiple ports to achieve multiple beams were proposed [21][22][23][24], but they can only achieve one beam at the same time, and the real-time switching of different ports is relatively complex. The antenna proposed in [25] achieved flexible beam control by placing programmable metasurface elements in SIW antennas, but the size of the antenna is relatively large. And as the array expands, the bias network of PIN diodes becomes complex. Our proposed multibeam switching antenna features flexible beam switching, a small size and a simple feeding network with a low cost. The idea of hybrid units can be applied to the design of a large aperture metasurface.

Performance Comparison
To further illustrate the advantages of the proposed multibeam switching antenna, the comparisons with the multibeam antennas based on metasurfaces are listed in Table 2. The antenna presented in [15] with a smaller size achieves two beams of radiation, but it cannot switch beams. Antennas using multiple ports to achieve multiple beams were proposed [21][22][23][24], but they can only achieve one beam at the same time, and the real-time switching of different ports is relatively complex. The antenna proposed in [25] achieved flexible beam control by placing programmable metasurface elements in SIW antennas, but the size of the antenna is relatively large. And as the array expands, the bias network of PIN diodes becomes complex. Our proposed multibeam switching antenna features flexible beam switching, a small size and a simple feeding network with a low cost. The idea of hybrid units can be applied to the design of a large aperture metasurface.

Conclusions
A novel reconfigurable hybrid metasurface antenna achieving multibeam switching is presented in this paper. As an elaborate arrangement of three types of metasurface elements, fewer PIN diodes are used, and the bias network is simple. The measured results are consistent with simulations, indicating that the antenna can realize dynamic switching between two beams and four beams with a compact structure. The proposed method is suggested for beamforming systems. And the possibility of flexible real-time control of more beams in a broad bandwidth will be investigated in future work.
Author Contributions: All authors have significantly contributed to the research presented in this manuscript; L.S. presented the main idea and wrote the manuscript; G.N. and Y.L. edited the manuscript; W.C. and L.M. reviewed and revised the manuscript. All authors have read and agreed to the published version of the manuscript.