- Photonics Research
- Vol. 12, Issue 1, 61 (2024)
Abstract
1. INTRODUCTION
A metasurface is a special metamaterial two-dimensional plane with great potential for electromagnetic manipulation [1–6]. Its main role is to artificially manipulate several parameters such as amplitude, phase, frequency, and polarization of incident electromagnetic waves. In 2014, the concept of programmable metasurfaces was introduced by Cui
2. THEORY OF THE GERCHBERG-SAXTON ALGORITHM
The key to multi-focus stochastic coding of the metasurface is to determine the phase pattern of each cell, which can be equated to the near-field superposition problem [32–37]. Simply put, for a given desired field amplitude, we determine the phase of the traces traced. We use the parallel vector Green’s function as the kernel of operation to calculate the electric field distribution of the metasurface in space:
is excited by the surface current , from which the electric field strength in space can be calculated as
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
We define the source as discrete N points and match them to the near field of each of the M points. is denoted as the current vector, and is denoted as the vector field. Based on the above equations, Eq. (3) can be simplified as
Writing the scalars and for all and in the vectors and , Eq. (4) can be simplified to
The above method transforms the design of a metasurface into an arrangement of equivalent dipole arrays and uses a vector Green’s function as the propagation kernel, which can be used in both the far and near fields. We consider the metasurface as an array of equivalent dipoles and use a rigorous concatenated vector Green’s function as the propagation kernel. It breaks the limitations of Fresnel and Fraunhofer integrals, which are commonly used in the design of far-field metasurface holograms. Importantly, the method can be used to generate complex and high-resolution field distributions in regions adjacent to the metasurface [38–40]. A particular feature of the method is that, when evaluating fields far from the source, a simple expression for the field quantity can be used, where R is much larger than the typical dimension of the source. The vector Green’s function degenerates as
The specific design implementation process is as follows: first read in the initial map, adjust it to the same size as the resolution of the modulator, and normalize it; then generate random phase and initial complex amplitude distributions as the initial distributions of the sum GS algorithm; next, iterate to adjust the algorithm’s convergence speed and accuracy by setting the maximum number of iterations and feedback parameter, and each iteration calculates the inverse adjustment according to the error to update the distribution.
3. METASURFACE CELL DESIGN
We first design the 1-bit encoded metasurface [41–51]. We design a three-layer cell structure, with the top and bottom layers consisting of 0.018 mm copper metal. The bottom metal is a whole piece of copper to ensure the reflection efficiency of the metasurface. The top metal is a piece of specially designed topology that has excellent maneuverability for cell phase control. The intermediate medium is a high-frequency sheet Teflon ZYF255DA (, ), which has a thickness of 1.52 mm. Figures 1(a) and 1(b) show all the details of the cell structure. According to the need, we designed a rectangular array cycle with one PIN diode (SMP1320) referenced to each cell, totaling 480 pieces. As shown in Fig. 1(c), when the diode is off, the reflection amplitude tends to be , and when the diode is on, the reflection amplitude is . The difference in the reflection amplitude between the two states is within 1.5 dB, which satisfies the GS theory of near-field random coding. Figure 1(d) shows the phase of the Y-polarized light under the incidence of the Y-polarized light at the resonance frequency , with a phase difference of 180° between the on and off states of the diode. From the surface electric field distribution, it is seen that, under the action of this topology, the electric field energy will be mainly concentrated in the middle two metal gaps, whose action can be equated to a capacitor with the ability of energy aggregation, as in Fig. 1(e). Finally, we calculate the RLC equivalent model of the diode (SMP1320) in the ADS software, with details shown in Fig. 1(f).
Figure 1.(a), (b) Schematic of the cell structure, (c) reflection amplitude, (d) reflection phase, (e) surface electric field distribution, and (f) RLC model of the diode in the “ON” and “OFF” states.
To verify the feasibility of the design theory, we calculate the array arrangement of the target field with the aid of MATLAB tools. Based on the consideration of algorithm accuracy, we first verify the multifocal coding control. According to the above theoretical derivation, we first initialize the elements of and , defining their amplitude as the formulated value and phase as the phase at the frequency of 4.611 GHz in the Fig. 1(d) diagram. Letting , we calculate the forward propagation [Eq. (5)], and after keeping the phase value of the element of , we replace the corresponding magnitude with the magnitude of . Next, the backward propagation [Eq. (6)] is computed, keeping the phase of the elements of unchanged and replacing the corresponding magnitude with the magnitude of .
4. ENCODING METASURFACE NEAR-FIELD MULTIFOCAL FOCUSING
As shown in Fig. 2(a), we process the to-be-focused focus map into a grayscale image before operating. After several iterations, we can obtain the phase distribution image corresponding to the focus in Fig. 2(b). Yellow represents the on state of the diode, purple represents the off state of the diode, and the two colors indicate two different states with constant amplitude.
Figure 2.(a) Grayscale map of each focal point; (b) phase alignment map obtained from GS algorithm calculation.
We utilize the finite integration method to create arrays of size . And we encode the array with the calculation results of Fig. 2(b). We utilize numerical simulation to verify the hologram properties. When the Y-polarized plane wave is excited, a focused spot is formed at the position from the metasurface, as in Fig. 3. The spot is displayed at the corresponding place of the grayscale map just like the theoretical result, which further verifies the feasibility of the theory and algorithm.
Figure 3.Near-field focusing characteristics of the focal spot for each phase arrangement.
5. CONVOLUTIONAL CODING THEORY
In the next step, in order to flexibly regulate the focusing spot, we will introduce the Fourier convolution principle to the phase arrangement of the encoded metasurface. First of all, we need to clarify the Fourier transform relationship between the electric field distribution on the encoded metasurface and the far-field direction map, as shown in Eq. (8), which can be calculated to derive its encoded electric field distribution as
The Fourier transform principle is a method of converting a time-domain signal into a frequency-domain signal, which allows the original function to be represented in the form of multiple trigonometric functions so that parameters such as frequency can be more easily observed. Using this principle, we can relate the coded pattern in the coded metasurface to the far-field directional map, which allows us to understand the principle of coded metasurfaces from the perspective of signal processing. In addition, the convolution theorem of the Fourier transform can also be expressed in this way below:
We try to express in Eq. (11) in terms of and in terms of to obtain the following equation:
We could get further,
In the following, we combine the convolution theory with the convolution operation on the focus C. First, according to the pitch angle formula of generalized Snell’s law [40],
Figure 4.(a) Coding S1 with “focus C” convolution. (b) Coding S2 with “focus C” convolution. (c) Encoding S3 with “focus C” convolution operation.
After obtaining the above encoding map, we use the finite integration method to import it into the numerical simulation method. A suitable monitor is set up to measure the near-field electric field distribution. As shown in Fig. 5(a), there are two symmetric beams at a distance of 100 mm from the super-surface, which is due to the addition of the encoding of Fig. 4(a). Figure 5(b) is the Y-axis section, and we can see that the beams are uniformly spread out in two beams at 43°, which is 3° different from the theoretical value of 46°. Figures 5(d) and 5(e) are the S2 encoding of Figs. 5(d) and 5(e). The result after convolution with focus C has a beam splitting angle of 31°, which is by the theoretical calculation. Similarly, Fig. 5(g) is the corresponding encoding sequence calculated by Fig. 4(c), and its beam splitting angle is 25°; all three simulation results are within the acceptable range of the theoretical calculation error. The reason for the error we believe is that the metasurface array is so small.
Figure 5.(a)–(c) Near-field Z, Y, 3D maps obtained by S1 encoding convolution. (d)–(f) Near-field Z, Y, 3D plots obtained by S2 coded convolution. (g)–(i) Near-field Z, Y, 3D maps obtained by S3 encoded convolution.
6. NEAR-FIELD ALPHABET HOLOGRAPHY
The phase distribution map is obtained by iterating using the GS algorithm, which can theoretically obtain any shape graph in the near-field range. Its resolution varies with the number of metasurface array units and is related to the frequency of the band in which it is located concerning the near-field focal length. Due to the limitation of the resolution of the metasurface, only the four letters “C,” “J,” “L,” and “U” are used here. “C,” “J,” “L,” and “U”, which are four letters with relatively simple English character outlines, are used as examples to verify the capability of the metasurface reconstruction. As shown in Fig. 6(a), the phase distribution of the letter “C” is calculated by the GS algorithm. The yellow color block represents the diode on state, and the purple color represents the diode off state. The obtained phase distribution and the corresponding near-field electric field hologram are displayed in Fig. 6.
Figure 6.(a) Phase distribution and near-field hologram corresponding to the letter C, (b) phase distribution and near-field hologram corresponding to the letter J, (c) phase distribution and near-field hologram corresponding to the letter L, and (d) phase distribution and near-field hologram corresponding to the letter U.
7. EXPERIMENTAL PREPARATION AND TESTING
We recorded the specific coding sequence into the FPGA, increased the number of pins through the pcf8575 expansion module, and finally accessed it to the metasurface. Figure 7(a) shows the flowchart of the whole test. It is worth mentioning that, to solve the problem of insufficient FPGA pins, we made an IO port expansion module with the chip model pcf8575, which works by utilizing the IIC bus protocol and has a 7-bit address bit input signal to ensure that each diode can be controlled individually. The bottom of the front of the ultra-surface is designed with two sets of in-line pins, which are connected to each ultra-surface in turn through wires and finally connected to the ground terminal of the FPGA to form an electrically closed loop. We connect a vector network analyzer (VNA), model Agilent E5071C, to the PC and control it with MATLAB software. Ports S1 and S2 are connected to an excitation horn (HD-70SGAH10) and a waveguide (HD-48WOEWPN), respectively. The other port of the PC is connected to the FPGA for burning the coded sequences. pcf8575’s expansion pins are connected to the DC port of the super surface. Finally, we use a vector network analyzer to transmit microwave signals to the metasurface and receive them from the other end speaker. Figure 7(b) shows the connection of some of the active devices, where the voltage source, the metasurface, the FPGA, and the pcf8575 share a common VCC and GND. The voltage source serves to help the FPGA to stabilize the output of the 3.3 V pin voltage. SCL and SDA are the clock and data lines of the IIC bus, respectively, to help the FPGA. To facilitate power feeding, the metasurface consists of five small sub-blocks arranged in a array. The experimental overview diagram is shown in Fig. 7(c). When a plane wave polarized in the Y direction is shot into the metasurface, a holographic image is generated within about 100 mm of the front surface of the metasurface. The PCB is drawn by Altium Designer software as shown in Fig. 7(d).
Figure 7.(a) Block diagram of the experimental flow, (b) electrical connection diagram, (c) experimental overview diagram, and (d) 3D rendering of a single PCB board on the super surface.
The waveguide is fixed to a three-axis bracket consisting of root aluminum tubes, which enables three-axis control, as shown in Figs. 8(a)–8(c). Each aluminum rod is equipped with a slide rail and is connected to a servo motor. The three aluminum rods are orthogonal to each other, which allows the waveguide horn to move spatially in three axes XYZ. We then prepare the metasurface samples by standard PCB techniques. Figure 8(d) shows the metasurface sample with 0.018 mm of copper metal laminated on the front side, and the thickness of the substrate is 1.52 mm. ZYF255DA (, ) is used as the substrate material. Each unit is fed by a copper post with a diameter of 0.4 mm, and copper posts with a diameter of 0.6 mm are punched at the top and bottom of the super surface for fixing the acrylic board on the backside. Six 16-pin SMD PCF connectors are lined up at the top of the back side of the PCB, which is connected to the expansion module through the flexible rows. During the measurement, the probe was scanned over the observation area, and the S21 parameters as well as the probe position were recorded synchronously. The measured field is plotted directly with the electric field received by the probe in Fig. 8(a). Figure 8(e) shows an overview view of the measurement apparatus with a 3D model.
Figure 8.Diagram of the experimental setup.
Considering the width of the caliber of the waveguide horn, we decided to use a stepping for the measurement. As represented in Fig. 8(a), the servo motor stepping program was written in MATLAB, and the horizontal X axis was first moved in the positive direction 40 times, with a stepping distance of 1 cm each time. Subsequently, 1 cm was stepped in the negative direction of the Z axis, and the waveguide horn was stepped in the negative direction of the X axis another 40 times. This is repeated 11 times, and the waveguide records the electric field distribution at 100 mm of the whole metasurface according to the S-measurement. Therefore, the whole experiment ends up, getting matrix data as shown in Fig. 9. For example, as shown in Fig. 9(a), there is a focused spot in the upper left corner at the approximate matrix position (5, 11), which is consistent with the “focus A” result in the Fig. 3 simulation. Similarly, Fig. 9(b) also shows a peak in the electric field at (33, 11) out, corresponding to the second plot simulated in Fig. 3, and so on. The small positional deviation is due to the selection error of the initial position of the waveguide. Although there are some deviations, the measured field agrees very well with the simulation results. Figure 9(f) shows the beam splitting image at focus C and , which is denoted by Convolution 6 here. As shown, we can see that the spot is separated in the Y direction, in agreement with the simulation in Fig. 5(a). The last two plots show the electric field experiments for Convolution 8 and Convolution 10, respectively. However, the image is not clear enough because the energy is too dispersed and receives spurious interference.
Figure 9.Eight-group multifocal experimental validation plot.
8. CONCLUSION
We have proposed a class of “1-bit” metasurface units, which have a strong manipulation of unidirectional polarized light. In addition, we propose a GS algorithm with a near-field scalar Green’s function to construct holograms in the microwave band. The corresponding phase alignments are calculated to generate multifocal holograms and also have some near-field holographic imaging effects for simple letters. Eight sets of multifocal holograms are verified in simulations. The measurement results are in good agreement with the simulation results. In the experiments, the FPGA is clocked at 50 MHz, and a set of 16-bit data is encoded at a rate of roughly 39 ns (with waiting time reserved between each module). Our proposed encoding array is , 480 bits of data. The FPGA encodes the metasurface once at a rate of 1170 ns. It is worth mentioning that the reverse recovery time of a PIN diode is only a few nanoseconds, which means that the near-field image can be dynamically changed at a high rate by using the FPGA to frequently switch the encoding sequence of the metasurface. If the near-field scanning speed is improved, the breakthrough from “near-field hologram” to “near-field holographic animation” may be realized in the future.
References
Set citation alerts for the article
Please enter your email address