Independent control of harmonic amplitudes and phases via a time-domain digital coding metasurface

  • Light: Science & Applications  7, Article number: 90 (2018)
More Information
  • Corresponding author:
    Qiang Cheng ( Jun Cui (
  • Received: 16 May 2018
    Revised: 23 October 2018
    Accepted: 03 November 2018
    Published online: 21 November 2018


  • Harmonic manipulations are important for applications such as wireless communications, radar detection and biological monitoring. A general approach to tailor the harmonics involves the use of additional amplifiers and phase shifters for the precise control of harmonic amplitudes and phases after the mixing process; however, this approach leads to issues of high cost and system integration. Metasurfaces composed of a periodic array of subwavelength resonators provide additional degrees of freedom to realize customized responses to incident light and highlight the possibility for nonlinear control by taking advantage of time-domain properties. Here, we designed and experimentally characterized a reflective time-domain digital coding metasurface, with independent control of the harmonic amplitude and phase. As the reflection coefficient is dynamically modulated in a predefined way, a large conversion rate is observed from the carrier signal to the harmonic components, with magnitudes and phases that can be accurately and separately engineered. In addition, by encoding the reflection phases of the meta-atoms, beam scanning for multiple harmonics can be implemented via different digital coding sequences, thus removing the need for intricate phase-shift networks. This work paves the way for efficient harmonic control for applications in communications, radar, and related areas.
  • 加载中
  • [1] Ni, X. J., Kildishev, A. V. & Shalaev, V. M. Metasurface holograms for visible light. Nat. Commun. 4, 2807 (2013). doi:  10.1038/ncomms3807
    [2] Huang, L. L. et al. Three-dimensional optical holography using a plasmonic metasurface. Nat. Commun. 4, 2808 (2013). doi:  10.1038/ncomms3808
    [3] Liu, L. X. et al. Broadband metasurfaces with simultaneous control of phase and amplitude. Adv. Mater. 26, 5031-5036 (2014). doi:  10.1002/adma.201401484
    [4] Zheng, G. X. et al. Metasurface holograms reaching 80% efficiency. Nat. Nanotechnol. 10, 308-312 (2015). doi:  10.1038/nnano.2015.2
    [5] Landy, N. I., Sajuyigbe, S., Mock, J. J., Smith, D. R. & Padilla, W. J. Perfect metamaterial absorber. Phys. Rev. Lett. 100, 207402 (2008). doi:  10.1103/PhysRevLett.100.207402
    [6] Liu, N., Mesch, M., Weiss, T., Hentschel, M. & Giessen, H. Infrared perfect absorber and its application as plasmonic sensor. Nano Lett. 10, 2342-2348 (2010). doi:  10.1021/nl9041033
    [7] Mei, J. et al. Dark acoustic metamaterials as super absorbers for low-frequency sound. Nat. Commun. 3, 756 (2012). doi:  10.1038/ncomms1758
    [8] Watts, C. M., Liu, X. L. & Padilla, W. J. Metamaterial electromagnetic wave absorbers. Adv. Mater. 24, OP98-OP120 (2012).
    [9] Akselrod, G. M. et al. Large‐area metasurface perfect absorbers from visible to near‐infrared. Adv. Mater. 27, 8028-8034 (2015). doi:  10.1002/adma.201503281
    [10] Sounas, D. L., Caloz, C. & Alù, A. Giant non-reciprocity at the subwavelength scale using angular momentum-biased metamaterials. Nat. Commun. 4, 2407 (2013). doi:  10.1038/ncomms3407
    [11] Popa, B. I. & Cummer, S. A. Non-reciprocal and highly nonlinear active acoustic metamaterials. Nat. Commun. 5, 3398 (2014). doi:  10.1038/ncomms4398
    [12] Sievenpiper, D. F., Schaffner, J. H., Song, H. J., Loo, R. Y. & Tangonan, H. Two-dimensional beam steering using an electrically tunable impedance surface. IEEE Trans. Antennas Propag. 51, 2713-2722 (2003). doi:  10.1109/TAP.2003.817558
    [13] Yao, Y. et al. Electrically tunable metasurface perfect absorbers for ultrathin mid-infrared optical modulators. Nano Lett. 14, 6526-6532 (2014). doi:  10.1021/nl503104n
    [14] Dabidian, N. et al. Electrical switching of infrared light using graphene integration with plasmonic Fano resonant metasurfaces. ACS Photonics 2, 216-227 (2015). doi:  10.1021/ph5003279
    [15] Chen, K. et al. A reconfigurable active Huygens' metalens. Adv. Mater. 29, 1606422 (2017). doi:  10.1002/adma.201606422
    [16] Chen, J. N. et al. Optical nano-imaging of gate-tunable graphene plasmons. Nature 487, 77-81 (2012). doi:  10.1038/nature11254
    [17] Rodrigo, D. et al. Mid-infrared plasmonic biosensing with graphene. Science 349, 165-168 (2015). doi:  10.1126/science.aab2051
    [18] Lee, J. et al. Ultrafast electrically tunable polaritonic metasurfaces. Adv. Opt. Mater. 2, 1057-1063 (2014). doi:  10.1002/adom.201400185
    [19] 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). doi:  10.1038/lsa.2014.99
    [20] Della Giovampaola, C. & Engheta, N. Digital metamaterials. Nat. Mater. 13, 1115-1121 (2014). doi:  10.1038/nmat4082
    [21] Cui, T. J., Liu, S., & Li, L. L. Information entropy of coding metasurface. Light Sci. Appl. 5, e16172 (2016).
    [22] Cui, T. J. Microwave metamaterials. Natl. Sci. Rev. 5, 134-136 (2018). doi:  10.1093/nsr/nwx133
    [23] Cui, T. J., Liu, S. & Zhang, L. Information metamaterials and metasurfaces. J. Mater. Chem. C 5, 3644-3668 (2017). doi:  10.1039/C7TC00548B
    [24] Cui, T. J. Microwave metamaterials-from passive to digital and programmable controls of electromagnetic waves. J. Opt. 19, 084004 (2017). doi:  10.1088/2040-8986/aa7009
    [25] Debogovic, T. & Perruisseau-Carrier, J. Low loss MEMS-reconfigurable 1-bit reflectarray cell with dual-linear polarization. IEEE Trans. Antennas Propag. 62, 5055-5060 (2014). doi:  10.1109/TAP.2014.2344100
    [26] Wan, X., Qi, M. Q., Chen, T. Y. & Cui, T. J. Field-programmable beam reconfiguring based on digitally-controlled coding metasurface. Sci. Rep. 6, 20663 (2016). doi:  10.1038/srep20663
    [27] Li, Y. B. et al. Transmission-type 2-bit programmable metasurface for single-sensor and single-frequency microwave imaging. Sci. Rep. 6, 23731 (2016). doi:  10.1038/srep23731
    [28] Huang, C. et al. Dynamical beam manipulation based on 2-bit digitally-controlled coding metasurface. Sci. Rep. 7, 42302 (2017). doi:  10.1038/srep42302
    [29] Li, L. L. et al. Electromagnetic reprogrammable coding-metasurface holograms. Nat. Commun. 8, 197 (2017). doi:  10.1038/s41467-017-00164-9
    [30] Hadad, Y., Sounas, D. L. & Alu, A. Space-time gradient metasurfaces. Phys. Rev. B 92, 100304 (2015). doi:  10.1103/PhysRevB.92.100304
    [31] Shaltout, A., Kildishev, A. & Shalaev, V. Time-varying metasurfaces and Lorentz non-reciprocity. Opt. Mater. Exp. 5, 2459-2467 (2015). doi:  10.1364/OME.5.002459
    [32] Liu, Z. Z., Li, Z. Y. & Aydin, K. Time-varying metasurfaces based on graphene microribbon arrays. ACS Photonics 3, 2035-2039 (2016). doi:  10.1021/acsphotonics.6b00653
    [33] Stewart, S. A., Smy, T. J. & Gupta, S. Finite-difference time-domain modeling of space-time-modulated metasurfaces. IEEE Trans. Antennas Propag. 66, 281-292 (2017). doi:  10.1109/TAP.2017.2772045
    [34] Li, G. X. et al. Continuous control of the nonlinearity phase for harmonic generations. Nat. Mater. 14, 607-612 (2015). doi:  10.1038/nmat4267
    [35] Rose, A., Huang, D. & Smith, D. R. Controlling the second harmonic in a phase-matched negative-index metamaterial. Phys. Rev. Lett. 107, 063902 (2011). doi:  10.1103/PhysRevLett.107.063902
    [36] Lee, J. et al. Giant nonlinear response from plasmonic metasurfaces coupled to intersubband transitions. Nature 511, 65-69 (2014). doi:  10.1038/nature13455
    [37] Klein, M. W., Enkrich, C., Wegener, M. & Linden, S. Second-harmonic generation from magnetic metamaterials. Science 313, 502-504 (2006). doi:  10.1126/science.1129198
    [38] Lapine, M., Shadrivov, I. V., Powell, D. A. & Kivshar, Y. S. Magnetoelastic metamaterials. Nat. Mater. 11, 30-33 (2011). doi:  10.1038/nmat3168
    [39] Zhao, J. et al. Controlling spectral energies of all harmonics in programmable way using time-domain digital coding metasurface. Preprint at Accessed June 2018.
    [40] Rocca, P., Zhu, Q. J., Bekele, E. T., Yang, S. W. & Massa, A. 4-D arrays as enabling technology for cognitive radio systems. IEEE Trans. Antennas Propag. 62, 1102-1116 (2014). doi:  10.1109/TAP.2013.2288109
    [41] Secmen M., Demir S., Hizal A., & Eker T. Frequency diverse array antenna with periodic time modulated pattern in range and angle. In Proc. of 2007 IEEE Radar Conference 427-430 (IEEE, Boston, MA, USA, 2007).
    [42] Poli, L., Rocca, P., Manica, L. & Massa, A. Handling sideband radiations in time-modulated arrays through particle swarm optimization. IEEE Trans. Antennas Propag. 58, 1408-1411 (2010). doi:  10.1109/TAP.2010.2041165
    [43] Poli, L., Rocca, P., Oliveri, G. & Massa, A. Harmonic beamforming in time-modulated linear arrays. IEEE Trans. Antennas Propag. 59, 2538-2545 (2011). doi:  10.1109/TAP.2011.2152323
    [44] Poli, L., Rocca, P., Manica, L. & Massa, A. Pattern synthesis in time-modulated linear arrays through pulse shifting. IET Microw. Antennas Propag. 4, 1157-1164 (2010). doi:  10.1049/iet-map.2009.0042
    [45] Zhao, J. et al. A tunable metamaterial absorber using varactor diodes. New J. Phys. 15, 043049 (2013). doi:  10.1088/1367-2630/15/4/043049
    [46] Liu, S. et al. Convolution operations on coding metasurface to reach flexible and continuous controls of terahertz beams. Adv. Sci. 3, 1600156 (2016). doi:  10.1002/advs.201600156
    [47] Wang, D., Liu, Z. G., Zhao, J., Cheng, Q. & Cui, T. J. Accurate design of low backscattering metasurface using iterative fourier transform algorithm. Sci. Rep. 7, 11346 (2017). doi:  10.1038/s41598-017-11719-7
通讯作者: 陈斌,
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(7) / Tables(3)

Article Metrics

Article views(3) PDF downloads(0) Citation(0) Citation counts are provided from Web of Science. The counts may vary by service, and are reliant on the availability of their data.

Independent control of harmonic amplitudes and phases via a time-domain digital coding metasurface

  • 1. State Key Laboratory of Millimeter Waves, Southeast University, 210096, Nanjing, China
  • 2. Synergetic Innovation Center of Wireless Communication Technology, Southeast University, 210096, Nanjing, China
  • Corresponding author: Qiang Cheng,
    Corresponding author: Tie Jun Cui,


Abstract: Harmonic manipulations are important for applications such as wireless communications, radar detection and biological monitoring. A general approach to tailor the harmonics involves the use of additional amplifiers and phase shifters for the precise control of harmonic amplitudes and phases after the mixing process; however, this approach leads to issues of high cost and system integration. Metasurfaces composed of a periodic array of subwavelength resonators provide additional degrees of freedom to realize customized responses to incident light and highlight the possibility for nonlinear control by taking advantage of time-domain properties. Here, we designed and experimentally characterized a reflective time-domain digital coding metasurface, with independent control of the harmonic amplitude and phase. As the reflection coefficient is dynamically modulated in a predefined way, a large conversion rate is observed from the carrier signal to the harmonic components, with magnitudes and phases that can be accurately and separately engineered. In addition, by encoding the reflection phases of the meta-atoms, beam scanning for multiple harmonics can be implemented via different digital coding sequences, thus removing the need for intricate phase-shift networks. This work paves the way for efficient harmonic control for applications in communications, radar, and related areas.

  • Metasurfaces provide an unprecedented route to control light propagation and engineer light-matter interaction using an array of resonators with carefully designed geometries. With the introduction of an abrupt phase discontinuity along the interface, the electromagnetic (EM) properties of incident light can be tailored in a controlled manner as governed by the generalized Snell's law, resulting in full control of the propagation direction, polarization, and wavefront shape over subwavelength distances. These distinctive features have driven a rich variety of new technologies such as holographic imaging, perfect absorption and nonreciprocal transmission1-11. With the rapid advent of various tuning technologies, a number of tunable metasurfaces have been shown to provide in situ dynamic control of electromagnetic waves, which opens a new route for a variety of applications such as beam shaping, active focusing and biosensing12-18.

    However, one challenging task for the construction of a tunable metasurface involves the design of an intricate feeding network that provides various excitation intensities and phases to each element with high precision. In this regard, an alternative method has been explored, namely, coding the metasurface and digital metasurface, which shapes the phase front of the electromagnetic waves using binary meta-atoms with carefully designed geometries19-24. Such a configuration is extremely helpful to simplify the feeding network since only the discrete ON or OFF states for the modulation are required to manipulate the radiation or scattering properties of electromagnetic waves25-29.

    In addition to engineering the spectral response with a space-gradient metasurface, the time discontinuities of the element phase provide an additional degree of freedom to control the normal momentum component at the interface of the time-varying metasurface, leading to a break in Lorentz reciprocity during the light-matter interactions without the use of nonlinear materials30-39. Due to the rapid growth of communication technologies, the time-modulated antenna was developed to radiate harmonics with tight control of wave behavior so that information could be delivered through multichannels40-44. However, such an antenna lacks the ability to modulate incoming waves with efficient harmonic manipulation.

    Here, we report on a scheme to achieve individual control of the amplitude and phase for harmonics via a reflective time-domain digital coding metasurface. As the reflection phase of the meta-atom is periodically switched between two states, a series of harmonics emerges, with their intensities determined by the phase difference within the switching function. Meanwhile, extra time delays are imposed on the modulation signal of the meta-atoms to create a metasurface with nonuniform harmonic phase distributions, thereby enabling wavefront reshaping for high-order harmonics. When discrete phase states of the meta-atoms are encoded into the binary sequences, control of multiple harmonics including the beam directions and intensities is realized using elaborate coding strategies. All measured signals agree closely with the analytical predictions, with these unique features likely to have key roles in future radar and communication systems.

  • We designed a reflective time-domain digital coding metasurface (Fig. 1) with each unit loaded with varactor diodes. The phase response of the meta-atom can be accurately tailored within a wide phase range (~270°) by tuning the biasing voltage of the diode. This can be easily implemented when the DC feed of the element is connected to the peripheral circuit responsible for the generation of the desired control signal level, as will be discussed later. A three-dimensional schematic view of the meta-atom is shown in Fig. 2a. The side and central rectangular patches of the proposed element are bridged by four varactor diodes in total. A slotted copper layer lies at the bottom of the element to act as the DC feeding network, where a biasing voltage is applied to the diode via holes. A more detailed element description is provided in Fig. 2b, where Px = 37 mm, Py = 33 mm, H = 4 mm, M = 16 mm, N = 8.45 mm, L = 7.45 mm, g = 0.5 mm, d = 2.1 mm, t = 4.6 mm, Φ (diameter of the hole) = 1 mm, and s = 0.3 mm. The top and bottom layers are spaced by a substrate (F4B) with a dielectric constant and loss tangent of 2.65 and 0.001, respectively. To prohibit wave transmission through the slots of the feeding layer, an extra copper layer is placed at the bottom of the metasurface, eliminating the possibility of energy leakage. In general, the varactor diode (SMV-2019, Skyworks, Inc.) can be regarded as an RLC model at the operation frequency with the equivalent circuit parameters listed in ref. 45. To grasp a complete picture of the reflection spectra for the meta-atom, a full wave simulation was performed using a commercial electromagnetic solver (CST Microwave Studio 2016) under different biasing voltages. The reflection minima was observed to redshift (Fig. 2c) as the biasing voltage was changed from 19 to 0 V. Good phase linearity and a large phase range near 3.7 GHz was observed, as shown in Fig. 2d. This enables beam steering based on the relative phase shift of adjacent elements similar to the principle of the phase array antenna. The layout of the metasurface consists of 7×8 elements (Fig. 2). For simplicity, all the varactor diodes in each column share the same biasing voltage. Targeting the discrete phases represented by the codes in our design, a series of biasing voltages, i.e., 0, 3, 6, 9, 12, 15, 18, and 21 V, were adopted to meet the phase demand in the experiment. The mapping relationship between the reflection phase and the biasing voltage for the meta-atom in the simulation is presented in Table 1, which provides a basis for the DC level on the diodes for each element in the following experiment.

    Fig. 1  Schematic diagram of the time-domain digital coding metasurface, which is able to generate nonlinear harmonics under excitation by a monochromatic wave as well as control their direction and magnitude

    Fig. 2  Illustration of the time-domain digital coding meta-atom.

    a, b Perspective and front views of the time-domain digital coding meta-atom. c, d The simulated reflection amplitude (c) and phase (d) spectra under different bias voltages. The dimensions of the meta-atom are Px = 37 mm, Py = 33 mm, H = 4 mm, M = 16 mm, N = 8.45 mm, L = 7.45 mm, g = 0.5 mm, d = 2.1 mm, t = 4.6 mm, Φ = 1 mm, s = 0.3 mm
    Voltage (V)036912151821
    Phase (Deg.)0103090180210250270

    Table 1.  The mapping relationship between the reflection phase of the meta-atom and biasing voltage at 3.7 GHz in the simulation

  • During the experiment, two horn antennae were used to illuminate the sample and receive the reflected signal, respectively. Both the microwave signal generator and the spectrum analyzer were utilized to monitor the nonlinear properties of the metasurface over a broad spectrum when connected to the horn antennae separately through phase stable cables. The sample was fabricated using standard PCB technology, with all components described in Fig. 3a and b. A FPGA controller, a digital-analog conversion (DAC) module and an analog amplifier module were employed to ensure a high output swing in accordance with the large range of regulation voltages used for the selected varactors. The operation frequency f0 was 3.7 GHz and the control circuit was programmed to generate square wave control signals with different voltages, periods and time delays.

    Fig. 3  Schematic and photograph of the time-domain digital coding metasurface.

    a Systematic description of the time-domain digital coding metasurface. b Photograph of the fabricated sample (on the left) in the experiment. The control system is presented in the right of b, consisting of a DC power source, FPGA platform (NI Corp.), DAC module and an external amplifier module

    Prior to measurement of the scattering patterns of the harmonics, it was important to observe the nonlinear generation capability of the sample through experiments (see theory in Materials and methods). Using Table 1, which relates the biasing voltage to the reflection phase, it is easy to identify the combination of biasing voltages V1/V2 corresponding to the required phases φ1 and φ2 in Eq. 4. During the experiment, we used a spectrum analyzer to directly measure the harmonic intensities. The harmonic phases were obtained through Fast Fourier Transformation of the echo signal from the metasurface, which was recorded by a software-defined radio reconfigurable device (USRP-2943, National Instruments Corp.). Figure 4 illustrates the relationship between the measured harmonic intensities/phases and the biasing voltage V1/V2 as well as the modulation period T when all the varactor diodes of the metasurface are biased with the same control signal.

    Fig. 4  The measured harmonic intensities/phases distribution of the time-domain digital coding metasurface at 3.7 GHz, where different biasing voltages V1/V2 and modulation periods T are considered.

    The same modulation period T = 6.4 μs (a-c), T = 3.2 μs (d-f) and T = 1.6 μs (g-i) are applied in each column, corresponding to a frequency interval of 156.25 kHz, 312.5 kHz and 625 kHz, while the same biasing voltage V1/V2 = 12 V/21 V (a-g), V1/V2 = 9 V/18 V (b-h), V1/V2 = 0 V/12 V (c-i) is adopted in each row, respectively

    Although the measured results in Fig. 4 indicate that the mapping relationship between the biasing voltage and reflection phase is not entirely consistent with the designs in Table 1, the nonlinear phenomenon still occurs. Consistent with theoretical predictions, the measured results clearly reproduce the trend of growing harmonics as φ1φ2 approaches π. In contrast, the intensity of the synchronous component experiences a significant drop as most of the reflection energy is shifted into the harmonic channels, which is expected to be totally eradicated in the case of opposite phases for φ1 and φ2 from Eq. (8). Such a discrepancy may be ascribed to phase deviations of the meta-atoms from their ideal values under different biasing voltages, largely stemming from the parasitic parameters of the varactor diode, the fabrication tolerance and distortion of the control signal. The last issue is caused by the small switching time of the biasing signal that is beyond the bandwidth of the control circuit, which in turn damages the edges of the square waveform.

    Figure 4 reveals that the modulation period has nearly no influence on the harmonic intensity but plays a vital role in determining the positions of the spectral lines. With a decrease of the modulation period from 6.4 to 1.6 μs, the observed frequency gap between adjacent harmonics is increased by up to 625 kHz in the experiment, proportional to 1/T. In addition to controlling the spectral position of the first harmonic, the possibility of amplitude modulation with the metasurface is shown in Fig. 4. For example, the amplitude of the +1st order harmonic is closely associated with the biasing voltage combination V1 and V2 from each column in Fig. 4. Therefore, the task of amplitude modulation can be easily implemented by tuning the external biasing voltages.

    The designed metasurface offers a wide dynamic range of up to 25 dB for the +1st order harmonic amplitude with respect to different voltage combinations. This is ideally suited for accurate tuning of the harmonic intensity in practice. For better demonstration in the following experiment, three sets of parameters with 0 V/12 V (A1), 9 V/18 V (A2), and 12 V/21 V (A3) were picked up manually to meet the requirement for amplitude attenuation by 0, 5, and 10 dB for the +1st order harmonic when reflected from the sample.

    We then experimentally investigated the modification of the scattering pattern of the +1st order harmonic via the proposed metasurface. The modulation period was T = 6.4 μs, indicating that the +1st order harmonic operates at a frequency of 3.70015625 GHz. Following the derivation in Eq. (8), the control signal shift in time t0 can lead to an extra phase lag ω0t0 for the +1st order harmonic without altering its magnitude (see Materials and methods). Thus, the time delay 0(0 μs) and T/2(3.2 μs) can be employed respectively to represent the coded elements '0′ and '1' with opposite phase. Consequently, the reflection phases of all the columns in the metasurface (Fig. 3a) can be described by the coding sequence of '0' and '1' bits. In case of the metasurface encoded by '00000000′ with no phase difference for all the columns, a directive scattering beam was measured for the +1st order harmonic along the normal direction, as shown in Fig. 5b. The biasing voltages also yielded magnitude control for the scattering pattern without changing its profile. By varying the voltage combination from A1 to A3, a decay rate of 0, 5, and 10 dB was observed for the scattering magnitude from the red, green and blue lines in Fig. 5b. Further characterization of the scattering patterns for the codes '00001111' and '00110011' can be found in Fig. 5d-f, in which the main scattering lobe was beamed toward different angular directions via the interference from the columns, enabling a large range of magnitude adjustment via control of the biasing voltages similar to that in Fig. 5b. By comparing Fig. 5a-e and Fig. 5b-f, one can observe a good correspondence between the theoretically calculated and experimentally measured results. Such agreement suggests potential for the beam forming the +1st order harmonic via the time-domain digital coding metasurface with sufficient accuracy.

    Fig. 5  Calculation and measurement results for the +1st order harmonic of the metasurface.

    The calculated (a, c, e) and measured (b, d, f) E-plane scattering patterns of the +1st order harmonic under different coding sequences of the metasurface. The first to third columns correspond to the coding sequences of '00000000', '00001111' and '00110011', respectively. The red, green and blue lines in (b, d, f) demonstrate variance of the scattering magnitude due to changing the voltage combination from A1 to A3

    Finally, we explored the feasibility of multiharmonic control with the same strategy. As mentioned above, the time delay t0 of the control signal introduced an additional phase shift 0t0 to the harmonic of kth order. This interesting property hints toward a novel route for controlling the scattering patterns of multiple harmonics at the same time. Here, we aimed to simultaneously tailor the scattering features of −3rd, −1st, +1st, and +3rd order harmonics in experiment. To achieve a more flexible control of the harmonics, the 3-bit meta-atoms with eight phase states were adopted by applying various time delays t0 to the square wave functions. The codes '0', '1', '2', '3', '4', '5', '6', '7' were used to represent 0, π/4, π/2, 3π/4, π, 5π/4, 3π/2, 7π/4 for the +1st order harmonic when t0 equals 0 (0 μs), T/8 (0.8 μs), T/4 (1.6 μs), 3 T/8 (2.4 μs), T/2 (3.2 μs), 5 T/8 (4 μs), 3 T/4 (4.8 μs), 7 T/8 (5.6 μs). However, the coding implications led to a great difference with changing harmonic order, which corresponds to (0, −π/4, −π/2, −3π/4, −π, −5π/4, −3π/2, −7π/4), (0, 3π/4, 3π/2, 9π/4, 3π, 15π/4, 9π/2, 21π/4) and (0, −3π/4, −3π/2, −9π/4, −3π, −15π/4, −9π/2, −21π/4) for −1st, +3rd, and −3rd order harmonics, respectively. This change means that these harmonics have different scattering orientations under the same coding sequences. The −3rd, −1st, +1st, and +3rd order harmonics were set to operate at 3.69953125, 3.69984375, 3.70015625, and 3.70046875 GHz, respectively. We investigated the waveform modulation phenomena over the aforementioned four harmonics when switching among the predesigned coding sequences. Figure 6 demonstrates the measured E-plane scattering patterns of the −3rd (purple line), −1st (blue line), +1st (red line) and +3rd (green line) order harmonics monitored by the spectrum analyzer, along with the corresponding simulation results. Five coding schemes were applied to the eight columns of the metasurface from left to right (Fig. 3), including '00000000', '00112334', '00111122', '33221100', and '76543210'. The calculated results (dashed lines) show how the variance of the coding sequence affects the scattering response of the four harmonics under normal illumination by plane waves at 3.7 GHz, based on the Fourier transformation between the surface currents and the far fields46, 47. Due to the linear relationship between the phase shift and the harmonic order, it was possible to achieve a larger deflection angle for high-order harmonics. This can be confirmed by comparing the scattering patterns with the same coding sequence in each line of Fig. 6. Moreover, symmetric scattering patterns could be found for the harmonic pairs (−1st, +1st) and (−3rd, +3rd) (Fig. 6) since opposite phase distributions were obtained with the same coding sequence, which thus gave rise to reversed in-plane momentum added to that of the incident wave. The close agreement between the measured and calculated scattering angles suggests that the time-domain digital coding metasurface can serve as a good candidate to realize a spatial scan with multiple harmonics, e.g., one can simultaneously scan the angular region from 0° to 90° by the +1st order harmonic and the region from −90° to 0° by the −1st order harmonic, and thus reduce the time cost for measurement by half.

    Fig. 6  The calculation and experimental results for multi order harmonics of the metasurface.

    Calculated (dashed line) and experimental (solid line) E-plane scattering patterns for the −3rd (purple lines), −1st (blue lines), +1st (red lines) and +3rd (green lines) order harmonics under the different coding sequences shown on the left
  • In this study, we created a time-domain digital coding metasurface to generate harmonics under the illumination of electromagnetic waves with a high conversion rate, which enables independent control of the harmonic amplitude and phase. By controlling the biasing voltages of the varactor diodes incorporated into the meta-atoms, it was possible to acquire a periodic time-variant reflection coefficient responsible for regulating the harmonic intensity. Meanwhile, the time delay of the switching functions for the meta-atoms can be carefully selected to modify the harmonic reflection phase, so that the metasurface can be employed to reshape the wavefront for multiple harmonics. Although our work was experimentally confirmed at microwave frequencies, it can be further extended into the THz and light regime in combination with advanced modulation technologies32, 33. Furthermore, the transmission type metasurface can also be employed to realize harmonic manipulation when it is fed by a patch antenna array from behind.

Materials and methods
  • We started by considering the reflection problem when EM waves are illuminated normally from the upper free space upon a reflective time-domain digital coding metasurface. The temporal expression of the reflected wave can be written as ${{E}}_r(t) = \Gamma (t) \cdot {{E}}_i(t)$, where ${{E}}_r(t),{{E}}_i(t)$ and $\Gamma (t)$ denote the reflected wave, incident wave, and reflection coefficient, respectively. Specifically, in the case of monochromic incidence at frequency fc where ${{E}}_i(t) = e^{ - j\omega _ct}$, the reflected wave in the frequency domain reads

    $$ {{E}}_r(\omega ) = \Gamma (\omega ) \ast \left[ {\delta \left( {\omega - \omega _c} \right)} \right] = \Gamma (\omega - \omega _c) $$ (1)

    where ωc the is angular frequency and δ(ωωc) represents the Dirac delta function at ω = ωc. From Eq. 1, it is straightforward to find that for a time-invariant metasurface with constant Γ(t), the reflected signal only contains information for the carrier frequency ωc. However, if the reflection coefficient Γ(t) turns out to be a periodic signal, it can be expressed as a linear combination of harmonically related complex exponentials

    $$ \Gamma (t) = \mathop {\sum}\limits_{k = - \infty }^{ + \infty } {a_ke^{jk\omega _0t}} = \mathop {\sum}\limits_{k = - \infty }^{ + \infty } {a_ke^{jk\frac{{2\pi }}{T}t}} $$ (2)

    with the corresponding spectral expression for the reflected wave given as follows

    $$ {{E}}_r(\omega ) = 2\pi \mathop {\sum}\limits_{k = - \infty }^{ + \infty } {a_k{{E}}_i(\omega - k\omega _0)} $$ (3)

    in which ω = 2π/T represents the angular frequency determined by the function period T and ak is the coefficient of the kth harmonic component. Due to strong wave-matter interactions, a nonlinear phenomenon occurs as characterized by the emergence of numerous harmonics beside the fundamental components. Here, we focus on a specific scenario of a lossless metasurface with a reflection phase response described by the periodic square wave

    $$ \Gamma (t) = Ae^{j\left\{ {\varphi _1 + (\varphi _2 - \varphi _1)\mathop {\sum}\limits_{n = - \infty }^{ + \infty } {\left[ {\varepsilon \left( {t - nT} \right) - \varepsilon \left( {t - \frac{T}{2} - nT} \right)} \right]} } \right\}} $$ (4)

    where A is the constant reflection amplitude, φ1, φ2 are the two phase states of the square wave and ε(t nT) represents the unit step function shifted by nT. From Eqs. (2)-(4), the harmonic coefficient is given by:

    $$ a_k = \left\{ {\begin{array}{*{20}{l}} {A{\kern 1pt} {{cos}}\frac{{\varphi _2 - \varphi _1}}{2}e^{j\frac{{\varphi _2 + \varphi _1}}{2}}} \hfill & {k = 0} \hfill \\ {\frac{{2A}}{{k\pi }}\sin \frac{{\varphi _2 - \varphi _1}}{2}e^{j\frac{{\varphi _2 + \varphi _1}}{2}}} \hfill & {k = \pm 1, \pm 3, \pm 5 \cdots} \hfill \\ 0 \hfill & {k = \pm 2, \pm 4, \pm 6 \cdots} \hfill \end{array}} \right. $$ (5)

    The Fourier transform of Eq. (5) yields

    $$ \begin{array}{*{20}{l}} {{{E}}_r(\omega )} \hfill & = \hfill & {2\pi A{\kern 1pt} {{cos}}\frac{{\varphi _2 - \varphi _1}}{2}e^{j\frac{{\varphi _2 + \varphi _1}}{2}}{{E}}_i(\omega )} \hfill \\ {} \hfill & {} \hfill & \hskip -20pt{ + {\kern 1pt} \mathop {\sum}\limits_{m = - \infty }^{ + \infty } {\frac{{4A}}{{2m - 1}}{\kern 1pt} {{sin}}{\kern 1pt} \frac{{\varphi _2 - \varphi _1}}{2}e^{j\frac{{\varphi _2 + \varphi _1}}{2}}{{E}}_i\left[ {\omega - (2m - 1)\omega _0} \right]}} \hfill \end{array} $$ (6)

    Only the synchronous component (k = 0) and odd harmonics are found to survive in the reflected waves, as can be understood in terms of the Fourier Transform of square waves. The envelope of the reflection spectra is highly associated with the phase pair (φ1, φ2), which in turn offers the opportunity to dynamically control the amplitude and phase distributions of high-order harmonics. To further illustrate this principle, the dependence of the reflection spectra on φ1 and φ2 is illustrated in Fig. 7. Three sets of phase combinations are taken into account with φ1/φ2 = 0°/36.87°, 0°/68.44° and 0°/180°, respectively. As the phase difference Δφ = φ1φ2 gradually approaches 180°, the intensity of the synchronous component tends to be significantly suppressed, as shown in Fig. 7a-c. This phenomenon is consistent with the mathematical predictions in Eq. 6, where the first term that represents the synchronous component is proportional to the cosine function of Δφ/2. Meanwhile, as a consequence, the harmonic energies are augmented, resulting in efficient frequency conversion. When the opposite phase requirement between φ1 and φ2 is satisfied, the synchronous components are totally cancelled under this circumstance, whereas the harmonics of the ±1st order reach the maximum amplitude with a1 = 0.6366. The calculated data in Fig. 7 were used to determine the contribution of the harmonics to overall reflection, which is dominated by the components k = ±1, accounting for 81.05 percent of the overall reflection energy in Fig. 7c. Such a proportion is especially favored for mixing applications, leading to new opportunities for up and down conversion with satisfactory efficiency in free space.

    Fig. 7  Theoretical harmonic amplitude and phase distributions.

    a-c The calculated harmonic amplitude and phase distributions under different combinations of φ12, which were set to 0°/36.87°, 0°/68.44° and 0°/180°, respectively
  • Until now, we have only considered the feasibility of harmonic generation through the time-domain digital coding metasurface. One drawback of the proposed scheme is the strong correlation between the harmonic amplitude and phase, thereby greatly hindering independent control of both parameters in practice. A possible recipe to overcome such limitation lies in the introduction of an additional time delay t0 to the time-varying reflection coefficient. Following the time shift property of a Fourier transform with $\Gamma (t - t_0) = {\cal F}^{ - 1}[e^{ - j\omega t_0}\Gamma (j\omega )]$, the time delay leads to an additional phase shift $e^{ - jk\omega _0t_0}$ for the kth order harmonic, while maintaining an unchanged amplitude. After a more formal derivation, we obtain the spectral expression for the reflected wave in this case

    $$ {{{E}}_r(\omega )} = {2\pi A{\kern 1pt} {{cos}}\frac{{\varphi _2 - \varphi _1}}{2}e^{j\frac{{\varphi _2 + \varphi _1}}{2}}{{E}}_i(\omega )} \\ \hskip 10pt + {\kern 1pt} \mathop {\sum}\limits_{m = - \infty }^{ + \infty } \frac{{4A}}{{2m - 1}}{\kern 1pt} {{sin}}{\kern 1pt} \frac{{\varphi _2 - \varphi _1}}{2}e^{j\left[ {\frac{{\varphi _2 + \varphi _1}}{2} - (2m - 1)\omega _0t_0} \right]} \\ \hskip 20pt {{E}}_i\left[ {\omega - (2m - 1)\omega _0} \right] $$ (7)

    Specifically, the amplitude and phase of the kth order harmonic are rewritten as

    $$ {{{E}}_r(k\omega _0 + \omega _c) = \left| {{A}} \right|\angle \Phi } \\ { = \left\{ {\begin{array}{*{20}{l}} {\left| {A{\kern 1pt} {{cos}}{\kern 1pt} \frac{{\varphi _2 - \varphi _1}}{2}} \right|e^{j\left\{ {\frac{{\varphi _2 + \varphi _1}}{2} + \pi \left[ {\varepsilon \left( {\frac{{\varphi _2 - \varphi _1 - \pi }}{2}} \right) + \varepsilon \left( {\frac{{\varphi _1 - \varphi _2 - \pi }}{2}} \right)} \right]} \right\}}} \hfill & {k = 0} \hfill \\ {\left| {\frac{{2A}}{{k\pi }}{\kern 1pt} {{sin}}{\kern 1pt} \frac{{\varphi _2 - \varphi _1}}{2}} \right|e^{j\left\{ {\frac{{\varphi _2 + \varphi _1}}{2} - \pi \left[ {\varepsilon (\varphi _2 - \varphi _1) + \varepsilon (k)} \right] - k\omega _0t_0} \right\}}} \hfill & {k = \pm 1, \pm 3, \pm 5 \cdots} \hfill \\ 0 \hfill & {k = \pm 2, \pm 4, \pm 6 \cdots} \hfill \end{array}} \right.} $$ (8)

    The extra phase shift provides new degrees of freedom to control the propagation features of high-order harmonics by breaking the constraint between their amplitude and phase, making it possible to realize independent regulation of both parameters, which is hard to achieve using conventional methods. In fact, Eq. (8) offers a promising route to tailor the response of reflected harmonics via the combination of three parameters φ1, φ2, and t0, consequently pushing the intriguing functionalities accomplished by traditional metasurfaces into the nonlinear regime.

    The concept of a digital coding metasurface arises from a simple analogy with digital circuits, where a finite number of element types are employed to manipulate electromagnetic waves. For a 1-bit metasurface, only two kinds of elements (named as '0' and '1' element) with opposite phase responses are utilized to constitute the entire structure. By alternating the spatial alignment of both elements, the metasurface demonstrates a powerful capability to reach a broad range of functionalities such as anomalous beam reflections and random diffusion. Such ideas can be easily extended to design of the metasurface in the time domain. For example, given the phase states φ1, φ2 in Eq. (8), it is easy to construct a '0' and '1' element by inserting a time delay t0 for the kth order harmonic. Here, we choose t0 = 0 and T/2 respectively, in which T is the period of the square function in Eq. (4). As a quantitative illustration of nonlinear control, we considered a scenario where a plane wave is illuminated normally at fc = 3.7GHz upon the 1-bit time-domain digital coding metasurface. In this case, we use eight columns of meta-atoms in total with a unit period P = 33 mm. The repeating frequency of the square wave was f0 = 156.25 kHz. Figure 5a, c, e shows the calculated E-plane scattering pattern for the +1st order harmonic under different coding sequences '0000…', '00001111…', and '00110011…'. From the antenna theory, the scattering pattern is shown to hold as a Fourier transform pair with the phase distributions within the metasurface. Therefore, the scattering pattern can be directly inferred from the coding sequences of the binary elements through mathematical calculations. The backward scattering beam is split into two symmetric beams directed along either side, with the tilting angle gradually increased, as shown in Figure 5a, c, e. Such scattering behavior is closely associated with the introduction of the parallel wave-vector that is imposed onto the reflected wave with the change of spatial arrangement for the '0' and '1' elements.

    In addition, the red, blue and green lines in Fig. 5a, c, e demonstrate the dependence of the scattering magnitude of the +1st order harmonic under different phase combinations of φ1/φ2 as 0°/180°, 0°/68.44° and 0°/36.87°, respectively. The scattering pattern undergoes a uniform attenuation in its amplitude by 0, 5, and 10 dB with the change in phase pair.

    One important concern related to the metasurface is the conversion efficiency from the fundamental harmonic to high-order harmonics. To address this issue, we compare the measured and theoretical conversion efficiency of the −4th to +4th order harmonics in Tables 2 and 3, with V1/V2(φ1/φ2)= V/12 V(0°/180°), T = 6.4 μs and V1/V2(φ1/φ2)= V/18 V (0°/68.44°), T = 6.4 μs, respectively. The efficiency is obtained from the energy ratio between the harmonic and incident wave. From Table 2, when φ1 and φ2 are out of phase, the measured fundamental wave occupies 3.49% of the total energy, indicating that 96.51% of the incident power is converted to high-order harmonics. The majority of the energy is assigned to ±1st harmonics, as expected. In general, the theoretical and measured results for each harmonic are in good accordance. In the measurement, when φ1φ2 = 68.44°, only 30% of the incident energy is converted into high-order harmonics, consistent with the theoretical prediction in Table 3.

    Harmonic order Measured V1/V2 = 0 V/12 V Theoretical φ1/φ2 = 0°/180° Error

    Table 2.  Measured and theoretical conversion efficiency when V1/V212) = 0 V/12 V (0°/180°) with T = 6.4 μs

    Harmonic order Measured V1/V2 = 9 V/18 V Theoretical φ1/φ2 = 0°/68.44° Error

    Table 3.  Measured and theoretical conversion efficiency when V1/V212) = 9 V/18 V (0°/68.44°) with T = 6.4 μs

  • This work is supported by the National Key Research and Development Program of China under Grant Nos. 2017YFA0700201, 2017YFA0700202, 2017YFA0700203, the National Science Foundation of China (61631007, 61571117, 61138001, 61371035, 61722106, 61731010, 11227904), and the 111 Project (111-2-05).

Author contributions
  • J.Y.D., J.Z., Q.C., and T.J.C. contributed equally to this work. J.Y.D. carried out the analytical modeling, numerical simulations, sample fabrication, and measurements together with J.Z. Q.C. and T.J.C. conceived the idea, suggested the designs, planned, coordinated, and supervised the work. All authors discussed the theoretical and numerical aspects and interpreted the results. All authors contributed to the preparation and writing of the manuscript.

Conflict of interest
  • The authors declare that they have no conflict of interest.

  • Reference (47)



      DownLoad:  Full-Size Img PowerPoint