1 Introduction
Metamaterials are artificial materials [
1,
2] that have a variety of special physics phenomena, including negative refraction [
3], reversals of both Doppler shift and Cherenkov radiation [
4], perfect absorption [
5], and super lenses [
6]. These phenomena cannot be found in natural materials. Owing to these properties, metamaterials have an amount of applications in civilian and military products such as electromagnetic compatibility (EMC) [
7], frequency selective surface (FSS) [
8], sensors [
9], transformation optics [
10,
11], and invisibility cloaks [
12]. The electromagnetic response in metamaterials is referred to as the electromagnetic multipole. The electric and magnetic multipoles play significant roles in the electromagnetic multipole family. However, as a third family, the toroidal multipole [
13] is negligible because of the other multipoles’ strong excitation. The toroidal dipole was first proposed by Zel’dovich [
14]. The formation of a toroidal dipole can be introduced by a current flow on a toroidal surface along its meridian. A distinct feature of toroidal metamaterials, anapole excitation, is of nontrivial, nonradiating nature [
15]. Recently, research surrounding toroidal dipoles in metamaterials has become popular and has shown many unusual properties, including polarization conversion [
16], unconventional optical activity [
17] and resonance transmission [
18].
Electromagnetically induced transparency (EIT) [
19,
20] is a three-level quantum interference process in an atomic system. When the different energy level coupling, the transmission spectrum forms a sharp transparency peak at a specific frequency. However, achieving the EIT phenomena needs a series of strict environmental conditions (a low-temperature, high-intensity optical pump) in experimental circumstances. Metamaterials provide convenience for researchers to realize EIT without difficult setups. The EIT phenomena has supported by considerable research from microwave, optics, near-infrared [
21], etc. The EIT metamaterial has many potential applications, for example, the sensing of dielectric media and biomolecules [
22], refractive index sensors [
23] and optical switches [
24]. The bright-dark mode coupling [
25,
26] is the most basic method to acquire a classical EIT phenomenon. Generally, the bright mode can be directly excited by an incident electromagnetic wave and has a low quality factor (
Q-factor). Therefore, the bright mode can be seen as the electric dipole response. By contrast, the dark mode is weakly radiated by an electromagnetic wave and has a high
Q-factor, which serves as the magnetic dipole response.
In recent years, many researchers have paid more attention to EIT and Fano resonance [
27] in metamaterials with the toroidal dipole response. In 2013, Fedotov
et al. [
28] proposed four dumbbell-shaped apertures in a thin metal plate, constituting a compass-shaped cross-section element, which created destructive interference between the toroidal and electric dipoles and acquired symmetric Lorentzian transparency lines. Guo
et al. proposed an electric toroidal dipole that etched a sun-like aperture metal pattern at Rogers 4003C in 2014 [
29]. In 2015, Han
et al. [
30] observed Fano resonance and EIT in a three-dimensional toroidal metamaterial that consisted of three identical combine split resonance ring by sharing a central connecting bridge, forming a hexagon. The first demonstration of toroidal excitation in 2D metamaterials was shown by Gupta
et al. [
31] in 2016. The setup consisted of two joint metallic loops with two capacitive gaps in each loop. Zhu
et al. [
32] employed a type-I cut wire and two spiral ring resonators to achieve a planar toroidal metamaterial in 2017. This work indicated that EIT could be attributed to the coupling between the electric and toroidal dipoles. In 2018, the classical Ising model was used to verify Fano resonance with strong magnetic and toroidal dipole responses, which was proposed by Cong
et al. [
33].
A planar dual-band EIT E– metamaterial is proposed, named so as it consists of E- and -shaped structures, with transparency peaks located at 7.60 and 10.27 GHz. The simulation and experimental results are seen to be consistent. The calculation results of the radiated power multipoles show that the toroidal dipole response is weaker than the electric dipole response at 7.60 GHz and stronger than the electric dipole response at 10.27 GHz. This phenomenon could be explained by the surface current distribution and the magnetic field. With increasing incident and polarizing angles, the transmission spectra change, and the Q-factors of the two transparency peaks increase. The proposed E– metamaterial can be used in refractive index sensors.
2 Design, simulation, and experiment
The unit cell structure is made up of two parts: E-shaped and -shaped metal patterns. The CST Microwave Studio simulation software is used to obtain simulation results, in which the x- and y-axes are set as unit cell boundaries and the z-direction is an open space (add). The geometrical parameters of the two different unit cell structures are shown in Fig. 1(a). The dielectric slab is 1 mm-thick F4B with permittivity and loss tangents of 2.2 and 0.001. The metal pattern is formed of 0.035 mm-thick copper with a conductivity of S/m. Fig. 1(b) shows the transmission spectra of the E-shaped pattern, in which the dip of the transmission coefficient was dB at 8.07 GHz. The transmission spectra of the -shaped pattern are shown in Fig. 1(c), where the frequency of the transmission dip was 10.34 GHz with a dB transmission coefficient. The Q-factors were 3.34 for the E-shaped pattern at 8.07 GHz and 6.01 for the -shaped pattern at 10.34 GHz, which were solved by using the formula , where f is the resonance frequency and FWHM is the full width at half maximum bandwidth. The transmission spectra of the E-shape show a line at 10.34 GHz, and the transmission coefficient approximates to 0 at 8.07 GHz for the -shape. Therefore, the E-shaped structure shows bright-mode coupling for the low frequency and dark-mode coupling for the high frequency.
In order to achieve EIT, the unit cell comprises two mirrored and interlocked E- or -shaped patterns. However, the transmission spectra only appear one transparency peak in two different metamaterials, as shown in Fig. 2. The transparency peak was located at 8.59 GHz, with a transmission coefficient of dB, as shown in Fig. 2(a), for the mirrored E-shaped pattern. The transmission spectrum of the mirrored asymmetric -shapes is shown in Fig. 2(b), with the transparency peak appearing at 13.35 GHz with a transmission coefficient of dB.
Using the above simulation results, the dual-band EIT metamaterial can be realized by combining the E-shape and the -shape in a single unit cell, as shown in Fig. 3(a). The geometrical parameters of this setup are the same as those in Fig. 1. The 200 mm × 200 mm sample was fabricated using a printed circuit board, and its image is shown in Fig. 3(b). The experimental results are measured using a vector network analyzer (Agilent PNA E8362B) and a pair of horn antennas in free space.
Figure 4 shows the transmission spectra of the E– metamaterial from the simulation and experiment results. In the simulation, there were two transparency peaks located at 7.60 and 10.27 GHz, with corresponding transmission coefficients of and dB, as well as three transmission dips at 7.22, 9.66, and 11.83 GHz, with corresponding transmission coefficients of , , and dB. In the experiment, the resonance frequencies of two transparency peaks were 7.67 and 10.25 GHz, for which the transmission coefficients were and dB. This agrees well with the simulation results. However, the three transmission dips were found to be at 7.30, 9.73, and 11.76 GHz, which deviates from the values found via simulation. Their transmission coefficients were , , and dB, respectively. There are two main factors causing the disparity between the simulation and experiment. Firstly, the permittivity or loss tangle of the dielectric slab in the experiment has a little discrepancy compared to that in the simulation, and second is the presence of strong electromagnetic wave diffraction at a higher frequency.
3 Results and discussion
To further probe the characteristics of the EIT E– metamaterial, we verify a variety of transmission spectra with different incident and polarizing angles by numerical simulation and experimental measurement. Fig. 5(a) shows the simulation results for the different incident angles , which changed from 0° to 60° with a step of 15°. When the incident angle rose from 0° to 45°, the four curves of transmission spectra became almost identical. The width of the two transmission windows decreases at an incident angle of 60°. However, the resonance frequencies of the two transparency peaks remain at 7.60 and 10.27 GHz. The Q-factors of the two transparency peaks are presented in Table 1 and are shown to increase with the incident angle. The experimental results are shown in Fig. 5(b). The resonance frequencies of the two transparency peaks perfectly coincide with the simulation results, which are 7.67 and 10.25 GHz. When the incident angles change from 0° to 60°, the variation of the transmission windows is nearly identical, and the width of the transmission windows decreases.
The simulated and measured results with different polarizing angles are presented in Fig. 6. The polarizing angles ranged from 0° to 45° with a step of 15°. The amplitude of the transmission dips gradually increased, and the resonance frequencies became blue-shifted, as shown in Fig. 6(a). For instance, in the second transmission dip, the transmission coefficient was dB and the resonance frequency was 9.66 GHz when φ=0°. When φ=45°, the resonance frequency moved to 9.84 GHz, and the corresponding transmission coefficient was dB. The first and third transmission dips have the same variation tendency and resonance frequencies, and the transmission coefficients of two transparency peaks remain unchanged. In Fig. 6(b), the measured results are in agreement with the simulated results. The resonance frequencies of the two transparency peaks do not vary far from around 7.67 and 10.25 GHz. The amplitudes of the three dips are higher than those for the simulation results at different polarizing angles. However, the transmission spectra show an identical frequency variation trend to that in the simulation.
To elaborate on the coupling process of the proposed E– metamaterial, the surface current distribution and the magnetic field are depicted in Figs. 7 and 8, respectively. We employed the same amplitude scale separately that make all results accomplishing a correct comparison in Figs. 7 and 8. The surface current distributions chosen correspond to the two transparency peaks and three transmission dips: 7.60 and 10.27 GHz, and 7.27, 9.66 and 11.85 GHz, respectively. At 7.60 GHz, the surface current mainly concentrates on the ring around -shaped and E-shaped patterns, taking the form of a magnetic dipole. Meanwhile, the surface current in two transverse lines flows in the opposite direction. Figure 7(b) shows the distribution of the surface current at 10.27 GHz. In the -shaped and E-shaped patterns, it forms two looped currents in opposite directions, which is equivalent to two magnetic dipoles and interconnected head-to-tail circle. In this way, a toroidal dipole is formed. In Fig. 7(c), the surface current forms a ring, like Fig. 7(a), and the current intensity at the transparency dip is lower than that at the transparency peak. There is an almost inexistent current at 9.66 GHz, which is shown in Fig. 7(d). The surface current distribution of 11.85 GHz is depicted in Fig. 7(e), shown to concentrate on the -shaped part of the pattern, and is equivalent to an electric dipole.
The magnetic fields of the two transparency peaks are shown in Fig. 8. In Fig. 8(a), the magnetic field is flowing from the back of the dielectric substrate to the front metal, which can be viewed as a magnetic dipole. At 10.27 GHz, most of the magnetic field is restricted to flowing inside the dielectric substrate and forms a closed circular magnetic current, which is shown in Fig. 8(b). It is considered as a toroidal dipole.
The radiated power of multipoles is an important method in the quantitative analysis of the EIT-like metamaterial coupling mechanism. The calculation formulas of the radiated power of multipoles are based on the conducting current density j and are given as
where
j is the induced current density,
c is the speed of light in vacuum,
ω is the frequency and
r is the distance vector from the origin to point
in a Cartesian coordinate system
. Therefore, the decomposed far-field scatter power by these multipole moments can be calculated using the following equations [
34]:
Eqs. (5)–(8) give the radiated power of the electric, magnetic, and toroidal dipoles and magnetic quadrupole, respectively. The calculation results for the radiated power of the multipoles are exhibited in Fig. 9. We can observe that the radiated powers of the electric and magnetic dipoles decline at 7.60 and 10.27 GHz. When it comes to the magnetic quadrupole, it significantly increases in the two transparency peaks. However, the change of toroidal dipole is very clear: firstly, there are two deep dips at 7.27 and 11.85 GHz, and additionally, two strong toroidal dipole responses appear in the two transparency peaks. The toroidal dipole is lower at 7.60 GHz and is larger at 10.27 GHz, as shown in Fig. 9. At 9.66 GHz, the radiated power of the toroidal dipole drops to the lowest point between 7.60 and 10.27 GHz. This manifests that coupling between electric and toroidal dipole can generate an EIT-like material.
Refractive index sensors are a common application for EIT-like metamaterials and have potential applications in environment sensing. Sensitivity (S) is a basic parameter to indicate the performance of a sensor. It is calculated by dividing the shift frequency by the refractive index unit (RIU) [
35]. The simulation results of various refractive indexes are presented in Fig. 10. Here, the refractive index ranges from 1.0 to 1.5 with a step of 0.1. With an increasing refractive index, the resonance frequencies of the two transparency peaks are both red-shifted. Ultimately, formulas are utilized to determine the refractive index that affects the transparency peak’s resonance frequency and are based on the resonance frequency (the first transparency peak
f1 and the second transparency peak
f2) and refractive index
n, as shown below:
The resonance frequencies of different refractive indexes and transparency peaks of the simulation results can be observed in Fig. 11. At the first transparency peak, the simulation points are in agreement with the fitting curve in Fig. 11(a), in which S is 1.57 GHz/RIU. The simulation points of Fig. 11(b) show a little deviation from the fitting curve. The S of the second transparency peak is 2.10 GHz/RIU.
4 Conclusions
In conclusion, an E– metamaterial with dual-band EIT is proposed, where the destructive interference between the electric and toroidal dipoles can generate EIT. The distribution of the surface current and the magnetic field verify the electric and toroidal dipoles. The radiated power of the multipoles reveals that the toroidal dipole response is weaker than the electric dipole response in the first transparency peak and stronger than the electric dipole response in the second transparency peak. The experimental results are in agreement with the simulation results. The proposed E– metamaterial can be used to achieve refractive index sensing of the surrounding medium, which will have potential applications in various civilian and military products.
Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature 2020