We report a new design of triple-band electric metamaterial resonator (TBEMR) based on integration of open delta loops within square ring resonator. This metamaterial resonator has three distinct -negative regions (ENG) over C, X frequency bands. The transmission and reflection response of the proposed subwavelength resonator is analyzed using full-wave electromagnetic solver Ansys HFSS to demonstrate the presence of electrical resonances within frequencies 4–12 GHz. Effective medium parameters permittivity and permeability are extracted from simulated complex scattering parameters to verify existence of ENG regions. The investigations are also carried out regarding dependence of resonant frequencies on design parameters of the TBEMR unit cell. The effective medium ratio for proposed subwavelength resonator is compared with various other metamaterial resonators to indicate its compact nature.

1. Introduction

Metamaterials are artificially structured composites made of subwavelength resonators that exhibit unique electromagnetic properties not found in individual constituent components. The classification of metamaterials is carried as -negative (ENG), -negative (MNG), or double-negative (DNG) based on negative values for either permittivity, permeability, or both [1]. These single negative metamaterials including ENG and MNG have exotic electromagnetic properties different from most naturally occurring materials. The first artificial realization of metamaterial was made of two different types of subwavelength structures, thin wire and split ring resonator (SRR), to get negative response for permittivity and permeability over common frequency range [2]. Metamaterials with negative refractive index can be realized by combining the single negative ENG and MNG structures within same unit cell [3].

These single negative metamaterials are important constituent in the formation of negative index LH metamaterials. Additionally, these single negative metamaterials also find applications in gain enhancement for antennas [4], microwave absorber [5], material characterization [6], and microwave sensors [7, 8]. Applications in microwave devises such as gain enhancement or beam forming for antennas for multiband operations require multiband characteristics for metamaterials [9]. To solve these problems several single negative (SNG) metamaterial structures having multiband responses were reported in the literature [1013]. A butterfly pattern resonator reported in [10] has triple-band response with two magnetic and one electric resonance. A four-arm square spiral metamaterial reported in [11] has dual-band magnetic and electric response in different frequency band. A square single loop resonator (square-SLR) [14] is reported to have triple-band response in X-band. This square-SLR has two distinct MNG regions and one very weak ENG region. A hybrid-SLR is reported to have quad-band response in S,C bands with two MNG and two ENG responses [12]. However, the above reported metamaterials have mixed response with MNG and ENG regions in nonoverlapping frequency bands. Conversely, the nested U-shaped ring resonators reported in [13] have triple-band response with only MNG characteristics in X-band region. A summary on conventional ENG metamaterials and recent developments in it are reviewed in the following paragraphs.

Conventional, single negative ENG metamaterials like thin wire structure were shown to have dielectric function to be of the form like plasma [15]. A planar complementary SRR reported in [16] shows ENG response and acts as an alternative to thin wire structure. However, thin wire and CSRR structure due to their interplane connectivity impose difficulty in fabricating countered objects like lenses. As a shortfall of this, Schurig et al. [17] proposed electrically coupled LC resonator for obtaining ENG response. Prime advantages of E-LC resonator are planar design, pure-electric response due to mirror symmetry, no need of intercell and interplane connectivity, and improved ratio. The E-LC resonator was miniaturized further by replacing capacitive gap with interdigitated (IDC) capacitance [18]. The maximum value for effective medium ratio achieved was 10.2. It has advantage of improved homogeneity and electrical size reduction without any fabrication limitations. However, these advantages were achieved at the cost of reduced coupling to the external field as indicated by low value (−15 dB) of transmission minima. Kolb et al. [19] proposed other notable low resonant frequency structures derived from E-LC resonator for further improvement in ratio. These structures were inspired from original E-LC resonator. It includes geometries such as I-shaped resonator, thin square loop, meander-shaped resonator, and its variations. A meander-shaped resonator presents electrical size reduction with lowest resonant frequency. However, it occurs at the cost of cross-polarization effect and reduced coupling to the external field. A planar Z-shaped resonator [20] for negative permittivity was shown to have lower resonant frequency than conventional E-LC resonator with same physical size. The double S-shaped resonators realized in [21, 22] give overlapped response for negative permittivity and permeability. However, above reported metamaterial resonators were having single-band response for negative permittivity property.

In recent years, research in this particular area includes electrical size reduction, improvement in effective medium ratio [23], and achieving dual-band characteristics [2426]. Size reduction is achieved by a common idea of raising overall inductance and capacitance in the unit cell, thus achieving lower resonant frequency. The approaches include changing geometrical feature dimensions in the structure, that is, shortening the gaps or lengthening the metallic strips [23], and use of surface mounted capacitor [27] and high permittivity substrate as host material [25]. However, above approaches have limitations due to difficulty and complexity in the fabrication and cost and availability of high permittivity substrates. The dual-band characteristic is obtained by combining different resonating structures and optimizing them to get response in desired frequency band [2426]. A pair of E-LCs resonators with different sized inductive loops and capacitive gaps were shown to have two distinct electric resonances [24, 25]. However, shortfalls in these designs are poor coupling to external field at higher resonant frequency [24] and use of high permittivity material as substrate [25]. A dual-band negative permittivity metamaterial design based on integration of the delta loop within closed square ring resonator was proposed in [26]. However, meander loaded delta increases complexity of the design.

Metamaterial design proposed in this paper has purely epsilon-negative response with triple-band characteristics. The multiband responses for this structure are obtained by integrating open delta loops with square ring resonator. Moreover effect of design parameters such as arm length of delta loops, position of shorting stubs, gaps formed due to open delta loops, metallic strip width, and substrate properties are investigated in detail. Also, proposed resonator shows notable improvements over other multiband structures in terms of pure-electric response, polarization stable response, use of low cost FR-4 as host material, and design simplicity for ease in the fabrication at higher frequencies.

2. Design and Simulations of Triple-Band Electric Metamaterial Resonator (TBEMR)

In this communication, triple-band electric metamaterial resonator (TBEMR) is proposed and investigated which couples strongly to electric component of the incident field. Design of this metamaterial resonator is based on integrating open delta loops within square ring resonator. The geometry of proposed resonator is described in Figure 1 together with design parameters and incident electromagnetic wave excitation. We started with square loop resonator that already appeared in the literature [19]. The square loop resonator that couples strongly to electric component of electromagnetic wave is widely used for getting -negative response for normal wave incidence as shown in [19, 29]. Although square loop resonator has advantage in terms of polarization insensitive behaviour for normally incident wave [29], the disadvantage is regarding single-band response and electrically large size as indicated by very high value of . The open delta loop resonator (ODLR) proposed by Sabah [30] is shown to have single-band magnetic response for parallelly incident wave, whereas Sarkhel [26] has recently shown electric response for meander-line loaded closed delta loop resonator (CDLR) for normally incident wave. However disadvantage of both designs is that they are highly polarization sensitive. In case of our design reported as shown in Figure 1(a), we have used open delta loop resonator (ODLR) formed by metallic strips , , , and . Two such ODLRs are used in the design of TBEMR. The open vertex of the each ODLR is connected with the diagonally opposite end of square ring resonator by shorting stubs and . The gaps formed between the base portions of the adjacent ODLR are responsible for creating capacitances. Due to diagonal symmetry of the structure, we get exactly same triple-band response for vertically (-polarized) and horizontally (-polarized) polarized electromagnetic wave when structure is rotated by angle 90° around -axis.

The proposed TBEMR is to be printed on single face copper cladded FR-4 epoxy dielectric spacer having permittivity , loss tangent () of 0.02, and a thickness of 1.6 mm. The unit cell of the proposed resonator has physical dimensions of 6 mm × 6 mm. All metallic strips in the design of TBEMR are made of copper having conductivity of  S/m and thickness of 35 m. The copper line width considered in simulation is 0.2 mm for fabrication tolerances. The electromagnetic response from the proposed resonator is investigated using finite element method based Ansys HFSS solver version 14. Simulations are performed on a unit cell of proposed resonator by applying perfect electric (PEC) and perfect magnetic (PMC) boundary conditions. The PEC and PMC boundaries along with two wave ports model source of plane wave excitation and indicate that unit cells are repeated infinitely along and direction. Perfect electric and perfect magnetic boundaries are assigned to the planes parallel to and . This indicates that electric field is oriented along -axis and magnetic field is oriented along -axis as depicted in Figure 1(b). We obtained complex scattering parameters and from unit cell simulation in the frequency range of 3–12 GHz. The complex S-matrix data is deembedded appropriately at the edges of the cubic unit cell. At a higher value of frequency sweep, that is, 12 GHz, the ratio has value 4.16 which is in accordance with the limit of 4. This ensures that unit cell dimensions are much smaller than operating wavelength and it conforms to the effective medium theory without introducing undesired effects of spatial dispersion [31].

Numerically simulated magnitude and phase of transmission () and reflection () spectra against frequency for proposed TBEMR is presented as shown in Figures 2(a) and 2(b). The dip in transmission spectra indicates stopband response and is an indication of resonance for single negative metamaterial structure. The transmission minima () is displayed at resonant frequencies 4.32 GHz, 7.55 GHz, and 9.76 GHz as shown in Table 1. This fact is also confirmed from the pulse like transition seen in transmission phase spectra at respective resonant frequencies as shown in Figure 2(b). The absorbance spectra for proposed resonator can be calculated from expression as shown in Figure 3. It peaks at frequencies 4.65 GHz and 8.38 GHz, where reflection spectra has minimum value. The absorptivity has maximum value of 0.36 and 0.31 for proposed resonator and it can be minimized to some extent by using high conductivity metal and high quality substrate. Although loss characteristic is drawback for many applications but it can be turned into an advantage for application like wave absorber [32]. The ideal value for loss characteristics () is near unity. It can be achieved by matching surface impedance of metamaterial slab with wave impedance and minimizing wave transmission in manner as shown in [32].

3. Parameter Extractions, Surface Current, and Electric Field Distribution

The effective medium parameters are extracted from the simulated complex reflection and transmission coefficients of wave incident on the metamaterial unit cell. A computer program is written to retrieve effective impedance, refractive index, effective permittivity, and permeability based on the formulation that appeared in [33]. For normally incident plane wave on the slab of a metamaterial unit cell, constitutive parameters such as effective impedance and refractive index can be computed from Fresnel’s relation by inverting S-matrix data.

The sign of the wave impedance in (1) is determined by imposing the conditions of passive material, that is, real part of the effective impedance and imaginary part of refractive index , to be greater than zero.

The complex refractive index can be determined fromwhere is an integer denoting the branch index of complex logarithm, is the thickness of the unit cell of metamaterial sample, is the wavenumber, and depends on wave impedance . The imaginary part of in (3) is not affected by the branches of logarithm function and it can be determined accurately. However, real part has multiple solutions due to branch ambiguity. It can be determined by the method based on the Kramers-Kronig relations to ensure the uniqueness of the effective parameters [34]. The effective permittivity and permeability can be determined from the effective impedance and refractive index as shown below

The effective medium parameters permittivity, permeability, and refractive index are obtained using the above formulation for the proposed TBEMR metamaterial resonators as shown in Figures 4(a), 4(b), and 4(c), respectively. For triple-band metamaterial (TBEMR), the electric resonance is clearly seen from the real part of permittivity which varies from high positive to negative values near narrow range around resonant frequency as shown in Figure 4(a). The real part of permittivity has negative values over frequencies 4.24–4.5 GHz, 6.6–8.28 GHz, and 8.95–11.28 GHz for first, second, and third band respectively. The proposed TBEMR shows only electric response and no magnetic response as observed from Figures 4(a) and 4(b). The plot of effective permeability has near unity value. However, antiresonant response is observed in corresponding negative permittivity frequency band. The antiresonant behaviour is an intrinsic property of such resonators due to finite spatial periodicity [35]. The antiresonant response is also seen from the imaginary part of the effective permeability which is negative at resonant frequency showing unphysical behaviour.

In quasistatic limit, when size of the resonator is very small compared to operating wavelength, it can be approximated with the LC resonant circuit. The resonant frequency is inversely proportional to square root of inductance and capacitance (). The current induced in the metallic strips offers inductance, whereas gaps within the unit cell or edges between adjacent unit cells form capacitance. Thus, inductance and capacitance are influenced mainly by geometrical design parameters of the unit cell. The surface current density and electric field distributions can be plotted at respective resonant frequencies to understand inductive and capacitive response from the structure (refer to Figures 5 and 6, resp.). At lower resonant frequency, that is, GHz, it is observed that the surface current flows over top right and bottom left portion of the TBEMR resonator (refer to Figure 5(a)). The surface current flows in clockwise and anticlockwise manner in arrow-like portion of the resonator. This, clockwise and anticlockwise flow of the time varying surface current in the upper right and lower left metallic portion of the resonator cancels any magnetic dipolar activity along -axis near the resonance. The magnetic field from the external plane wave excitation is oriented along -axis and it lies in the plane - of substrate without threading the metallic loops of the TBEMR resonator. This fact explains that resonator does not couple to magnetic field component of the plane wave electromagnetic excitation. The strong intensity of the localized electric field observed near upper right and lower left inclined strips of the unit cell as shown in Figure 6(a) indicates capacitive effects. The two delta loops are coupled by mean of capacitance strongly along inclined strips . Capacitive effect occurs due to spacing formed between inclined metallic strips of delta loops as shown in Figure 6(a). The metallic strips printed on substrate are planar and not fully 3D. The accurate expressions for capacitances for coupled parallel strips deposited on dielectric substrate can be found from reference [36, 37].

The equivalent circuit is drawn at resonant frequencies from the plots of surface currents and electric field intensities (refer to Figure 7). The equivalent circuit presents qualitative response from the structure and helps to understand resonance mechanism at respective resonant frequencies. The application of vertically polarized plane wave induces voltages across parallel branches of the resonator and forms capacitances. The capacitor discharges and current flows via inductive metal path along the length of arrow-like geometrical portion as shown in Figure 7(a). The self-inductance of each metal strip is calculated by using (5) from the path where current flows [36].The terms , , and are strip width, strip thickness, and strip length, respectively, in microns. The mutual inductance between strips is calculated by usingwhere is the effective length of two metal strips and is the mean distance between the strips. Dimensions , , , , and are in microns. The capacitance between coplanar strips is calculated by using the following [36, 37]:where effective dielectric constant and filling factor is defined as follows: , , , , , , , , where represents a complete elliptical function of first kind with being its complementary function [36]. and are width and spacing between coplanar metal strips. The computation of resonant frequency requires calculation of total inductance and capacitance in the structure at respective resonant mode. At first resonant frequency, that is, GHz, self-inductance () for the current path as shown in Figure 7(a) has value of 18.5 nH calculated using (5). Mutual inductance between the metal strips having antiparallel current has value of 1.5 nH calculated using (6). Then, the total inductance is calculated as  nH. The capacitance between coupled parallel strips is  pF. The resonant frequency is calculated as  GHz which is very much close to simulated resonant frequency of 4.32 GHz. It is observed that, at this frequency, the effect of capacitor is dominant compared to intercell capacitance . This fact is verified by observing the effect of closing the gaps and variations in gap length over resonant frequency (refer to Figures 10 and 8). Note that we have neglected ohmic losses occurring due to finite conductivity of metal strips.

At second resonant frequency, that is, GHz, surface current is excited mainly along top and bottom horizontal strips of outer square loop resonator, left and right portions of delta loops as shown by arrows in Figure 7(b). The total inductance for the current path is  nH calculated from (5). The localized electric fields can be observed at the top right and bottom left horizontal branches of the square loop. These electric fields are due to capacitances formed between gaps of the adjacent unit cells denoted as intercell capacitance . The capacitance is also formed at the center within the unit cell . This capacitance is relatively weak compared to that of intercell capacitance (0.18 pF) and has pronounced effect on determining resonant frequency as they act in series. This fact is substantially similar to that of resonator discussed in [20]. Determining the value of capacitor is difficult as it is formed between nonadjacent diagonally opposite arms of delta loop. However, in calculation of , it is approximated as 1/6th the value of . The equivalent capacitance is then calculated as  pF. The predicted resonant frequency has value 7.51 GHz close to simulated value.

At third resonant mode , total self-inductance for the current path is  nH calculated from (5) and mutual inductance between the metal strips having antiparallel current has value of 1.5 nH calculated using (6). Then the total inductance is calculated as  nH. Capacitor is relatively weaker than and it is approximated as 1/8th that of . The equivalent capacitance is then calculated as  pF. The predicted resonant frequency has value 9.8 GHz close to simulated value. The increased resonant frequency as against is due to decreased value of inductance and capacitor .

4. Effect of Design Parameters on Resonance Frequency

Parametric analysis of the TBEMR is carried out in order to point out importance of different geometrical parameters on the LC resonant frequency and guidelines to design such material. Series of numerical simulations are performed concerning variations in the arm length of inner delta loops (, ), positions of shorting stub , (, ), closing the gaps formed due to arms of open delta loops, metal strip width , and substrate properties like permittivity. The transmission spectra under normal incidence are indicted and analyzed for each configuration as presented in subsequent subsection.

4.1. Effect of Changing Arm Length of Inner Delta Loops

The arm length for inner delta loops is varied from  mm to  mm in steps of 0.4 mm. For optimized design as shown in Figure 8, we have considered length  mm and three resonances occur at 4.32 GHz, 7.55 GHz, and 9.76 GHz, respectively. For decrement in value of from optimized value, we find that resonant frequency increases to 5.21 GHz (for  mm) from earlier values of 4.32 GHz (for  mm). There is almost a shift of 20% at resonant frequency . From simulated surface current density and electric field distribution, it is observed that resonant mode at is preserved when arm length is decreased. The increase in resonant frequency is due to decrease in capacitance , as decreases, whereas second resonant frequency shifts upward from value of 7.5 GHz (for  mm) to 8.62 GHz (for ). For further decrements in value of , the second and third resonance merge into single value at 9.22 GHz. It is found that resonating modes at and are preserved up to value of  mm. The shifts in values of and are probably due to weaker capacitance and . However for lower values of  mm and  mm modes at and merge and new resonating mode is formed.

4.2. Effect of Variations in Position of Shorting Stubs ,

The lengths and of both the ODLRs are varied to see the effect for the change in position of shorting stubs and as shown in the inset of the Figure 9. Parameters ( and ) are varied from optimized value of 3.6 mm to minimum value of 2.0 mm in steps of 0.4 mm. It is observed that, for optimized value of 3.6 mm, three resonances are obtained as shown in Table 1. For decreased value of 3.2 mm, resonances at and shift to higher frequencies. The shift at frequency is very small compared to shift at . For lowest value of 2.0 mm resonance at disappears and two closely spaced resonances and are observed at 5.4 GHz and 7.32 GHz, respectively.

4.3. Effect of Gaps on Formation of Resonance

To understand the effects of gaps formed due to open delta loops, transmission spectra is computed by closing gaps between and within delta loops as shown in the inset of Figure 10. When gaps between inclined arms of open delta loops are closed, it is found that electric resonance at totally disappears. This fact confirms the role of capacitance in the formation of resonance at . However, the other two resonances at and are still present and they are shifted to lower frequencies by 6.13% and 5.1%. The resonant frequencies at and are preserved because capacitances and are preserved as they are formed between diagonally opposite arms of the delta loops. These facts are confirmed from the plots of the surface current distributions plotted at resonant frequencies (not shown here). The resonances at and are present, when gaps within the delta loops are closed in a manner as shown in inset of Figure 10. The resonance at frequency disappears totally. The resonance frequency at is lowered by 3.33% and at is lowered by 5.1%. These results simply imply that resonant frequencies at and are directly related to the gaps between and within delta loops, while is slightly affected by closing the gaps in either way. The ENG response at first and second band can be tuned by closing corresponding gaps. Moreover, this study also paves the way for designing dual and triple-band response and switching the bands at first and second resonance.

4.4. Effect of Changing Substrate Properties

The material properties such as permittivity of the substrate affect resonant frequency and effective medium parameters such as relative permittivity , permeability , and refractive index . To study this effect, we have taken different low and high permittivity substrates like Rogers Duroid 5880 , FR-4 , Rogers TMM4 , Alumina , and Taconic CER-10 . Considering normal wave incidence setup, it can be found that, without any geometrical and dimensional changes in the structure, a high permittivity substrate presents a lower resonant frequency. The downward shift in resonant frequency can be clearly observed for high permittivity Alumina, Taconic CER-10 substrate as shown in Figure 11. This happens because the capacitance formed between coplanar metallic strips in the structure increases with increase in the permittivity of the substrate, thereby lowering resonant frequency. Another advantage of using high permittivity substrate is that we get improvement in at respective resonant frequency (refer to Table 2) which is figure of merit for miniaturization and homogenization. We also find that more than three resonances, that is, four, come into picture for high permittivity Alumina and Taconic CER-10. Although high permittivity substrates are useful in getting resonant response towards lower frequency, its cost and availability have to be considered. The above study can be employed to optimize a structure for various applications and to get a desired electromagnetic response. Moreover, it is also observed that electric field coupling (dip in ) improves for low loss substrate. Low loss Rogers TMM4 has improved coupling compared to that of high loss FR-4 at respective resonant frequencies. Similar observations can be made while comparing Alumina and Taconic CER-10.

4.5. Effect of Changing Metallic Strip Width

For optimized design, we have considered width for metallic strips as  mm. For parametric analysis, width is varied from 0.18 mm to 0.28 mm in steps of 0.02 mm. As seen from plots shown in Figure 12, the variations in resonant frequencies at , , and are less than 2%, 3.5%, and 4%, respectively. It can be inferred that resonant frequencies do not vary much. However slight increase observed at and can be attributed to decrease in self-inductance with increment in metallic strip width [13]. Moreover, it is reported that for very thin metal width conductor losses increase [36].

Further, performance comparison of the proposed structure with previously reported metamaterial structures is carried out as shown in Table 3. For comparison, we have considered unit cells operating in microwave regime. From Table 3, it is indicated that, compared to [17, 18, 20, 26], the proposed TBEMR structure exhibits triple-band response and highest value of ratio at lower resonant frequency of 4.32 GHz. Moreover at other two resonant frequencies, its value more than 4 is maintained. The proposed structure has advantage of having purely electric response as against [10, 28]. Letters E and M in Table 3 denote type of resonance, namely, electric and magnetic, respectively. An improved ratio, miniaturization, and triple-band response are the main characteristics of the proposed structure.

5. Conclusion

In conclusion, we present new compact triple-band microwave metamaterials evolved by integrating open delta loops within square ring resonator. The scattering parameters and effective medium parameters of the proposed metamaterial are investigated through full-wave simulation over the frequency range 3 to 12 GHz. Proposed geometry has purely electric response as indicated by negative permittivity property at respective resonant frequencies. The effect of various geometrical design parameters on resonant frequencies is investigated and it serves as a guideline to get dual, triple, and tunable response from the proposed structure. Moreover, proposed TBEMR resonator is compact in size as indicated by improved wavelength to unit cell ratio having value of 11.5 and 5.1 at lower and higher resonant frequency, respectively. The proposed metamaterial presents improvement over conventional single-band ENG structures such as thin wire and complementary SRR. The proposed TBEMR is expected to find applications in realization of multiband devices such as antennas, filters, and wave absorbers at microwave and higher frequency regime.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.