## 1.

## Introduction

Approximately 40 years ago, the curiosity of a Russian scientist, Veselago, began the negative index metamaterials (MTMs) process.^{1} Thirty years later (around the 1990s), a British scientist, Sir John Pendry, presented a proposal which soon gained concreteness.^{2}

MTM structures that may be produced by some artificial methods do not exist in nature. The good features of these structures are the negative index of refraction and the fact that a wave and energy in the moving axis have opposite directions. One of the most interesting features of MTMs that attracts the researchers is that these artificially manufactured structures can provide the opportunity to obtain a negative refractive index. The pioneering study in this area was theoretically proposed by Pendry.^{2}

The presence of these materials was demonstrated and was experimentally manufactured by Shelby et al.^{3} Their work consists of a copper wire and a split-ring resonator system. Dielectric photonic crystals were then investigated by Ozbay et al. for the first time.^{4}

Application of MTMs in the terahertz (THz) range has recently gained great importance.^{5}6.7.8.9.10.^{–}^{11} THz MTMs are the new class of artificial composite materials under development that work in the THz frequency regime. Generally, the THz frequency range used in the research of materials is defined in the frequency spectrum from 0.1 to 10 THz.

MTMs are a promising new concept to be used in the design and development of many new technological areas such as electronics, microwave, optical circuits, antenna miniaturization, antenna design, and so on. Working in the THz frequency range leads us to develop useful devices to be used in security screening, medical imaging, wireless communication systems, nondestructive evaluation, and development of the applications of chemical identification.

Chirality is a very important concept of MTMs for which a negative refractive index can be obtained around the resonance point through both the right- and left-hand circularly polarized waves for different transmission values.^{12}

In this study, we theoretically and experimentally examined and analyzed the asymmetric transmission phenomenon for linearly polarized electromagnetic (EM) waves in both microwave and THz frequency regimes. The asymmetric characteristics of the chiral MTMs are investigated and evaluated in detail. The proposed dynamic structures that provide asymmetric transmission can easily be scaled down for any desired frequency range in order to fit any desired application.

## 2.

## Theoretical Analysis

Chirality, one of the most interesting features that MTMs have, can be used to easily obtain a negative refractive index for various transmission values around any desired resonance point compared to other procedures for both right and left circularly polarized waves. Asymmetric transmission properties allow a certain difference between the polarization states when a wave is applied to opposite sides. In this study, for theoretical analysis, an incoming plane wave propagating in the ($+z$) direction with a time dependence ${e}^{-iwt}$ of is considered^{12}13.14.15.^{–}^{16}

This medium is assumed to be a dielectric medium that is symmetrically embedded (e.g., vacuum). The $T$-matrix relates the complex amplitudes of the incident field to that of the transmitted field

## (3)

$$\left(\begin{array}{c}{T}_{x}\\ {T}_{y}\end{array}\right)=\left(\begin{array}{cc}{T}_{xx}& {T}_{xy}\\ {T}_{yx}& {T}_{yy}\end{array}\right)\left(\begin{array}{c}{E}_{x}\\ {E}_{y}\end{array}\right)={\widehat{T}}_{\mathrm{lin}}^{f}\left(\begin{array}{c}{E}_{x}\\ {E}_{y}\end{array}\right).$$^{12}13.

^{–}

^{14}

## (4)

$${T}_{\mathrm{circ}}^{f}=\left(\begin{array}{cc}{T}_{++}& {T}_{+-}\\ {T}_{-+}& {T}_{--}\end{array}\right)\phantom{\rule[-0.0ex]{2em}{0.0ex}}\text{and}\phantom{\rule[-0.0ex]{2em}{0.0ex}}{T}_{\mathrm{circ}}^{b}=\left(\begin{array}{cc}{T}_{++}& {T}_{-+}\\ {T}_{+-}& {T}_{--}\end{array}\right).$$## (5)

$${T}_{\mathrm{circ}}^{f}=\left(\begin{array}{cc}{T}_{++}& {T}_{+-}\\ {T}_{-+}& {T}_{--}\end{array}\right)=\frac{1}{2}\left[\begin{array}{cc}{T}_{xx}+{T}_{yy}+i({T}_{xy}-{T}_{yx})& {T}_{xx}-{T}_{yy}+i({T}_{xy}+{T}_{yx})\\ {T}_{xx}-{T}_{yy}-i({T}_{xy}+{T}_{yx})& {T}_{xx}+{T}_{yy}-i({T}_{xy}-{T}_{yx})\end{array}\right].$$The linearly and circularly polarized waves of the asymmetric transmission are commonly characterized by the parameter $\mathrm{\Delta}$. This parameter shows the difference between transmittances in the ($+z$) and ($-z$) directions. The asymmetric transmission parameter $\mathrm{\Delta}$ is defined as

## (6)

$${\mathrm{\Delta}}_{\mathrm{lin}}^{(x)}={|{T}_{yx}|}^{2}-{|{T}_{xy}|}^{2}=-{\mathrm{\Delta}}_{\mathrm{lin}}^{(y)},$$## (7)

$${\mathrm{\Delta}}_{\mathrm{cir}}^{(+)}={|{T}_{-+}|}^{2}-{|{T}_{+-}|}^{2}=-{\mathrm{\Delta}}_{\mathrm{circ}}^{(-)},$$^{12}13.14.15.

^{–}

^{16}

However, some studies have shown that the mirror symmetry is getting lost in the propagation direction for hybridizing MTM structures.^{5}^{,}^{17} In addition, the reason for the chirality can be explained by the fact that MTM structures that have a polarization conversion feature that the polarization state of the wave will change when it propagates through the structure. This case is called as optical activity of chiral structures. Depending on the handedness of a material, the azimuth of the polarization plane rotates clockwise or counter-clockwise when it passes through a chiral medium. This polarization conversion feature may lead to design polarization sensitive devices in any desired frequency range.

## 3.

## Numerical and Experimental Setup and Design

We numerically and experimentally investigated the asymmetric transmission properties of the proposed chiral MTM structures. The designed structures have resonator-metallic parts placed on the front and the back side of the dielectric substrate in order to provide asymmetric transmission for linearly polarized EM waves in both microwave and THz frequency regions. Furthermore, the proposed structures can also provide asymmetry by changing the parameters of the structure such as the angle, gap, and width of the metallic strips.

For a microwave frequency regime, FR4 is chosen as the dielectric substrate and the metallic patterns are modeled as copper sheets with an electrical conductivity of $5.8001\times {10}^{7}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{S}/\mathrm{m}$ and a thickness of 0.036 mm. The thickness, loss tangent, and relative permittivity of FR4 are 1.6 mm, 0.02, and 4.2, respectively. Figure 1 shows the schematic diagram of the unit cell. The unit cell dimensions of the proposed resonator are shown in Table 1.

## Table 1

Unit cell dimensions of the proposed resonator (Fig. 1).

Parameters | Microwave frequency regime | Terahertz frequency regime |
---|---|---|

l1 | 40 mm | 900 μm |

l2 | 24 mm | 600 μm |

l3 | 20 mm | 500 μm |

g | 2 mm | 40 μm |

d | 1.6 mm | 450 μm |

Angle | 15 deg | 15 deg |

On the other hand, for a THz frequency regime, the proposed resonators on the top and bottom layers are separated by a Quartz (fused) dielectric substrate which is selected as the substrate with the thickness of 450 *μ*m. The loss tangent and relative permittivity of the Quartz are 0.0004 and 3.75, respectively. The metallic structures on the top and bottom layers of the substrate are chosen as silver sheets with an electrical conductivity of $6.3\times {10}^{7}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{S}/\mathrm{m}$ and a thickness of 0.2 *μ*m. The reason for this is that silver is a soft, white, lustrous transition metal, and possesses extremely low resistivity. The unit cell dimensions of the proposed resonator for this frequency range are also shown in Table 1.

A commercial full wave EM solver (CST Microwave Studio, Darmstadt, Germany) based on a finite integration technique has been used for numerical studies. Unit cell and open add space boundary conditions are used in both gigahertz (GHz) and THz regions. For determining the reflection and transmission properties of the proposed structure, the numerical simulations have been implemented in both GHz and THz ranges. Transmission matrices of the proposed chiral MTMs are calculated from Eq. (8). ${T}_{xx}$ and ${T}_{yy}$ represent copolarized transmissions, and ${T}_{xy}$ and ${T}_{yx}$ define cross-polarized transmissions. These values are obtained by ${T}_{ij}={E}_{j}^{\text{out}}/{E}_{i}^{\text{in}}$, relating the incident (${E}_{i}^{\text{in}}$) and transmitted (${E}_{j}^{\text{out}}$) waves, for example, ${T}_{xy}={E}_{y}^{\text{out}}/{E}_{x}^{\text{in}}$, where ${E}_{x}^{\text{in}}$ is the incident $x$-polarized electric field and ${E}_{y}^{\text{out}}$ is the transmitted $y$-polarized electric field. The transmitted wave of the electric field can be written as ${E}^{\text{out}}(r,t)=|{E}_{x}^{\text{out}}|\mathrm{cos}\text{\hspace{0.17em}}(kz-\omega t)$, $|{E}_{y}^{\text{out}}|\mathrm{cos}\text{\hspace{0.17em}}(kz-\omega t+\varphi )$ for both copolarized and cross-polarized transmissions. Transmission matrices can be found by using the following equation:^{18}19.^{–}^{20}

## (8)

$$\left(\begin{array}{cc}{T}_{++}& {T}_{+-}\\ {T}_{-+}& {T}_{--}\end{array}\right)\phantom{\rule{0ex}{0ex}}=\frac{1}{2}\left[\begin{array}{cc}{T}_{xx}+{T}_{yy}+i({T}_{xy}-{T}_{yx})& {T}_{xx}-{T}_{yy}+i({T}_{xy}+{T}_{yx})\\ {T}_{xx}-{T}_{yy}-i({T}_{xy}+{T}_{yx})& {T}_{xx}+{T}_{yy}-i({T}_{xy}-{T}_{yx})\end{array}\right],$$## 4.

## Numerical and Experimental Results

For investigating the frequency responses of the proposed chiral MTMs, a CWS structure has been designed and fabricated in the microwave frequency regime with respect to the plane waves propagating along the ($+z$) and ($-z$) directions. Figures 2 (numerical and experimental results in GHz) and 3 (numerical results in THz) show the results of the four transmission matrix elements (${T}_{xx}$, ${T}_{xy}$, ${T}_{yx}$, ${T}_{yy}$) of the slab for propagations in the ($+z$) and ($-z$) directions, respectively. One can see from these figures that the numerical results verify the measurement results.

The resulting numerical and experimental magnitudes of ${T}_{xx}$, ${T}_{xy}$, ${T}_{yx}$, and ${T}_{yy}$ are shown in Figs. 2 and 3, in order. The obtained results show asymmetric transmission for linear polarization in the resonance frequencies which may be acquired when the propagation direction is the same. This inequality confirms the presence of a circular dichroism ${T}_{+-}\ne {T}_{-+}$.

For the microwave frequency regime, the copolarization parameter of the transmission coefficient ${T}_{xx}$ reaches its maximum value of 0.86 at the resonance frequency of $f=5.99\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{GHz}$. With this regard, the simulation and experimental results show that ${T}_{xx}$ and ${T}_{yy}$ are almost equal to each other. Cross-polarization parameters of the transmissions ${T}_{xy}$ and ${T}_{yx}$ reach a maximum value of approximately 0.38 at the resonance frequency of $f=4.43\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{GHz}$.

For the THz frequency regime, the copolarization parameter of transmission coefficient ${T}_{xx}$ reaches its maximum value of around 0.86 at the resonance frequency of $f=0.595\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{THz}$. In this regard, the simulation and experimental results show that ${T}_{xx}$ and ${T}_{yy}$ are almost equal. In addition, cross-polarization parameters of the transmissions ${T}_{xy}$ and ${T}_{yx}$ reach to a maximum value of approximately 0.41 at the resonance frequency of $f=0.583\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{THz}$.

For both linear and circular polarizations, the asymmetric transmission parameter, $\mathrm{\Delta}$, is calculated and normalized using theoretical analysis via Eqs. (6) and (7), as shown in Fig. 4. It can be seen that the suggested model exhibits asymmetric transmission at two different resonance frequencies, 0.583 THz and 0.587 THz, for linear polarization. No significant change is observed in ($\mathrm{\Delta}\approx 0$) for the circular polarization.

For the next investigation, we numerically retrieved the circular dichroism, theta degree, and chirality values of the proposed chiral MTM with constitutive relations. First, the circular dichroism is retrieved and realized, which can be defined as the differences between RCP and LCP waves of two polarizations propagating through the medium as shown in Fig. 5. The circular dichroism is theoretically characterized by the ellipticity^{18}19.^{–}^{20}

## (9)

$$\eta =\frac{1}{2}\text{\hspace{0.17em}}{\mathrm{sin}}^{-1}\left(\frac{{|{T}_{++}|}^{2}-{|{T}_{--}|}^{2}}{{|{T}_{++}|}^{2}+{|{T}_{--}|}^{2}}\right).$$Second, optical activity is retrieved and realized which could rotate the polarization plane of a linearly polarized wave propagating through, as shown in Fig. 5. A linearly polarized plane wave will rotate when it passes through a chiral medium. This rotation is called optical activity and is characterized by the polarization azimuth rotation angle of an elliptically polarized wave. The azimuth rotation is calculated by the copolar and crosspolar transmission data^{18}19.^{–}^{20}

We realized a retrieval procedure to calculate the chirality parameters of chiral MTMs by using Eq. (9). Simulated chirality values calculated from copolar and crosspolar transmission data for the proposed structure between 0.58 and 0.6 THz are shown in Fig. 6 (Refs. 1819.–20).

The proposed structure does not show chiral MTM characteristics between 0.592 THz and 0.594 THz. The chirality value is extremely low and near zero, but the structure shows a very large chirality value at the resonance frequencies of 0.583 THz and 0.587 THz.

Furthermore, the effects of the changes in the dimensions of the proposed resonator on the chirality are realized and evaluated. Obtained numerical results are shown in Fig. 7. The proposed geometry consists of a square-shaped resonator with gaps in the unit cell on one side and the rotated (15 deg) version of the same geometry on the other side. We have selected this shape because it is very simple, flexible, and easy to manufacture, which will allow us to adjust the structure easily by simply tuning the geometrical dimensions in order to work in any desired frequency regime. This property provides a shiftable frequency band and large chirality. The length of the wire ($l2$) is increased from 595 um to 650 um and the resonator at the back is rotated from 15 deg to 30 deg for parametric examination. It can be seen from the parametric study that the suggested structure has reconfigurable chirality values. In addition, the best value for the angle of the back resonator is found to be 15 deg. Therefore, this value is used in the design for the asymmetric transmission calculations.

Moreover, up- and down-shifts in the frequency response of the chirality value can be achieved when the dimensions of the resonator vary in any desired range. This feature provides mechanical tunability and flexibility to obtain any desired optical activity. Furthermore, this design possesses a very strong and dynamic optical activity due to the features previously mentioned.

Figures 8 and 9 show the electric field and surface current distributions which verify the characteristics of the resonances of the proposed design. The electric field and surface current distributions are obtained and evaluated separately for the proposed chiral MTM at the resonance frequencies of 0.583 THz and 0.587 THz. It can be seen from the surface current distributions shown in Fig. 8 that there are parallel and antiparallel currents on the top and the bottom layers of the resonator. This means that the asymmetric resonance is directly related to the antiparallel currents excited by the magnetic dipoles and the symmetric resonance mode is realized by the parallel currents excited by the electric field (Fig. 9). Moreover, the physical mechanism of the resonance mode for 0.583 THz is different than the one for 0.587 THz.

## 5.

## Conclusion

In conclusion, the proposed chiral MTM consisting of a dielectric substrate and a square-shaped resonator with gaps are analyzed theoretically, experimentally, and numerically for both microwave and THz frequency regimes. In addition, obtained numerical results are confirmed with our experimental results in GHz and with the simulated results in THz with very good agreement. For designing novel devices in the THz frequency band, different frequency regimes can be studied based on the proposed GHz study.

Numerical simulations showed that it is possible to adjust the optical activity, for instance by actively controlling the angle of polarization of light via integrating semiconductor layers into the chiral MTM design. Moreover, a tunability feature of the proposed model provides an additional degree of freedom for polarization control. In other words, the tunable chiral MTM design enables applications such as antireflection coatings, new THz polarization rotators, polarization control devices, security screening, medical imaging, nondestructive evaluation and the development of chemical identification applications, etc., in the THz-wave frequency range.

## Acknowledgments

This work is supported by the Scientific and Technological Research Council of Turkey (TUBITAK) under Project No. 113E290. One of the authors M. K. also acknowledges partial support from the Turkish Academy of Sciences. Moreover, the help of Professor C. Sabah is greatly appreciated with respect to the installation of the experimental setup and authors would like to thank the editors and anonymous reviewers for their suggestions to improve the paper.

## References

## Biography

**Zafer Ozer** received his PhD degree from the Physics Department of the University of Cukurova, Adana, Turkey. He is now working in Mustafa Kemal University, Hatay, Turkey. His research interests are photonics, EM waves, metamaterials, and waveguides.

**Furkan Dincer** received the BSc and MSc degrees in electrical and electronics engineering from Sutcu Imam and Yuzuncu Yil Universities, Turkey. He is now a PhD student in Mustafa Kemal University, Hatay, Turkey. His research interests are related to functional microwave structures and metamaterials.

**Muharrem Karaaslan** received his PhD degree from the Physics Department of the University of Cukurova, Adana, Turkey, in 2009. He is the coauthor of about 20 scientific contributions published in journals and conference proceedings. His research interests are applications of metamaterials, analysis and synthesis of antennas, and waveguides.

**Oguzhan Akgol** received the BSc, MSc, and PhD degrees in electrical and electronics engineering from Inonu University, Turkey; Polytechnic University, Brooklyn, NY, USA; and University of Illinois at Chicago, (UIC), Chicago, IL, USA, respectively. He is now working in Mustafa Kemal University, Hatay, Turkey. His research interests are EM scattering, antennas, and DNG materials.