# Broadside-Coupled Microstrip Lines as Low Loss Metamaterial for Microwave Circuit Design

Authors: Shaofang Gong; Xin Xu; Magnus Karlsson
Journal: Wireless Communications and Mobile Computing
 (2019)

---

## 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 broadside-coupled 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 time-space shrink to a singularity as seen from the propagating wave perspective.


---

## Body

## 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 time-invariant 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 so-called 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 mm-wave 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 mm-wave 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 man-made 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 free-space 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 broadside-coupled 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 right-handed 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 time-space 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 broadside-coupled 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].(1)∇×H^=J^+ε∂E^∂t,(2)∇×E^=-μ∂H^∂t,(3)∇·εE^=ρ,(4)∇·μH^=0,where H^ is the magnetic field intensity, E^ the electric field intensity, ρ the electric charge density, ε permittivity of a medium, μ permeability of a medium, and t the time variable. From (1) and (2),(5)∇×H^=J^,ifε=εrε0=0,(6)∇×E^=0,ifμ=μrμ0=0,where εr is the relative permittivity (dielectric constant) and ε0 the permittivity of free-space, whereas μr is the relative permeability and μ0 the permeability of free-space. Using (5),(7)∇·J^=∇·∇×H^≡0,which is equivalent to the following equation in the integral form according to the divergence theorem [14]. (8)∭v∇·J^dv=∯sJ^·ds^=∑kIk=0atanode,where Ik 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]: (9)∬s∇×E^·ds^=∮cE^·dl^=∑kVk=0alongaclosedloop,where Vk 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,∂E^/∂t and ∂H^/∂t in (1) and (2) can be nonzero, which means that the above derivations are valid under time-variant 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 fast-changing field conditions as in a microwave or a mm-wave circuitry.


### 2.2. Broadside-Coupled Microstrip Line Model

Figures1(a) and 1(b) show oblique-top and side views of our proposed metamaterial structure with two series broadside-coupled microstrip lines and one shunt short-circuit 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 ist1 = 0.17 mm andt2 = 1.52 mm, respectively. The substrate material used is Rogers RO4350B with a dielectric constant εr= 3.66 ± 0.05 and a loss factor of 0.0031 @ 2.5 GHz (0.0037 @ 10 GHz).Table 1
Unit dimensions.
f 0 (GHz) w a (mm) l a (mm) w b (mm) l b (mm) 1 3.67 9.80 0.25 12.3 2 3.67 4.67 0.25 7.05 3 3.67 3.20 0.25 4.46Figure 1
Metamaterial with broadside-coupled microstrip lines and short-circuit stubs with multiple vias. (a) Oblique-top view of a unit cell. (b) Side view of the unit cell. (c) Two or more unit cells are cascaded.
(a) (b) (c)Figure2(a) depicts a lossless equivalent model of the metamaterial unit cell shown in Figures 1(a) and 1(b), whereCp is a series capacitance of the broadside-coupled line segment, either to the left or to the right of the middle short-circuit stab having a shunt inductance ofLs. The shunt capacitance of the broadside-coupled line segment to ground isCpg. 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 right-handed 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 broadside-coupled microstrip lines.Figure 2
Lossless equivalent model of (a) the unit cell shown in Figure1(a), and (b) the middle unit section marked with ∆z shown in Figure 1(c).
(a) (b)Using the model shown in Figure2(b), the series impedanceZ and shunt admittanceY of a unit section of ∆z length are described below. (10)Z=jωLR′△z-1ωCL=jωLR′-1ωCL△z△z,(11)Y=jωCR′△z-1ωLL=jωCR′-1ωLL△z△z,where LR′=2Lp/△z and CR′=2Cpg/△z are unit length series inductance and unit length shunt capacitance of the broadside-coupled line, whereas CL=Cp/2 is the total series capacitance of the broadside-coupled line with a length of∆z. LL = LS is the total shunt inductance of the short-circuit stub, and ω=2πf is an angular frequency, whereasf is a frequency. Note that in (10) and (11), we have used the total inductance and total capacitance forLL andCL, instead of unit length capacitance and unit length inductance. The reason is that, unlike LR′ and CR′,LL andCL 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.,la in Figure 1(a), is dependent on the shunt stub lengthlb. As seen in Table 1, whenlb = 12.3 mm, the unit cell lengthla= ∆z =9.80 mm at 1 GHz that has a wavelengthλ of 300 mm in free-space and 157 mm when εr = 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 (12)γ=α+jβ≈jβ,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 quasi-TEM (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(13)γ=Z∆zY∆z.Substituting (10), (11), and (12) into (13), the following equation is obtained: (14)β=ωLR′-1ωCL∆zωCR′-1ωLL∆z=ω2LRCL-1ω2CRLL-1ω∆zLLCL,where LR=LR′∆z=2Lp is the total series inductance and CR=CR′∆z=2Cpg 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: (15)Z0=Z/∆zY/∆z=LLω2LRCL-1CLω2CRLL-1.Under the following designated condition,(16)LRCL=LLCR,the phase constant of (14) and the characteristic impedance of (15) are simplified, respectively:(17)β=ω2LRCL-1ω∆zLLCL=ω2LLCR-1ω∆zLLCL=ω2/ω02-1ωLω∆z,(18)Z0=LLCL=LRCR,where (19)ω0=1LRCL=1LLCR,(20)ωL=1CLLL.From (17), it is seen that (21)β<0(left-handed),whenω<ω0,(22)β=0,whenω=ω0,(23)β>0(right-handed),whenω>ω0.


### 2.3. Derivation of Phase and Group Velocities

The phase velocityvp and the group velocity vg of a TEM-mode wave or a sinusoidal wave on a transmission line are [14] (24)vp=ωβ,(25)vg=∂β∂ω-1.Using (17), the phase and group velocities of (24) and (25) are derived:(26)vp=∆zω02ω2ω2-ω02ωL,(27)vg=∆zω02ω2ω2+ω02ωL.From (26) and (27), and with relations of ω=ωL=1/LLCL for ω≪ω0, and ω=ωR=1/LRCR for ω≫ω0, the following properties are obtained:(28)vp=-vg=-∆zωL=-∆zLLCL,whenω≪ω0,(29)vp=∝,whenω=ω0,(30)vp=vg=∆zωR=∆zLRCR,whenω≫ω0,(31)vg=∆zω022ωL,whenω=ω0,(32)vg>0,whenω≠0.The above theoretical derivations indicate that the phase velocity of a sinusoidal wave on the metamaterial shown in Figure1 can be either negative (left-handed) or positive (right-handed), and it approaches infinity around the transition region of ω0. However, the group velocity is always positive, no matter ω<ω0 or ω>ω0. This means that when ω<ω0, the phase and group velocities on the metamaterial are antiparallel; i.e., the two vectors of vp and vg point to the opposite directions. However, when ω>ω0, the phase and group velocities on the metamaterial are parallel; i.e., the two vectors of vpandvg point to the same direction.


### 2.4. Realization of Time-Space Singularity

In order to make a detailed analysis around the transition frequency, i.e.,ω≈ω0, a linear Taylor expansion to the phase constant of (17) leads to the following simplification:(33)β=βω0+∂βω∂ωω=ω0∆ω≈2ωL∆zω02ω-ω0.The phase velocity of (26) is then expressed as(34)vp≈∆zω022ωL1-ω0/ω.If the structure shown in Figure 1(a) is cascadedm times, as shown in Figure 1(c), the length of a broadside-coupled microstrip line is(35)lm=m∆z,m=1,2,3…Thus, the phase change of a sinusoidal wave propagating along a distance lm is the following, when (33) is used.(36)θ=βlm=m2ωLω02ω-ω0.Using (34) and (35), the phase delay (or ahead) time is(37)td=lmvp=2mωL1-ω0/ωω02=mfL1-f0/fπf02.From (36) and (37), it is seen that when ω=ω0, both phase changeθ and delay timetd become zero, independent ofm, i.e., line length. This means that the time-space 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 ω0, along any line length described by (35). As common sense, it sounds unlikely that time-space 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):(38)vg=∆z21LRCR.Under a condition of ω≫ω0, i.e., for right-handed traveling wave, the phase velocity is described as the following [17]: (39)vp=1εμ=1εrε0μrμ0=c0εrμr=c0r,where ε0 and μ0 are the permittivity and the permeability of free-space, and εr and μr are the relative permittivity and permeability in the substrate, whereas c0 and c0r are light velocities in free-space and in a medium, respectively. Comparing (39) to (30), we get the following equation: (40)∆zLRCR=c0r.Substituting (40) into (38), we get(41)vg=c0r2.It is seen from (41) that at the transition frequency, i.e., ω=ω0, 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 TEM-mode or sinusoidal wave with a constant frequency, e.g., a single frequency carrier for communication transceivers, the singularity of time-space, i.e., vanishing of time-space, is realized when the phase velocityvp approaches infinity at ω=ω0. However, a deviation from this sinusoidal wave, e.g., a frequency- or phase-modulated wave, smears the singularity, as described by (36) and (37).


## 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].(1)∇×H^=J^+ε∂E^∂t,(2)∇×E^=-μ∂H^∂t,(3)∇·εE^=ρ,(4)∇·μH^=0,where H^ is the magnetic field intensity, E^ the electric field intensity, ρ the electric charge density, ε permittivity of a medium, μ permeability of a medium, and t the time variable. From (1) and (2),(5)∇×H^=J^,ifε=εrε0=0,(6)∇×E^=0,ifμ=μrμ0=0,where εr is the relative permittivity (dielectric constant) and ε0 the permittivity of free-space, whereas μr is the relative permeability and μ0 the permeability of free-space. Using (5),(7)∇·J^=∇·∇×H^≡0,which is equivalent to the following equation in the integral form according to the divergence theorem [14]. (8)∭v∇·J^dv=∯sJ^·ds^=∑kIk=0atanode,where Ik 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]: (9)∬s∇×E^·ds^=∮cE^·dl^=∑kVk=0alongaclosedloop,where Vk 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,∂E^/∂t and ∂H^/∂t in (1) and (2) can be nonzero, which means that the above derivations are valid under time-variant 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 fast-changing field conditions as in a microwave or a mm-wave circuitry.


## 2.2. Broadside-Coupled Microstrip Line Model

Figures1(a) and 1(b) show oblique-top and side views of our proposed metamaterial structure with two series broadside-coupled microstrip lines and one shunt short-circuit 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 ist1 = 0.17 mm andt2 = 1.52 mm, respectively. The substrate material used is Rogers RO4350B with a dielectric constant εr= 3.66 ± 0.05 and a loss factor of 0.0031 @ 2.5 GHz (0.0037 @ 10 GHz).Table 1
Unit dimensions.
f 0 (GHz) w a (mm) l a (mm) w b (mm) l b (mm) 1 3.67 9.80 0.25 12.3 2 3.67 4.67 0.25 7.05 3 3.67 3.20 0.25 4.46Figure 1
Metamaterial with broadside-coupled microstrip lines and short-circuit stubs with multiple vias. (a) Oblique-top view of a unit cell. (b) Side view of the unit cell. (c) Two or more unit cells are cascaded.
(a) (b) (c)Figure2(a) depicts a lossless equivalent model of the metamaterial unit cell shown in Figures 1(a) and 1(b), whereCp is a series capacitance of the broadside-coupled line segment, either to the left or to the right of the middle short-circuit stab having a shunt inductance ofLs. The shunt capacitance of the broadside-coupled line segment to ground isCpg. 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 right-handed 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 broadside-coupled microstrip lines.Figure 2
Lossless equivalent model of (a) the unit cell shown in Figure1(a), and (b) the middle unit section marked with ∆z shown in Figure 1(c).
(a) (b)Using the model shown in Figure2(b), the series impedanceZ and shunt admittanceY of a unit section of ∆z length are described below. (10)Z=jωLR′△z-1ωCL=jωLR′-1ωCL△z△z,(11)Y=jωCR′△z-1ωLL=jωCR′-1ωLL△z△z,where LR′=2Lp/△z and CR′=2Cpg/△z are unit length series inductance and unit length shunt capacitance of the broadside-coupled line, whereas CL=Cp/2 is the total series capacitance of the broadside-coupled line with a length of∆z. LL = LS is the total shunt inductance of the short-circuit stub, and ω=2πf is an angular frequency, whereasf is a frequency. Note that in (10) and (11), we have used the total inductance and total capacitance forLL andCL, instead of unit length capacitance and unit length inductance. The reason is that, unlike LR′ and CR′,LL andCL 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.,la in Figure 1(a), is dependent on the shunt stub lengthlb. As seen in Table 1, whenlb = 12.3 mm, the unit cell lengthla= ∆z =9.80 mm at 1 GHz that has a wavelengthλ of 300 mm in free-space and 157 mm when εr = 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 (12)γ=α+jβ≈jβ,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 quasi-TEM (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(13)γ=Z∆zY∆z.Substituting (10), (11), and (12) into (13), the following equation is obtained: (14)β=ωLR′-1ωCL∆zωCR′-1ωLL∆z=ω2LRCL-1ω2CRLL-1ω∆zLLCL,where LR=LR′∆z=2Lp is the total series inductance and CR=CR′∆z=2Cpg 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: (15)Z0=Z/∆zY/∆z=LLω2LRCL-1CLω2CRLL-1.Under the following designated condition,(16)LRCL=LLCR,the phase constant of (14) and the characteristic impedance of (15) are simplified, respectively:(17)β=ω2LRCL-1ω∆zLLCL=ω2LLCR-1ω∆zLLCL=ω2/ω02-1ωLω∆z,(18)Z0=LLCL=LRCR,where (19)ω0=1LRCL=1LLCR,(20)ωL=1CLLL.From (17), it is seen that (21)β<0(left-handed),whenω<ω0,(22)β=0,whenω=ω0,(23)β>0(right-handed),whenω>ω0.


## 2.3. Derivation of Phase and Group Velocities

The phase velocityvp and the group velocity vg of a TEM-mode wave or a sinusoidal wave on a transmission line are [14] (24)vp=ωβ,(25)vg=∂β∂ω-1.Using (17), the phase and group velocities of (24) and (25) are derived:(26)vp=∆zω02ω2ω2-ω02ωL,(27)vg=∆zω02ω2ω2+ω02ωL.From (26) and (27), and with relations of ω=ωL=1/LLCL for ω≪ω0, and ω=ωR=1/LRCR for ω≫ω0, the following properties are obtained:(28)vp=-vg=-∆zωL=-∆zLLCL,whenω≪ω0,(29)vp=∝,whenω=ω0,(30)vp=vg=∆zωR=∆zLRCR,whenω≫ω0,(31)vg=∆zω022ωL,whenω=ω0,(32)vg>0,whenω≠0.The above theoretical derivations indicate that the phase velocity of a sinusoidal wave on the metamaterial shown in Figure1 can be either negative (left-handed) or positive (right-handed), and it approaches infinity around the transition region of ω0. However, the group velocity is always positive, no matter ω<ω0 or ω>ω0. This means that when ω<ω0, the phase and group velocities on the metamaterial are antiparallel; i.e., the two vectors of vp and vg point to the opposite directions. However, when ω>ω0, the phase and group velocities on the metamaterial are parallel; i.e., the two vectors of vpandvg point to the same direction.


## 2.4. Realization of Time-Space Singularity

In order to make a detailed analysis around the transition frequency, i.e.,ω≈ω0, a linear Taylor expansion to the phase constant of (17) leads to the following simplification:(33)β=βω0+∂βω∂ωω=ω0∆ω≈2ωL∆zω02ω-ω0.The phase velocity of (26) is then expressed as(34)vp≈∆zω022ωL1-ω0/ω.If the structure shown in Figure 1(a) is cascadedm times, as shown in Figure 1(c), the length of a broadside-coupled microstrip line is(35)lm=m∆z,m=1,2,3…Thus, the phase change of a sinusoidal wave propagating along a distance lm is the following, when (33) is used.(36)θ=βlm=m2ωLω02ω-ω0.Using (34) and (35), the phase delay (or ahead) time is(37)td=lmvp=2mωL1-ω0/ωω02=mfL1-f0/fπf02.From (36) and (37), it is seen that when ω=ω0, both phase changeθ and delay timetd become zero, independent ofm, i.e., line length. This means that the time-space 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 ω0, along any line length described by (35). As common sense, it sounds unlikely that time-space 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):(38)vg=∆z21LRCR.Under a condition of ω≫ω0, i.e., for right-handed traveling wave, the phase velocity is described as the following [17]: (39)vp=1εμ=1εrε0μrμ0=c0εrμr=c0r,where ε0 and μ0 are the permittivity and the permeability of free-space, and εr and μr are the relative permittivity and permeability in the substrate, whereas c0 and c0r are light velocities in free-space and in a medium, respectively. Comparing (39) to (30), we get the following equation: (40)∆zLRCR=c0r.Substituting (40) into (38), we get(41)vg=c0r2.It is seen from (41) that at the transition frequency, i.e., ω=ω0, 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 TEM-mode or sinusoidal wave with a constant frequency, e.g., a single frequency carrier for communication transceivers, the singularity of time-space, i.e., vanishing of time-space, is realized when the phase velocityvp approaches infinity at ω=ω0. However, a deviation from this sinusoidal wave, e.g., a frequency- or phase-modulated 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 Table1 are used for all simulations. The substrate material used is Rogers RO4350B with a dielectric constant εr= 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

Figure3 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. Form = 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 3
Illustration of time domain simulations, when a layout component (m = 1) is used. Form = 2 and 3, the layout component is cascaded twice or three times.Figure4 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 timetd is negative, indicating (see (37)) f<f0, 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 timetd is zero and thusf = f0 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 timetd is positive, indicating (see (37))f > f0, i.e., the phase velocity (see (34)) is positive.Figure 4
Time domain simulation when the section numberm = 1, and the frequency of the input signal is chosen to be (a) 0.9 GHz, (b) 1.0 GHz, and (c) 1.2 GHz.
(a) (b) (c)Obviously, the results shown in Figures4(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 (m=1) used in Figure 3 has a length of 9.80 mm (see Table 1) that is much shorter than the wavelength (300 mm in free-space and 157 mm when εr = 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 f≤f0. Whenf > f0, 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.Figure5 depicts simulation results when the frequency is chosen to be f=f0=1GHz, and the layout component shown in Figure 3 is cascaded such thatm (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 whenm = 2, 3, and 4, since the line lengths, i.e.,lm = 19.6, 29.4, and 39.2 mm, are still much shorter than a wavelength of 300 mm in free-space or 157 mm in the substrate with a relative permittivity ofεr = 3.66.Figure 5
Time domain simulation of cascaded sections at 1 GHz, whenm = 2, 3, and 4, respectively.Simulation results similar to Figure4(b) were also observed, when m is chosen to be 3 and f0=2 and 3 GHz, respectively. Obviously, zero phase delay time, i.e., td=0, 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 m∆z, according to (35).


### 3.2. Frequency Domain

Due to symmetry of the line structure (see Figure1), 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 f=f0=1 GHz, is clearly seen on the S21 phase curve in Figure 6(b).Figure 6
Simulation of S-parameters whenm = 1: (a) S21 and S11 magnitudes (dB), and (b) S21 phase response; zero phase delay at 1 GHz is seen.
(a) (b)Figure7 depicts an S21 diagram using the same structure as that in Figure 6, butm 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 f=f0=1 GHz, remains the same, even though the tangents of the phase curves have changed.Figure 7
Simulation of S-parameters,f=f0=1 GHz, andm = 2 and 3, respectively: (a) S11 and S21 magnitudes (dB), and (b) S21 phase response.
(a) (b)Figure8 depicts an S21 diagram like that in Figure 6, butf0= 2 and 3 GHz, respectively. Similar properties of low loss and zero phase delays are observed, as compared to the case atf0= 1 GHz.Figure 8
Simulation of S-parameters,m = 1: (a) S11 and S21 magnitudes (dB), and (b) S21 phase response.
(a) (b)Figure9 depicts an S21 diagram like that in Figure 6, butm = 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 atf0= 1 GHz, even though the line length has changed such that it is random with respect to the wavelength.Figure 9
Simulation of S-parameters,m = 10, 20, and 30, respectively: (a) S11 and S21 magnitudes (dB), and (b) S21 phase response.
(a) (b)


## 3.1. Time Domain

Figure3 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. Form = 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 3
Illustration of time domain simulations, when a layout component (m = 1) is used. Form = 2 and 3, the layout component is cascaded twice or three times.Figure4 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 timetd is negative, indicating (see (37)) f<f0, 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 timetd is zero and thusf = f0 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 timetd is positive, indicating (see (37))f > f0, i.e., the phase velocity (see (34)) is positive.Figure 4
Time domain simulation when the section numberm = 1, and the frequency of the input signal is chosen to be (a) 0.9 GHz, (b) 1.0 GHz, and (c) 1.2 GHz.
(a) (b) (c)Obviously, the results shown in Figures4(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 (m=1) used in Figure 3 has a length of 9.80 mm (see Table 1) that is much shorter than the wavelength (300 mm in free-space and 157 mm when εr = 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 f≤f0. Whenf > f0, 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.Figure5 depicts simulation results when the frequency is chosen to be f=f0=1GHz, and the layout component shown in Figure 3 is cascaded such thatm (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 whenm = 2, 3, and 4, since the line lengths, i.e.,lm = 19.6, 29.4, and 39.2 mm, are still much shorter than a wavelength of 300 mm in free-space or 157 mm in the substrate with a relative permittivity ofεr = 3.66.Figure 5
Time domain simulation of cascaded sections at 1 GHz, whenm = 2, 3, and 4, respectively.Simulation results similar to Figure4(b) were also observed, when m is chosen to be 3 and f0=2 and 3 GHz, respectively. Obviously, zero phase delay time, i.e., td=0, 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 m∆z, according to (35).


## 3.2. Frequency Domain

Due to symmetry of the line structure (see Figure1), 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 f=f0=1 GHz, is clearly seen on the S21 phase curve in Figure 6(b).Figure 6
Simulation of S-parameters whenm = 1: (a) S21 and S11 magnitudes (dB), and (b) S21 phase response; zero phase delay at 1 GHz is seen.
(a) (b)Figure7 depicts an S21 diagram using the same structure as that in Figure 6, butm 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 f=f0=1 GHz, remains the same, even though the tangents of the phase curves have changed.Figure 7
Simulation of S-parameters,f=f0=1 GHz, andm = 2 and 3, respectively: (a) S11 and S21 magnitudes (dB), and (b) S21 phase response.
(a) (b)Figure8 depicts an S21 diagram like that in Figure 6, butf0= 2 and 3 GHz, respectively. Similar properties of low loss and zero phase delays are observed, as compared to the case atf0= 1 GHz.Figure 8
Simulation of S-parameters,m = 1: (a) S11 and S21 magnitudes (dB), and (b) S21 phase response.
(a) (b)Figure9 depicts an S21 diagram like that in Figure 6, butm = 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 atf0= 1 GHz, even though the line length has changed such that it is random with respect to the wavelength.Figure 9
Simulation of S-parameters,m = 10, 20, and 30, respectively: (a) S11 and S21 magnitudes (dB), and (b) S21 phase response.
(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 Zi-B 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 Figure1 and dimension parameters listed in Table 1 have been designed and fabricated. Figure 10 shows a photo of the fabricated samples, at f0=1, 2, and 3 GHz, respectively.Figure 10
Photo of fabricated samples.


### 4.2. Time Domain Measurement

Figure11 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.Figure 11
Time domain measurement whenm = 1, and the frequency of the input signal is chosen to be 1.0 GHz.


### 4.3. Frequency Domain Measurement

Figure12 shows an S21 diagram using the same structure as that in Figure 7. It is seen that at the transition point of f0=1 GHz, zero phase delay is realized, when the section numberm = 2 and 3, respectively. Figure 13 depicts S-parameter diagrams like that in Figure 12, but f0=2 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.Figure 12
Measurement of S-parameters,f0=1 GHz, andm = 2 and 3, respectively: (a) S11 and S21 magnitudes (dB), and (b) S21 phase response.
(a) (b)Figure 13
Measurement of S-parameters,m = 1: (a) S11 and S21 magnitudes (dB), and (b) S21 phase response.
(a) (b)


## 4.1. Sample Design and Fabrication

Samples with the structure shown in Figure1 and dimension parameters listed in Table 1 have been designed and fabricated. Figure 10 shows a photo of the fabricated samples, at f0=1, 2, and 3 GHz, respectively.Figure 10
Photo of fabricated samples.


## 4.2. Time Domain Measurement

Figure11 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.Figure 11
Time domain measurement whenm = 1, and the frequency of the input signal is chosen to be 1.0 GHz.


## 4.3. Frequency Domain Measurement

Figure12 shows an S21 diagram using the same structure as that in Figure 7. It is seen that at the transition point of f0=1 GHz, zero phase delay is realized, when the section numberm = 2 and 3, respectively. Figure 13 depicts S-parameter diagrams like that in Figure 12, but f0=2 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.Figure 12
Measurement of S-parameters,f0=1 GHz, andm = 2 and 3, respectively: (a) S11 and S21 magnitudes (dB), and (b) S21 phase response.
(a) (b)Figure 13
Measurement of S-parameters,m = 1: (a) S11 and S21 magnitudes (dB), and (b) S21 phase response.
(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ω0 on a broadside-coupled transmission line with short-circuit 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 (c0=3×108 m/s in free-space), even though the phase velocity can approach infinity.


### 5.1. Realization of Zero Permittivity and Permeability

Equation (29) can be rewritten as(42)vp=1εω0μω0=∝,whenω=ω0,where εω0and μω0 are permittivity and permeability of the metamaterial shown in Figure 1, at a designated frequency of ω0. 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 have(43)Z=jω0LR′△z-1ω0CL=jω0LR-1ω0CL=jω0m′μω0(44)Y=jω0CR′△z-1ω0LL=jω0CR-1ω0LL=jω0n′εω0where m′ and n′ are two nonzero constants, and m′μω0 and n′εω0are equivalent inductance and capacitance of the metamaterial, respectively. Inserting (19) into (43) and (44), we obtainZ = 0 andY = 0, and thus μω0=0 and εω0=0, 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 ω0 of the metamaterial of any length according to (35).


### 5.2. Traveling Wave or Standing Wave

One probable doubt can be that the observed time-space 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 Figures6(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 frequencyf0. 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∆z 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 frequencyf0 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 aroundf0. Therefore, a certain bandwidth can be utilized, in which both permittivity and permeability are near zero, but not exactly zero.With a near-zero 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 mm-wave 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.Figure 14
Illustration of time-space reduces to a singularity of the metamaterial shown in Figure1(a): (a) the unit cell is connected in both series and parallel; (b) its equivalent circuit where the input and output ports merge to the singularity around the designated frequency where the phase velocity approaches infinity.


### 5.4. Phase, Group, and Propagation Velocities

It is commonly understood that in a homogenous and right-handed 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 off0, 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 free-space; see Figure 4 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., near-zero delay, through the metamaterial; see Figure 5 of [8].Within a small bandwidth∆f aroundf0, a phase velocity on our presented metamaterial still approaches infinity. Thus, narrowband signals, e.g., a frequency-modulated signal with a bandwidth of∆f aroundf0 or a phase-modulated signal off0, can be transferred with near-zero 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 off0 propagates through the metamaterial with a phase velocity approaching infinity.


## 5.1. Realization of Zero Permittivity and Permeability

Equation (29) can be rewritten as(42)vp=1εω0μω0=∝,whenω=ω0,where εω0and μω0 are permittivity and permeability of the metamaterial shown in Figure 1, at a designated frequency of ω0. 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 have(43)Z=jω0LR′△z-1ω0CL=jω0LR-1ω0CL=jω0m′μω0(44)Y=jω0CR′△z-1ω0LL=jω0CR-1ω0LL=jω0n′εω0where m′ and n′ are two nonzero constants, and m′μω0 and n′εω0are equivalent inductance and capacitance of the metamaterial, respectively. Inserting (19) into (43) and (44), we obtainZ = 0 andY = 0, and thus μω0=0 and εω0=0, 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 ω0 of the metamaterial of any length according to (35).


## 5.2. Traveling Wave or Standing Wave

One probable doubt can be that the observed time-space 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 Figures6(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 frequencyf0. 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∆z 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 frequencyf0 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 aroundf0. Therefore, a certain bandwidth can be utilized, in which both permittivity and permeability are near zero, but not exactly zero.With a near-zero 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 mm-wave 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.Figure 14
Illustration of time-space reduces to a singularity of the metamaterial shown in Figure1(a): (a) the unit cell is connected in both series and parallel; (b) its equivalent circuit where the input and output ports merge to the singularity around the designated frequency where the phase velocity approaches infinity.


## 5.4. Phase, Group, and Propagation Velocities

It is commonly understood that in a homogenous and right-handed 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 off0, 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 free-space; see Figure 4 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., near-zero delay, through the metamaterial; see Figure 5 of [8].Within a small bandwidth∆f aroundf0, a phase velocity on our presented metamaterial still approaches infinity. Thus, narrowband signals, e.g., a frequency-modulated signal with a bandwidth of∆f aroundf0 or a phase-modulated signal off0, can be transferred with near-zero 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 off0 propagates through the metamaterial with a phase velocity approaching infinity.


## 6. Conclusion

Using our invented low loss metamaterial of broadside-coupled transmission lines with short-circuit 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 mm-wave 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.



---

*Source: 9249352-2019-07-04.xml*