https://www.nature.com/articles/s41467-023-40717-9 Skip to main content Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript. Advertisement Advertisement Nature Communications * View all journals * Search * Log in * Explore content * About the journal * Publish with us * Sign up for alerts * RSS feed 1. nature 2. nature communications 3. articles 4. article A universal metasurface antenna to manipulate all fundamental characteristics of electromagnetic waves Download PDF Download PDF * Article * Open access * Published: 24 August 2023 A universal metasurface antenna to manipulate all fundamental characteristics of electromagnetic waves * Geng-Bo Wu ORCID: orcid.org/0000-0002-8226-1663^1,2^ na1, * Jun Yan Dai ORCID: orcid.org/0000-0003-1589-1232^3,4,5^ na1, * Kam Man Shum ORCID: orcid.org/0000-0002-8750-8581^1, * Ka Fai Chan ORCID: orcid.org/0000-0002-5719-6855^1, * Qiang Cheng ORCID: orcid.org/0000-0002-2442-8357^3,4,5, * Tie Jun Cui ORCID: orcid.org/0000-0002-5862-1497^3,4,5 & * ... * Chi Hou Chan ORCID: orcid.org/0000-0002-9718-0831^1,2,6 Show authors Nature Communications volume 14, Article number: 5155 (2023) Cite this article * 10k Accesses * 49 Altmetric * Metrics details Subjects * Electronic devices * Metamaterials * Nonlinear optics * Optical manipulation and tweezers Abstract Metasurfaces have promising potential to revolutionize a variety of photonic and electronic device technologies. However, metasurfaces that can simultaneously and independently control all electromagnetics (EM) waves' properties, including amplitude, phase, frequency, polarization, and momentum, with high integrability and programmability, are challenging and have not been successfully attempted. Here, we propose and demonstrate a microwave universal metasurface antenna (UMA) capable of dynamically, simultaneously, independently, and precisely manipulating all the constitutive properties of EM waves in a software-defined manner. Our UMA further facilitates the spatial- and time-varying wave properties, leading to more complicated waveform generation, beamforming, and direct information manipulations. In particular, the UMA can directly generate the modulated waveforms carrying digital information that can fundamentally simplify the architecture of information transmitter systems. The proposed UMA with unparalleled EM wave and information manipulation capabilities will spark a surge of applications from next-generation wireless systems, cognitive sensing, and imaging to quantum optics and quantum information science. Similar content being viewed by others [41467_2023] Pulse-driven self-reconfigurable meta-antennas Article Open access 06 February 2023 Daiju Ushikoshi, Riku Higashiura, ... Hiroki Wakatsuchi [41928_2022] Sideband-free space-time-coding metasurface antennas Article 03 November 2022 Geng-Bo Wu, Jun Yan Dai, ... Chi Hou Chan [41467_2023] Arbitrary aperture synthesis with nonlocal leaky-wave metasurface antennas Article Open access 20 July 2023 Gengyu Xu, Adam Overvig, ... Andrea Alu Introduction Electromagnetic (EM) waves, ranging from microwave and terahertz waves to visible light, are the bases of various disciplines ranging from optics^1, telecommunications^2, and material engineering^3,4 to quantum systems^5. Plane-wave mode with the electric field of \ ({{{{{\boldsymbol{E}}}}}}(\mathop{r}\limits^{ \rightharpoonup },t)=\ hat{e}A\cos (2\pi {ft}+\mathop{k}\limits^{ \rightharpoonup }\cdot \ mathop{r}\limits^{ \rightharpoonup }+\varphi )\) is the eigensolution of Maxwell's equation in free space^6,7. Classical EM waves can be completely characterized by five intrinsic properties: polarization \ (\hat{e}\), amplitude A, frequency f, momentum k, and initial phase ph . These constitutive properties can also be the function of position \(\mathop{r}\limits^{ \rightharpoonup }\) and time t, leading to a higher-dimensional structured wave. Wave science and technology development is intimately related to flexibly controlling and fully utilizing these fundamental properties. Conventional EM wave manipulations primarily rely on the accumulated propagation effect in naturally existing dielectric materials such as lenses, optical modulators, and waveplates^8,9. These optical components are bulky and generally have curved shapes, making them unsuitable for modern integrated electronic and photonic systems. The advance of metasurfaces^10,11, allowing wave-matter interactions within an ultrathin artificial surface, provides a paradigm shift for EM-wave manipulations^12,13,14,15,16,17,18. Since then, the science and engineering communities have long sought-after a dream 'universal' metasurface that enables simultaneous and independent controls over all the fundamental properties of EM waves. Although some effects have been made (Table 1), a single metasurface component that can manipulate all EM waves' properties remains elusive due to the following challenges. Firstly, most metasurfaces are passive, whose functionalities are set in stone and cannot be altered post-fabrication. However, many modern wave-empowered applications such as communications, holographic displays, and light detection and ranging demand dynamic and active controls for environment adapting and/or information processing. Some tunable metasurfaces integrated with active functional materials have been explored to achieve dynamic wave controls upon the external stimuli, including the electrical bias^19,20,21, mechanical deformation^22,23,24, optical pumping^25,26,27,28, and thermal excitation^29,30,31. However, most of the tunable metasurfaces hitherto can control only one or two wave properties due to insufficient degrees of freedom in the element's geometrical parameters and external control variables to support the regulation of all wave properties^32. Another grand challenge is the difficulty of independent wave property manipulations, as the controls over these properties are generally coupled with each other. For instance, unique geometric structures and co-optimization with active materials should be adopted to decouple the amplitude and phase regulations^33,34,35, but the design complexity and insertion loss are amplified geometrically with the increase of the control degrees of freedom. Table 1 Comparison of state-of-the-art metasurfaces in the literature Full size table In recent years, spatiotemporally modulated metasurfaces have become an emerging technique to engineer EM waves in both space and time. Spatiotemporally modulated metasurfaces have introduced an additional dimension, time, into the conventional metasurface design, enabling intriguing physical phenomena^36,37,38,39,40,41 and wave manipulations in frequency-momentum spaces^42,43,44,45,46,47,48,49,50 . In particular, the waveguide-integrated metasurface antenna allows high-efficiency frequency controls without sideband pollution^44,51, 52. However, only the frequency and momentum controls are validated for most spatiotemporally modulated metasurfaces (Table 1). Whether it is possible to control all of the fundamental properties of EM waves with a single radiation aperture remains an outstanding question. The situation is far more challenging to realize simultaneously and independently programmable radiation characteristics. In this article, we report a microwave universal metasurface antenna (UMA) capable of dynamically, simultaneously, independently, and precisely manipulating all fundamental properties of EM waves, including amplitude, phase, polarization, frequency, and momentum. This is enabled by the additional degree of freedom provided by the spatiotemporal modulation and the unique waveguide-integrated structure of the UMA. We show that our UMA facilitates space-varying wave properties for more complicated wave manipulation by demonstrating the Airy beam and focusing beam generation. Our UMA can also enable time-varying wave properties for information manipulation. For validation, we build a wireless communications link in which our UMA directly transmits the modulated waveforms with different modulation schemes. We further reveal our UMA's unique direction information modulation phenomenon, making it an ideal candidate for eavesdropper-proof communications. The complicated wave manipulations and information modulations are achieved via spatiotemporally switching the meta-atoms' ON-OFF (1-bit) coding states. Our UMA offers a universal EM platform for wave and information manipulations, making it a promising enabler for next-generation wireless communications, integrated photonics, and quantum information science. Results Full manipulations of all fundamental EM-wave properties Our UMA consists of an array of subwavelength anisotropic meta-atoms on top of a waveguiding structure (Fig. 1a). In this work, we consider the metasurface antenna operating at microwave frequencies and utilize the positive-intrinsic-negative (PIN) diode and substrate-integrated waveguide (SIW) as the active element and waveguide, respectively (Fig. 1b). Each meta-atom consists of two +-45 ^o-inclined slot openings to radiate two orthogonal eigen-polarization states |u> (+45^o linear polarization) and |v> (-45^o linear polarization) in free space. Each slot opening can be independently switched between the radiating ('1') and non-radiating ('0') states in real time by the PIN diode. Assume that the radiating state of the meta-atom is temporally modulated with a time cycle T[M] = 1/f[M], and the frequency of the injected monochromatic wave is f [0], subject to f[0] [?] f[M]. Applying independent time-coding sequences to all meta-atoms forms a two-dimensional (2D) "0/1" space-time-coding (STC) matrix^42,43,44 (see Supplementary Fig. S1b), which is controlled by a field programmable gate array (FPGA). Our STC metasurface antenna enables extracting and converting the in-plane guided wave (GW) into the out-of-plane propagating wave (PW) with arbitrary wave properties. Fig. 1: The UMA for independent controls of all fundamental wave properties. figure 1 a Schematic of the UMA. Slot-opening meta-atoms are positioned on top of the waveguide to convert the GW into PW with software-defined properties. By applying 1-bit '0/1' space-time-coding sequences to switch the meta-atom between the radiating ('1') and non-radiating ('0') states, we can control all fundamental properties of the radiated wave. Only the momentum of the target harmonic frequency (m = + 1 harmonic in this case) matches with that of free space. In contrast, other unwanted higher-order harmonics are highly suppressed in both the waveguide and free space without phase-matching. b Configuration of the anisotropic meta-atom, consisting of two +- 45^ o-inclined slot openings loaded with PIN diodes. c-g Independent control of wave properties for the m = +1 harmonic frequency radiation. c Measured frequency shifting of the extracted free-space PW with different modulation frequencies f[M]. The momentum control ( d), phase control (e), amplitude control (f), and polarization control (g) are achieved by tuning the time gradient \(\partial {t}_ {i}(x)/\partial x\), reference time shift t[i=1], duty cycle t, and time gradient and reference time shift for the |u> , |v> polarizations, respectively. h-l Independent control of wave properties for the m = 0 fundamental frequency radiation. h Measured spectrum of the radiated PW for m = 0 fundamental frequency radiation. The momentum control (i), phase control (j), amplitude control (k), and polarization control (l) are achieved by tuning the equivalent spatial period L, space shift [?]x of the amplitude envelope, modulation depth M, and spatial period and space shift for the |u> , |v> polarizations, respectively. The input frequency and modulation frequency in (d-l) are f[0] = 23.5 GHz and f[M] = 1.2 MHz, respectively. Here, CW continuous wave, SIW substrate-integrated waveguide, PIN positive-intrinsic-negative. Full size image We first demonstrate the independent controllability of our UMA over all of the fundamental wave properties. Frequency manipulation is generally challenging due to the requirement to change a photon's energy. The typical approach based on nonlinear bulk media suffers from the weak nonlinear effect and stringent phase-matching conditions. Here, the functions of our STC meta-atom are twofold: (i) temporal modulation for producing nonlinear effects; and (ii) simultaneous space-time modulations for converting the newly generated waves into free space to mitigate the phase matching condition (Fig. 1a). Specifically, the periodic switching of the meta-atoms' radiation states leads to an infinite number of harmonic frequencies f[0] + mf[M], where m is an arbitrary integer. We apply an identical rectangular time sequence yet with a position-dependent time shift^42,45 t[i](x) to the meta-atoms (Supplementary Fig. S1c). According to the time-shifting property of the Fourier Transform, a shift in the time domain t[i] corresponds to a linear phase shift -2 pmf[M]t[i] in the frequency domain. Therefore, the total phase shift of the radiated PW consists of two parts (see Supplementary Fig. S1d ): (i) the phase accumulation from the propagation of GW ph[GW] = -x[ GW]x[i] (where x[GW] is the wavenumber inside the waveguide); and (ii) the abrupt phase shift induced by the spatiotemporal modulation ph[ST] = -2pmf[M]t[i](x). The corresponding linear momentum of the radiated PW along the x direction is k[x](f[0] + mf[M]) = x[GW] + k[ ST], where k[ST] = 2pmf[M] [?]t[i](x)/[?]x is the additional momentum imparted by the spatiotemporal modulation. In this illustrative example (Fig. 1a), we leverage this momentum to compensate for the momentum mismatch between the waveguide and free space at the target m = +1 harmonic frequency, i.e., \({-{\xi }_{m=+ 1} \, < \, \xi }_ {{GW}}+{k}_{{ST}} \, < \, {\xi }_{m=+ 1}\), where x[m=+1] is the free-space wavenumber at the m = +1 harmonic frequency. Other unwanted harmonics are not supported and are highly suppressed in both free space and the waveguide due to the tremendous momentum mismatch (see Supplementary Fig. S2 for the measured spectrum distribution). In this manner, our UMA can simultaneously achieve a nearly perfect GW-to-PW conversion and frequency shifting. Figure 1c presents the extracted free-space wave's measured spectra with different frequency values shifting from -1.8 MHz to +1.8 MHz through changing the modulation frequency f[M] by the FPGA. Next, we show the momentum control of the up-converted PW (Fig. 1d). The linear momentum imparted by the spatiotemporal modulation \({k}_ {{ST}}=2\pi m{f}_{M}\partial {t}_{i}(x)/\partial x\) is proportional to the applied time gradient [?]t[i](x)/[?]x. Supposing [?]t[i](x)/[?]x is a constant, the radiated PW has a well-defined radiation angle \({\ theta }_{r}={\sin }^{-1}({k}_{x}/{\xi }_{m})\) (see Methods). Therefore, the momentum and the corresponding output angle of the radiated PW can be easily tuned by changing the applied time gradient ^42 (Fig. 1d). Moreover, since the abrupt phase shift imparted by the spatiotemporal modulation satisfies \({\varphi }_{{ST}}=-2\pi m{t}_ {i}/{T}_{M}\), the initial phase of the extracted PW can be tuned from 0^o to 360^o by altering the reference time shift t[i=1] (the time shift for the 1st meta-atom) while fixing the time gradient [?]t[i ](x)/[?]x (Fig. 1e). Furthermore, the harmonic frequencies' power distribution depends on the time sequence's coding context (see Methods). The amplitude of the extracted PW can be tuned by varying the duty cycle t of the rectangular time sequence for all meta-atoms (Fig. 1f). Finally, it is known that arbitrary polarizations can be decomposed into a linear combination of two complete orthogonal polarization bases, e.g., |u> and |v> , and vice versa. The anisotropic meta-atom consists of a pair of +-45^o-inclined elliptical slot openings with a large length-to-width ratio, whose radiated electric field polarization is perpendicular to the long side (Fig. 1b and Supplementary Fig. S1a). Applying independent STC matrixes to the +-45^o-inclined slots in each meta-atom allows independent control over the amplitude and phase contents of the extracted |u> and |v> components, resulting in arbitrary polarization generation by our UMA. Specifically, an identical time gradient [?]t[i](x)/[?]x is adopted for both the +-45^o-inclined slot openings such that the |u> and |v> components share the same momentum and output direction. We apply different initial time delays and duty cycles to the +-45^o-inclined slot openings to control the amplitude ratio and phase difference between the radiated |u> and |v> components, respectively. Six representative polarizations, including |x> , |y> , |u> , |v> , |LCP > (left-hand circular polarization), and |RCP> (right-hand circular polarization) are given as illustrative examples to show the polarization controllability (see Fig. 1g and Supplementary Table S1 ). The coupling effects between the two +-45^o-inclined slot openings for polarization control are given in Supplementary Note 5. The required STC matrixes for different polarizations are given in Supplementary Fig. S9a. Our UMA can also independently manipulate all the properties of the radiated PW without shifting the frequency (Fig. 1h-l and Methods). We leverage the time-average effect of the spatiotemporal modulation to form an equivalent sinusoidal amplitude distribution at the fundamental frequency (see Methods). In this manner, the n = -1 space harmonic becomes a fast wave and converts to a PW in free space (Fig. 1h). The momentum and the corresponding output angle of the extracted PW can be tuned by varying the spatial period L through changing the applied STC matrix (Fig. 1i). Next, we utilize the space-shifting property of the Fourier Transform to control the phase of the extracted PW. In analogy with the time-shifting property but in the space domain, a shift in the space domain [?]x results in a phase shift 2pn[?]x/L for the nth-order space harmonic. We apply different space shifts [?]x to the sinusoidal amplitude envelope, resulting in different phases to the radiated n = -1 space harmonic (Fig. 1j). The extracted power for each meta-atom from the waveguide is proportional to the modulation efficiency of the sinusoidal amplitude modulation \(\eta=\frac{{(\frac{M}{2})}^{2}}{1+{2(\frac{M} {2})}^{2}}=\frac{{M}^{2}}{4+2{M}^{2}}\). Therefore, the amplitude of the extracted PW of the UMA can be tuned by changing the modulation depth M (Fig. 1k). Finally, our UMA can generate arbitrary polarization by controlling the amplitude ratio and phase difference of the extracted |u> and |v> components. To this end, we apply equivalent sinusoidal amplitude distributions with an identical spatial period L to both the +-45^o-inclined slot openings such that the extracted |u> and |v> components share the same output angle in free space. We apply different space translations and modulation depths to the +-45^o-inclined slot openings to control the amplitude ratio and phase difference of the |u> and |v> components, respectively. As proof-of-concept examples, the UMA is loaded with different STC matrixes (see Supplementary Fig. S9b) to generate six representative polarizations, including |x> , |y> , |u> , |v> , |LCP > , and |RCP> (Fig. 1l and Supplementary Table S1). We further show the simultaneous controls over EM wave's properties by investigating the far-field radiation patterns of the UMA. In each panel in Fig. 2, three radiating cases with different magnitudes are given to demonstrate the independent amplitude controllability. The panels in the same row in Fig. 2 represent the radiated PWs with identical wave properties except with different linear momentums, whose corresponding beam angle steers from -40^o, 0^o to +40^o, verifying the simultaneous amplitude and momentum controls. We change the polarization from |x> in the first-row panels to |u> in the second-row panels, leading to simultaneous amplitude, momentum, and polarization controls. We further change the output frequency from the fundamental frequency (m = 0) in the 2nd-row panels to m = +1 harmonic frequency in the 3rd-row panels while fixing all other wave properties. In this case, simultaneous amplitude, momentum, polarization, and frequency controls have been verified. Fig. 2: The UMA for simultaneous manipulations of the radiated EM wave properties. figure 2 a-c Measured far-filed radiation patterns of the UMA at the fundamental frequency with the |x> polarization, whose main beam scans to -40^o (a), 0^o (b) and 40^o (c), respectively. d-f Same as ( a-c), except that the polarization is changed to |u> . g-i Same as (d -f), except that the frequency of the extracted PW is changed to the m = +1 harmonic. Three curves in each panel demonstrate the amplitude controls of the radiated PWs. The input frequency and modulation frequency are f[0] = 23.5 GHz and f[M] = 1.2 MHz, respectively. Full size image Space-varying wave properties for complicated wave manipulation We have demonstrated the full-dimensional tailoring of the EM waves yet with spatial-invariant (uniform) wave properties. Here, we show that our UMA can generate more complicated EM waves with space-varying wave properties for more complicated waveform generation and beamforming. As the first illustrative example, our metasurface antenna can generate an Airy beam with phase-varying properties (Fig. 3a). The required phase profile along the metasurface aperture should fulfill the following equation to generate a parabolic trajectory x [?] z^2 ^53 $$\varphi \left(x\right)=\frac{4}{3}{a}^{\frac{1}{2}}{\xi }_{m}{(-x)} ^{3/2}$$ (1) where a is the acceleration factor and x[m] is the free-space wavenumber at the target mth-order harmonic frequency. For the UMA, the phase distribution of the extracted wave at the mth-order harmonic frequency is the sum of the accumulated phase shift of the guided wave and the equivalent phase shift imparted by the spatiotemporal modulation, given by \({\varphi }_{{GW}}+{\varphi }_ {{ST}}={-\xi }_{{gw}}x-2\pi m{f}_{M}{t}_{i}(x)\) (Supplementary Fig. S1d). Combining with Eq. (1), the required normalized time shift of the UMA is given by $$\frac{{t}_{i}\left(x\right)}{{T}_{M}}=\frac{\frac{4}{3}{a}^{\frac {1}{2}}_{\xi_{m}}{(-x)}^{3/2}+{\xi }_{{gw}}x}{-2\pi m}$$ (2) Fig. 3: The UMA for complicated beam shaping with space-varying wave properties. figure 3 a-c UAM for Airy beam generation. a Schematic illustration of the UMA for generating the Airy beam at the m = -1 harmonic frequency in free space. The Airy beam requires a spatial-varying phase profile as \(\ frac{4}{3}{a}^{1/2}{\xi }_{m}{(-x)}^{3/2}\). b The required aperture phase profile of the Airy wavefront and the corresponding normalized time shifting for the meta-atoms at different positions. c The measured E-field intensities in the xz-plane at different harmonic frequencies. The green dashed line represents the theoretical parabolic trajectory of the Airy beam x = az^2, where the acceleration factor a = 0.003 in this case. Only the target m = -1 harmonic frequency possesses high field intensity, whereas other harmonic frequencies are all highly suppressed in free space. d-h UAM for focused beam generation. d Schematic illustration of the UMA for wave focusing at the m = -1 harmonic frequency. In this case, the local momentum or output angle of the extracted PWs from all meta-atoms are different. e The required normalized local momentum k[ x]/x[m] as a function of the meta-atom position for different focal spot positions from F[1] to F[5], where F[1] = (-60, 0, 60) mm, F[2] = (-30, 0, 60) mm, F[3] = (0, 0, 60) mm, F[4] = (30, 0, 60) mm, and F [5] = (60, 0, 60) mm. f Measured cross-section intensity distributions at the m = -1 harmonic frequency in different transversal planes for the designed focal point F[3] (0, 0, 60) mm. g The measured E-field intensities in the xz-plane at the m = -1 harmonic frequency for different focus spots. h The corresponding 1-D normalized intensity distributions on the focal plane z = 60 mm. The input frequency and modulation frequency are f[0] = 23.5 GHz and f[M] = 1.2 MHz, respectively. Full size image The calculated time shift for generating the Airy beam with the acceleration factor a = 0.003 at the m = -1 harmonic frequency is given in Fig. 3b. Figure 3c presents the measured electric field intensity distributions at different harmonic frequencies. We observe that the UMA can extract and mould the waves into an accelerating and non-diffractive Airy beam with a well-defined parabolic trajectory. Moreover, other undesired harmonic frequencies are highly suppressed in free space. Our metasurface antenna can also mould the extracted waves with space-varying momentum properties k[x](x) for wave-focusing applications (Fig. 3d). To focus the extracted wave into a desired focal point F = (x[F], z[F]), the linear momentum of the extracted PW should fulfill the following formula based on the geometric ray and the geometrical relationship (Supplementary Fig. S4a) $${k}_{x}(x)={\xi }_{m}\frac{{x}_{F}-x}{\sqrt{{\left({x}_{F}-x\ right)}^{2}+{({z}_{F}-z)}^{2}}}$$ (3) The corresponding required time gradient to achieve such momentum of the extracted PW is $$\partial {t}_{i}(x)/\partial x=\left[{\xi }_{m}\frac{{x}_{F}-x}{\ sqrt{{\left({x}_{F}-x\right)}^{2}+{({z}_{F}-z)}^{2}}}-{\xi }_{{gw}}\ right]/(2\pi m{f}_{M})$$ (4) As illustrative examples, we consider the UMA generates different intended focal points from F[1] to F[5] at the m = -1 harmonic frequency, whose required spatial-varying linear momentums are presented in Fig. 3e. The corresponding 2-D field intensities distributions at the m = -1 harmonic frequency are presented in Fig. 3g, respectively. We observe that the metasurface antenna can extract and mould the PW wave into intended focal spots, whose positions are software-defined according to the applied STC matrix. Figure 3f shows the measured cross-section intensity distribution in different transversal planes for the designed focal point F[3] (0, 0, 60) mm. We can observe that the E-field intensity increases as close to the focal plane. Due to the free-space wave spreading effect, the density peak is located at around z = 50 mm, which is between the designed focal point and antenna's aperture^54. Nevertheless, the strongest field on the interested focal plane z = 60 mm is exactly located at the designed focal position for all the focusing scenarios (Fig. 3h), making it suitable for microwave real-time imaging systems^55. The effects of the 1D and 2D radiating aperture on the fading distance for the Airy beam and focusing beam are investigated in Supplementary Note 6. The wave focusing effect of the UMA is weaker as the designed focus moves away from the radiating aperture, as shown in Supplementary Fig. S15. Our UMA's flexible and agile beam-shaping capability holds promising potential in sensing, imaging, and wireless power transfer applications. Time-varying wave properties for information manipulation We further demonstrate that our UMA enables information manipulation by generating time-varying wave properties. Figure 4a presents the digital baseband information and the corresponding transmitted modulated PWs in free space for some popular modulation schemes in modern wireless communications systems. The amplitude-shift keying (ASK), phase-shift keying (PSK), and quadrature amplitude modulation (QAM) schemes map to the free-space PWs with time-varying amplitude, time-varying phase, and time-varying amplitude and phase properties, respectively. Generation of such modulated waveforms in free space by conventional radiofrequency transmitters relies on a heterodyne architecture consisting of a series of active/passive modules, including digital to analog converters, modulators, mixers, filters, phase shifters, and antenna arrays to mould the PW with time-varying amplitude and phase (information-carried) properties. In contrast, our UMA--a single component only-- can directly generate identical modulated waveforms with time-varying amplitude and phase properties in a single step. This can be achieved by loading the dynamic STC matrices, mapping to the desired digital information stream to be transmitted (Fig. 4a). This provides a radically new communication paradigm in the physical layer with the advantages of a much simpler structure, higher integration, lower cost, and lower power consumption. Figure 4c, d presents the measured decoded constellation diagrams and eye diagrams at the receiver end when the information is carried at the m = +1 and m = 0 harmonic frequencies, respectively. The information manipulation functionality of our UMA is validated for different information modulation schemes, including QPSK, 8PSK, 16 amplitude and phase-shift keying (16APSK), and 16QAM with transmission data rates of 1 Mbps, 1.5 Mbps, 2 Mbps, and 2 Mbps, respectively. Fig. 4: The UMA for information modulation with time-varying wave properties. figure 4 a One-to-one mapping relationship among the transmitted digital information stream, the radiated information-carried waveforms in the free space (2ASK, BPSK and 16QAM), and the applied dynamic STC matrixes for communications at the m = +1 harmonic frequency. The required dynamic STC matrixes for generating modulated waveform at the fundamental frequency are presented in Supplementary Fig. S5b. b The UMA is used to directly generate the information-carried EM waveforms in free space. c, d The measured decoded constellation diagrams and eye diagrams at the receiver end for communications at the m = +1 (c) and m = 0 (d) frequencies, respectively. CW continuous wave, Mod modulated, Ant. antenna, VSA vector signal analyzer, EVM error vector magnitude. Full size image Moreover, the control over other properties of EM waves (momentum, frequency, and polarization) by our UMA opens additional opportunities to achieve space-division multiplexing (SDM), frequency-division multiplexing (FDM), and polarization-division multiplexing (PDM), which can establish multiple independent channels to improve the communications capacity. We set up a joint SDM-FDM-PDM data transmission link (Figs. 5a and 6c), in which our UMA simultaneously transmits two independent channels with different beam directions (-30^o and +30^o for SDM), different polarizations (|u> and |v> for PDM), different frequencies (m = -1 and m = -2 harmonics for FDM), and different information modulation schemes (QPSK and 8PSK). To this end, we divide the UMA into two interwoven sub-metasurfaces (shown in Fig. 5a with two different colors). Each sub-metasurface, loaded with independent dynamic STC matrixes, generates one specific modulated waveform corresponding to one communication channel. Specifically, we apply two different time gradients to the two sub-metasurfaces; the equivalent momentums imparted by the spatiotemporal modulation push the m = -1 and m = -2 harmonic frequencies to radiate into the free space with output angles of -30^o and 30^o for the two sub-metasurfaces, respectively. Moreover, the amplitude and phase contents (carried information) of the two extracted PWs can be independently controlled by the applied duty cycle and reference time shift according to the transmitted binary bit streams of the two channels. The polarization states of the two radiated beams can be controlled by applying different space-time coding sequences to the +-45^o-inclined slot openings in each sub-metasurface. Each sub-metasurface's bit stream is first mapped to the corresponding dynamic STC matrixes according to the desired conversion harmonic frequency, polarization, and signal modulation format. The final dynamic STC matrixes for dual-channel multiplexing is the sum of the two dynamic STC matrixes for the two interwoven sub-metasurfaces (see Supplementary Fig. S6). The lattice of meta-atoms in the sub-metasurface is 0.31l[0] (where l[0] is the free-space wavelength at 23.5 GHz). The sub-metasurface is free of higher-order diffractions in the free space (Fig. 5b). Moreover, we can flexibly change the radiation directions of CH1 and CH2 (Fig. 5c ). This can be achieved by changing the applied time gradients of the two sub-metasurfaces to vary the momentum property of the extracted PWs. The measured radiation patterns (Fig. 5b, c) and the measured decoded constellation diagrams of the two channels at the receiver end (right inset of Fig. 5b, c) verify that our UMA can realize simultaneous and independent multiplexing successfully. Fig. 5: The UMA for multi-channel and physical-layer security communications. figure 5 a The concept illustration of the UMA for directly generating two independently modulated waveforms carrying different digital data streams. Each sub-metasurface is responsible for one modulated waveform generation. Two designated users at different locations can simultaneously and independently receive the information from the UMA. b, c Measured radiation patterns of the UMA at the m = -1 (CH1) and m = -2 (CH2) harmonic frequencies. The insets on the right present the measured decoded constellation diagrams of the two channels at the receiver end. The beam directions for the CH1 and CH2 are (-30^o, 30^o) and (0^o, 30^o) in (b) and (c), respectively. d The concept illustration of the UMA with the inherent direction modulation property at the fundamental frequency. Only the target user in the main beam direction can successfully decode the information, whereas eavesdroppers at other positions totally lose the information. e Measured radiation pattern of the UMA at the fundamental frequency, and the measured decoded constellation diagrams when the receiver is located in different directions. The physical-layer security communication link for other main-beam directions is validated in Supplementary Fig. S8a, b. CW continuous wave, Mod modulated, Ant. antenna, VSA vector signal analyzer, CH channel, EVM error vector magnitude. Full size image Fig. 6: Prototype characterization. figure 6 a, b Photograph and schematic of the measurement setup for radiation pattern measurement of the UMA. c Dual-channel wireless communications link testbed, in which the UMA directly and simultaneously generates two independent modulated waveforms carrying different information. Two horn antennas connected to two vector signal analyzers (VSAs) are used to receive and demodulate the signals. DUT device under test, VNA vector network analyzer, Rx receiver, Sig. Gen. signal generator. Full size image In conventional transmitter architectures, the radiated EM waves in different directions hold identical time-varying wave properties (information). Therefore, an eavesdropper can recover the information even in the sidelobe region by using a sufficiently sensitive receiver (Supplementary Fig. S8c, d). Here, we reveal an intriguing inherent directional modulation (IDM) phenomenon of our UMA at the fundamental frequency. This unique IDM property is attributed to the direction-dependent phase control of our UMA based on spatial modulation (see Methods and Supplementary Fig. S7 for a detailed discussion). For validation, we establish a communication link, where our UMA directly emits the modulated waveforms carrying the 8PSK signals at the fundamental frequency (Fig. 5d). Figure 5e presents the corresponding measured radiation pattern and the decoded constellation diagrams when the receiver is located in different directions. We clearly note that only the receiver in the main-beam direction of the UMA can successfully decode the information, while the eavesdroppers outside the main beam lose the information entirely. This IDM effect can effectively mitigate malicious eavesdropper attacks in different directions, establishing physical-layer security on top of wireless communications. The IDM is an intrinsic property of our UMA and is free of any optimization and performance trade-off generally incurred for conventional approaches to achieve physically secured links^43,56. Discussion In summary, we proposed and demonstrated a UMA capable of dynamically, simultaneously, independently, and precisely manipulating all of the fundamental properties of EM waves. We further demonstrated that our UMA could facilitate complicated wave and information manipulations. The UMA can directly generate the modulated waveforms with dynamic wave properties, leading to a paradigm shift for new information-transmitting architectures. The UMA concept can be extended to a 2-D aperture by periodically repeating the 1-D metasurface antenna along the y-axis, fed by a power-dividing network (See Supplementary Note 4 for the detailed metasurface configuration and 2D wavefront engineering). Our UMA uniquely combines a number of significant merits, including full-dimensional wave controllability, inherent IDM, simplified coding scheme (1-bit), free of sideband pollution, and potential on-chip integration, making it an appealing enabler for the next-generation large-capacity and high-security information systems. Although demonstrated at the microwave band, the proposed concept can be extended to terahertz frequencies based on CMOS chips^21, monolayer molybdenum disulfide switches^57, and high-electron-mobility InAIN/GaN-based metadevices^58. In the optical band, the nanoscale switch can be implemented by phase-changed materials, such as germanium antimony telluride^59,60. These technologies greatly extend the applications of the UMA in augmented reality, holography, integrated sensing and communications for 6 G, quantum optics, and quantum information science. Methods Wave manipulations with frequency shifting Here we extend the theoretical model of the waveguide-integrated metasurface antenna^44 to a more generalized one by taking the polarization into account. In contrast to the continuous field modeling in ref. ^44, we describe the metasurface by an array of subwavelength scatterers with the discrete nature of the metasurface. The UMA consists of an array of slot-opening meta-atoms etched on the top conductive layer of the SIW waveguide along the x-axis (Supplementary Fig. S1a). The anisotropic slot-opening meta-atom is equivalent to two orthogonally orientated waveguide-fed magnetic dipoles to extract energy from the waveguide and radiate to free space. The wave properties of the extracted |u> and |v> components can be independently controlled by applying independent control voltages to the PIN diodes incorporated in the meta-atom (Supplementary Fig. S1b). The excited polarized magnetic dipole moments \({\bar{\bar{m}}}_{i}=\left[\begin{array}{c}{m}_{i}^{u}\\ {m} _{i}^{v}\end{array}\right]\) with |u> and |v> components for the i th meta-atom at the position x[i] is $${\bar{\bar{m}}}_{i}={\bar{\bar{P}}}_{i}H_{i}={\bar{\bar{P}}}_{i}(t) {H}_{0}{e}^{-j{\xi }_{{gw}}{x}_{i}}{e}\,^{j2\pi {f}_{0}t}$$ (5) where \({\bar{\bar{P}}}_{i}=\left[\begin{array}{c}{P}_{i}^{u}\left(t\ right)\\ {P}_{i}^{v}\left(t\right)\end{array}\right]\) is the magnetic polarizability Jones matrix at instant t, and H[i] is the magnetic field of the reference guided wave inside the waveguide. Since the radiating state of the meta-atom is periodically ON-OFF switched with a time cycle T[M], the magnetic polarizability is a periodic function of time, satisfying \({\bar{\bar{P}}}_{i}\left(t\ right)={\bar{\bar{P}}}_{i}(t+{T}_{M})\), which can be decomposed into a Fourier series $${\bar{\bar{P}}}_{i}\left(t\right)=\mathop{\sum }\limits_{m=-{{\ infty }}}^{{{\infty }}}{\bar{\bar{p}}}_{i,m}{e}\,^{j2\pi {{mf}}_{M}t} $$ (6) where f[M] = 1/T[M] is the modulation frequency. The Fourier coefficients \({\bar{\bar{p}}}_{i,m}=\left[\begin{array}{c}{p}_{i,m}^ {u}\\ {p}_{i,m}^{v}\end{array}\right]\) can be calculated by $${\bar{\bar{p}}}_{i,m}=\frac{1}{{T}_{M}}{\int }_{0}^{{T}_{M}}{\bar{\ bar{P}}}_{i}\left(t\right){e}^{-j2\pi m{f}_{M}t}{dt}$$ (7) We substitute Eq. (6) into Eq. (5) to yield: $${\bar{\bar{m}}}_{i}={H}_{0}\mathop{\sum }\limits_{m=-{{\infty }}}^ {{{\infty }}}{e}\,^{j2\pi \left({f}_{0}+m{f}_{M}\right)t} \, {\bar{\ bar{p}}}_{i,m}{e}^{-j{\xi }_{{gw}}{x}_{i}}$$ (8) From Eq. (8), one can obtain that the periodic ON-OFF switching of the meta-atom will generate an infinite number of harmonic frequencies with a frequency interval f[M]. This nonlinear effect opens the opportunity to control the frequency property of EM waves. Moreover, the spatiotemporal modulation introduces an equivalent magnetic polarizability \({\bar{\bar{p}}}_{i,m}\left(\,{f}_{0}+{{mf}} _{M}\right)\) at the mth-order harmonic frequency, providing additional degrees of freedom to control the other properties of EM waves. Once the excited magnetic polarizability of each meta-atom is known, the radiation pattern of the UMA in free space can be obtained (Supplementary Note 1). Firstly, we adopt the most straightforward identical rectangular time sequence (Supplementary Fig. S1c with [?]t = 0) to all meta-atoms. The corresponding magnetic polarizability can be written as $${\bar{\bar{P}}}_{i}\left(t\right)={P}_{0}\left\{\begin{array}{cc}{\ bar{\bar{1}}} & {-{\bar{\bar{\tau}}}}_{i}/2\le {\bar{\bar{t}}}/{T}_ {M}\le {{\bar{\bar{\tau}}}}_{i}/2 \\ {\bar{\bar{0}}} & {others}\end {array}\right.$$ (9) where P[0] is the constant magnetic polarizability as the meta-atom is in the radiating state. \({{\bar{\bar{\tau}}}}_{i}=\left[\begin {array}{c}{\tau }_{i}^{u}\\ {\tau }_{i}^{v}\end{array}\right]\) is the duty cycle for the |u> and |v> components, which is defined as the ratio of the time a meta-atom is in the coupling state ('1') over a modulation time cycle. Substitute Eq. (9) into Eq. (7), the equivalent magnetic polarizability \({\bar{\bar{p}}}_{i,m}\) for the ith meta-atom at the mth-order harmonic frequency reads (See Supplementary Note 2 for the detailed derivation) $${\bar{\bar{p}}}_{i,m}={{P}_{0}{\bar{\bar{\tau}}}}_{i}{sinc}(\pi m {{\bar{\bar{\tau}}}}_{i})$$ (10) From Eq. (10), the amplitude of the magnetic polarizability at each harmonic is a function of the employed duty cycle \({{\bar{\bar{\ tau}}}}_{i}\). To further enable phase control, we introduce a time delay \({{\bar{\bar{t}}}}_{i}=\left[\begin{array}{c}{t}_{i}^{u}\\ {t} _{i}^{v}\end{array}\right]\) to the meta-atom (Supplementary Fig. S1c ), and the equivalent magnetic polarizability \({\bar{\bar{p}}}_{i,m} \) becomes $${\bar{\bar{p}}}_{i,m}={P}_{0}{{\bar{\bar{\tau}}}}_{i}{sinc}(\pi m {{\bar{\bar{\tau}}}}_{i}){e}^{-j2\pi m{f}_{M}{{\bar{\bar{t}}}}_{i}}$$ (11) Comparing Eqs. (10) and (11), an additional phase term \({\varphi }_ {{ST}}=-\!2\pi m{f}_{M}{{\bar{\bar{t}}}}_{i}\), depending on the applied time delay \({{\bar{\bar{t}}}}_{i}\), is introduced to the equivalent magnetic polarizability at the mth-order harmonic frequency (Supplementary Fig. S1d). In this manner, the amplitude and phase of the excited magnetic polarizability can be decoupled and controlled independently. Wave manipulations without frequency shifting For the wave manipulations at the fundamental frequency (m = 0), the equivalent magnetic polarizability due to the spatiotemporal modulation in Eq. (7) can be simplified as $${\bar{\bar{p}}}_{i,m=0}=\frac{1}{{T}_{M}}{\int }_{0}^{{T}_{M}}{\bar {\bar{P}}}_{i}\left(t\right){dt}$$ (12) We can observe that the equivalent magnetic polarizability at the fundamental frequency is the time-average magnetic polarizability for one modulation cycle. Specifically, when a rectangular time sequence is adopted, the equivalent magnetic polarizability at the fundamental frequency in Eq. (11) simplifies as \({\bar{\bar{p}}}_{i,m=0}={{P}_ {0}{\bar{\bar{\tau}}}}_{i}\). The spatiotemporal modulation generates an equivalent amplitude of the magnetic polarizability without introducing any phase shifts or momentum for the fundamental frequency^44. To facilitate GW to PW transformation at the fundamental frequency, we leverage the equivalent magnetic polarizability to form a sinusoidal amplitude modulation along the length of the waveguide such that the n = -1 space harmonic becomes fast and radiates into free space (Supplementary Fig. S3a)^61 $${\bar{\bar{A}}}(x)={A}_{0}\left[1+{\bar{\bar{M}}}{\cos}\left(\frac {2\pi }{\Lambda }x\right)\right]$$ (13) where \({\bar{\bar{M}}}=\left[\begin{array}{c}{M}^{u}\\ {M}^{v}\end {array}\right]\) are the modulation depths for the |u> and |v> components, and L is the spatial period of the sinusoidal amplitude envelope. The radiation of the equivalent magnetic dipole can be viewed as spatially sampling the reference guided wave at each meta-atom position. The momentum of the n = -1 space harmonic along the x direction is k[ x] = x[gw] - 2p/L, which matches that of the free space providing that -1 < k[x]/x[0] < 1. The corresponding output angle of the n = -1 space harmonic is \({\theta }_{r}={\sin }^{-1}(\frac{{\xi }_{{gw}}-2\ pi /\Lambda }{{\xi }_{0}})\)^61. Furthermore, to suppress the higher-order harmonic frequencies, we randomize the applied time delays and the time sequence yet maintain the equivalent sinusoidal modulation to the meta-atoms^44 (Supplementary Fig. S3b). Since none of the higher-order harmonic frequencies' momentums match that of the free space and waveguide, all the undesired higher-order harmonic frequencies (m [?] 0) are highly suppressed. Inherent directional modulation of the UMA The direction-dependent phase property at the fundamental frequency can be elucidated from the perspective of spatial Fourier transform. The aperture field distribution of the metasurface and its spatial frequency spectrum (far-field radiation pattern in free space) F(th) fulfills the Fourier transform relationship \(A\left(x\right){e}^{-j {\xi }_{{gw}}x}\mathop{\longrightarrow}\limits^{FT}F(\theta )\), where A(x) is the equivalent spatial amplitude envelope in Eq. (13) imparted by the spatiotemporal modulation. Now we introduce a space translation [?]x to the spatial amplitude envelope (see Supplementary Fig. S7a). The aperture field distribution becomes \(A\left(x-\Delta x\right){e}^{-j{\xi }_{{gw}}x}\), whose far-field radiation pattern is (see Supplementary Note 3 for the detailed derivation) $${F}^{\prime} \left(\theta \right)=e^{-j({{\xi }_{{gw}}-\xi }_{0}\ sin \theta )\Delta x}F(\theta )$$ (14) We can observe from Eq. (14) that a space translation [?]x of the amplitude envelope introduces an additional phase shift \(\Delta \phi (\theta )=-\!({{\xi }_{{gw}}-\xi }_{0}\sin \theta )\Delta x\) to the phase pattern without affecting the power pattern of the UMA. Moreover, the introduced phase shift Dph is a function of the observation direction th, meaning that different observation directions possess different phase shifts for a fixed space translation [?]x of the amplitude envelope. As a special case, the phase shift in the main beam direction th[r] is \(\Delta \phi \left({\ theta }_{r}\right)=-\frac{2\pi }{\Lambda }\Delta x\), considering the main beam direction is \({\theta }_{r}={\sin }^{-1}\left(\frac{{\xi } _{{gw}}-2\pi /\Lambda }{{\xi }_{0}}\right)\)^61. For clear illustration purposes, here we take the simplest BPSK modulation scheme as an example to show the directional information modulation of our UMA. Supplementary Fig. S7a shows the two equivalent sinusoidal amplitude envelopes to generate high-directivity beams with the same amplitude yet 180^o phase difference at the main-beam direction (th[r] = 0^o), mapping to the digital information '0' and '1', respectively. To this end, the two equivalent sinusoidal amplitude envelopes share an identical spatial period L yet with a space shift [?]x = L/2. The phase difference between the two radiating cases as a function of the observation direction th is presented in Supplementary Fig. S7b based on Eq. (14). We observe that a desired 180^o phase difference is generated at the main beam direction (th[r] = 0^o in this case), while deviating significantly from off-broadside directions. To directly generate BPSK modulated waveform, we apply the dynamic STC matrixes to the UMA according to the desired transmitted '0/1' digital stream (Supplementary Fig. S7c). Supplementary Fig. S7d further shows the theoretical radiation power pattern of the UMA (Supplementary Note 1) and the calculated error vector magnitude (EVM) at different observation directions in a noiseless environment. We observe that eavesdroppers receive a much weaker power at off-main-beam directions. Most importantly, off-angle eavesdroppers sense incorrect phase variation (information) compared to the target user in the main beam direction, leading to corrupted constellations. We further calculate the bit error rates (BERs) for the receivers at different directions to validate the direction modulation of the UMA. The closed-form BER equation for BPSK signal in an additive white Gaussian noise (AWGN) channel can be expressed as^62 $${P}_{b}=Q(\sqrt{2{\gamma }_{b}})$$ (15) where g[b] is the signal-to-noise power ratio (SNR) per bit. Q(t) is the cumulative distribution function of the standard Gaussian random variable $$Q\left(t\right)={\int }_{t}^{{{\infty }}}\frac{1}{\sqrt{2\pi }}{e}^ {-\frac{{x}^{2}}{2}{dx}}$$ (16) The Monte Carlo simulation is widely adopted to calculate the signal BER in digital communications systems. Supplementary Fig. S7e shows the theoretical (black dot) and simulated (black line) BERs as a function of SNR. We can observe that excellent agreement between the two results, verifying the Monte Carlo simulation as an effective approach to calculating the BER. Due to the lack of a closed-form expression in BER calculation for the UMA, here we adopt the Monte Carlo simulation to calculate the BER performance in different directions. In the simulation, the power of the AWGN in different directions is identical. The simulated BERs versus SNR for receivers at 0^o to 50^o are shown in Supplementary Fig. S7e. Moreover, Supplementary Fig. S7f presents the simulated BER as a function of the observation direction when the SNR at the main beam direction is set as 20 dB. We can observe that receivers in the main-beam direction enjoy a much smaller BER than that of receivers at the sidelobe regions. Therefore, eavesdroppers at the sidelobe region have an extremely low probability of intercepting the information. The inherent directional modulation of the UMA is a kind of static direction modulation^63. An eavesdropper may still exploit the correlation between the main beam and sidelobe. However, the BER at the sidelobe region is significantly larger than that of the main beam region. An eavesdropper can easily intercept the information for conventional transmitters without any direction modulation due to the broadcast nature of wireless communications. Our UMA can significantly reduce the probability of successful data interception by eavesdroppers. The IDM property of the UMA is equally applied to other higher-order modulation formats, such as 8PSK and QAM. Prototype design The configurations of the anisotropic meta-atom are illustrated in Fig. 1b. For easy integration with other components, here we adopt the SIW as the waveguiding structure, which uses two parallel rows of metallic via holes and the thin dielectric substrate to realize the rectangular waveguide in planar form. Two +-45^o inclined elliptical slot openings are etched on the top metallic surface of the SIW. The lattice size of the unit cell along the x-direction is 2 mm, corresponding to 0.158l[0]. Each slot opening works as a polarizable magnetic dipole whose extracted electric field polarization is perpendicular to the long side of the slot. The slot meta-atom is intentionally designed at an off-resonance state by carefully tuning the geometric parameters of the slot opening. Four PIN diodes (MACOMMADP-000907-14020x) are placed across the capacitive gaps of the slot openings in each meta-atom (Fig. 1b). A DC bias circuit is integrated into the meta-atom design. The circuit consists of a fan-shaped bias line (for radiofrequency choking) on the bottom biasing circuitry and a control via connecting the top meta-atom to the bottom biasing circuitry. We directly utilize the vias fences that existed in the SIW as the control vias to alleviate the perturbation of the bias network on the guided wave. Two PIN diodes in the same slot opening are biased in the same state, while PIN diodes in different inclined slot openings are biased and controlled independently. The UMA consists of 41 meta-atoms with 164 PIN diodes in this proof-of-concept example. A low-cost FPGA control board (ALTERA Cyclone IV) generates 82 independent control signals (two independent control signals for each anisotropic meta-atom) to control the ON-OFF states of the PIN diodes and the radiating states of the meta-atoms. The proposed UMA is based on the 1-bit ON-OFF spatiotemporal modulation of the meta-atoms. Therefore, the limitation for the EM wave's property control is largely determined by the ON-OFF switch speed of the meta-atom. The adopted PIN diode (MACOMMADP-000907-14020x) can support the fastest switching speed of 2-3 ns. Although the adopted FPGA (ALTERA Cyclone IV) allows a maximum internal clock speed of 1.3 GHz using a phase-locked loop, the relatively low-speed interconnection (Supplementary Fig. S13a) between the FPGA board and metasurface antenna board limits the maximum control speed of the whole metasurface system. Supplementary Fig. S13b-d shows the measured control signal waveforms with different frequencies from the FPGA using an oscilloscope. We can observe that the signal waveform with 50 MHz distorts from a perfect rectangular waveform due to the parasitic inductance of the interconnection. Therefore, the maximum frequency shift for the current metasurface system is around 20 MHz and the transmission data rate of 2 Mbps for information transmission. Nevertheless, the limitation of the control speed of the metasurface system can be largely mitigated by adopting, e.g., (1) high-speed interconnect technologies (such as FMC connector, PCI Express), (2) packaging the FPGA and metasurface on the same PCB board, and (3) using high-speed PCB routing. The radiation characteristics of the meta-atom are modeled and simulated using the commercially available ANSYS HFSS numerical simulator based on the finite element method. The PIN diode was modeled as a series of resistance R = 8 O and inductance L = 30 pH for forward biased (State '0') and a series of capacitance C = 0.052 pF and inductance L = 30 pH for non-bias (State '1'), respectively. Wave ports are adopted to excite the fundamental TE[10] mode of the SIW waveguide. Prototype fabrication and characterization A commercial multilayer printed circuit board technology is used to fabricate the UMA. The SIW and bias circuitry are made of two Rogers 5880 substrates that are 1.575 and 0.787 mm-thick, respectively. These substrates have a relative dielectric constant of 2.2 and a loss tangent of 0.0009. They are bonded by a thin Rogers 4450 F film. After this, the 164 PIN diodes are mounted on the gaps of the slot openings through reflow soldering. The radiation pattern of the UMA was measured in a microwave anechoic chamber using a reconfigurable robotic measurement system (Fig. 6a, b). A signal generator (Agilent E8267D) launches monochromatic waves at 23.5 GHz to feed the UMA. A horn antenna is connected to a vector network analyzer (VNA, Keysight N9041B) to detect the radiated EM waves from the metasurface antenna. The robotic arm holds the linearly polarized horn antenna, which can rotate along a user-defined circular path from -90^o to 90^o to measure the radiation pattern of the UMA. Different polarization components of the EM waves can be measured by physically rotating the linearly polarized horn antenna. We built up an indoor wireless communication experiment to illustrate the information manipulations of our UMA (Fig. 6c). For the transmitter side, a microwave signal generator (Agilent E8267D) is adopted to generate a monochromatic wave with a frequency f[0] = 23.5 GHz to feed the UMA. The modulated waves with temporal-variant amplitude and phase properties (information) are directly generated and launched into free space by the UMA. Random binary bit streams with different modulation formats (QPSK, 8PSK, and 16QAM) are generated and mapped to FPGA's corresponding dynamic STC matrixes. For the receiver side, a linearly polarized diagonal horn antenna connected to a vector signal analyzer (VSA, Keysight N9041B) is adopted to receive and demodulate the PW radiated from the UMA. The received digital messages can be recovered by the VSA, which provides the performance of the real-time constellation diagram, eye diagram, SNR, and EVM. The distance between the UMA and the received horn antenna is ~1.2 m. The modulation frequency f[M] is set as 0.5 MHz, and the switching speed of the PIN diodes is 10 MHz. For the multiplexing wireless communication experiment, the UMA generates two independent channels in different directions (Fig. 6c). At the receiver side, two horn antennas with |u> and |u> polarizations are located at the two main beam directions, respectively. The horns are connected to two VSAs to decode the digital message of the two channels. The radiation efficiency, defined as the radiated power into free space to the incident power of the antenna, can be calculated by \(\ eta=\frac{G}{D}\)^64. The measured realized gain G of the 1-D universal metasurface antenna is 9.3 dBi at 23.5 GHz based on the gain-comparison method, which uses a standard gain horn to determine the gain of the antenna under test. Due to the unavailable cylindrical near-field scanning system for 1-D antenna directivity measurement, we utilize the simulated directivity D = 13.1 dB here. The calculated radiation efficiency of the metasurface antenna is 41.7%, corresponding to the insertion loss of -3.8 dB. According to the Friis transmission equation \({P}_{r}={P}_{t}{(\frac{\lambda }{4\ pi R})}^{2}{G}_{t}{G}_{r}\)^64, the received power for the detector depends on various factors, including the input power, the free-space loss related to the distance between the transmitter and receiver, and the gains of the transmitting and receiving antennas. In the communication link setup in Fig. 4, the measured CW received power at the detector is -33.8 dBm with a distance of 1.2 m and input power of 0 dBm for the UMA. The measured received power is -37.7 dBm in each channel for the dual-user communication link in Fig. 5. Data availability The data that support the findings of this study are presented in the paper and the Supplementary Information file. Code availability The codes that support the theoretical modeling of the UMA are available from the corresponding authors upon request. References 1. Saleh, B. E. A. & Teich, M. C. Fundamentals of Photonics (John Wiley & Sons, 2019). 2. Tse, D. & Viswanath, P. Fundamentals of Wireless Communication (Cambridge University Press, 2005). 3. Li, P. et al. Infrared hyperbolic metasurface based on nanostructured van der Waals materials. Science 359, 892-896 (2018). ADS CAS PubMed Google Scholar 4. Hu, G. et al. Topological polaritons and photonic magic angles in twisted a-MoO3 bilayers. Nature 582, 209-213 (2020). ADS CAS PubMed Google Scholar 5. Gerry, C., Knight, P. & Knight, P. L. Introductory Quantum Optics . (Cambridge university press, 2005). 6. Harrington, R. F. Time-harmonic Electromagnetic Fields (McGraw-Hill, 1961). 7. He, C., Shen, Y. & Forbes, A. Towards higher-dimensional structured light. Light Sci. Appl. 11, 1-17 (2022). Google Scholar 8. Lalanne, P., Astilean, S., Chavel, P., Cambril, E. & Launois, H. Blazed binary subwavelength gratings with efficiencies larger than those of conventional echelette gratings. Optics Lett. 23, 1081-1083 (1998). ADS CAS Google Scholar 9. Bomzon, Z. E., Biener, G., Kleiner, V. & Hasman, E. Radially and azimuthally polarized beams generated by space-variant dielectric subwavelength gratings. Optics Lett. 27, 285-287 (2002). ADS Google Scholar 10. Yu, N. et al. Light propagation with phase discontinuities: generalized laws of reflection and refraction. Science 334, 333-337 (2011). ADS CAS PubMed Google Scholar 11. Ni, X., Emani, N. K., Kildishev, A. V., Boltasseva, A. & Shalaev, V. M. Broadband light bending with plasmonic nanoantennas. Science 335, 427 (2012). ADS CAS PubMed Google Scholar 12. Arbabi, A., Horie, Y., Bagheri, M. & Faraon, A. Dielectric metasurfaces for complete control of phase and polarization with subwavelength spatial resolution and high transmission. Nat. Nanotechnol. 10, 937-943 (2015). ADS CAS PubMed Google Scholar 13. Sun, S. et al. Gradient-index meta-surfaces as a bridge linking propagating waves and surface waves. Nat. Mater. 11, 426-431 (2012). ADS CAS PubMed Google Scholar 14. Pfeiffer, C. & Grbic, A. Metamaterial Huygens' surfaces: tailoring wave fronts with reflectionless sheets. Phys. Rev. Lett. 110, 197401 (2013). ADS PubMed Google Scholar 15. Della Giovampaola, C. & Engheta, N. Digital metamaterials. Nat. Mater. 13, 1115-1121 (2014). ADS CAS PubMed Google Scholar 16. Cui, T. J., Qi, M. Q., Wan, X., Zhao, J. & Cheng, Q. Coding metamaterials, digital metamaterials and programmable metamaterials. Light Sci. Appl. 3, e218 (2014). ADS Google Scholar 17. Van de Groep, J. et al. Exciton resonance tuning of an atomically thin lens. Nat. Photonics 14, 426-430 (2020). ADS Google Scholar 18. Zhang, J. C. et al. A 6G meta-device for 3D varifocal. Sci. Adv. 9, eadf8478 (2023). ADS PubMed PubMed Central Google Scholar 19. Sherrott, M. C. et al. Experimental demonstration of> 230 phase modulation in gate-tunable graphene-gold reconfigurable mid-infrared metasurfaces. Nano Lett. 17, 3027-3034 (2017). ADS CAS PubMed Google Scholar 20. Park, J., Kang, J.-H., Kim, S. J., Liu, X. & Brongersma, M. L. Dynamic reflection phase and polarization control in metasurfaces. Nano Lett. 17, 407-413 (2017). ADS CAS PubMed Google Scholar 21. Venkatesh, S., Lu, X., Saeidi, H. & Sengupta, K. A high-speed programmable and scalable terahertz holographic metasurface based on tiled CMOS chips. Nat. Electron. 3, 785-793 (2020). CAS Google Scholar 22. Holsteen, A. L., Raza, S., Fan, P., Kik, P. G. & Brongersma, M. L. Purcell effect for active tuning of light scattering from semiconductor optical antennas. Science 358, 1407-1410 (2017). ADS MathSciNet CAS PubMed MATH Google Scholar 23. Malek, S. C., Ee, H.-S. & Agarwal, R. Strain multiplexed metasurface holograms on a stretchable substrate. Nano Lett. 17, 3641-3645 (2017). ADS CAS PubMed Google Scholar 24. Arbabi, E. et al. MEMS-tunable dielectric metasurface lens. Nat. Commun. 9, 1-9 (2018). CAS Google Scholar 25. Shcherbakov, M. R. et al. Ultrafast all-optical tuning of direct-gap semiconductor metasurfaces. Nat. Commun. 8, 1-6 (2017). CAS Google Scholar 26. Guo, X., Ding, Y., Duan, Y. & Ni, X. Nonreciprocal metasurface with space-time phase modulation. Light Sci. Appl. 8, 1-9 (2019). CAS Google Scholar 27. Alam, M. Z., De Leon, I. & Boyd, R. W. Large optical nonlinearity of indium tin oxide in its epsilon-near-zero region. Science 352, 795-797 (2016). ADS CAS PubMed Google Scholar 28. Clerici, M. et al. Controlling hybrid nonlinearities in transparent conducting oxides via two-colour excitation. Nat. Commun. 8, 1-7 (2017). Google Scholar 29. Sautter, J. et al. Active tuning of all-dielectric metasurfaces. ACS Nano 9, 4308-4315 (2015). CAS PubMed Google Scholar 30. Lewi, T., Evans, H. A., Butakov, N. A. & Schuller, J. A. Ultrawide thermo-optic tuning of PbTe meta-atoms. Nano Lett. 17, 3940-3945 (2017). ADS CAS PubMed Google Scholar 31. He, J. et al. Terahertz tunable metasurface lens based on vanadium dioxide phase transition. Plasmonics 11, 1285-1290 (2016). CAS Google Scholar 32. Thureja, P. et al. Toward a universal metasurface for optical imaging, communication, and computation. Nanophotonics 11, 3745-3768 (2022). CAS Google Scholar 33. Kafaie Shirmanesh, G., Sokhoyan, R., Pala, R. A. & Atwater, H. A. Dual-gated active metasurface at 1550 nm with wide (>300) phase tunability. Nano Lett. 18, 2957-2963 (2018). ADS CAS PubMed Google Scholar 34. Han, S. et al. Complete complex amplitude modulation with electronically tunable graphene plasmonic metamolecules. ACS Nano 14, 1166-1175 (2020). CAS PubMed Google Scholar 35. Park, J. et al. All-solid-state spatial light modulator with independent phase and amplitude control for three-dimensional LiDAR applications. Nat. Nanotechnol. 16, 69-76 (2021). ADS MathSciNet CAS PubMed Google Scholar 36. Hadad, Y., Soric, J. C. & Alu, A. Breaking temporal symmetries for emission and absorption. Proc. Natl Acad.Sci. USA 113, 3471-3475 (2016). ADS CAS PubMed PubMed Central Google Scholar 37. Cardin, A. E. et al. Surface-wave-assisted nonreciprocity in spatio-temporally modulated metasurfaces. Nat. Commun. 11, 1-9 (2020). ADS Google Scholar 38. Ramaccia, D., Sounas, D. L., Alu, A., Toscano, A. & Bilotti, F. Doppler cloak restores invisibility to objects in relativistic motion. Phys. Rev. B 95, 075113 (2017). ADS Google Scholar 39. Dai, J. Y. et al. Simultaneous in situ direction finding and field manipulation based on space-time-coding digital metasurface. IEEE Trans. Antennas Propag. 70, 4774-4783 (2022). ADS Google Scholar 40. Dai, J. Y. et al. Realization of multi-modulation schemes for wireless communication by time-domain digital coding metasurface. IEEE Trans. Antennas Propag. 68, 1618-1627 (2019). ADS Google Scholar 41. Dai, J. Y. et al. High-efficiency synthesizer for spatial waves based on space-time-coding digital metasurface. Laser Photonics Rev. 14, 1900133 (2020). ADS Google Scholar 42. Zhang, L. et al. Space-time-coding digital metasurfaces. Nat. Commun. 9, 1-11 (2018). ADS Google Scholar 43. Zhang, L. et al. A wireless communication scheme based on space-and frequency-division multiplexing using digital metasurfaces. Nat. Electron. 4, 218-227 (2021). Google Scholar 44. Wu, G.-B., Dai, J. Y., Cheng, Q., Cui, T. J. & Chan, C. H. Sideband-free space-time-coding metasurface antennas. Nat. Electron. 5, 808-819 (2022). 45. Dai, J. Y., Zhao, J., Cheng, Q. & Cui, T. J. Independent control of harmonic amplitudes and phases via a time-domain digital coding metasurface. Light Sci. Appl. 7, 1-10 (2018). Google Scholar 46. Dai, J. Y. et al. Wireless communications through a simplified architecture based on time-domain digital coding metasurface. Adv. Mater. Tech. 4, 1900044 (2019). Google Scholar 47. Wang, S. R., Chen, M. Z., Ke, J. C., Cheng, Q. & Cui, T. J. Asynchronous space-time-coding digital metasurface. Adv. Sci. 9, 2200106 (2022). CAS Google Scholar 48. Ke, J. C. et al. Frequency-modulated continuous waves controlled by space-time-coding metasurface with nonlinearly periodic phases. Light Sci. Appl. 11, 1-11 (2022). Google Scholar 49. Ke, J. C. et al. Linear and nonlinear polarization syntheses and their programmable controls based on anisotropic time-domain digital coding metasurface. Small Struct. 2, 2000060 (2021). Google Scholar 50. Ke, J. C. et al. Space-frequency-polarization-division multiplexed wireless communication system using anisotropic space-time-coding digital metasurface. Natl. Sci. Rev. https:// doi.org/10.1093/nsr/nwac225 (2022). 51. Wu, G. B., Dai, J. Y. & Chan, C. H. Four-dimensional (4-D) leaky-wave antenna for frequency translation. Asia-Pacific Microwave Conference (APMC), 267-269 (2022). 52. Wu, G. B., Dai, J. Y. & Chan, C. H. Time-modulated leaky-wave Antenna. IEEE International Symposium on Antennas and Propagation and USNC-URSI Radio Science Meeting (AP-S/URSI), 812-813 (2022). 53. Kaganovsky, Y. & Heyman, E. Wave analysis of Airy beams. Opt. Express 18, 8440-8452 (2010). ADS CAS PubMed Google Scholar 54. Buffi, A., Nepa, P. & Manara, G. Design criteria for near-field-focused planar arrays. IEEE Antennas Propag. Mag. 54, 40-50 (2012). ADS Google Scholar 55. Mutai, K. K., Sato, H. & Chen, Q. Active millimeter wave imaging using leaky-wave focusing antenna. IEEE Trans. Antennas Propag. 70, 3789-3798 (2021). ADS Google Scholar 56. Zheng, Y. et al. Metasurface-assisted wireless communication with physical level information encryption. Adv. Sci. 9, 2204558 (2022). MathSciNet Google Scholar 57. Kim, M. et al. Monolayer molybdenum disulfide switches for 6G communication systems. Nat. Electron. 5, 367-373 (2022). CAS Google Scholar 58. Samizadeh Nikoo, M. & Matioli, E. Electronic metadevices for terahertz applications. Nature 614, 451-455 (2023). ADS CAS PubMed Google Scholar 59. Zhang, Y. et al. Electrically reconfigurable non-volatile metasurface using low-loss optical phase-change material. Nat. Nanotechnol. 16, 661-666 (2021). ADS CAS PubMed Google Scholar 60. Wang, Y. et al. Electrical tuning of phase-change antennas and metasurfaces. Nat. Nanotechnol. 16, 667-672 (2021). ADS CAS PubMed Google Scholar 61. Wu, G.-B., Zhang, Q.-L., Chan, K. F., Chen, B.-J. & Chan, C. H. Amplitude-modulated leaky-wave antennas. IEEE Trans. Antennas Propag. 69, 3664-3676 (2020). ADS Google Scholar 62. Goldsmith, A. Wireless Communications (Cambridge University Press, 2005). 63. Ding, Y. & Fusco, V. F. Establishing metrics for assessing the performance of directional modulation systems. IEEE Trans. Antennas Propag. 62, 2745-2755 (2014). ADS Google Scholar 64. Balanis, C. A. Antenna Theory: Analysis and Design (John Wiley & Sons, 2015). Download references Acknowledgements This work was supported by the Hong Kong Research Grants Council of the Hong Kong SAR (T42-103/16-N, C.H.C.), the Guangdong Provincial Department of Science and Technology, China (2020B1212030002, C.H.C.), the Program of Song Shan Laboratory (Included in the management of Major Science and Technology Program of Henan Province) (221100211300-02, Q.C., 221100211300-03, J.Y.D.), the National Key Research and Development Program of China (2018YFA0701904, Q.C.), the National Natural Science Foundation of China (62288101, T.J.C. and Q.C., 62201139, J.Y.D.), the 111 Project (111-2-05, T.J.C.), the Jiangsu Province Frontier Leading Technology Basic Research Project (BK20212002, T.J.C.), the Fundamental Research Funds for the Central Universities (2242022k60003, Q.C.), the National Science Foundation (NSFC) for Distinguished Young Scholars of China (62225108, Q.C.), the Southeast University - China Mobile Research Institute Joint Innovation Center (R207010101125D9, Q.C.), and the IEEE Antennas and Propagation Society (AP-S) Fellowship (G.B.W.). Author information Author notes 1. These authors contributed equally: Geng-Bo Wu, Jun Yan Dai. Authors and Affiliations 1. State Key Laboratory of Terahertz and Millimeter Waves, City University of Hong Kong, Hong Kong, 999077, China Geng-Bo Wu, Kam Man Shum, Ka Fai Chan & Chi Hou Chan 2. Department of Electrical Engineering, City University of Hong Kong, Hong Kong, 999077, China Geng-Bo Wu & Chi Hou Chan 3. State Key Laboratory of Millimeter Waves, Southeast University, Nanjing, 210096, China Jun Yan Dai, Qiang Cheng & Tie Jun Cui 4. Institute of Electromagnetic Space, Southeast University, Nanjing, 210096, China Jun Yan Dai, Qiang Cheng & Tie Jun Cui 5. Frontiers Science Center for Mobile Information Communication and Security, Southeast University, Nanjing, 210096, China Jun Yan Dai, Qiang Cheng & Tie Jun Cui 6. Guangdong-Hong Kong Joint Laboratory for Big Data Imaging and Communication, Shenzhen, 518048, China Chi Hou Chan Authors 1. Geng-Bo Wu View author publications You can also search for this author in PubMed Google Scholar 2. Jun Yan Dai View author publications You can also search for this author in PubMed Google Scholar 3. Kam Man Shum View author publications You can also search for this author in PubMed Google Scholar 4. Ka Fai Chan View author publications You can also search for this author in PubMed Google Scholar 5. Qiang Cheng View author publications You can also search for this author in PubMed Google Scholar 6. Tie Jun Cui View author publications You can also search for this author in PubMed Google Scholar 7. Chi Hou Chan View author publications You can also search for this author in PubMed Google Scholar Contributions Q.C., T.J.C., and C.H.C. suggested the designs, planned and supervised the entire study, and led the project. G.B.W. and J.Y.D. conceived the idea of this work and designed the metasurface. G.B.W., J.Y.D., K.M.S., K.F.C. carried out the measurements and data analysis. All authors contributed to the writing of the paper. All authors discussed the theoretical modeling and numerical simulations and reviewed the manuscript. Corresponding authors Correspondence to Jun Yan Dai, Qiang Cheng, Tie Jun Cui or Chi Hou Chan. Ethics declarations Competing interests The authors declare no competing interests. Peer review Peer review information Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available. Additional information Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary information Supplementary Information Peer Review File Rights and permissions Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Reprints and Permissions About this article Check for updates. Verify currency and authenticity via CrossMark Cite this article Wu, GB., Dai, J.Y., Shum, K.M. et al. A universal metasurface antenna to manipulate all fundamental characteristics of electromagnetic waves. Nat Commun 14, 5155 (2023). https://doi.org/10.1038/ s41467-023-40717-9 Download citation * Received: 27 January 2023 * Accepted: 01 August 2023 * Published: 24 August 2023 * DOI: https://doi.org/10.1038/s41467-023-40717-9 Share this article Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative Comments By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Download PDF Advertisement Advertisement Explore content * Research articles * Reviews & Analysis * News & Comment * Videos * Collections * Subjects * Follow us on Facebook * Follow us on Twitter * Sign up for alerts * RSS feed About the journal * Aims & Scope * Editors * Journal Information * Open Access Fees and Funding * Calls for Papers * Editorial Values Statement * Journal Metrics * Editors' Highlights * Contact * Editorial policies * Top Articles Publish with us * For authors * For Reviewers * Language editing services * Submit manuscript Search Search articles by subject, keyword or author [ ] Show results from [All journals] Search Advanced search Quick links * Explore articles by subject * Find a job * Guide to authors * Editorial policies Nature Communications (Nat Commun) ISSN 2041-1723 (online) nature.com sitemap About Nature Portfolio * About us * Press releases * Press office * Contact us Discover content * Journals A-Z * Articles by subject * Protocol Exchange * Nature Index Publishing policies * Nature portfolio policies * Open access Author & Researcher services * Reprints & permissions * Research data * Language editing * Scientific editing * Nature Masterclasses * Live Expert Trainer-led workshops * Research Solutions Libraries & institutions * Librarian service & tools * Librarian portal * Open research * Recommend to library Advertising & partnerships * Advertising * Partnerships & Services * Media kits * Branded content Professional development * Nature Careers * Nature Conferences Regional websites * Nature Africa * Nature China * Nature India * Nature Italy * Nature Japan * Nature Korea * Nature Middle East * Privacy Policy * Use of cookies * Your privacy choices/Manage cookies * Legal notice * Accessibility statement * Terms & Conditions * Your US state privacy rights Springer Nature (c) 2023 Springer Nature Limited Close Nature Briefing Sign up for the Nature Briefing newsletter -- what matters in science, free to your inbox daily. Email address [ ] Sign up [ ] I agree my information will be processed in accordance with the Nature and Springer Nature Limited Privacy Policy. Close Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing * *