Research Article  Open Access
BroadsideCoupled Microstrip Lines as Low Loss Metamaterial for Microwave Circuit Design
Abstract
The entire microwave theory is based on Maxwell’s equations, whereas the entire electronic circuit theory is based on Kirchhoff’s electrical current and voltage laws. In this paper, we show that the traditional microwave design methodology can be simplified based on a broadsidecoupled microstrip line as a low loss metamaterial. That is, Kirchhoff’s laws are still valid in the microwave spectrum for narrowband signals around various designated frequencies. The invented low loss metamaterial has been theoretically analyzed, simulated, and experimentally verified in both time and frequency domains. It is shown that the phase velocity of a sinusoidal wave propagating on the low loss metamaterial can approach infinity, resulting in timespace shrink to a singularity as seen from the propagating wave perspective.
1. Introduction
Microwave theory and techniques, based on Maxwell’s equations, and electronic circuit theory and techniques, based on Kirchhoff’ electrical current and voltage laws, are traditionally two distinguishable disciplines in both university education and industrial research and development. It is known that Kirchhoff’s electrical current and voltage laws can be derived from Maxwell’s equations under a timeinvariant or slowly changing field condition [1]. Therefore, to use Kirchhoff’s laws in analog circuit design is truly valid at a relatively low frequency, e.g., below 300 MHz. At a higher frequency than 300 MHz, the socalled radio frequency (RF) design methodology has been utilized, which takes the inaccuracy of Kirchhoff’s laws into account. At a frequency of a few gigahertz and above, microwave theory and techniques are employed.
Now a challenging question is, can we find a way to accurately use Kirchhoff’s laws in the microwave and even mmwave spectrums? If the answer is yes, microwave designs that are difficult for many electronic engineers can be simplified. Moreover, new circuitry architectures and topologies may be worked out for future communication technologies. Theoretically, it is possible to converge Maxwell’s equations and Kirchhoff’s laws for microwave and even mmwave circuit designs, as shown in the next section. The precondition is to let both permittivity () and permeability (μ) be zero in Maxwell’s equations.
Metamaterial has become a research interest since early 2000s, when its existence was experimentally verified with manmade structures [2–4], even though a theoretical concept of metamaterial had been published much earlier in the literature [5]. The properties of metamaterial are characterized by its negative and zero permittivity and/or permeability. These properties do not exist with material in nature, since the permittivity and the permeability of freespace are positive values, i.e., ε_{0} = 8.854 x 10^{−12} F/m and μ_{0} = 4π x 10^{−7} H/m, whereas the permittivity and permeability values in a medium are larger than these values. To date, various metamaterial structures in one, two, and three dimensions have been proposed and studied [6–13]. However, previous published metamaterial structures have some common problems, e.g., high loss for electromagnetic wave propagation and difficulty in manufacturing metamaterial structures with precision at a high frequency.
In this paper we present a low loss metamaterial structure based on broadsidecoupled microstrip lines, which avoids the aforementioned problems. Moreover, it is compatible with traditional transmission line design, e.g., microstrip lines on a planar structure. Based on this low loss metamaterial structure, a detailed study around the transition region from left to righthanded wave propagations on the metamaterial has been done. It is shown that at various designated frequencies, the phase velocity of a sinusoidal wave approaches infinity, leading to timespace shrink to a singularity and resulting in zero permittivity and zero permeability of the metamaterial. Consequently, a simplified design methodology based on Kirchhoff’s laws can be utilized for microwave circuitry analyses, when the low loss metamaterial is utilized at various designated frequencies above 1 GHz.
2. Model and Theory
We start from the first principle, i.e., from Maxwell’s equations, and let both permittivity and permeability be zero, and then we study our proposed metamaterial structure with broadsidecoupled microstrip lines.
2.1. Maxwell’s Equations When Both ε and μ Are Zero
The four Maxwell’s equations in differential form under the condition of null magnetic current source are described below [14].where is the magnetic field intensity, the electric field intensity, the electric charge density, permittivity of a medium, μ permeability of a medium, and the time variable. From (1) and (2),where is the relative permittivity (dielectric constant) and the permittivity of freespace, whereas is the relative permeability and the permeability of freespace. Using (5),which is equivalent to the following equation in the integral form according to the divergence theorem [14]. where is an electrical current, flowing either in or out at a node.
Similarly, from (6) the following form can be obtained according to the Stokes’ theorem [14]: where is the electrical voltage between two points on a loop.
It is apparent that Kirchhoff’s current law, i.e., (8), can be derived from Maxwell’s equations under the precondition of ε = ε_{r}ε_{o} = 0. Similarly, Kirchhoff’s voltage law, i.e., (9), is derived under the precondition of μ = μ_{r}μ_{o} = 0. Note that (3) and (4) are not needed to obtain (8) and (9). However, from (3) ρ = 0 when ε = ε_{r}ε_{o} = 0, indicating no electric charge accumulation. When ε deviates from null, ρ becomes nonzero.
Note that in the above derivations, and in (1) and (2) can be nonzero, which means that the above derivations are valid under timevariant electrical and magnetic field conditions. This means that, when both ε and μ are zero, Kirchhoff’s laws are valid not only under a static or slowly changing field condition as in the standard circuit theory for a direct current (dc) or an alternate current (ac) circuitry, but also under fastchanging field conditions as in a microwave or a mmwave circuitry.
2.2. BroadsideCoupled Microstrip Line Model
Figures 1(a) and 1(b) show obliquetop and side views of our proposed metamaterial structure with two series broadsidecoupled microstrip lines and one shunt shortcircuit stub with multiple vias for good grounding. Figure 1(c) illustrates that two or more unit cells are cascaded. Dimensions of the structure designated for 1, 2, and 3 GHz are listed in Table 1. The substrate thickness is t_{1} = 0.17 mm and t_{2} = 1.52 mm, respectively. The substrate material used is Rogers RO4350B with a dielectric constant = 3.66 ± 0.05 and a loss factor of 0.0031 @ 2.5 GHz (0.0037 @ 10 GHz).

(a)
(b)
(c)
Figure 2(a) depicts a lossless equivalent model of the metamaterial unit cell shown in Figures 1(a) and 1(b), where C_{p} is a series capacitance of the broadsidecoupled line segment, either to the left or to the right of the middle shortcircuit stab having a shunt inductance of L_{s}. The shunt capacitance of the broadsidecoupled line segment to ground is C_{pg}. To simplify this model for analysis, we can choose the middle unit section between two cascaded unit cells with a length of ∆z, as illustrated in Figure 1(c). The equivalent model is then simplified to Figure 2(b). Apparently, Figure 2(b) shows a typical lossless metamaterial model with a composed left and righthanded transmission line. In our case, we have used only microstrip lines to realize the metamaterial, avoiding any discrete capacitor or inductor, e.g., interdigital capacitor used in [4]. This can be an advantage, since the design of the metamaterial is then scalable with frequency, similar to traditional microwave circuit design using microstrip lines. Moreover, our previous study [15, 16] has shown that low loss and broadband properties can be achieved with broadsidecoupled microstrip lines.
(a)
(b)
Using the model shown in Figure 2(b), the series impedance Z and shunt admittance Y of a unit section of length are described below. where and are unit length series inductance and unit length shunt capacitance of the broadsidecoupled line, whereas is the total series capacitance of the broadsidecoupled line with a length of ∆z. = is the total shunt inductance of the shortcircuit stub, and is an angular frequency, whereas f is a frequency. Note that in (10) and (11), we have used the total inductance and total capacitance for L_{L} and C_{L}, instead of unit length capacitance and unit length inductance. The reason is that, unlike and , L_{L} and C_{L} are not linearly scalable with ∆z that is the minimum length of a metamaterial section at a designated frequency.
It should be pointed out that, in a conventional transmission line theory [17], ∆z should be infinitesimally small. However, for the metamaterial structure shown in Figure 1, ∆z, i.e., l_{a} in Figure 1(a), is dependent on the shunt stub length l_{b}. As seen in Table 1, when l_{b} = 12.3 mm, the unit cell length l_{a} = ∆z = 9.80 mm at 1 GHz that has a wavelength λ of 300 mm in freespace and 157 mm when = 3.66 is considered. Thus, the ration of ∆z/λ is 0.033 to 0.062, indicating ∆z is reasonably small as compared to the wavelength.
The propagation constant is expressed as where is the attenuation constant that is ignored in our case, since the metamaterial shown in Figure 1 has low loss, which has been verified in our experiments shown later in this paper. The other term β is the phase constant, i.e., wave number.
For a microstrip line, the quasiTEM (transverse electromagnetic wave) mode can be utilized [17], so when ∆z is small as compared to the wavelength, the propagation constant can also be expressed as
Substituting (10), (11), and (12) into (13), the following equation is obtained: where is the total series inductance and is the total shunt capacitance. Other variables in (14) have the same definitions as in (10) and (11). Similarly, the characteristic impedance can be derived:
Under the following designated condition, the phase constant of (14) and the characteristic impedance of (15) are simplified, respectively:where From (17), it is seen that
2.3. Derivation of Phase and Group Velocities
The phase velocity and the group velocity of a TEMmode wave or a sinusoidal wave on a transmission line are [14] Using (17), the phase and group velocities of (24) and (25) are derived:From (26) and (27), and with relations of for , and for , the following properties are obtained:
The above theoretical derivations indicate that the phase velocity of a sinusoidal wave on the metamaterial shown in Figure 1 can be either negative (lefthanded) or positive (righthanded), and it approaches infinity around the transition region of . However, the group velocity is always positive, no matter or . This means that when , the phase and group velocities on the metamaterial are antiparallel; i.e., the two vectors of and point to the opposite directions. However, when , the phase and group velocities on the metamaterial are parallel; i.e., the two vectors of point to the same direction.
2.4. Realization of TimeSpace Singularity
In order to make a detailed analysis around the transition frequency, i.e., , a linear Taylor expansion to the phase constant of (17) leads to the following simplification:The phase velocity of (26) is then expressed asIf the structure shown in Figure 1(a) is cascaded m times, as shown in Figure 1(c), the length of a broadsidecoupled microstrip line isThus, the phase change of a sinusoidal wave propagating along a distance is the following, when (33) is used.Using (34) and (35), the phase delay (or ahead) time isFrom (36) and (37), it is seen that when , both phase change θ and delay time t_{d} become zero, independent of m, i.e., line length. This means that the timespace reduces to a singularity from the sinusoidal wave perspective. That is, time stops and space reduces to a single point, as experienced by the propagating wave of , along any line length described by (35). As common sense, it sounds unlikely that timespace reduces to a singularity here, but this result from the theoretical analysis will later on be verified with simulations and experiments.
The group velocity on the metamaterial described by (31) has a limited value. In order to find the relation between this group velocity and the light velocity, (31) is rewritten using the expression of (19) and (20):Under a condition of , i.e., for righthanded traveling wave, the phase velocity is described as the following [17]: where and are the permittivity and the permeability of freespace, and and are the relative permittivity and permeability in the substrate, whereas and are light velocities in freespace and in a medium, respectively. Comparing (39) to (30), we get the following equation: Substituting (40) into (38), we getIt is seen from (41) that at the transition frequency, i.e., , the group velocity is only half of the light velocity in the substrate medium, even though the phase velocity approaches infinity.
The above derivations indicate that for a TEMmode or sinusoidal wave with a constant frequency, e.g., a single frequency carrier for communication transceivers, the singularity of timespace, i.e., vanishing of timespace, is realized when the phase velocity approaches infinity at . However, a deviation from this sinusoidal wave, e.g., a frequency or phasemodulated wave, smears the singularity, as described by (36) and (37).
3. Simulation
To verify the above analytical results, simulations in both time and frequency domains are done. The simulator used is Advanced Design System (ADS) version 2017 from Keysight. Parameters listed in Table 1 are used for all simulations. The substrate material used is Rogers RO4350B with a dielectric constant = 3.66 ± 0.05 and a loss factor of 0.0031 @ 2.5 GHz (0.0037 @ 10 GHz). Metal layers are Cu with a resistivity of 1.68×10^{−8} Ωm.
3.1. Time Domain
Figure 3 shows an illustration for simulations in the time domain, with a unit cell having a layout component, i.e., m = 1 in (35), generated from a simulation within ADS Momentum. For m = 2 and 3, the layout component is cascaded twice or three times (see Figure 1(c)). The source at the input generates a transient sinusoidal voltage wave. The load is a 50Ω resistive termination.
Figure 4 shows simulation results when the frequency is chosen to be 0.9, 1.0, and 1.2 GHz, respectively. The solid line curve depicts the wave at the input, whereas the dashed line curve depicts the wave at the output. It is seen in Figure 4(a) that the dashed line curve is on the left side of the solid line curve, i.e., the phase delay time t_{d} is negative, indicating (see (37)) , i.e., the phase velocity (see (34)) is negative. In Figure 4(b), the solid and dashed lines overlap with each other, i.e., the phase delay time t_{d} is zero and thus f = f_{0} according to (37). In Figure 4(c), the dashed line curve is on the right side of the solid line curve, i.e., the phase delay time t_{d} is positive, indicating (see (37)) f > f_{0}, i.e., the phase velocity (see (34)) is positive.
(a)
(b)
(c)
Obviously, the results shown in Figures 4(a) and 4(b) are abnormal as compared to those from a conventional microstrip line. On the one hand, one may argue that this can happen if the output curve is one period (360°) after the input curve on a conventional microstrip line. The fact is that the unit cell () used in Figure 3 has a length of 9.80 mm (see Table 1) that is much shorter than the wavelength (300 mm in freespace and 157 mm when = 3.66 is considered) at 1 GHz, so this argument cannot be true. On the other hand, the negative or zero phase delay time agrees very well with the formula shown in (37), when . When f > f_{0}, the phase delay time shown in (37) is positive as shown in Figure 4(c), which is similar to that from a conventional microstrip line.
Figure 5 depicts simulation results when the frequency is chosen to be , and the layout component shown in Figure 3 is cascaded such that m (see (35)) is chosen to be 2, 3, and 4, respectively. It is seen that all output curves overlap with the input curve, similar to that in Figure 4(b). Similar to the explanation for Figure 4(b), the overlap of the input and output curves in Figure 5 cannot be interpreted with repeated wave periods when m = 2, 3, and 4, since the line lengths, i.e., l_{m} = 19.6, 29.4, and 39.2 mm, are still much shorter than a wavelength of 300 mm in freespace or 157 mm in the substrate with a relative permittivity of ε_{r} = 3.66.
Simulation results similar to Figure 4(b) were also observed, when is chosen to be 3 and and 3 GHz, respectively. Obviously, zero phase delay time, i.e., , has been realized at all the designated frequencies of 1, 2, and 3 GHz, which is independent of line length, when it is chosen to be , according to (35).
3.2. Frequency Domain
Due to symmetry of the line structure (see Figure 1), only scatter parameters of S11 and S21 are presented in the paper. Figure 6 depicts simulation results, when the layout component used is the same as that in Figure 4. At 0.9, 1.0, and 1.2 GHz, the S21 phase is +11.53°, 0°, and 14.67°, respectively. Low loss is observed in the S21 and S11 amplitude curves up to 4 GHz in Figure 6(a). Moreover, the transition point, from positive to negative phase when GHz, is clearly seen on the S21 phase curve in Figure 6(b).
(a)
(b)
Figure 7 depicts an S21 diagram using the same structure as that in Figure 6, but m is chosen to be 2 and 3, respectively. Again, low loss up to 4 GHz is observed in Figure 7(a), and it is also seen in Figure 7(b) that the transition point, from positive to negative S21 phase at GHz, remains the same, even though the tangents of the phase curves have changed.
(a)
(b)
Figure 8 depicts an S21 diagram like that in Figure 6, but f_{0} = 2 and 3 GHz, respectively. Similar properties of low loss and zero phase delays are observed, as compared to the case at f_{0 }= 1 GHz.
(a)
(b)
Figure 9 depicts an S21 diagram like that in Figure 6, but m = 10, 20 and 30, respectively. That is, the line length is 1.14 λ, 2.29 λ, and 3.44 λ, respectively, where λ is the wavelength at 1 GHz. It is seen in Figure 9(b) that the zero phase delay remains at f_{0 }= 1 GHz, even though the line length has changed such that it is random with respect to the wavelength.
(a)
(b)
4. Experiment
For verification purposes, experimental samples were designed and fabricated in our printed circuit laboratory at Linköping University. Measurements in the time domain were done in our test and measurement laboratory with an oscilloscope WaveMaster/SDA/DDA 8 ZiB from LeCroy, whereas measurements in the frequency domain were done with a vector network analyzer ZVM 20 GHz from Rohde & Schwarz.
4.1. Sample Design and Fabrication
Samples with the structure shown in Figure 1 and dimension parameters listed in Table 1 have been designed and fabricated. Figure 10 shows a photo of the fabricated samples, at , 2, and 3 GHz, respectively.
4.2. Time Domain Measurement
Figure 11 shows that the measurement result corresponds to Figure 4(b). Other measurement results similar to simulation results shown in Figures 4(a) and 4(c), as well as in Figure 5, are all observed with the samples shown in Figure 10.
4.3. Frequency Domain Measurement
Figure 12 shows an S21 diagram using the same structure as that in Figure 7. It is seen that at the transition point of GHz, zero phase delay is realized, when the section number m = 2 and 3, respectively. Figure 13 depicts Sparameter diagrams like that in Figure 12, but and 3 GHz, respectively. By comparing Figure 12 with Figure 7, and Figure 13 with Figure 8, it is apparent that the simulation and experiment results agree well with each other. The larger loss from experimental samples shown in Figures 12 and 13 as compared to those in Figures 7 and 8 is due to discontinuities from SMA connectors and some misalignment of metal layers when processing experimental samples in our lab.
(a)
(b)
(a)
(b)
5. Results and Discussion
Both the simulation and experimental results have verified that the phase velocity of a sinusoidal wave can approach infinity at a designated angular frequency on a broadsidecoupled transmission line with shortcircuit stubs, as theoretically described by (26) and (29). However, according to (31) the group velocity has still a limited value, i.e., half of the light velocity according to (41), indicating that the group velocity can never exceed the light velocity ( m/s in freespace), even though the phase velocity can approach infinity.
5.1. Realization of Zero Permittivity and Permeability
Equation (29) can be rewritten aswhere and are permittivity and permeability of the metamaterial shown in Figure 1, at a designated frequency of . From (42), we still cannot answer the question as to whether we have realized both zero permittivity and zero permeability, or only one of them. However, from (10) and (11) we havewhere and are two nonzero constants, and and are equivalent inductance and capacitance of the metamaterial, respectively. Inserting (19) into (43) and (44), we obtain Z = 0 and Y = 0, and thus and , respectively. Thus, we can conclude that both zero permittivity and zero permeability are realized when the phase velocity approaches infinity at the designated angular frequency of the metamaterial of any length according to (35).
5.2. Traveling Wave or Standing Wave
One probable doubt can be that the observed timespace singularity at a designated frequency in this study is just standing waves between the input and output ports such that the phase difference is zero between the two ports. However, the simulated S21 curves in Figures 6(a), 7(a), and 8(a) and the measured S21 curves in Figures 12(a) and 13(a) have clearly shown traveling wave properties, no matter how the section numbers or frequencies are changed. Moreover, as shown in Figure 9, when the line length is randomly changed with respect to the wavelength , the zero phase delay remains at the designated frequency f_{0}. All those results also support the analytical results of (34)(37) derived for a traveling wave, not a standing wave between the input and output ports.
5.3. Bandwidth and Application in Microwave Design
This study has started from the first principle of Maxwell’s equations and the transmission line theory. The only two approximations made are, first, the minimum unit cell length that is not infinitesimally small in (10) and (11) and, second, the Taylor expansion of (33). Thus, the analytical results should provide fundamental insights. For instance, the correctness of (34)(37) has been verified with both simulation and experimental results. One of the remaining questions is that regarding the bandwidth of the metamaterial at a designated frequency f_{0} when both permittivity and permeability are near zero. This can be analyzed further with (34), (36), and (37). As seen in Figures 12(b) and 13(b), we have a rather smooth positive to negative phase transition around f_{0}. Therefore, a certain bandwidth can be utilized, in which both permittivity and permeability are near zero, but not exactly zero.
With a nearzero permittivity and permeability realized with the presented low loss metamaterial, microwave designs can be simplified. According to (8) and (9), Kirchhoff’s current and voltage laws are still valid at a microwave or mmwave frequency, when circuitries are designed with the metamaterial around that frequency. As illustrated in Figure 14, no matter how the unit cell shown in Figure 1(a) is connected in series or parallel, the input and output ports merge to a singularity at the designated frequency where the phase velocity approaches infinity. This means that the microwave design methodology in terms of impedance matching, utilization of Smith Chart, and consideration of phase delay is not needed. Instead, the design methodology for analog circuits in terms of voltage loop and current divergence can be used directly. Moreover, some new circuit topology with zero phase delay of interconnects of arbitrary lengths can be worked out.
5.4. Phase, Group, and Propagation Velocities
It is commonly understood that in a homogenous and righthanded material, a harmonic electromagnetic wave propagates with a phase velocity, but a modulated wave from a harmonic electromagnetic wave propagates with a group velocity. However, one must notice that phase, group, and propagation velocities are different concepts; the propagation velocity of a wave can be different from either a phase or a group velocity.
This study has shown that a phase velocity on our presented metamaterial can approach infinity at a designated frequency of f_{0}, but the group velocity is only half of the light velocity in the same substrate of the metamaterial. This result agrees well with that from an electromagnetic field study [8], in which the author shows that a dispersive wave front propagates in the metamaterial, in the beginning, with a group velocity that is half of the propagation, i.e., light velocity in freespace; see Figure of [8]. However, at late times when a steady state is reached the incoming and outcoming waves are “in lockstep with each other”, resulting in an infinitively large propagation velocity, i.e., nearzero delay, through the metamaterial; see Figure of [8].
Within a small bandwidth ∆f around f_{0}, a phase velocity on our presented metamaterial still approaches infinity. Thus, narrowband signals, e.g., a frequencymodulated signal with a bandwidth of ∆f around f_{0} or a phasemodulated signal of f_{0}, can be transferred with nearzero time delay between the input and output ports, according to the principle illustrated in Figure 14, even though the group velocity on the metamaterial before reaching a steady state is only half of the light velocity in the substrate. This is simply due to the fact that the group velocity according to (41) describes the speed of wave dispersion fronts before they reach the output ports shown in Figure 14. After reaching an output port, the input wave of f_{0} propagates through the metamaterial with a phase velocity approaching infinity.
6. Conclusion
Using our invented low loss metamaterial of broadsidecoupled transmission lines with shortcircuit stubs, it is shown that a phase velocity of a traveling sinusoidal wave can approach infinity at various designated frequencies, resulting in zero permittivity and zero permeability. This means that the traveling sinusoidal wave experiences time stop and space shrink to a singularity, independent of its size. This property has been derived from theory, analyzed with simulation, and verified with experimental results.
With this low loss metamaterial, the traditional microwave theory and techniques can be simplified for narrowband signals around various designated frequencies where the phase delay is near zero. That is, the traditional electrical circuit theory based on Kirchhoff’s laws is still valid in the microwave or mmwave spectrum for narrowband signals, utilizing a low loss metamaterial. Consequently, the microwave design methodology in terms of impedance matching, utilization of Smith Chart, and consideration of phase delay is not needed. Instead, the design methodology for analog circuits in terms of voltage loop and current divergence can be used directly. Moreover, some new circuit topology with zero phase delay of interconnects of arbitrary lengths can be worked out.
Data Availability
The simulation and experimental data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors thank Gustav Knutsson at Linköping University for fabrication of experimental samples. This work has been financially supported by the Faculty of Science and Engineering at Linköping University, Sweden. The second author also acknowledges the support by the China Scholarship Council (File No. 201706075057).
References
 S. Gong and M. Karlsson, “Pushing the wireless data rate to the internet speed,” IEEE Access, vol. 4, pp. 8787–8792, 2016. View at: Publisher Site  Google Scholar
 D. R. Smith, W. J. Padilla, D. C. Vier, S. C. NematNasser, and S. Schultz, “Composite medium with simultaneously negative permeability and permittivity,” Physical Review Letters, vol. 84, no. 18, pp. 4184–4187, 2000. View at: Publisher Site  Google Scholar
 R. A. Shelby, D. R. Smith, S. C. NematNasser, and S. Schultz, “Microwave transmission through a twodimensional, isotropic, lefthanded metamaterial,” Applied Physics Letters, vol. 78, no. 4, pp. 489–491, 2001. View at: Publisher Site  Google Scholar
 A. Lai, C. Caloz, and T. Itoh, “Composite right/lefthanded transmission line metamaterials,” IEEE Microwave Magazine, vol. 5, no. 3, pp. 34–50, 2004. View at: Publisher Site  Google Scholar
 V. G. Veselago, “The electrodynamics of substances with simultaneously negative values of ε and μ,” Soviet Physics Uspekhi, vol. 10, no. 4, pp. 509–514, 1968. View at: Publisher Site  Google Scholar
 G. V. Eleftheriades, A. K. Iyer, and P. C. Kremer, “Planar negative refractive index media using periodically LC loaded transmission lines,” IEEE Transactions on Microwave Theory and Techniques, vol. 50, no. 12, pp. 2702–2712, 2002. View at: Publisher Site  Google Scholar
 C. Caloz and T. Itoh, “Transmission line approach of lefthanded (LH) materials and microstrip implementation of an artificial LH transmission line,” IEEE Transactions on Antennas and Propagation, vol. 52, no. 5, pp. 1159–1166, 2004. View at: Publisher Site  Google Scholar
 R. W. Ziolkowski, “Propagation in and scattering from a matched metamaterial having a zero index of refraction,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 70, no. 4, Article ID 046608, 2004. View at: Publisher Site  Google Scholar
 D. Eliecer, X. Cheng, and Y. Yoon, “Highly compact surface micromachined metamaterial circuits using multilayers of lowloss Benzocyclobutene for microwave and millimeter wave applications,” in Proceedings of the 2012 IEEE 62nd Electronic Components and Technology Conference (ECTC), pp. 2062–2069, San Diego, CA, USA, May 2012. View at: Publisher Site  Google Scholar
 A. Rennings, T. Liebig, C. Caloz, and P. Waldow, “MIM CRLH series mode zeroth order resonant antenna (ZORA) implemented in LTCC technology,” in Proceedings of the 2007 AsiaPacific Microwave Conference (APMC '07), pp. 1–4, Bangkok, Thailand, December 2007. View at: Publisher Site  Google Scholar
 Y. Horii, “Development of LTCCbased supercompact multilayered CRLH transmission lines and broadband applications,” in Proceedings of the 2012 AsiaPacific Symposium on Electromagnetic Compatibility (APEMC), pp. 149–152, Singapore, May 2012. View at: Publisher Site  Google Scholar
 A. M. Mahmoud and N. Engheta, “Wavematter interactions in epsilonandmunearzero structures,” Nature Communications, vol. 5, pp. 1–7, 2014. View at: Google Scholar
 D. Schuette, I. Irfanullah, S. Nariyal, S. Khattak, and B. D. Braaten, “Strong coupling (crosstalk) between printed microstrip and complementary split ring resonator (CSRR) loaded transmission lines in multilayer printed circuit boards (pcb),” in Proceedings of the 2016 IEEE International Conference on Electro Information Technology (EIT), pp. 0196–0199, Grand Forks, ND, USA, May 2016. View at: Publisher Site  Google Scholar
 D. M. Pozar, Microwave Engineering, Wiley, 3rd edition, ISBN 0471448788, pp. 69, 151154, 682683, 2005.
 M. Karlsson and S. Gong, “Circular dipole antenna for mode 1 UWB radio with integrated balun utilizing a flexrigid structure,” IEEE Transactions on Antennas Propagation, vol. 57, no. 10, pp. 2967–2971, 2009. View at: Google Scholar
 M. Karlsson, P. Hakansson, and S. Gong, “A frequency triplexer for ultrawideband systems utilizing combined broadside and edgecoupled filters,” IEEE Transactions on Advanced Packaging, vol. 31, no. 4, pp. 794–801, 2008. View at: Publisher Site  Google Scholar
 R. Ludwig and P. Bretchko, RF Circuit Design: Theory and Applications, PrenticeHall, Inc., 2nd edition, 2009.
Copyright
Copyright © 2019 Shaofang Gong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.