Modified Coptic Cross Shaped Split-Ring Resonator Based Negative Permittivity Metamaterial for Quad Band Satellite Applications with High Effective Medium Ratio

This research article describes a modified Coptic cross shaped split ring resonator (SRR) based metamaterial that exhibits a negative permittivity and refractive index with a permeability of nearly zero. The metamaterial unit cell consists of an SRR and modified Coptic cross shaped resonator providing quadruple resonance frequency at 2.02, 6.985, 9.985 and 14.425 GHz with the magnitude of −29.45, −25.44, −19.05, and −24.45 dB, respectively. The unit cell that was fabricated on a FR-4 substrate with a thickness of 1.6 mm has an electrical dimension of 0.074λ × 0.074λ; the wavelength (λ) is computed at the frequency of 2.02 GHz. The computer simulation technology (CST) microwave studio was employed to determine the scattering parameters and their effective medium properties, i.e., permittivity, permeability and refractive index, also calculated based on NRW (Nicolson–Ross–Weir) method through the implementation of MATLAB code. The frequency range of 2.02–2.995 GHz, 6.985–7.945 GHz, 9.985–10.6 GHz, and 14.425–15.445 GHz has been found for negative permittivity. An effective medium ratio (EMR) of 13.50 at 2.02 GHz shows that the proposed unit cell is compact and effective. The lumped component based equivalent circuit model is used to validate with simulation results. The proposed unit cell and its array were fabricated for experimental verification. The results show that the simulation result using CST and high-frequency structure simulator (HFSS) simulator, equivalent circuit model result using advanced design system (ADS) simulator and measurement results match each other better. Its near zero permeability, negative permittivity, negative refractive index, high EMR and simple unit cell design allow the proposed metamaterial to be used for S-, C-, X- and Ku-band satellite applications.


Introduction
The beginning of the metamaterial-based research using theoretical and experimental demonstration was considered by John Pendry and David R. Smith in 2000, after ignoring Veselago's idea in a similar case for more than 30 years [1][2][3]. Metamaterials are outlined as a man-made periodic structure designed to alter the electromagnetic behavior of materials to achieve excellent properties not found in nature. Metamaterials can have different types, for example single negative (ε or µ), double negative (ε and µ), and a double positive. Metamaterials have become an advanced technology due to their ability to obscure their identities and manipulate light. With the crystallization of this new technology, further applications have originated in radar and satellite communications [4,5], sensing applications [6,7], microwave and terahertz absorber [8][9][10][11], superlens [12], cloaking devices [13,14], seismic protection [15], and specific absorption rate (SAR) detection [16,17].
A double negative (DNG) metamaterial was proposed in [18] based on the modified split H-shape that exhibits single resonance frequency of 6.296 GHz and an EMR of 6.7 fabricated on Rogers RO 3010 substrate material. In [19], an epsilon negative (ENG) metamaterial based on crossed line SRR printed on Rogers RT 5880 dielectric material was introduced that presents triple band resonance frequency covering C-, X-and Ku-band microwave applications. This metamaterial's unit cell size is 10 × 10 mm 2 and an EMR of 4.5. In [20], a double C shaped metamaterial for microwave applications was introduced that displays four resonance frequencies at 2.37, 6.12, 8.16, and 14.51 GHz. This metamaterial shows a good EMR of 10.55. A DNG metamaterial based on modified Z-shaped was presented in [21] that is suitable for wideband microwave applications, showing a dual resonance frequency at 7.32 GHz and 11.84 GHz with poor EMR of 4.09. In [22], a DNG metamaterial was proposed based on double C shape that exhibits multiband resonance covering S-, C-, and X-band with EMR of 7.44. A DNG metamaterial based on the double Z-shape was outlined in [23] that presents a triple resonance frequency at 7.3, 8.1, and 9.4 GHz. This metamaterial unit cell size is 8.5 × 8.5 mm 2 showing an EMR of 4.8. The negative permittivity metamaterial based on double dumbbell shaped SRR was suggested in [24], that is fabricated on Rogers RT 6002 substrate material. This metamaterial exhibits five resonance frequency covering the S-, X-, and Ku-bands and also shows an excellent EMR of 11.51. The electrical dimension of this metamaterial is 0.087λ × 0.087λ, where λ is calculated at the first resonance frequency. Aziz et al. reported the single negative (SNG) metamaterial structure based on concentric ring crossed lines that employed Kuand Ka-bands, fabricated on FR-4 substrate with an EMR of 4.44 [25]. In another study [26], a DNG metamaterial with size of 5.5 × 5.5 mm 2 was introduced that covered C-band microwave applications. This metamaterial exhibits a moderate EMR of 8. Another study demonstrated a metamaterial cloaking operation based on the epsilon negative, in [27]. This cloaking operates in a C-band microwave region and hides the metal cylinder. A hexagonal split-ring resonator with a coupled aperture is proposed as a microwave absorber in [28]. This structure exhibits a dual resonance frequency at 3.56 GHz and 11.67 GHz that is fabricated on a FR-4 substrate covering C and X-band microwave applications. An EMR of this resonator is 8.4. A triple band double negative metamaterial based on the double L-shape was proposed in [29]. This resonator shows a quadruple resonance frequency at 7.69, 8.47, 13.14, and 12.04 GHz and also presents an EMR of 3.9. In addition, [30] presented a double H-shaped SRR based metamaterial that exhibits quadruple resonance for satellite application. The dimension of the unit cell is 9 × 9 mm 2 , showing an EMR of 10.75. In other research, a cross-coupled interlinked based negative permittivity metamaterial was introduced in [31]. This metamaterial displayed triple resonance at 4.15 GHz, 10.38 GHz, and 14.93 GHz that also shows an EMR of 8.03. This work was also validated through the equivalent circuit model and measurement results. The effect of the dielectric properties of biological materials on the propagation of electromagnetic waves was presented in [32]. This work also delineates the technique of implementing the negative permittivity-based metamaterial structure on a microwave applicator that is enhanced by the control of the electromagnetic field in microwave therapy. In recent years, great progress has been made in the field of SPP. Much SPP research has been conducted, especially for metamaterials and photonic crystals, as they can dramatically adjust the performance of SPPs [33]. The above literature review is explained through the explanation of metamaterial behavior with multiband characteristics. However, high EMR with multiband resonance exhibiting metamaterial structure design still needs to be improved.
Since metamaterials consider their EMR to be inadequate in many of the applications discussed above, this research outlines a modified Coptic cross shaped split ring resonator with negative permittivity metamaterial for quad band satellite applications. This metama-terial is fabricated on low cost FR-4 substrate material showing resonances at 2.02, 6.985, 9.985, and 14.425 GHz that cover S-, C-, X-and Ku-bands with an EMR of 13.50. The negative permittivity was found to be 2.02-2.995 GHz, 6.985-7.945 GHz, 9.985-10.6 GHz, and 14.425-15.445 GHz, respectively. The novelty of this work is that the proposed design can provide a unique shape with miniaturized size showing a dimension of 0.074λ ×0.074λ. The proposed metamaterial produced a quadruple resonance covering with four bands including S-, C-, X-, and Ku-band satellite applications. Moreover, it can provide wider bandwidth, i.e., 0.385 GHz (1.815-2.2 GHz), 0.92 GHz (6.41-7.33 GHz), 0.56 GHz (9.66-10.22), and 1.39 GHz (13.57-14.96 GHz). The developed modified Coptic cross shape exhibited a high effective medium ratio (EMR) of 13.50 that showed compactness, which is relatively high in this regard compared to the recently proposed metamaterials. The array results are almost like a unit cell due to low mutual coupling effect. A unique characteristic of this structure is that the resonance frequency can be tuned by modifying the sizes of the Coptic width. A widely used NRW method is used to achieve the effective medium properties. E-field, H-field, and surface current were also analyzed to understand the behavior of the metamaterial. In addition, this work describes the equivalent circuit modeling using ADS simulator and also analyzes the measurement results through the fabrication of unit cell and array of unit cells.

Metamaterial Structure Design and Simulation Geometry
The proposed metamaterial unit cell is represented on an inexpensive FR-4 dielectric material (PCBWay, Xiacheng, China) with a thickness of 1.6 mm, a physical size of 11 × 11 mm 2 , and electrical dimension of 0.074λ × 0.074λ. The FR-4 substrate has a dielectric constant of 4.4 with a dissipation factor of 0.02. The resonator of the unit cell is made of annealed copper with a thickness of 0.035 mm and a conductivity of 5.96 × 10 7 S/m. The metamaterial resonator consists of a split ring resonator and modified Coptic cross shaped resonator, as indicated in Figure 1. The top view of the unit cell represents a resonating patch of the metamaterial, as shown in Figure 1a; on the other hand, the isometric view shows the annealed copper resonating patch with FR-4 substrate, as presented in Figure 1b. The length and width of the SRR are 10.4 mm and 1 mm, respectively. The split gap of the SRR is 0.2 mm, which is fixed by the trial and error method. The split gap of the modified Coptic cross is 0.4 mm. The design specification of the metamaterial is presented in Table 1. The time domain solver with hexahedral meshing system is used to derive the scattering parameters. The simulation scenario is displayed in Figure 2, where a transverse electromagnetic wave (TEM) is usually propagated through the proposed unit cell. As shown in Figure 2, PEC-PMC boundary geometry is used for simulation. The metamaterial unit cell and arrays of unit cells are placed between the two waveguide ports (one is the transmitting end, the other is the receiving end).

Effective Medium Parameters Extraction Method
The simulation approach of the proposed modified Coptic cross shape-based unit cell was executed using a time domain solver with hexahedral mesh in the frequency span of 1 to 16 GHz. The NRW technique was used to obtain the relative permittivity, relative permeability, refractive index, and impedance [34]. The equations were as follows to determine the effective parameters: where, V 1 and V 2 are the composite terms, which is derived using transmission and reflection coefficient of the unit cell; k 0 = 2π f c , V 1 = |S 21 | + |S 11 |, V 2 = |S 21 | − |S 11 |, and d is the substrate thickness; and ε r , µ r , η r , and Z are the relative permittivity, relative permeability, refractive index, and normalized impedance, respectively.

Design Procedure of the Metamaterial Unit Cell
The proposed design was chosen to analyze different unit cell design approaches based on a trial and error method, as shown in Figure 3, and the S 21 and S 11 plot for the different design approaches is shown in Figure 4a,b. Multiple resonance frequencies require more factors, i.e., the size and thickness of the substrate, width and split gap of the resonating patch, and mutual coupling of the resonator. The variation in capacitance and inductance of the unit cell is responsible for creating a multiple resonance frequency. The unit cell design was initiated with the insertion of a split ring resonator through the split gap of 0.2 mm. Three resonance frequencies were obtained from design 1, the resonance frequencies of reflection coefficients (S 11 ) are 3.325, 9.05, and 11.5 GHz with the magnitude of −27.78, −9.88, and −11.04 dB, respectively, whilst the resonance frequencies of transmission coefficients (S 21 ) are 2.05, 7.19, and 9.47 GHz with the magnitude of −29.42, −28.92, and −11.25 dB, respectively. In the next step, a Coptic shape is added to the previous design showing triple resonances. This inclusion improves the inductance value and thus creates a similar resonance with the previous design, but increased the magnitude by slightly changing the first and second resonance frequencies to a lower frequency range. The resonance frequencies of S 11 are at 2.995, 7.39, and 15.97 GHz with the magnitude of −27.75, −12.15, and −3.57 dB; the S 21 is at 2.02, 6.79, and 12.91 GHz with the magnitude of −28.48, −22.41, and −34 dB, respectively. According to design 3 of Figure 3, the same resonance frequency exists as the earlier design by inserting four cross shapes into design 2, because of only improving the inductive property. The magnitudes of the resonance frequencies of S 11 are −27.52, −12.06, and −4.05 dB at 2.98, 7.39, and 15.8 GHz, respectively, whereas, the resonance frequencies of S 21 are 2.02, 6.79, and 12.88 GHz with the magnitude of −28.97, −22.45, and −34 dB, respectively. The design is completed by inserting two split gaps in Coptic shape, as shown in Figure 3. The inserting split gap in design 3 creates additional resonance frequency due to generating a capacitance property. According to Figure 4, the resonance frequencies of S 11 are 3.085, 8      The S11 plot on different design approaches.

Parametric Study
The performance of the metamaterial unit cell depends on several factors such as the split gap of the ring, the width of the metamaterial resonator, substrate and its thickness, and boundary condition of the simulation environment. These factors directly hit the variation of inductance and capacitance property, which results in a change in resonance frequency of the metamaterial. Numerous parametric studies have been executed to investigate the effects of the split gap of the SRR and Coptic cross resonator, the width of the Coptic cross resonator, effect of different conductor on metamaterial resonator, effect of substrate material, and boundary condition of the simulation geometry.
Effect of split gap of SRR on S21 results: The S21 plot for the various split gap of SRR is introduced in Figure 5a. The split gap of the SRR initiates a capacitance variation that influences the resonance frequency of the metamaterial. The capacitance value decreases with the increasing split gap of the ring, which affects the resonance frequency to a higher value. The initial step introduces a split gap of 0.2 mm, exhibiting quadruple resonances with excellent magnitude. Then, the split gap was changed to 0.4 mm and 0.6 mm, showing the same response as the result of 0.2

Parametric Study
The performance of the metamaterial unit cell depends on several factors such as the split gap of the ring, the width of the metamaterial resonator, substrate and its thickness, and boundary condition of the simulation environment. These factors directly hit the variation of inductance and capacitance property, which results in a change in resonance frequency of the metamaterial. Numerous parametric studies have been executed to investigate the effects of the split gap of the SRR and Coptic cross resonator, the width of the Coptic cross resonator, effect of different conductor on metamaterial resonator, effect of substrate material, and boundary condition of the simulation geometry.

Effect of split gap of SRR on S 21 results:
The S 21 plot for the various split gap of SRR is introduced in Figure 5a. The split gap of the SRR initiates a capacitance variation that influences the resonance frequency of the metamaterial. The capacitance value decreases with the increasing split gap of the ring, which affects the resonance frequency to a higher value. The initial step introduces a split gap of 0.2 mm, exhibiting quadruple resonances with excellent magnitude. Then, the split gap was changed to 0.4 mm and 0.6 mm, showing the same response as the result of 0.2 mm split gap of SRR, but the resonance frequency moves to higher value. Finally, the no split gap is inserted inside the ring, which affects the resonance frequency, where the first resonance is moved to 4.14 GHz, and the second resonance frequency is almost diminished. It can be seen from Figure 5a that the resonance frequency for 0.2 mm split gap is better than other specified values. Finally, a 0.2 mm split gap is selected for the proposed design.

Effect of Coptic shape resonator width on S 21 results:
The S 21 plot for the different widths of the Coptic shape resonator is outlined in Figure 5c. As the Coptic shape resonator's width increases, the overall inductance property decreases and the resonance frequency increases, because the resonance frequency is inversely proportional to the inductance. The width of the Coptic shape resonator was taken in four sizes, such as 0.4, 0.6, 0.8, and 1 mm. The effect of width analysis found that the first and second resonance frequencies were almost the same for different widths, but the third and fourth resonance frequencies increased slightly. In this analysis, the width was set to 0.6 mm for the proposed design. mm split gap of SRR, but the resonance frequency moves to higher value. Finally, the no split gap is inserted inside the ring, which affects the resonance frequency, where the first resonance is moved to 4.14 GHz, and the second resonance frequency is almost diminished. It can be seen from Figure 5a that the resonance frequency for 0.2 mm split gap is better than other specified values. Finally, a 0.2 mm split gap is selected for the proposed design.

Effect of split gap of Coptic shape resonator on S21 results:
The S21 plot for the various split gaps of the Coptic shape resonator is shown in Figure  5b. The capacitance property depends on the split size of the resonator. The first design starts with a 0.2 mm split gap of the Coptic shape resonator showing quadruple resonance frequency with excellent magnitude. The split gap was uniformly increased by 0.2 mm. An increase in the value of the resonance frequency was also observed. Lastly, no split gap was inserted into the Coptic shape resonator that affects the resonance frequency, where the third resonance frequency was shifted from 10.3 to 12.895 GHz and the fourth resonance disappeared. As shown in Figure 5b, it can be realized that the resonance frequency with magnitude is better than the other specified values for a Coptic shape resonator split gap of 0.4 mm. Finally, a split gap of 0.4 mm was chosen for the proposed design.

Effect of split gap of Coptic shape resonator on S 21 results:
The S 21 plot for the various split gaps of the Coptic shape resonator is shown in Figure 5b. The capacitance property depends on the split size of the resonator. The first design starts with a 0.2 mm split gap of the Coptic shape resonator showing quadruple resonance frequency with excellent magnitude. The split gap was uniformly increased by 0.2 mm. An increase in the value of the resonance frequency was also observed. Lastly, no split gap was inserted into the Coptic shape resonator that affects the resonance frequency, where the third resonance frequency was shifted from 10.3 to 12.895 GHz and the fourth resonance disappeared. As shown in Figure 5b, it can be realized that the resonance frequency with magnitude is better than the other specified values for a Coptic shape resonator split gap of 0.4 mm. Finally, a split gap of 0.4 mm was chosen for the proposed design.
Effect of different resonating conductor on S 21 results: The effect of different resonating conductors on the resonance frequency of S 21 is described in Figure 5d. The different metals such as gold, nickel, platinum, tantalum, and copper are used to test the performance of the metamaterial. First, copper is used as a patch to indicate the quadruple resonance frequency with excellent bandwidth. The effect of a resonating conductor (gold, platinum, tantalum) indicated that the response of S 21 was the same as copper, whereas nickel used as a patch showed the same response with lower magnitude. Since gold, platinum, and tantalum are very expensive metals, annealed copper was used as a resonating conductor for the proposed design.

Effect of substrate material on S 21 response:
Three different substrate materials were used for substrate selection. The three substrates are FR-4 (lossy and loss-free), Rogers RT-5880 (lossy), and Rogers RT-6002 (lossy). All substrate properties are different. In the first analysis, the dielectric constant, dielectric loss, and thickness (mm) of FR-4 (lossy) are 4.3, 0.025, and 1.6, respectively. Figure 6a shows that this material has demonstrated quadruple resonance frequency with the magnitudes of −29.45, −25.44, −19.05, and −24.45 dB; on the other hand, FR-4 (loss-free) material exhibits the same resonance with the magnitudes of −39.35, −44.55, −38.94, and −49.02 dB. The dielectric constant, dielectric loss, and thickness (mm) of Rogers RT-5880 (lossy) and RT-6002 (lossy) are 2.2, 0.0009, and 1.575 and 2.94, 0.0012 and 1.524, respectively. As shown in Figure 6a, both substrate materials exhibit triple resonance frequency with good magnitude. Compared to the FR-4 substrate, the resonance frequency of both Rogers materials moves to a higher value. However, the performance of the FR-4 substrate material is better than the other substrate materials. Therefore, the FR-4 substrate material was selected solely because it was consistent with the purpose of this study.
patch to indicate the quadruple resonance frequency with excellent bandwidth. The effect of a resonating conductor (gold, platinum, tantalum) indicated that the response of S21 was the same as copper, whereas nickel used as a patch showed the same response with lower magnitude. Since gold, platinum, and tantalum are very expensive metals, annealed copper was used as a resonating conductor for the proposed design.
Effect of substrate material on S21 response: Three different substrate materials were used for substrate selection. The three substrates are FR-4 (lossy and loss-free), Rogers RT-5880 (lossy), and Rogers RT-6002 (lossy). All substrate properties are different. In the first analysis, the dielectric constant, dielectric loss, and thickness (mm) of FR-4 (lossy) are 4.3, 0.025, and 1.6, respectively. Figure 6a shows that this material has demonstrated quadruple resonance frequency with the magnitudes of −29.45, −25.44, −19.05, and −24.45 dB; on the other hand, FR-4 (loss-free) material exhibits the same resonance with the magnitudes of −39.35, −44.55, −38.94, and −49.02 dB. The dielectric constant, dielectric loss, and thickness (mm) of Rogers RT-5880 (lossy) and RT-6002 (lossy) are 2.2, 0.0009, and 1.575 and 2.94, 0.0012 and 1.524, respectively. As shown in Figure 6a, both substrate materials exhibit triple resonance frequency with good magnitude. Compared to the FR-4 substrate, the resonance frequency of both Rogers materials moves to a higher value. However, the performance of the FR-4 substrate material is better than the other substrate materials. Therefore, the FR-4 substrate material was selected solely because it was consistent with the purpose of this study.

Effect of substrate length on S21 response:
The design starts with an initial selection of basic dimensions and considers the same general subwavelength criteria as λ/10. The S21 plot for different substrate lengths is exhibited in Figure 6b. Since the substrate length and resonance frequency are inversely proportional to each other, we examined the performance of four different unit cell substrate length, such as 14, 13, 12, and 11 mm, in order to obtain the lowest resonance frequency

Effect of substrate length on S 21 response:
The design starts with an initial selection of basic dimensions and considers the same general subwavelength criteria as λ/10. The S 21 plot for different substrate lengths is exhibited in Figure 6b. Since the substrate length and resonance frequency are inversely proportional to each other, we examined the performance of four different unit cell substrate length, such as 14, 13, 12, and 11 mm, in order to obtain the lowest resonance frequency in the lower dimension as the higher EMR. For the effect of substrate length on S 21 response, it was found that first and second resonance frequency shifted to higher value but third and fourth resonance frequency decreased slightly. The unit cell with four different substrate lengths exhibits the same resonance frequency. Finally, the proposed metamaterial resonator with a substrate length of 11 mm on each side is selected to achieve quadruple resonance frequency covering for S-, C-, X-, and Ku-band.
Effect of different boundary condition on S 21 response: The aim of this study is to elucidate the effects of different boundary conditions on modified Coptic cross shape-based metamaterial in the frequency range of 1 to 16 GHz. Five different boundary conditions, perfect electric conductor (PEC)-perfect magnetic conductor (PMC), PMC-PEC, PEC-PEC, PMC-PMC, and Unit Cell, are used to describe the metamaterial performance. The proposed metamaterial is used to demonstrate that the use of different boundary conditions can lead to completely different electromagnetic responses, as shown in Figure 6c. In the first study, PEC-PMC (E-field in x-axis and H-field in y-axis) boundary conditions are used in modified Coptic cross shape-based metamaterial; there are quadruple resonance frequencies in a given frequency range, representing 2.02, 6.985, 9.985, and 14.425 GHz. In the second study, PMC-PEC (H-field in x-axis and Efield in y-axis) boundary conditions are utilized in the proposed metamaterial exhibiting five resonance frequencies with smaller magnitude compared to PEC-PMC conditions. When PEC-PEC (E-field in x-axis and y-axis) conditions are employed in the proposed metamaterial, three non-uniform resonance frequencies arise in a given frequency range. In the orientation of PMC-PMC (H-field in x-axis and y-axis), conditions apply in both x and y directions. As seen from Figure 6c, these boundary conditions have a completely different response from the other boundary conditions. Lastly, unit cell boundary conditions (unit cells repeat periodically in x-axis and y-axis) are used to see the response of metamaterial. This boundary condition shows only one resonance peak. However, only the PEC-PMC (Efield in x-axis and H-field in y-axis) boundary condition was preserved because it fulfilled the purpose of this study.

Equivalent Circuit of the Proposed Unit Cell
The equivalent circuit model of the proposed unit cell is made up of passive components, such as inductor (L) due to metallic resonator patch, and capacitor (C) due to split gap of the ring resonator. The split gap of the ring with metallic resonator creates LC resonance circuits that supply the resonance frequency at a certain frequency, which is changed by the variation in the split gap and metallic strip of the metamaterial resonator [35].
The overall inductance was computed according to the principle of the transmission line [36]: Furthermore, the overall capacitance was calculated using: where ε 0 = 8.854 × 10 −12 F/m, µ 0 = 4π × 10 −7 H/m, w = metallic resonating patch width, h = substrate thickness, t = metallic resonating patch thickness, and l = length. The optimized equivalent circuit of the proposed metamaterial unit cell using the passive element (L, C) responsible for generating the quadruple resonance frequency is shown in Figure 7a. The component values are adjusted using an ADS simulator. By adjusting the component values, the S 21 results are computed so that the ADS results are very similar to the CST results. A comparison of S 21 is shown in Figure 7b. It was observed that the equivalent circuit results were expressed identically to those observed in unit cell results using CST in case of first to third resonance frequency. The metamaterial resonator consists of a split ring resonator and two split based Coptic cross shape resonator. The split ring resonator can be represented by L1, C1, L2, and C2, which is responsible for the creation of resonance frequency of 2.02 and 6.985 GHz. The Coptic shape with two split rings contributes to the creation of two resonance frequencies at 9.985 and 14.425 GHz through the representative L3, C3, L4, and C4 component of proper tuning. In the equivalent circuit model, L5 and L6 boost the magnitude of resonances representing the equivalent of the four cross shaped resonators. C5 represents the coupling capacitor.

Results and Discussion
The numerical simulation of the modified Coptic cross shape-based metamaterial unit cell was executed using a Finite Integration Technique (FIT) based CST simulator. The scattering parameters (S11 and S21) graph of the proposed unit cell is delineated in Figure 8a. The S11 results of the proposed unit cell provide quadruple resonances with frequencies of 3.

Results and Discussion
The numerical simulation of the modified Coptic cross shape-based metamaterial unit cell was executed using a Finite Integration Technique (FIT) based CST simulator. The scattering parameters (S 11 and S 21 ) graph of the proposed unit cell is delineated in Figure 8a Figure 9d shows that the measured transmission coefficient (S21) results (second, third, and fourth) are slightly inconsistent with the simulation results due to vector network analyzer (VNA) calibration errors, waveguide coupling effects, and cable losses. The first resonance frequency of measured results is more varied with simulated results. The proposed metamaterial structure has explored the tunability properties. Smaller Coptic resonator width sizes result in lower annealed copper inclusion, resulting in lower ohmic losses. Meanwhile, a large Coptic width size has several advantages in the attempt to obtain an infinite medium at higher frequency. This indicates that by changing the Coptic width, the metamaterial structure can be tailored for a variety of applications of interest. The influence of the Coptic width of the metamaterial on the transmission coefficient has already been investigated in Figure  5c. The Coptic width size ranges from 0.4 mm to 1 mm. The resonance frequencies shift higher as the Coptic width size increases. A unique characteristic of this structure is that the resonance frequency can be tuned by modifying the sizes of the Coptic width. Selecting the resonance frequency through a tuning process based on the size of the Coptic width makes it more suitable for satellite applications.  Figure 9d shows that the measured transmission coefficient (S 21 ) results (second, third, and fourth) are slightly inconsistent with the simulation results due to vector network analyzer (VNA) calibration errors, waveguide coupling effects, and cable losses. The first resonance frequency of measured results is more varied with simulated results. The proposed metamaterial structure has explored the tunability properties. Smaller Coptic resonator width sizes result in lower annealed copper inclusion, resulting in lower ohmic losses. Meanwhile, a large Coptic width size has several advantages in the attempt to obtain an infinite medium at higher frequency. This indicates that by changing the Coptic width, the metamaterial structure can be tailored for a variety of applications of interest. The influence of the Coptic width of the metamaterial on the transmission coefficient has already been investigated in Figure 5c. The Coptic width size ranges from 0.4 mm to 1 mm. The resonance frequencies shift higher as the Coptic width size increases. A unique characteristic of this structure is that the resonance frequency can be tuned by modifying the sizes of the Coptic width. Selecting the resonance frequency through a tuning process based on the size of the Coptic width makes it more suitable for satellite applications. The effective medium parameters (ε, µ, η, and Z) of modified Coptic cross shapebased negative permittivity metamaterial are extracted from the real and imaginary value of S11 and S21 using the NRW method, considered by Equations (1)-(4). The NRW-based MATLAB code is used to achieve the real and imaginary value of effective medium properties, as shown in Figure 10. Figure 10a shows a plot of relative permittivity versus frequency with negative permittivity (real) in the frequency range of 2.02 to 2.995, 6.985 to 7.945, 9.985 to 10.6, and 14.425 to 15.475. The imaginary values of effective permittivity are positive, and thus fulfill the criteria for metamaterial behavior. The effective permeability (real) values were near to zero at each resonance, with an amplitude of 0.048, 0.084, 0.16, and 0.12, respectively, as indicated in Figure 10b. The refractive index versus frequency plot is shown in Figure 10c, which shows the negative refractive index (real) in the 2.065 to 2.86 GHz, 7.09 to 7.765 GHz, 10.18 to 10.315, and 14.725 to15.04 GHz frequency ranges. On the other hand, the imaginary value of the refractive index exhibits a positive scenario in the same frequency range that reduces the incident wave in the structure. The impedance versus frequency plot is outlined in Figure 10d. It shows positive impedance (real) values over the entire frequency range. Impedance is observed close to zero, such as 0.006, The effective medium parameters (ε, µ, η, and Z) of modified Coptic cross shape-based negative permittivity metamaterial are extracted from the real and imaginary value of S 11 and S 21 using the NRW method, considered by Equations (1)-(4). The NRW-based MATLAB code is used to achieve the real and imaginary value of effective medium properties, as shown in Figure 10. Figure 10a shows a plot of relative permittivity versus frequency with negative permittivity (real) in the frequency range of 2.02 to 2.995, 6.985 to 7.945, 9.985 to 10.6, and 14.425 to 15.475. The imaginary values of effective permittivity are positive, and thus fulfill the criteria for metamaterial behavior. The effective permeability (real) values were near to zero at each resonance, with an amplitude of 0.048, 0.084, 0.16, and 0.12, respectively, as indicated in Figure 10b. The refractive index versus frequency plot is shown in Figure 10c, which shows the negative refractive index (real) in the 2.065 to 2.86 GHz, 7.09 to 7.765 GHz, 10.18 to 10.315, and 14.725 to15.04 GHz frequency ranges. On the other hand, the imaginary value of the refractive index exhibits a positive scenario in the same frequency range that reduces the incident wave in the structure. The impedance versus frequency plot is outlined in Figure 10d. It shows positive impedance (real) values over the entire frequency range. Impedance is observed close to zero, such as 0.006, 0.016, 0.04, and 0.029, at the resonance frequencies of 2.02, 6.985, 9.985, and 14.425 GHz, respectively, that represent the proposed metamaterial performing as a passive medium. The proposed metamaterial structure is effective or homogenized due to fulfilling the criteria for metamaterial behavior. The impedance (real), permittivity (imaginary), permeability (imaginary), and refractive index (imaginary) of the proposed metamaterial structure exhibits positive value, as shown in Figure 10. Moreover, the proposed metamaterial structure is used in S-, C-, X-, and Ku-band satellite applications. The C-band is used in satellite communication between the ground station and satellite. This band is also used in weather radar systems, terrestrial microwave links, and Wi-Fi devices. The S-and X-band are used in the weather radar system and communication satellite systems. Meanwhile, the Ku-band is commonly employed in satellite TV networks. Therefore, the metamaterial unit cell design outlined in this work can be effectively used in satellite. However, using this metamaterial structure offers various microwave products in S-, C-, X-and Ku-bands for satellite applications. This proposed metamaterial design can be used to make a variety of devices including antennas, bandpass filters, power dividers, and couplers. This metamaterial structure essentially enhances the performance of the satellite antenna. 0.016, 0.04, and 0.029, at the resonance frequencies of 2.02, 6.985, 9.985, and 14.425 GHz, respectively, that represent the proposed metamaterial performing as a passive medium. The proposed metamaterial structure is effective or homogenized due to fulfilling the criteria for metamaterial behavior. The impedance (real), permittivity (imaginary), permeability (imaginary), and refractive index (imaginary) of the proposed metamaterial structure exhibits positive value, as shown in Figure 10. Moreover, the proposed metamaterial structure is used in S-, C-, X-, and Ku-band satellite applications. The C-band is used in satellite communication between the ground station and satellite. This band is also used in weather radar systems, terrestrial microwave links, and Wi-Fi devices. The S-and Xband are used in the weather radar system and communication satellite systems. Meanwhile, the Ku-band is commonly employed in satellite TV networks. Therefore, the metamaterial unit cell design outlined in this work can be effectively used in satellite. However, using this metamaterial structure offers various microwave products in S-, C-, Xand Ku-bands for satellite applications. This proposed metamaterial design can be used to make a variety of devices including antennas, bandpass filters, power dividers, and couplers. This metamaterial structure essentially enhances the performance of the satellite antenna.

Surface Current, Magnetic Field, and Electric Field Analysis
In a metallic resonator, the surface current delineates the actual current generated by the time varying electromagnetic field. This actual current produces a magnetic field, so the changing magnetic field produces an electromotive force. Maxwell and conventional

Surface Current, Magnetic Field, and Electric Field Analysis
In a metallic resonator, the surface current delineates the actual current generated by the time varying electromagnetic field. This actual current produces a magnetic field, so the changing magnetic field produces an electromotive force. Maxwell and conventional law have proposed several equations to elucidate the relationship between surface current, E-field, and H-field through the following equations [37]: where E, H, D, B, ε, µ, and J are E-field intensities, H-field intensities, electric flux densities, magnetic flux densities, electric permittivity, magnetic permeability, and time-varying electric current density in a medium, respectively. Figure 11 displays the surface current, H-field, and E-field distribution of the proposed unit cell at the frequencies of 2.02, 6.985, 9.985, and 14.425 GHz.
At a frequency of 2.02 GHz, the outer SRR provides a low impedance path, thereby creating a strong surface current, as shown in Figure 11(ai). In Figure 11(bi), a good understanding was observed between surface current and H-field distribution, which was confirmed by Ampere's law ( B.dl = µ 0 I). As the current increases, the intensity of the H-field also increases. Strong surface current creates a higher H-field in same place. According to Maxwell's law, the greater the value of the H-field in SRR, the greater value of the E-field on the other side of the SRR due to capacitive effect is shown in Figure 11(ci).
At the resonance frequency of 6.985 GHz, the outer SRR produces a little bit of current. In the other portion, the surface current was evenly distributed, as shown in Figure 11(aii). The Coptic shape resonator also contributed less to the surface current flow due to increasing the capacitive response for split gap. On the other hand, the SRR shows a small strength of H-field; the Coptic cross shape resonator has almost no current due to the split gap and the H-field approaches zero, as shown in Figure 11(bii). Meanwhile, the split gap of the SRR and Coptic shape resonator create an additional electric field due to forming a capacitor, as indicated in Figure 11(cii).
At a frequency of 9.985 GHz, the entire unit cell resonator, except the split gap of the SRR and Coptic shape, four cross shaped resonators, contributed most of the surface current, as indicated in Figure 11(aiii). This higher surface current creates a strong H-field through the same place of the resonator, as shown in Figure 11(biii). As shown in Figure 11(ciii), an increased E-field is observed throughout the structure at this resonance frequency, except in the resonator where the H-field is larger.
At the resonance frequency of 14.425 GHz, lower current densities were recorded across the unit cell as the impedance increased with higher frequencies, as shown in Figure 11(aiv). From the Figure 11(biv), it can be seen that the concentration of H-field remains the same where the surface current is higher. Meanwhile, the E-field strength is high where the H-field is low as a result of creating a mutual capacitance in a metamaterial unit cell, as shown in Figure 11(civ). The scale and axis views are shown in Figure 11v for the distribution of surface current, H-field, and E-field.
where E, H, D, B, ɛ, µ, and J are E-field intensities, H-field intensities, electric flux densities, magnetic flux densities, electric permittivity, magnetic permeability, and time-varying electric current density in a medium, respectively. Figure 11 displays the surface current, H-field, and E-field distribution of the proposed unit cell at the frequencies of 2.02, 6.985, 9.985, and 14.425 GHz.  At a frequency of 2.02 GHz, the outer SRR provides a low impedance path, thereby creating a strong surface current, as shown in Figure 11(ai). In Figure 11(bi), a good understanding was observed between surface current and H-field distribution, which was confirmed by Ampere's law (

. =
). As the current increases, the intensity of the H-field also increases. Strong surface current creates a higher H-field in same place. According to Maxwell's law, the greater the value of the H-field in SRR, the greater value of the E-field on the other side of the SRR due to capacitive effect is shown in Figure 11(ci).
At the resonance frequency of 6.985 GHz, the outer SRR produces a little bit of current. In the other portion, the surface current was evenly distributed, as shown in Figure  11(aii). The Coptic shape resonator also contributed less to the surface current flow due to

Array Metamaterial Results
The metamaterial unit cell does not act alone to exhibit the correct exotic electromagnetic properties; some series of array designs were simulated. The array performance was analyzed by choosing a few array types, such as 1 × 2 array (11 × 22 mm 2 ), 2 × 2 array (22 × 22 mm 2 ), and 4 × 4 array (44 × 44 mm 2 ). This selected array was numerically simulated, using CST simulator as the same method of unit cell. Figure 12 shows the results of the scattering parameters for selected types of array. From Figure 12a, it can be seen that the S 11 result of all arrays shows the same response only for the first and second resonance frequency. The S 11 results of the unit cell and 1 × 2 array reveal the same response. The third and fourth resonance frequency of the 2 × 2 array and 4 × 4 array are shifted to the lower frequency due to the mutual coupling effect. Figure 12(b) presents the S 21 results of changing at each resonance frequency for selected types of array. The S 21 results of the unit cell and 1 × 2 arrays exhibit the same resonance frequency at 2.02, 6.985, 9.985, and 14.425 GHz with the magnitude of −29.45, −25.44, −19.05, and −24.45 dB, respectively. The third and fourth resonance frequencies shift lower for the 2 × 2 array and 4 × 4 array compared to the unit cell. The first and second resonance frequency for all the arrays exhibits the same response. The mutual coupling effect is responsible for this discrepancy between the resonance frequencies.

Array Metamaterial Results
The metamaterial unit cell does not act alone to exhibit the correct exotic electromagnetic properties; some series of array designs were simulated. The array performance was analyzed by choosing a few array types, such as 1 × 2 array (11×22 mm 2 ), 2 × 2 array (22×22 mm 2 ), and 4 × 4 array (44×44 mm 2 ). This selected array was numerically simulated, using CST simulator as the same method of unit cell. Figure 12 shows the results of the scattering parameters for selected types of array. From Figure 12a, it can be seen that the S11 result of all arrays shows the same response only for the first and second resonance frequency. The S11 results of the unit cell and 1 × 2 array reveal the same response. The third and fourth resonance frequency of the 2 × 2 array and 4 × 4 array are shifted to the lower frequency due to the mutual coupling effect. Figure 12(b) presents the S21 results of changing at each resonance frequency for selected types of array. The S21 results of the unit cell and 1 × 2 arrays exhibit the same resonance frequency at 2.02, 6.985, 9.985, and 14.425 GHz with the magnitude of −29.45, −25.44, −19.05, and −24.45 dB, respectively. The third and fourth resonance frequencies shift lower for the 2 × 2 array and 4 × 4 array compared to the unit cell. The first and second resonance frequency for all the arrays exhibits the same response. The mutual coupling effect is responsible for this discrepancy between the resonance frequencies.

Effective Medium Ratio and Comparison with Recent Study
The effective medium ratio (EMR) is the fundamental aspect of metamaterial related applications. EMR stands for the smallness and usefulness of metamaterials. Many devices operate at low resonant frequencies, but the small size of the metamaterial makes it difficult to achieve low resonant frequencies. A high EMR represents the integrity and achievement of metamaterial design standards. If the EMR is less than 4, the subwavelength criteria of the metamaterial are not met. The operating frequency and size of metamaterials are inversely related to each other. Therefore, the EMR should be considered to design the metamaterial. EMR is the ratio of wavelength to unit cell length and is expressed by Equation (13).
The modified Coptic cross shape-based metamaterial unit cell exhibits an EMR of 13.50, 3.90, 2.73, and 1.89 at the frequency of 2.02, 6.985, 9.985, and 14.425 GHz. An EMR

Effective Medium Ratio and Comparison with Recent Study
The effective medium ratio (EMR) is the fundamental aspect of metamaterial related applications. EMR stands for the smallness and usefulness of metamaterials. Many devices operate at low resonant frequencies, but the small size of the metamaterial makes it difficult to achieve low resonant frequencies. A high EMR represents the integrity and achievement of metamaterial design standards. If the EMR is less than 4, the subwavelength criteria of the metamaterial are not met. The operating frequency and size of metamaterials are inversely related to each other. Therefore, the EMR should be considered to design the metamaterial. EMR is the ratio of wavelength to unit cell length and is expressed by Equation (13).

EMR =
Wavelength, λ Unit cell length, L The modified Coptic cross shape-based metamaterial unit cell exhibits an EMR of 13.50, 3.90, 2.73, and 1.89 at the frequency of 2.02, 6.985, 9.985, and 14.425 GHz. An EMR of 13.50 at 2.02 GHz shows that the proposed unit cell is compact and effective. The main advantages of high EMR are: (1) improved characteristic uniformity, (2) shorter electrical dimension without fabrication limitations, and (3) reduced transmission coupling effects. The unit cell reveals an electrical dimension of 0.074λ × 0.074λ; λ is calculated at the frequency of 2.02 GHz which fulfills the criteria of high EMR condition.
The performance of the designed metamaterial is compared with recently published work considering unit cell size (physical and electrical), a range of resonance frequencies with covered bands, and the effective medium ratio of the metamaterial. The comparison is shown in Table 2. The EMR is an important aspect for metamaterial. The proposed metamaterial exhibits an excellent EMR of 13.50, compared to other works listed in Table 2. The table shows that Refs. [24,30,31,38,39] have the highest electrical dimensions and lower EMRs compared to the proposed work; the proposed work is, therefore, the most effective. The proposed metamaterial contains moderate size and quadruple resonance with higher EMR, which makes this design more fulfilling than others shown in the comparison Table 2.

Conclusions
In this article, numerical and measurement representation of a modified Coptic cross shape-based metamaterial is presented in which the unit cell and array structure was designed, simulated, and investigated using CST and HFSS simulator. The proposed unit is fabricated on a FR-4 substrate with an electrical dimension of 0.074λ × 0.074λ, and λ is calculated at the frequency of 2.02 GHz. The simulated resonance frequencies using CST and HFSS are 2.02, 6.985, 9.985, and 14.425 GHz; 2.05, 6.97, 9.98, and 14.215 GHz, respectively. The measured resonance frequencies using Agilent N5227 PNA Microwave Network Analyzer are 1.77, 7.181, 9.787, and 14.425 GHz. The simulated results cover the S-, C-, X-, and Ku-band; additionally, measured results incorporate the L-, C-, X-, and Ku-band. In addition, the effective permittivity shows as negative at all resonance frequencies, whilst effective permeability exhibits near zero characteristics that treat the metamaterial as an SNG.
Comparing simulation and measurement results shows strong agreement. The equivalent circuit model is sketched in the ADS simulator and the response of S 21 agrees very well with the CST and HFSS simulation results. This study was also rigorously performed using surface current, H-field, and E-field analysis, and the overall performance is superior and better than other work described in this article. The inexpensive FR-4 substrate and small electrical size provide the desired resonance frequency and exhibit a high EMR of 13.50. Therefore, the proposed modified Coptic cross shape-based metamaterial can be widely used for S-, C-, X-, and Ku-band satellite applications. S-band is employed for Mobile Satellite Service (MSS) networks, C-band can be used in commercial telecommunications via satellites, X-band satellite communication is utilized by military forces, and Ku-band can also be applied to satellite TV networks and satellite communication systems.