Achieving a multi-band metamaterial perfect absorber via a hexagonal ring dielectric resonator*
Li Li-Yang, Wang Jun†, Du Hong-Liang, Wang Jia-Fu, Qu Shao-Bo‡
College of Science, Air Force Engineering University, Xi'an 710051, China

Corresponding author. E-mail: wangjun563@163.com

Corresponding author. E-mail: qushaobo@mail.xjtu.edu.cn

*Project supported by the National Natural Science Foundation of China (Grant Nos. 61331005, 11204378, 11274389, 11304393, and 61302023), the Aviation Science Foundation of China (Grant Nos. 20132796018 and 20123196015), the Natural Science Foundation for Post-Doctoral Scientists of China (Grant Nos. 2013M532131 and 2013M532221), the Natural Science Foundation of Shaanxi Province, China (Grant No. 2013JM6005), and the Special Funds for Authors of Annual Excellent Doctoral Degree Dissertations of China (Grant No. 201242).

Abstract

A multi-band absorber composed of high-permittivity hexagonal ring dielectric resonators and a metallic ground plate is designed in the microwave band. Near-unity absorptions around 9.785 GHz, 11.525 GHz, and 12.37 GHz are observed for this metamaterial absorber. The dielectric hexagonal ring resonator is made of microwave ceramics with high permittivity and low loss. The mechanism for the near-unity absorption is investigated via the dielectric resonator theory. It is found that the absorption results from electric and magnetic resonances where enhanced electromagnetic fields are excited inside the dielectric resonator. In addition, the resonance modes of the hexagonal resonator are similar to those of standard rectangle resonators and can be used for analyzing hexagonal absorbers. Our work provides a new research method as well as a solid foundation for designing and analyzing dielectric metamaterial absorbers with complex shapes.

Keyword: 42.25.Bs; 81.05.Xj; 78.20.Bh; 78.20.Ci; multi-band; perfect absorber; all-dielectric; metamaterial
1. Introduction

In 2008, Landy et al.[1] put forward a perfect absorber based on metallic structure metamaterial, which contains an electric resonator coupling with incident electric field and a magnetic resonator coupling with incident magnetic field. Owing to strong resonance of the metallic structure, the absorptivity can reach nearly unity. Afterward, the perfect metamaterial absorber became one of the focal fields, which aroused the interest of researchers all over the world. Metamaterial absorbers based on metallic patterns developed rapidly in recent years, from polarization-independent, [2] wide angle-of-incidence absorbers[36] to multi-band, [6, 7] broadband absorbers, [8, 9] from microwave[1, 4, 6, 8] to terahertz, [2, 3, 7, 9] infrared [1013] and optical absorbers.[5, 1416]

The studies of perfect absorbers have achieved remarkable success in improving the absorbing performance. However, most of the above-mentioned metamaterial absorbers are based on metallic structures. The metallic structures have some advantages. Nevertheless, metallic structures are prone to be oxidized and corroded under high temperatures. Alternatively, dielectrics can make up these deficiencies and are advantageous in high-power and high-temperature applications, so it is meaningful to investigate dielectric absorber design. In recent years, remarkable progress has been achieved for dielectric metamaterials in theory, [1722] experiments, [17, 2325] and applications.[2628] This provides a good foundation for the study of dielectric absorbers. However, at present, few relevant studies are reported, giving no applicable examples. In 2013, Liu et al.[29] proposed a perfect absorber via high-permittivity dielectric cubes on metal plates, and the effects on absorption were studied by changing the edge length and the permittivity of the dielectric material. The absorber achieves single frequency absorption in X-band.

In this paper, a multi-band dielectric perfect absorber in microwave regime is proposed. A hexagonal ring structure made of high-permittivity dielectrics is employed to fabricate the absorber. Considering the fact that Mie scattering theory[29] and dielectric resonator theory[22] are essentially equivalent in revealing resonant modes of high permittivity dielectrics, the dielectric resonator theory (which is more convenient especially for resonators with complex shapes and boundaries) is used to analyze the multi-band perfect absorption. It is found that the resonant modes of the complex hexagonal structure are similar to those of standard rectangle resonators. This approximate methodology provides new ideas and methods for future studies on dielectric metamaterial absorbers with complex structures.

2. Design and simulation

The dielectric hexagonal ring structure is shown in Fig.  1. The structure is made of SrTiO3 (the relative permittivity is 300 and the loss tangent is 0.005). As shown schematically in Fig.  1(a), the inner side length of the hexagonal ring is 2.5  mm and the height of the hexagonal ring is 1.0  mm. The thickness values for the two upper walls and the two lower walls are all 0.433  mm, while the two vertical walls are both 0.2165-mm thick. The hexagonal ring is placed on a copper backing plate and the unit cell is 14  mm× 14  mm in area as shown in Fig.  1(b).

Fig.  1. Description of the dielectric metamaterial absorber: (a) the hexagonal ring structure with geometrical parameters; (b) unit cell and periodic structures.

The absorptivity of an absorber can be calculated from the transmission (T) and reflection (R), and is given by the following formula:

This absorber adopts a metal backing plate, so the transmission T is 0.

The commercial software CST MICROWAVE STUDIO 2011 is used for this simulation. A full-wave simulation based on the frequency domain solver is implemented in 9.0  GHz– 14.0  GHz. As shown in Fig.  2, the hexagonal ring structure has three resonant dips below – 10  dB in 9.0  GHz– 14.0  GHz. The three resonant dips are at 9.785  GHz, 11.525  GHz, 12.37  GHz, respectively. The corresponding reflections are -15.057  dB, -12.446  dB, and -25.94  dB, corresponding to absorptions about 96.88% , 94.31% , and 99.75% , respectively.

Fig.  2. Reflection (a) and absorption (b) properties for the absorber.

3. Theoretical analysis

The electromagnetic responses at each absorption peak are analyzed to reveal the mechanism of the multi-band perfect absorption. We illustrate the main distribution features of electric and magnetic fields, and point out the mechanism of the perfect absorption.

The hexagonal high-permittivity inclusions are treated as dielectric resonators. The resonance modes are analyzed based on the dielectric resonator theory. Analysis of a rectangular dielectric resonator is conducted using the method of mixed magnetic wall. This is helpful in providing the inner field distributions and the formula for calculating resonant modes under different boundary conditions.[21, 22] Take mode for example. When the surfaces of the dielectric resonators do not contact magnetic conductors (namely, equivalent magnetic walls), the normal component of the electric field vector vanishes on the magnetic wall. Therefore, electric field lines form loops inside the resonator as shown in Fig.  3(a). When one surface of the dielectric resonator in z direction contacts electric conductors, it is changed into a perfect electric wall. The tangential component of the electric field vector vanishes at this surface, and the electric field lines are vertical to this surface. Since the remaining five surfaces do not contact the conductor, the electric field lines  are tangential to these surfaces. The field distributions are shown in Fig.  3(b). Since a metal backing plate is used in the absorber, the bottom faces of the hexagonal dielectric ring structure become electric walls. Thus, the field distributions are similar to those in Fig.  3(b).

Fig.  3. Electric distributions for modes under different boundary conditions. (a) Six sides are magnetic walls. From left to right, the three panels show the distributions of TEδ 11, TEδ 31, and TEδ 51. (b) One side wall is an electric wall, and the others are magnetic walls.

We first analyze the field distribution features at the first resonant frequency 9.785  GHz. The reflection dip is -15.057  dB, with absorption about 96.88% . As shown in Fig.  4, electromagnetic fields inside the hexagonal ring are focused on the two upper walls and the two lower walls, while there is almost no field in the two vertical walls. Figure  4 shows the electric and magnetic field distributions on xy plane at 9.785  GHz. As shown in Fig.  4, the electric fields inside the upper and lower walls resonate along the walls back and forth. In order to describe the resonance process in detail, we monitor the field distributions at different excitation phases (the phase on the plane where the incident port lies). The electric field resonances are stronger when the excitation phases are 135° and 315° than when at the other phases. In addition, the magnetic fields in the walls are vertical to the electric fields. The magnetic field distribution can be regarded as a magnetic dipole. The magnetic field is more intense when the excitation phases are 45° and 225° .

Fig.  4. Electric and magnetic fields at 9.785  GHz. (a) Electric field distribution in xoy plane for the excitation phase is 135° . (b) Electric field distribution in xoy plane for the excitation phase of 315° . (c) Magnetic field distribution in xoy plane for the excitation phase of 45° . (d) Magnetic field distribution in xoy plane for the excitation phase of 225° .

Figure  5 shows the electric fields of the four upper/lower sides further simulated by Matlab programming. Denote xoy plane by plane A and the inner ring side by plane B. In Fig.  5, the bottom surface of the hexagonal ring contacts the copper plate. Thus, the bottom surface can be treated as the electric wall. The electric fields are perpendicular to the electric wall. The top surface of the hexagonal ring can be considered as a magnetic wall, on which the electric field is tangential. Therefore, semicircular electric field lines are formed on plane B. Due to the mirror effect of the copper plate, equivalent closed loops of electric field lines are formed on the four upper/lower sides. The electric field component perpendicular to the loops remains zero. The resonant mode can be seen as TE mode. Using Ploss ∝ | Eloc| 2(1 − e− 2α x) (where α represents the attenuation coefficient, x represents the thickness of the absorber), we can see that the power loss is proportional to the squared localized field intensity in the dielectric resonator. It is the resonance that leads to intensified fields and thus a strong absorption effect.

Fig.  5. Electric field simulation for the four upper/lower sides.

A simple analogy between the four upper/lower sides and standard rectangular dielectric resonators is carried out to clarify the resonant modes. As shown in Fig.  6, the standard rectangular resonator has basically the same sizes of the four upper/lower sides. The geometrical dimensions of the rectangular resonator are 2.5  mm× 0.433  mm× 1.0  mm with a tilted angle of 30° with respect to x direction. The resonant mode of the standard rectangular resonator is similar to those of the four upper/lower sides of the hexagonal ring dielectric resonator.

Fig.  6. Standard rectangular resonator extracted from the hexagonal dielectric resonator.

As shown in Fig.  7, the magnetic field and the electric field of the rectangular resonator are similar to those of the four upper/lower hexagonal ring sides at the resonant frequency. Figures  7(a) and 7(d) show the electric field distributions in plane A for the excitation phases of 135° and 315° respectively. The electric field distribution is almost the same as that of the hexagonal ring. Figures  7(b) and 7(e) show the magnetic field distributions in plane A for the excitation phases of 45° and 225° respectively. The magnetic fields are perpendicular to the electric fields, equivalent to a magnetic dipole. This resonant mode is TE mode. Figures  7(c) and 7(f) show the electric field distributions in plane B for the excitation phases of 135° and 315° respectively. It can be shown that the electric field lines are also in the shape of a semicircle, similar to those in Fig.  5.

Fig.  7. Electric and magnetic field distributions for the rectangular resonator. (a) Electric field distribution in the xoy plane for the excitation phase of 135° . (b) Magnetic field distribution in the xoy plane for the excitation phase of 45° . (c) Electric field distribution in plane B for the excitation phase of 135° . (d) Electric field distribution in the xoy plane for the excitation phase of 315° . (e) Magnetic field distribution in the xoy plane for the excitation phase of 225° . (f) Electric field distribution in plane B for the excitation phase of 315° .

The second resonant frequency is located at 11.525  GHz. The reflection dip is -12.446  dB, with absorption about 94.31% . As shown in Fig.  8, strong electromagnetic responses are excited in all the six sides of the hexagon ring. Figures  8(a) and 8(d) show electric field distributions in the xoy plane for the excitation phases of 0° and 180° respectively. Electric fields on the two vertical sides resonate in the y direction. In this case, each vertical side is equivalent to an electric dipole. Figures  8(b) and 8(e) show magnetic field distributions in the xoy plane for the excitation phases of 90° and 270° respectively. The magnetic field is perpendicular to the electric field vector. For each of the four upper/lower sides, the magnetic fields form two loops in the xoy plane, which are opposite in direction. In this magnetic resonant response, the magnetic fields vanish in the z direction (Hz = 0), while none of the electric fields in the x, y, and z directions is zero. Thus the corresponding resonance mode is TM mode. On each vertical side, the magnetic fields form one loop in the xoz plane. The magnetic field in the y direction is zero (Hy = 0) as shown in Figs.  8(c) and 8(f), indicating that the resonant mode on the vertical side is TM mode. Figure  8(c) shows that the magnetic field forms an asymmetric ring, leading to a composite magnetic field in the positive x direction for the excitation phase of 90° . Similarly, the composite magnetic field direction is along the negative x direction for the excitation phase of 270° as shown in Fig.  8(f). From Ploss ∝ | Eloc| 2(1 − e− 2α x), we can see that the power loss is proportional to the squared field intensity. It is also the TM mode resonance on the six sides that leads to the enhanced localization field and produces a perfect absorption at the resonant frequency.

Fig.  8. Electric and magnetic field distributions at 11.525 GHZ. (a) Electric field distribution in the xoy plane for the excitation phase of 0° . (b) Magnetic field distribution in the xoy plane for the excitation phase of 90° . (c) Magnetic field distribution in the xoz plane for the excitation phase of 90° . (d) Electric field distribution in the xoy plane for the excitation phase of 180° . (e) Magnetic field distribution in the xoy plane for the excitation phase of 270° . (f) Magnetic field distribution in the xoz plane for the excitation phase of 270° .

The third resonant frequency is located at 12.37  GHz, with a dip of -25.94  dB. The absorption is about 99.75% . As shown in Fig.  9, a strong electromagnetic response happens on all six sides of the hexagonal ring resonator, and the response is more intense on the two vertical sides. Figures  9(a) and 9(d) show the electric field distributions in the xoy plane when the excitation phases are 0° and 180° , when the electric fields reach their highest intensities. Each of the electric fields on the vertical sides resonates in the y direction, leading to an electric dipole. Figures  9(b) and 9(e) show the magnetic field distributions in the xoy plane when the excitation phases are 90° and 270° respectively. The magnetic field vector is perpendicular to the electric field vector. For each of the four upper/lower sides, the magnetic fields form three loops in the xoy plane, which are opposite in direction alternatively one by one. In this case, the magnetic fields are zero in the z direction (Hz = 0), while none of the electric fields in the x, y, and z directions is zero. Thus the corresponding resonance mode is TM mode. For the vertical sides, the magnetic fields form one loop in the xoz plane. The magnetic fields in the y direction are all zero (Hy = 0) as shown in Figs.  9(c) and 9(f), indicating that the resonant modes on the vertical sides are all TM mode. Figure  9(c) shows that the magnetic fields form an asymmetric ring, leading to a composite magnetic field in the positive x direction when the excitation phase is 90° . Similarly, the composite magnetic field direction is along the negative x direction when the excitation phase is 270° as shown in Fig.  9(f). The power loss is proportional to squared field intensity. It is the TM resonant modes on the four upper/lower sides and the two vertical sides that lead to enhanced field intensity. This produces a perfect absorption at the resonant frequency.

Fig.  9. Electric and magnetic field distributions at 12.37  GHZ. (a) Electric field distribution in the xoy plane for the excitation phase of 0° . (b) Magnetic field distribution in the xoy plane for the excitation phase of 90° . (c) Magnetic field distribution in the xoz plane for the excitation phase of 90° . (d) Electric field distribution in the xoy plane for the excitation phase of 180° . (e) Magnetic field distribution in the xoy plane for the excitation phase of 270° . (f) Magnetic field distribution in the xoz plane for the excitation phase of 270° .

The multi-band perfect absorption is investigated by analyzing the electric and magnetic field distributions at each absorption peak frequency. It is shown that the multi-band perfect absorptions originate from the enhanced field intensity. The enhanced localized fields, which are excited in the high dielectric resonators, lead to the enhanced absorption loss for incident waves.

We investigate the variation of absorption peaks with the whole size of the dielectric resonator, the metal backing plate, and the dielectric permittivity. We set the scale of the original absorber as S = 1, then the scale factor varies from 0.7 to 1.2 (under ɛ = 300). Then, we change the dielectric permittivity from 285 to 315 (under S = 1). It can be seen from Fig. 10 that if we change the size or permittivity of the dielectric resonator, the absorption peak of the absorber will change, which is in good agreement with theoretical predictions. Figure  10(a) shows that the simulated absorption peak frequency decreases as the whole size of the absorber increases from 0.7 to 1.2 (under ɛ = 300). Figure  10(b) indicates that the absorption peak frequency decreases as the ɛ increases from 285 to 315 (under S = 1).

Fig.  10. Absorption properties for absorbers with different geometrical scales (a) and permittivies (b).

4. Conclusions

In this paper, a multi-band absorber in the microwave regime is proposed using a hexagonal ring dielectric resonator. Unlike the metallic pattern-based absorbers, this absorber is composed of dielectric resonators and a metallic ground plate. Such an absorber is superior in inoxidizability, anticorrosion, and high-temperature resistance and high-power absorption, compared with their metallic counterpart. The mechanism for the near-unity absorption is caused by the electric and magnetic resonances where enhanced localized fields are excited in the dielectric resonator. In addition, the research offers an analogical method of simplifying a complex geometrical dielectric resonator into an equivalent simple dielectric resonator with similar resonance modes. By analyzing its resonant modes and adjusting resonant frequencies, broadband dielectric perfect absorbers can be achieved.

Reference
1 Land y N I, Sajuyigbe S, Mock J J, Smith D R and Padilla W J 2008 Phys. Rev. Lett. 100 207402 DOI:10.1103/PhysRevLett.100.207402 [Cited within:2] [JCR: 7.943]
2 Land y N I, Bingham C M, Tyler T, Jokerst N, Smith D R and Padilla W J 2009 Phys. Rev. B 79 125104 DOI:10.1103/PhysRevB.79.125104 [Cited within:2]
3 Tao H, Bingham C M, Strikwerda A C, Pilon D, Shrekenhamer D, Land y N I, Fan K, Zhang X, Padilla W J and Averitt R D 2008 Phys. Rev. B 78 241103 DOI:10.1103/PhysRevB.78.241103 [Cited within:2]
4 Cheng Y Z, Yang H L, Cheng Z Z and Wu N 2011 Appl. Phys. A 102 99 DOI:10.1007/s00339-010-6022-4 [Cited within:1] [JCR: 1.545]
5 Hao J M, Wang J, Liu X L, Padilla W J, Zhou L and Qiu M 2010 Appl. Phys. Lett. 96 251104 DOI:10.1063/1.3442904 [Cited within:1] [JCR: 3.794]
6 Li L, Yang Y and Liang C 2011 J. Appl. Phys. 110 063702 DOI:10.1063/1.3638118 [Cited within:3] [JCR: 0.71]
7 Tao H, Bingham C M, Pilon D, Fan K, Sorikwerda A C, Shrekenhamer D, Padilla W J, Zhang X and Averitt R D 2010 J. Phys. D: Appl. Phys. 43 225102 DOI:10.1088/0022-3727/43/22/225102 [Cited within:2] [JCR: 2.528]
8 Bao S, Luo C R, Zhang Y P and Zhao X P 2010 Acta Phys. Sin. 59 3187(in Chinese) 118.145.16.217/magsci/article/article?id=173916542010 [Cited within:2] [JCR: 1.016] [CJCR: 1.691]
9 Grant J, Ma Y, Saha S, Khalid A and Cumming David R S 2011 Opt. Lett. 36 3476 DOI:10.1364/OL.36.003476 [Cited within:2] [JCR: 3.385]
10 Liu X, Starr T, Starr A F and Padilla W J 2010 Phys. Rev. Lett. 104 207403 DOI:10.1103/PhysRevLett.104.207403 [Cited within:1] [JCR: 7.943]
11 Liu X, Tyler T, Starr T, Star A F, Jokerst N M and Padilla W J 2011 Phys. Rev. Lett. 107 045901 DOI:10.1103/PhysRevLett.107.045901 [Cited within:1] [JCR: 7.943]
12 Mason J A, Smith S and Wasserman D 2011 Appl. Phys. Lett. 98 241105 DOI:10.1063/1.3600779 [Cited within:1] [JCR: 3.794]
13 Cui Y, Xu J, Fung K H, Jin Y, Kumar A, He S and Fang N X 2011 Appl. Phys. Lett. 99 253101 DOI:10.1063/1.3672002 [Cited within:1] [JCR: 3.794]
14 Aydin K, Ferry V E, Briggs R N and Atwater A 2011 Nat. Commun. 2 517 DOI:10.1038/ncomms1528 [Cited within:1] [JCR: 10.015]
15 Wu C, III B N and Shvets G 2011 Phys. Rev. B 84 075102 DOI:10.1103/PhysRevB.84.075102 [Cited within:1]
16 Wu C and Shvets G 2012 Opt. Lett. 37 308 DOI:10.1364/OL.37.000308 [Cited within:1] [JCR: 3.385]
17 Peng L, Ran L, Chen H, Zhang H, Kong J A and Grzegorczyk T M 2007 Phys. Rev. Lett. 98 157403 DOI:10.1103/PhysRevLett.98.157403 [Cited within:2] [JCR: 7.943]
18 Lai Y J, Chen C K and Yen T J 2009 Opt. Express 17 12960 DOI:10.1364/OE.17.012960 [Cited within:1] [JCR: 3.546]
19 Kim J, Yoo H S and Gopinath A 2010 Electron. Lett. 46 1590 DOI:10.1049/el.2010.8795 [Cited within:1]
20 Lepetit T, Akmansoy É and Ganne J 2011 J. Appl. Phys. 109 023115 DOI:10.1063/1.3537920 [Cited within:1] [JCR: 0.71]
21 Wang J, Xu Z, Yu Z H, Wei X Y, Yang Y M, Wang J F and Qu S B 2011 J. Appl. Phys. 109 084918 DOI:10.1063/1.3575326 [Cited within:1] [JCR: 0.71]
22 Du B, Wang J, Xu Z, Xia S, Wang J F and Qu S B 2014 J. Appl. Phys. 115 234104 DOI:10.1063/1.4883962 [Cited within:3] [JCR: 0.71]
23 Zhang W, Potts A, Bagnall D M and Davidson B R 2006 J. Vac. Sci. Technol. B 24 1455 [Cited within:1] [JCR: 1.267]
24 Zhao Q, Kang L, Du B, Zhao H, Xie Q, Huang X, Li B, Zhou J and Li L 2008 Phys. Rev. Lett. 101 027402 DOI:10.1103/PhysRevLett.101.027402 [Cited within:1] [JCR: 7.943]
25 Schuller J A, Zia R, Taubner T and Brongersma M L 2007 Phys. Rev. Lett. 99 10 [Cited within:1] [JCR: 7.943]
26 Ueda T, Michishita N, Akiyama M and Itoh T 2008 IEEE Trans. Microwawe Theor. Technol. 56 2259 DOI:10.1109/TMTT.2008.2003524 [Cited within:1]
27 Kim J and Gopinath A 2009 IEEE Anten. Technol. 2009 1 DOI:10.1109/IWAT.2009.4906878 [Cited within:1]
28 Gaillot D P, Croenne C and Lippens D 2008 Opt. Express 16 3986 DOI:10.1364/OE.16.003986 [Cited within:1] [JCR: 3.546]
29 Liu X, Zhao Q, Lan C and Zhou J 2013 Appl. Phys. Lett. 103 031910 DOI:10.1063/1.4813914 [Cited within:2] [JCR: 3.794]