## Abstract

In this paper we present the wave propagation method for the retrieving of effective properties of media with circularly polarized eigenwaves, in particularly for chiral metamaterials. The method is applied for thick slabs and provides bulk effective parameters. Its strong sides are the absence of artificial branches of the refractive index and simplicity in implementation. We prove the validity of the method on three case studies of homogeneous magnetized plasma, bi-cross and U-shaped metamaterials.

© 2010 OSA

## 1. Introduction

Chiral metamaterials (MTMs) attract a lot of attention nowadays due to their promising properties. They can be used to achieve giant optical activity [1–3], large or/and broadband circular dichroism [4–6] and even negative index of refraction [7–12]. Effective parameters (EPs) introduction is a common approach to describe the electromagnetic properties of MTMs; therefore it is of crucial importance to have a proper tool for correct EPs restoration.

Usually effective parameters of chiral metamaterials are retrieved with a method based on inversion of reflection and transmission spectra [13–15], so we will refer to it as the standard method (SM). The main problem of the SM is multi-branching and it arises when the real part of the refractive index *n* is retrieved. Different branches of *n* can be very close to each other thus creating difficulties in choosing the correct refractive index. The SM is applicable only to thin MTM slabs consisting of few MTM monolayers as it relies on the transmission simulations. In a thick slab transmission can be at the noise level that distorts the restored parameters severely. At the same time few layers of MTM often do not show bulk effective parameters [16]. In conclusion, the SM is ambiguous since effective parameters and reflection/transmission spectra are not in biunique correspondence.

Recently, we have proposed the wave propagation retrieval method (WPRM) [17] for unambiguous restoration of EPs of bulk MTMs with linearly polarized eigenwaves. In this work we extend the WPRM on the case of media whose eigenwaves are circularly polarized. We therefore show the possibility of using this method for an important class of media, the chiral MTMs. To distinguish this extension from the conventional WPRM we will refer to it as WPRMC. As the problem of ambiguity is not really an issue for the impedance restoration but rather for the refractive indices, in this paper we will present the refractive indices retrieval only. Impedance retrieval is a conventional routine procedure identical to that described in [17]. We use case studies of a homogeneous gyrotropic medium, bi-cross MTM and U-shaped MTM, all having eigenmodes as right-circular polarized (RCP) and left-circular polarized (LCP) waves. These examples serve to illustrate the WPRMC capabilities.

## 2. Wave propagation retrieval methodology

The main idea behind the WPRM is to retrieve the EPs from the cell-averaged electric field E dependence on the coordinate for the wave propagation in a thick (quasi-semi-infinite) metamaterial slab [17]. In this case we can neglect the reflection of the wave from the rear interface and consider only the wave propagating in the forward direction.

Let us consider a linearly *x*-polarized wave propagating along the *z*-axis in vacuum:

*k*

_{0}is the vacuum wavenumber,

*n*

_{R}and

*n*

_{L}are the RCP and LCP refractive indices correspondingly. The linear polarized wave can be decomposed in a sum of RCP and LCP waves as read in formula (1).

The RCP and LCP have equal Fresnel’s transmission coefficients *t*
_{R} = *t*
_{L} = *t*, so the electric field of the wave inside the chiral material is:

These formulas are essential for the WPRMC. In practice, the electric field is highly inhomogeneous inside the MTM unit cell, so we should use the volume averaged field. Bloch impedance *Z*
_{B} (see the relevant discussion in review [18]) is retrieved unambiguously from the single reflection *R* on the vacuum–chiral MTM interface ${Z}_{\text{B}}=(1+R)/(1-R).$ Other EPs can be restored in the following conventional way: refractive index $n=({n}_{\text{R}}+{n}_{\text{L}})/2,$ chirality $\kappa =({n}_{\text{R}}-{n}_{\text{L}})/2.$

Strictly speaking, the complex logarithm in Eqs. (6) and (7) is a multivalued function, so we can find the values of $i\mathrm{ln}({E}_{\text{x}}\pm {E}_{\text{y}})$ with accuracy up to ± 2π*m*, where *m* is an integer. In order to apply a linear fit while using Eqs. (6) and (7), the dependence of $\mathrm{Re}\left[i\mathrm{ln}({E}_{\text{x}}\pm {E}_{\text{y}})\right]$on the unit cells coordinates *z*
_{UC} should be made continuous. If the metamaterial is quasi-homogeneous the difference of the wave phase between neighboring unit cells should be less than π/2, thus aforementioned continuity is easy to persuade.

The constraints of the WPRMC are the same as of the WPRM for linearly polarized waves [17]. First, a metamaterial slab should be thick enough to prevent the reflection from the rear interface and standing wave formation. We found less than 1% error for the WPRM when $-45<n\text{'}{k}_{0}d<45$ and $5<n"{k}_{0}d<56$, where *d* is the slab thickness. For the case of low absorption, one can always increase *d*. Another option is to use a time-domain method and to terminate the simulation when the excitation pulse reaches the rear interface of the slab. Second, the accuracy of the simulations should be sufficient to observe the linear dependence of $i\mathrm{ln}({E}_{x}\pm {E}_{y})$ on the coordinate *z*, for that the electric field doesn’t have to drop down to the noise level within less than three unit cells.

## 3. Case studies

To check the validity of the WPRMC and illustrate its application we consider 3 case studies: 1) homogeneous gyrotropic medium, 2) bi-cross MTM and 3) U-shaped MTM. Commercial software CST Microwave Studio [19] was used in simulations. Simulations were performed in the frequency domain with open boundary conditions in the direction of propagation and periodic boundary conditions in the lateral directions.

#### 3.1 Homogeneous magnetized plasma

The main criteria for the WPRMC employment is that eigenwaves of the material are RCP and LCP waves, so this means that WPRMC is applicable to any material having circularly polarized eigenwaves. As a first example we chose homogeneous plasma magnetized along the direction of the wave propagation *z*. We should note that such medium is not a chiral medium, but its effective refractive index can be calculated analytically, so it is a reliable reference. Its permittivity tensor is $\epsilon =\left[\left({\epsilon}_{1},i{\epsilon}_{2},0\right),\left(-i{\epsilon}_{2},{\epsilon}_{1},0\right),\left(0,0,{\epsilon}_{3}\right)\right]$, where ${\epsilon}_{1}={\epsilon}_{\infty}-\frac{{\omega}_{\text{P}}^{2}(\omega +i{\omega}_{\text{COL}})/\omega}{{(\omega +i{\omega}_{\text{COL}})}^{2}-{\omega}_{\text{C}}^{2}},$
${\epsilon}_{2}=\frac{-{\omega}_{\text{P}}^{2}{\omega}_{\text{C}}/\omega}{{(\omega +i{\omega}_{\text{COL}})}^{2}-{\omega}_{\text{C}}^{2}},$ and ${\epsilon}_{3}={\epsilon}_{\infty}-\frac{{\omega}_{\text{P}}^{2}}{\omega (\omega +i{\omega}_{\text{COL}})}.$ The eigenwaves are RCP and LCP with refractive indices ${n}_{\text{R}}={\left({\epsilon}_{1}+{\epsilon}_{2}\right)}^{1/2}$ and ${n}_{\text{L}}={\left({\epsilon}_{1}-{\epsilon}_{2}\right)}^{1/2}$ [20]. For the simulations we chose plasma frequency ${\omega}_{\text{P}}=\text{1 .0}\times {\text{10}}^{\text{15}}{s}^{-1}$, cyclotron and collision frequencies ${\omega}_{\text{C}}=\text{\hspace{0.17em}}{\omega}_{COL}=\text{3 .8}\times {\text{10}}^{\text{14}}{s}^{-1}$and ε_{∞} = 3. For the simulation the slab of the plasma medium was divided into 100 layers. The WPRMC restores EPs in perfect agreement with the analytical ones (Fig. 1
).

#### 3.2 Bi-cross metamaterial

As the second example we investigated a slightly modified bi-cross metamaterial [21], that is a chiral metamaterial for the optical range. The geometrical parameters are indicated in Fig. 2
. The crosses in the bi-cross are connected in the middle with a cylindrical conducting rod. The bi-cross material is silver, which is approximated as the Drude metal with plasma frequency ${\omega}_{\text{P}}=1.37\times {10}^{16}{\text{s}}^{\text{-1}}$ and collision frequency ${\omega}_{\text{COL}}=8.\text{5}\times {\text{10}}^{\text{13}}{s}^{-1}$ [22]. The silver structure is embedded into silica (refractive index *n* = 1.44).

In this case we lack the analytical results. Therefore WPRMC results are challenged against the ones, obtained by the standard method [15] as a reference. As we have already emphasized in the Introduction the EPs of few MTM layers are not the same as the bulk EPs so we do not anticipate the coincidence of the SM and WPRMC results. We have in mind rather the convergence of the SM results to WPRMC with the increasing thickness of the slab.

Indeed, the effective refractive indices converge very fast with the increase of the number of MTM layers (slab thickness) (Fig. 3 ). The effective properties of 3 layers are similar to the bulk ones restored with WPRMC from simulations of the 36 layers-thick structure.

To further illustrate the restoration procedure we plotted the logarithm of the electric field amplitude dependence on the distance (layer number) inside the MTM slab (Fig. 4 ) for two different cases. At frequency 110 THz (Fig. 4a) the RCP and LCP waves have small absorption, but the MTM exhibits strong circular dichroism. We clearly see the standing wave after the layer number 20, so for retrieving we should use data points before that. At frequency 150 THz (Fig. 4b) absorption of both modes is high (and equal). We observe linear behavior of the field until the 19th layer then the wave amplitude is reaching the noise level. In this case, the last part should be excluded when applying Eqs. (6) and (7) for restoration.

#### 3.3 U-shaped metamaterial

To test our method in the microwave range we chose the recently published U-shaped metamaterial design [23]. The original paper presents the structure consisting of only two layers of MTM, so we stacked multiple layers for the WPRMC application (Fig. 5a ). All the geometrical and material parameters were taken as in [23]. U-shaped resonators in each layer are 90 degrees rotated around their centers with respect to the preceding layer.

As we see from Fig. 5b and 5c, the effective properties change dramatically with the increasing thickness. One bi-layer of the U-shaped MTM shows pronounced negative values of the refractive index real part. For two bi-layers the refractive index is just slightly negative and further it becomes progressively positive for a thick structure (24 bi-layers). We explain this behavior as the result of strong interaction between meta-atoms in the adjacent layers, similar to the one described recently in [24]. Such interaction leads to mutual annihilation of total response finishing with a rather mitigated resonance. It is encouraging that correspondingly the absorption is much smaller than that of a single bi-layer (Fig. 5d and 5e). We suspect that the strong coupling between meta-atoms can represent a problem for the construction of bulk MTM on the basis of most of the proposed chiral MTMs designs. However, this is worth a special investigation that is beyond the scope of this paper.

## 4. Conclusion

In this paper we presented the wave propagation retrieval method for chiral metamaterials. However, the WPRMC is suitable to retrieve bulk effective parameters of any kind of materials with RCP and LCP eigenwaves. The method’s validity is verified on the case studies of homogeneous magnetized plasma and bi-cross chiral MTMs. We also analyzed the case of U-shaped MTMs possessing extremely poor convergence of effective parameters with the slab thickness. In such case the WPRMC provides important information on the bulk effective properties which are difficult or even impossible to obtain by the standard method with the few layers structures. Our method is straightforward in implementation and restores effective properties unambiguously. All this make the WPRMC a useful tool in chiral MTMs research.

## Acknowledgments

The authors gratefully acknowledge partial support from the Danish Research Council for Technology and Production Sciences via the NIMbus project and thank Prof. C. Simovski for pointing on the problems in the effective impedance restoration.

## References and links

**1. **A. V. Rogacheva, V. A. Fedotov, A. S. Schwanecke, and N. I. Zheludev, “Giant gyrotropy due to electromagnetic-field coupling in a bilayered chiral structure,” Phys. Rev. Lett. **97**(17), 177401 (2006). [CrossRef]

**2. **E. Plum, V. A. Fedotov, A. S. Schwanecke, N. I. Zheludev, and Y. Chen, “Giant optical gyrotropy due to electromagnetic coupling,” Appl. Phys. Lett. **90**(22), 223113 (2007). [CrossRef]

**3. **M. Decker, M. Ruther, C. E. Kriegler, J. Zhou, C. M. Soukoulis, S. Linden, and M. Wegener, “Strong optical activity from twisted-cross photonic metamaterials,” Opt. Lett. **34**(16), 2501–2503 (2009). [CrossRef]

**4. **M. Decker, M. W. Klein, M. Wegener, and S. Linden, “Circular dichroism of planar chiral magnetic metamaterials,” Opt. Lett. **32**(7), 856–858 (2007). [CrossRef]

**5. **D. H. Kwon, P. L. Werner, and D. H. Werner, “Optical planar chiral metamaterial designs for strong circular dichroism and polarization rotation,” Opt. Express **16**(16), 11802–11807 (2008). [CrossRef]

**6. **J. K. Gansel, M. Thiel, M. S. Rill, M. Decker, K. Bade, V. Saile, G. von Freymann, S. Linden, and M. Wegener, “Gold helix photonic metamaterial as broadband circular polarizer,” Science **325**(5947), 1513–1515 (2009). [CrossRef]

**7. **S. Tretyakov, I. Nefedov, A. Sihvola, S. Maslovski, and C. Simovski, “Waves and Energy in Chiral Nihility,” J. Electromagn. Waves Appl. **17**(5), 695–706 (2003). [CrossRef]

**8. **J. B. Pendry, “A chiral route to negative refraction,” Science **306**(5700), 1353–1355 (2004). [CrossRef]

**9. **S. Zhang, Y. S. Park, J. Li, X. Lu, W. Zhang, and X. Zhang, “Negative refractive index in chiral metamaterials,” Phys. Rev. Lett. **102**(2), 023901 (2009). [CrossRef]

**10. **E. Plum, J. Zhou, J. Dong, V. A. Fedotov, T. Koschny, C. M. Soukoulis, and N. I. Zheludev, “Metamaterial with negative index due to chirality,” Phys. Rev. B **79**(3), 035407 (2009). [CrossRef]

**11. **B. Wang, J. Zhou, T. Koschny, and C. M. Soukoulis, “Nonplanar chiral metamaterials with negative index,” Appl. Phys. Lett. **94**(15), 151112 (2009). [CrossRef]

**12. **J. Zhou, J. Dong, B. Wang, T. Koschny, M. Kafesaki, and C. Soukoulis, “Negative refractive index due to chirality,” Phys. Rev. B **79**(12), 121104 (2009). [CrossRef]

**13. **C. Menzel, C. Rockstuhl, T. Paul, and F. Lederer, “Retrieving effective parameters for quasiplanar chiral metamaterials,” Appl. Phys. Lett. **93**(23), 233106 (2008). [CrossRef]

**14. **D. H. Kwon, D. H. Werner, A. V. Kildishev, and V. M. Shalaev, “Material parameter retrieval procedure for general bi-isotropic metamaterials and its application to optical chiral negative-index metamaterial design,” Opt. Express **16**(16), 11822–11829 (2008). [CrossRef]

**15. **E. Plum, J. Zhou, J. Dong, V. A. Fedotov, T. Koschny, C. M. Soukoulis, and N. I. Zheludev, “Metamaterial with negative index due to chirality,” Phys. Rev. B **79**(3), 035407 (2009). [CrossRef]

**16. **C. R. Simovski, “Bloch material parameters of magneto-dielectric metamaterials and the concept of Bloch lattices,” Metamaterials (Amst.) **1**(2), 62–80 (2007). [CrossRef]

**17. **A. Andryieuski, R. Malureanu, and A. V. Lavrinenko, “Wave propagation retrieval method for metamaterials: Unambiguous restoration of effective parameters,” Phys. Rev. B **80**(19), 193101 (2009). [CrossRef]

**18. **C. R. Simovski, “Material parameters of metamaterials,” Opt. Spectrosc. **107**(5), 726–753 (2009). [CrossRef]

**19. **“CST Computer Simulation Technology AG,” http://cst.com/

**20. **A. Akhiezer, *Plasma Electrodynamics, Volume One: Linear Theory (Monographs in Natural Philosophy)*, 1st ed. (Pergamon, 1975).

**21. **J. Dong, J. Zhou, T. Koschny, and C. Soukoulis, “Bi-layer cross chiral structure with strong optical activity and negative refractive index,” Opt. Express **17**(16), 14172–14179 (2009). [CrossRef]

**22. **G. Dolling, C. Enkrich, M. Wegener, C. M. Soukoulis, and S. Linden, “Low-loss negative-index metamaterial at telecommunication wavelengths,” Opt. Lett. **31**(12), 1800–1802 (2006). [CrossRef]

**23. **X. Xiong, W. Sun, Y. Bao, M. Wang, R. Peng, C. Sun, X. Lu, J. Shao, Z. Li, and N. Ming, “Construction of a chiral metamaterial with a U-shaped resonator assembly,” Phys. Rev. B **81**(7), 075119 (2010). [CrossRef]

**24. **Z. Li, H. Caglayan, E. Colak, J. Zhou, C. M. Soukoulis, and E. Ozbay, “Coupling effect between two adjacent chiral structure layers,” Opt. Express **18**(6), 5375–5383 (2010). [CrossRef]