Issue
EPJ Appl. Metamat.
Volume 8, 2021
Metamaterials Research and Development in China
Article Number 15
Number of page(s) 4
DOI https://doi.org/10.1051/epjam/2021006
Published online 11 June 2021

© S. Yang et al., Published by EDP Sciences, 2021

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Metamaterials, normally composed by artificially periodic structures with sub-wavelength scales, made it possible to design multifunction coding metamaterials [13]. Various of extraordinary EM phenomena have been found in microwave, terahertz and optical fields [47]. Consequently, coding metamaterials have tremendous application potentials, numerous devices like invisibility cloak, planar lens, and perfect absorber were made by using coding metamaterials [810]. Most of aforementioned metamaterials have a specific function, therefore the manipulation of EM wave is fixed once the design is completed. In recent years, much attention has been paid on active coding metamaterials whose EM response can be dynamically controlled [11,12]. While a variety of functionalities have been successfully presented with different metamaterials, the way to get coding matrix and simultaneous demonstration of different functionalities on a certain metamaterial has not been fully reported.

The presented study focuses on the optimization of the coding matrix based on machine learning. Initially, we presented two designs of metamaterial unit cells with π phase difference and we simulated the reflection performances using full-wave simulation software. Subsequently we arranged the metamaterial particles periodically to construct a metasurface. Finally, we get certain scattering properties by programming a coding matrix with the model which is constructed under the open-source machine learning framework of TensorFlow [13].

2 Construction of semianalytical model and tensor model

First we propose a purpose architecture to automatically model and optimize coding matrix. Figure 1 illustrates the flowchart of the Machine learning based optimization. In what follows, we study the scattering properties of our coding metamaterials for plane-wave illumination. As a meaningful observable to quantify the scattering properties, we consider the bistatic RCS: σ ( θ , φ ) = lim x 4 π r 2 | E s | 2 | E i | 2 (1)

Accordingly, we define the RCS ratio γ = max θ , φ σ M M ( θ , φ ) max θ , φ σ P E C ( θ , φ ) (2)where Ei and Es denote the incident and scattered electric fields respectively, (θ, φ) denote the observation direction in the far fields, and the subscripts “MM” and “PEC” indetify the metamaterials and PEC targets. In equation (2) requires the full-wave numerical solution of two 3D scattering problems. Nevertheless, the Semianalytical approximation that was developed in “Coding metamaterials, digital metamaterials and programmable metamaterials” [1]. And the directivity of metamaterials can be described as D ( θ , φ ) = 4 π | f ( θ , φ ) | 2 | F | 2 0 2 π d φ 0 π 2 d θ | f ( θ , φ ) | 2 | F | 2 (3)where F = F (k sin θ cos φ, k sin θ sin φ). The expression in equation (3) closely resembles the directivity of antenna arrays. So, we considered every unitcells as an antenna.

The proposed machine-learning model for designing coding metamaterials is schematically depicted in Figure 2. As indicated by the orange zone, the model contains tensor module and semianalytical module, where both modules allowing data to flow towards to coding matrix and design pattern. The dataflow in the machine-learning is denoted by the arrows.

thumbnail Fig. 1

Flowchart of optimization based generic algorithm.

thumbnail Fig. 2

Structure of the machine-learning model for designing coding metamaterials.

3 Tailoring the scattering properties

After analysis by semi-analytical model, we designed the unitcells to verify the model. The unitcells are composed of three-layer dielectric substrate with double metallic split-rings structure and full copper ground plane on the bottom. The dielectric substrate can be described as ϵr = 2.65, tan δ = 0.001, the radius of each ring is 2.5, 3.8, 4.1, 4.6 mm, respectively. The phase differences are approximately equal to π in the range of 6–18 GHz as shown in Figure 3. So, we can use the unitcells to build a metamaterial according to the coding matrix T.

In order to achieve the various functionalities, we used machine-learning model that shown in Figure 2 to optimize the coding matrix. We combined the scattering properties of unit cell and the semianalytical modeling and get the relationship between scattering properties and coding matrix without using full-wave simulation.

For a reflective metamaterial with 8 × 8 lattices, each lattice has 2 × 2 unitcells, the different coding matrices with different information entropy and different corresponding scattering patterns are shown in Figure 4. First we define the coding matrix T as: T 1 = [ 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 ] (4)

Its corresponding scattering pattern shown in Figure 4a, we can find that the reflection of EM wave gathers in one direction.

Then we try to get the maximum of information entropy. Thus, T2: T 2 = [ 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 ] (5)

In T2 the prior probability of “0” equals to the prior probability of “1” (p = 0.5), Its corresponding scattering pattern shown in Figure 4b, we can find that the beam of reflection divided into three direction. In other words, while the larger information entropy, the smaller radar cross section.

Furthermore, by using machine-learning to optimize the coding matrix, the coding matrix staggered alone the x-direction, then we get the third coding matrix: T 3 = [ 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 ] (6)

Its corresponding scattering pattern shown in Figure 4c, the number of “0”' lattice in T2 is same as in T3, So the information entropy also gets the maximum value, the only difference is that the reflection is more similar to diffuse reflection in x-direction.

Finally, we optimized the coding matrix in both x-direction and y-direction, and we got the fourth coding matrix: T 4 = [ 1 0 1 0 1 0 1 0 0 1 0 1 0 1 0 1 1 0 1 0 1 0 1 0 0 1 0 1 0 1 0 1 1 0 1 0 1 0 1 0 0 1 0 1 0 1 0 1 1 0 1 0 1 0 1 0 0 1 0 1 0 1 0 1 ] (7)

Its corresponding scattering pattern shown in Figure 4d, compared to the results above, we can find that is a metamaterial with diffuse scattering. By using machine-learning, we encode and initialize a population, and calculated its information entropy and calculated the pattern based on the semianalytical model. Compared to the full-wave numerical simulation, the semianalytical method is easier to handle.

We used information entropy theory to describe the metamaterial. We defined H ( X ) = H ( p 0 , p 1 ) = i = 0 1 p i l o g p i where pi is probability of coding matrix X takes xi. According to the entropy theory we can find that the H (T1) = 0.8113, H (T2,3,4) = 1. Figure 4 presents the simulation results of scattering properties. We can describe the RCS reduction qualitatively based on the number of beam dispersions.

thumbnail Fig. 3

Geometry of the unitcells and their reflective parameters. (a) Geometry of “0” lattice; (b) geometry of “1” lattice; (c) S-parameters of unitcells; (d) phase difference between “0” and “1” unitcells.

thumbnail Fig. 4

Simulated scattering pattern of different coding matrices.

4 Conclusions

In conclusion, we have investigated the design of coding metamaterials for tailor scattering properties by optimizing the coding matrix. The above results and predictions rely on semianalytical modeling and full-wave numerical simulations. They provide a way to design multifunction coding metamaterials, this may enable to approach the absolute lower RCS or shaping the EM wave. the information entropy theory and genetic optimization algorithm were used for optimized the numbers and combinational format of unitcells. From the comparison, we can tailor scattering properties by coding matrix. In addition to the reduction of RCS this method may also find more general applications to complex field.

5 Implications and Influences

Numerical simulation is employed to generate the scattering properties data by CST Microwave Studio. The scattering properties of interest are set in the microwave region from 6 to 18 GHz. We have collected scattering properties for training and testing. The model is constructed under the open-source machine learning framework of TensorFlow.

References

  1. T.J. Cui, M.Q. Qi, X. Wan, J. Zhao, Q. Cheng, Coding metamaterials, digital metamaterials and programmable metamaterials, Light Sci. Appl. 3, e218 (2014) [CrossRef] [Google Scholar]
  2. Q. Zheng, Y. Li, J. Zhang, H. Ma, J. Wang, Y. Pang, Y. Han, S. Sui, Y. Shen, H. Chen, S. Qu, Wideband, wide-angle coding phase gradient metasurfaces based on pancharatnamberry phase, Sci. Rep. 7, 43543 (2017) [Google Scholar]
  3. Q. Zhang, X. Wan, S. Liu, J.Y. Yin, L. Zhang, T.J. Cui, Shaping electromagnetic waves using software-automatically-designed metasurfaces, Sci. Rep. 7, 3588 (2017) [Google Scholar]
  4. H.T. Wu, S. Liu, X. Wan, L. Zhang, D. Wang, L.L. Li, T.J. Cui, Controlling energy radiations of electromagnetic waves via frequency coding metamaterials, Adv. Sci. 4, 12 (2017) [Google Scholar]
  5. J. Su, Y. Lu, H. Zhang, Z. Li, Y. Yang, Y. Che, K. Qi, Ultra-wideband, wide angle and polarization-insensitive specular reflection reduction by metasurface based on parameter radjustable meta-atoms, Sci. Rep. 7, 42283 (2017) [Google Scholar]
  6. Y. Shen, J. Zhang, Y. Pang, Y. Li, Q. Zheng, J. Wang, H. Ma, S. Qu, Broadband reflectionless metamaterials with customizable absorption-transmission-integrated performance, Appl. Phys. A 123, 530 (2017) [Google Scholar]
  7. M. Moccia, S. Liu, R.Y. Wu, G. Castaldi, A. Andreone, T.J. Cui, V. Galdi, Coding metasurfaces for diffuse scattering: Scaling laws, bounds, and suboptimal design, Adv. Opt. Mater. 5, 1700455 (2017) [CrossRef] [Google Scholar]
  8. T.J. Cui, S. Liu, L. Zhang, Information metamaterials and metasurfaces, J. Mater. Chem. C 5, 3644 (2017) [CrossRef] [Google Scholar]
  9. K. Chen, Y. Feng, F. Monticone, J. Zhao, B. Zhu, T. Jiang, L. Zhang, Y. Kim, X. Ding, S. Zhang, A. Alu, C.-W. Qiu, A reconfigurable active Huygens' metalens, Adv. Mater. 29, 1606422 (2017) [Google Scholar]
  10. K. Chen, L. Cui, Y. Feng, J. Zhao, T. Jiang, B. Zhu, Coding metasurface for broadband microwave scattering reduction with optical transparency, Opt. Express 25, 5571 (2017) [CrossRef] [Google Scholar]
  11. H. Yang, X. Cao, F. Yang, J. Gao, S. Xu, M. Li, X. Chen, Y. Zhao, Y. Zheng, S. Li, A programmable metasurface with dynamic polarization, scattering and focusing control, Sci. Rep. 6, 35692 (2016) [CrossRef] [Google Scholar]
  12. C. Huang, W. Pan, X. Ma, X. Luo, Multispectral metasurface for different functional control of reflection waves, Sci. Rep. 6, 23291 (2016) [Google Scholar]
  13. K. Chen, Y. Feng, Z. Yang, L. Cui, J. Zhao, B. Zhu, T. Jiang, Geometric phase coded metasurface: from polarization dependent directive electromagnetic wave scattering to diffusion-like scattering, Sci. Rep. 6, 35968 (2016) [CrossRef] [Google Scholar]

Cite this article as: Shuai Yang, Kuang Zhang, Xumin Ding, Guohui Yang, Qun Wu, Tailoring the scattering properties of coding metamaterials based on machine learning, EPJ Appl. Metamat. 8, 15 (2021)

All Figures

thumbnail Fig. 1

Flowchart of optimization based generic algorithm.

In the text
thumbnail Fig. 2

Structure of the machine-learning model for designing coding metamaterials.

In the text
thumbnail Fig. 3

Geometry of the unitcells and their reflective parameters. (a) Geometry of “0” lattice; (b) geometry of “1” lattice; (c) S-parameters of unitcells; (d) phase difference between “0” and “1” unitcells.

In the text
thumbnail Fig. 4

Simulated scattering pattern of different coding matrices.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.