## Abstract

In this paper, the dispersive behavior around the plasma frequency in a magnetically uniaxial metamaterial is experimentally investigated. We show by theoretical analysis, parameter retrieval and experiment that when material loss is considered, while the plasma frequency is defined by the frequency where the real part of permeability approaches zero, ultra fast phase velocity actually appears at a slightly lower frequency, due to the change of the dispersion diagram. Both parameter retrieval and experimental data show that within a narrow frequency band to the left of the plasma frequency, the inherent loss keeps finite and is much less than that in the corresponding resonant region. In a real metamaterial sample, an ultra fast phase velocity of 24,440 times the speed of light in free space is measured, and negative phase propagation due to the only negative permeability is observed. The existence of such ultra fast phase velocity with finite loss perfectly explains how the highly directivity antennas based on near-zero refractive index metamaterial work, and can be further used in other applications such as in-phase wave divider and coherent wave sources.

©2010 Optical Society of America

## 1. Introduction

Dispersion describes the electromagnetic (EM) characteristics of a medium at different frequencies, and the phase velocity (PV) of an EM wave is the most important physical quantity characterizing the dispersion [1]. For an EM wave propagating in an isotropic medium is the speed of light in free space. Although PV is always lower than *c* in naturally occurring media that possess refractive indices greater than unity, we also know that it can be faster than *c* without violating any known physical laws, for example, in a plasma medium with a relative permittivity of $1-{\omega}_{p}^{2}/{\omega}^{2}$, where${\omega}_{p}$is known as the plasma frequency (PF). In an extreme and interesting case where the frequency of the EM wave overlaps with the PF, its PV will diverge due to a vanishing refractive index. An infinitely fast PV means that there is no phase delay along the wave propagation direction, which can be used in a number of technologically important applications, such as in-phase power divider, coherent EM wave sources, and highly directivity antennas. However, naturally occurring plasma media, for example ionized gases, are usually difficult to handle in laboratory and adapt to real applications. Fortunately, metamaterials provide a substitution of a man-made plasma medium that can be easily tuned and manipulated in laboratory.

Metamaterials are structured composites possessing unusual electromagnetic properties, such as permittivity and permeability not commonly seen in natural homogeneous materials. Many of these properties have enabled new applications and excited the imaginations of both physicists and engineers in the past years [2, 3]. In this paper, we focus on an anisotropic metamaterial, whose permeability along its optical axis exhibits a plasma-like frequency dependence, at a certain microwave frequency band. We show by theoretical analysis, parameter retrieval and experiment that when material loss is considered, while the PF is defined by the frequency where the real part of permeability approaches zero, ultra fast PV actually appears at a slightly lower frequency, due to the dispersion diagram change. Both parameter retrieval and experimental data show that within a narrow frequency band to the left of the PF, the inherent loss still keeps finite and is much less than that in the corresponding resonant region. In the real metamaterial sample, an ultra fast phase velocity of 24,440 times the speed of light in free space is measured, and negative phase propagation due to the only negative permeability is observed. The existence of such ultra fast phase velocity with finite loss perfectly explains how the highly directivity antennas based on near-zero refractive index metamaterial work [4–9], why the ε-near-zero material can be used for tunneling and squeezing electromagnetic energy through subwavelength narrow channels and waveguides [10], and can be further used in other applications such as in-phase wave divider [11] and coherent wave sources.

## 2. Theoretical analysis

We begin with a magnetically uniaxial effective medium with a scalar permittivity ${\epsilon}_{r}$ and a permeability tensor defined in an (*x, y, z*) Cartesian coordinates of

*y*direction and a wave vector$\overline{k}=\widehat{x}{k}_{x}+\widehat{z}{k}_{z}$, the dispersion relation isFor simplicity and having explicit definition, we consider the PV along the

*x*direction, which means ${k}_{z}=0$. When the loss is neglected (${\gamma}_{e}$ = 0, ${\epsilon}_{r}$ has no imaginary part and ${\mu}_{rz}$ is positive and approaches to zero), ${k}_{x}=\sqrt{{\mu}_{rz}({\epsilon}_{r}{(\omega /c)}^{2})}$ also approaches zero at PF. The PV of such wave, calculated from ${\nu}_{p}=\omega /{k}_{x}$, is theoretically infinite at ${\omega}_{p}$ without transmitting loss. However, when the material loss is taken into account (${\epsilon}_{r}$has a positive imaginary part and ${\gamma}_{e}$ is no longer zero), there would be two complex roots for${k}_{x}$, in which we choose the one with a positive imaginary part to satisfy the energy conservation law. Although its real part still approaches zero at a frequency slightly lower than${\omega}_{p}$, we can still expect to observe ultra-fast PVs along the

*x*direction at frequencies where ${k}_{x}$ approaches zero. To illustrate, Fig. 1 (a) and (b) show the dispersion

*k*-surface of ${k}_{x}$when ${\epsilon}_{r}$ and ${\mu}_{rz}$ are lossless (when ${\gamma}_{e}$ = 0 and ${\epsilon}_{r}$ = 1) and lossy (when ${\gamma}_{e}=0.002{\omega}_{p}$ and ${\epsilon}_{r}=1+0.05\text{i}$), respectively. We see that in lossy case, the real part of ${k}_{x}$ approaches zero at a frequency around 0.98${\omega}_{p}$, no longer exactly at ${\omega}_{p}$, and at the adjacent frequencies lower than 0.98${\omega}_{p}$, the real part of ${k}_{x}$ becomes negative, implying negative phase propagations, albeit with a non-zero imaginary part. Similar dispersion diagram can also be found in Ref [12–14]. Later we will show that in the metamaterial sample we investigated, the loss caused by such imaginary part of ${k}_{x}$ can be actually finite and much smaller than that in the resonant region.

## 3. Parameter retrieval

In the following simulation and experiment, we used a slab-shaped metamaterial sample consisting of multilayered printed circuit boards (PCBs) carrying unit cells made of metallic split-ring resonators (SRRs). Similar metamaterial slabs have been used in Ref [9]. to obtain highly directivity radiation around the PF. This kind of metallic pattern was first proposed in Ref [15]. for obtaining metamaterial showing effective negative permeability. In our realization illustrated in Fig. 2(a)
, the metallic (copper) SRRs are printed on 1-mm-thick FR4 substrates, with a relative permittivity around 4.6 at microwave frequencies. Each board contains 82 unit cells along the *x* direction with a periodicity of 6 mm and 15 unit cells in the *y* direction with a periodicity of 10 mm. Eight layers of such PCBs are aligned along the z direction spaced by 16 mm to form a slab-shaped sample with a width of 500 mm along the *x* direction and a thickness of 105 mm in the *z* direction. In such a slab, for an incidence with z-polarized magnetic field, magnetic resonances can be induced at X-band microwave frequencies inside the metallic rings, yielding a complex dispersion with a negative band of permeability along the *z* direction. Meanwhile, there is no magnetic response in the *x* and *y* directions, allowing${\mu}_{rx}$ and ${\mu}_{ry}$to be regarded as unity. By utilizing the homogenization approach proposed in Ref [16, 17], the dispersion of ${\mu}_{rz}$can be retrieved from the scattering parameters obtained from a full-wave FDTD simulation for the real configuration of the slab, and the result is shown in Fig. 2(b). We see that when the frequency is between 9.3 GHz and 9.5 GHz (marked in gray), ${\mu}_{rz}$exhibits a dispersion very similar to a typical plasma medium with a PF at 9.43 GHz and a finite loss that is much smaller than that in the region centered around 9.2 GHz, which is the resonant frequency of the SRRs.

## 4. Experimental investigation

To experimentally verify the aforementioned analysis and conclusion, we measure the PV within a frequency band near the PF. In the measurement, a small coaxial connector (the inset in Fig. 3(a)
) with a protruded 15-mm-long internal conductor serving as a monopole antenna is placed in the center of the slab, and another identical connector is used as a receiving antenna. To measure the PV inside the metamaterial along the *x* direction (such that ${k}_{z}=0$), three points designated as a, b and c, which are equally spaced by 12 mm on the *x*-axis, are chosen. The phase at each point is measured using an Agilent 8722ES vector network analyzer, and the phase velocity along the *x* direction is calculated by ${v}_{p}=\omega /k=\omega \cdot \Delta d/\Delta \varphi $, where $\Delta d$ is the distance between two adjacent measurement points, i.e., 12 mm, and $\Delta \varphi $ is the corresponding phase difference.

Figure 3(b) shows the phases measured ranging from 9.2 GHz to 9.6 GHz. We see that at the frequencies far away from a frequency around 9.32GHz, the phase difference between the points a, b and c increase. However, at the frequencies close to 9.32GHz, the phase differences become negligible and less frequency dependent. Particularly, at 9.3235 GHz, the phase difference almost vanishes, implying an existence of ultra fast PV. This 9.3235 GHz frequency is 0.9887 times of the retrieved PF at 9.43 GHz in Fig. 2(b), and thus fits quite well with previous analysis. It is also seen that in the adjacent narrow bands lower and higher than the 9.3235 GHz, the phase difference has different signs, which is also in accordance with the previous theoretical analysis, i.e., negative phase propagation could occur at the frequencies to the left of the 9.3235 GHz.

In addition to the phases, Fig. 3(c) shows the amplitude of the electrical field at the same measuring points. As expected, as the frequency is leaving 9.4 GHz (We can assume that the PF of the metamaterial is actually around this frequency) to higher frequencies, there behaves a good transmission due to a real and positive refractive index, while as the frequency is moving to the left, the amplitude is decaying due to the increasing imaginary part of the ${\mu}_{rz}$. The thing interesting is that in the region between 9.3 and 9.4 GHz, the decaying is actually very slow, and the amplitudes measured at points a, b and c are still very much higher than the sensitivity of the instrument, while in the region lower than 9.3 GHz, the decay is drastically accelerated. This trend well matches the variation of the imaginary part of ${\mu}_{rz}$ in Fig. 2(b) between 9.2 and 9.4 GHz, again implying that that in the small band containing the 9.3235 GHz frequency, the loss is finite and therefore the ultra fast PV phenomenon can be used in real, loss-insensitive applications.

In Fig. 4
, the PV curve calculated from the experimental data (in dashed red) is shown, where the *y*-axis represents $sign({\nu}_{p})\times \mathrm{lg}(|{\nu}_{p}|/c)$. We see that at the frequency of 9.3235 GHz, a PV of 24,440*c*, where *c* is the speed of light in free space is experimentally obtained, and with the frequency moves to higher frequencies, the PV drastically falls. This variation trend of the PV versus frequency clearly reflects a plasma-like dispersion. As predicted before, when frequency moves from 9.3235 GHz toward lower frequencies, while the PV also drastically decreases, its direction is negative, meaning that in such lossy area, negative phase propagation, and therefore negative refraction, can be expected to occur without requiring simultaneously negative permittivity and permeability. Although the loss is substantial, the negative refraction at this frequency band is of interest and needs to be further investigated.

For comparison, the PV calculated from the retrieved data is also shown in sold black. Under a sweeping interval of 0.01 MHz (versus the 0.1 MHz interval provided by Agilent 8722ES), the fastest PV of 4,300,000*c* appears at 9.3379 GHz, which is 0.014 GHz higher than the 9.3235 GHz in experimental case. Considering the unavoidable difference of the material and geometry parameters between the simulation and experiment, as well as the difference of 10 times the sampling intervals, the diversity is actually reasonable.

## 5. Conclusion

In short, we discussed the dispersion of a magnetically uniaxial metamaterial and measured experimentally a metamaterial sample consisting of periodically arranged SRRs exhibiting plasma-like dispersion in microwave frequencies. An ultra fast PV of 24,440*c* is observed and the PV can be negative in an adjacent band. In the mean time, the loss is inherent but is also finite. Experimental results fit quite well with theoretically analysis and parameter retrieval. The ultra fast phase velocity along the x direction can be used to perfectly explain how the highly directivity antennas based on such uniaxial metamaterial work, and can be used in other applications such as in-phase wave divider and coherent wave sources. Moreover, at optical frequencies, similar dispersive behaviors can be found in the band diagrams of photonic crystals, indicating that ultra fast PV can also be found in optics, and therefore can be used in optical applications.

## Acknowledgements

The authors sincerely thank Dr. Z. Wang, who is with the Dept. of Physics, MIT, for his valuable suggestions and kind helps in text revisions. This work is sponsored by the ZJNSF (No. Y1080715), 863 Project (No. 2009AA01Z227), NSFC (No. 61071063, 60701007), NCET-07-0750, and the National Key Laboratory Foundation (No. 9140C5304020901).

## References and links

**1. **L. Brillouin, *Wave Propagation and Group Velocity* (Academic, 1960).

**2. **V. G. Veselago, “The electrodynamics of substances with simultaneously negative values of ε and μ,” Sov. Phys. Usp. **10**(4), 509–514 (1968). [CrossRef]

**3. **D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, “Composite medium with simultaneously negative permeability and permittivity,” Phys. Rev. Lett. **84**(18), 4184–4187 (2000). [CrossRef] [PubMed]

**4. **R. W. Ziolkowski, “Propagation in and scattering from a matched metamaterial having a zero index of refraction,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. **70**(4), 046608 (2004). [CrossRef] [PubMed]

**5. **F. Zhang, S. Potet, and J. Caobonell, “Negative-Zero-Positive Refractive Index in a Prism-Like Omega-Type Metamaterial,” IEEE Trans. Microw. Theory Tech. **56**(11), 2566–2573 (2008). [CrossRef]

**6. **A. Alù, M. G. Silveirinha, and N. Engheta, “Transmission-line analysis of epsilon -near-zero-filled narrow channels,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. **78**(1), 016604 (2008). [CrossRef] [PubMed]

**7. **B. Edwards, A. Alù, M. E. Young, M. Silveirinha, and N. Engheta, “Experimental verification of epsilon-near-zero metamaterial coupling and energy squeezing using a microwave waveguide,” Phys. Rev. Lett. **100**(3), 033903 (2008). [CrossRef] [PubMed]

**8. **S. Enoch, G. Tayeb, P. Sabouroux, N. Guérin, and P. Vincent, “A metamaterial for directive emission,” Phys. Rev. Lett. **89**(21), 213902 (2002). [CrossRef] [PubMed]

**9. **T. Jiang, Y. Luo, Z. Wang, L. Peng, J. Huangfu, W. Cui, W. Ma, H. Chen, and L. Ran, “Rainbow-like radiation from an omni-directional source placed in a uniaxial metamaterial slab,” Opt. Express **17**(9), 7068–7073 (2009). [CrossRef] [PubMed]

**10. **M. Silveirinha and N. Engheta, “Tunneling of electromagnetic energy through subwavelength channels and bends using ε-near-zero materials,” Phys. Rev. Lett. **97**(15), 157403 (2006). [CrossRef] [PubMed]

**11. **M. Bozzi, L. Perregrini, D. Deslandes, K. Wu, and G. Conciaurol, “A compact, wideband, phase-equalized waveguide divider/combiner for power amplification,” Microwave Conference, 33rd European (2003).

**12. **G. A. Zheng, “Abrupt change of reflectivity from the strongly anisotropic metamaterial,” Opt. Commun. **281**(8), 1941–1944 (2008). [CrossRef]

**13. **S. Qiao, G. A. Zheng, H. Zhang, and L. X. Ran, “Transition behavior of k-surface: from hyperbola to ellipse,” Prog. Electromagn. Res. **81**, 267–277 (2008). [CrossRef]

**14. **S. Qiao, G. A. Zheng, W. Ren, and L. X. Ran, “Possible abnormal group velocity in the normal dispersive anisotropic media,” J. Electromagn. Waves Appl. **22**(10), 1309–1317 (2008). [CrossRef]

**15. **J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, “Magnetism from conductors and enhanced nonlinear phenomena,” IEEE Trans. Microw. Theory Tech. **47**(11), 2075–2084 (1999). [CrossRef]

**16. **D. R. Smith, S. Schultz, P. Markoš, and C. M. Soukoulis, “Determination of effective permittivity and permeability of metamaterials from reflection and transmission coefficients,” Phys. Rev. B **65**(19), 195104 (2002). [CrossRef]

**17. **X. D. Chen, T. M. Grzegorczyk, B.-I. Wu, J. Pacheco Jr, and J. A. Kong, “Robust method to retrieve the constitutive effective parameters of metamaterials,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. **70**(1), 016608 (2004). [CrossRef] [PubMed]