Broadband Microwave Coding Metamaterial Absorbers
Broadband Microwave Coding Metamaterial Absorbers
Broadband Microwave Coding Metamaterial Absorbers
com/scientificreports
Research involving metamaterial has been massively focused by many groups worldwide in recent decades1.
Metamaterial is an artificially structured material, generally made up of periodic metal circuits on a dielectric
substrate. This material was first theoretically predicted by Veselago in his paper on negative refractive index2.
But, almost four decades later, based on experimental evidence, the importance and applicability of metamate-
rial have been confirmed3,4. This material has a variety of electromagnetic and optical properties different from
ordinary materials, leading to many special applications such as super-lens, cloaking, wireless power transfer,
high-performance antenna systems1,4,5, and most recently the electromagnetic perfect absorber that the energy
of the incident wave is mostly absorbed2,6,7. This dernier special feature allows many potential applications such
as increased solar cell efficiency, military radar camouflage, super-sensitive sensor design and so on to be real-
ized8–18. For achieving efficient applications, current research focuses on two main issues: controlling the absorp-
tion band and finding ways to extend the working frequency range. Different studies have been carried out for
perfect absorber in the frequency domain from microwave to visible18–34. Various methods have been proposed
and yielded positive results such as multilayered, asymmetric, super-cell structure, all-metal or all-dielectric
structure, and hybrid structure34–46. Even now research has been geared up towards digital or coding metamate-
rial absorbers, which is related to computational metamaterial, to increase the applicability and meet the develop-
ment of digital technology47–59. However, the challenge is that experimental structures often have high complexity,
low flexibility, and therefore cause limitations in application.
Next, the study of the metamaterial absorber is performed generally from one unit cell with appropriate
boundary conditions in the simulation, where the structure is assumed to be infinite. With such a simulation
model, controlled defects in the structure due to its periodicity cannot be investigated. In this study, we consider
a full-sized structure, which resembles the structure in a realistic size. Therefore, the introduction of controlled
defects for absorption manipulation is feasible and realistic.
1
Faculty of Physics, Hanoi National University of Education, 136 Xuan Thuy, Cau Giay, Hanoi, Vietnam. 2Institute of
Physics, Vietnam Academy of Science and Technology, 18 Hoang Quoc Viet, Hanoi, Vietnam. 3Institute of Materials
Science, Vietnam Academy of Science and Technology, 18 Hoang Quoc Viet, Hanoi, Vietnam. 4Department of
Physics, Hanoi University of Mining and Geology, 18 Pho Vien, Bac Tu Liem, Hanoi, Vietnam. 5Graduate University of
Science and Technology, Vietnam Academy of Science and Technology, 18 Hoang Quoc Viet, Hanoi, Vietnam. *email:
[email protected]; [email protected]
Figure 1. (a) An elementary unit cell with its dimension, (b) A full-sized structure of 144 UCs.
Studying the application of defects in the full-sized absorber for broadening the working frequency band has
been described in our previous reports30,31. It can be seen that the method has many advantages such as simple
structure, high flexibility, easy to integrate into the electronic or optical circuits. The challenge to this method is
that the optimization process is quite complex and takes a lot of work. This work integrates the digital concept
into the metamaterial absorber system, namely coding metamaterials. Coding and digital metamaterials are a
new field of research and attract a lot of attention48,49. Researches often focus on investigating the ability to control
the characteristics of reflected beams with different coding configurations53–59. Here, in this study we focus on
investigating the electromagnetic wave absorption performance of the structure with different coding metasur-
faces through the reflection and transmission parameters. Furthermore, in this study we place more emphasis on
the ability to improve the performance of the meta-structure through a search algorithm using a large number of
different coding configurations.
For the calculation of the absorption, 200 samples are used in each approach, we obtain the optimal config-
uration with a broadband up to 17 GHz, the results then theoretically analyzed and experimentally confirmed.
This approach is expected opening up a new way to study and design broadband or multiband absorber metama-
terial. Moreover, this method promises to fit very well with studies of the actual size structure, which requiring
pre-analysis of electromagnetic response via simulation.
Combination process. In the first approach, named SRS (simple random sample) method, to ensure the
rationality of the full-sized structure and adapt to our calculation system, the configuration of 144 UCs (12 × 12
UCs) is investigated. The idea is that the defect positions or bits 0 will be generated by randomly removing the top
copper layer of the UCs at an arbitrary position on the configuration, the locations with a normal UC are assigned
as bits 1. In principle, the structure of n × n UCs could generate 2n×n configurations. The structure of 144 UCs
can generate therefore 2144 configurations and the possible number of simulations is very large. Each generated
configuration will be transferred to the environment with appropriate boundary conditions to simulate the elec-
tromagnetic absorption response. The final results are aggregated and analyzed to provide optimal configuration.
It is worth mentioning that a simple random sample is a subset of a statistical population in which each mem-
ber of the subset has an equal probability of being chosen. As an example, we draw multiple samples consisting
of 200 configurations each from a set of possible 2144 configuration (configurational space). Here, the variable s is
used to represent the size of the sample; thus configuration space size N = 2144 and sample size s = 200). By ran-
domizing the selection procedure, any member of this configuration space has an equal chance of being selected
as part of this first sample, and an equal chance of being selected for the next sample of the same size (and so on).
After the simulation process, the absorption spectra of 5 selective cases (SRS 12 × 12–1 to SRS 12 × 12–5) are
shown for a convenience view in Fig. 2.
Figure 2. Absorption spectrum of SRS configurations and of the full-sized structure of 144 UCs. Right panels
are different configurations SRS 12 × 12–2 to SRS 12 × 12–5.
One can see that a wide absorption band of ~17 GHz from 16 to 33 GHz (with the average absorptivity is about
85%) is obtained for almost of cases. We also present the absorption of the non-defect full-sized absorber with
144 UCs for the comparison (Fig. 2). The absorption range is improved for the case of applying defect generation
algorithms in the structure.
This is a favorable result and can be used to confirm that by randomly generate the defects on the structure, a
remarkable broad absorption frequency band can be obtained.
However, it is easy to point out that the absorption band for most random configurations has a deep recess
around the 19 GHz and 28 GHz region. This suggests us to find out a more optimal strategy that can improve the
absorption strength in the range. The next section will discuss this.
The second approach is named CFM (combination of fundamental meta-block), based on the combination
process of the basic meta-blocks (MBs). MB is a set of basic blocks with the number of UC greater than 1, in
this study, MBs are created from blocks including 2 × 2, 3 × 3, 4 × 4 and 6 × 6 UCs, which are the divisors of 12.
Derived from the idea that relevant absorption peaks can overlap and form a broader absorption band, first we
calculate the absorption of the different MBs, then the structures can combine to have the broadband selected to
build a full-sized structure (of 144 UCs).
The first case discussed is MB of 4 UCs (called MB2 × 2), for the 2 × 2 block, there are 16 distinct configu-
rations corresponding to different combinations of 0 and 1 bits. The absorption results of 5 cases (MB2 × 2–1
to MB2 × 2–5) are shown in Fig. 3(a) for the illustration and two best selective cases (in bold curves) are cho-
sen for the combinational process to generate the full-sized structure. The corresponding two best configurations
(MB2 × 2–1 and MB2 × 2–2) are also illustrated in the inset of Fig. 3(a). In the same scenario, the simulation results
of different types of MB as MB3 × 3, MB4 × 4, and MB6 × 6 are presented in Fig. 3(b–d). For each type of MB, we
choose two best configurations (which own two best absorption curves in bold) to combine the full-sized struc-
ture. E.g. for MB3 × 3 type, MB3 × 3–1 and MB3 × 3–2 are chosen for the combination of a full-sized structure.
From the coding metamaterial perspective, a meta-block could be treated as a ‘super-bit’, where the bit 0 or 1 is
named by two best configurations of MB, for example, MB2 × 2–1 is set as bit 2–1 (for ON state) and MB2 × 2–2
is set as bit 2–0 (for OFF state), that is why in the insets of Fig. 3, we set the name ‘bit MB’ for the configurations.
By choosing the suitable ON/OFF state of this super-bit on a surface of 144 UCs, we can optimize the structure to
the widest absorption band by combining the MB or bits in the most convenient position. This is also the way that
our algorithm performs in the simulation process to combine MBs (or ‘super-bit’) to find the optimal full-sized
absorber structure. According to the same principle, we call the MB3 × 3–1 and MB3 × 3–2 are bit 3–1 and 3–0;
MB4 × 4–1 and MB4 × 4–2 are bit 4–1 and 4–0; MB6 × 6–1 and MB6 × 6–2 are bit 6–1 and 6–0, respectively.
These bit digits are illustrated in the inset of Fig. 3. When applied in practice, by integrating the ON/OFF circuit
into the location of each MB, we can easily implement the desired full-sized configurations through a designed
digital controller.
In the next step, each type of MB will be used to combine the full-sized configuration of 144 UCs. With
MB2 × 2 type, 36 MBs are employed to combine the full-sized structure. The number of possibilities of simula-
tions for this type is 236. As an illustration of a straightforward approach and similar to the discussion in the SRS
method, we calculate the absorption spectrum for a number of different structures, i.e. 200 structures and the
absorptions of 5 selective cases are shown in Fig. 4(a). The best full-sized configuration, in this case, is assigned
as CFM2 × 2–1 and is illustrated in the inset of Fig. 4(a). We see that this combination process of MB2 × 2–1 and
MB2 × 2–2 does not greatly improve the absorption of the structure.
Figure 3. Absorption spectrum of different fundamental MB types: (a) MB2 × 2, (b) MB3 × 3, (c) MB4 × 4 and
(d) MB6 × 6.
With MB3 × 3 type, we need 16 MBs to combine the full-sized structure of 144 UCs. The number of inves-
tigated simulation is 216. As in the previous case, we perform 200 simulation rounds and 5 selective cases are
shown in Fig. 4(b). The MB3 × 3-based full-sized structures show that there exists an optimal configuration with
a superior result with the average absorptivity in the frequency range of 16–33 GHz is over 90% and the maximal
absorption reaches 98%, this full-sized configuration is assigned as CFM3 × 3–1 and is illustrated in the inset of
Fig. 4(b). In this case, the absorption of the full-sized configuration has been significantly improved compared to
the SRS method.
The scenario is similar for the combination using the MB4 × 4 and MB6 × 6. The structure of 4 × 4 type needs
9 MBs then 29 simulations could be performed. The 6 × 6 type requires 4 MBs so only 24 simulations are needed
for the investigation. The simulations are then performed and the absorption of five selective cases of each type
is reported in Fig. 4(c,d). The inset of Fig. 4(c,d) also shows two best full-sized configurations CFM4 × 4–1 and
CFM6 × 6–1, respectively.
Thus, the structures show that they can be optimized for a wide absorption frequency range from 16–33 GHz
with all case of MB, but the MB3 × 3 types (configuration CFM3 × 3–1) shows the most effective improvement
of the CFM method compared to the SRS method. A more detailed description of the simulation method and
additional results can be found in Supporting Information.
Figure 4. Absorption spectrum of full-sized absorbers after the combination of: (a) MB2 × 2 type, (b) MB3 × 3
type, (c) MB4 × 4 type and (d) MB6 × 6 type. The two types of colors in each full-sized configuration illustrate
two different types of bits.
Experiment Verification
After the simulation, a full-sized non-defect structure, an optimized full-sized structure by SRS method (SRS
12 × 12–1), four most optimal absorber structures by CFM method (CFM 2 × 2–1, CFM 3 × 3–1, CFM 4 × 4–1,
CFM 6 × 6–1) are fabricated and measured the absorption (Fig. 10). All samples have a size of 12 × 12 UC equiv-
alent to an area of ~11 cm2. In the measurement, a vector network analyzer (Rohde & Schwarz ZNB20) with the
working frequency range from 100 kHz to 20 GHz is used. Two standard horn antennas working at the same fre-
quency range are applied for the electromagnetic emission and reception. The absorption rate is deduced from the
reflection coefficient. Due to the limitation of our equipment, we perform the measurement only in the range of
13–20 GHz and the comparisons are also in this frequency range. It is clear that from 15 to 20 GHz, the measured
results are in good agreement with the simulated ones for all cases (Fig. 9). From the tendency of the experimental
results, we believe that the measurements will adapt to the simulations if equipment with a wider measuring range
is applied. This also confirms the correctness of our method.
Figure 5. (a) The simulated power loss energy density on SRS and CFM structures at 25 GHz, which is in the
absorption frequency range of the full-sized absorbers; (b) The simulated electric field density on SRS and CFM
structures at 25 GHz.
Figure 6. Simulation of the surface current of the unit cells at the absorption frequency of CFM3 × 3–1.
Additionally, the measured absorption spectra for different incident electric field polarization angles of all
prototypes are shown in the frequency range from 13–20 GHz. As seen in Fig. 11, the absorption spectra of SRS
and CFM structure depend weakly on the polarization of the incident wave until 45 degrees.
Conclusions
This report presents a new study on broadband metamaterial absorber based on defect optimization coding of
the full-sized metamaterial structure. Two approaches are presented for considering the broadband character-
istic of the material. First, a full-sized structure is used to investigate the absorption when randomly removing
the UC in the configuration. The second strategy uses the combination of meta-blocks of 2 × 2, 3 × 3, 4 × 4 and
Figure 7. EMT model used for the calculation of the theory absorption.
Figure 8. Relative impedance of the full-sized absorber types: (a) Full 144 UCs, (b) SRS12 × 12–1 and (c)
CFM3 × 3–1.
6 × 6 UC to control and improve the absorption rate in the working frequency range. By integrating the digital
concept into the metamaterial absorber system, namely coding metamaterials, the calculation, and analysis of
the electromagnetic responses are significantly improved. A broadband absorber could be obtained in all cases
when the optimal defect induced. The EMT model is then applied for understanding the absorption mechanism.
Samples are fabricated for validating experimentally the proposed structures. The results show that the optimized
absorbers work steadily in the 16–33 GHz range with high absorptivity up to 98%. The full-sized metamaterial
absorbers designed have great potential in the community of anti-detector for the radar system, electromagnetic
energy harvesting and so on. The presented method could pave a new road for the study and realization of the
broadband coding metamaterial absorbers.
Methods
Simulation. Different configurations are generated using our computational algorithm which associated with
the simulator. Each configuration corresponds to a two-dimensional logic frame of n × n unit cells where bit “0”
represents the cell without metal plate, namely defect state, and bit “1” represents the cell with the metal plate, n
is the number of row or column of the full-sized structure. Simulations are then performed using the commercial
simulator CST Studio Suite. We use the open boundary conditions in the investigation process. A waveguide
Figure 9. Comparison the simulation, theory and measured absorption spectrum of (a) full 144 UCs structure,
(b) SRS12 × 12–1 structure and (c) CFM3 × 3–1 structure.
Figure 10. Image of the fabricated samples using for the measurement.
Figure 11. Measured angular dependence of the absorption of different full-sized absorbers under incident
electromagnetic waves.
port is used for the excitation; the absorption is calculated from the magnitude of S11 parameters. All results are
collected and analyzed afterward.
Theory calculation. As in30, the final absorptivity is calculated through the total reflection (Rtot) as A=1−Rtot
while the total reflection is given from the real part and imaginary part of S11 and S21 (the transmission equals
zero by using a metal plate at the back of the structure), which are extracted from the numerical simulation. Our
Matlab code was used for data processing.
Measurements. Measurements have been realized using a vector network analyzer R&S ZNB20 and a pair of
standard horn antennas (emission and reception), working in the range of 100 kHz–20 GHz. The reflection coeffi-
cient is normalized using a full metal reflector at the same sample position. The angle of incidence wave is altered
by changing the position of the horn antennas around the sample surface normal. The experimental results are
compared with the simulations.
References
1. Engheta, N. and Ziolkowski, R. W. Metamaterials: Physics and Engineering Explorations, Wiley-IEEE Press (2006).
2. Veselago, V. G. The electrodynamics of substances with simultaneously negative values of ε and μ. Sov. Phys. Usp. 105, 509 (1968).
3. Watts, C. M., Liu, X. & Padilla, W. J. Metamaterial Electromagnetic Wave Absorbers. Adv. Mater. 24, 23 (2012).
4. Tran, C.-M. et al. High impedance surfaces based antennas for high data rate communications at 40 GHz. Prog. Electromagn. Res. 13,
217–229 (2010).
5. Kundtz, N. & Smith, D. R. Extreme-angle broadband metamaterial lens. Nat. Mater. 9, 129–132 (2010).
6. Cuong, T. M., Thuy, N. T. & Tuan, L. A. Influence of Structural Parameters on a Novel Metamaterial Absorber Structure at K-band
Frequency. J. Electron. Mater. 45(5), 2591–2596 (2016).
7. Mehmet, E. et al. Investigation of electromagnetic and ultraviolet properties of nano-metal-coated textile surfaces. Appl. Nanosci.
(2019). https://doi.org/10.1007/s13204-019-01122-1.
8. Cui, Y. et al. Ultrabroadband Light Absorption by a Sawtooth Anisotropic Metamaterial Slab. Nano Lett. 12(3), 1443 (2012).
9. Chen, X. et al. Dual-polarity plasmonic metalens for visible light. Nat. Commun. 3, 1198 (2012).
10. Aydin, K., Ferry, V. E., Briggs, R. M. & Atwater, H. A. Broadband polarization-independent resonant light absorption using ultrathin
plasmonic super absorbers. Nat. Commun. 2, 517, https://doi.org/10.1038/ncomms1528 (2011).
11. Gao, Q., Yin, Y., Yan, B. & Yuan, N. C. Application of metamaterials to ultra-thin radar-absorbing material design. Electron. Lett.
41(17), 936–937 (2005).
12. Lee, Y., Kim, S.-J., Park, H. I. D. & Lee, B. Metamaterials and Metasurfaces for Sensor Applications. Sensors 17, 1726, https://doi.
org/10.3390/s17081726 (2017).
13. Mehmet, B. et al. A Comprehensive Study on Fuel Adulteration Sensing by Using Triple Ring Resonator Type. J. Electrochem. Soc.
166(12), B1044–B1052 (2019).
14. Abdulkarim, Y.I., Deng, L.W., Karaaslan, M., Unal, E. Determination of the liquid chemicals depending on the electrical
characteristics by using metamaterial absorber based sensor. Chem. Phys. Lett, 732, Article 136655 (2019).
15. Fatih, O. A. et al. Enhancement of image quality by using metamaterial inspired energy harvester. Phys. Lett. A 384(1), 126041
(2020).
16. Murat, O. et al. Experimental work on mechanical, electromagnetic and microwave shielding effectiveness properties of mortar
containing electric arc furnace slag. Constr. Build. Mater 165, 58–63 (2018).
17. Liu, N. et al. Infrared Perfect Absorber and Its Application As Plasmonic Sensor. Nano Lett. 10, 2342 (2010).
18. Khuyen, B. X. et al. Ultra-subwavelength thickness for dual/triple-band metamaterial absorber at very low frequency. Sci. Rep 8,
11632 (2018).
19. Heijun, J. & Sungjoon, L. Broadband frequency reconfigurable metamaterial absorber using switchable ground plane. Sci. Rep 8,
9226 (2018).
20. Hai, L. D. et al. Conductive polymer for ultra-broadband, wide-angle, and polarization-insensitive metamaterial perfect absorber.
Opt. Express 26, 33253 (2018).
21. Cheng, Y. Z. & Yang, H. L. Design, simulation, and measurement of metamaterial absorber. J. Appl. Phys. 108, 034906 (2010).
22. Yao, G., Ling, F., Yue, J., Luo, C. & Ji, J. and Yao, Dual-band tunable perfect metamaterial absorber in the THz range. J. Opt. Express.
24(2), 1518, https://doi.org/10.1364/OE.24.001518 (2016).
23. Long, C. et al. Broadening the absorption bandwidth of metamaterial absorbers by transverse magnetic harmonics of 210 mode. Sci.
Rep. 6, 21431 (2016).
24. Jianfei, Z. et al. Ultra-broadband terahertz metamaterial absorber. Appl. Phys. Lett. 105, 021102 (2014).
25. Dang, H. L. et al. Broadband metamaterial perfect absorber obtained by coupling effect. J. Nonlinear Opt. Phys. Mater 26, 3 (2017).
26. Yuping, Z. et al. Independently tunable dualband perfect absorber based on graphene at mid-infrared frequencies. Sci. Rep. 25,
18463 (2015).
27. Lin, C. H., Chern, R. L. & Lin, H. Y. Polarization-independent broad-band nearly perfect absorbers in the visible regime. Opt.
Express 19(2), 415–424 (2011).
28. Avitzour, Y., Urzhumov, Y. A. & Shvets, G. Wide-angle infrared absorber based on a negative-index plasmonic metamaterial. Phys.
Rev. B 79(4), 045131 (2009).
29. Chen, X. et al. Robust method to retrieve the constitutive effective parameters of metamaterials. Phys. Rev. E. 70, 16608 (2004).
30. Cuong, T. M. et al. Controlled Defect Based Ultra Broadband Full-sized Metamaterial Absorber. Sci. Rep. 8, 9523 (2018).
31. Cuong, T. M. et al. Creating a Multiband Perfect Metamaterial Absorber at K Frequency Band Using Defects in the Structure. J.
Electron. Mater. 46(1), 413 (2017).
32. Landy, N. I. et al. Perfect Metamaterial Absorber. Phys. Rev. Lett. 100, 207402 (2008).
33. Landy, N. I. et al. Design, theory, and measurement of a polarization-insensitive absorber for terahertz imaging. Phys. Rev. B
2009(79), 125104 (2008).
34. Hao, J. M. et al. High performance optical absorber based on a plasmonic metamaterial. Appl.Phys. Lett. 96, 251104 (2010).
35. Cuong, T. M. et al. Creating Multiband and Broadband Metamaterial Absorber by Multiporous Square Layer Structure. Plasmonics,
https://doi.org/10.1007/s11468-019-00953-6. (2019).
36. Gong, C. et al. Broadband terahertz metamateriavl absorber based on sectional asymmetric structures. Sci. Rep. 6, 32466 (2016).
37. Kenney, M. et al. Octave-Spanning Broadband Absorption of Terahertz Light Using Metasurface Fractal-Cross. ACS Photonics 4(10),
2604 (2017).
38. Chen, K. et al. Broadband microwave metamaterial absorber made of randomly distributed metallic loops, IMWS-AMP, 7588420.
https://doi.org/10.1109/IMWS-AMP.2016.7588420. (2016).
39. Peng, Y. et al. Broadband Metamaterial Absorbers. Opt. Mater 7, 1800995 (2019).
40. Lobet, M., Lard, M., Sarrazin, M., Deparis, O. & Henrard, L. Plasmon hybridization in pyramidal metamaterials: a route towards
ultra-broadband absorption. Opt. Express 22, 12678 (2014).
41. Xin, W., Binzhen, Z., Wanjun, W., Junlin, W. & Junping, D. Fabrication and Characterization of a Flexible Dual-band Metamaterial
Absorber. IEEE Photonics J. 9, 4600213 (2017).
42. Tao, H. et al. Highly flexible wide angle of incident terahertz metamaterial absorber: Design, fabrication, and characterization. Phys.
Rev. B 78(24), 241103 (2008).
43. Mattiucci et al. Impedance matched thin metamaterials make metals absorbing. Sci. Rep. 3, 3203, https://doi.org/10.1038/srep03203
(2013).
44. Zheludev, N. I. & Eric, P. Reconfigurable nanomechanical photonic metamaterials. Nat. Nanotech. 11, 16 (2016).
45. Xiong, X. et al. Structured Metal Film as a Perfect Absorber. Adv. Mater. 25, 3994, https://doi.org/10.1002/adma.201300223. (2013).
46. Liu, C. & Qi, L. and Mingjing. Broadband infrared metamaterial absorber with polarization-independent and wide-angle
absorption. Opt. Mater. Express 8, 82439 (2018).
47. Zhang, L. et al. Space-time-coding digital metasurfaces. Nat. Commun. 9, 4334 (2018).
48. Silva, A. et al. Performing mathematical operations with metamaterials. Science 343, 160–3 (2014).
49. Cui, T. J. et al. Coding metamaterials, digital metamaterials and programmable metamaterials. Light. Sci. Appl. 3, e218, https://doi.
org/10.1038/lsa.2014.99 (2014).
50. Chen, P. Y. et al. Synthesis design of artificial magnetic metamaterials using a genetic algorithm. Opt. Express. 16(17), 12806 (2008).
51. Zhao, Y. X. et al. Prog. Electromagn. Res. Lett. 22, 95 (2011).
52. Shaw, L. A. et al. Computationally efficient design of directionally compliant metamaterials. Nat. Commun. 10, 291 (2019).
53. Gao, L. et al. Broadband diffusion of terahertz waves by multi-bit coding metasurfaces. Light Sci Appl 4, e324, https://doi.
org/10.1038/lsa.2015.97 (2015).
54. Zhang, L. et al. Realization of low scattering for a high-gain Fabry–Perot antenna using coding metasurface. IEEE Trans. Antennas
Propag. 65(7), 3374–3383 (2017).
55. Massimo, M. et al. Coding Metasurfaces for Diffuse Scattering: Scaling Laws, Bounds, and Suboptimal Design. Adv. Opt. Mater. 5,
19 (2017).
56. Jie, Z. et al. Controlling the Bandwidth of Terahertz Low‐Scattering Metasurfaces. Adv. Opt. Mater 4(11), 1773–1779 (2016).
57. Zhang, L. et al. Transmission‐Reflection‐Integrated Multifunctional Coding Metasurface for Full‐Space Controls of Electromagnetic
Waves. Adv. Funct. Mater 28, 33 (2018).
58. Zhang, L. et al. Breaking Reciprocity with Space‐Time‐Coding Digital Metasurfaces. Adv. Mate 31, 41 (2019).
59. Lei, Z., Shuo, L., Lianlin, L. & Tie, J. C. Spin-Controlled Multiple Pencil Beams and Vortex Beams with Different Polarizations
Generated by Pancharatnam-Berry Coding Metasurfaces. ACS Appl. Mater. Interfaces 9, 41 (2017).
Acknowledgements
This research was funded by the Vietnam National Foundation for Science and Technology Development (Grant
No. 103.99-2017.26).
Author contributions
M.C.T. and V.H.P. conceived the idea. M.C.T. and V.H.P. performed the simulation and carried out the theoretical
explanation. The fabrication and measurements were implemented by H.T.D., D.T.L., T.H.H.,T.T.N., X.K.B.,
S.T.B., T.L.P. The article text was reviewed and wrote by D.L.V., M.C.T. and V.H.P. All of the authors discussed and
commented on the final manuscript.
Competing interests
The authors declare no competing interests.
Additional information
Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-58774-1.
Correspondence and requests for materials should be addressed to M.C.T. or V.H.P.
Reprints and permissions information is available at www.nature.com/reprints.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and
institutional affiliations.
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 Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article’s Creative Commons license and your intended use is not per-
mitted 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 license, visit http://creativecommons.org/licenses/by/4.0/.