## 1. Introduction

Metamaterials are artificial materials composed by the periodic arrangement of a unit cell; among the different materials in the unit cell, stratified or layered media involving metallic layers have been extensively studied and throughout this chapter, metallic material is thought in the far infrared, thus in a frequency range where metal behaves as an opaque medium, at the boundaries of which Neumann boundary condition applies. Such metallic arrays enter in the design of metallo-dielectric structures, as the artificial magnetic conductors used in the design of antenna and there is currently renewed interest in such array since they are the key piece of so-called metasurfaces (with a typical resonant behavior as in the mushroom structure of Sievenpiper [1]). In addition to be able to produce unexpected scattering of electromagnetic waves, these structures can support guided waves and it is for this property that they have been proposed by Pendry et al. [2]. In this context, these guided waves have been called “spoof plasmons” since they mimic, in the far infrared regime, the behavior of plasmons observed in the visible range (plasmons are the wave guided at the flat interface between air and metal and this requires a negative permittivity).

Because of their periodic subwavelength structuration, homogenization approaches are ideal tools to predict within a rigorous mathematical framework the scattering properties of these devices. In this chapter, we will present results coming from homogenization techniques. Our starting point is the simplest homogenization; owing to the small parameter *η = kh*, with *k* the wavenumber and *h* the typical spacing of the periodic structuration, we call simplest homogenization the homogenization performed at the leading order in *η*, which is the 0 order. Such homogenization predicts that a stratified medium can be replaced by an equivalent homogeneous and anisotropic medium associated to the usual conditions of continuity of the electromagnetic fields at the boundaries between two media; it has been used for the practical realization of several metamaterial devices, notably in the context of cloaking [3, 4]. Starting from this classical homogenization, which regards the effect of wave propagation in the bulk of the stratified medium, we will focus on two refinements that may be needed to get accurate results: (i) the case of a metallic array of small thickness for which boundary layer effects can become significant or even dominant when compared to the effect of propagation in the bulk; this case requires to conduce the homogenization up to order 1, which does not affect the equation of propagation in the bulk but makes new conditions at the boundaries to appear (these new conditions are called jump conditions), (ii) the case where the metamaterial device is composed of a succession of metallic elements arranged in a three-dimensional geometry. Although homogenization techniques could be used considering the whole device, the resolution may become tricky. Rather, a heuristic extension of the result coming from the classical homogenization is proposed, based on the analysis of the dispersion relations in each part of the whole device.

These two aspects are presented theoretically and the main theoretical results (in terms of the scattering properties or in terms of the dispersion relations of guided waves) are supported by experimental results. As a practical application, the realization of a demultiplexer is proposed, allowing for the frequency selection of an incoming large band signal into different “colored” channels.

## 2. The starting point: homogenization of metallic arrays at the leading order (order 0)

The present derivation considers first a stratified medium, with the layers made of a transparent media, afterwards the limiting case of metallic arrays in the far infrared is considered (thus, with metallic layers being associated to Neumann boundary conditions at its boundaries). Also, we restrict ourselves to polarized TM waves, for which the transverse magnetic field **e**_{3} and the whole system is invariant by translation along *x*_{3} (thus, *H* is

For a succession of layers made of transparent media, the wave equation reads

with *ω* the frequency). In Eq. (1), *H* and of

Next, Eq. (1) can be written in an equivalent form, introducing the field *C*

with boundary conditions being the continuity of *H* and of ** C**.

**n**at the boundaries between two layers (note that the field

**is linked to the electric field, see [5]).**

*C*### 2.1. The asymptotic analysis

The idea is to define three regions where different asymptotic expansions will be used, Eq. (3), with respect to the small parameter *h* the periodic spacing of the layers). The inner region contains the boundary between the stratified medium and the air. The two outer regions for *x*_{1} > 0 and *x*_{1} < 0 are the regions far enough the interface, where the evanescent field can be neglected, while the inner region contains the evanescent field. Next, the inner region and the outer regions are connected using so-called matching conditions, which will constitute the boundary conditions for the outer solutions (see **Figure 1**).

#### 2.1.1. The outer and inner expansions

The asymptotic expansions are thought with spatial dependences on a *macroscopic* coordinate **x** associated with low variations of the fields (with the typical scale 1/*k* of the wave) and a *microscopic* coordinate **y** associated with rapid variations (with the typical scale *h* of the layers) and in each region, we keep the coordinates that are relevant to describe the variations of the field. To do so and with *H*, ** C**) can be expanded by using the following asymptotic expansions

(3) |

with the outer terms (*H*^{n}, *C*^{n}) for *x*_{1} < 0 and the inner terms (*h*^{n}, *c*^{n}) being periodic *w.r.t*. *y*_{2}. The differential operator reads, in the different regions, as

#### 2.1.2. The matching conditions

The inner and outer problems have to be associated with boundary conditions or radiation conditions which ensure that the problems are well-posed. Since the outer expansions hold true only far away from the interface, the outer terms do not have to satisfy the continuity conditions at *x*_{1} = 0. Reversely, the conditions at infinity satisfied by the inner terms are unknown *a priori*. These missing conditions are provided by the matching conditions, which read at leading order

### 2.2. The homogenized wave equation at the leading order

We want to establish the wave equation satisfied by the mean fields *x*_{1} < 0, where we have defined the average over *x*_{1} < 0, only. For *x*_{1} > 0, the wave equation is obviously

being the same at each order and the fields equal their averages. Eq. (2), at leading order (1/*η*), read

Now, we establish the relation between *H*^{0} (this latter being equal to its average). Eq. (2) at order *η*^{0} in the outer problems *x*_{1} < 0 give

and

Averaging both equations, with *H*^{1} and of *y*_{2} (thus,

### 2.3. The continuity relations at the leading order

To the homogenized wave equation (10), we have to associate continuity (or discontinuity) conditions at the interface *x*_{1} = 0. To that aim, we have to consider the inner solution and its matching with the two outer solutions. We are looking for the quantities *η*), that *h*^{0} does not depend on **y**. From the previous section, we already know that *H*^{0}(**x**) does not depend on *y*_{2}, from which the matching conditions, Eqs. (5a) and (5b), give

Next, Eq. (2) in the inner region gives also, at the leading order,

(we have used the periodicity of *c*^{0} w.r.t. *y*_{2}). From the matching conditions Eqs. (5c) and (5d) integrated over *Y*, we get

At the first order, the usual continuities of the electromagnetic fields are found.

### 2.4. The homogenized problem for metallic layers in the far infra red

For rigid layers in the air and denoting *φ* the filling fraction of air within the stratified medium, we have: in the air

(14) |

It is worth noting that the above continuity relation means notably that (i) at the interface with the air,

## 3. Weakness of the first-order homogenization for small thickness devices

In this section, we inspect the validity of the homogenized problem, Eq. (14). This is done in the configuration of **Figure 2**, where a metallic array is placed in the air and illuminated by a plane wave at oblique incidence *θ*. We measured experimentally the transmission coefficient *T*^{exp} for the array of thickness *e* (and the array structuration is characterized by *φ* and *h*) and compared *T*^{exp} to the transmission coefficient *T* obtained in the homogenized problem. To anticipate, we shall see that the leading order homogenization may become unsatisfactory for small thicknesses of the array and going up to the homogenization at order 1 is necessary. Thus, the derivation of the scattering coefficients in these homogenized problems is presented first, afterwards comparisons with those measured experimentally are presented.

### 3.1. The scattering coefficients in the homogenized problems

Let us start with the first-order homogenization, for which the explicit solution of Eq. (14) read

(15) |

We used the dispersion relations coming from the wave equations

where *z*_{i},

and *H* and ** C.n** as boundaries conditions, jump conditions are obtained. Specifically, the homogenization at order 1 read

In the above equations, it appears that both *H* and ** C**.

**n**are now discontinuous (the obtained conditions are jump conditions at the boundaries of the stratified medium). This is why

*H*

^{±}(same for

**) are defined, as the limit values of**

*C**H*at the boundary. In this case, the expressions of (

*R*,

*T*) in Eq. (16) are still valid, but we get

(19) |

The parameters

In principle, (*H*, ** C**) in the homogenization at order 0, Eq. (14), approximate the solution of the actual problem up to

*O*(

*η*) and (

*H*,

**) in the homogenization at order 1, Eq. (18), approximate the solution of the actual problem up to**

*C**O*(

*η*

^{2}). Thus, we could expect that the difference between both remains incidental. We will see that it is not the case. The reason is that the jump conditions obtained at order 1 encapsulate the effect of the boundary layers at the entrance and at the exit of the stratified slab and these effects may become dominant compared to the effect of the propagation in the bulk of the slab. This is what we inspect below.

### 3.2. Experimental measurements of the scattering coefficients for varying slab thicknesses e

To test the ability of the leading order homogenization to capture the scattering properties of a metallic array, we realized six arrays of different thicknesses *e* = 30 μm and 0.25, 1, 4, 14 and 20 mm. Otherwise, *h* = 6 mm and *T*^{exp} using two X-band frequency horn antennas at both extremities of an electromagnetic chamber (**Figure 3**). The arrays were amounted on a plate able to rotate in order to realize varying incidence angles *θ* (and the incidence wave is polarized in TM polarization).

We start by reporting in **Figure 4** the spectra *e* = 30 μm and *e* = 20 mm (right panels). For frequencies in the range [8,12] GHz, the thinnest array realizes *e* = 30 μm, the leading order homogenization predict *T*| given by the homogenization at order 1 and this latter appears to be in good agreement with *f* and *θ*).

For the thicker array, the transmission predicted by the homogenization at order 0 is closer to the measured transmission; in this case, although the spectra obtained using the homogenization at order 1 reproduces better the form of the measured spectra, the relative errors are in both cases about 15% (these highest discrepancies are due to highest relative errors for transmissions close to 0).

Next, we inspect the variations of the transmission coefficient as a function of *ke* for the six arrays, **Figure 5**. Each plot corresponds to a fixed frequency (*f* = 8, 10 and 12 GHz) and we reported

It is visible that thin arrays are not correctly described by the homogenization at order 0. More specifically, it largely underestimates the scattering strength of thin arrays (small *e* produces systematic large errors in the prediction) and it becomes accurate only when *ke* > 1. To the contrary, the homogenization at order 1 is able to describe the scattering strength of thin and thick arrays.

In conclusion, the homogenization at leading order is valid for *ke* < 1 and *kh* < 1 and thus helpful to predict the behavior of metallic arrays as used in the design of many metamaterial devices (see also [6]). However, care has been taken when using arrays with vanishing thicknesses, typically *ke* < 1. In such cases, the homogenization at order 1 has to be considered. As previously said, this is because the leading order homogenization does not account for the scattering effects of the wave at the entrance and at the exit of the array (these are boundary layer effects) which are correctly accounted for in the homogenization at order 1 through the—boundary—parameters

## 4. Use of the results coming for the first-order homogenization for metamaterials with complex geometry

In this section, we investigate another configuration of wave propagation where predictions provided by the homogenization are useful. We consider the ability of periodic structures to present a band structure, with frequency ranges where the wave propagation is forbidden (band gaps) and frequency ranges where the propagation of guided waves is allowed (passbands). We study a structure based on such principle. First, the dispersion relation of guided waves within a periodic media (made of metallic plates or metallic rods) is discussed and it is shown that the leading order homogenization is able to reproduce the real dispersion relation. Next, a system of waveguides is thought, based on the band structures of the periodic media, in the waveguide and in the surrounding medium. The application of demultiplexing RF signals is proposed and experimentally validated.

### 4.1. Dispersion relation of waves guided in a periodic medium

Let us start with the derivation of the wave guided in a 2D stratified medium (the famous “spoof plasmon,” **Figure 6a**). In the homogenized problem, Eq. (14), this wave corresponds to the solution of the homogeneous problem (solution in the absence of source), for which the solution reads

In the above expression, the field *H* for

with *φ* the filling fraction of air and *e* the length of the layers in the stratified medium. This dispersion relation has been established previously using approximate modal method [2] and it is easily obtained by considering the equivalent homogenized problem.

It is worth noting that such guided wave propagates in the homogenized medium described by the wave equation

(according to Eq. (14)) and basically, the wave equation (23) tells us that the wave inside the grooves can only propagate along one direction (the *x*_{1}-direction). While in principle the upper frequency

Now, we want to go toward a 3D structuration (**Figure 6b**), where rods are considered, with radius *r*, periodic spacing *h* along *x*_{2} and *x*_{3} and height *e*. The extension of the homogenization results (Eq. (14)) to three dimensions is easy and we find that a periodic structuration of rods produce an equivalent transverse isotropic medium, with the axis of anisotropy along *e*_{1} (the two other directions are equivalent, now *φ* the volume fraction of air in the rods, whence *e*_{s} (the surrounding rods form the surrounding medium SM) and imagine that we work in a frequency range such that

that is in the band gap of the surrounding medium, from which *H* = 0 at *e*. Looking for the existence of guided wave in the waveguide, we extend the homogenization result to this 3D configuration assuming a solution of the form

In the simple form thought above, we added heuristically dependence in the *x*_{3} direction which accounts for Dirichlet boundary conditions at

Otherwise, Eq. (26) accounts for the continuity of *H* at the interface with the air and for the Neumann boundary condition at the ground plane; as previously, the additional condition

In Eq. (27), the band gap above

The existence of the resulting finite pass band operating in the frequency range

this frequency range and we refer to *n* an integer, a waveguide obtained when *n* lines of rods have been shortened (with resulting height *e*) with respect to the rods of height *e*_{s} forming the surrounding medium (SM); the *w*_{n} waveguide has a width

### 4.2. Experimental validation of the homogenized dispersion relation

To begin with, we validate experimentally the existence of the pass band and check the validity of our predictions of the FR with bounds *h* = 7 mm, *r* = 3 mm and *e*_{s} = 30 mm. The shortened rods have *e* = 17 mm.

First, we report measurements of the electric field in the range [2.1– 4] GHz (**Figure 7** for *w*_{1}). This has been done in a semianechoic chamber using an Agilent 8722ES network analyzer; an *S*-band coaxial-to-waveguide transition has been used as an excitation source and an electric near-field probe mounted on a motorized two-dimensional scanning system has been used to measure the field distribution at a distance of about 1 mm above the structure. The wave guided within the waveguide *w*_{1} is visible at *f* = 3.7, 3.8 and 3.9 GHz, as expected from the frequency range FR(*w*_{1}) = [3.6–4] GHz (see **Table 1**) .

To go further, we report the transmission in the *w*_{1} and *w*_{3} waveguides as a function of the frequency (**Figure 8**). This has been done by placing a second coaxial-to-waveguide transition at the end of the waveguide and by implementing a normalization to the free air transmission between the emitter and the receptor. In both cases, the existence of a finite pass band is confirmed (the waveguide is called colored) and the observed bounds of the FR are in good agreement with our predictions (the theoretical dispersion relations are reported, with FR(*w*_{3}) = [3.1–4]). One can notice here the importance of the attenuation for the spoof-plasmon modes providing the smallest wavelengths. These wavelengths correspond to the highest frequencies in the transmission bands. This phenomenon, known in plasmonics, happens due to intrinsic losses in the considered materials for the high-wave-vector components.

### 4.3. Application to the design of a demultiplexer

This validation being performed, a multichannel demultiplexer is easy to design; the principle of the demultiplexing is shown in **Figure 5**. A main waveguide, called white guide, is built in order that the FR(white) covers the working frequency range; this is done by choosing: (i) *w* large enough to produce a small enough *e* small enough to produce a large enough *e* = 15 mm, we expect FR(white) [2.5–4.5] GHz.

Next, three-colored waveguides (red, green and blue) are thought in order to support guided mode propagation in three different frequency ranges with no overlapping (see **Figure 9**). Again from Eqs. (23) and (27), thin FR are obtained for small *w* and we set *L* = 17, 19 and 21 mm, we expect this condition to be fulfilled with

The efficiency of the demultiplexer has been tested experimentally and it is illustrated in **Figure 10**. The frequency selection of the colored channel are visible, with the red channel being active for

## 5. Conclusion

The design of devices made of subwavelength periodic structure can be helped using the various theoretical predictions provided by homogenization approaches. Among the different homogenization technique, the one presented in this chapter has the advantage to be developed within a rigorous mathematical framework and it yields predictions in a deductive way, that is without any adjustable parameters. We illustrate the ability of such techniques for the scattering properties of metamaterial devices (underlying the limit of the simplest homogenization) and for the ability of certain metamaterial devices to support guided waves of “spoof plasmon” type.

As indicated throughout this chapter, alternative forms of homogenizations can be used, which are more adapted to thin devices. These theoretical tools can be used in order to realize the control of light propagation in a desired way.