## 1. Introduction

The electromagnetic characterization of materials is a fundamental problem in many research areas of Electromagnetism. The goal is to determine the electrical permittivity and magnetic permeability, which make up the usual constitutive relations:

During recent decades, a great variety of novel and complex materials have been designed with promising practical applications. It has been found, however, that in many cases the relations between electric and magnetic fields cannot be described by standard constitutive equations. Among these new materials we may mention, for example, chiral, nonreciprocal, gyrotropic and negative refractive index media. Bi-isotropic media are the most general linear, homogeneous and isotropic materials, and they respond to electromagnetic excitation according to the following relations:

where

Electromagnetic activity at optical frequencies was first observed in the XIX century, [2, 3]. Although nature presents many examples of chiral media at optical frequencies, there is, to the best of our knowledge, no natural medium able to produce electromagnetic activity at microwave frequencies. At these frequencies, the first man-made chirals were manufactured by embedding a random distribution of conducting helices with the same handedness into a host medium. The first experimental studies were developed using basic free-wave systems with the aim of verifying electromagnetic activity at microwave frequencies. Lindman [4] and later Tinoco and Freeman [5] found a critical wavelength for copper helix at

The ordered distribution of resonant structures in a periodic lattice enhances the electromagnetic effects and is able to produce larger electromagnetic activity and circular dichroism. These enhanced chirals are often known as chiral metamaterials. It has been found that a material with very high electromagnetic activity can also possess a negative refractive index [13]. During the last decade, a large variety of resonant structures have been designed and analyzed in order to produce high values of optical activity, circular dichroism and negative refractive index; for example, see [14-15] and references therein.

The presence of the chirality parameter in the constitutive equations prevents the use of standard characterization procedures. Furthermore, the existence of resonant structures, which finally produce strong variations in the characteristic parameters, makes the retrieval algorithms more complicated. In this paper, the main techniques for chiral characterization are reviewed. The second section is devoted to free-wave methods, which are the most widely used systems for chiral characterization. The guided methods presented in the third section, usually provide more accuracy than the free-wave methods, although the search for electromagnetic rotatory dispersion limits the guided systems to those of axial symmetry. In section four, resonator techniques are considered. Finally, the accuracy and feasibility of each technique is discussed.

## 2. Electromagnetic fields in chiral media

The constitutive relations (1) for reciprocal media (

where

It is possible to find a linear combination of the E and H fields, called normal modes and written *,* that satisfy independent rotational Maxwell equations so that:

These combinations of electromagnetic fields propagate through the chiral material as through an isotropic non-chiral medium, with equivalent constitutive parameters:

The normal modes satisfying eq. (4) are:

the medium parameters being:

where

### 2.1. Plane wave propagation in chiral media

A propagating wave can always be expressed as a linear combination of the propagation of normal modes. The propagation constant for each mode

If the wave is propagating in the

where

It can be deduced, from eq. (4, 5), that the vectors

## 3. Free-wave

### 3.1. Experimental systems

For regular materials where _{11} and S_{21} are sufficient to solve the problem. But there remains an extra parameter to be determined, the chirality

The usual free-wave experimental setup consists of a pair of transmitting and receiving antennas, a sample holder, a vector network analyzer (to determine magnitude and phase of the scattering parameters) and a calibration kit. In order to minimize the diffraction effects at the edges of the sample it is necessary to decrease the size of the beam or place the sample close to the antennas. But if the antenna is close to the sample, the conditions of plane wave are not accomplished and the far-field retrieval algorithms cannot be applied. Moreover, the effects of multiple reflections between the antenna and the sample are enhanced, which introduces extra errors into the measurements. The size of the beam can be reduced by placing focusing devices, such us spot focusing lens in the antennas [17] or an ellipsoidal reflector mirror after the transmitting antenna [18]. Behind the focus, the wave is not plane, but takes the configuration of a Gaussian beam. In this case, the plane wave assumption does not introduce serious errors in the measurements [19].

Figure 1 shows a schematic diagram of the experimental system developed by Varadan et al. [6, 17, 20]. It makes use of a lens in each antenna to focus the wave; the sample is located in the common focal plane. The transmitting and receiving antennas are mounted on a carriage and the distance between them can be adjusted for calibration. This setup was initially developed by Ghodgaonkar [21] to measure permittivity and permeability, and later adapted to determine rotation angle, ellipticity and power absorption. Making use of this system, Guire et al. [20] measured the normal incidence reflection of linearly polarized waves of metal-backed chiral composite samples at microwave frequencies. Ro et al. [6], determined the ellipticity, rotation angle and power absorption for samples containing helix and observed that, at resonance, the angle of rotation rapidly increases or decreases or even changes its sign, a phenomenon also known as the Cotton effect. They also reported that this phenomenon is observed in the frequency region where maximum electromagnetic power absorption occurs.

Figures 2 and 3 show a picture and a schematic diagram, respectively, of the free-wave experimental system developed by the authors. The setup is based on a previous one for permittivity and permeability measurements [22], and adapted to measure electromagnetic activity [10]. This system has been used to characterize chiral materials based on random distribution of elements and periodical lattices [10-11, 23-25].

The incident beam is focused by an ellipsoidal concave mirror, with the transmitting antenna placed at one of the mirror foci and the sample at the other. The position of the receiving antenna is not crucial in this configuration. Figure 2 shows a diagram where the focusing mirror and sample holder are not shown for clarity. The incident wave is linearly polarized in the x-axis direction. The rotation on the polarization angle,

The experimental system presented by Wang et al. [14] to characterize chiral metamaterials with high electromagnetic activity makes use of a pair of standard gain horn antennas, which are located close to the sample, without any focusing device. This configuration leads to multiple reflections and border diffractions, which affect the experimental reflection and transmission coefficients.

The scattering parameters are measured using a vector network analyzer. First, a two port calibration should be made, for which purpose the TRL (Transmission, Reflection, Line) technique is one of the most frequently used. The basic TRL calibration process is a two port calibration widely used in non-coaxial environments, such as waveguide and free-wave systems[1] -. TRL uses a 12-term error model with three calibration standards: THRU, REFLECT and LOAD. For the THRU step, the test ports are connected and the four S-parameters are measured in both directions. For the REFLECT step, a highly reflective device is connected to each port and the reflection coefficients are measured. In the LINE step, a short transmission line is inserted and, again, the four S-parameters are measured.

The reference planes of the calibration can be located at the sample surface [17] or elsewhere, for example, at the two waveguide flanges where the antennas are connected [22]. In the first case, Ghodgaonkar et al. [17, 21] implemented the TRL calibration for free-wave systems: The THRU standard was configured by keeping the distance between the two antennas equal to twice the focal distance, without the sample. The LINE standard was achieved by separating the focal plane of transmitting and receiving antennas by a distance equal to a quarter of the free space wavelength at the center of the band. The REFLECT standards for both ports were obtained by placing a metal plate in the focal planes of the transmitting and receiving antenna, respectively.

The second option is the one used by the authors. In this case, it is possible to use standards for calibration in rectangular waveguide: a short-circuit is located at the antenna flanges for the REFLECT standard and a rectangular guide of λ/4 length for the LINE standard. A second calibration is necessary to shift the reference planes to the sample surface. To do that, the reflection coefficient, with a metal plate in sample position, and the transmission coefficient without the sample were used as reference measurements. Then, the ratio between the measurement with and without sample provides the magnitude of the scattering parameters of the sample. The phase of the reflection coefficient must be corrected by a π factor introduced by metal plate. For the phase of the transmission coefficient, it is necessary to subtract a quantity equals to

At first sight, the second TRL calibration option seems more complicated; however, it has the advantage of using LINE and THRU standards in rectangular waveguide, which are more accurate and easier to use than standards obtained by moving the antennas.

A time domain (TD) transform can be used to filter out mismatches from the antennas, edge diffraction effects and unwanted reflections from antennas, sample, mirrors or elsewhere. The network analyzer makes measurements in the frequency domain, and the time domain response is obtained by applying an inverse Fourier transform. The unwanted signals, which have longer paths than the direct signal, are separated from the main signal in the TD response and then can be removed by gating the signal. In order to obtain good filtering, it is necessary to clearly identify the sample response in the time domain and to choose the most suitable filter parameters [26]. The Fourier transform of the gated time-domain response provides the filtered signal. As will be shown later, this gating process greatly improves the experimental results.

### 3.2. Retrieval algorithm

The standard retrieval procedure for a normal medium relates the transmission and reflection coefficients with the electromagnetic properties of the medium. As mentioned above, for a chiral material a new experimental quantity is necessary. Here reflection coefficient, *R*, and the transmission coefficients parallel (co-polar,

Consider a chiral slab of thickness *L*, which is illuminated by a normally incident linearly polarized plane wave

The incident wave may be decomposed as

The reflected wave in z=*0* is *L* are

The reflected and transmitted electric fields are,

where the subscript *R* (*L*) refers to right-hand (left-hand) circular polarization and subscript r to the reflected wave. The magnetic field can be expressed in a similar way in terms of wave impedance.

Applying the condition of continuity of tangential electric and magnetic fields at z=*0* and z=*L*, the reflection coefficients for RCP and LCP waves can be expressed as:

The experimental reflection coefficient is

These are related with the experimental transmission coefficients

It is useful to consider the following additional parameters:

We are now ready to compute the electromagnetic parameters. The relative impedance

The refractive index is obtained from equation (21) and is:

where *m* is an integer related to the branch index of the logarithm function. Only the imaginary part of *n* is unequivocally determined, but ambiguity arises for the real part because an infinite number of solutions exist for different choices of *m*. The value of *m* can be determined by applying Kramers-Kronig relations to the wavenumber [28] or by expanding equation (22) in a Taylor series [29].

From

The chirality parameter is obtained from

where

The retrieval procedure may fail and problems arise when the thickness *L* of the effective slab is not estimated accurately [30], when S_{11} and S_{21} are very small in magnitude [31] or when the first boundary of the effective homogeneous slab is not well determined [29]

The rotation angle

### 3.3. Results

As already mentioned, periodical structures provide certain advantages over random distributions: for example, they are easier to manufacture and are more homogeneous and they present extreme values for the rotation angle and chirality parameter. Figure 5 shows the reflection and transmission coefficients, and the retrieved parameters

A clear resonance at around 11 GHz is present in all the plots, both measured and simulated. Below the resonance, the retrieved parameters correspond to the blank FR4 board, which does not present electromagnetic activity (

There is good agreement between the retrieved parameters obtained from experiments with TD filtering and simulated data. It is important to highlight that the correct selection of TD parameters is crucial for obtaining good results in the retrieval process. The values of the time width (span) of the gating and the window shape affect the final results and, in this respect, the use of experimental and numerical data improves the calculations. Two propagation modes with two different velocities inside the chiral slab can be observed in the signal received by the receiving antenna when TD is applied, specially when the TCR coefficient is measured. Two Gaussian curves arrive at two different times, which implies a time span greater than that used for standard dielectric material.

## 4. Waveguides

### 4.1. Experimental systems

Waveguide techniques for characterizing chiral media are based on the same principles as free-wave techniques, and it is necessary to measure three complex quantities. The experimental setup, therefore, must allow measurements in, at least, two different polarizations for the transmission wave. All the experimental systems reported in the literature [7-9, 33-35] use circular waveguides because in rectangular waveguides the effects related to rotatory dispersion cannot be observed. Lubkowski [36] proposed a square waveguide for measurements with SRR (Split Ring Resonators). There, the orthogonal modes TE_{10} and TE_{01} are degenerated and, therefore, allow propagation for co-polar and cross-polar polarization. This could be a good alternative for measurement involving chiral materials but, to the best of our knowledge, such a system has not been implemented experimentally.

Figures 6 and 7 show a picture and a schematic diagram, respectively, of the experimental setup using a circular waveguide developed by the authors. It is similar to that used by Brewitt-Taylor [7] et al. and Liu et al. [34]. The circular measurement cell containing the chiral sample (2 in Figures 6, 7) is fed with a linearly polarized wave through a rectangular-to-circular waveguide transition (1 in Figures 6, 7). The transmission wave passes to the rectangular waveguide through another transition which can rotate around the longitudinal axis (3 in Figures 6, 7), enabling measurements in different polarizations. The waveguide transition contains a section of resistive film to absorb any field cross-polarized from the fields being measured, which prevents them from being reflected from the waveguide transition back on to the sample.

S-parameters are measured using a network analyzer and the TRL technique is used for calibration as described for free-wave measurements. The REFLECT standard is achieved by shorting with a metal plate the flanges where the measurement cell is connected, which defines the reference planes. The TRHU standard is achieved simply by connecting directly these flanges without the measurement cell. Finally, the LINE standard is a section of the circular waveguide quarter wavelength at midband.

Hollinger et al. [33] measured some effects of chiral materials including rotatory dispersion and circular dichroism, but they did not study the inversion procedure. Liu et al. [9, 34] used a slotted line instead of the network analyzer to measure the reflection coefficients of short-circuited and open-circuited samples, completing the experimental measurements with the rotation angle and axial ratio. Reinert et al. [8, 35] used coaxial-to-circular waveguide transition to excite the fundamental mode in circular waveguide TE_{11}, avoiding rectangular-to-circular transitions.

Two aspects that can produce important errors must be taken into account in waveguide measurements: the presence of local density variations in the sample, which causes fluctuations to the measurements, depending on sample orientation, and the existence of gaps between the sample and the waveguide walls. In order to improve the accuracy of the results, the samples are rotated and the measured scattering parameters for several positions are averaged in an attempt to overcome the sample variations and inhomogeneities [8, 9, 33, 34]. The influence of the gap depends on sample manufacturing and can be overcome, for instance, by locating the chiral particles in non-conducting spheres [8], since the multiple air gaps between the spheres are part of the sample, or by using an elastic material as dielectric host for distributing the chiral particles [12].

### 4.2. Modes in circular waveguides

The propagation of electromagnetic fields in circular waveguides filled with chiral media, also known as chirowaveguides, can be better understood by using the normal modes formulation. Assuming, as usual, solutions of the type

where

Once

Similar expressions for the magnetic components can be deduced. Assuming an azimuthal dependence as

where classification of the normal modes is described by

where

Figure 8 shows the dispersion diagram of a circular waveguide for X band (*(ϵr=3.0, μr=1.0,*

### 4.3. Retrieval algorithm

When several modes propagate in the chirowaveguide it is not possible to obtain inverse equations. Only if just one mode propagates can a simple retrieval algorithm be developed. However, this is only valid for low values of the chirality parameter.

Liu et al. [9, 34] assumed that only the fundamental mode, HE_{±11} (_{11} mode, the fundamental mode in the empty waveguide. In this case, the procedure for obtaining the inversion equations is similar to the one used in the free-wave technique. The above mentioned authors used the following experimental quantities:

The characteristic parameters of the material are calculated from the measured quantities by:

where _{11} in the empty waveguide.

These equations are relatively easy to use but are only valid for media with small values of the characteristic parameters. For example, Figure 8 shows two propagating modes in the X-Band for a chirality relative constant _{.} Using the expression for waveguide modes in the empty and filled regions and starting with approximate values for permittivity, permeability and chirality, the material properties are obtained by iterative fitting of the calculated S parameters to the experimental ones.

Brewitt-Taylor et al. [7] followed a similar procedure to calculate the characteristic parameters. More specifically, they used four modes and eight different points to perform the modal analysis in the interfaces. In this case, mode matching was performed at the two interfaces simultaneously. The number of iterations required varied from 5 to 30, although the authors admit that high chirality samples could require more iterations.

Both Reinert et al. and Brewitt-Taylor et al. established the number of modes empirically, this number increasing with the values of the characteristic parameters. For instance, the samples used by Brewitt-Taylor et al. presented maximum values of 4 for the permittivity, 0.5 for the chirality and around 1 for the permeability, while the corresponding rotation angle was always less than 20º. Samples manufactured by the authors present values of up to 7 for permittivity, and 2 for chirality and permeability (see Figure 5), finding rotation angles in the range of 90º (Figure 9). The number of modes which should be used is, probably, more than five, increasing the computation requirements.

## 5. Resonators

Some authors propose using resonator techniques to characterize chiral media [38, 39]. Such techniques are usually the most suitable method for obtaining high precision measurements. However, there are several reasons why they are difficult to use with chiral media. In particular, resonators work at discrete frequencies this is not useful to deal with dispersive media; only two real quantities (resonant frequency and quality factor) can be measured and, as mentioned above, three complex quantities have to be measured to extract permittivity, permeability and chirality. Measurements using different resonant modes and samples in several positions is a possibility, but, in this case, the experimental procedure becomes very complicated.

Tretyakov et al. [38, 39] proposed a perturbation technique, although this gives rise to a new problem. It is well known that perturbation techniques are only valid for small samples. But small samples of chiral media will have a small number of inclusions, perhaps only one, and so it makes no sense to treate them as a homogeneous material.

In our opinion, resonant techniques are not suitable for characterizing chiral media. In fact, although Tretyakov et al. studied the proposed technique theoretically, none was implemented experimentally.

## 6. Discussion and conclusions

In this work, microwave techniques for characterizing chiral media have been revised, and free-wave, waveguide and resonator techniques have been analyzed. Since chiral materials are characterized by three complex parameters, permittivity, permeability and chirality, at least, three experimental quantities have to be measured. All the techniques considered present some difficulties, while free-wave techniques are, in our opinion, the most suitable: i) They do not need any special component. ii) Quasi-plane wave propagation can be assumed in the sample, so that inversion equations for extracting characteristic parameters can be developed without difficulty. Usually, the measured quantities are the reflection coefficient, the co-polar and cross-polar transmission coefficients.

Free-wave methods present an inherent problem due to the finite size of the sample, which should be large enough to avoid diffraction at the borders. Several techniques found in the literature overcome this problem by using focusing systems, lenses or ellipsoidal mirrors. They concentrate the propagated wave at a tiny spot, where the sample is situated and so reasonably sized samples can be used. Many authors use time domain transformations to remove unwanted reflections from the experimental measurements.

Waveguide measurements may be considered a possible alternative. The measuring principle is similar to that of free-wave techniques, but they present the advantage of avoiding many kinds of external interference and use only small samples, although good machining is needed. However, they present one serious problem: usually several modes are propagated in the sample and inversion equations cannot be developed. All the techniques found in the literature use a circular waveguide because axial symmetry is necessary for measurements in different polarizations (only waves with polarization parallel to the short side can be propagated in rectangular waveguide). Unfortunately many circular waveguide components are not commercially available. A common solution is to use a rectangular waveguide and to feed the measurement cell, made in the circular waveguide, by means of two rectangular to circular waveguide transitions. The transition contains a resistive film to absorb any fields cross-polarized from those being measured. Square waveguides could, in principle, be another alternative for chiral characterization, although we have not found any theoretical nor experimental study in the current literature.

Finally, it should be emphasized that resonator techniques are not suitable for characterizing chiral media. Only one research group has proposed a resonator technique in a theoretical study without experimental implementation. In fact, resonators work at a fixed frequency and are therefore not useful for measurements in dispersive media.

## Acknowledgments

This work was supported by the Dirección General de Investigación (TEC 2010-21496-C03-02) of the Spanish Ministry of Education and Science and by Fundación Séneca (11844/PI/09) Región de Murcia (Spain).

## Notes

- Agilent Technologies Inc. Applying Error Correction to Network Analyzer Measurements, Application Note 2002, 1287-3, 5965-7709E.