Miura origami based reconfigurable polarization converter for multifunctional control of electromagnetic waves

Zhibiao Zhu; Yongfeng Li; Zhe Qin; Lixin Jiang; Wenjie Wang; Hongya Chen; Jiafu Wang; Yongqiang Pang; Shaobo Qu

doi:10.1364/PRJ.504027

Note: This section is automatically generated by AI . The website and platform operators shall not be liable for any commercial or legal consequences arising from your use of AI generated content on this website. Please be aware of this.

Abstract

Polarization is one of the basic characteristics of electromagnetic (EM) waves, and its flexible control is very important in many practical applications. At present, most of the multifunction polarization metasurfaces are electrically tunable based on PIN and varactor diodes, which are easy to operate and have strong real-time performance. However, there are still some problems in them, such as few degrees of freedom of planar structure control, complex circuit, bulky sample, and high cost. In view of these shortcomings, this paper proposes a Miura origami based reconfigurable polarization conversion metasurface for multifunctional control of EM waves. The interaction between the electric dipoles is changed by adjusting the folding angle

θ

, thereby tuning the operating frequency of the polarization conversion and the polarization state of the reflected wave. This mechanical control method brings more degrees of freedom to manipulate EM waves. And the processed sample is with lightweight and low cost. To verify the performance of the proposed origami polarization converter, a Miura origami structure loaded with metal split rings is designed and fabricated. The operating frequency of the structure can be tuned in different folding states. In addition, by controlling the folding angle

θ

, linear-to-linear and linear-to-circular polarization converters can be realized at different folding states. The proposed Miura origami polarization conversion metasurface provides a new idea for reconfigurable linear polarization conversion and multifunctional devices.

1. INTRODUCTION

Recently, metasurfaces have developed into a new research field. Metasurfaces can be designed to be materials with arbitrary permittivity and permeability, manipulating electromagnetic (EM) waves in ways that have never been possible before. In addition, metasurfaces are gradually developing new structural forms, which can achieve more and more diverse functions. Examples include holography, vortex beams for wireless communication systems [1], radar cross section (RCS) reduction [2], focusing [3 –6], and polarization conversion [7 –17]. Polarization is one of the important characteristics of EM waves, which can be applied to communication and radar systems [18]. Polarization clearly has significant research value. Polarization conversion metasurfaces experienced the development process of passive and active polarization conversion metasurfaces. They are classified into two types: reflection type and transmission type polarization conversion metasurfaces. These polarization conversion metasurfaces enable functions such as linear to linear [19], linear to circular [7], and circular to circular [20] conversions. For the active polarization conversion metasurface, the switching between polarization conversion functions can be realized.

With the development of intelligent technology in recent years, metasurfaces are required to flexibly control EM waves by sensing environmental changes and reconfiguring structural elements. The active metasurface uses a varactor diode, PIN diode, and variable resistance diode to electrically control the parameters. For example, Tian et al. [21] used PIN diodes as switches to realize real-time and fast conversion of linear-to-linear polarization and linear-to-circular polarization in 2019. The method of active metasurface to realize related functions is simple, but the limitations of complex circuits and high cost make it unable to be better utilized in a large area in practical application scenarios. In addition, the control degrees of freedom (DOF) of active metasurfaces are related to the number and type of active devices. When the control DOF of active metasurfaces increases, the complexity of the structure and the number of devices increase. Thus, the cost increases significantly.

The emergence of paper-cut/origami metamaterials [22 –25] has solved the above problems and achieved lightweight structure while achieving multi-DOF control. Origami is a very promising technique for realizing reconfigurable metamaterials. This is the ancient art of folding a piece of paper into a complex and delicate 3D object. In recent years, origami technology has attracted great interest from scientists and engineers due to its ability to create materials with various shape variations. Origami has excellent mechanical properties and provides more control freedom. In recent years, origami has been gradually applied to the EM field. By changing the folding state of origami and loading the origami structure of metallic patterns, different functions are realized. For example, Li et al. [24] proposed a multifunctional metamaterial based on Miura origami, which realized three functions of the absorber, mirror reflector, and negative reflector by changing the folding angle. Le et al. [25] proposed a multifunctional metamaterial based on paper-cut metamaterials, which changed the folding state of paper folding through mechanical stimulation, so as to realize three functions: a bandpass filter, absorber, and reflector. Wang et al. [26] reported a metamaterial based on Miura origami, which results in a strong chiral response by deforming the origami structure. However, few scientists have reported in detail how origami modulates the polarization of EM waves.

Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you！Sign up now

In this paper, a reconfigurable polarization conversion metasurface based on Miura origami is designed for multifunctional control of EM waves. Among the origami types, Miura origami is chosen for the following reasons: the reconfigurability of the continuous state range and the large surface reconfigurable property of folding to unfolded state. The fabricated Miura origami substrate is loaded with metallic split rings. The shape of Miura origami is regulated by mechanical stimulation. After changing the folding angle of the origami structure, the interaction between the electric dipoles on the surface of the structure changes accordingly. Thus, the multifunctional control of EM waves is realized by the structure, as shown in Fig. 1. One function of the origami polarization conversion metasurface is that the operating frequency of polarization conversion can shift when the folding angle is changed. In addition, the origami polarization conversion metasurface can also control the polarization state of the reflected EM wave. When the folding angle is $θ = 15 °$ , the origami polarization conversion metasurface can be regarded as a linear-to-circular polarization converter. When the folding angle is $θ = 75 °$ , this origami polarization conversion metasurface can be regarded as a linear-to-linear polarization converter. The proposed polarization conversion metasurface achieves multifunctional control of EM waves under the premise of lightweight structure and low processing cost. This provides reference for the development of multifunctional metasurfaces in the future.

Figure 1.Polarization converter based on Miura origami can achieve multiple functional control of EM waves. First, in different folding states of origami, the polarization converter operates at different frequencies. In addition, by adjusting the folding angle of the origami structure, linear-to-linear and linear-to-circular polarization converters can be realized under the incidence of linear polarization waves, respectively.

2. ORIGAMI STRUCTURE DESIGN AND SIMULATION

CST full-wave simulation software is used to simulate the EM response of the origami polarization converter. The simulation process is as follows. (i) CST solver is set as frequency domain solver. (ii) The frequency range is set as 0–5 GHz. (iii) We set two Floquet port modes, and both ports are linearly polarized wave excitation in different folding states of the structure. (iv) Periodic boundaries are set in the $x$ and $y$ directions, and open (added space) boundaries are set in the $Z_{\max}$ and $Z_{\min}$ directions. After setting the simulation conditions, the model is established according to the structural relations of Miura origami. In order to reduce the computational cost of simulation, we redefined the unit cell, as shown in Figs. 2(a) and 2(b). We intercept the period of a rectangle on consecutive adjacent Miura origami structures. As the structure is designed to work under reflective conditions, there is a metal back plate on the bottom of the origami. The air gap exists between the metal back plate and the origami polarization converter, with a thickness of $d_{1} = 2 mm$ . In the experimental stage, this layer of air gap is replaced with foam.

Figure 2.(a) Simulation setting of planar origami polarization converter ( $θ = 0 °$ ). (b) Simulation setting of nonplanar origami polarization converter ( $θ \neq 0 °$ ). (c) Structural parameters of Miura origami substrate. (d) Dimension parameters of metallic split ring.

The unit cell of the origami polarization converter is composed of four symmetrical parallelograms arranged periodically. For convenience, the elements of each parallelogram are called facets. The facet design parameters are shown in Figs. 2(c) and 2(d). The geometric parameters of the facet are $a = b = 40 mm$ , $α = 60 °$ . The outside diameter of the split ring is $r_{1} = 13 mm$ . The inside diameter of the ring is $r_{2} = 12 mm$ . The width of opening is $g = 5 mm$ . The folding state of the reconfigurable polarization converter is determined by the folding angle $θ$ . When $θ = 0 °$ , the origami polarization converter is in the flat state. When $θ$ gradually increases, the lengths of the parameters $l_{x}$ and $l_{y}$ change accordingly. In the following formula, $l_{x}$ and $l_{y}$ can be represented by $θ$ as [22] ${\begin{matrix} l_{x} = 2 b (\sin δ \cos α + \cos δ \cos θ \sin α) \\ l_{y} = 2 a \cos δ \end{matrix},$ (1)where $δ = \arctan (\frac{1}{\tan α \cos θ})$ . $l_{x}$ and $l_{y}$ are the length parameters of the unit mentioned above. By changing the folding angle $θ$ , $l_{x}$ and $l_{y}$ change in size. A foldable polarization converter is investigated without tearing or stretching the substrate. It is worth noting that there are several parameters that determine the folding state of origami. In order to ensure that the desired frequency converges at the desired location, there can only be one variable ( $θ$ ) that describes the folded state of the origami polarization converter. The frequency shift caused by change of the $θ$ is calculated as follows: $\frac{Δ f}{f_{mid}} = \frac{2 | f_{θ 1} - f_{θ 2} |}{| f_{θ 1} + f_{θ 2} |},$ (2)where $Δ f$ is the frequency shift produced by $θ$ changing to $θ_{2}$ from $θ_{1}$ , and $f_{mid}$ is the midpoint frequency between $f_{θ 1}$ and $f_{θ 2}$ .

As shown in Fig. 3(a), as the folding angle of origami polarization converter increases, the working frequency gradually moves to low frequency ( $f_{1} = 3.4 GHz$ , $f_{2} = 3.8 GHz$ , and $f_{3} = 4.1 GHz$ ). In addition, the physical mechanism of operating frequency migration is discussed by analyzing the surface current distribution of metallic patterns, as illustrated in Figs. 3(b)–3(d). When the origami structure is in any folded state, the surface currents on both arms of the split ring oscillate back and forth along the metallic arms. Remarkably, the current oscillates in the same direction on both arms. Thus, surface currents on each metallic spilt ring can be equivalent to the electric dipole. Four metallic spilt rings are equivalent to four electric dipoles ( $p_{1}$ , $p_{2}$ , $p_{3}$ , and $p_{4}$ ). For example, for $p_{1}$ and $p_{2}$ , the interaction energy between them is $V = (p_{1} \cdot p_{2}) / 4 π ε_{0} r^{3}$ . As the folding angle $θ$ increases, the distance between the electric dipoles $p_{1}$ and $p_{2}$ decreases. In this case, the strength of the surface current on the metallic rings increases. When the folding angle is $θ = 15 °$ , the surface currents on the split rings are weak and the attraction between the dipoles is weak, as shown in Fig. 3(b). The origami polarization converter works in the high frequency band ( $f_{3}$ ). As the folding angle $θ$ increases, the surface currents on the split rings gradually strengthen [Figs. 3(c) and 3(d)]. The two dipoles on the origami surface become more attractive to each other, making the system more stable. As a result, the operating frequency is shifted to lower frequencies ( $f_{2}$ and $f_{3}$ ). The analysis of other electric dipole groups is similar.

Figure 3.(a) Simulated reflection spectra of the origami polarization converter. The surface current distributions on metallic parts under $y$ -polarization incidence: (b) $θ = 15 °$ , (c) $θ = 35 °$ , and (d) $θ = 55 °$ .

3. ANALYSIS AND DISCUSSION

Polarization is the phenomenon that the electric field vector of EM wave vibrates according to a certain law. Based on EM field theory, this paper analyzes the basic principle of circular polarization conversion and resonant frequency regulation after the EM wave incident on the origami structure. This provides the basis and guidance for the final design. Aiming at the polarization conversion efficiency of the origami circular polarization converter, we assume that a $y$ -polarized EM wave is incident on the surface of the structure from the $+ z$ direction. Hence, the electric field of the incident EM wave can be expressed as $E_{i} = E_{y} = \hat{y} E_{0} e^{i k_{z} z}$ . Furthermore, the electric field of the reflected wave can be expressed as [27] $E_{r} = E_{y} = \hat{y} R_{y y} E_{r} + \hat{x} R_{x y} E_{r},$ (3)where $R_{y y} = r_{y y} e^{i φ_{y y}}$ and $R_{x y} = r_{x y} e^{i φ_{x y}}$ , and $r_{y y}$ and $r_{x y}$ are the amplitudes of reflection cofficients $R_{y y}$ and $R_{x y}$ . Additionally, $φ_{x x}$ and $φ_{x y}$ are the phases of reflection waves. Then, $R_{LCP}$ and $R_{RCP}$ are used to denote the reflection efficiency of linear-to-circular polarization. $r_{LCP}$ and $r_{RCP}$ represent the corresponding reflection amplitude. $φ_{LCP}$ and $φ_{RCP}$ represent the corresponding phase of reflection. Meanwhile, the reflected electric field of LCP and RCP waves can be expressed as $E_{LCP}^{r} = \frac{\sqrt{2}}{2} R_{LCP} (E_{i} e^{i \frac{π}{2}} \hat{x} + E_{i} \hat{y}),$ (4) $E_{LCP}^{r} = \frac{\sqrt{2}}{2} R_{LCP} (E_{i} \hat{x} + E_{i} e^{i \frac{π}{2}} \hat{y}) .$ (5)

It is worth noting that $E_{r} = E_{LCP}^{r} + E_{RCP}^{r}$ . Combined with the above formula, the following formulas are deduced: $R_{LCP} = \frac{\sqrt{2}}{2} (R_{y y} - i R_{x y}),$ (6) $R_{RCP} = \frac{\sqrt{2}}{2} (R_{x y} - i R_{y y}) .$ (7)

It can be found that the amplitudes of reflection coefficients of co-polarization and cross-polarization should be equal and as high as possible to achieve an effective conversion from linear polarization to circular polarization. At the same time, the phase difference between the co- and cross-polarization reflection coefficients is $n π / 2$ ( $n$ is an integer). To study the polarization of reflected waves, four Stokes parameters are introduced as follows: ${\begin{matrix} I = {| r_{y y} |}^{2} + {| r_{x y} |}^{2} \\ Q = {| r_{x y} |}^{2} - {| r_{y y} |}^{2} \\ U = 2 | r_{y y} | | r_{x y} | \cos Δ φ \\ V = 2 | r_{y y} | | r_{x y} | \sin Δ φ \end{matrix} .$ (8)

As shown in Fig. 4(a), polarization azimuth $α$ and ellipticity $β$ describe the orientation of the principal axis of the ellipse and the shape of the ellipse. They can be obtained by the following formula: $\tan 2 α = \frac{U}{Q},$ (9) $\sin 2 β = \frac{V}{I} .$ (10)

Figure 4.(a) Sketch of the polarization azimuth angle $α$ and the ellipticity angle $β$ . (b) Co- and cross-polarized reflection amplitudes and phase difference between them under $y$ - polarized wave incidence. The folding angle $θ$ changes from 15° to 75° ( $f = 3.1 GHz$ , $ψ = 45 °$ ). (c) Calculated azimuth angle $α$ and the ellipticity angle $β$ at different folding angles. (d) Polarization ellipses of the reflected wave at different folding angles (75°, 65°, 55°, 45°, 30°, and 15°).

Figure 4(b) describes the changes of $r_{x y}$ and $r_{y y}$ amplitudes when the folding angle $θ$ changes from 15° to 75° ( $f = 3.1 GHz$ , $ψ = 45 °$ ). Meanwhile, the phase difference between the two reflected waves remains the same, which is 90°. According to Eqs. (9) and (10), the result shown in Fig. 4(c) can be obtained. The above derivation shows that when the folding angle of the origami structure changes, the polarization ellipticity angle at $f = 3.1 GHz$ can be adjusted from 0° to 45° ( $ψ = 45 °$ ). The azimuth angle $α$ of polarization remains constant when the folding angle $θ$ changes. Therefore, when the origami structure is dynamically regulated, the control from linear polarization to circular polarization can be realized under the $y$ -polarized wave incidence. Figure 4(d) shows the polarization states of the polarization ellipse under different folding angles. When the folding angle $θ = 75 °$ , the reflected wave is a linearly polarized wave. When the folding angle $θ = 15 °$ , the reflected wave is a circularly polarized wave. When $θ \in (15 °, 17 °)$ , the ellipticity angle $β$ changes with the fold angle $θ$ .

By analyzing the surface current distribution, the physical mechanism of polarization conversion under a different folding state of the structure is discussed. Decompose the electric field vector $E$ ( $ψ = 45 °$ ) into the vertical directions $x (E_{x})$ and $y (E_{y})$ , as shown in Fig. 4(c). Then, the commercial simulation software CST is used to obtain the simulation results of the surface current distribution under the incidence of $x$ -polarized and $y$ -polarized waves. Figure 5 shows the surface current distribution of the metallic spilt ring under two folding angles ( $θ = 15 °$ and $θ = 75 °$ ) at $f = 3.1 GHz$ . Under the incident of $x$ or $y$ polarized wave, strong surface currents are generated on the surface of the metallic split ring. And the surface currents on the two arms oscillate back and forth in the same direction. Therefore, the surface currents on the four split rings can be equivalent to four electric dipoles $p_{i}$ ( $i = 1$ , 2, 3, 4). According to Ref. [26], the total effective electric dipole is ${\bar{p}}_{eff} = \sum {\bar{p}}_{i} = p_{eff, x} \hat{x} + p_{eff, y} \hat{y} + p_{eff, z} \hat{z}$ . In addition, according to the calculation formula in Ref. [28], the total effective electric dipoles of the four electric dipoles have zero components in the $y$ and $z$ directions. It can be concluded that ${\bar{p}}_{eff} = p_{eff, x} \hat{x}$ . Therefore, the total effective electric dipole can be equivalent to an electric dipole in the $x$ direction.

Figure 5.(a) When the folding angle $θ = 15 °$ , the co-polarized reflection amplitudes and phases at 3.1 GHz under $x$ - and $y$ -polarized wave incidence. Current distributions on the surface of the metallic split ring at folding angle 15° ( $f = 3.1 GHz$ ), under two linearly polarized incident waves: (b) $x$ -polarization incidence, (c) $y$ -polarization incidence. (d) When the folding angle $θ = 75 °$ , the co-polarized reflection amplitudes and phases at 3.1 GHz under $x$ - and $y$ -polarized wave incidence. Current distributions on the surface of the metallic split ring at folding angle 75° ( $f = 3.1 GHz$ ), under two linearly polarized incident waves: (e) $x$ -polarization incidence, (f) $y$ -polarization incidence.

When the folding angle $θ = 15 °$ , the two total effective dipoles $p_{eff}^{a}$ and $p_{eff}^{b}$ in the $x$ direction radiate $x$ and $y$ polarized waves, respectively, with a phase difference $Δ φ$ of 90° and an equal amplitude ( $r_{x x} = r_{y y} = 1$ ). As shown in Figs. 5(b) and 5(c), in this case, the surface currents are strong under both $x$ - and $y$ -polarized waves. Polarization conversion occurs for both incident linearly polarized waves. Therefore, both $x$ and $y$ incident polarized waves undergo polarization conversion of 90°, and a circularly polarized wave is synthesized after the polarization conversion. As shown in Fig. 5(d), the total effective electric dipoles $p_{eff}^{c}$ and $p_{eff}^{d}$ emit $x$ - and $y$ -polarized waves, respectively. The phase difference $Δ φ$ of the two linearly polarized waves is 180°, and the amplitude is equal ( $r_{x x} = r_{y y}$ ). As shown in Figs. 5(e) and 5(f), when the folding angle $θ = 75 °$ , the surface currents are strong under $x$ -polarized wave irradiation, and the incident $x$ -polarized wave undergoes 90° polarization conversion. When illuminated by $y$ -polarized waves, the surface current is weak and no polarization conversion occurs for the incident wave. It can be concluded that at the folding angle $θ = 75 °$ , the incident $x$ -polarized wave undergoes polarization conversion. And it combines with another outgoing $y$ -polarized wave into a linearly polarized wave in the same direction.

4. FABRICATION AND MEASUREMENT

To verify the feasibility of the reconfigurable origami polarization converter, a sample is fabricated and measured, as shown in Fig. 6(a). In the sample processing link, the metallic split rings are printed on the cut polyimide films. Subsequently, the numerical control is used to etch creases on the polypropylene (PP) material substrate with a thickness of 0.05 mm. Finally, the polyimide films are pasted on the PP material substrate. The experimental setup consists of a pair of linearly polarized wideband horn antennas operating in the frequency range 2–18 GHz. In this pair of antennas, one is used as the transmitting antenna and the other is used as the receiving antenna. Meanwhile, the two antennas are connected to the two interfaces of the vector network analyzer to obtain more accurate analysis results. In addition, to minimize the impact of test noise, the measurements are chosen to perform in a microwave anechoic chamber. When setting up the antenna, the transmitting and receiving antennas are placed in a symmetric manner. Moreover, this set of antennas should have some space from the origami polarization converter under test. Setting the transmitting and receiving antennas at the same height allows the EM waves emitted by the transmitting antenna to be more effectively reflected by the origami device. At the same time, it ensures that the receiving antenna can fully receive the reflected EM waves, thus increasing the accuracy of the test results.

Figure 6.(a) Measurement experimental setup in microwave anechoic chamber. (b) Measured and simulated reflection spectra of the origami polarization converter under different folding states ( $θ = 15 °$ , 30°, and 55°). (c) When the folding angle $θ = 15 °$ , the measured co-polarized reflection amplitudes and phases at 3.1 GHz under $x$ - and $y$ -polarized wave incidence. (d) When the folding angle $θ = 75 °$ , the measured co-polarized reflection amplitudes and phases at 3.1 GHz under $x$ - and $y$ -polarized wave incidence.

As shown in Fig. 6(b), during the test, as the folding angle of the origami polarization converter continues to increase, the working frequency of the structure continues to move to low frequency. In addition, Fig. 6(c) shows that when the folding angle of the origination is 15°, the phase difference between the $x$ - and the $y$ -polarized waves remains 90° around 3.1 GHz. The reflection amplitude of both linearly polarized waves is close to one. Therefore, in this folded state, the reflected wave is circularly polarized. Figure 6(d) shows that when the folding angle of the origination is 75°, the phase difference between the $x$ - and the $y$ -polarized waves remains 180° around 3.1 GHz. The reflection amplitude of both linearly polarized waves is close to 1. Therefore, in this folded state, the reflected wave is linearly polarized. Since a flexible material is chosen as the substrate, the structure is not strictly rigid. And the deformation between different folded Miura origami structural elements may be slightly different. In addition, the experimental environment can also lead to differences between the experimental and the simulation results.

5. CONCLUSION

In summary, we propose a reconfigurable polarization conversion metasurface based on Miura origami. The surface current distribution simulation results show that the interaction between the electric dipoles can be changed by adjusting the folding angle $θ$ , so as to tune the operating frequency of polarization conversion and the polarization state of the reflected wave. In addition, by controlling the folding angle $θ$ , linear-to-linear and linear-to-circular polarization converters can be realized in different folding states. Compared with other polarization converters, this mechanical control method brings more DOFs to manipulate the EM wave. To verify the performance of the proposed origami polarization converter, a Miura origami polarization converter loaded with metal split rings is designed and fabricated. The processed samples are with lightweight and low cost. The proposed Miura origami polarization conversion metasurface provides a new idea for the design of reconfigurable antennas and multifunctional devices in the future.

References

[1] Q. Wu, W. Fan, C. Qin. Generation and superposition of perfect vortex beams in terahertz region via single-layer all-dielectric metasurface. Nanomaterials, 12, 3010(2022).

[2] Y. Gou, Y. Chen, S. Yang. Radar cross section reduction of wideband Vivaldi antenna arrays with array-level scattering cancellation. IEEE Trans. Antennas Propag., 70, 6740-6750(2022).

[3] L. Li, P. Zhang, F. Cheng. An optically transparent near-field focusing metasurface. IEEE Trans. Microwave Theory Tech., 69, 2015-2027(2021).

[4] M. Boyarsky, T. Sleasman, M. F. Imani. Electronically steered metasurface antenna. Sci. Rep., 11, 4693(2021).

[5] S. Vellucci, M. Longhi, A. Monti. Phase-gradient Huygens’ metasurface coatings for dynamic beamforming in linear antennas. IEEE Trans. Antennas Propag., 71, 7752-7765(2023).

[6] E. Martini, S. Maci. Theory, analysis, and design of metasurfaces for smart radio environments. Proc. IEEE, 110, 1227-1243(2022).

[7] H. L. Zhu, S. W. Cheung, K. L. Chung. Linear-to-circular polarization conversion using metasurface. IEEE Trans. Antennas Propag., 61, 4615-4623(2013).

[8] S. Teng, Q. Zhang, H. Wang. Conversion between polarization states based on a metasurface. Photonics Res., 7, 246-250(2019).

[9] X. Wu, Y. Meng, L. Wang. Anisotropic metasurface with near-unity circular polarization conversion. Appl. Phys. Lett., 108, 183502(2016).

[10] J. Lončar, A. Grbic, S. Hrabar. A reflective polarization converting metasurface at X-band frequencies. IEEE Trans. Antennas Propag., 66, 3213-3218(2018).

[11] P. C. Wu, R. Sokhoyan, G. K. Shirmanesh. Near-infrared active metasurface for dynamic polarization conversion. Adv. Opt. Mater., 9, 2100230(2021).