Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Reflectionless zero refractive index metasurface in the terahertz waveband

Open Access Open Access

Abstract

Zero refractive index materials behave electromagnetically as single points despite the finite dimensions because a propagating electromagnetic wave passes with an infinite phase velocity. However, the composition of naturally occurring materials cannot produce a zero refractive index material because any effect of both of the dielectric and magnetic properties would have to vanish (be near zero). In this report, we demonstrate a zero refractive index metasurface with a refractive index of 0.16 + j0.09, the reflectance of 0.7%, and transmittance of 97.3% at 0.505 THz. The measured relative permittivity and relative permeability are 0.18 − j0.10 and 0.004 + j0.16 at 0.505 THz, respectively. Both the relative permittivity and relative permeability simultaneously approach zero at the same frequency, and the dielectric and magnetic properties appear to be absent (vanish) in the artificial material. The zero refractive index metasurface can offer a material platform for terahertz applications with unprecedented functionalities for 6G (beyond 5G) wireless communications, imaging, and security.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

The on-demand control of light propagation is one of the ultimate goals in optics, and naturally occurring materials with fundamental phenomena such as refraction, reflection, and transmission provide an accessible platform for conventional optical components. The control of high-frequency electric and magnetic fields can artificially provide a design for unprecedented phenomena of electromagnetic waves in a two-dimensional material consisting of smartly engineered sub-wavelength structures [13]. The sub-wavelength structures are termed meta-atoms, and the two-dimensional artificial material is called a metasurface. Artificial materials consisting of sub-wavelength building blocks are metamaterials, and metamaterials have made it possible to produce extraordinary phenomena including negative refraction, a milestone in artificial material development [4]. The control of electromagnetic waves based on a metasurface is called flat optics [5], and conventional bulky optical components can be replaced with optically thin planar components. Planar optical components in the terahertz waveband are desired for cutting-edge applications such as in 6G (beyond 5G) [6] wireless communications [7,8], imaging [9,10], and security [11]. Terahertz continuous-wave (CW) sources with the oscillation of a single frequency with narrow bandwidths at room temperature such as quantum cascade lasers (QCLs) [1216] and resonant-tunneling diodes (RTDs) [17,18] have been powerful and useful tools among terahertz applications. However, most terahertz optical components are bulky compared with terahertz CW sources because there are not many naturally occurring materials for applications to the terahertz waveband [1921]. The development of materials in the terahertz waveband remains a significant problem, and metasurfaces can offer a platform for artificial materials enabling a compact integration of planer terahertz components with terahertz CW sources. The metasurfaces allow arbitrary design of material properties due to the direct control of electric and magnetic fields in the artificial material. A refractive index is described as the root of the product of the relative permittivity and relative permeability ($n = \sqrt {{\varepsilon _\textrm{r}}{\mu _\textrm{r}}} $), and Fresnel reflection is related to the relative wave impedance, the root of the relative permeability over the relative permittivity (${Z_\textrm{r}} = \sqrt {{\mu _\textrm{r}}/{\varepsilon _\textrm{r}}} $).

A wide range of refractive indices from high [2224] to negative [2527] values are presented in the terahertz waveband by the design approach of metamaterials due to the direct control of relative permittivity and relative permeability. Material properties with high and negative refractive indices produce slow and fast wave effects, respectively. Phase manipulation based on metasurfaces is applied to terahertz flat optic as meta lenses [2830]. The control of the relative permeability in naturally occurring materials is commonly difficult to achieve at high frequencies [31]. Conventional metasurfaces face the problem of high Fresnel reflectance because it is not straightforward to design both relative permittivity and relative permeability of similar values. The work in [32] and [33] has recently reported reflectionless terahertz metasurfaces with high and negative refractive indices, respectively, where the relative permeability is close to the relative permittivity enabled by engineered meta-atoms. However, it remains to be elucidated whether the relative permittivity and also the relative permeability simultaneously vanish in an optically thin metasurface. Epsilon-near-zero (ENZ) photonics [34], where the relative permittivity vanishes to approximately zero, is a useful platform for optical components. The work in [3539] has recently reported ENZ structures for the phase manipulation of terahertz waves. Furthermore, ENZ structures can offer attractive wave properties such as stretching, supercoupling, and squeezing for the development of optical components [40]. Unfortunately, Fresnel reflections in ENZ structures are unavoidable because ENZ structures control only the relative permittivity, but not the relative permeability. Near-zero refractive index photonics [41], where both the relative permittivity and relative permeability simultaneously vanish to approximately zero, is a unique platform substituted for ENZ platforms. The work in [42] simulates that zero refractive index artificial materials electromagnetically behave as single points due to an infinite phase velocity and infinite wavelength despite the finite structures. The phase velocity in a material is described as the phase velocity in vacuum over a refractive index. The phase velocity and wavelength are infinite in a material with a refractive index of zero resulting in no spatial phase shift [43]. Measurements in the optical region have confirmed the existence of zero refractive index metamaterials consisting of sub-wavelength building blocks [44,45]. Measurements in the microwave band have also confirmed a zero refractive index metamaterial consisting of three-dimensional structures [46]. To the best of our knowledge, only the research in [47] has shown a zero refractive index metasurface enabled by the simultaneous control of the relative permittivity and permeability in the terahertz waveband achieved by simulations. The work in [48] has shown a periodic array of split-ring resonator holes with a refractive index of zero due to localized waveguide resonance. The work in [49] has shown a nanosphere dispersed liquid crystal in the mid-infrared region.

In this article, we experimentally demonstrate a reflectionless metasurface with a refractive index of zero where both the relative permittivity and the relative permeability vanish to reach zero in the terahertz waveband. The zero refractive index metasurface consists of 10,080 pairs (126 × 80) of meta-atoms, symmetrically aligned paired cut metal wires, on both the front and back of a dielectric substrate. Terahertz time-domain spectroscopy (THz-TDS) confirms a refractive index of 0.16 + j0.09, extremely low reflectance of 0.7%, and high transmittance of 97.3% at 0.505 THz. The measured relative permittivity and permeability are 0.18 − j0.10 and 0.004 + j0.16 at 0.505 THz, respectively. The reflectionless metasurface also has a low refractive index |n| of less than 1.0 and high transmittance characteristics above 80% from 0.49 to 0.545 THz resulting from both the permittivity and permeability exhibiting near-zero properties. The positive and negative signs for the frequency characteristics of the relative permittivity are very similar to those of the relative permeability, and electromagnetic waves propagate without attenuation in the metasurface. At 0.485 THz the metasurface has a refractive index of −1.21 + j0.28 and transmittance of 82.4%, and at 0.57 THz the metasurface has a refractive index of 1.29 + j0.04 and transmittance of 83.0%. Summing up, the metasurface has a transmittance above 80% from 0.485 to 0.57 THz. This zero refractive index metasurface would enable building a platform of artificial materials for planar terahertz components with unprecedented functionalities such as gradient-refractive-index (GRIN) metalenses. The CW sources such as QCLs [1216] and RTDs [17,18] have extremely narrow bandwidths. Further, the laminated structure of the zero refractive index metasurface could be applied to a three-dimensional structure with a refractive index of zero.

2. Zero refractive index metasurface

Figure 1 shows the full model of the zero refractive index metasurface consisting of symmetrically aligned paired cut metal wires on both front and back of a dielectric substrate. The double-sided paired cut metal wires perform as meta-atoms in the metasurface and simultaneously control both the dielectric and magnetic properties. The dielectric property arises from the electric field of the propagating terahertz waves parallel to the cut metal wires. The relative permittivity describing the dielectric properties resonates and changes from negative to positive values around a resonance frequency. The relative permittivity is zero at the state between the negative permittivity and the positive permittivity. The magnetic property is caused by the magnetic field of propagating terahertz waves perpendicular to the cut metal wires on both the front and back of the dielectric substrate. The phenomenon is explained by Faraday’s law, and the double-sided paired cut metal wires perform as a microscale induction coil. The relative permeability describing the magnetic properties resonates and changes from negative to positive values around a resonance frequency. The relative permeability is zero between the negative and positive permeabilities. The relative permittivity of zero and the relative permeability of zero result in a refractive index of zero.

 figure: Fig. 1.

Fig. 1. Zero refractive index metasurface consisting of symmetrically aligned paired cut metal wires on both the front and back of a dielectric substrate.

Download Full Size | PDF

An equivalent circuit adequately explains the dielectric and magnetic properties in the zero refractive index metasurface. Figures 2(a) and (b) show the equivalent circuits of symmetrically aligned paired cut metal wires with the dielectric and magnetic properties, respectively. The dielectric property in the metasurface is mainly determined by the inductance component of a cut wire and the capacitance component at a gap between cut wires along the y-axis as in a series circuit. The series circuit behaves as a resonance circuit with the inductance and capacitance components at a resonant frequency. The resonance of the series circuit represents the resonance of the dielectric property. The magnetic property in the metasurface is mainly determined by the inductance and capacitance components of the cut wires on the front and back substrate as in a parallel circuit. The parallel circuit behaves as a resonance circuit with the inductance and capacitance components at a resonant frequency. The resonance of the parallel circuit expresses the resonance of the magnetic properties.

 figure: Fig. 2.

Fig. 2. Equivalent circuits of the double-sided paired cut metal wires with (a) dielectric and (b) magnetic properties.

Download Full Size | PDF

3. Design of zero refractive index metasurface

Figure 3 shows a single meta-atom extracted from the zero refractive index metasurface in Fig. 1. The one-unit cell model with periodic boundary conditions efficiently designs the properties of the zero refractive index metasurface because the large full model of Fig. 1 is periodic along the x- and y-axes. The simulation of a large full model compared with the wavelength is a very time-consuming process. The control of length l and gap g of the cut metal wires enables the control of both the relative permittivity and relative permeability and the design of a zero refractive index property. Figures 4(a)–(f) show simulated contour maps for the real and imaginary parts of refractive indices, the real parts of the relative permittivity and permeability, reflectance, and transmittance, respectively, with varied gap g and length l of the cut metal wires at 0.50 THz. Other parameters are set as s = 361 µm, w = 120 µm, d = 50 µm, and h = 0.5 µm. The metal of the cut wires is copper, with a conductivity σ = 5.8 × 107 S/m. The dielectric substrate is a cyclo-olefin polymer with a measured refractive index of 1.53 + j0.0012 at 0.5 THz and has low loss properties in the terahertz waveband. The simulations are performed by the frequency domain solver of a finite element method simulator: ANSYS HFSS. As the convergence criteria in HFSS, the maximum number of passes and the maximum Delta S per pass are set to 15 and 0.005, respectively. The simulations use Floquet ports on the incident plane and output plane. The scattering matrices S11 and S21 from the simulations of the one-unit cell model derive effective optical constants such as the refractive index, relative permittivity, and relative permeability, from the following equations [50],

$${n_{\textrm{eff}}} = \frac{{\textrm{Im}[{\ln ({{e^{jn{k_0}({d + 2h} )}}} )} ]+ 2\textrm{m}\pi - j\textrm{Re}[{\ln ({{e^{jn{k_0}({d + 2h} )}}} )} ]}}{{{k_0}({d + 2h} )}}, $$
$${e^{jn{k_0}({d + 2h} )}} = \frac{{{S_{21}}}}{{1 - {S_{11}}\frac{{{Z_\textrm{r}} - 1}}{{{Z_\textrm{r}} + 1}}}}, $$
$${Z_\textrm{r}} ={\pm} \sqrt {\frac{{{{({1 + {S_{11}}} )}^2} - S_{21}^2}}{{{{({1 - {S_{11}}} )}^2} - S_{21}^2}}} , $$
where k0 is the wave number in free space, m is an integer, and Zr is the relative wave impedance divided by the wave impedance in free space. A propagating electromagnetic wave is defined by exp(jk0zjωt), where ω is the angular frequency, and t is time. The thickness of the metasurface is defined as d + 2h in the above derivation. The transmittance is the square of S21, and the reflectance is the square of S11. The black dots in Fig. 4 are the fabricated parameters of the zero refractive index metasurface described in Sec. 4 and have a refractive index of −0.12 + j0.03, reflectance of 0.3%, and high transmittance of 97.0% with l = 202 µm and g = 106 µm. The black dots are the optimized parameters and were chosen to design the reflectionless zero refractive index metasurface. Figure 4(a) shows that the real parts of the refractive indices are changed with a gap g and are slowly changed with a length l. Figure 4(b) shows that the imaginary parts of the refractive indices are changed with a gap g and length l, and the imaginary part of a refractive index increases with longer lengths l. The gray dots with a gap g of 210 µm and length l of 150 µm have a refractive index of 0.007 + j0.33. Figures 4(c) and (d) show the real parts of the relative permittivity and relative permeability, respectively. Dielectric and magnetic properties appear to vanish simultaneously around the parameters with both the relative permittivity and permeability close to 0 in Figs. 4(c) and (d). The black dots show the real part of the relative permittivity of −0.26 and that of the relative permeability of −0.06 at 0.50 THz. The gray dots show that the real parts of the relative permittivity and relative permeability are −0.43 and 0.26, respectively. Figures 4(e) and (f) show the reflectance and transmittance, respectively. The gray dots have a reflectance of 3.2% and a transmittance of 94.2%, respectively. The contour maps effectively derive the optimized parameters of the reflectionless zero refractive index metasurface.

 figure: Fig. 3.

Fig. 3. One-unit cell model of a single meta-atom extracted from a zero refractive index metasurface for the design of dielectric and magnetic properties.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. Contour maps of (a) the real and (b) imaginary parts of refractive indices, real parts of (c) relative permittivity and (d) relative permeability, (e) reflectance, and (f) transmittance at 0.50 THz.

Download Full Size | PDF

Figure 4(a) shows that the refractive indices change from −0.98 to 0.26 with the gap g and length l. The reflectance and transmittance displayed in Figs. 4(e) and (f) are from 0 to 5.0% and from 90 to 100%, respectively, and the material properties in Fig. 4(e) and (f) show low reflectance and high transmittance. Figures 4(a), (e), and (f) show that controlling the gap g and length l enables the design of a reflectionless metasurface with a low refractive index |n| of less than 1.0 and with high transmittance. The material properties with low refractive indices |n| of less than 1.0, reflectionless, and high transmittance characteristics can be directly applied to a wide variety of applications in near-zero refractive index photonics [41]. The control of refractive indices |n| of less than 1.0 with reflectionless and high transmittance characteristics can also fill the gap between high [32] and negative [33] refractive indices and help build a platform for reflectionless artificial materials. Further, Figs. 4(a), (e), and (f) show the only area of meta-atoms where it is possible to control refractive indices when gap g and length l are changed and with the other parameters of s, d, w, and h fixed in the contour maps. The area of meta-atoms would be able to allow the design of a wide range of refractive indices from extremely high to negative values with reflectionless and high transmittance characteristics. A challenge as a yet unresolved problem is to solve the relationship between the area and material properties using formulas. The design concept of material properties with identical meta-atoms developed here could offer an unprecedented platform for reflectionless artificial materials with a full range of refractive indices from extremely high to negative values, including the zero refractive index.

Figures 5(a)–(f) show simulated contour maps for the real and imaginary parts of refractive indices, the real parts of the relative permittivity and permeability, reflectance, and transmittance, respectively, with varied width w and spacing s of the cut metal wires at 0.50 THz. Other parameters are set as g = 106 µm, l = 202 µm, d = 50 µm, and h = 0.5 µm. The design of the reflectionless zero refractive index metasurface also needs the control of the width w and the spacing s. Figure 5(a) shows that the control of the width w and spacing s designs the real parts of the refractive indices from −0.73 to 0.19. Figures 5(e) and (f) show that the reflectance changes from 0% to 50% and the transmittance changes from 50% to 100%, respectively. Figures 5(c) and (d) show that the change of relative permeability is more stable than the change of relative permittivity with varied width w and spacing s. Figure 5(c) shows that the relative permittivity changes from −4.0 to 1.0 with varied width w and spacing s. Figure 5(d) shows that the relative permeability changes from −0.16 to 0.04 with varied width w and spacing s. With w = 100 µm and s = 340 µm, the difference between a relative permittivity of −3.63 and a relative permeability of −0.15 is a local maximal value, resulting in a refractive index of −0.73 + j0.05, reflectance of 43.2%, and transmittance of 52.4% due to the impedance mismatching.

 figure: Fig. 5.

Fig. 5. Contour maps of (a) the real and (b) imaginary parts of refractive indices, real parts of (c) relative permittivity and (d) relative permeability, (e) reflectance, and (f) transmittance at 0.50 THz.

Download Full Size | PDF

4. Measurements and discussion

Figure 6(a) shows a fabricated zero refractive index metasurface consisting of 10,080 pairs (126 × 80) of meta-atoms, symmetrically aligned paired cut metal wires. The fabricated metasurface is a square 4 cm on a side. The wavelength is 600 µm at 0.50 THz. The dimensions of the zero refractive index metasurface are large enough compared with the wavelength and can be simply applied to terahertz flat optics. Figure 6(b) shows a laser microscope image of the fabricated meta-atoms. A low-loss cyclo-olefin polymer film coated with copper layers on both front and back is etched for the fabrication of the zero refractive index metasurface. Additional layers to ensure adhesion are usually used to ensure adhesion between the dielectric substrate and metal, but layers ensuring adhesion (chrome and titanium) will increase conductor loss. The cyclo-olefin polymer is processed using reverse sputtering technique before the sputtering procedure to ensure adhesion between the dielectric substrate and metal. The cyclo-olefin polymer is directly coated with copper without the additional layers to reduce conductor loss, using sputtering technique. Metals with high conductivity such as copper, gold, and silver are acceptable to reduce conductor loss. Copper layers need to coat both the front and back of the thin cyclo-olefin polymer before the wet etching because the double-sided meta-atoms produce the material property of a zero refractive index. The meta-atoms on the front side are fabricated by the wet etching, and those on the back side are also fabricated by the wet etching. The parameters of the fabricated meta-atoms are l = 202 µm, g = 106 µm, s = 361 µm, and w = 120 µm. The thickness of the dielectric substrate is d = 50 µm. The thickness of the copper is h = 0.5 µm, sufficiently thick considering the skin depth of copper of approximately 0.1 µm.

 figure: Fig. 6.

Fig. 6. (a) Fabricated zero refractive index metasurface consisting of symmetrically aligned paired cut metal wires on both the front and back of the cyclo-olefin polymer film. (b) Laser microscope image of meta-atoms.

Download Full Size | PDF

Figures 7(a)–(e) show the frequency characteristics of the optical constants for measurements and simulations. Experiments by THz-TDS TOPTICA with a frequency resolution of 0.005 THz were conducted to measure the transmittance and reflectance of the fabricated metasurface. In the transmittance measurements, pulsed terahertz waves radiated from a transmitting photoconductive antenna pass through the fabricated metasurface and arrive at a receiving photoconductive antenna. The terahertz waves in the transmittance measurements are collimated with a diameter of approximately 18 mm. In reflectance measurements, vertically incident terahertz waves are vertically reflected on the fabricated metasurface and arrive at a receiving photoconductive antenna. The terahertz waves in the reflectance measurements are focused at the fabricated metasurface, and the calculated focusing spot and focal depth are approximately 4.2 mm and 47 mm at 0.50 THz, respectively. Scattering matrices from the measurements by THz-TDS derive the effective optical constants from Eqs. (1), (2), and (3). Deflection of the fabricated metasurface is unavoidable in the measurements and is a cause of the difference in optical lengths between the fabricated metasurface and a reference mirror in the reflectance measurements. The reflectance measurements are compensated for on the condition that the measured phases of the reflectance S11 are the same as the simulated phases. The compensation predicts that the fabricated metasurface is concave with a depth of 112.5 µm along the incident terahertz waves. Figures 7(a) and (b) show that at 0.505 THz the measured refractive index ($n = \sqrt {{\varepsilon _\textrm{r}}{\mu _\textrm{r}}} $), reflectance, and transmittance are 0.16 + j0.09, 0.7%, and 97.3%, respectively. The measurements by THz-TDS demonstrate material properties with a refractive index of zero, reflectionless characteristics, and high transmittance. Figures 7(c) and (d) show the frequency characteristics of the measured relative permittivity and relative permeability, respectively. The measurements agree well with the relative permittivity of 0.18 − j0.10 and relative permeability of 0.004 + j0.16 at 0.505 THz. Both the relative permittivity and the relative permeability approach zero at 0.50 THz, and the dielectric and magnetic properties appear to vanish in the fabricated zero refractive index metasurface. Figure 7(e) shows the real and imaginary parts of the relative wave impedance (${Z_\textrm{r}} = \sqrt {{\mu _\textrm{r}}/{\varepsilon _\textrm{r}}} $). The measurements confirm a relative wave impedance Zr of 0.45 + j0.75 at 0.505 THz. The frequency characteristics of the relative permittivity are similar to those of the relative permeability. The high transmittance, above 80%, is obtained from 0.485 THz to 0.570 THz due to the impedance matching between the fabricated metasurface and free space.

 figure: Fig. 7.

Fig. 7. Measurements and simulations of the frequency characteristics for (a) the refractive index, (b) transmittance and reflectance, (c) relative permittivity, (d) relative permeability, and (e) relative wave impedance.

Download Full Size | PDF

The power absorption in the zero refractive index metasurface is expressed by the following equation [51]

$$Q = \frac{{\omega {\varepsilon _0}{{|E |}^2}}}{{2|{{\mu_\textrm{r}}} |}}[{|{{\mu_\textrm{r}}} |\textrm{Im}({{\varepsilon_\textrm{r}}} )+ |{{\varepsilon_\textrm{r}}} |\textrm{Im}({{\mu_\textrm{r}}} )} ], $$
where ω is the angular frequency, ε0 is the permittivity in vacuum, and E is the electric field. Equation (4) consists of the sum of the dielectric energy loss |µr|Im(εr) and magnetic energy loss |εr|Im(µr). The energy is not amplified in naturally occurring materials, and the imaginary parts of the relative permittivity and relative permeability are positive for a propagating wave expressed as exp(jk0zjωt). Equation (4) shows that the conservation of energy is satisfied in materials when the sum of the dielectric energy loss |µr|Im(εr) and magnetic energy loss |εr|Im(µr) is positive, and the energy is not amplified. Equation (4) should always be positive in the metasurface even when either the imaginary part of the relative permittivity or the imaginary part of the relative permeability is negative [52]. The relative permittivity and relative permeability of the metasurface are calculated effective values [50]. Figure 8(a) shows the measured and simulated frequency characteristics of the dielectric energy loss |µr|Im(εr) and the magnetic energy loss |εr|Im(µr). The dielectric energy loss |µr|Im(εr) in the measurements has negative values, and the magnetic energy loss |εr|Im(µr) in the measurements has positive values. The sum of the dielectric and magnetic losses should be positive to satisfy the conservation of energy, to ensure that the energy is not amplified in the zero refractive index metasurface. Figure 8(b) shows the measured and simulated frequency characteristics of the sum of the dielectric and magnetic energy losses. The sum of the dielectric and magnetic energy losses is almost only positive in the measurements and is everywhere positive in the simulations from 0.45 to 0.55 THz. The errors between the measurements and simulations could be caused by experimental errors and misalignments in optical systems for the transmittance and reflectance measurements. Figure 8(b) shows that the conservation of energy is satisfied in the zero refractive index metasurface. Figure 9 shows the measured and simulated frequency characteristics of the power loss calculated from the measured transmittance and reflectance by THz-TDS. The measured power loss is 2.0% at 0.505 THz. The power loss is positive everywhere in the simulations from 0.45 to 0.55 THz and agrees well with the measurements. Figure 9 shows that the energy is not amplified in the zero refractive index metasurface.

 figure: Fig. 8.

Fig. 8. Measurements and simulations of the frequency characteristics of (a) dielectric energy loss |µr|Im(εr) and magnetic energy loss |εr|Im(µr) and (b) the sum of dielectric energy loss |µr|Im(εr) and magnetic energy loss |εr|Im(µr).

Download Full Size | PDF

 figure: Fig. 9.

Fig. 9. Measurements and simulations of the power loss.

Download Full Size | PDF

5. Summary

We experimentally demonstrate a reflectionless zero refractive index metasurface in the terahertz waveband. The zero refractive index metasurface consists of 10,080 pairs (126 × 80) of double-sided meta-atoms on a dielectric substrate. Measurements by THz-TDS show a refractive index of 0.16 + j0.09, a reflectance of 0.7%, and a transmittance of 97.3% at 0.505 THz. The measurements also yield a relative permittivity of 0.18 − j0.10 and relative permeability of 0.004 + j0.16 at 0.505 THz, and show that the relative permittivity and relative permeability simultaneously approach zero at 0.50 THz. Both dielectric and magnetic properties effectively vanish in the zero refractive index metasurface while the dielectric and magnetic properties cannot vanish in naturally occurring materials. The control of both the relative permittivity and permeability with similar values enable the reflectionless properties, avoiding the high Fresnel reflectance in conventional metasurfaces. The sum of the dielectric and magnetic energy losses is positive in the zero refractive index metasurface, and the conservation of energy is satisfied. The reflectionless zero refractive index metasurface would open the doors to a variety of optically thin planar terahertz components with attractive functionalities such as focusing, collimating, wavefront control, and light vortexes. Further, the path-breaking metasurface helps build a platform of artificial materials for 6G (beyond 5G) wireless communications and terahertz imaging.

Funding

Japan Society for the Promotion of Science (18K04970); Precursory Research for Embryonic Science and Technology (JPMJPR1815); Inamori Foundation; TEPCO Memorial Foundation.

Acknowledgments

The authors wish to thank Tatsuya Kimura, Junichi Yasuda, and Tatsuya Sato in performing the simulations and measurements.

Disclosures

The authors declare no conflicts of interest.

References

1. D. Lin, P. Fan, E. Hasman, and M. L. Brongersma, “Dielectric gradient metasurface optical elements,” Science 345(6194), 298–302 (2014). [CrossRef]  

2. A. V. Kildishev, A. Boltasseva, and V. M. Shalaev, “Planar Photonics with Metasurfaces,” Science 339(6125), 1232009 (2013). [CrossRef]  

3. N. Yu, P. Genevet, M. A. Kats, F. Aieta, J.-P. Tetienne, F. Capasso, and Z. Gaburro, “Light propagation with phase discontinuities: generalized laws of reflection and refraction,” Science 334(6054), 333–337 (2011). [CrossRef]  

4. 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]  

5. N. Yu and F. Capasso, “Flat optics with designer metasurfaces,” Nat. Mater. 13(2), 139–150 (2014). [CrossRef]  

6. P. Yang, Y. Xiao, M. Xiao, and S. Li, “6G Wireless Communications: Vision and Potential Techniques,” IEEE Netw. 33(4), 70–75 (2019). [CrossRef]  

7. T. Harter, S. Ummethala, M. Blaicher, S. Muehlbrandt, S. Wolf, M. Weber, M. M. H. Adib, J. N. Kemal, M. Merboldt, F. Boes, S. Nellen, A. Tessmann, M. Walther, B. Globisch, T. Zwick, W. Freude, S. Randel, and C. Koos, “Wireless THz link with optoelectronic transmitter and receiver,” Optica 6(8), 1063–1070 (2019). [CrossRef]  

8. T. Nagatsuma, G. Ducournau, and C. C. Renaud, “Advances in terahertz communications accelerated by photonics,” Nat. Photonics 10(6), 371–379 (2016). [CrossRef]  

9. H. Guerboukha, K. Nallappan, and M. Skorobogatiy, “Toward real-time terahertz imaging,” Adv. Opt. Photonics 10(4), 843–938 (2018). [CrossRef]  

10. D. M. Mittleman, “Twenty years of terahertz imaging [Invited],” Opt. Express 26(8), 9417–9431 (2018). [CrossRef]  

11. J. Ma, R. Shrestha, J. Adelberg, C.-Y. Yeh, Z. Hossain, E. Knightly, J. M. Jornet, and D. M. Mittleman, “Security and eavesdropping in terahertz wireless links,” Nature 563(7729), 89–93 (2018). [CrossRef]  

12. K. Vijayraghavan, Y. Jiang, M. Jang, A. Jiang, K. Choutagunta, A. Vizbaras, F. Demmerle, G. Boehm, M. C. Amann, and M. A. Belkin, “Broadly tunable terahertz generation in mid-infrared quantum cascade lasers,” Nat. Commun. 4(1), 2021 (2013). [CrossRef]  

13. M. Razeghi, Q. Y. Lu, N. Bandyopadhyay, W. Zhou, D. Heydari, Y. Bai, and S. Slivken, “Quantum cascade lasers: from tool to product,” Opt. Express 23(7), 8462–8475 (2015). [CrossRef]  

14. K. Fujita, M. Hitaka, A. Ito, T. Edamura, M. Yamanishi, S. Jung, and M. A. Belkin, “Terahertz generation in mid-infrared quantum cascade lasers with a dual-upper-state active region,” Appl. Phys. Lett. 106(25), 251104 (2015). [CrossRef]  

15. K. Fujita, M. Hitaka, A. Ito, M. Yamanishi, T. Dougakiuchi, and T. Edamura, “Ultra-broadband room-temperature terahertz quantum cascade laser sources based on difference frequency generation,” Opt. Express 24(15), 16357–16365 (2016). [CrossRef]  

16. K. Fujita, A. Ito, M. Hitaka, T. Dougakiuchi, and T. Edamura, “Low-threshold room-temperature continuous-wave operation of a terahertz difference-frequency quantum cascade laser source,” Appl. Phys. Express 10(8), 082102 (2017). [CrossRef]  

17. H. Kanaya, T. Maekawa, S. Suzuki, and M. Asada, “Structure dependence of oscillation characteristics of resonant-tunneling-diode terahertz oscillators associated with intrinsic and extrinsic delay times,” Jpn. J. Appl. Phys. 54(9), 094103 (2015). [CrossRef]  

18. T. Maekawa, H. Kanaya, S. Suzuki, and M. Asada, “Oscillation up to 1.92 THz in resonant tunneling diode by reduced conduction loss,” Appl. Phys. Express 9(2), 024101 (2016). [CrossRef]  

19. A. Podzorov and G. Gallot, “Low-loss polymers for terahertz applications,” Appl. Opt. 47(18), 3254–3257 (2008). [CrossRef]  

20. M. Wichmann, A. S. Mondol, N. Kocic, S. Lippert, T. Probst, M. Schwerdtfeger, S. Schumann, T. Hochrein, P. Heidemeyer, M. Bastian, G. Bastian, and M. Koch, “Terahertz plastic compound lenses,” Appl. Opt. 52(18), 4186–4191 (2013). [CrossRef]  

21. D. Han, K. Lee, J. Lim, S. S. Hong, Y. K. Kim, and J. Ahn, “Terahertz lens made out of natural stone,” Appl. Opt. 52(36), 8670–8675 (2013). [CrossRef]  

22. M. Choi, S. H. Lee, Y. Kim, S. B. Kang, J. Shin, M. H. Kwak, K.-Y. Kang, Y.-H. Lee, N. Park, and B. Min, “A terahertz metamaterial with unnaturally high refractive index,” Nature 470(7334), 369–373 (2011). [CrossRef]  

23. H. T. Yudistira, A. P. Tenggara, V. D. Nguyen, T. T. Kim, F. D. Prasetyo, C.-G. Choi, M. Choi, and D. Byun, “Fabrication of terahertz metamaterial with high refractive index using high-resolution electrohydrodynamic jet printing,” Appl. Phys. Lett. 103(21), 211106 (2013). [CrossRef]  

24. S. Tan, F. Yan, L. Singh, W. Cao, N. Xu, X. Hu, R. Singh, M. Wang, and W. Zhang, “Terahertz metasurfaces with a high refractive index enhanced by the strong nearest neighbor coupling,” Opt. Express 23(22), 29222–29230 (2015). [CrossRef]  

25. M. Awad, M. Nagel, and H. Kurz, “Negative-index metamaterial with polymer-embedded wire-pair structures at terahertz frequencies,” Opt. Lett. 33(22), 2683–2685 (2008). [CrossRef]  

26. O. Paul, C. Imhof, B. Reinhard, R. Zengerle, and R. Beigang, “Negative index bulk metamaterial at terahertz frequencies,” Opt. Express 16(9), 6736–6744 (2008). [CrossRef]  

27. P. Weis, O. Paul, C. Imhof, R. Beigang, and M. Rahm, “Strongly birefringent metamaterials as negative index terahertz wave plates,” Appl. Phys. Lett. 95(17), 171104 (2009). [CrossRef]  

28. S. Indrišiūnas, H. Richter, I. Grigelionis, V. Janonis, L. Minkevičius, G. Valušis, G. Račiukaitis, T. Hagelschuer, H.-W. Hübers, and I. Kašalynas, “Laser-processed diffractive lenses for the frequency range of 4.7 THz,” Opt. Lett. 44(5), 1210–1213 (2019). [CrossRef]  

29. T. Zhou, J. Du, Y. Liu, and X. Zang, “Helicity multiplexed terahertz multi-foci metalens,” Opt. Lett. 45(2), 463–466 (2020). [CrossRef]  

30. N. Ullah, W. Liu, G. Wang, Z. Wang, A. U. R. Khalid, B. Hu, J. Liu, and Y. Zhang, “Gate-controlled terahertz focusing based on graphene-loaded metasurface,” Opt. Express 28(3), 2789–2798 (2020). [CrossRef]  

31. T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz magnetic response from artificial materials,” Science 303(5663), 1494–1496 (2004). [CrossRef]  

32. K. Ishihara and T. Suzuki, “Metamaterial demonstrates both a high refractive index and extremely low reflection in the 0.3-THz band,” J. Infrared, Millimeter, Terahertz Waves 38(9), 1130–1139 (2017). [CrossRef]  

33. T. Suzuki, M. Sekiya, T. Sato, and Y. Takebayashi, “Negative refractive index metamaterial with high transmission, low reflection, and low loss in the terahertz waveband,” Opt. Express 26(7), 8314–8324 (2018). [CrossRef]  

34. X. Niu, X. Hu, S. Chu, and Q. Gong, “Epsilon-near-zero photonics: A new platform for integrated devices,” Adv. Opt. Mater. 6(10), 1701292 (2018). [CrossRef]  

35. T. Suzuki, T. Sato, M. Sekiya, and J. C. Young, “Epsilon-near-zero three-dimensional metamaterial for manipulation of terahertz beams,” Appl. Opt. 58(11), 3029–3035 (2019). [CrossRef]  

36. A. I. Hernandez-Serrano, R. Mendis, K. S. Reichel, W. Zhang, E. Castro-Camus, and D. M. Mittleman, “Artificial dielectric stepped-refractive-index lens for the terahertz region,” Opt. Express 26(3), 3702–3708 (2018). [CrossRef]  

37. V. Pacheco-Peña, N. Engheta, S. Kuznetsov, A. Gentselev, and M. Beruete, “Experimental realization of an epsilon-near-zero graded-index metalens at terahertz frequencies,” Phys. Rev. Appl. 8(3), 034036 (2017). [CrossRef]  

38. V. P. -Peña, V. Torres, B. Orazbayev, M. Beruete, M. N. -Cía, M. Sorolla, and N. Engheta, “Mechanical 144 GHz beam steering with all-metallic epsilon-near-zero lens antenna,” Appl. Phys. Lett. 105(24), 243503 (2014). [CrossRef]  

39. V. Torres, V. P. Peña, P. R. -Ulibarri, M. N. -Cía, M. Beruete, M. Sorolla, and N. Engheta, “Terahertz epsilon-near-zero graded-index lens,” Opt. Express 21(7), 9156–9166 (2013). [CrossRef]  

40. N. Engheta, “Pursuing near-zero response,” Science 340(6130), 286–287 (2013). [CrossRef]  

41. I. Liberal and N. Engheta, “Near-zero refractive index photonics,” Nat. Photonics 11(3), 149–158 (2017). [CrossRef]  

42. A. M. Mahmoud and N. Engheta, “Wave-matter interactions in epsilon-and-mu-near-zero structures,” Nat. Commun. 5(1), 5638 (2014). [CrossRef]  

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

44. P. Moitra, Y. Yang, Z. Anderson, I. I. Kravchenko, D. P. Briggs, and J. Valentine, “Realization of an all-dielectric zero-index optical metamaterial,” Nat. Photonics 7(10), 791–795 (2013). [CrossRef]  

45. Y. Li, S. Kita, P. Muñoz, O. Reshef, D. I. Vulis, M. Yin, M. Lončar, and E. Mazur, “On-chip zero-index metamaterials,” Nat. Photonics 9(11), 738–742 (2015). [CrossRef]  

46. X. Huang, Y. Lai, Z. H. Hang, H. Zheng, and C. T. Chan, “Dirac cones induced by accidental degeneracy in photonic crystals and zero-refractive-index materials,” Nat. Mater. 10(8), 582–586 (2011). [CrossRef]  

47. Q.-L. Zhang, L.-M. Si, Y. Huang, X. Lv, and W. Zhu, “Low-index-metamaterial for gain enhancement of planar terahertz antenna,” AIP Adv. 4(3), 037103 (2014). [CrossRef]  

48. J.-K. Yang, C. Kang, I. Sohn, and C.-S. Kee, “Effective description of THz localized waveguide resonance through metal film with split ring resonator holes: zero refractive index,” Opt. Express 18(24), 25371–25378 (2010). [CrossRef]  

49. I. C. Khoo, D. H. Werner, X. Liang, A. Diaz, and B. Weiner, “Nanosphere dispersed liquid crystals for tunable negative-zero-positive index of refraction in the optical and terahertz regimes,” Opt. Lett. 31(17), 2592–2594 (2006). [CrossRef]  

50. X. 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 70(1), 016608 (2004). [CrossRef]  

51. Y. Minowa, T. Fujii, M. Nagai, T. Ochiai, K. Sakoda, K. Hirao, and K. Tanaka, “Evaluation of effective electric permittivity and magnetic permeability in metamaterial slabs by terahertz time-domain spectroscopy,” Opt. Express 16(7), 4785–4796 (2008). [CrossRef]  

52. T. Koschny, P. Markoš, D. R. Smith, and C. M. Soukoulis, “Resonant and antiresonant frequency dependence of the effective parameters of metamaterials,” Phys,” Rev. E 68(6), 065602 (2003). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (9)

Fig. 1.
Fig. 1. Zero refractive index metasurface consisting of symmetrically aligned paired cut metal wires on both the front and back of a dielectric substrate.
Fig. 2.
Fig. 2. Equivalent circuits of the double-sided paired cut metal wires with (a) dielectric and (b) magnetic properties.
Fig. 3.
Fig. 3. One-unit cell model of a single meta-atom extracted from a zero refractive index metasurface for the design of dielectric and magnetic properties.
Fig. 4.
Fig. 4. Contour maps of (a) the real and (b) imaginary parts of refractive indices, real parts of (c) relative permittivity and (d) relative permeability, (e) reflectance, and (f) transmittance at 0.50 THz.
Fig. 5.
Fig. 5. Contour maps of (a) the real and (b) imaginary parts of refractive indices, real parts of (c) relative permittivity and (d) relative permeability, (e) reflectance, and (f) transmittance at 0.50 THz.
Fig. 6.
Fig. 6. (a) Fabricated zero refractive index metasurface consisting of symmetrically aligned paired cut metal wires on both the front and back of the cyclo-olefin polymer film. (b) Laser microscope image of meta-atoms.
Fig. 7.
Fig. 7. Measurements and simulations of the frequency characteristics for (a) the refractive index, (b) transmittance and reflectance, (c) relative permittivity, (d) relative permeability, and (e) relative wave impedance.
Fig. 8.
Fig. 8. Measurements and simulations of the frequency characteristics of (a) dielectric energy loss |µr|Im(εr) and magnetic energy loss |εr|Im(µr) and (b) the sum of dielectric energy loss |µr|Im(εr) and magnetic energy loss |εr|Im(µr).
Fig. 9.
Fig. 9. Measurements and simulations of the power loss.

Equations (4)

Equations on this page are rendered with MathJax. Learn more.

n eff = Im [ ln ( e j n k 0 ( d + 2 h ) ) ] + 2 m π j Re [ ln ( e j n k 0 ( d + 2 h ) ) ] k 0 ( d + 2 h ) ,
e j n k 0 ( d + 2 h ) = S 21 1 S 11 Z r 1 Z r + 1 ,
Z r = ± ( 1 + S 11 ) 2 S 21 2 ( 1 S 11 ) 2 S 21 2 ,
Q = ω ε 0 | E | 2 2 | μ r | [ | μ r | Im ( ε r ) + | ε r | Im ( μ r ) ] ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.