As one of the most rapidly expanding frontiers of modern photonics, metasurfaces hold promise for a wide range of applications due to their compact physical size and unconventional optical functionalities.1^–6 These flat and miniaturized devices usually consist of an array of optically thin meta-atoms (plasmonic or dielectric structures) with spatially varying geometrical parameters and subwavelength separations aimed at imprinting on the incoming optical wavefront a functionality7^–12 of choice through tailored local interactions of the phase,13^–15 amplitude,16^–18 and polarization.19^–21 The further development of this field entails not only endeavors to optimize the static metasurface performance for certain functionalities22^–24 but also most importantly techniques for tunable/reprogrammable metasurfaces that can be flexibly controlled on demand and in real time.25^–27 In the terahertz and optical regimes, dynamic metasurfaces have mainly been achieved by utilizing the response of tunable materials to external stimuli, such as optical pump,28^–31 heating,32^–34 and bias voltage.35^–38 Nevertheless, limited by the current techniques of fabrication and modulation, these metasurfaces have been limited to simultaneous tuning of a large portion of the meta-atoms, and thus can only be switched among a limited number of operations.

In the microwave regime, PIN diodes and varactor diodes have been recently used to design reprogrammable metasurfaces with individually reconfigured meta-atoms,39^–54 due to their miniaturized sizes compatible with passive metallic resonators. By judicious designs, several dynamic functionalities including holographic imaging,39 special beam steering,41^,43^,48^,53 and topological surface state control50 have been demonstrated. Moreover, digital and real-time features promote several system-level applications that can hardly be achieved by static metasurfaces, such as space-time-coding,42^,44^,54 intelligent autonomous self-adaptive systems,40^,45^,46 and dynamic optimization of wireless communication channels.47^,49^,51 Despite this progress, the degrees of freedom of PIN diode-based meta-atoms have been limited to two or four phase levels per meta-atom, i.e., the reprogrammability was limited to a locally controllable abrupt phase by a step of π or π/2 at best. This limitation causes a roughly discontinuous phase distribution that unavoidably introduces diffraction losses. Meanwhile, it also limits the amount of information that can be processed by a metasurface of given size.55 Programmable metasurfaces utilizing varactor diodes can enable phase levels larger than four;41^,46^,53 however, their operation was based on the dynamic phase around a resonance dip,38 resulting in the significant amplitude variation when adjusting the phase response. To eliminate such amplitude-phase correlation-induced distortion in specific wavefront control, due consideration should be paid into inverse design optimization algorithms.56^,57 Moreover, each meta-atom in electrically or optically reprogrammable metasurfaces requires at least one PIN diode, photodiode, or varactor diode,39^–54 whose typical power dissipation is about several hundred milliwatts, and they need a continuous power supply to maintain the functionality. Such energy consumption entails a trade-off between the size/number of meta-atoms and the overall metasurface size, hindering large-scale applications.

Notably, the shortcomings of the above-discussed microwave metasurfaces, including the meta-atom capability and the energy consumption, are mainly resulting from the limitations of the adopted voltage-driven elements. Improving the electrical metric of PIN diodes and varactor diodes, no doubt can optimize the metasurface performance, but it is an even harder challenge in electronic engineering, especially for higher frequencies.58 As thus, researchers have also carried out studies on the other physical mechanisms for achieving reconfigurable metasurfaces. For example, it has been demonstrated that an array of cavities can be filled with liquid mercury in a controllable manner using microfluidic technology in which the overall electromagnetic response can be tailored for different functionalities.59 In addition, origami60 and mechanical lifting61 approaches have also been reported for achieving reconfigurable metasurfaces where each of the meta-atoms is spatially folded/unfolded or lifted/dropped to tailor the overall metasurface functionality. Such approaches,59^–61 free from voltage-driven elements inside meta-atoms, hold an advantage of nonvolatility that the maintaining of designated functionality does not consume electricity. These works have played an important role in exploring the reconfigurable metasurfaces; however, the wavefront control capability remains insufficient.

In this article, we propose and demonstrate a reprogrammable metasurface platform based on mechanical control for quasicontinuous Pancharatnam–Berry (PB) phase tunability operating at microwave frequencies. PB phase,7^,8^,62^,63 also known as geometric phase, is a robust control method for incident circularly polarized waves, which is determined by the rotation angle of meta-atoms and it is therefore decoupled from amplitude control. Figure 1(a) schematically shows our reprogrammable PB metasurface platform, which consists of 20×20 supercells covering an area of 870mm×870mm. Each supercell, as shown in Fig. 1(b), is composed of a stepping motor, a set of transmission gears, and an array of 4×4 PB meta-atoms. The reprogrammable function is enabled by mechanically rotating each PB meta-atom to achieve the desired phase control, enabling a continuous and arbitrary phase control pattern over the entire surface with a high-efficiency and uniform amplitude. In principle, this metasurface platform can be employed for on-demand generation of a large variety of functionalities in real time. As proof-of-concept experiments, we demonstrate metalensing, focused vortex generation, and holographic imaging using the same metasurface, exhibiting superior performance and real-time tunability.

To empower the metasurface with well-controlled programmability, 400 commercial stepping motors are utilized to respectively control the 20×20 supercells. Each stepping motor is equipped with an addressed circuit for the power supply (terminal voltage of 5 V) and rotation control. The step angle is 5.625 deg, corresponding to 64 steps per turn. More details about the stepping motor can be found in Supplementary Material 1. The rotation step number (clockwise or counterclockwise) of each motor is controlled by wireless signals from a host computer with full addressability. As shown in Fig. 1(b), a three-layered gear set is utilized to transmit the torque from the motor to the PB meta-atoms: the first layer consists of the main gear attached to the output shaft of the stepping motor; the middle layer consists of four gears that connect the first and last layers; and the bottom layer consists of 16 gears connecting with the PB meta-atoms one by one. The modules of all gears are 0.5. The yellow-green gear has 30 teeth, the light-blue gear has 30 and 16 teeth in the bottom and top layers, respectively, and the pink gear has 14 teeth. In this condition, when the stepping motor rotates a turn, it drives the PB meta-atoms to rotate 8/7 turns. Accordingly, the rotation of each meta-atom consists of 56 steps per turn. For reference, Supplemental Video S1 (MP4, 30 MB [URL: https://doi.org/10.1117/1.AP.4.1.016002.1]) is provided to show the controllable rotations. The photos of the metasurface template and more details about each component can be found in Supplementary Material 2.

Different from traditional static metasurfaces, the meta-atoms of the proposed reprogrammable metasurface are discretely arranged to avoid possible obstructions during their rotation. With the overall performance taken into consideration, a circular-shaped metal–insulator–metal structure, as schematically shown in Fig. 2(a), is chosen as the meta-atom to manipulate the PB phase in reflection at an operating frequency of 7 GHz (see Supplementary Material 3 for simulation details). The PB meta-atom consists of two Archimedean spirals with the same geometric parameters, which are C2 rotation-symmetric with each other along the normal axis. Both the spirals and the bottom metallic film are made from 35-μm thick copper. They are fabricated using standard printed circuit board technology over a 3-mm-thick FR4 substrate, with relative permittivity of 4.2 and a loss tangent of 0.025. The radius of the meta-atom is R=4.8mm and the center distance between adjacent meta-atoms along the x and y directions is about 10.7 mm, leaving a minimum distance of about 1.1 mm to ensure smooth rotation of each meta-atom. Two important features characterize the chosen PB meta-atom in its interaction with the incident light: it provides strong anisotropy to achieve efficient polarization conversion, and it can strongly confine the resonance fields inside the meta-atom to minimize near-field coupling between neighboring meta-atoms (see Supplementary Material 3 for details). By carefully optimizing the meta-atom dimensions, such meta-atoms could efficiently convert the incident right- and left-handed circular polarizations (RCP and LCP) to reflected RCP and LCP, respectively. In principle, when this meta-atom is rotated by an angle θ, the phases of reflected Rrr and Rll (the subscripts r and l stand for RCP and LCP) are expected to follow a phase gradient of 2θ and −2θ, respectively, where the sign depends on the spin direction of incident polarization. Such geometric-dependent phase behavior is also known as the PB phase response.7^,8^,62^,63

Since the PB meta-atom can be rotated by 56 steps per turn, our metasurface enables quasicontinuous PB phase control over 28 levels. As schematically shown by the pink and sky-blue dials in Fig. 2(a), the reflectance Rrr and Rll have a phase control resolution of π/14 and opposite phase gradients. To experimentally characterize the metasurface, an experimental setup is established by a vector network analyzer (VNA, Agilent N5230C) and two broadband horn antennas (see Appendix). Figures 2(b) and 2(c) show the normalized amplitude of measured Rrr and Rll, respectively, where the amplitude peak around the target working frequency of 7 GHz can be obtained in both spectra. |Rrr| has a maximum amplitude of 0.91 at 6.925 GHz and |Rll| has a maximum amplitude of 0.7 at 6.55 GHz. Notably, when light impinges onto a chiral meta-atom that lacks mirror-symmetry, photons of different CPs will undergo different absorption, an effect known as circular dichroism or circular conversion dichroism.24 The chosen PB meta-atom shows a chiral response, thus we expect a difference between Rrr and Rll, as shown in Fig. 2. For different applications, one can either design a PB meta-atom with in-plane mirror-symmetry to eliminate circular conversion dichroism7 or design chiral PB meta-atoms with stronger symmetry breaking to enhance circular conversion dichroism.64 The red and blue circles in Figs. 2(d) and 2(e) illustrate the measured phase response of Rrr and Rll, respectively, at 7 GHz, as we simultaneously rotate all meta-atoms from 0 to π. Clearly, the phase variation of Rrr and Rll has a positive and negative gradient versus the rotation angle, respectively. The solid lines in Figs. 2(d) and 2(e) show the linearly fitted results, where the fitting gradient of Rrr and Rll is 1.9624 and −2.0665, respectively, which agree well with the theoretical value. In addition, Fig. S4 in the Supplementary Material shows the measured phase responses at 6, 6.7, and 7.4 GHz, which also comply well with the PB phase response. These results experimentally verify the PB phase controllability of the proposed metasurface platform in a broad band.

Since each supercell can be individually controlled, the metasurface yields a variety of optical functionalities by designing the rotation distribution over the surface. We define the metasurface as the xy-plane and its center the origin of our coordinate system. To experimentally characterize the performance, a broadband antenna is coaxially mounted at (0, 0, 2100 mm) as a feeding source, and a waveguide probe is used to scan the electric filed distribution (see Appendix). Since the wave emitted by the antenna is not a plane wave as it reaches the metasurface, we should compensate its phase to make it collimated (see Supplementary Material 5 for details). Notably, due to the quasicontinuous phase controllability, our metasurface can flexibly work under various incident wavefronts, since the functionality distortion produced by a nonplanar incident wavefront can be easily compensated. Such flexibility is quite important in practical applications, especially in scenarios in which the feeding antenna can hardly be placed in the Fraunhofer region of the metasurface.39^,42 Once the desired phase distribution for certain target functionality is obtained, the corresponding rotation distribution can be calculated. In the following, we choose RCP as target operating polarization to demonstrate several representative functionalities, including metalensing, focused vortex beam generation, and holographic imaging.

We first investigate the functionality of a metalens using the proposed metasurface. Figure 3(a) shows the required rotation profile for a metasurface focusing the impinging wave at (0, 0, 600 mm). The measured electric field intensity (|Rrr|2) in the xy-plane at z=600mm is shown in Fig. 3(b), where a highly symmetric focal spot can be obtained. The horizontal cut of the focal spot is shown in Fig. 3(c) showing a full-width at half-maximum (FWHM) ∼42mm. Notably, the numerical aperture of the demonstrated metalens is about 0.587, corresponding to a diffraction-limited FWHM of ∼36.5mm at 7 GHz. We remark that the generated focal spot is quite close to the diffraction limit, indicating the excellent performance of our quasicontinuous phase control. The programmability of our metasurface enables scanning the focus off-axis. Figures 3(d) and 3(g) show the desired rotation profiles for focusing at (−80, 0, 600 mm) and (120, 0, 600 mm), respectively. The measured electric field intensities and the horizontal cuts of the focal spots are shown in Figs. 3(e), 3(h) and 3(f), 3(i), respectively, to perform very well in terms of their focusing and scanning capabilities. In addition, lensing with continuously varying focal lengths can also be achieved. For example, Fig. S6 in the Supplementary Material shows experimentally measured focusing at (0, 0, 450 mm), (0, 0, 600 mm), and (0, 0, 750 mm), respectively. These high-quality and compact focal spots experimentally verify the accuracy of the phase realization of our proposed metasurface.

Electromagnetic fields with a phase profile eilφ carry orbital angular momentum (OAM), where φ is the azimuthal coordinate in the transverse plane and l is the topological charge. To corroborate the generality of the proposed metasurface, we show reprogrammable focused vortex generations in Fig. 4. These vortices are designed to coaxially focus at (0, 0, 600 mm), thus their target rotation angle distribution is the one in Fig. 3(a) plus l/2 times the azimuthal angle. Figures 4(a), 4(e), 4(i), and 4(m) show the required rotation distribution across the surface to generate focused vortex beams of topological charges l=1,2,3,4, respectively. Figures 4(b), 4(f), 4(j), and 4(n) show the measured electric intensity (|Rrr|2) distributions of the vortex beams, respectively, at the focal plane z=600mm. As expected, these intensity distributions exhibit doughnut shapes, and the central dark area of the vortices is larger for increased topological charge. Figures 4(c), 4(g), 4(k), and 4(o) show the corresponding measured phase distributions, where the azimuthal angle dependence clearly reveals the vortex nature of the focused beams. To quantitatively analyze the quality of the vortex beams, we extracted and integrated the complex amplitudes to get the amplitude |Sn| of each OAM (see Supplementary Material 7 for details). Figures 4(d), 4(h), 4(l), and 4(p) show the corresponding normalized |Sn| as a function of l from −6 to 6, respectively. Clearly, the target OAM component of each measurement is the strongest, whereas the other OAM components are quite weak, indicating the high purity of the generated vortex beams. The existence of unwanted OAM beams can be attributed to the square shape of the metasurface and the experimental errors (see Supplementary Material 8 for details). These results experimentally illustrate the large reprogrammability and at the same time high efficiency of our metasurface in tailoring the impinging fields at will.

Metasurface holography has become a well-established approach to address inverse engineering problems for electromagnetic waves. To further explore the wavefront control capability of our programmable metasurface, we apply it to construct computer-generated holographic images in Fig. 5. In our demonstration, the Chinese characters “天津大学” (Tianjin University) and “大同云冈” (Datong Yungang) were generated in the image plane at z=600mm. The required rotation profile for each word, as shown in Figs. 5(a)–5(d) and 5(i)–5(l), was calculated by a modified Gerchberg–Saxton algorithm,64^,65 in which the conventionally used Fresnel diffraction formula is replaced by the Rayleigh–Sommerfeld diffraction formula (see Supplementary Material 9 for details). To confirm the functionality, we calculated the holography imaging by considering the meta-atoms as point sources with ideally designated hologram phase profiles. The calculated results are shown in Fig. S9 in the Supplementary Materials, which comply well with the target characters. To experimentally demonstrate these holographic images, we rotated the meta-atoms according to the required rotation profile and measured the corresponding electric intensity (|Rrr|2) distribution. Although the holograms have only 20×20 pixelated phase points, the generated holographic images, as shown in Figs. 5(e)–5(h) and 5(m)–5(p), are of very high-quality and in good agreement with calculations, indicating the excellent wavefront control capability of the proposed reprogrammable metasurface platform.

Polarization sensitivity is an important feature of metasurfaces. Since the proposed metasurface is based on the PB phase concept, once all the meta-atoms are arranged to achieve the desired phase profile Φ for specific functionality under LCP (RCP) incidence, the corresponding phase profile under RCP (LCP) incidence turns to Φ′=−Φ. In most cases, the phase profile of Φ′ cannot maintain the same functionality as that of Φ. For instance, Fig. S10 in the Supplementary Material shows the measurements in the case of the metasurface are designated to generate a holographic image of the word “PB” under an LCP incidence, where the measured |Rll|2 successfully forms the holographic image while the measured |Rrr|2 is very noisy. Similarly, the demonstrated functionalities shown in Figs. 3–5 are also polarization-sensitive and limited to work for RCP. Inspired by a recent work that utilized binary geometric phase control to achieve polarization-insensitive metalensing in the visible range,66 we find that the proposed metasurface can also maintain the same functionality for multipolarizations. For instance, Fig. S11 in the Supplementary Material shows the rotation profile of meta-atoms to achieve metalensing at (0, 0, 600 mm), where the orientation angles are binarily arranged with −π/4 and π/4. In this condition, the reflected intensity profiles of |Rrr|2,|Rll|2,|Rxy|2, and |Ryx|2 can focus at the same spot, as shown in Fig. S11f-I in the Supplementary Material.

In addition, the polarization-dependence and wavelength-independence of PB phase control can be further exploited to achieve polarization and frequency multiplexed functionalities.67 For achieving functionality under RCP for wavelength λ1, the desired phase distribution can be calculated as ΦR,1. Similarly, for achieving other functionality under LCP for wavelength λ2, the desired phase distribution can be calculated as ΦL,2. To multiplex these phase distributions into one set of meta-atoms, one can arrange the rotation profile as θ=arg[exp(iΦR,1)+exp(−iΦL,2)]/2. In this condition, under RCP incidence the metasurface will perform the phase profiles of ΦR,1 and −ΦL,2, and similarly under LCP incidence the metasurface will perform the phase profiles of −ΦR,1 and ΦL,2. If the target functionalities are metalensings, their corresponding negative phase distributions (−ΦL,2 and −ΦR,1) correspond to divergences and thus have slight influences on the focusing spots achieved by positive phase distributions (ΦR,1 and ΦL,2). For instance, Fig. S12 in the Supplementary Material experimentally shows that the proposed metasurface can focus the |Rrr|2 of 6.3 GHz and |Rll|2 of 7 GHz at different positions. For multiplexing further complicated functionalities, the crosstalk from unwanted parts becomes serious. If only the intensity profile of the target functionality is important (such as metalensing and holography), the phase profiles of ΦR,1 and ΦL,2 can be carefully optimized to reduce the crosstalk.67 Nevertheless, to delink the multiplexed functionalities for RCP and LCP channels, due consideration must be paid into the spin-independent design freedom, such as propagation phase (see Supplementary Material 13 for details).

We have demonstrated a reprogrammable mechanical PB metasurface platform and showcased its superior dynamic and efficient control of the impinging wavefront by reconfiguring in real time its operation across a number of functionalities. From the diffraction optics perspective, the performance of our optical metasurface gradually converges to the ideal response as the phase control level increases. Supplementary Material 14 uses metalensing as an example and numerically compares the performance of 2-level, 4-level, quasicontinuous (28-level), and ideal phase control schemes, showing that the focus intensity of our quasicontinuous scheme is nearly identical to the ideal one, about 2.61 and 1.25 larger than the case of 2- and 4-level phase control schemes, respectively. In the current design, the meta-atoms in the same supercell are designed to have the same orientation angle, where the gear set identically transmits the torque from the stepping motor to each meta-atom. Since the mechanically reprogrammable modules and meta-atoms can be flexibly attached/detached from the metasurface, the proposed metasurface platform can conveniently incorporate different gear sets and different meta-atoms, which may empower multidimensional manipulation of electromagnetic waves with programmability (see Supplementary Material 13 for details).

Modulation frequency is also an important metric of reprogrammable metasurfaces. Among the voltage-driven elements-based metasurfaces, the reported fastest modulation frequency could achieve 10 MHz.54 Microfluidic59 and mechanical61 tuning is generally slow. With our metasurface, the rotation speed of the PB meta-atom is determined by the rotating speed of the step motor, the rotation ratio of the gear set, and the geometric phase dependence coefficient. In the current design, the rotation speed of stepping motors is set as ∼2s per turn, the rotation ratio of the gear set is 8/7, and the PB phase dependence coefficient is ±2, corresponding to a modulation frequency of about 8/7 Hz. It should be emphasized that by adopting step motors with a higher speed, gear sets with a larger rotation ratio, and meta-atoms with high-fold rotational symmetries,63 the PB phase modulation frequency can be much improved. Such speed, although it is still slower than what can be achieved using voltage-driven diodes, can meet the requirement of some important applications, such as dynamic optimization of wireless communication channels,47^,49^,51 intelligent imager/recognizer,40^,45 and self-adaptive cloak,46 since they aim to reconfigure functionality according to variations of weather, temperature, and the motion of pedestrian or stealth vehicles. On the other hand, the proposed mechanical metasurface features low power dissipation (lower than 37.5 milliwatts at the meta-atom level) and nonvolatility, and therefore it is energy-saving and environment-friendly for large-scale and long-term applications.

For the applications that directly launch designated microwave beams, say the military radar system that steers a beam quickly across the sky to detect planes, the well-established phased arrays are more suitable. In contrast, the reprogrammable metasurfaces, either diodes-based metasurfaces or the proposed mechanical metasurface, aim to interact with incident beams and then tailor the out beams, playing the role as a part of the communication channels or surroundings. Therefore, from the function aspect, reprogrammable metasurfaces are more suitable for applications that are designed to control the wavefront of an external source of radiation, such as dynamic optimization of wireless communication channels47^,49^,51 or self-adaptive cloak.46 With our metasurface, the low power dissipation and the nonvolatility are also key for the cloak: the low power dissipation leads to low heat generation, which reduces the risk of detection by infrared technology; the nonvolatility makes the stealth vehicles possible for being hidden a long time and on standby with a relative low energy consumption. It also should be emphasized that each kind of technology has different performances in multiple dimensions. Such diversity ensures each technology has its own place in future applications, and on the other end will encourage continuous endeavors to optimize each specific technology for most suitable scenario.

In summary, we have demonstrated a versatile approach to implement a mechanically reprogrammable PB metasurface with superior performance. In principle, this metasurface can be reconfigured to yield various functionalities of choice. High-performance functionalities can be dynamically switched by this cost-effective metasurface, which is of great interest in various microwave applications. Precise and reversible control of micro/nanoelements in space and time has long been pursued in the context of micro/nanotechnologies.68^–71 We envision that our work may further boost endeavors into the interdisciplinary fields between micro/nanotechnologies to push these concepts to higher frequencies, enabling precise rotation control of micro/nanoscale PB meta-atoms with full addressability and programmability to compress optical reprogrammable metasurfaces into a single chip.

In experiments, a VNA was used to acquire the response data by measuring the transmission coefficients (S21). During measurements, the feed horn antenna and receive antenna/probe are connected to the output and input ports, respectively, of the VNA through two 3-m-long 50Ω coaxial cables. The excitation source is a linearly polarized wave emitted by the feeding antenna with a working bandwidth from 2 to 18 GHz. The experimental results shown in Fig. 2 are obtained by adopting a broadband antenna (the same as the feeding antenna) as the receiver. By respectively rotating the feed and receive antenna, the reflection coefficients for different linearly polarized waves, i.e., Rxx,Rxy,Ryx, and Ryy, can be obtained, then the reflection coefficients for circularly polarized waves can be calculated using the transformation (RrrRrlRlrRll)=12(Rxx−Ryy+i(Rxy+Ryx)Rxx+Ryy−i(Rxy−Ryx)Rxx+Ryy+i(Rxy−Ryx)Rxx−Ryy−i(Rxy+Ryx)).(1)

The experimental results shown in Figs. 3–5 are obtained by adopting a waveguide probe as the receiver. Driven by an electric-controlled translation stage, the waveguide probe raster scans the linearly polarized waves over a large area, and the field distributions of (|Rrr|2) are also calculated by Eq. (1). For reference, Fig. S17 in the Supplementary Material shows the photos and schematics of the experimental setups.

Quan Xu received his BEng degree in communication engineering and his MEng degree in optical engineering from Harbin Engineering University, China, in 2011 and 2014, respectively, and his PhD in optical engineering from Tianjin University, China, in 2019. He is currently an assistant professor at the Center for Terahertz Waves, Tianjin University, Tianjin, China. His research interests focus on metamaterial/metasurface devices, surface plasmons, and terahertz photonics.

Xiaoqiang Su received his BEng degree in physics from Yuncheng College, China, in 2009, his MEng degree in optics from Zhejiang Normal University, China, in 2012, and his PhD in optical engineering from Tianjin University, China, in 2016. He is currently an associate professor at the Institute of Solid State Physics, Shanxi Datong University, Datong, China. His research interests focus on reconfigurable intelligence surfaces and spoof surface plasmons.

Xueqian Zhang is an associate professor at Tianjin University, China. He received his BEng degree in electronic science and technology and his PhD in optical engineering from Tianjin University, China, in 2010 and 2016, respectively. He was a visiting scholar at the University of Birmingham, United Kingdom, and at King Abdullah University of Science and Technology, Saudi Arabia, from 2013 and 2014. His main research interests are terahertz metasurfaces, surface plasmonic waves, and nonlinear plasmonics.

Lijuan Dong received her BSc degree in physics from Yanbei Normal College, China, in 2000, her MSc degree in physics from Hebei University of Technology, China, in 2003, and her PhD in physics from Tongji University, China, in 2009. She is a full professor at the Institute of Solid State Physics, Shanxi Datong University, Datong, China. Her research interests focus on functional electromagnetic devices and wireless energy transfer technology.

Lifeng Liu received his BEng degree in physics from Shanxi Datong University, China, in 2006, his MEng degree in condensed matter physics from Chongqing University, China, in 2009, and his PhD in optical engineering from Chongqing University, China, in 2012. He is currently an associate professor at the Institute of Solid State Physics, Shanxi Datong University, Datong, China. His research interests focus on active metamaterial devices and spoof surface plasmon polaritons.

Yunlong Shi is currently a full professor at Shanxi Datong University. He received his PhD in condensed matter physics from Tongji University, China, in 1999. He was a visiting scholar at the Cavendish Laboratory of Cambridge University from 2003 to 2004. His research interests include photonic crystal, strongly correlated systems, and metamaterials.

Qiu Wang received his BEng degree in optoelectronic science and technology and his PhD in optical engineering from Tianjin University, China, in 2014 and 2020, respectively. He is currently an associate professor at the School of Information Engineering, Wuhan University of Technology, Wuhan, China. His current research mainly includes applications of the metasurfaces and optical calculation algorithms.

Jiaguang Han received his BSc degree in material physics from Beijing Normal University, China, in 2000, and his PhD in applied physics from Shanghai Institute of Applied Physics, Chinese Academy of Sciences, China, in 2006. From 2006 to 2007, he was a postdoctoral researcher at Oklahoma State University, United States. In 2007, he joined the National University of Singapore, Singapore, as a Lee Kuan Yew research fellow. He is a full professor and his research focuses on surface plasmon polaritons, metamaterials, and material studies in the terahertz regime.

Weili Zhang is a full professor at Oklahoma State University. He received his PhD in optical engineering from Tianjin University, China, in 1993. He was a postdoctoral research associate at the Hong Kong University of Science and Technology from 1993 to 1995. His research interests include terahertz optoelectronics, nano- and microstructured materials optics, and ultrafast phenomena. He is a fellow of OSA.

Biographies of the authors are not available.